Volume 525, January 2011
|Number of page(s)||8|
|Published online||10 December 2010|
Modeling of the Vela complex including the Vela supernova remnant, the binary system γ2 Velorum, and the Gum nebula
Humboldt Universität zu Berlin, Institut für Physik,
2 National Taras Shevchenko University of Kyiv, Department of Physics, Kyiv, Ukraine
3 National Taras Shevchenko University of Kyiv, Astronomical Observatory, Kyiv, Ukraine
4 ISDC, Versoix, Switzerland
Accepted: 18 October 2010
We study the geometry and dynamics of the Vela complex including the Vela supernova remnant (SNR), the binary system γ2 Velorum and the Gum nebula. We show that the Vela SNR belongs to a subclass of non-Sedov adiabatic remnants in a cloudy interstellar medium (ISM), the dynamics of which is determined by the heating and evaporation of ISM clouds. We explain observable characteristics of the Vela SNR with a SN explosion with energy 1.4 × 1050 erg near the step-like boundary of the ISM with low intercloud densities (~10-3 cm-3) and with a volume-averaged density of clouds evaporated by shock in the north-east (NE) part about four times higher than the one in the south-west (SW) part. The observed asymmetry between the NE and SW parts of the Vela SNR could be explained by the presence of a stellar wind bubble (SWB) blown by the nearest-to-the Earth Wolf-Rayet (WR) star in the γ2 Velorum system. We show that the size and kinematics of γ2 Velorum SWB agree with predictions of numerical calculations for the evolution of the SWB of Mini = 35 M⊙ star. The low initial mass of the WR star in γ2 Velorum implies that the luminosity of the nuclear line of 26Al, produced by γ2 Velorum, is below the sensitivity of existing gamma-ray telescopes.
Key words: ISM: supernova remnants / ISM: clouds / ISM: bubbles / ISM: individual objects: Vela SNR / shock waves / stars: Wolf-Rayet
© ESO, 2010
The Vela complex is one of the most interesting regions in the galactic plane. The observable flux from the Vela region ranges from radio to TeV energies. It consists of many objects, including the Gum nebula, the Vela supernova remnant (SNR) and Vela Jr. SNR (SNR RX J0852.0-4622 superposed on Vela), the binary system γ2 Velorum, the IRAS Vela Shell, an OB-association, ζ Pup, etc. Some of them are shown in Fig. 2. Because of the distances to these objects it is possible that some of them can intersect not only in projection, but also physically. We propose a scenario of joint evolutionary interaction of the Vela SNR, the binary system γ2 Velorum, the IRAS Vela Shell, the Vela OB2-association, and the Gum nebula.
The Vela SNR is one of the closest supernova remnants to us. Different estimates of the distance to the Vela SNR suffer from a large uncertainty: from 250 ± 30 pc (Cha et al. 1999) to 350 pc (Dubner et al. 1998 and references therein). Hubble Space telescope parallax observations of the Vela pulsar give the distance to the pulsar of pc (Caraveo et al. 2001), and the best estimate is from the VLBI parallax measure (Dodson et al. 2003): (1)Similarly uncertain are the estimates of the Vela SNR age, which range from a few thousand years (Stothers 1980) to tSNR ≃ 2.9 × 104 yr (Aschenbach et al. 1995). The most commonly cited estimate is tSNR = tpulsar ≃ 1.14 × 104 yr, where tpulsar is the age of the Vela pulsar (Reichley et al. 1970). The total 0.1−2.4 keV X-ray luminosity from Vela SNR in erg/s is (Lu & Aschenbach 2000) (2)The main peculiarity of Vela SNR is the difference in the X-ray brightness and radius of its south-west (SW) and north-east (NE) parts. The ROSAT All-Sky Survey image (Fig. 1) of the Vela SNR reveals a shell with a diameter of about 8° (Aschenbach et al. 1995), which implies a mean linear diameter of (3)The SW part of the shell appears to have a radius larger by a factor of (4)than the NE part.
ROSAT All-Sky Survey image (0.1–2.4 KeV) of the Vela SNR (Aschenbach et al. 1995). A–F are extended features outside the boundary of the remnant (“bullets”). Light blue to white contrast represents a contrast in surface brightness of a factor of 500 (Aschenbach et al. 1995). Blue curves show the NE and SW hemispheres of the Vela SNR. The yellow curve shows the contour of the SWB of γ2 Velorum.
Apart from the difference of the radii, the spatially-resolved spectroscopic analysis by Lu & Aschenbach (2000) shows that the SW part of the shell appears to be hotter than the NE one. The shell is bright only on the NE side, while the SW side appears to be dim and is apparently more extended (see Fig. 1). The boundary between the bright and the dim part of shell is quite sharp (Lu & Aschenbach 2000). The change in the properties of the shell at different sides indicates that the characteristics of the ISM in the NE part differ from those in the SW part. The sharpness of the boundary shows that the change in the ISM properties is abrupt rather than gradual.
Lu & Aschenbach (2000) estimate the contrast change of the two faint regions, each 1.5° × 1.5°, in the NE and SW parts, to be a factor of ≃11 in brightness and ~6 in emission measure. Assuming the same estimate for the emission measure contrast change for the entire NE and SW parts of the SNR and taking into account the difference of the radii of the two parts of the SNR, one can find that the luminosity of the NE part Lx,NE is a factor of ~3.7 higher than the luminosity Lx,SW of the SW part.
Apart from the difference in the overall luminosity and temperature, the characteristics of the X-ray emission from the NE and SW parts of the SNR also reveal a difference in the column density of the absorbing material along the line of sight (Lu & Aschenbach 2000). The absorption column densities NH range from 5.0 × 1019 cm-2 in the NE part to 6.0 × 1020 cm-2 in the SW part.
Finally, the Vela SNR is peculiar in still another aspect: the main shock of the SNR is not observed. Instead, the bulk of the X-ray emission is distributed all over the SNR volume. Such an observational appearance can be caused by the SNR expanding into a highly inhomogeneous (“cloudy”) medium. In this case, the main shock advances through a low-density interstellar medium (ISM), leaving behind denser clouds, which are subsequently heated and partially evaporated by thermal conductivity and transmitted shocks. This results in the appearance of a distributed emission throughout the entire volume of the SNR, instead of from a thin shell at the interface of the main shock with the low-density ISM. The main X-ray emitters in the remnant are the two (cool and hot) phases (components) of heated cloud matter (Lu & Aschenbach 2000; Miceli et al. 2005, 2006).
Despite extensive investigations, there is no common agreement about the evolutionary status of the Vela SNR. In early studies, a distance of 500 pc was assumed, which corresponds to a rather large size of 70 pc in diameter. Together with this assumption, the absence of the X-ray limb-brightening effect and of the clear signature of the high-speed shock in the SNR suggested a radiative stage of the Vela SNR evolution. The appearance of filament structures in optic and radio waves as well as 100 km s-1 absorption lines in the spectra of background stars additionally supported this scenario (see, e.g. Gvaramadze (1999) and references therein). But recent observational data strongly support the hypothesis of the Vela SNR being in the adiabatic stage. New results for the distance place the Vela SNR at D ≈ 290 pc, implying that the mean radius of the Vela SNR is R ≈ 20 pc, and the dynamics of bullets (ejecta fragments) outside the SNR boundary suggest the expansion of the SNR in the low-density ISM. The presence of the Vela pulsar indicates that the Vela SN progenitor was of Mini ≤ 25 M⊙. For Mini = 11−25 M⊙ presupernova masses are Mfin = 10.6−16.6 M⊙ and masses of the ejecta Mej = 9−15 M⊙. (Limongi & Chieffi 2006; Eldridge & Tout 2004; Kasen & Woosley et al. 2009). Therefore, the interaction (merger) of the massive SN ejecta with a velocity of over 1000 km s-1 with the red supergiant (RSG) shell occurs in the adiabatic regime, and the Vela SNR with a radius of 20 pc and an age of about 10 000 years should be in the adiabatic stage of evolution with characteristic velocities of about mRSNR/tSNR ≈ 1000 km s-1 (m ≤ 1 in the free expansion case and m = 0,4 in the adiabatic one) without forming a thin dense radiative shell. Another observational confirmation of the absence of the 100 km s-1 shell follows from the studies of absorption lines in the spectra of background stars (Cha & Sembach 2000), where all stars in the Vela SNR direction at distances smaller than 350 pc do not show evidences of a 100 km s-1 absorption line. The maximum broadening corresponds to <50 km s-1. Only stars with distances exceeding 500 pc show 100 km s-1 features. Radio and optical shells (shell-like filaments) can be naturally explained by the emission of filamentary structures, exited by the SNR shock (akin to the Cyg Loop case, where filaments coexist with the main fast adiabatic shock).
The hydrodynamical model of the Vela SNR evolution cannot be described directly by the Sedov solution (Sedov 1959), because the expansion proceeds in a cloudy rather than homogeneous ISM. Instead, to describe the Vela SNR evolution one can use the White & Long solution, which describes the evolution of a supernova remnant expanding into a cloudy ISM (White & Long 1991).
As explained above, the observed asymmetry of the Vela SNR is most probably due to the difference in the properties of the ISM on different sides of the remnant, in particular, by the different densities of clouds or different density contrast between the clouds and the IC medium. Below we assume that before the SN explosion the average number density (concentration) and, therefore, volume filling factor of clouds in NE part was larger than that in the SW part. This naturally explains the smaller radius and the lower temperature of the NE part of the remnant. Once more, interaction of shock wave with clouds results in numerous filamentary structures of disrupted cloud material, visible in radio band, and the radio image of the Vela SNR really shows more numerous shell structures in the NE part than in SW one (Bock et al. 1998). The NE-SW asymmetry proposed here can also explain a pulsar wind nebula (PWN) displacement to the SW according to the pulsar position. Namely, the reverse shock from more dense NE part reaches and destroys PWN earlier (Blondin et al. 2001; LaMassa et al. 2008)
The conjecture that the densities on the two sides of the SNR are different is supported also by the dynamics of the shrapnel (bullets) – high-velocity clumps of the SN ejecta with an overabundance of heavy elements (Aschenbach et al. 1995; Miceli et al. 2008; LaMassa et al. 2008; Yamaguchi & Katsuda 2009). X-ray spectra of shrapnel argue a low density of the ISM around the Vela SNR, whereas a difference in the distances traveled by the protruding shrapnel on the NE (shrapnel pieces A-D/D’) and on SW (shrapnel pieces E and F) sides of the SNR, as it is illustrated for the shrapnel pieces D and E shown in Fig. 1, is consistent with a corresponding difference in mean ambient density, assuming that each shrapnel piece has approximately the same density and initial velocity.
The solution for the shock radius rs of the remnant as a function of the age t has the form (White & Long 1991) (5)where E is the explosion energy, ρic is the density of the intercloud medium and γ = 5/3 is the adiabatic index. The solution (5) differs from the standard Sedov solution without evaporating clouds by a different choice of a phenomenological constant K.
The constant K depends on two parameters: on the ratio C of the fraction of the ISM mass initially contained in clouds and evaporated behind the shock wave to the ISM mass in the intercloud medium (6)where ⟨ρc⟩ is the volume-averaged density of the clouds, and on the ratio of the cloud evaporation time scale tev to the SNR age t(7)To express the dependence of K on these two parameters we can use a simple analytical approximation (8)and K/KS ~ C-1 for C ≪ 1,τ ≫ 1, where KS = 1.528 is the value of constant K for the Sedov solution (when C = 0). From Eq. (5) follows the value of the shock velocity (9)which allows to calculate the temperature at the shock (10)(here μ = 16/27 is the average mass per particle in hydrogen mass units mH for helium/hydrogen abundance ratio nHe = (1/12)nH and k is the Boltzmann constant) for the White & Long (1991) solution.
The dynamics of expansion of the SNR is fully characterized by a set of four parameters: the explosion energy E, the preshock intercloud ISM density ρic, the cloud/intercloud density ratio C and the evaporation time/SNR age ratio τ. The values of these parameters can be derived from the set of the observed characteristics of Vela SNR, such as the shock radii, characteristics of the X-ray radiation for the both NE and SW parts, etc.
For the NE and SW shock radii RNE ≃ 18 pc and RSW ≃ 23 pc and SNR age tVela ≃ 1.14 × 104 yr from Eq. (9) we obtain corresponding shock velocities (11)(12)and from Eq. (10) shock temperatures , . Average (emission measure weighted) temperature of plasma inside the SNR is about twice as high (White & Long (1991)): (13)(14)According to the results of Lu & Aschenbach (2000) X-ray radiating plasma in both parts of the Vela SNR consists of two phases: a hot one (T ~ 0.5−1.2 keV) and a cool one (T ~ 0.09−0.25 keV), and the cool one dominates the X-ray luminosity of the SNR. The temperatures found above correspond to the hot component in the two-temperature Raymond-Smith thermal plasma model used in Lu & Aschenbach 2000). The hot evaporated gas component with the volume filling factor fhot ≃ 1 dominates the shock dynamics, while the cool one with fcool = 1 − fhot ≪ 1 dominates in X-ray radiation. The role of the initial intercloud interstellar gas is negligible in both shock dynamics and X-ray radiation. It means that (15)where Chot = ⟨ ρc,hot ⟩ /ρic, and hereafter we take τ ≪ 1 and use the approximation . In this case, Eq. (5) for shock radius is reduced to (16)where nhot = ρc,hot/mH is the nucleon number density of the hot component. As we can see, Eq. (16) is the Sedov solution, in which a mean density of intercloud plasma inside the remnant ρic is replaced by the mean density of the evaporated clouds (of the dominant hot component ρc,hot in our case).
The nucleon number density of the hot component nhot can be estimated from its X-ray radiation. The X-ray luminosity of the SNR is an integral over SNR volume for X-ray emissivity ϵX = nenHΛX(T) of plasma with temperature T, electron (hydrogen) number density ne(nH) and cooling function ΛX(T) (17)Within the model of White & Long (1991), LX is estimated in erg/s as (18)where nic = ρic/mH is the intercloud ISM nucleon number density (nucleon number density n = (4/3)nH )), Λ-22 is the cooling function in units 10-22 erg cm3 s-1, Q is a number on the order of one and depends on C and τ (for C ≫ 1 and τ ≪ 1 Q ≃ 1).
Unfortunately, we do not know the X-ray luminosity of the hot component of Vela SNR with satisfactory accuracy. Therefore, we use the more accurate data about the emission measure of the hot plasma for the calculation of the nucleon number density of the hot component. According to the observation of Lu & Aschenbach (2000), the mean value of the emission measure A of ϕpixel × ϕpixel = 8.75′ × 8.75′ pixel (19)of the hot plasma in NE region is . Meanwhile, from Eq. (19) and conditions nenH = 0.66n2, it follows that (20)and, taking into account that , (21)Now, from Eq. (16), we can find nucleon number density of hot component in SW part of Vela SNR (22)and the energy of the Vela SNR explosion (23)The cool component of the X-ray radiating plasma in the Vela SNR does not influence the shock dynamics, but dominates in the X-ray luminosity. The mean value of the emission measure A and the temperature of the cool plasma in the NE region is , and in the SW region , (Lu & Aschenbach 2000). From Eq. (20) for the NE part of the Vela SNR follows (24)The filling factors of hot and cool plasma in the NE part can be estimated from the equality of the pressure P ∝ nT of both components: (ncoolTcool)NE ≃ (nhotThot)NE(25)or , . And, finally, from Eq. (20) follow nucleon number densities , .
Similarly, for the SW part of the Vela SNR , , , and , .
Table 1 gives a summary of the parameters of the NE and SW parts of the SNR derived from the X-ray data.
As a test for the self-consistency of our model we can calculate the predicted X-ray luminosities of different parts/components of Vela SNR in erg/s, using Eq. (17): (26)where we use the approximation validated by White & Long (1991) Λ-22(T) = 1. Similarly, (27)(28)(29)The total model luminosity LX, tot ≃ 3.2 × 1035Λ-22 erg s-1 is close to the observed one of Eq. (2)
Characteristics of Vela SNR.
Estimates of the physical parameters of the Vela SNR derived in the previous section show that it is evolving in the inhomogeneous medium whose density changes in a step-like manner. In this section we explore the possibility that this step-like change of the ISM parameters can be related to the presence of the boundary of the stellar wind bubble (SWB) around a Wolf-Rayet (WR) star in the γ2 Velorum system, which is situated in the vicinity of the Vela SNR.
γ2 Velorum is a WC8+O8-8.5III binary system. The WR component in this binary system is the closest to Earth WR star (WR11). The estimates of the distance to γ2 Velorum have evolved over the years and are still controversial. An early estimate of the distance based on HIPPARCOS parallax was pc (Schaerer et al. 1997). It was recently revised by several scientific groups. Millour et al. (2007) give an interferometric estimate of the distance pc. North et al. (2007) estimate the distance to γ2 Velorum to be based on the orbital solution for the γ2 Velorum binary obtained from the interferometric data. Finally, a revision of the analysis of HIPPARCOS data gives a distance of van Leeuwen (2007). Therefore, in the following we take the distance to γ2 Velorum to be Dγ2Vel ≃ 330pc.
The current mass estimate of the WR star is MWR = 9.0 ± 0.6 M⊙ (North et al. 2007). The mass of the O star is MO = (28.5 ± 1.1) M⊙ (North et al. 2007). The γ2 Velorum system is an important source from the viewpoint of nuclear gamma-ray astronomy. Because it is the nearest WR star, this is the only source that can potentially be detected as a point source of 1.8 MeV gamma-ray line emission from the radioactive 26Al with current generation instruments. Previous observations by COMPTEL put an upper limit at the level of 1.1 × 10-5 γ cm-2s-1 on the line flux from the source (Oberlack et al. 2000), which is, apparently, below the typical predictions of the models for stars with the initial mass Mini ~ 60 M⊙. One should note, however, that modeling of 26Al production suffers from uncertainties of the nuclear reaction cross-sections, stellar parameters (such as rotation and metallicity), etc. (Limongi & Chieffi 2006; Palacios et al. 2005).
The estimate of the initial mass of the WR star, Mini ~ (57 ± 15) M⊙, is obtained from the evolutionary models of isolated rather than binary stars (Schaerer et al. 1997). It is possible that within a binary system, the mass transfer between the companions can change the stellar structure of both components, so that, for example, stars with initial masses as low as 20M⊙ can become WR stars in binary systems (Vanbeveren 1991). Modeling the binary evolution of γ2 Velorum leads to a lower limit on the initial mass of the WR star: MWR,i ≥ 38 M⊙ in (Vanbeveren et al. 1998) and MWR,i ≃ 35 M⊙, MO,i ≃ 31.5 M⊙ in (Eldridge 2009).
WR stars are expected to be surrounded by the multi-parsec scale bubbles blown by the strong stellar wind inside a photoionized HII region. The size of the HII region and bubble depends on the (time-dependent) joint action of the photon luminosity in the Lyman continuum LLyC and the mass-loss rate of the star, Ṁw, on its age, t, on the wind velocity vw and on the density of the ambient medium into which the bubble expands, ρ0. Qualitatively, the relation between the bubble radius Rbub and the above parameters can be found using the analytical calculation of Weaver et al. (1977)(30)where κ is a numerical coefficient which, in the simple case considered by Weaver et al. (1977) was κ = [125/(154π)]1/5, while in a realistic case it can be found from numerical modeling based on a stellar evolution model. For example, a numerical model considered by Arthur (2007), which studied the evolution of a star with initial mass 40 M⊙ in a medium of density n0 = 15 cm-3, predicts the final radius of the shock in the ambient interstellar medium Rsh ≃ 36 pc at the end of WR phase before SN explosion. Using the Eq. (30) one can re-scale the numerical simulations for the particular values of ISM density to find that for the typical value n0 = 0.1 cm-3 in the ISM the radius of the SWB can reach 100 pc. Meanwhile, as we will see later, γ2 Velorum system was born in a molecular cloud and its bubble should be considerably smaller.
Comparing the distances of Vela SNR (~290 pc) and γ2 Velorum system (~330 pc) and taking into account the proximity of the two objects on the sky, one can notice that if the SWB around the γ2 Velorum system has indeed a radius of 30−70 pc (≥5° on the sky), the Vela SNR, which itself has a size of ~40 pc, is expected to physically intersect with the γ2 Velorum bubble. In view of this geometrical argument it is natural to ascribe the observed step-like change in the parameters of the ISM at the location of the Vela SNR to the boundary of the SWB of γ2 Velorum.
The hypothesis of intersection between the Vela SNR and γ2 Velorum SWB is further supported by the simple geometrical form of the boundary between the bright and the dim part of the SNR shell. Indeed, the boundary roughly follows the contour of an ellipse whose major axis is perpendicular to the direction from the center of Vela SNR toward the γ2 Velorum, so that the minor axis is aligned with the direction toward γ2 Velorum, see Fig. 1. The projected distance between the Vela SNR and γ2 Velorum is D ′ = 5.2°. Assuming the distances DVela ≃ 290pc and Dγ2Vel ≃ 330pc one finds that the physical distance between the two objects is (31)which we adopt as an estimate for the radius of the SWB around γ2 Velorum.
Within the geometrical model discussed above, the observed difference in the absorption column density NH between the NE and SW parts of the Vela SNR can be used to estimate the total mass of the ISM swept up by the stellar wind of γ2 Velorum over the entire lifetime of the SWB. Taking the difference between the measured NH values in the NE part and SW parts (Lu & Aschenbach 2000) (32)one can find the total mass of the γ2 Velorum SWB: (33)Assuming that the bubble has expanded into a homogeneous ISM over the entire expansion history, one would estimate the ISM density around γ2 Velorum as
(34)One could notice that this estimate of the density of the ISM is much higher than the estimates of both the intercloud and of the volume-averaged cloud density around Vela SNR (Table 1). We come back to the discussion of the origin of this discrepancy below.
The angular size of the γ2 Velorum SWB, found from our interpretation of the asymmetry of the Vela SNR, i.e., from the ratio of Rγ2Vel ≃ 44pc and the distance 330pc as ≃7.6° coincides with the angular size of a large circular arc like structure, centered on (l,b) = (263°, −7°) with the radius RIVS ≃ 7.5°, which is visible in the infrared band, known as the “IRAS Vela shell” (IVS) (Sahu 1992). This structure surrounds the Vela OB2 association including γ2 Velorum and ζ Puppis. A recent study of the spatial distribution of the neutral hydrogen and radio continuum emission at 1420 MHz of the IVS by Testori et al. (2006) provides a new estimate of the coordinates of the centroid of IVS from the observed IR emission (l,b) = (259.9°, −8.3°) and the radius of the neutral hydrogen shell in the SW sector between position angles ~162° and ~265°RIVS ≃ 5.7°. Testori et al. (2006) also estimate the mass of ionized and atomic components of the shell, assuming a distance of 400 pc: MIVS ≃ 9.1 × 104 M⊙ (or MIVS ≃ 6.0 × 104 M⊙ for our distance 330pc). It is expected that the amount of molecular gas in the IVS is about ~105 M⊙(Rajagopal & Srinivasan 1998).
Adopting this interpretation, one can find that the estimates of the parameters of the SWB of γ2 Velorum derived above from the analysis of X-ray data on the Vela SNR, agree well with those found from the analysis of the IVS data. In particular, adopting an estimate of the distance toward γ2 Velorum for the center of the IVS, one finds the radius (35)and the total mass of the shell (36)comparable to the estimates found in Eqs. (31) and (33), respectively.
The initial mass of the WR11 star Mini is equal the sum of the current mass of the star MWR = 9.0 M⊙ and the mass which was blown by the wind: (37)where ṀWR is the mass loss rate and tWR is the time interval of the WR wind. The sum in the above equation includes mass loss via different types of winds ejected by the star at different stages of stellar evolution: the main sequence stage (MSS) wind, the red supergiant stage (RSGS) (for stars with initial mass Mini ≤ 40 M⊙) or the luminous blue variable stage (LBVS) (for Mini ≥ 40 M⊙) wind, and the continuing WR wind. Numerical calculations from Freyer et al. (2003), Freyer et al. (2006), van Marle et al. (2005), Arthur (2007), van Marle et al. (2007), Perez-Rendon et al. (2009) show that the MSS and RSGS/LBVS dominate the mass loss with typical values of mass loss before WR stage (21−26) M⊙ for Mini ~ (30−60) M⊙. Meanwhile, the WR stage dominates in the kinetic energy of the wind, injected into the wind bubble.
For the γ2 Velorum binary one can calculate the total stellar wind mass loss ΔM ≃ 29 M⊙ as the difference between the initial and contemporary masses of stars in binary systems, which are MWR,ini ≃ 35 M⊙, MO,ini ≃ 31.5 M⊙ and MWR ≃ 9.0 M⊙, MO ≃ 28.5 M⊙ respectively, according to Eldridge (2009). This includes the mass-loss rate from the slow red supergiant stage wind ΔMRSGS ≃ 19 M⊙ (Freyer et al. 2006). Therefore, the remaining 10M⊙ of the hot intercloud gas in the γ2 Velorum SWB correspond to fast MSS and WR winds. Clumps, created by the interaction of the fast WR wind and the slow RSG wind, together with ISM clouds (Gum nebula interior, see below), survived the passage through the expanded IVS, are the main sources of X-ray emitting plasma in SW part of Vela SNR. As follows from Table 1, the lower limit on the total mass of gas in clumps and clouds inside the SWB is (38)Let us estimate the energy budget of IVS. Assuming that the IVS is a boundary of the SWB of γ2 Velorum we can estimate its kinetic energy, assuming an expansion velocity Vexp ≃ 13km s-1 (Testori et al. 2006) (39)
Thermal energy inside of the SWB of γ2 Velorum is limited by conditions of strong Vela SNR shock and the dominance of the mean density of evaporated clouds in postshock plasma.
Strong shock condition means that the ratio of pressures inside the remnant and bubble and equivalently the ratio of energy densities should be high or (40)The dominance of evaporated cloud material in the postshock region means that and we take hereafter (i.e. C = 10) and as reasonable parameters of the hot (intercloud) gas inside the SWB enclosed in the IVS. The total mass and thermal energy of this hot gas inside γ2 Velorum is and . Meanwhile, the thermal energy of clumps and clouds with a reasonable value of temperature Tcl ≤ 104K is only .
To summarize, in our hydrodynamical model of IVS the total (thermal and kinetic) energy of the γ2 Velorum SWB is (41)with evident dominance of kinetic energy of the massive (~105 M⊙) IVS in the total sum. Inside the SWB we predict about 10 M⊙ of hot intercloud gas and about 260 M⊙ of immersed warm/cool clumps/clouds (in case of their complete evaporation inside Vela SNR). It is worth noting that temperature and density of plasma inside IVS are not well restricted by observations, and here we use approximate values, which are consistent with the adopted limits.
We can compare these estimates with the results of numerical simulations of stellar wind bubbles around WR stars. Numerical simulations of the evolution of a star with an initial mass Mini = 35 M⊙, as proposed by Eldridge (2009) for WR11, were made by Freyer et al. (2006) for an environment with a density of n0 = 20cm-3 and a temperature of T0 = 200K. They show that at the end of the calculations (before the SN explosion) the hot gas bubble has a mean radius of 34pc and shell-like HII and HI regions of the swept up ambient gas extend out to a distance of 43−44pc, the total mass is 1.5 × 105 M⊙, the kinetic energy is 4.9 × 1049 erg, the thermal energy of the hot gas is 1.1 × 1050 erg, of the warm gas 4.3 × 1049 erg, i.e., the radius, total mass and total energy (kinetic and thermal, 2 × 1050 erg) are surprisingly close to our estimate for γ2 Velorum. Nevertheless, two important differences should be clarified for γ2 Velorum namely, the kinetic energy dominance and the low density of the ISM.
The three estimates for the “typical” ISM density derived above provide widely different values. The estimate based on the dynamics of expansion of Vela SNR suggests a low value for the ISM density nISM ≤ 0.01 cm-3, while the estimate based on the total mass of the SWB around γ2 Velorum (nISM ~ 16 cm-3) and the estimate based on the dynamics of expansion of the SWB (nISM ~ 20 cm-3) suggest a much higher density. This points to the fact that the distribution of the ISM in the direction of the Vela region, in the distance range ~300−400 pc, is highly inhomogeneous. This is, in principle, not surprising, because the region is known to contain several stellar formations with different properties.
First, γ2 Velorum belongs to the γVel association, which is a subcluster of the OB-association Vela OB2 (Jeffries et al. 2008). The density of the ISM in the OB association soon after γ2 Velorum’s formation is expected to have been much higher because of the presence of a parent molecular cloud. The initial expansion of the γ2 Velorum bubble into a dense molecular cloud (the progenitor of Vela OB2) can explain the large mass of the swept up ISM in a shell around the SWB. In this scenario, the stellar wind and the radiation of γ2 Velorum destroyed the parent molecular cloud and swept up its gas. For some time the stellar wind of γ2 Velorum was practically trapped inside a ~105 M⊙ cloud, and only when the energy accumulated in the SWB and HII region exceeded the gravitationally bound energy of the cloud (~1050 erg s for Rcl = 10pc cloud), the dense shell of swept-up cloud material began to be accelerated by the thermal pressure of the SWB and the HII region gas without considerable additional mass loading and counter pressure of hot rarefied gas of the ISM. At this evolutionary stage, the thermal energy of the system converts into the kinetic energy of the shell, resulting in the atypical dominance of the kinetic energy of the IVS in the total energy balance of the γ2 Velorum SWB/HII region.
Locations of the Vela SNR (Vela pulsar is shown as a cross), γ2 Velorum (shown as a circle), IRAS Vela Shell (IVS) bubble and Gum nebula (center is shown as a square) in Galactic coordinate system.
Next, both the Vela SNR and the γ2 Velorum SWB could be interacting with a still larger scale SNR, known as the Gum nebula. This nebula is a very large region of the ionized gas about 36° in diameter, centered approximately at (l,b) = (258°, −2°) (Gum 1952), shown in Fig. 2. Brandt et al. (1971) explained the Gum nebula as a fossil Stromgren sphere of the Vela X supernova and estimated the distance to nebula to be (42)Reynolds (1976) suggested that the Gum nebula is an 1 Myr old SNR, which is now heated and ionized by the two very hot stars zeta Puppis and gamma Velorum within it. Later Woermann et al. (2001) showed that the Gum nebula can possibly be a SNR of the zeta Puppis companion. They showed that runaway O-star zeta Puppis was within <0.5 deg of the expansion center of the Gum nebula about 1.5 Myr ago, which is evidence of the relation between the SN explosion of the binary companion and the Gum nebula expansion. Assuming a distance of 400 pc the radius of Gum nebula is about (43)This means that both the Vela SNR and γ2 Velorum are situated inside the cavity formed by the expansion of the SNR associated to the Gum nebula.
The Gum nebula as a very old (0.9−2.0 Myr) SNR should be at the late radiative stage of evolution, which can be modeled with the solution of Cioffi et al. (1988) for the shock dynamics. Assuming the the current radius 124 pc, age ~1.5 Myr and the explosion energy of 1051E51 erg, we can find out the density of the ISM in which the Gum nebula expands: (44)where ζ is the metallicity factor, equal to 1 for solar abundances.
The estimate of the average density of the ISM around Gum nebula shows that actually the nebula expands into the medium with a much lower density than the one suggested by the estimate of the total mass of the SWB around γ2 Velorum. This supports the hypothesis that the enhancement of the density of the ISM was locally present around the γ2 Velorum system at moment of its birth, most probably because of an OB association.
We developed a model for the interaction of the Vela SNR and the γ2 Velorum SWB, which explains the observed NE/SW asymmetry of the Vela SNR.
Adopting a model of the expansion of the Vela SNR into a “cloudy” ISM, we showed that the volume-averaged density of the shock-evaporated clouds in the NE part of the SNR has to be about four times higher than in the SW part. We noticed that a plausible explanation for the observed density contrast is that the Vela SNR exploded at the boundary of the SWB around a nearby Wolf-Rayet star in the γ2 Velorum system, which is situated at approximately the same distance as the Vela SNR.
Within our model of the interaction of the Vela SNR and the γ2 Velorum SWB, the difference of the spectral characteristics of the X-ray emission from the NE and SW parts of the remnant can be used for an estimate of the parameters of the γ2 Velorum bubble. We showed that the measurement of the change of the column density of the neutral hydrogen gives an estimate of the total mass of the SWB, ~1 × 105 M⊙.
On the basis of modeling the dynamics of expansion of the bubble around γ2 Velorum, we confirmed the initial mass of the Wolf-Rayet star in the γ2 Velorum system suggested by Eldridge (2009) to be ≃35 M⊙. This estimate is lower than the previous estimates used for the derivations of predictions of the flux of a γ-ray spectral line at 1.8 MeV, expected from the decays of 26Al in this source. Taking into account the revised estimate of the initial mass of the Wolf-Rayet star, the 26Al line flux flux from γ2 Velorum is expected to be much below the COMPTEL limit.
We would like to thank the referee, John Dickel, for many useful comments and suggestions, which appreciably inproved the paper. I.S. acknowledges support from Erasmus Mundus, “External Cooperation Window”.
- Arthur, S. J. 2007, in Diffuse Matter from Star Forming Regions to Active Galaxies, ed. T.W. Hartquist, J. M. Pittard, & S. A. E. G. Falle (Dordrecht: Springer) 183 [Google Scholar]
- Aschenbach, B., Egger, R., & Trümper, J. 1995, Nature, 373, 587 [NASA ADS] [CrossRef] [Google Scholar]
- Blondin, J. M., Chevalier, R. A., & Frierson, D. M. 2001, ApJ, 563, 806 [NASA ADS] [CrossRef] [Google Scholar]
- Bock D. C.-J., Turtle A. J., & Green A. J. 1998, AJ, 116, 1886 [NASA ADS] [CrossRef] [Google Scholar]
- Brandt, J. C., Stecher, T. P., Crawford, D. L., & Maran, S. P. 1971, ApJ, 163, L99 [NASA ADS] [CrossRef] [Google Scholar]
- Caraveo, P. A., De Luca, A., Mignani, R. P., & Bignami, G. F. 2001, ApJ, 561, 930 [NASA ADS] [CrossRef] [Google Scholar]
- Cha, A. N., & Sembach, K. R. 2000, ApJS, 126, 399 [NASA ADS] [CrossRef] [Google Scholar]
- Cha, A. N., Sembach, K. R., & Danks, A. C. 1999, ApJ, 515, 25 [Google Scholar]
- Cioffi, D. F., McKee, C. F., & Bertschinger, E. 1988, ApJ, 334, 252 [NASA ADS] [CrossRef] [Google Scholar]
- De Marco, O., Schmutz, W., Crowther, P. A., et al. 2000, A&A, 358, 187 [NASA ADS] [Google Scholar]
- Diehl, R., Bennett, K., Bloemen, H., et al. 1995, A&A, 298, L25 [NASA ADS] [Google Scholar]
- Dodson, R., Legge, D., Reynolds, J. E., & McCulloch, P. M. 2003, ApJ, 596, 1137 [NASA ADS] [CrossRef] [Google Scholar]
- Dubner, G. M., Green, A. J., Goss, W. M., Bock, D. C.-J., & Giacani, E. 1998, AJ, 116, 813 [NASA ADS] [CrossRef] [Google Scholar]
- Eldridge, J. J. 2009, MNRAS, 400, L20 [Google Scholar]
- Eldridge, J. J., & Tout, C. A. 2004, MNRAS, 353, 87 [NASA ADS] [CrossRef] [Google Scholar]
- Freyer, T., Hensler, G., & Yorke, H.W. 2003, ApJ, 594, 888 [NASA ADS] [CrossRef] [Google Scholar]
- Freyer, T., Hensler, G., & Yorke, H.W. 2006, ApJ, 638, 262 [NASA ADS] [CrossRef] [Google Scholar]
- Gum, C. S. 1952, Observatory, 72, 151 [Google Scholar]
- Gvaramadze, V. 1999, A&A, 352, 712 [NASA ADS] [Google Scholar]
- Ignace, R., Bessey, R., & Price, C. S. 2009, MNRAS, 395, 962 [NASA ADS] [CrossRef] [Google Scholar]
- Jeffries, R. D., Naylor, T., Walter, F. M., Pozzo, M. P., & Devey, C. R. 2009, MNRAS, 393, 538 [NASA ADS] [CrossRef] [Google Scholar]
- Kasen, D., & Woosley, S. E. 2009, ApJ, 703, 2205 [NASA ADS] [CrossRef] [Google Scholar]
- LaMassa, S. M., Slane P. O., & de Jager, O. C. 2008, ApJ, 689, L121 [NASA ADS] [CrossRef] [Google Scholar]
- Limongi, M., & Chieffi, A. 2006, ApJ, 647, 483 [NASA ADS] [CrossRef] [Google Scholar]
- Lu, F. J., & Aschenbach, B. 2000, A&A, 362, 1083 [NASA ADS] [Google Scholar]
- Miceli, M., Bocchino, F., Maggio, A., & Reale, F. 2005, A&A, 442, 513 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Miceli M., Bocchino F., & Reale F. 2008, ApJ, 676, 1064 [NASA ADS] [CrossRef] [Google Scholar]
- Miceli, M., Reale, F., Orlando, S., & Bocchino, F. 2006, A&A, 458, 213 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Millour, F., Petrov, R. G., Chesneau, O., et al. 2007, A&A, 464, 107 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- North, J. R., Tuthill, P. G., Tango, W. J., & Davis, J. 2007, MNRAS, 377, 415 [NASA ADS] [CrossRef] [Google Scholar]
- Oberlack, U., Wessolowski, U., Diehl, R., et al. 2000, A&A, 353, 715 [NASA ADS] [Google Scholar]
- Palacios, A., Meynet, G., Vuissouz, C., et al. 2005, A&A, 429, 613 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Perez-Rendon, B., Garcia-Segura, G., & Langer, N. 2009, A&A, 506, 1249 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Rajagopal, J., & Srinivasan, G. 1998, JApA, 19, 79 [NASA ADS] [CrossRef] [Google Scholar]
- Reichley, P. E., Downs, G. S., & Morris, G. A. 1970, ApJ, 159, 35 [Google Scholar]
- Reynolds R. J. 1976, ApJ, 206, 679 [NASA ADS] [CrossRef] [Google Scholar]
- Sahu, M.S. 1992, Ph.D. Thesis, Rijksuniversiteit Groningen [Google Scholar]
- Schaerer D.,Schmutz, W., & Grenon, M. 1997, ApJ, 484, 153 [Google Scholar]
- Schanne, S., Sizun, P., Maurin, D., et al. 2007, Proceedings of the VI INTEGRAL Workshop. July 2–8, 2006, Moscow, Russian Federation, ed. S. Grebenev, R. Sunyaev, C. Winkler, ESA SP-622 (Noordwijk: ESA Publication Division), 113 [Google Scholar]
- Sedov, L. I. 1959, Similarity and Dimensional Methods in Mechanics (New York: Academic Press) [Google Scholar]
- Stothers, R. 1980, PASP, 92, 145 [NASA ADS] [CrossRef] [Google Scholar]
- Testori, J.,C., Arnal, E. M., Morras, R., et al. 2006, A&A, 458, 163 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Vanbeveren, D. 1991, A&A, 252, 159 [NASA ADS] [Google Scholar]
- Vanbeveren, D., De Donder, E., van Bever, J., van Rensbergen, W., & De Loore, C. 1998, New Astron. 3, 443 [NASA ADS] [CrossRef] [Google Scholar]
- van Leeuwen, F. 2007, Hipparcos, the new reduction of the raw data (Dordrecht: Springer) [Google Scholar]
- van Marle, A., Langer, N., & Garcia-Segura, G. 2005, A&A, 444, 837 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- van Marle, A. J., Langer, N., & Garcia-Segura G. 2007, A&A, 469, 941 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Weaver, R., McCray, R., Castor, J., Shapiro, P., & Moore, R. 1977, ApJ, 218, 377 [NASA ADS] [CrossRef] [Google Scholar]
- White, R. L., & Long, K. S. 1991, ApJ, 373, 543 [NASA ADS] [CrossRef] [Google Scholar]
- Woermann, B., Gaylard, M. J., & Otrupcek, R. 2001, MNRAS, 325, 1213 [NASA ADS] [CrossRef] [Google Scholar]
- Yamaguchi, H., & Katsuda, S. 2009, ApJ, 696, 1548 [NASA ADS] [CrossRef] [Google Scholar]
ROSAT All-Sky Survey image (0.1–2.4 KeV) of the Vela SNR (Aschenbach et al. 1995). A–F are extended features outside the boundary of the remnant (“bullets”). Light blue to white contrast represents a contrast in surface brightness of a factor of 500 (Aschenbach et al. 1995). Blue curves show the NE and SW hemispheres of the Vela SNR. The yellow curve shows the contour of the SWB of γ2 Velorum.
|In the text|
Current usage metrics show cumulative count of Article Views (full-text article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.
Data correspond to usage on the plateform after 2015. The current usage metrics is available 48-96 hours after online publication and is updated daily on week days.
Initial download of the metrics may take a while.