Ionized X-ray winds in the radio galaxy Mrk 6
Abstract
Active galactic nucleus (AGN) outflows including jets and ionized winds have been key phenomena such as jet collimation and AGN feedback to the host galaxy in astrophysics. Radio galaxies, a type of AGN with misaligned jets, have provided valuable insights into the properties and relationships of these outflows. However, several aspects regarding AGN outflows remain unclarified, such as the relationship between jets and ultrafast outflows (UFOs) and the differences between the properties of radio-loud AGN disk winds and radio-quiet AGN ionized winds. To clarify these aspects, radio galaxies containing UFOs and warm absorbers (WAs) must be studied. Currently, both UFOs and WAs have been reported in only two radio galaxies, 3C 390.3 and 4C 74.26. To enhance our understanding of the connection between the jets and ionized winds, we conducted a study on Mrk 6, a potential candidate for the third case of a radio galaxy displaying these characteristics. Using X-ray spectra obtained from XMM-Newton, we performed photoionization modeling using the SPEX code. The best-fit model analysis results revealed the presence of a UFO component with a relatively low ionization parameter (Fe xix/xviii lines blueshifted by ) and a WA component with an outflow velocity of . To further confirm the nature of the UFO and WA in Mrk 6, we simulated the X-ray imagining and spectroscopy mission spectra.
keywords:
galaxies: active — X-rays: galaxies1 Introduction
Active galactic nuclei (AGNs) accrete mass onto supermassive black holes (SMBHs) while rejecting outflow as jets or ionized winds. Ionized winds are characterized by their ionization parameter , hydrogen column density , and outflow velocity, and are classified as ultrafast outflows (UFOs, Tombesi et al., 2010b), warm absorbers (WAs, Crenshaw et al., 2003; Halpern, 1984) and obscures (e.g. Kaastra et al., 2014; Mehdipour et al., 2017; Mao et al., 2022). Ionization is defined as , where ( Ryd) denotes the ionizing luminosity, represents the disk wind hydrogen number density, and indicates the distance from the source (Tarter et al., 1969). Generally, UFOs are observed mainly as blueshifted absorption lines at keV (Tombesi et al., 2010b). The UFO ionization parameter is , the hydrogen column density is , the outflow velocity is (Laha et al., 2021), and the kinetic power of the mass outflow is of Eddington luminosity (e.g. Gofford et al., 2015; Tombesi et al., 2012a). As such, UFOs are estimated to exist in AGNs (Tombesi et al., 2010b, 2011; Gofford et al., 2013). Based on simulations (e.g. Nomura et al., 2016),UFOs are launched at times the Schwarzschild radius , and from observations (e.g. Gofford et al., 2015), UFOs originate at . Although prior research suggested that UFOs may originate from the jet base (Ghisellini et al., 2004), the launch radius of UFOs remains unclear and the relationship between the jet and UFO has not yet been characterized. In general, WAs are detected mainly as blueshifted absorption lines or edges in the soft X-ray spectrum. The WAs ionization parameter is , the hydrogen column density is , the outflow rate is (Sako et al., 2001; Laha et al., 2016). WAs exist at pc, and their hydrogen number density has been estimated to be (Laha et al., 2021, and references therein). Recently, the possibility of a new classification of entrained ultra-fast outflows (E-UFO, e.g. Serafinelli et al., 2019, and references therein), which have fast outflow velocity, low-ionization, and low-density, has been suggested.
AGNs can be classified by radio intensity (or radio loudness). The radio intensity is defined based on Kellermann et al. (1989) as using the radio flux () and the optical flux (Schmidt & Green, 1983). For radio-loud AGN, , and for radio-quiet AGN, . Although the spectral energy distribution (SED) of radio-loud and radio-quiet AGNs are similar across the infrared to X-ray spectrum, their radio and gamma radiations differ owing to the strong relativistic jets emitted by radio-loud AGNs (e.g. Urry & Padovani, 1995; Kayanoki & Fukazawa, 2022). In principle, radio galaxies are radio-loud AGNs with jets misaligned to the line of sight, which facilitate the observation of the jets as well as the central core in such galaxies. Therefore, radio galaxies are deemed as the most suitable objects for investigating the environment surrounding SMBHs and characterizing the relationship between the jet and accretion disk and that between the jet and disk winds. Recently, several scholars have conducted statistical studies on radio galaxies, e.g., X-ray absorption in radio galaxies in Kayanoki & Fukazawa (2022), UFOs in radio galaxies in Tombesi et al. (2010a), Tombesi et al. (2012b), Gofford et al. (2013), Tombesi et al. (2013) and Tombesi et al. (2014), and WAs in radio galaxies in Reeves et al. (2009), Torresi et al. (2010), Torresi et al. (2012), Tombesi et al. (2016), and Mehdipour & Costantini (2019). Notably, WAs are seldom observed in radio galaxies and might be related to the power of the jets (Mehdipour & Costantini, 2019); however, a direct physical relationship has not yet been established.
Jets are presumed to be magnetically powered; however, the collimation of jets in the vicinity of SMBHs is not understood. A quasi-analytical study suggested that disk winds may promote the collimation of jets near SMBHs (Fukumura et al., 2014; Globus & Levinson, 2016; Blandford et al., 2019). Moreover, as reported by Mehdipour & Costantini (2019), WAs are related to the jet intensity (radio-loudness). However, several questions regarding the AGN outflows remain unsolved, such as the possible variations in the properties of disk winds in radio-loud AGNs and radio-quiet AGNs, including the relationship between jets and UFOs. To resolve these questions, the AGN with observed jets, UFOs, and WAs, i.e., radio galaxies containing both UFOs and WAs, provide vital insights. However, only two radio galaxies have been reported to contain both UFOs and WAs, 3C 390.3 and 4C 74.26 (Tombesi et al. (2014) for UFOs, Mehdipour & Costantini (2019) for WAs). Therefore, understanding the physical state of the disk winds in radio galaxies is crucial for resolving these problems. Mrk 6 is the Seyfert 1.5-type AGN in a lenticular galaxy (Kharb et al., 2014). -band magnitude and 6 cm radio flux () of Mrk 6 is and Jy (Véron-Cetty & Véron, 2010), Mrk 6 can be classified as a radio-loud galaxy (). Specifically, it has a black hole mass of (Doroshenko et al., 2012) with (de Vaucouleurs et al., 1991). In Mehdipour & Costantini (2019), the X-ray data analysis of the European Photon Imaging Camera (EPIC, Strüder et al., 2001; Turner et al., 2001) and Reflection Grating Spectrometer (RGS, den Herder et al., 2001) detectors on the XMM-Newton satellite for 16 radio galaxies reported WAs, among which Mrk 6 exhibited the highest density of WAs. Furthermore, Schurch et al. (2006) reported that Mrk 6 contains disk winds or an ionized skin or an atmosphere of the molecular torus. However, as the UFO of the Mrk 6 has not yet been discussed, this study analyzed the Mrk6 disk winds in more detail. The cosmological parameters (, , ) used in this study are , , .
2 Observation and data reduction
Mrk 6 was observed with XMM-Newton in March 2001, April 2003, and October 2005 (Table 1). In this study, we used the RGS data for the soft X-ray spectra and the positive-negative junction (pn) charged-coupled device (CCD) camera (Strüder et al., 2001) of the EPIC for the hard X-ray spectra, because EPIC/MOS data does not significantly increase the statistics. The XMM-Newton data were reprocessed and analyzed using Science Analysis System version 19.1.0 and the calibration data from current calibration files (CCF) were updated on November 18, 2021. After eliminating the data at instances unsuitable for analysis because of background flares, the exposure durations of the EPIC detector were shorter than those of the RGS detector data (Table 1). Moreover, the observations (IDs 0061540101 and 0305600501) were unsuitable for disk wind analysis, because their extremely short exposure durations produce excessively large errors in the spectral data points. Herein, the analysis was performed on the data with the longest EPIC/pn exposure time, i.e., observation ID 0144230101.
Observation ID | Net Exposure (ks) | Date | ||
---|---|---|---|---|
RGS1 | RGS2 | EPIC/pn | ||
0061540101 | 33.8 | 24.9 | 10.5 | 2001-03-27 |
0144230101 | 34.4 | 34.4 | 31.5 | 2003-04-26 |
0305600501 | 21.0 | 21.0 | 4.49 | 2005-10-27 |
For the EPIC/pn data reduction, the epchain pipeline task was used. In EPIC/pn (), periods with high flaring were filtered out by applying the #XMMEA_EP filters. The EPIC/pn spectra were extracted from a circular region of the center coordinate (06:52:11.6, +74:25:35.1) with a radius of 35.750 arcsec. Moreover, the background was extracted from a neighboring source-free rectangle region of center coordinates (06:52:05.0, +74:22:05.9), length and width of 96 arcsec and 262 arcsec, respectively, and a tilt of 342.29 in the same CCD. Overall, the pile-up was evaluated to be negligible. Thereafter, the response matrices were generated for the EPIC/pn spectra of each observation using the rmfgen and arfgen tasks. For the RGS data, the rgsproc pipeline task was used for data reduction, and the source and background spectra were extracted to generate the response matrices. Moreover, the time intervals with background count rates were filtered out with the Chandra Interactive Analysis of Observations (CIAO, Fruscione et al., 2006) version 4.14 lc_sigma_clip task. To rectify the cross-calibration between the EPIC/pn and RGS spectra, the RGS1 and RGS2 spectra were rescaled by a factor of 0.91 and 0.89, respectively, to correspond to the flux level in the keV energy range.
3 Spectral Fitting and Result
In this study, the spectral analysis was performed using the SPEX code v3.06.01 (Kaastra et al., 2020). For spectral analysis, EPIC/pn spectra were used in the energy range of keV, and the obin command of the SPEX code was applied to perform optimal binning according to Kaastra & Bleeker (2016). The two RGS spectra were utilized for the Å wavelength range and binned by a factor of 2 to ensure that the bin size was approximately 3rd of the RGS energy resolution. In particular, the spectral fitting was performed on the RGS and EPIC/pn spectra using the C-statistic (Kaastra, 2017). In all the fittings, performed hereinafter, all models were multiplied by the neutral absorption model (hot model) to account for the absorption by our Galaxy. The column density of our Galaxy was set to (HI4PI Collaboration, 2016). For the reflection and absorption models, we considered the protosolar abundances of Lodders et al. (2009). The error for each parameter refers to the 68.3% confidence level.
3.1 Search for emission and absorption lines
First, to inspect the emission and absorption lines in the spectra, we analyzed the EPIC/pn spectra in the keV energy band using the power-law (pow) model, the power-law (pow) + emission line model, and the power-law (pow) + emission line (gaus) + absorption line (gaus) model (Table 2). The EPIC/pn spectra in the keV energy band and the power-law (pow) + emission line (gaus) + absorption line model (gaus) are illustrated in Figure 1. The Fe emission line at 6.4 keV and absorption line at keV were observed in Figure 1. In addition to Table 2, the model with the absorption lines displayed a more accurate C-statistic (, this would be a marginal detection) in Figure 1, thereby suggesting the existence of UFOs. The emission line at 6.4 keV was considered to be fluorescence X-rays owing to the neutral or low ionized iron in AGN, whereas the keV absorption lines were considered to be caused by a highly ionized gas.
power-law | power-law | power-law | |
+ emission line | + emission line | ||
+ absorption line | |||
pow | |||
pow Norm | |||
gaus e | - | ||
gaus norm | - | ||
gaus FWHM | - | ||
gaus e | - | - | |
gaus norm | - | - | |
gaus FWHM | - | - | |
Total C-stat / d.o.f. | 149.9 / 83 | 82.4 / 80 | 76.5 / 77 |
Total Expected C-stat | 85.113.1 | 85.113.1 | 85.113.1 |
C-stat | - | -67.5 | -73.4 |
-
•
Fit the EPIC/pn hard X-ray data (5-10 keV). C-stat: Difference in C-statistics between the simple power-law model. gauss model used for obtaining emission and absorption lines. (1) Photon index of the X-ray power-law component. (2) Normalization of the X-ray power-law component in at 1 eV. (3) Line energy in keV. (4) Normalization of Gaussian in . (5) Line FWHM in keV.
![Refer to caption](x1.png)
3.2 Continuum
Furthermore, we investigated the features (photon index and flux) of the continuum component of the X-ray spectra. As the high-energy band is considered to be weakly affected by absorption, the EPIC/pn data were fitted with a model in which the power law was multiplied by the absorption (hot) by the neutral gas of the host galaxy. In addition, as the keV Fe emission lines possibly existed based on the results of the previous subsection, the power-law (pow) + reflection component (refl) multiplied by the absorption (hot) of the neutral gas of host galaxy was fitted to the EPIC/pn data. Moreover, the RGS and EPIC/pn data were fitted with a model that multiplied the modified blackbody component (mbb) + power-law (pow) + reflection component (refl) by the absorption (hot) of the neutral gas in the host galaxy.
Compared to the power-law model with the power-law (pow) + reflection component (refl) model, the C-statistic for the total spectrum (EPIC/pn + RGS spectrum) is more suitable for the power-law model without the reflection component ( C = 100.5). Nonetheless, the C-statistic for the EPIC/pn data spectrum is more appropriate for the model with the reflection component. A comparison of the C-statistics of the EPIC/pn data spectra indicated the considerably superior performance of the model incorporating the reflective component. This is because the reflection model can adequately explain the keV Fe emission lines. However, the power-law (pow) + reflection component (refl) model does not adequately represent the RGS (soft X-ray spectrum). Accordingly, the modified blackbody component (mbb) was added to the continuum. As listed in Table 3, the C-statistic for the total spectrum (EPIC/pn + RGS spectrum) was 2286.8, which is a relatively more accurate result. Based on the aforementioned findings, we adopted the modified blackbody model (mbb) + power-law (pow) + reflection (refl) as the continuum component. In the following subsection, we investigated the suitability of the neutral gas hot model (de Plaa et al., 2004; Steenbrugge et al., 2005) and the ionized gas pion model (Miller et al., 2015; Mehdipour et al., 2016) as the absorbed component.
3.3 Absorption component
We investigated whether the neutral gas (hot) model or the ionized gas (pion) model (disk wind model) is appropriate as the absorber. As discussed in Section 3.2, the M1 model multiplying the neutral gas (hot) model to the continuum did not explain the keV absorption line. Therefore, we adopted the photoionization model pion. The photoionization equilibrium of the pion model is calculated self-consistently using the available plasma routines of SPEX 111 https://spex-xray.github.io/spex-help/models/pion.html#sec-pion . Therefore, the pion model could calculate the SED and SED filtering (the second PION sees a filtered SED) with the photoionization equilibrium simultaneously. However, the photoionization equilibrium of the other photoionization model (e.g. the xabs model) is pre-calculated for a grid of -values by an external code, for instance, Cloudy or XSTAR. For the pion model, the hydrogen density were set to 1 and the broadening velocity were set to 100 . These two values are the default values of the pion model.
First, we added a single pion component to the M1 model (M2 model). Compared to the M1 model, the M2 model improved the C-statistic. However, the ionizing absorption component of the M2 model exhibited soft X-ray absorption and did not represent the keV absorption line. Subsequently, the neutral gas (hot) component in the M2 model was modified to an ionized gas (pion) component (M3 model). Initially, the M3 model was assumed to contain one neutral gas (hot) and two ionized gases. As the column density for the hot model was obtained as zero from the fitting results, the M3 model potentially did not require a neutral gas (hot). Compared to the M2 model, the M3 model improved the C-statistic by 61.7. The M3 model exhibited the soft X-ray absorption (Figure 4) and keV absorption lines (however, the keV absorption lines were weak (Figure 5). Although the 7 keV absorption lines of UFOs are typically considered to be Fe xxv/xxvi absorption lines (e.g. Tombesi et al., 2014). However, the M3 model has a low ionization parameter and it ascribed the keV absorption lines to Fe xix/xviii. Therefore, the M3 model has weak absorption lines. Thus, the parameters of the M3 model were modified and fitted again. The obtained result may be a local minimum (M4 model). Although the M4 model more accurately represented the keV absorption lines (Figure 5), its C-statistic of 36.6 was inferior to that obtained with the M3 model. However, the keV absorption lines can be possibly attributed to the Fe xxv/xxvi absorption lines, the spectra obtained herein suggested that the keV absorption lines are Fe xix/xviii. As such, the data from the X-ray imagining and spectroscopy mission (XRISM) can be utilized to verify whether the keV absorption lines in Mrk 6 are caused by Fe xxv/xxvi or Fe xix/xviii. The fitting results for the M1, M2, M3, and M4 models are summarized in Table 3.
At this point, we tried thawing the covering factor and Fe abundance of the UFO component for both the M3 model and the M4 model. First, we thawed the covering factor () of the UFO component for the M3 model and the M4 model. For the M3 model, the covering factor changed to = 0.7 and the C-stat improved significantly (). On the other hand, for the M4 model, the covering factor did not change significantly. Therefore, we thawed the covering factor model for the M3 model. Next, we thawed the Fe abundance of the UFO component for the M3 model and the M4 model. For the M3 model, the Fe abundance did not change significantly. On the other hand, for the M4 model, the Fe abundance changed to 3.3 solar and the C-stat improved significantly (). As a result, the M3 model has the lowest C-stat value (Table 3) among the above 4 cases. Therefore, the M3 model with a covering factor of 0.7 was deemed to be the best model.
To investigate whether the M3 model and M4 model are not a local minimum, we varied the value of outflow velocity in the range while performing the fitting and checked the C-stat. In this test, both the M3 model and the M4 model look not local minimum. (We could not determine whether the M4 model is just a local minimum or not.) In addition, we made contour maps for the M3 model of the UFO component (Fig 2).
![Refer to caption](x2.png)
![Refer to caption](x3.png)
![Refer to caption](x4.png)
![Refer to caption](x5.png)
We did the likelihood ratio test to estimate the statistical confidence level of the M3 model with PN spectra. The likelihood ratio is (PN C-stat of the M2 model) (PN C-stat of the M3 model) = 183.0 166.8 = 16.2. PN spectra have 162 bins, and then we decided that we could use the likelihood ratio test for PN spectra and approximate the test following the distribution. Since the difference in degrees of freedom between the M2 model and the M3 model is 3, the -value of the M3 model should be a . From the above, the confidence level of the M3 model is .
In addition, to investigate the number of layers of the ionized gas, the M3 model was multiplied by an additional component of the ionized gas (pion) model (M5 model). Similar to the M3 model, the column density of the neutral absorption component was zero in this model. Thus, the neutral absorption (hot) component was not required in the M5 model. Compared to the M3 model, the M5 model improved the C-statistic by only 7.9 with the 3 degrees of freedom. Thus, the third ionized gas (pion) was regarded as insignificant. (Even if the ionized gas (pion) is included within the M4 model was not significant.)
According to the results obtained by the best model (M3 model), an ionized gas exhibited high ionization energy of , high column density of , and an extremely fast outflow velocity of . Thus, this ionized gas component was identified to be a UFO. The other ionized gas has a lower ionization parameter of , lower column density of and a lower outflow velocity of . Thus, this ionized gas component was deemed as a WA.
We evaluated the absorption lines of UFOs and WAs. The UFO absorption lines above 6.0 keV with and soft X-ray WA absorption lines with are illustrated in Figure 7 and Figure 8, respectively. In addition, the information on UFO absorption is likely better identified with presented in Table 4.
M1 | M2 | M3 | M4 | |
pow | ||||
pow Norm | ||||
refl Scal | ||||
hot | - | - | ||
mbb Norm | ||||
mbb t | ||||
pion | - | |||
pion | - | |||
pion | - | |||
pion | - | - | - | |
pion | - | - | - | |
pion | - | - | ||
pion | - | - | ||
pion | - | - | ||
Total C-stat / d.o.f. | 2286.8 / 1727 | 2218.8 / 1724 | 2157.1 / 1722 | 2193.7 / 1722 |
RGS1 C-stat / bins | 1085.0 / 803 | 1054.9 / 803 | 1020.6 / 803 | 1015.6 / 803 |
RGS2 C-stat / bins | 982.7 / 768 | 981.0 / 768 | 969.7 / 768 | 997.8 / 768 |
EPIC/pn C-stat/bins | 220.2 / 162 | 183.0 / 162 | 166.8 / 162 | 180.3 / 162 |
Total Expected C-stat | ||||
RGS1 Expected C-stat | 886.1 | 892.8 | 896.2 | 895.5 |
RGS2 Expected C-stat | 862.6 | 863.2 | 859.0 | 854.5 |
EPIC/pn Expected C-stat | 162.1 | 162.1 | 162.2 | 162.2 |
C-stat |
-
•
Fit the total spectra (0.34-10 keV) with a continuum (modified blackbody model power-law reflection) and the absorption component. The cut-off energy of the power law and reflection was set to 120 keV. C-stat: Compare the C-stat with M1 model. M1: With hot model. M2: With one pion model and one hot model. M3, M4: With two pion models. (1) Photon index of the X-ray power-law component. (2) Normalization of the X-ray power-law component in at 1 keV. (3) Scale parameter of the reflection component. (4) Column density of neutral ISM gas component in the host galaxy of the AGN in cm. (5) Normalization of the modified blackbody in . (6) Temperature of the modified blackbody in eV. (7) Logarithm of the ionization parameter of the ionized wind components in . (8) Column density of the ionized wind components in cm. (9) Velocity of the outflowing ionized wind components in . (10) Covering factor of the pion component. (11) Fe abundance of the pion component.
![Refer to caption](x6.png)
![Refer to caption](x7.png)
![Refer to caption](x8.png)
![Refer to caption](x9.png)
![Refer to caption](x10.png)
![Refer to caption](x11.png)
Absorption line | Rest-frame Energy (keV) | (eV) | Observed Energy (keV) |
---|---|---|---|
M3 model () | |||
Fe xix | 6.4656 | 3.5 | 7.0794 |
Fe xviii | 6.4347 | 3.4 | 7.0456 |
M4 model () | |||
Fe xxvi | 6.9732 | 23.4 | 7.1180 |
Fe xxvi | 6.9517 | 20.4 | 7.0960 |
Fe xxv | 6.7004 | 18.2 | 6.8395 |
-
•
These lines are based on the M3 model () and M4 model () for each. The red dashed line and blue dashed-dotted line in Fig. 7 are the absorption lines shown here.
4 Discussion
4.1 UFO
As reported in prior research, 3C 390.3 and 4C 74.26 are the only two galaxies where both UFOs and WAs have been observed, with Mrk 6 being the third case considered in this study (Table 5). Compared to the Mrk 6 UFOs with 3C 390.3 and 4C 74.26, the ionization was lower, and the column density and outflow velocity are comparable (Table 5). Although, both 3C 390.3 and 4C 74.26 have lower limits for the column density (Tombesi et al., 2014).
Mrk 6 | 3C 390.3 | 4C 74.26 | |
---|---|---|---|
-
•
*: Tombesi et al. (2014). The xabs model was used for fitting. (1) Logarithm of the ionization parameter of the ionized wind components in . (2) Column density of the ionized wind components in cm. (3) Velocity of the outflowing ionized wind components .
The ionization parameter is defined as (Tarter et al., 1969) using the hydrogen density , disk wind distance , and luminosity of the light source . Thus, the distance of the disk wind was calculated as
(1) |
Based on the M3 model result, the UFO ionization parameter was and the non-filtered (before getting absorbed by the UFO pion component) ionized luminosity was calculated to . Therefore, the distance of the UFO from the SMBH could be estimated as
(2) | |||||
This is equivalent to , where the Schwarzschild radius . Assuming that the wind has a higher velocity than the escape velocity. Therefore, the lower limit of the distance of the disk wind was evaluated as
(3) |
using the gravitational constant , SMBH mass , and outflow velocity . Based on the analysis results, the outflow velocity of the UFO was , and therefore, the lower limit of the UFO distance was calculated as . Assuming that the disk wind width did not exceed the distance to the SMBH, the upper limit of the disk wind distance was computed as
(4) |
considering the covering factor of the disk wind (Wang et al., 2022). M3 model assumed a covering factor of , such that the upper limit of the UFO distance was , considering a column density of . Therefore, the UFOs of Mrk 6 originated at , and the range of the UFOs is . Based on observational studies, UFOs are considered to originate in the (e.g. Gofford et al., 2015), and according to simulation research, they originate in the (e.g. Nomura et al., 2016). Given the observational and model uncertainties, here you can say that the range is consistent with what expected from disk winds.
Furthermore, we evaluated the range of hydrogen density from that of UFOs. Based on the definition of ionization, the density of the disk wind can be evaluated as
(5) |
The range of the UFOs was , and thus, the range of the hydrogen density was calculated as .
4.2 WA
WAs are often reported to contain multiple components, e.g., NGC 3516 has 3 WAs (Mehdipour et al., 2010), NGC 7469 has 4 WAs (Mehdipour et al., 2018), NGC 5548 has 6 WAs (Mao et al., 2017), NGC 3783 has 9 WAs (Mao et al., 2019) and NGC 3227 has 4 WAs (Wang et al., 2022). Both 3C 390.3 and 4C 74.26 have two WAs components each (Mehdipour & Costantini, 2019). However, for Mrk 6, the second WA component is not significant. Thus, Mehdipour & Costantini (2019) used only one WA component along with a neutral absorption component. Compared to the Mrk 6 WAs with 3C 390.3 and 4C 74.26, the ionization was low and the outflow velocity was high, however, the column densities did not vary significantly (Table 6).
Mrk 6 | 3C 390.3 | 4C 74.26 | |
---|---|---|---|
, | , | ||
, | , | ||
, | , |
-
•
*: Mehdipour & Costantini (2019). The pion model was used for fitting. (1) Logarithm of the ionization parameter of the ionized wind components in . (2) Column density of the ionized wind components in cm. (3) Velocity of the outflowing ionized wind components .
Based on the results of the analysis with the M3 model, the ionization parameter of WA was and the filtered (after getting absorbed by the UFO pion component and before getting absorbed by the WA pion component) ionized luminosity was calculated to . Therefore, with Equation (1), the distance of WA from SMBH was
(6) | |||||
which is equivalent to . The fitting results exhibited that the outflow velocity of the WA was , and the lower limit of the WA distance was ) using Equation (3). As the column density of the WA was according to the analysis and the coverage factor was assumed to be , the upper limit of the WA distance was calculated to be using Equation (4). Therefore, the WAs of Mrk 6 originated at , and the range of the WAs was . WAs are considered to originate in the pc (Laha et al., 2021), and the present findings are consistent with this result.
Similar to the case with UFOs, the hydrogen density of the WA was calculated as well. The WA distance range was calculated to be , such that the hydrogen density of WA was calculated as using Equation (5). The hydrogen number density of WAs was estimated to be (Laha et al., 2021, and references therein), and the present result was wider than the range.
4.3 XRISM simulation
As the EPIC/pn instrument of XMM-Newton exhibited a low energy resolution, the absorption lines could not be clearly observed (Figure 3) and the ionized wind parameters (ionization parameter, hydrogen column density, and outflow velocity) yielded large errors (Table 3). However, XRISM/Resolve (Tashiro et al., 2020; XRISM Science Team, 2022) provides an ultrahigh energy resolution, which facilitates the prominent observation of the absorption line and reduces the uncertainties of the parameters and XRISM to verify whether Mrk 6 contained only one WA component. Thereafter, we simulated the observation in XRISM/Resolve based on the M3 model, which was determined to be the best-fit model. The XRISM/Resolve response file was used with an energy resolution of 7 eV version 2019a222https://xrism.isas.jaxa.jp/research/proposer/obsplan/response/index.html. For the simulations and fitting, the SPEX code was used, and the exposure time was set to 150 ks. In the analysis, the energy range was set to keV. For comparison, we simulated the M4 model to obtain results with Fe xxv/xxvi absorption lines at keV.
Based on the findings listed in Table 7, the errors of the disk wind parameters in both the M3 and M4 models were smaller than those reported in Table 3. In particular, the error in the column density , which is a vital parameter for the significant detection of the disk wind, was reduced. The keV absorption lines can be more prominently observed (Figure 9). In addition, the soft X-ray spectra error bar was more accurate (Figure 10). Therefore, XRISM/Resolve can more accurately detect the disk wind. As observed in Figure 9, the Fe xix/xviii absorption lines are weak, and the Fe xxv/xxvi absorption lines are strong. Furthermore, the Fe xxv/xxvi and Fe xix/xviii absorption lines appeared prominently at various energies. Thus, we can employ the XRISM/Resolve to decide whether the keV absorption lines in Mrk 6 can be attributed to Fe xxv/xxvi or Fe xix/xviii.
M3 | M4 | |
pow | ||
pow Norm | ||
refl Scal | ||
mbb Norm | ||
mbb t | ||
pion | ||
pion | ||
pion | ||
pion | - | |
pion | - | |
pion | ||
pion | ||
pion | ||
Total C-stat / d.o.f. | 4095.6 / 3935 | 3920.9 / 3935 |
Total Expected C-stat |
-
•
Fit the total spectra (0.34-10 keV) with a continuum (modified blackbody model power-law reflection) and the absorption component. The cut-off energy of the power law and reflection was set to 120 keV. (1) Photon index of the X-ray power-law component. (2) Normalization of the X-ray power-law component in at 1 keV. (3) Scale parameter of the reflection component. (4) Normalization of the modified blackbody in . (5) Temperature of the modified blackbody in eV. (6) Logarithm of the ionization parameter of the ionized wind components in . (7) Column density of the ionized wind components in cm. (8) Velocity of the outflowing ionized wind components) in . (9) Covering factor of the pion component. (10) Fe abundance of the pion component.
![Refer to caption](x12.png)
![Refer to caption](x13.png)
5 Summary
This study analyzed the disk wind of the radio galaxy Mrk 6 using XMM-Newton archival data with the SPEX code. Based on the results of the M3 model, the existence of both UFOs and WAs was suggested in Mrk 6. In addition to the radio galaxies 3C 390.3, 4C +74.26, Mrk 6 is the third radio galaxy in which both UFOs and WAs have been observed. Compared to the Mrk 6 disk wind with 3C 390.3 and 4C 74.26, the UFOs exhibited low ionization, however, the UFO and WA parameters are broadly consistent with the ones reported for the other two radio galaxies. The distance from SMBH of UFOs and WAs are and , respectively. These distances were consistent with the typical ones. The M3 model attributes the 7 keV absorption lines to Fe xix/xviii. For comparison, the M3 model was re-fitted with varying parameters to obtain a model that attributed the 7 keV absorption lines to Fe xxv/xxvi (M4 model) and it may represent a local minimum of the C-stat. Although the keV absorption lines can be attributed to Fe xxv/xxvi absorption lines, the spectra derived herein suggested that the keV absorption lines are actually Fe xix/xviii absorption lines. Finally, we performed simulations of observations with the XRISM satellite and determined that the errors in the pion model parameters were extremely small using XRISM. Therefore, XRISM can be used to detect significant disk winds. In addition, based on XRISM simulations we investigated whether the keV absorption lines in Mrk 6 were caused by Fe xxv/xxvi or Fe xix/xviii.
Acknowledgements
The authors thank the anonymous referee for their useful comments to improve the quality of this work. This work is based on observations obtained with XMM-Newton. This work was supported by JST, the establishment of university fellowships toward the creation of science technology innovation, Grant Number JPMJFS2129.
Data Availability
A reproduction package is available at Zenodo DOI: 10.5281/zenodo.10394008. This package includes data and scripts used to reproduce the fitting results and figures presented here.
References
- Blandford et al. (2019) Blandford, R., Meier, D., & Readhead, A. 2019, ARA&A, 57, 467. doi:10.1146/annurev-astro-081817-051948
- Crenshaw et al. (2003) Crenshaw, D. M., Kraemer, S. B., & George, I. M. 2003, ARA&A, 41, 117. doi:10.1146/annurev.astro.41.082801.100328
- Kaastra et al. (2020) Kaastra, J. S., Raassen, A. J. J., de Plaa, J., et al. 2020, Zenodo
- de Plaa et al. (2004) de Plaa, J., Kaastra, J. S., Tamura, T., et al. 2004, A&A, 423, 49. doi:10.1051/0004-6361:20047170
- den Herder et al. (2001) den Herder, J. W., Brinkman, A. C., Kahn, S. M., et al. 2001, A&A, 365, L7. doi:10.1051/0004-6361:20000058
- de Vaucouleurs et al. (1991) de Vaucouleurs G., de Vaucouleurs A., Corwin H. G., Buta R. J., Paturel G., Fouque P., 1991, rc3..book
- Doroshenko et al. (2012) Doroshenko, V. T., Sergeev, S. G., Klimanov, S. A., et al. 2012, MNRAS, 426, 416. doi:10.1111/j.1365-2966.2012.20843.x
- Fruscione et al. (2006) Fruscione, A., McDowell, J. C., Allen, G. E., et al. 2006, Proc. SPIE, 6270, 62701V. doi:10.1117/12.671760
- Fukumura et al. (2014) Fukumura, K., Tombesi, F., Kazanas, D., et al. 2014, ApJ, 780, 120. doi:10.1088/0004-637X/780/2/120
- Ghisellini et al. (2004) Ghisellini, G., Haardt, F., & Matt, G. 2004, A&A, 413, 535. doi:10.1051/0004-6361:20031562
- Globus & Levinson (2016) Globus, N. & Levinson, A. 2016, MNRAS, 461, 2605. doi:10.1093/mnras/stw1474
- Gofford et al. (2015) Gofford, J., Reeves, J. N., McLaughlin, D. E., et al. 2015, MNRAS, 451, 4169. doi:10.1093/mnras/stv1207
- Gofford et al. (2013) Gofford, J., Reeves, J. N., Tombesi, F., et al. 2013, MNRAS, 430, 60. doi:10.1093/mnras/sts481
- Kaastra (2017) Kaastra, J. S. 2017, A&A, 605, A51. doi:10.1051/0004-6361/20162931910.48550/arXiv.1707.09202
- Kaastra & Bleeker (2016) Kaastra, J. S. & Bleeker, J. A. M. 2016, A&A, 587, A151. doi:10.1051/0004-6361/201527395
- Kaastra et al. (2014) Kaastra, J. S., Kriss, G. A., Cappi, M., et al. 2014, Science, 345, 64. doi:10.1126/science.1253787
- Kayanoki & Fukazawa (2022) Kayanoki, T. & Fukazawa, Y. 2022, PASJ, 74, 791. doi:10.1093/pasj/psac036
- Kellermann et al. (1989) Kellermann, K. I., Sramek, R., Schmidt, M., et al. 1989, AJ, 98, 1195. doi:10.1086/115207
- Kharb et al. (2014) Kharb, P., O’Dea, C. P., Baum, S. A., et al. 2014, MNRAS, 440, 2976. doi:10.1093/mnras/stu421
- Halpern (1984) Halpern, J. P. 1984, ApJ, 281, 90. doi:10.1086/162077
- HI4PI Collaboration (2016) HI4PI Collaboration, Ben Bekhti, N., Flöer, L., et al. 2016, A&A, 594, A116. doi:10.1051/0004-6361/201629178
- Laha et al. (2021) Laha, S., Reynolds, C. S., Reeves, J., et al. 2021, Nature Astronomy, 5, 13. doi:10.1038/s41550-020-01255-2
- Laha et al. (2016) Laha, S., Guainazzi, M., Chakravorty, S., et al. 2016, MNRAS, 457, 3896. doi:10.1093/mnras/stw211
- Lodders et al. (2009) Lodders, K., Palme, H., & Gail, H.-P. 2009, Landolt Börnstein, 4B, 712. doi:10.1007/978-3-540-88055-4_34
- Mao et al. (2022) Mao, J., Kriss, G. A., Landt, H., et al. 2022, ApJ, 940, 41. doi:10.3847/1538-4357/ac99de
- Mao et al. (2019) Mao, J., Mehdipour, M., Kaastra, J. S., et al. 2019, A&A, 621, A99. doi:10.1051/0004-6361/201833191
- Mao et al. (2017) Mao, J., Kaastra, J. S., Mehdipour, M., et al. 2017, A&A, 607, A100. doi:10.1051/0004-6361/201731378
- Mehdipour & Costantini (2019) Mehdipour, M. & Costantini, E. 2019, A&A, 625, A25. doi:10.1051/0004-6361/201935205
- Mehdipour et al. (2018) Mehdipour, M., Kaastra, J. S., Costantini, E., et al. 2018, A&A, 615, A72. doi:10.1051/0004-6361/201832604
- Mehdipour et al. (2017) Mehdipour, M., Kaastra, J. S., Kriss, G. A., et al. 2017, A&A, 607, A28. doi:10.1051/0004-6361/201731175
- Mehdipour et al. (2016) Mehdipour, M., Kaastra, J. S., & Kallman, T. 2016, A&A, 596, A65. doi:10.1051/0004-6361/201628721
- Mehdipour et al. (2010) Mehdipour, M., Branduardi-Raymont, G., & Page, M. J. 2010, A&A, 514, A100. doi:10.1051/0004-6361/200913049
- Miller et al. (2015) Miller, J. M., Kaastra, J. S., Miller, M. C., et al. 2015, Nature, 526, 542. doi:10.1038/nature15708
- Nomura et al. (2016) Nomura, M., Ohsuga, K., Takahashi, H. R., et al. 2016, PASJ, 68, 16. doi:10.1093/pasj/psv124
- Reeves et al. (2009) Reeves, J. N., Sambruna, R. M., Braito, V., et al. 2009, ApJ, 702, L187. doi:10.1088/0004-637X/702/2/L187
- Sako et al. (2001) Sako, M., Kahn, S. M., Behar, E., et al. 2001, A&A, 365, L168. doi:10.1051/0004-6361:20000081
- Schmidt & Green (1983) Schmidt, M. & Green, R. F. 1983, ApJ, 269, 352. doi:10.1086/161048
- Schurch et al. (2006) Schurch, N. J., Griffiths, R. E., & Warwick, R. S. 2006, MNRAS, 371, 211. doi:10.1111/j.1365-2966.2006.10581.x
- Serafinelli et al. (2019) Serafinelli, R., Tombesi, F., Vagnetti, F., et al. 2019, A&A, 627, A121. doi:10.1051/0004-6361/201935275
- Steenbrugge et al. (2005) Steenbrugge, K. C., Kaastra, J. S., Crenshaw, D. M., et al. 2005, A&A, 434, 569. doi:10.1051/0004-6361:20047138
- Strüder et al. (2001) Strüder, L., Briel, U., Dennerl, K., et al. 2001, A&A, 365, L18. doi:10.1051/0004-6361:20000066
- Tarter et al. (1969) Tarter, C. B., Tucker, W. H., & Salpeter, E. E. 1969, ApJ, 156, 943. doi:10.1086/150026
- Tashiro et al. (2020) Tashiro, M., Maejima, H., Toda, K., et al. 2020, Proc. SPIE, 11444, 1144422. doi:10.1117/12.2565812
- Tombesi et al. (2016) Tombesi, F., Reeves, J. N., Kallman, T., et al. 2016, ApJ, 830, 98. doi:10.3847/0004-637X/830/2/98
- Tombesi et al. (2014) Tombesi, F., Tazaki, F., Mushotzky, R. F., et al. 2014, MNRAS, 443, 2154. doi:10.1093/mnras/stu1297
- Tombesi et al. (2013) Tombesi, F., Reeves, J. N., Reynolds, C. S., et al. 2013, MNRAS, 434, 2707. doi:10.1093/mnras/stt1213
- Tombesi et al. (2012b) Tombesi, F., Sambruna, R. M., Marscher, A. P., et al. 2012, MNRAS, 424, 754. doi:10.1111/j.1365-2966.2012.21266.x
- Tombesi et al. (2012a) Tombesi, F., Cappi, M., Reeves, J. N., et al. 2012, MNRAS, 422, L1. doi:10.1111/j.1745-3933.2012.01221.x
- Tombesi et al. (2011) Tombesi, F., Cappi, M., Reeves, J. N., et al. 2011, ApJ, 742, 44. doi:10.1088/0004-637X/742/1/44
- Tombesi et al. (2010b) Tombesi, F., Cappi, M., Reeves, J. N., et al. 2010, A&A, 521, A57. doi:10.1051/0004-6361/200913440
- Tombesi et al. (2010a) Tombesi, F., Sambruna, R. M., Reeves, J. N., et al. 2010, ApJ, 719, 700. doi:10.1088/0004-637X/719/1/700
- Turner et al. (2001) Turner, M. J. L., Abbey, A., Arnaud, M., et al. 2001, A&A, 365, L27. doi:10.1051/0004-6361:20000087
- Torresi et al. (2010) Torresi, E., Grandi, P., Longinotti, A. L., et al. 2010, MNRAS, 401, L10. doi:10.1111/j.1745-3933.2009.00773.x
- Torresi et al. (2012) Torresi, E., Grandi, P., Costantini, E., et al. 2012, MNRAS, 419, 321. doi:10.1111/j.1365-2966.2011.19694.x
- Urry & Padovani (1995) Urry, C. M. & Padovani, P. 1995, PASP, 107, 803. doi:10.1086/133630
- Véron-Cetty & Véron (2010) Véron-Cetty, M.-P. & Véron, P. 2010, A&A, 518, A10. doi:10.1051/0004-6361/201014188
- XRISM Science Team (2022) XRISM Science Team 2022, arXiv:2202.05399. doi:10.48550/arXiv.2202.05399
- Wang et al. (2022) Wang, Y., Kaastra, J., Mehdipour, M., et al. 2022, A&A, 657, A77. doi:10.1051/0004-6361/202141599