{'fraction_non_alphanumeric': 0.03802872694728503, 'fraction_numerical': 0.03848354582666999, 'mean_word_length': 4.13755935780508, 'pattern_counts': {'":': 0, '<': 8, '<?xml version=': 0, '>': 1, 'https://': 0, 'lorem ipsum': 0, 'www.': 0, 'xml': 0}, 'pii_count': 0, 'substrings_counts': 0, 'word_list_counts': {'cursed_substrings.json': 2, 'profanity_word_list.json': 0, 'sexual_word_list.json': 0, 'zh_pornsignals.json': 0}}
XMM-Newton observations of the supernova remnant RX J1713.7-3946 and its central source

We present new results coming from the observations of the supernova remnant (SNR) RX J1713.7--3946 (also G347.3--0.5) performed in five distinct pointings with the EPIC instrument on board the satellite XMM-Newton. RX J1713.7--3946 is a shell-type SNR dominated by synchrotron radiation in the X-rays. Its emission (emission measure and photon index) as well as the absorption along the line-of-sight has been characterized over the entire SNR. The X-ray mapping of the absorbing column density has revealed strong well constrained variations (0.4 10^22 cm^-2<NH<1.1 10^22 cm^-2) and, particularly, a strong absorption in the southwest. Moreover, there are several clues indicating that the shock front of RX J1713.7--3946 is impacting the clouds responsible for the absorption as revealed for instance by the positive correlation between X-ray absorption and X-ray brightness along the western rims. The CO and HI observations show that the inferred cumulative absorbing column densities are in excellent agreement with the X-ray findings in different places of the remnant on condition that the SNR lies at a distance of 1.3 +/- 0.4 kpc, probably in the Sagittarius galactic arm, instead of the commonly accepted value of 6 kpc. An excess in the CO emission is found in the southwest suggesting that the absorption is due to molecular clouds. A search for OH masers in the southwestern region has been unsuccessful, possibly due to the low density of the clouds. The X-ray mapping of the photon index has also revealed strong variations (1.8<Gamma<2.6). [See the article for the entire abstract]


Introduction
Supernova remnants (SNRs) are well known to be sources of radio synchrotron emission from GeV relativistic electrons. During the last 10 years, X-ray synchrotron emission from TeV relativistic electrons has been discovered in several shell-type SNRs among which are SN 1006 (Koyama et al. 1995), RX J1713.7-3946 (Koyama et al. 1997, Slane et al. 1999) and RX J0852.0-4622 (Slane et cm −3 , where D 6 is the distance in units of 6 kpc. Recently, this thermal component may have been detected in the interior of the SNR over large fields-of-view using ROSAT, ASCA and RXTE satellites (Pannuti et al. 2003). However, some possible local variations of absorbing column density and photon index can compromise this interpretation.
An intriguing X-ray point-source (1WGA J1713.4-3949) with no optical and radio counterparts was also detected at the center of the remnant, being possibly a neutron star or an extra-galactic object (SL99). Recently, a combined Chandra and XMM-Newton spectral analysis of this source has shown that its X-ray properties (twocomponent spectrum, luminosity, absence of pulsations) resemble those of other compact central objects in SNRs and suggests that 1WGA J1713.4-3949 is the neutron star associated to RX J1713.7-3946 (Lazendic et al. 2003).
Neither the distance nor the age of RX J1713.7-3946 is well established. From the X-ray measurements of the column density toward this source, Koyama et al. (1997) derived a distance of 1 kpc, while SL99 proposed a larger distance of 6 kpc based on its probable association with three dense and massive molecular clouds and the Hii region G347.6+0.2 located northwest to it. However that may be, if the Sedov evolution is assumed with 6 kpc, an age of a few 10 4 years results, whereas an age of one order of magnitude below is compatible with 1 kpc. Hence, this nearer distance could be in agreement with the hypothesis proposed by Wang et al. (1997), based on historical records, that RX J1713.7-3946 is the remnant of the supernova (SN) that exploded in AD 393. Following studies adopted the largest distance of 6 kpc, until new XMM-Newton observations (Cassam-Chenaï et al. 2004b) that re-opened the whole distance question together with new high-resolution CO mm-wave observations (Fukui et al. 2003). These new results suggest possible indications of interaction between the SNR shock front and molecular gas located at 1 kpc in the northwest and southwest sides of the SNR.
In the radio, the emission arises from faint filaments aligned with the X-ray shell of RX J1713.7-3946, except for the brightest filament located in the northwestern rim which is perpendicular to the shock and coincident with part of the Hii region G347.6+0.2 (SL99, Ellison et al. 2001, Lazendic et al. 2004. GeV γ-ray emission was detected by the EGRET instrument to the northeast of the SNR (Hartman et al. 1999). This emission was interpreted as the decay of neutral pions attributed to the interaction of CR nuclei (accelerated at the shock in RX J1713.7-3946) with the massive and dense cloud assumed to be interacting with the remnant (Butt et al. 2001). By arguing on the photon spectral index of the GeV source (α = 2.3 ± 0.2 which is expected from the hadronic interactions of CR source) and on the non-detection in radio of the massive cloud Butt et al. (2001) ruled out the possibility that electrons are responsible for the GeV luminosity.
At TeV γ-rays, the CANGAROO imaging Cerenkov telescope detected emission in the northwest of the SNR, which was interpreted as Inverse Compton (IC) emission from accelerated electrons (Muraishi et al. 2000). The low matter densities in the ambient medium is unfavorable to the interpretation in terms of π 0 decay process. A second CANGAROO observation of RX J1713.7-3946 leads to a TeV γ-ray spectrum (between 400 GeV and 8 TeV) whose photon spectral index (α = 2.8 ± 0.2) was shown to be consistent with the π 0 decay process and no other mechanism (Enomoto et al. 2002). However, such a spectrum would exceed the EGRET observed emission by a factor three (Reimer & Pohl 2002, Butt et al. 2002. Several self-consistent models have been constructed using the broad-band spectrum from radio to γ-ray wavelengths, essentially for the northwest region of RX J1713.7-3946. Using a non-linear diffusive shock acceleration model, Ellison et al. (2001) were able to reproduce the broadband spectrum with synchrotron emission from shock-accelerated electrons in X-rays and IC emission in γ-rays. In this model, more than 25% of the shock kinetic energy is going into relativistic ions. Other attempts to reproduce the broadband spectrum were undertaken (Uchiyama et al. 2003;Pannuti et al. 2003;Lazendic et al. 2004).
In this paper, we provide for the first time a detailed description of the X-ray emission of RX J1713.7-3946 based on XMM-Newton observatory data. Thanks to the large field-of-view provided by the EPIC imaging spectrometers, we can cover the bulk of the SNR extent in a few pointings. The high sensitivity of XMM-Newton allows us to carry out a spectral analysis at medium-scale of the emission structures and then to produce for the first time a mapping of the spectral parameters of RX J1713.7-3946.
In Section 2, we present our observations and the different data reduction methods that we have used. In Section 3, we present the results obtained from the XMM-Newton data analysis. Section 4 contains a discussion on the potential association between the central point-like source 1WGA J1713.4-3949 and SNR RX J1713.7-3946 (Sect. 4.1), on the interaction between the remnant and clouds as well as on the SNR distance (Sect. 4.2), on the SNR age and energetics (Sect. 4.4) and finally on the overall picture that comes out from the results (Sect. 4.5). In Section 5, we summarize our results.

Observation
SNR RX J1713.7-3946 was observed with the three EPIC instruments (namely the MOS1, MOS2 and pn cameras) on board the XMM-Newton satellite in the course of the AO-1 program. The observations were performed in five distinct pointings, each of around 10 ks duration. Table 1 indicates the date of observation and the OBS ID for each pointing. The five pointings were chosen based on the ROSAT PSPC image (SL99) to allow a maximum coverage of the remnant and the possibility to extract a local X-ray background. The pointings correspond schematically to the center, northeast, northwest, southwest and southeast (hereafter CE, NE, NW, SW, SE, respectively) parts of RX J1713.7-3946. During the observation, the MOS and pn cameras were operated in Full Frame and Extended Full Frame providing a temporal resolution of 2.6 s and 200 ms, respectively. For all the pointings, the medium filter was used to avoid contamination by visible sources contained in the field of view and to limit the number of low energy photons (and hence pile-up) coming from the bright X-ray point source 1WGA J1713.4-3949.

Flare rejection
The Science Analysis System (SAS version 5.3) was used for data reduction.
First, we have generated calibrated events files with the SAS tasks emchain and epchain. To create light curves, spectra and images, we have selected the single, double, triple and quadruple events (pattern≤12) for the MOS cameras, and the singles and doubles (pattern≤4) for the pn camera (with FLAG=0 for all cameras). The data quality of the CE, SW, SE pointings is very good. However, it is not the case for the NE and NW pointings which are contamined by soft protons flares. The method used to clean the data from the flares consists in applying a threshold on the light curves at high energy, i.e.
where the EPIC cameras are almost insensitive to X-ray photons. We took a threshold of 18 counts per 100 seconds time interval in the 10-12 keV band for the MOS cameras (pattern≤12) and a threshold of 22 counts per 100 seconds time intervals in the 12-14 keV band for the pn camera (pattern=0). Table 1 indicates the remaining exposure time after flares cleaning for each EPIC camera.

Data reduction methods
In this section, we describe the different methods that have been used to create the images, the hardness ratio maps, the spectra and the mapping of the spectral parameters.

Images
All the mosaiced images of the X-ray emission presented here are background subtracted, exposure and vignetting corrected (Sect. 3.1). The background used for the subtraction is only instrumental and is estimated by using a reference observation (Read & Ponman 2003) which has been obtained with the same medium filter as in our observations. No correction for the local astrophysical background is made because we have no information on it in the center of the remnant and further more it strongly varies around the remnant which renders difficult any extrapolation towards the interior. To correct for the exposure and the vignetting, we divide the mosaiced background-subtracted count image by its mosaiced exposure map which includes the vignetting correction using the SAS task eexpmap. Each time, the mosaiced count image is adaptively smoothed with the SAS task asmooth and the resulting template is applied to smooth both the mosaiced background-subtracted count image and its associated exposure map before the division.
We have not used the pn data for the intensity maps as it introduces some problems of homogeneity in the correction for vignetting and exposure due to the lack of pn data in the NE pointing (see Table 1).

Hardness ratio maps
To create hardness ratio maps, we divide two background subtracted images made in different energy bands that are previously adaptively smoothed. For the smoothing, the template of the image with the lowest statistics is always applied to that with higher statistics. When we compute the ratio of the two maps, the statistical fluctuations can be amplified. To limit this effect, the MOS and pn data are used, but in return, the NE is missing in the hardness ratio maps (Sect. 3.3).

Spectra
All spectral analyses use the weight method described in Arnaud et al. (2001) to correct for the vignetting. As for the background subtraction, we always use the double subtraction method described in Arnaud et al. (2002) except for the mapping of the spectral parameters (see below, Sect. 2.3.4). This method allows to take into account the difference, at low energy, between the local astrophysical background of our observations and the one found in the reference observation. This is done for local spectral analyses which require a precise knowledge of the spectrum at low energy (Sect. 3.3). In our observations, there are two regions where the local astrophysical background can be extracted with enough statistics: the SW and the SE. Each time, we choose the nearest available local astrophysical background for each selected region.
The spectra are grouped so that each bin has a signalto-noise ratio greater than 5σ. To fit the spectra, we use XSPEC (version 11.2, Arnaud 1996). The chosen on-axis response files are m * medv9q20t5r6 all 15.rsp for MOS and epn ef20 sdY9 medium.rsp for pn.

Mapping of the spectral parameters
The X-ray spectrum of SNR RX J1713.7-3946 can be well fitted by a simple absorbed power-law model the spectral parameters of which are the absorbing column N H and the photon index Γ. Here, we describe how we map the spatial variations of these two parameters (Fig. 4).
First, we build a spatial grid which takes into account the statistics so that we have approximately the same number of counts in each grid pixel. The spatial grid is constructed from a mosaiced count image (MOS+pn) between 0.8-10 keV where all the point-like sources have been removed. The mean number of total counts in a pixel is fixed to 9000 counts from this image which corresponds approximately to 2500 and 4000 counts in most of the boxes of the MOS and pn cameras, respectively.
For each pixel of the grid, we extract the X-ray spectra (MOS, pn) and fit them with a simple absorbed power law. A map of these parameters (N H , Γ) is constructed from the best-fit values in each grid pixel.
In this spectral analysis, we use a single backgroundsubtracted method (as described for the images, Sect. 2.3.1) since it is not possible to have the true local astrophysical background in each grid pixel. Because of that, the spectral bins only above 0.8 keV are selected, i.e. where the reference background provided by Read & Ponman (2003) is almost one order of magnitude below the spectrum of our observations. We have checked that the absence of spectral bins below 0.8 keV did not change significantly the value of the spectral parameters in the X-ray bright regions.
There is no significant difference between the single and the double background-subtraction methods for the brightest regions. Nevertheless, in the X-ray faint regions, the effect of the double substraction is to increase N H and Γ by a factor ∼ 10 − 15% which is not enough to alter our results and interpretation.

X-ray results
Unlike the X-ray spectra of typical shell-like SNRs which exhibit emission lines, the spectrum of RX J1713.7-3946 is found to be featureless as shown for example in Fig.  1. This suggests a non-thermal origin and has been interpreted as synchrotron X-ray emission from TeV shock accelerated electrons (SL99). RX J1713.7-3946 is a second example, after SN 1006, of a shell-like SNR for which non-thermal emission dominates the X-ray flux, but more extreme in the sense that no thermal emission has been unambiguously detected yet. Unlike SN 1006, a symmetrical limb-brightened type Ia remnant located in a high-latitude environment, RX J1713.7-3946 has a quite distorted morphology which may be related to its complex surrounding environment in the Galactic plane and is likely to be the remnant of a type II explosion as suggested by the presence of a central point source detected at its center.
In section 3.1, we describe in detail the unusual Xray morphology of RX J1713.7-3946, and compare it to the radio continuum morphology. In section 3.2, we trace and accurately diagnose the spectral variations of the synchrotron emission. In section 3.3, we search for the thermal emission and constraints on its properties. Finally, in sec- Both images were adaptively smoothed to a signal-to-noise ratio of 10. The scaling is square root with the higher cut fixed to 30% of the SNR maximum X-ray brightness. The color bar numbers are in units of counts/s/pixel with a pixel size of 4 ′′ . The bad resolution in the NE is due to a very short exposure time (see Table 1). Moreover, it is possible that the smoothing has removed a few small-scale features in this NE pointing. tion 3.4, we intent to give better spectral constraints on the nature of the point-like source 1WGA J1713.4-3949 to provide arguments in favor or against a potential association with SNR RX J1713.7-3946.

X-ray morphology
The good spatial resolution and sensitivity of XMM-Newton combined to its large field-of-view allow us for the first time to both have an overall view of RX J1713.7-3946 and resolve some X-ray emission features not seen before in the ROSAT and ASCA images (SL99, Uchiyama et al. 2003). Figure 2 shows the 5 pointings mosaiced image of RX J1713.7-3946 in the 0.8-2 keV (left panel) and 2-10 keV (right panel) energy bands. However, the 5 pointings are not enough for a complete coverage of the SNR. The general shape of this distorted remnant resembles roughly an ellipse stretched out along the NW-SE axis over ∼ 70 ′ corresponding to ∼ 20 pc at 1 kpc or ∼ 120 pc at 6 kpc. An overall view of RX J1713.7-3946 obtained with ASCA (Uchiyama et al. 2003) shows that the boundaries of the SNR are irregular and suggests that the remnant extends farther to the north.
Except for the point-like sources, the comparison between the low energy and high energy images (see Fig. 2) does not show significant changes in the X-ray morphology at first sight, at least for the brightest X-ray features. From these images, the X-ray emission can be schemat-  ically classified in three distinct levels. The lower intensity one is due to the diffuse X-ray emission of the local astrophysical background. In both images, it can be found (in blue-black) on the outside of the NE, NW, SW and SE pointings. It stands out that the local astrophysical background intensity varies around the remnant. This fact is expected since RX J1713.7-3946 lies just below the Galactic Ridge in an environment where the properties of the interstellar gas rapidly change with Galactic latitude (which is the direction where the SNR is the more stretched out). The medium intensity level corresponds to the diffuse and faint X-ray emission arising from the SNR which clearly dominates in the 0.8-2 keV energy band (in green and also in blue). The higher intensity level is made of a more structured emission which appears in both images but more clearly at high energy in the 2-10 keV band (in yellow and red). Its surface brightness is about twice higher than the diffuse emission associated with the SNR. Both the NW and SW limbs are very bright in X-rays. The X-ray emission in the western half is more structured than anywhere else in the SNR. At large scale, a doublering morphology stands out. The internal ring is more or less circular whereas the external one looks broken into a few arcs. At smaller scale, the Chandra image of the NW region shows that the brightest X-ray features correspond to thin filaments of 20 ′′ apparent width or hot spots embedded in a diffuse plateau emission (Uchiyama et al. 2003, Lazendic et al. 2004). This is confirmed by XMM-Newton in the NW and also in the SW.
In spite of its bad quality, the NE observation shows that the inner circular thin shell seen in the west seems to continue there (see Fig. 2, right panel). Nevertheless, the diffuse X-ray emission from the NE pointing seems to exhibit a concave border between the remnant and the ambient medium (green emission in Fig. 2, left panel). Such a shape is not naturally expected from an explosion into a homogeneous low density interstellar medium. That shape must be caused by some obstacles either interacting with the SNR or obstructing the X-ray emission arising from it.
The X-ray emission from the SE is weak and diffuse but we distinguish even so a possible border between the SNR and the ambient medium.
Finally, the central part of the SNR shows diffuse X-ray emission rather weaker than in the NW and SW but similar to what is observed in the SE. In addition, two point sources are observed. The first one is 1 WGA J1714.4-3945 located at the position α J2000 =17h14m30s, δ J2000 =-39 • 46 ′ 00 ′′ which disappears above 2 keV. The second one located in the middle of the CE pointing is 1WGA J1713.4-3949, a potential neutron star, which will be studied in detail in Section 3.4. Figure 3 shows a high-resolution radio image obtained at 1.4 GHz with ATCA (from Lazendic et al. 2004) overlaid with the X-ray contours. The radio morphology of RX J1713.7-3946 is quite complex, with the brightest emission corresponding to the Hii region G347.6+0.2, to the NW corner. A well defined ring of emission, about 15 ′ in size, is visible in the interior of the SNR. This radio feature matches the inner X-ray bright ring, particularly towards the northern, northwestern and western sides. Also, a very good correspondence is observed on the X-ray-bright outer rim to the SW. The observed similar morphology in the radio and in the X-rays strongly suggests that the synchrotron emission detected at the two spectral domains must originate in the same electrons population accelerated at the shocks.

Spatial and spectral characterization of the synchrotron emission
In this section, we present the maps of absorbing column density and photon index (method described in Sect. 2.3.4). Figure 4 (left panel) shows the mapping of absorbing column density in the line-of-sight at medium scale. The mean relative error on the absorbing column in each pixel grid is ∼ 9% with a maximum value of 16%. The variations of absorbing column density are strong with N H varying from ∼ 0.4 × 10 22 cm −2 to ∼ 1.1 × 10 22 cm −2 . In the SE and CE pointings, N H is low with a value ∼ 0.4−0.5×10 22 cm −2 whereas it is larger in the NW (N H ∼ 0.6−0.7×10 22 cm −2 ) and even more in the SW (N H ∼ 0.8 − 1.1 × 10 22 cm −2 ). Figure 5 (top panel) illustrates the variations of absorbing column from two spectra extracted in the SW region.
The comparison between the mapping of absorbing column density and the X-ray brightness contours (Fig. 4, left panel) suggests that there is a positive correlation them. Figure 6 shows the X-ray brightness at 1 keV versus N H in each grid pixel. We see that the regions of large (small) N H correspond to those of high (weak) X-ray brightness confirming the significance of the observed correlation. The fact that high absorbing column coincides with brighter X-ray emission, as observed to the SW (and NW), is unexpected. Indeed, a dense cloud located in front of an Xray source and not related to it must diminish the photon source flux in our direction (at low energy in particular). That we observe a positive correlation between the absorbing column density and the X-ray brightness means  photons/cm 2 /s/keV/arcmin 2 ). Each point corresponds to the value found for a particular pixel of the mapping grid (Sect. 2.3.4) and the error bars are given at a 90% confidence level. The points of the SW region are selected so that they are above 2.0 × 10 −4 photons/cm 2 /s/keV/arcmin 2 . Bottom panel : Correlation plots between the photon index Γ and the normalized brightness at 1 keV for the same set of points as above. The error bars are given at a 90% confidence level. that the increased ambient density somehow amplifies the X-ray brightness. This strongly suggests that the remnant is interacting in the brightest regions with part of the lineof-sight absorbing material. Figure 4 (right panel) shows the mapping of the photon spectral index. The mean relative error on the photon index value in each grid pixel is ∼ 4% with a maximum value of 7%. The variations of photon index are important seeing that Γ varies from ∼ 1.8 to ∼ 2.6. The spectrum is steep in the faint central region whereas it is a bit flatter in the NW and even more in the SE. The fact that the spectrum gets steeper as we move inwards away from the shock front is expected in SNRs. This effect is generally explained by the adiabatic expansion of the remnant and the synchrotron radiative losses (Reynolds 1998) as is also observed in SN 1006. In the SW, the spectrum is generally steep, but some regions exhibit strong variations of spectral index as is shown in Fig. 4 (right panel). For illus-  Table 2). The black triangle correspond to the power-law plus thermal equilibrium model (cf. Table 3). The error bars are given at a 90% confidence level. tration, Fig. 5 (bottom panel) shows the spectrum of two regions in the SW, for which the index is very different.
The synchrotron X-ray morphology, though quite complex, can be related to the absorption along the line-ofsight (or equivalently brightness). On the one hand, there is a weakly absorbed SE region which seems to extend far away from the center of the SNR (taken as the pointlike source 1WGA J1713.4-3949), and on the other hand, a strongly absorbed western region which is closer and where RX J1713.7-3946 is suspected to interact with the material that makes the absorption. Such an asymmetry in the morphology and the absorption must be the result of different ambient densities around the remnant. To investigate how particle acceleration is modified in different environment, we have selected two sets of points located at the presumed shock locations: the one in the SW external rim coincident with the radio emission (see Fig. 3) and in the SE. These points are shown in Fig. 7 in the absorbing column density-X-ray brightness plane (top panel) and in the photon index-X-ray brightness plane (bottom panel). It is obvious that two groups come out. The regions where the shock propagates in the SW have a steep spectrum (Γ ∼ 2.3 − 2.7) and those where the shock propagates in the SE have a flat spectrum (Γ ∼ 1.9 − 2.2) 1 . This point is discussed further in Sect.4.3.

Thermal emission
Determining the thermal properties of the hot gas is crucial to constrain the evolutionary state of the remnant. That the thermal emission is not clearly detected in SNR Region 952 (909) Table 2. Best fit parameters for the power-law model using the MOS and pn data all at once. The errors are in the range ∆χ 2 < 2.7 (90% confidence level) on one parameter. The regions shown in Fig. 8. N orm is given in units of photons/cm 2 /s/keV/arcmin 2 at 1 keV. Region 933 (907) Table 3. Best fit parameters for the power-law plus equilibrium model using the MOS and pn data all at once. The errors are in the range ∆χ 2 < 2.7 (90% confidence level) on one parameter. The regions are shown in Fig. 8. N orm is given in units of photons/cm 2 /s/keV/arcmin 2 at 1 keV.
RX J1713.7-3946 sets however strong limits on the temperature of the hot gas and, above all, on its mean density. On the one hand, the temperature determination is important because it gives an idea on the shock velocity (under a few assumptions on the heating of the electrons behind the shock). On the other hand, the estimation of the mean ambient medium density is crucial to provide an estimate of the swept-up mass (assuming that the distance to the SNR is known) which then constrains the phase of evolution of the remnant.
To characterize the synchrotron emission (Sect. 3.2), the spectra of SNR RX J1713.7-3946 were modelled by an absorbed phenomenological power-law which provides, at the first order, an excellent representation of the data. To search for the thermal emission, we add a second component, the thermal model, which will not necessarily improve the spectral fit but will be essential to constrain the properties of the thermal gas. For this task, it is important to look at spectrally homogeneous areas since the spectrum of a region with strong variations of N H and Γ modelled with an uniformly absorbed single power-law could simulate artificially the need for a thermal additional component instead of a second non-thermal component. Eight particular regions (point-like sources removed) with no or few spectral variations were chosen based on the low-energy and high-energy hardness ratio maps, and such that they have enough statistics. They are shown in Fig. 8. Although the hardness ratio maps do not give the value of Γ or N H , they permit to keep a more precise spatial information than in the photon index or absorbing column map which intrinsically depend on the shape and size of the pixels of the grid.
Like in previous studies where an attempt to investigate the thermal emission properties was done (SL99, Pannuti et al. 2003), we will first assume the plasma to be in an equilibrium ionization state. This hypothesis will be discussed further. In the following study, the elemental abundances are fixed to solar.
First, we investigate how our previous results (Sect. 3.2) might be modified by the addition of a thermal component in the model. Table 2 gives the best-fit parameters obtained for an absorbed power-law model in each of the eight regions selected for this analysis, while Table 3 gives the results obtained with the addition of a thermal equilibrium ionization component. It appears that the introduction of the thermal model improves very slightly the fit (∆χ 2 ∼ 10 − 30) and does not change qualitatively the results found for the non-thermal emission as is shown in Fig. 9. The derived parameters (N H , Γ) are similarly changed for all the eight regions: they always get slightly larger as we add the thermal component to the model.
The comparison between the fluxes of the thermal emission obtained in different regions of the SNR gives an idea on the spatial distribution of the shocked medium density. Table 3 gives the value of n e n H dl and the temperature found for the thermal emission (for the 8 regions). The range of temperatures kT e is always below 1 keV. The values of the square density integrated along the line-ofsight are rather homogeneous (∼ 1 × 10 16 cm −5 ). We note however that this quantity seems to increase from the SNR shock front (regions 1 and 8) to the interior (regions 2 and 6) when we look at a temperature around 0.5-0.6 keV.
The postshock electronic number densityn e ≡ n 2 e dl/ dl 1/2 is determined by (assuming that n e = where the angular distance between the center of the SNR (taken as the position of the central point-like source 1WGA J1713.4-3949) and the considered region is θ r while that of the blast-wave is θ s . Table 4 shows that the average postshock electronic density is approximately at the same level over the SNR (n e ∼ 2.1 +0.8 −0.9 × 10 −2 cm −3 for D = 1 kpc), being a bit larger in the interior (n e ∼ 2.4±0.5×10 −2 cm −3 ) than close to the shock (n e ∼ 1.9±0.5×10 −2 cm −3 ). In a Sedov model, the average postshock density is similar to the mean number density in the ambient medium in the inner regions of the SNR. Then, the postshock densities lead to a mean hydrogen number density n 0 of the ambient preshock medium between 0.010 and 0.024 D −1/2 1 cm −3 .
Nevertheless, because of the extremely low density in the ambient medium, the assumption of equilibrium ionization state does not hold any more. Indeed, since the ionization age must be lower than max(n e ) t s = 9.0 × 10 8 D 1 kpc 1/2 ts 10 3 yr s cm −3 where t s is the SNR age, even a remnant of 10 5 years old would be insufficient to reach the ionization equilibrium. Consequently, we have applied a non-equilibrium model instead of the equilibrium model while keeping the power-law model. This has been done by relating the ionization age τ to the emission measure as described in Cassam-Chenaï et al. (2004a, Sect. 6.3). It is found that this model does not improve the fit obtained with the power-law model alone. In the power-law plus non-equilibrium model, the power-law characteristics are the same as those obtained for a single power-law model ( Table 2). The temperature of the thermal gas is found to be around 3 keV but is poorly constrained. Moreover, an upper limit on the ionization age of 4.5 × 10 8 s cm −3 is obtained. Assuming t s = 1600 years (cf. Sect. 4.4), this leads to an upper limit of 2.6 × 10 −2 D −1/2 1 cm −3 for the average postshock electronic density and 0.02D −1/2 1 cm −3 for the mean hydrogen number density of the ambient preshock medium (D 1 is the distance in units of 1 kpc). Whatever the thermal model (equilibrium or nonequilibrium) added to the nonthermal model, we derive similar mean hydrogen number densities n 0 of the ambient preshock medium. In the following, we will adopt the upper limit n 0 ≤ 2 × 10 −2 D −1/2 1 cm −3 .

Central source 1WGA J1713.4-3949
The nature of the point-like source 1WGA J1713.4-3949 located at the center of SNR RX J1713.7-3946 is not clearly established: a neutron star or a background extragalactic source seem to be possible interpretations (SL99). However, the neutron star hypothesis is the most attractive because of the lack of radio emission and optical counterpart, and because the derived upper limit column density from the ASCA X-ray spectra is comparable with the total absorption through the Galaxy in this direction (SL99).
Here, we intent to give better spectral constraints on the nature of this source 1WGA J1713.4-3949 thanks to the high sensitivity of XMM-Newton. The spectral analysis of the source, as well as the knowledge of the absorbing column density map of the remnant (Sect. 3.2), must provide arguments in favor or against a potential association with the SNR RX J1713.7-3946.
We extracted the source counts from a 25 ′′ radius circle and the background was taken in an annulus of radii R min = 60 ′′ and R max = 200 ′′ around the point source center. The background subtracted source count rates obtained are 0.410 ± 0.006, 0.413 ± 0.006 and 1.148 ± 0.013 counts s −1 for MOS1, MOS2 and pn, respectively. We extracted the spectrum of the central source in the 0.3-12.0 keV band. As shown in Fig. 10, it does not extend higher than 6.5 keV. No thermal spectral emission lines are found in the EPIC data, nor in the RGS data around 1 keV.
We first considered one component simple models, the results of which are summarized in Table 5. Fitting the spectrum with a simple absorbed power law yields a large photon index Γ = 4.24 ± 0.06. Such a value is much steeper than what is generally found in other young neutron stars which are well known for their non-thermal emission from relativistic particles accelerated in the pulsar magnetosphere (for instance, Willingale et al. 2001 give Γ = 1.63 ± 0.09 for the Crab pulsar). Moreover, the fit overestimates the spectrum at high energy resulting in a bad χ 2 r = 1.58 (for 475 dof). The column density N H = 1.10 ± 0.03 × 10 22 cm −2 inferred exceeds the values found in the SNR surrounding medium.
If the source is not an active pulsar but a neutron star, thermal radiation from the star surface can be observed. A single blackbody model improves the fit (χ 2 r = 1.31 for 475 dof) but does not reproduce the hard emission. It yields a high temperature of kT = 0.40 ± 0.2 keV, a small radius R = 0.37±0.01D 1 km of emitting area and a luminosity of L = 0.49 ± 0.02 × 10 33 D 2 1 erg s −1 where D 1 is the distance in units of 1 kpc. The value of L is lower than predicted by standard cooling neutron star models for which the typical luminosity is ∼ 10 34 erg s −1 for a neutron star of a few thousand years (e.g. Tsurata 1998). Such a low luminosity could be explained by an accelerated cooling process. However the small size of emitting area along with the high temperature means that the cooling cannot take place from the entire neutron star surface. Note that the temperature, radius and luminosity of the central point source 1WGA J1713.4-3949 are almost identical to what is found for the central point source in SNR RX J0852.0-4622 (also G266.2-1.2 or "Vela Junior") adopting a distance of 1 kpc (see Becker & Aschenbach 2002, Kargaltsev et al. 2002. The column density N H = 0.35 ± 0.02 × 10 22 cm −2 is lower than the values found in the SNR. We note that our range of temperatures is in the error bars given by SL99 whereas both our column density and radius are lower (in SL99, N H = 0.52 +0.18 −0.16 × 10 22 cm −2 and R = 0.50 +0.17 −0.11 D 1 km). The latter point should be related to the weak flux reduction. Indeed, the unabsorbed flux in the 0.5 − 5.0 keV energy band obtained with ASCA (5.3×10 −12 erg cm −2 s −1 , SL99) is larger than the one obtained with XMM-Newton (3.8 × 10 −12 erg cm −2 s −1 ).
A simple thermal bremsstrahlung yields an excellent fit (χ 2 r = 1.05 for 475 dof) with a high temperature kT = 0.94 ± 0.2 keV and an absorbing column density of N H = 0.70 ± 0.02 × 10 22 cm −2 compatible with the value derived at the center of the SNR.
A more precise description of 1WGA J1713.4-3949 can be given by testing more elaborated models with two spec-   Table 6. Best fit parameters for different two components spectral models for 1WGA J1713.4-3949 using the MOS and pn data all at once. The errors are in the range ∆χ 2 < 2.7 (90% confidence level) on one parameter. D 1 is the distance in units of 1 kpc.
tral components. Table 6 summarizes the results obtained.
To allow for different emission mechanisms, we first introduce a soft thermal component arising from the neutron star surface with a blackbody model and a non-thermal component described by a power law due to relativistic particles heating the polar caps of the neutron star. The resulting fit is similar to that obtained with the thermal bremsstrahlung model. The photon index Γ = 4.24 ± 0.06 is steep in comparison to spectral index found in young neutron stars. Both the luminosity and the temperature are a bit lower than what is expected for a cooling neutron star of a few thousand years (Tsurata 1998). The column density N H = 0.91 ± 0.05 × 10 22 cm −2 inferred exceeds the values found in the SNR center.
A two-component blackbody model allows to represent the non-uniformity of the surface temperature distribution of the neutron star. The double blackbody is found to describe the data best (χ 2 r = 1.04 for 473 dof) resulting in temperatures of kT 1 = 0.57 +0.07 −0.05 keV and kT 2 = 0.32 +0.02 −0.03 keV with R 1 = 0.59 +0.10 −0.09 D 1 km and R 2 = 0.32 +0.02 −0.03 D 1 km for the radii of the emitting areas. It suggests that the first thermal component is emitted from a small fraction of the neutron star surface since the expected radius of a neutron star is ∼ 10 − 15 km. The second thermal component emission can be due to the heating of polar caps by relativistic particles which are accelerated in the magnetosphere or to anisotropic heat conductivity in the neutron star crust (e.g. Becker & Aschenbach 2002). The fit yields a column density N H = 0.47 ± 0.03 × 10 22 cm −2 that is compatible with the values found in the central part of the remnant.
These results are globally compatible with those of Lazendic et al. (2003).

On the nature of the central source
To say whether the neutron star is associated with SNR RX J1713.7-3946 or not, we must compare the results of our spectral analysis with the values found for other types of compact objects found in other SNRs on one hand. On the other hand, the absorbing column value derived from the best-fit model must be consistent with the absorbing column density found for the diffuse emission of the remnant in the central region.
The spectral properties (cf. Table 6) of the central point source 1WGA J1713.4-3949 are very similar to those of the Compact Central Objects (CCOs) in "Vela Junior" (SNR G266.1-1.2 or RX J0852-4622), in Cas A and in Puppis A (Becker & Aschenbach 2002, Pavlov et al. 2004 In addition, the absorbing column N H ∼ 0.5 × 10 22 cm −2 of the diffuse emission around central source 1WGA J1713.4-3949 (top panel of Fig. 4 and Table 7) is Fig. 11. Optical image between 2300 and 3000Å (Second Generation Digitized Sky Survey in red color at 450 − 600 THz) overlaid with the contour of the 2-10 keV X-ray image of SNR RX J1713.7-3946. The scaling is linear with the higher cut fixed at 5σ of the mean optical brightness value.
consistent with the 0.47 ± 0.03 10 22 cm −2 found for the two-component blackbody model ( Table 6). Note that this is also true considering the bremsstrahlung model which leads to a still acceptable N H ∼ 0.7 × 10 22 cm −2 (Table  5) but is weakened adopting the blackbody+power-law model which yields a too high (almost a factor 2) column density N H ∼ 0.9 × 10 22 cm −2 ( Table 6).
As a conclusion, it is highly probable that the point source 1WGA J1713.4-3949 is the compact relic of the supernova progenitor ranking RX J1713.7-3946 in the category of type II supernovae with strong particle acceleration.

Interaction with clouds and distance estimate
The mapping of the CO emission at 8 ′ resolution shown in SL99 has revealed the presence of three dense and massive clouds probably belonging to the same molecular complex of a Galactic arm at 6 kpc. An enhanced intensity ratio of the two lowest rotational transition of the CO molecule was found in one of these clouds located north to RX J1713.7-3946 suggesting that the SNR shock front is hitting this cloud. Based on this interaction, it was concluded that the remnant should be at a distance of 6 kpc.
On the other hand, XMM-Newton has revealed a positive correlation between the absorbing column density and the synchrotron X-ray brightness (see Fig. 6), which strongly suggests that the brightest regions of the SNR interact with part of the line-of-sight absorbing material (cf. Sect. 3.2). In particular, this is the case of the bright southwestern region which exhibits the strongest absorp- Fig. 12. Correlation plot between the absorbing column N H (in units of 10 22 cm −2 ) and the normalized optical brightness. To obtain this plot, the optical image shown in Fig. 11 has been cut following the absorbing column map grid. Each point corresponds to the value found for a particular pixel of the mapping grid (Sect. 2.3.4) and the error bars are given at a 90% confidence level. The symbols are defined as in Table 1. tion in the line-of-sight. At 6 kpc, however, the three dense and massive CO clouds are observed in the north of the remnant; none is observed in the SW (SL99) where the highest absorption and a clear interaction are found with XMM-Newton.
Independent information on the absorption in the lineof-sight of SNR RX J1713.7-3946 can be obtained from a map of integrated optical star light. The optical star light is sensitive to the whole line-of-sight, but it should somehow reflect column density variations. Figure 11 shows the red DSS optical image. The absorption morphology in the optical is quite different from the 6 kpc CO emission morphology presented in SL99 but similar to the X-ray absorption map. The optical brightness appears to be high in the X-ray faint CE and seems particularly weak in the X-ray bright SW region. If we extract the optical brightness and N H for each pixel of the mapping grid (Sect. 2.3.4), we obtain a one dimension plot of these two quantities. Figure  12 shows the very good inverse correlation obtained between our inferred column density and the optical brightness towards RX J1713.7-3946. This is consistent with the low and high absorptions found in the N H map in the CE and SW, respectively (Fig. 4, left panel). This last result provides additional support to the reliability of the N H distribution derived from the X-ray measurements.
The disagreement between the absorption morphology towards RX J1713.7-3946 and the molecular gas distribution at 6 kpc has an immediate consequence either on the distance of the SNR or on the nature of the material interacting with the SNR. If RX J1713.7-3946 is really located at 6 kpc, this implies that the remnant interacts with an interstellar component of different nature than molecular gas, as Hi gas for instance. In the other alternative, if RX J1713.7-3946 is really impacting molecular clouds, this im-plies that RX J1713.7-3946 is situated at a distance lower or larger than 6 kpc. Because there is no evidence that the absorption toward the SNR is strongest in the north contrary to what is observed in the CO map integrated up to 6 kpc (Bronfman et al. 1989, SL99), it indicates that the remnant must be located at a distance lower than 6 kpc.
An indirect estimate of the density of the absorbing material located in the SW can be obtained from the local variations of column density in the interacting regions. We assume that the angular size θ abs of the absorbing material is the same along the line-of-sight and in projection on the sky. If it is so, the density of the absorbing matter n abs is given by: (2) where D is the SNR distance and ∆N H is the variation of the absorbing column density between two regions separated by angular size θ abs . Typically in the SW (see left panel of Fig. 4), θ abs ∼ 15 ′ and ∆N H ≃ 0.4 × 10 22 cm −2 .
This gives n abs ≃ 300 D 1 kpc −1 cm −3 which corresponds to 50 cm −3 for a distance of the SNR of 6 kpc or 300 cm −3 for 1 kpc, a typical value for cold interstellar molecular clouds (e.g. Tatematsu et al. 1990, Wilner, Reynolds & Moffett 1998. The pending question is then to determine the nature and the distance of the matter interacting with RX J1713.7-3946. Since the strongly absorbed regions are mainly found at the limit or around the SNR (particularly in the NW and SW) and not inside, a possible scenario would be that the remnant is within a wind blown bubble created by the precursor star, and that the SN shock would be presently encountering the wall of such a bubble, as previously suggested (SL99, Ellison et al. 2001). However, the optical image shows that the thickness of the regions of lowest optical brightness found at the edges of the remnant (Fig. 11) is much larger than that expected from a stellar wind shell. Then, the most plausible picture is that of preexisting interstellar clouds, possibly pushed by the progenitor stellar wind, being struck by the SNR shock front. This can be tested from CO and Hi observations.
We have reanalyzed the behaviour of the absorbing material across the SNR RX J1713.7-3946 based on data from the CfA CO survey ) and from the IAR High Sensitivity Hi survey (Arnal et al. 2000). The angular resolution of the CO J:1-0 data is 8 ′ .8, the grid sampling of 7 ′ .5 and the velocity resolution of 1.3 km s −1 (note that although Dame et al.'s 2001 new survey improves the sampling interval to half a beamwidth, the region containing the SNR RX J1713.7-3946 in the fourth Galactic quadrant, is still based on the coarser sampling as published by Bronfman et al. 1989). For the Hi data the angular resolution and the sampling are 30 ′ and the velocity resolution of 1.03 km s −1 .   Table 7. Values of the accumulative absorbing column density towards the NW, SW and CE as obtained from CO and Hi observations for several LSR radial velocities or equivalently for several distances. The values of N H which are in bold font are those compatible with the Xrays. Figure 13 (top panel) displays the cumulative absorbing column density N H = N (Hi) + 2 × N (H 2 ) as a function of the Local Standard of Rest (LSR) radial velocity, towards the NW (continuous line), the SW (dotted line) and the center (dashed line). The plotted curves to the NW and SW were obtained from the average over ∼ 8 ′ areas centered near α J2000 =17h11m45s, δ J2000 =-39 • 32 ′ 45 ′′ and near α J2000 =17h11m50s, δ J2000 =-39 • 57 ′ 00 ′′ , respectively. The third N H curve was calculated for a line-of-  Table 7 gives the value of the cumulative absorbing column density towards the NW, SW and center for several LSR radial velocities or equivalently for several distances (only the closer distance is provided).
The line-of-sight towards l = 347.5 • basically crosses three galactic arms: Sagittarius, at an approximate distance of 1.4 kpc (V LSR ∼ −10 km s −1 ), Scutum-Crux at about 3.7 kpc (V LSR ∼ −34 km s −1 ), and Norma at about 5.3 kpc (V LSR ∼ −66 km s −1 ) (Georgelin & Georgelin 1976). Figure 14 displays two average spectra obtained towards the NW (continuous line) and towards the SW (dotted line), where, in addition to the local gas, the three important galactic components are evident.
From the X-ray absorption map (see Fig. 4, left panel), the absorbing column density was found to be higher in the SW (N H ∼ 0.8 − 1.1 × 10 22 cm −2 ) than in the NW (N H ∼ 0.6 − 0.7 × 10 22 cm −2 ). The N H curve obtained to the SW from CO and Hi observations, lies above the curve obtained for the NW region only for velocities below V LSR ∼ −23 km s −1 as shown in Fig. 13 (top panel). Although the limited angular resolution of the available Hi data "dilutes" small spatial variations and reduces the ∆N H between the two regions, this sets an upper limit on the SNR distance at ∼ 2.9 kpc using the galactic rotation curve (Fig. 13, bottom panel).
Besides, the only negative LSR velocities for which the CO emission is clearly higher in the SW than in NW are found between 0 and −10 km s −1 as revealed by the spectra displayed in Fig. 14. Indeed, it can be noticed the presence of a "shoulder" of higher emission obtained to the SW (dotted line), as compared to the NW spectra (solid line). This constitutes an evidence for the presence of molecular gas in the SW direction. This feature is consistent with the molecular cloud recently detected by Fukui et al. (2003) in the SW (between −3 and −11 km s −1 ) on the basis of 2.6 mm CO J:1-0 observations carried out with NANTEN telescope. The range of LSR velocities corresponding to the upper limit distance of 1.6 kpc, using the circular galactic rotation model (Fich et al. 1989), is complety compatible with the allowed distance range of ∼ 0.9 − 1.7 kpc.
A search for ground-state satellite line of the hydroxyl molecule at 1720.53 MHz was carried out towards two spots in G347.3-0.5 where both, X-rays and CO observations, suggested strong shock/molecular cloud interaction. The positive detection of OH (1720 MHz) masers permits: to confirm the existence of such interaction and to constrain the physical parameters of the shock-excited gas on the basis of the very specific conditions required to excite the 1720 MHz OH masers.
The results of the search for OH (1720 MHz) masers were negative. This fact is probably due either to the fact that the shocked molecular clouds have a volume density smaller than required to excite the masers, or the shock front was not of a non-dissociative C-type (Wardle 1999).
The X-ray analysis of SNR RX J1713.7-3946 with XMM-Newton has revealed a strong absorption in the SW, even stronger than in the NW. The absorption derived from the optical star brightness is in agreement with the X-ray derived map. This high absorption found in the SW was unexpected with regard to the morphology of the CO distribution at 6 kpc. Moreover, the positive correlation between X-ray absorption and the X-ray brightness is indicating that the shock front of RX J1713.7-3946 is impacting the molecular clouds responsible for the absorption. The CO and Hi observations show that the inferred cumulative absorbing column densities are in excellent agreement with the X-ray findings in different places of the remnant provided that the SNR lies at a distance of 1.3 ± 0.4 kpc, probably in the Sagittarius galactic arm. Other indications of an interaction between the SNR blast-wave and molecular clouds are provided by the existence of an excess of fast-moving Hi gas (Koo et al. 2004) and possibly, the existence of a broad CO profile in LSR velocity found in a coincident bright X-ray feature in the SW at 1 kpc (Fukui et al. 2003). Based on the absorption morphology and the probable association between the SNR and the molecular clouds, the distance of SNR RX J1713.7-3946 is revised to 1.3 ± 0.4 kpc instead of the commonly accepted value of 6 kpc.

Particle acceleration
The asymmetry in the synchrotron X-ray emission morphology is the result of the asymmetry in the ambient density around the remnant. In the strongly absorbed bright SW rim, we have established that the shock front of RX J1713.7-3946 is interacting with molecular clouds as suggested by the positive correlation between X-ray absorption and X-ray brightness. In the weakly absorbed faint SE region, the same correlation suggests that the shock front of RX J1713.7-3946 is propagating into a tenuous ambient medium which is consistent with the larger extent of the SNR than in the SW. In Sect. 3.2, we have looked at the spatial modifications (in the SW and SE) of particle acceleration at the shock in RX J1713.7-3946 through the photon spectral index changes. The regions where the shock impacts molecular clouds have a steeper spectrum (Γ ∼ 2.3−2.7) than those where the shock propagates into a low density medium (Γ ∼ 1.9 − 2.2).
Qualitatively, we can interpret such modifications in the synchrotron spectrum by considering a modification of the acceleration in the interaction region with dense molecular clouds. To investigate that effect, we have used a simple model of nonlinear diffusive shock acceleration (Berezhko & Ellison 1999 varying the upstream density n 0 , the shock speed V sk and the magnetic field B 0 for an injection (i.e. fraction of gas particles which are accelerated) set to 10 −3 . As the SNR blastwave runs into a high density medium such as molecular clouds, its velocity is expected to diminish according to the ram pressure conservation (n 0 V 2 sk constant) and the upstream magnetic field is expected to increase following the square root of the density (B 0 n −1/2 0 constant). The result of the density increase is to lead to a steepening of the synchrotron spectrum but no increase in brightness. Therefore, this cannot explain the observed difference in X-ray brightness between the SE and the SW (for the same magnetic field orientation). However, if the shock is preferentially tangential in the region of high density but normal in the region of lower density, the X-ray brightness can be strongly enhanced due to the compression of the downstream magnetic field on condition that the density contrast between the two regions is weak (not more than a factor 4). Then, it is possible to explain stronger X-ray brightness and steeper synchrotron spectrum in the SW than in the SE if the magnetic field is tangential in the SW but normal in the SE and if the acceleration process occured before the shock entered the dense molecular clouds.
A modeling of these effects is beyond the scope of this paper. Also, further smaller scale spectral studies with the Chandra satellite would help to investigate the origin of the steepening of the spectrum in the SW.

SNR age and energetics
The study of the thermal emission with XMM-Newton has lead to an upper limit on the density in the ambient medium of n 0 ≤ 2 × 10 −2 D −1/2 1 cm −3 where D 1 is the distance in units of 1 kpc (cf. Sect. 3.3). If the density derived from the emission measure is indeed representative of the density in the ambient medium (which is not the case in the radiative phase), an estimate of the swept-up mass is given by: where θ s is the angular size of the blast-wave radius. This leads to an upper limit on the swept-up mass of 1 M ⊙ for a distance of 1 kpc and θ s = 32 ′ . If RX J1713.7-3946 is the remnant of a type II explosion, the initial mass of the progenitor is expected to be above 5 − 8 M ⊙ (Kennicutt 1984) and then the mass of the ejected matter of a few solar masses. This is well above the derived swept-up mass leading to the conclusion that RX J1713.7-3946 must be in the free expansion phase.
The hypothesis proposed by Wang et al. (1997), based on historical records, that RX J1713.7-3946 is the remnant of the supernova that exploded in AD 393 would imply an age of t s = 1600 years. Assuming selfsimilar solutions (Chevalier 1982), we can estimate the shock velocity given by v s = 2.85 × 10 2 (n − 3)/(n − s) (D/1 kpc) (θ s /1 ′ ) t s /10 3 yr −1 km s −1 where n and s are the power-law indexes of the initial outer density profile in the ejecta and of the initial ambient medium density profile, respectively. n ∼ 9 − 12 is typical of SN II, s = 0 corresponds to an homogeneous ambient medium, s = 2 to an ambient medium modified by the progenitor wind. For n = 9, v s ≃ 3.8 × 10 3 km s −1 if s = 0 and v s ≃ 4.9 × 10 3 km s −1 if s = 2.
The determination of the kinetic energy of the explosion provides a strong argument for or against the hypothesis that RX J1713.7-3946 is the remnant of the AD 393 stellar explosion. In the framework of spherical self-similar model (Chevalier 1982), we can estimate the kinetic energy of the explosion as a function of the SNR age for different progenitor masses and ambient medium densities. The energy curves depend on three parameters among which n and s for a given progenitor mass (Decourchelle & Ballet 1994). The power-law index n of the ejecta density profile was fixed to 9. The power-law index s of the ambient density profile was fixed to 0 to simulate a flat ambient medium which seems to be the situation in the SE. The kinetic energy of the explosion is then determined by the relation: where r s is the shock radius set to 10 pc, ρ 0 = 1.4 m H n 0 and t s the SNR age. The parameters A and ǫ are determined by solving the self-similar hydrodynamical equations (Chevalier 1982). Energy curves are illustrated in Fig. 15 for supernova ejecta masses of M ej = 10, 15 and 20 M ⊙ representative of type II supernova and for different ambient medium densities between 1×10 −4 and 5×10 −3 cm −3 . If we adopt the standard 10 51 erg as the energy of the explosion, it shows that an ambient medium density of 10 −3 cm −3 is appropriate for ejecta masses between 10 M ⊙ and 15 M ⊙ if we consider a remnant of 1600 years old. At this age, higher density of the ambient medium (5 × 10 −3 cm −3 ) leads to a higher kinetic energy of the explosion, still possible, whereas below 10 −4 cm −3 , the kinetic explosion energy begins to be too weak If we consider the effect of efficient cosmic-ray acceleration on the structure of the blast-wave as described in Chevalier (1983), energy curves move a bit higher and would only very slightly widen the range of allowed ambient medium densities. Hence, these results based on energy argumentation are consistent with the picture of a young SNR expanding in a low density medium.

Evolutionary stage
Although the star formation process is complex, there is observational evidence that massive stars are born in the densest parts of molecular clouds. The presence of molecular clouds around the northern and southwestern parts of RX J1713.7-3946 strongly suggests that its progenitor star is indeed born in such conditions. Evolved SNRs interacting with molecular clouds are expected to be efficient electron accelerators (Bykov et al. 1999) as observed in RX J1713.7-3946. Nevertheless, these molecular clouds appear to be located at the periphery of the SNR at ∼ 7 pc away from the center, assuming a distance of 1 kpc (here and hereinafter). Massive stars are effectively known for initiating and maintaining powerful stellar winds during their lives so that they can move the molecular clouds from the vicinity of the progenitor star. These winds may have homogeneized the immediate environment of the progenitor star, explaining the extremely low mean density inside the SNR (McKee et al. 1984, Chevalier 1999. In the standard picture of a stellar wind interacting with the interstellar medium, a four-zone structure (Castor et al. 1975, Weaver et al. 1977) settles very rapidly, consisting of: the innermost hypersonic stellar wind characterized by a density profile ρ ∝ r −2 , a hot shocked stellar wind region, a shell of shocked interstellar medium, and the ambient interstellar gas. After a few thousand years, the swept-up shocked interstellar gas collapses into a thin cold (T ∼ 10 4 K) shell strongly affected by radiative losses while the low-density (< 1 cm −3 ) shocked stellar wind region occupying most of the bubble volume, remains very hot (T > 10 6 K).
A possible interpretation of our observation of RX J1713.7-3946 is that the SN event occured in a windblown bubble. In the SE, the shock front of the remnant is still propagating into the low-density hot shocked wind (n 0 ≤ 0.02 cm −3 ). In the SW, the blast-wave of RX J1713.7-3946 has already progagated into the shocked stellar wind and is impacting the cold and dense (∼ 300 cm −3 , see Eq. (2)) shell of shocked ambient material. The remnant could be in the free expansion phase in the SE whereas in the radiative phase in the SW.
In this scenario, the size of the bubble provides an estimate of the strength of the wind and hence of the progenitor mass. From Fig. 2, the position of the bright X-ray rims in the NW and SW (r b ∼ 6 − 9 pc, taking into account the projection effects) provides an estimate of the radius of the bubble r b , assuming that the molecular clouds have been removed from the vicinity of the star by the progenitor stellar wind.
wind-bubble size where the molecular clouds are located with those estimated from observations of stars of different masses. Table 8 shows two cases for which the bubble radii (5.3 and 10.0 pc) are close to the 6 − 9 pc values estimated in the NW and SW for RX J1713.7-3946. This implies that the mass of RX J1713.7-3946's progenitor must be between 12 and 16 M ⊙ in the main sequence evolution phase, which corresponds to a B-star whose lifetime on the main sequence is about 10 7 years. In Table 8, the mechanical wind luminosity L w ≡Ṁ loss v 2 w /2 is computed assuming an appropriate wind velocity v w = 700 km s −1 as chosen in Chevalier (1999). The maximum size of the wind bubble is given by (in pc): r b = 10.7 L w 10 34 erg s −1 1/3 p/k 10 5 K cm −3 −1/3 t ms 10 7 yr 1/3 (5) assuming that the surrounding interclump medium exerts a pressure p/k = 10 5 K cm −3 , typical of cloud (Chevalier 1999).
From the characteristics of B-stars, we can estimate the mean density in the wind before the supernova event and compare it to the density derived from the thermal X-rays (see Sect. 3.3). The mass lost during the lifetime of the progenitor diluted in the volume of the bubble gives an estimate of the mean number density n b of the windbubble (in cm −3 ): n b = 6.9 × 10 −4 Ṁ loss 10 −8 M ⊙ yr −1 t ms 10 7 yr r b 10 pc −3 (6) For M ⋆ = 12 and 16 M ⊙ (cf. Table 8), we obtain a mean density of 3.6 × 10 −3 and 3.7 × 10 −3 cm −3 , respectively. This is consistent with the ∼ 2 × 10 −2 cm −3 upper limit found for the thermal emission in the X-rays. Note that the mass lost in the red giant phase was not taken into account because it produces a local excess in density only in the vicinity of the progenitor star.
The proposed scenario in which the SNR blast wave expands into a low-density shocked wind medium in the SE and strikes molecular clouds in the SW appears to be a possible physical picture since the density derived from the thermal X-ray emission measures can be found to be consistent with the expected density in the pre-SN shocked stellar wind. However, it is hard to understand why the shell of molecular clouds shocked by the SNR blast-wave does not produce copious thermal X-ray emission since the X-ray brightness is proportional to the square of the density. The regions of enhanced X-ray emission as observed in the SW are those where there is a strong increase of the synchrotron X-ray emission (correlated with an increased of absorbing column density), but where there is no indication of any increase of the thermal X-ray emission.
However, if the SW rim of RX J1713.7-3946 is in the radiative phase, then we do not expect thermal X-ray emission from the shocked interstellar clouds. It would also explain why the thermal X-ray density does not follow the variations of the ambient density. A high density of 300 cm −3 and a temperature of 10 5 K in the shocked interstellar clouds yield a cooling time of ∼ 25 years.
However, in SNRs where the shock is radiative, as observed in the Cygnus Loop, bright optical filaments are generally observed (Levenson et al. 1996). The optical emission is attributed to the cooling of the gas behind the shock front. This cooling gas emits many optical lines including hydrogen Balmer lines with a temperature of ∼ 10 4 K. If RX J1713.7-3946 is indeed in the radiative phase where the blast wave has reached the cavity wall and is interacting with dense matter, such optical filaments might be visible, although the absorption along the lineof-sight is much stronger in RX J1713.7-3946 (N H ∼ 10 22 cm −2 ) than in the Cygnus Loop SNR (N H ∼ 4 × 10 20 cm −2 ).
Finally, we could expect an enhancement of the thermal X-ray emission in the shocked stellar wind just behind the shocked shell. Indeed, a blast wave impacting on a dense cloud gives rise to both a reflected shock propagating back and a transmitted shock propagating into the cloud (McKee & Cowie 1975). The shocks that are reflected off molecular clouds surfaces are expected to propagate back into the hot, previously shocked gas, further heating and compressing it (e.g., Levenson et al. 1996). Then, enhanced X-ray thermal emission is expected (Hester & Cox 1986). However, the density in the shocked gas is still small so that its contribution might be negligeable.

Conclusion
The XMM-Newton observation has yielded several new results on SNR RX J1713.7-3946. Observational facts and their astrophysical interpretations are summarized here: -The X-ray bright central point source 1WGA J1713.4-3949 detected at the center of SNR RX J1713.7-3946 shows spectral properties very similar to those of the Compact Central Objects found in SNRs (e.g., Vela Junior) and which are consistent with the absorbing column density of the central diffuse X-ray emission arising from the SNR. It is highly probable that the point source 1WGA J1713.4-3949 is the compact relic of RX J1713.7-3946's supernova progenitor.
-The X-ray mapping of the absorbing column density has revealed strong variations (0.4 × 10 22 cm −2 ≤ N H ≤ 1.1 × 10 22 cm −2 ) and, particularly, an unexpected strong absorption in the southwest. The column density variations are well reflected by the extinction in the map of integrated optical star light. The strong positive correlation between the X-ray absorption and the X-ray brightness along the western rims suggests that the shock front of RX J1713.7-3946 is impacting molecular clouds there. -CO and Hi observations show that the inferred cumulative absorbing column densities are in excellent agreement with the X-ray measurements in different places of the remnant only if the SNR is placed at a distance of 1.3 ± 0.4 kpc, probably in the Sagittarius galactic arm, and not at 6 kpc as previously claimed. An excess in the CO emission found in the southwest at ∼ 1 kpc strongly suggests that molecular clouds produce the enhancement in absorption. A search for OH masers in the southwestern region has been unsuccessful, possibly due to the low density of the clouds. -The X-ray mapping of the photon index has revealed strong variations (1.8 ≤ Γ ≤ 2.6). The spectrum is steep in the faint central regions and flat at the presumed shock locations, particularly in the southeast. However, the regions where the shock strikes molecular clouds have a steeper spectrum than those where the shock propagates into a low density medium. -The search for the thermal emission in RX J1713.7-3946 has been still unsuccessful leading to a number density upper limit of 2 × 10 −2 cm −3 in the ambient medium. This low density in the ambient medium yields to a reasonable kinetic energy of the explosion provided that the remnant is less than a few thousand years old. It can be reconciled with the high density in molecular clouds if the remnant is in the radiative phase where the SNR shock encounters a dense ambient medium whereas it is in the free expansion phase elsewhere. -RX J1713.7-3946's progenitor mass is estimated to lie between 12 and 16 M ⊙ based on a scenario involving the effect of stellar wind of the progenitor star.