License: arXiv.org perpetual non-exclusive license
arXiv:2312.11167v1 [astro-ph.HE] 18 Dec 2023

Ionized X-ray winds in the radio galaxy Mrk 6

Taishu Kayanoki,11{}^{1}start_FLOATSUPERSCRIPT 1 end_FLOATSUPERSCRIPT Junjie Mao,2,121{}^{2,1}start_FLOATSUPERSCRIPT 2 , 1 end_FLOATSUPERSCRIPT Yasushi Fukazawa11{}^{1}start_FLOATSUPERSCRIPT 1 end_FLOATSUPERSCRIPT
11{}^{1}start_FLOATSUPERSCRIPT 1 end_FLOATSUPERSCRIPTDepartment of Physics, Graduate School of Advanced Science and Engineering, Hiroshima University 1-3-1 Kagamiyama, Higashi-Hiroshima,
Hiroshima 739-8526, Japan
22{}^{2}start_FLOATSUPERSCRIPT 2 end_FLOATSUPERSCRIPTDepartment of Physics, Tsinghua University, Beijing 100084, China
E-mail: [email protected]
(Accepted XXX. Received YYY; in original form ZZZ)
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 34700200+400kms1subscriptsuperscript34700400200kmsuperscripts1-34700^{+400}_{-200}~{}{\rm km~{}s^{-1}}- 34700 start_POSTSUPERSCRIPT + 400 end_POSTSUPERSCRIPT start_POSTSUBSCRIPT - 200 end_POSTSUBSCRIPT roman_km roman_s start_POSTSUPERSCRIPT - 1 end_POSTSUPERSCRIPT) and a WA component with an outflow velocity of 7600±200kms1plus-or-minus7600200kmsuperscripts1-7600\pm 200~{}{\rm km~{}s}^{-1}- 7600 ± 200 roman_km roman_s start_POSTSUPERSCRIPT - 1 end_POSTSUPERSCRIPT. 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: galaxies
pubyear: 2023pagerange: Ionized X-ray winds in the radio galaxy Mrk 6Ionized X-ray winds in the radio galaxy Mrk 6

1 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 ξ𝜉\xiitalic_ξ, hydrogen column density NHsubscript𝑁HN_{\rm H}italic_N start_POSTSUBSCRIPT roman_H end_POSTSUBSCRIPT, 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 ξ𝜉\xiitalic_ξ is defined as ξ=Lion/nHr2𝜉subscript𝐿ionsubscript𝑛Hsuperscript𝑟2\xi=L_{\rm ion}/n_{\rm H}~{}r^{2}italic_ξ = italic_L start_POSTSUBSCRIPT roman_ion end_POSTSUBSCRIPT / italic_n start_POSTSUBSCRIPT roman_H end_POSTSUBSCRIPT italic_r start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT, where Lionsubscript𝐿ionL_{\rm ion}italic_L start_POSTSUBSCRIPT roman_ion end_POSTSUBSCRIPT (11000110001-10001 - 1000 Ryd) denotes the ionizing luminosity, nHsubscript𝑛Hn_{\rm H}italic_n start_POSTSUBSCRIPT roman_H end_POSTSUBSCRIPT represents the disk wind hydrogen number density, and r𝑟ritalic_r indicates the distance from the source (Tarter et al., 1969). Generally, UFOs are observed mainly as blueshifted absorption lines at 79797-97 - 9 keV (Tombesi et al., 2010b). The UFO ionization parameter is log[ξ(ergcms1)]=35𝜉ergcmsuperscripts135\log[\xi({\rm erg~{}cm~{}s^{-1}})]=3-5roman_log [ italic_ξ ( roman_erg roman_cm roman_s start_POSTSUPERSCRIPT - 1 end_POSTSUPERSCRIPT ) ] = 3 - 5, the hydrogen column density is log[NH(cm2)]=2223.5subscript𝑁Hsuperscriptcm22223.5\log[N_{\rm H}({\rm cm^{-2}})]=22-23.5roman_log [ italic_N start_POSTSUBSCRIPT roman_H end_POSTSUBSCRIPT ( roman_cm start_POSTSUPERSCRIPT - 2 end_POSTSUPERSCRIPT ) ] = 22 - 23.5, the outflow velocity is vout=10,00070,000kms1subscript𝑣out1000070000kmsuperscripts1v_{\rm out}=10,000-70,000~{}{\rm km~{}s^{-1}}italic_v start_POSTSUBSCRIPT roman_out end_POSTSUBSCRIPT = 10 , 000 - 70 , 000 roman_km roman_s start_POSTSUPERSCRIPT - 1 end_POSTSUPERSCRIPT (Laha et al., 2021), and the kinetic power of the mass outflow is 0.110%similar-toabsent0.1percent10\sim 0.1-10\%∼ 0.1 - 10 % of Eddington luminosity (e.g. Gofford et al., 2015; Tombesi et al., 2012a). As such, UFOs are estimated to exist in 40%similar-toabsentpercent40\sim 40\%∼ 40 % AGNs (Tombesi et al., 2010b, 2011; Gofford et al., 2013). Based on simulations (e.g. Nomura et al., 2016),UFOs are launched at 30similar-toabsent30\sim 30∼ 30 times the Schwarzschild radius RS=2GMBH/c2subscript𝑅S2𝐺subscript𝑀BHsuperscript𝑐2R_{\rm S}=2GM_{\rm BH}/c^{2}italic_R start_POSTSUBSCRIPT roman_S end_POSTSUBSCRIPT = 2 italic_G italic_M start_POSTSUBSCRIPT roman_BH end_POSTSUBSCRIPT / italic_c start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT, and from observations (e.g. Gofford et al., 2015), UFOs originate at 102104RSsuperscript102superscript104subscript𝑅S10^{2}-10^{4}~{}R_{\rm S}10 start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT - 10 start_POSTSUPERSCRIPT 4 end_POSTSUPERSCRIPT italic_R start_POSTSUBSCRIPT roman_S end_POSTSUBSCRIPT. 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 log[ξ(ergcms1)]=13𝜉ergcmsuperscripts113\log[\xi({\rm erg~{}cm~{}s^{-1}})]=-1-3roman_log [ italic_ξ ( roman_erg roman_cm roman_s start_POSTSUPERSCRIPT - 1 end_POSTSUPERSCRIPT ) ] = - 1 - 3, the hydrogen column density is log[NH(cm2)]=2122.5subscript𝑁Hsuperscriptcm22122.5\log[N_{\rm H}({\rm cm^{-2}})]=21-22.5roman_log [ italic_N start_POSTSUBSCRIPT roman_H end_POSTSUBSCRIPT ( roman_cm start_POSTSUPERSCRIPT - 2 end_POSTSUPERSCRIPT ) ] = 21 - 22.5, the outflow rate is vout=1002,000kms1subscript𝑣out1002000kmsuperscripts1v_{\rm out}=100-2,000~{}{\rm km~{}s^{-1}}italic_v start_POSTSUBSCRIPT roman_out end_POSTSUBSCRIPT = 100 - 2 , 000 roman_km roman_s start_POSTSUPERSCRIPT - 1 end_POSTSUPERSCRIPT (Sako et al., 2001; Laha et al., 2016). WAs exist at 0.110000.110000.1-10000.1 - 1000 pc, and their hydrogen number density nHsubscript𝑛Hn_{\rm H}italic_n start_POSTSUBSCRIPT roman_H end_POSTSUBSCRIPT has been estimated to be nH1041011cm3similar-tosubscript𝑛Hsuperscript104superscript1011superscriptcm3n_{\rm H}\sim 10^{4}-10^{11}~{}{\rm cm}^{-3}italic_n start_POSTSUBSCRIPT roman_H end_POSTSUBSCRIPT ∼ 10 start_POSTSUPERSCRIPT 4 end_POSTSUPERSCRIPT - 10 start_POSTSUPERSCRIPT 11 end_POSTSUPERSCRIPT roman_cm start_POSTSUPERSCRIPT - 3 end_POSTSUPERSCRIPT (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 R𝑅Ritalic_R is defined based on Kellermann et al. (1989) as R=F6cm/Fopt𝑅subscript𝐹6cmsubscript𝐹optR=F_{\rm 6~{}cm}/F_{\rm opt}italic_R = italic_F start_POSTSUBSCRIPT 6 roman_cm end_POSTSUBSCRIPT / italic_F start_POSTSUBSCRIPT roman_opt end_POSTSUBSCRIPT using the radio flux (F6cmsubscript𝐹6cmF_{\rm 6~{}cm}italic_F start_POSTSUBSCRIPT 6 roman_cm end_POSTSUBSCRIPT) and the optical flux logFopt=(48.36B)/2.5subscript𝐹opt48.36𝐵2.5\log F_{\rm opt}=(48.36-B)/2.5roman_log italic_F start_POSTSUBSCRIPT roman_opt end_POSTSUBSCRIPT = ( 48.36 - italic_B ) / 2.5 (Schmidt & Green, 1983). For radio-loud AGN, R10𝑅10R\geq 10italic_R ≥ 10, and for radio-quiet AGN, R<10𝑅10R<10italic_R < 10. 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). B𝐵Bitalic_B-band magnitude and 6 cm radio flux (F6cmsubscript𝐹6cmF_{\rm 6~{}cm}italic_F start_POSTSUBSCRIPT 6 roman_cm end_POSTSUBSCRIPT) of Mrk 6 is B=15.16𝐵15.16B=15.16italic_B = 15.16 and F6cm=0.115subscript𝐹6cm0.115F_{\rm 6~{}cm}=0.115italic_F start_POSTSUBSCRIPT 6 roman_cm end_POSTSUBSCRIPT = 0.115 Jy (Véron-Cetty & Véron, 2010), Mrk 6 can be classified as a radio-loud galaxy (R=30𝑅30R=30italic_R = 30). Specifically, it has a black hole mass of (1.8±0.2)×108Mplus-or-minus1.80.2superscript108subscript𝑀direct-product(1.8\pm 0.2)\times 10^{8}M_{\odot}( 1.8 ± 0.2 ) × 10 start_POSTSUPERSCRIPT 8 end_POSTSUPERSCRIPT italic_M start_POSTSUBSCRIPT ⊙ end_POSTSUBSCRIPT (Doroshenko et al., 2012) with z=0.019𝑧0.019z=0.019italic_z = 0.019 (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 (H0subscript𝐻0H_{0}italic_H start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPTΩmsubscriptΩ𝑚\Omega_{m}roman_Ω start_POSTSUBSCRIPT italic_m end_POSTSUBSCRIPTΩΛsubscriptΩΛ\Omega_{\Lambda}roman_Ω start_POSTSUBSCRIPT roman_Λ end_POSTSUBSCRIPT) used in this study are H0=70kms1Mpc1subscript𝐻070kmsuperscripts1superscriptMpc1H_{0}=70~{}\rm{km~{}s^{-1}Mpc^{-1}}italic_H start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT = 70 roman_km roman_s start_POSTSUPERSCRIPT - 1 end_POSTSUPERSCRIPT roman_Mpc start_POSTSUPERSCRIPT - 1 end_POSTSUPERSCRIPT, ΩM=0.3subscriptΩ𝑀0.3\Omega_{M}=0.3roman_Ω start_POSTSUBSCRIPT italic_M end_POSTSUBSCRIPT = 0.3, ΩΛ=0.7subscriptΩΛ0.7\Omega_{\Lambda}=0.7roman_Ω start_POSTSUBSCRIPT roman_Λ end_POSTSUBSCRIPT = 0.7.

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.

Table 1: XMM-Newton observation of Mrk 6.
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 (>0.50countss1absent0.50countssuperscripts1>0.50~{}{\rm counts~{}s^{-1}}> 0.50 roman_counts roman_s start_POSTSUPERSCRIPT - 1 end_POSTSUPERSCRIPT), 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{}^{\circ}start_FLOATSUPERSCRIPT ∘ end_FLOATSUPERSCRIPT 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 >3σabsent3𝜎>3\sigma> 3 italic_σ 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 1.21.31.21.31.2-1.31.2 - 1.3 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 1.7710.01.7710.01.77-10.01.77 - 10.0 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 7.027.07.027.07.0-27.07.0 - 27.0 Å 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 NH=7.63×1020cm2subscript𝑁H7.63superscript1020superscriptcm2N_{\rm H}=7.63\times 10^{20}~{}{\rm cm}^{-2}italic_N start_POSTSUBSCRIPT roman_H end_POSTSUBSCRIPT = 7.63 × 10 start_POSTSUPERSCRIPT 20 end_POSTSUPERSCRIPT roman_cm start_POSTSUPERSCRIPT - 2 end_POSTSUPERSCRIPT (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 5105105-105 - 10 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 5105105-105 - 10 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 7similar-toabsent7\sim 7∼ 7 keV were observed in Figure 1. In addition to Table 2, the model with the absorption lines displayed a more accurate C-statistic (ΔC=5.9Δ𝐶5.9\Delta C=-5.9roman_Δ italic_C = - 5.9, 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 7similar-toabsent7\sim 7∼ 7 keV absorption lines were considered to be caused by a highly ionized gas.

Table 2: Emission line and absorption line search
power-law power-law power-law
+ emission line + emission line
+ absorption line
pow ΓΓ\Gammaroman_Γ (1)1{}^{(1)}start_FLOATSUPERSCRIPT ( 1 ) end_FLOATSUPERSCRIPT 1.61±0.05plus-or-minus1.610.051.61\pm 0.051.61 ± 0.05 1.59±0.05plus-or-minus1.590.051.59\pm 0.051.59 ± 0.05 1.580.02+0.05subscriptsuperscript1.580.050.021.58^{+0.05}_{-0.02}1.58 start_POSTSUPERSCRIPT + 0.05 end_POSTSUPERSCRIPT start_POSTSUBSCRIPT - 0.02 end_POSTSUBSCRIPT
pow Norm (2)2{}^{(2)}start_FLOATSUPERSCRIPT ( 2 ) end_FLOATSUPERSCRIPT 2.90.2+0.3subscriptsuperscript2.90.30.22.9^{+0.3}_{-0.2}2.9 start_POSTSUPERSCRIPT + 0.3 end_POSTSUPERSCRIPT start_POSTSUBSCRIPT - 0.2 end_POSTSUBSCRIPT 2.70.2+0.3subscriptsuperscript2.70.30.22.7^{+0.3}_{-0.2}2.7 start_POSTSUPERSCRIPT + 0.3 end_POSTSUPERSCRIPT start_POSTSUBSCRIPT - 0.2 end_POSTSUBSCRIPT 2.70.2+0.1subscriptsuperscript2.70.10.22.7^{+0.1}_{-0.2}2.7 start_POSTSUPERSCRIPT + 0.1 end_POSTSUPERSCRIPT start_POSTSUBSCRIPT - 0.2 end_POSTSUBSCRIPT
gaus e (3)3{}^{(3)}start_FLOATSUPERSCRIPT ( 3 ) end_FLOATSUPERSCRIPT - 6.43±0.02plus-or-minus6.430.026.43\pm 0.026.43 ± 0.02 6.43±0.02plus-or-minus6.430.026.43\pm 0.026.43 ± 0.02
gaus norm (4)4{}^{(4)}start_FLOATSUPERSCRIPT ( 4 ) end_FLOATSUPERSCRIPT - 13±2plus-or-minus13213\pm 213 ± 2 12±2plus-or-minus12212\pm 212 ± 2
gaus FWHM (5)5{}^{(5)}start_FLOATSUPERSCRIPT ( 5 ) end_FLOATSUPERSCRIPT - 130±10plus-or-minus13010130\pm 10130 ± 10 120±10plus-or-minus12010120\pm 10120 ± 10
gaus e (3)3{}^{(3)}start_FLOATSUPERSCRIPT ( 3 ) end_FLOATSUPERSCRIPT - - 7.22±0.02plus-or-minus7.220.027.22\pm 0.027.22 ± 0.02
gaus norm (4)4{}^{(4)}start_FLOATSUPERSCRIPT ( 4 ) end_FLOATSUPERSCRIPT - - 3±1plus-or-minus31-3\pm 1- 3 ± 1
gaus FWHM (5)5{}^{(5)}start_FLOATSUPERSCRIPT ( 5 ) end_FLOATSUPERSCRIPT - - <160absent160<160< 160
Total C-stat / d.o.f. 149.9 / 83 82.4 / 80 76.5 / 77
Total Expected C-stat 85.1±plus-or-minus\pm±13.1 85.1±plus-or-minus\pm±13.1 85.1±plus-or-minus\pm±13.1
ΔΔ\Deltaroman_ΔC-stat - -67.5 -73.4
  • Fit the EPIC/pn hard X-ray data (5-10 keV). ΔΔ\Deltaroman_Δ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 1051photonss1keV1superscript1051photonssuperscripts1superscriptkeV110^{51}~{}{\rm photons~{}s^{-1}~{}keV^{-1}}10 start_POSTSUPERSCRIPT 51 end_POSTSUPERSCRIPT roman_photons roman_s start_POSTSUPERSCRIPT - 1 end_POSTSUPERSCRIPT roman_keV start_POSTSUPERSCRIPT - 1 end_POSTSUPERSCRIPT at 1 eV. (3) Line energy in keV. (4) Normalization of Gaussian in 1048photonss1superscript1048photonssuperscripts110^{48}~{}{\rm photons~{}s^{-1}}10 start_POSTSUPERSCRIPT 48 end_POSTSUPERSCRIPT roman_photons roman_s start_POSTSUPERSCRIPT - 1 end_POSTSUPERSCRIPT. (5) Line FWHM in keV.

Refer to caption
Figure 1: EPIC/pn spectra in 5105105-105 - 10 keV energy band and the power-law + emission line + absorption line model; black cross point and red line represent EPIC/pn spectra and power-law + gaussian (emission line) + gaussian (absorption line), respectively. For the EPIC/pn spectra we used the obin command for optimal binning.

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 6.46.46.46.4 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 (ΔΔ\Deltaroman_Δ 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 6.46.46.46.4 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 7similar-toabsent7\sim 7∼ 7 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 ξ𝜉\xiitalic_ξ-values by an external code, for instance, Cloudy or XSTAR. For the pion model, the hydrogen density were set to 1 cm3superscriptcm3{\rm cm}^{-3}roman_cm start_POSTSUPERSCRIPT - 3 end_POSTSUPERSCRIPT and the broadening velocity were set to 100 kms1kmsuperscripts1{\rm km~{}s^{-1}}roman_km roman_s start_POSTSUPERSCRIPT - 1 end_POSTSUPERSCRIPT. 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 7similar-toabsent7\sim~{}7∼ 7 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 7similar-toabsent7\sim~{}7∼ 7 keV absorption lines (however, the 7similar-toabsent7\sim~{}7∼ 7 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 7similar-toabsent7\sim 7∼ 7 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 7similar-toabsent7\sim 7∼ 7 keV absorption lines (Figure 5), its C-statistic of 36.6 was inferior to that obtained with the M3 model. However, the 7similar-toabsent7\sim 7∼ 7 keV absorption lines can be possibly attributed to the Fe xxv/xxvi absorption lines, the spectra obtained herein suggested that the 7similar-toabsent7\sim 7∼ 7 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 7similar-toabsent7\sim 7∼ 7 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 (CfsubscriptCf{\rm C_{f}}roman_C start_POSTSUBSCRIPT roman_f end_POSTSUBSCRIPT) of the UFO component for the M3 model and the M4 model. For the M3 model, the covering factor changed to CfsubscriptCf{\rm C_{f}}roman_C start_POSTSUBSCRIPT roman_f end_POSTSUBSCRIPT = 0.7 and the C-stat improved significantly (ΔC=4.3Δ𝐶4.3\Delta C=-4.3roman_Δ italic_C = - 4.3). 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 (ΔC=4.7Δ𝐶4.7\Delta C=-4.7roman_Δ italic_C = - 4.7). 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 00.35c00.35𝑐0-0.35~{}c0 - 0.35 italic_c 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
Refer to caption
Refer to caption
Refer to caption
Figure 2: Contour maps based on the M3 model (upper left panel: the ionization parameter vs. column density of UFO component. upper right panel: ionization parameter vs. outflow velocity of UFO component. bottom left panel: the ionization parameter vs. column density of WA component. bottom right panel: ionization parameter vs. outflow velocity of WA component.). Black cross point, black line, red line, green line, blue line, and cyan line show the best-fit parameter, 68.3% confidence level, 90% confidence level, 95.4% confidence level, 99% confidence level, and 99.99% confidence level, respectively.

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 χ2superscript𝜒2\chi^{2}italic_χ start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT distribution. Since the difference in degrees of freedom between the M2 model and the M3 model is 3, the p𝑝pitalic_p-value of the M3 model should be a p<0.001𝑝0.001p<0.001italic_p < 0.001. From the above, the confidence level of the M3 model is >99.9%absentpercent99.9>99.9\%> 99.9 %.

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 logξ(ergcms1)=2.10.1+0.5𝜉ergcmsuperscripts1subscriptsuperscript2.10.50.1\log\xi~{}({\rm erg~{}cm~{}s}^{-1})=2.1^{+0.5}_{-0.1}roman_log italic_ξ ( roman_erg roman_cm roman_s start_POSTSUPERSCRIPT - 1 end_POSTSUPERSCRIPT ) = 2.1 start_POSTSUPERSCRIPT + 0.5 end_POSTSUPERSCRIPT start_POSTSUBSCRIPT - 0.1 end_POSTSUBSCRIPT, high column density of NH=3.50.9+1.6×1022cm2subscript𝑁Hsubscriptsuperscript3.51.60.9superscript1022superscriptcm2N_{\rm H}=3.5^{+1.6}_{-0.9}\times 10^{22}~{}{\rm cm}^{-2}italic_N start_POSTSUBSCRIPT roman_H end_POSTSUBSCRIPT = 3.5 start_POSTSUPERSCRIPT + 1.6 end_POSTSUPERSCRIPT start_POSTSUBSCRIPT - 0.9 end_POSTSUBSCRIPT × 10 start_POSTSUPERSCRIPT 22 end_POSTSUPERSCRIPT roman_cm start_POSTSUPERSCRIPT - 2 end_POSTSUPERSCRIPT, and an extremely fast outflow velocity of vout=34700200+400kms1subscript𝑣outsubscriptsuperscript34700400200kmsuperscripts1v_{\rm out}=-34700^{+400}_{-200}~{}{\rm km~{}s}^{-1}italic_v start_POSTSUBSCRIPT roman_out end_POSTSUBSCRIPT = - 34700 start_POSTSUPERSCRIPT + 400 end_POSTSUPERSCRIPT start_POSTSUBSCRIPT - 200 end_POSTSUBSCRIPT roman_km roman_s start_POSTSUPERSCRIPT - 1 end_POSTSUPERSCRIPT. Thus, this ionized gas component was identified to be a UFO. The other ionized gas has a lower ionization parameter of logξ(ergcms1)=0.510.15+0.19𝜉ergcmsuperscripts1subscriptsuperscript0.510.190.15\log\xi~{}({\rm erg~{}cm~{}s}^{-1})=0.51^{+0.19}_{-0.15}roman_log italic_ξ ( roman_erg roman_cm roman_s start_POSTSUPERSCRIPT - 1 end_POSTSUPERSCRIPT ) = 0.51 start_POSTSUPERSCRIPT + 0.19 end_POSTSUPERSCRIPT start_POSTSUBSCRIPT - 0.15 end_POSTSUBSCRIPT, lower column density of NH=8.10.6+0.8×1021cm2subscript𝑁Hsubscriptsuperscript8.10.80.6superscript1021superscriptcm2N_{\rm H}=8.1^{+0.8}_{-0.6}\times 10^{21}~{}{\rm cm}^{-2}italic_N start_POSTSUBSCRIPT roman_H end_POSTSUBSCRIPT = 8.1 start_POSTSUPERSCRIPT + 0.8 end_POSTSUPERSCRIPT start_POSTSUBSCRIPT - 0.6 end_POSTSUBSCRIPT × 10 start_POSTSUPERSCRIPT 21 end_POSTSUPERSCRIPT roman_cm start_POSTSUPERSCRIPT - 2 end_POSTSUPERSCRIPT and a lower outflow velocity of vout=7600±200kms1subscript𝑣outplus-or-minus7600200kmsuperscripts1v_{\rm out}=-7600\pm 200~{}{\rm km~{}s}^{-1}italic_v start_POSTSUBSCRIPT roman_out end_POSTSUBSCRIPT = - 7600 ± 200 roman_km roman_s start_POSTSUPERSCRIPT - 1 end_POSTSUPERSCRIPT. 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 EW>2eV𝐸𝑊2eVEW>2~{}{\rm eV}italic_E italic_W > 2 roman_eV and soft X-ray WA absorption lines with EW>2eV𝐸𝑊2eVEW>2~{}{\rm eV}italic_E italic_W > 2 roman_eVare illustrated in Figure 7 and Figure 8, respectively. In addition, the information on UFO absorption is likely better identified with presented in Table 4.

Table 3: Absorber Search.
M1 M2 M3 M4
pow ΓΓ\Gammaroman_Γ (1)1{}^{(1)}start_FLOATSUPERSCRIPT ( 1 ) end_FLOATSUPERSCRIPT 1.35±0.02plus-or-minus1.350.021.35\pm 0.021.35 ± 0.02 1.48±0.03plus-or-minus1.480.031.48\pm 0.031.48 ± 0.03 1.570.04+0.03subscriptsuperscript1.570.030.041.57^{+0.03}_{-0.04}1.57 start_POSTSUPERSCRIPT + 0.03 end_POSTSUPERSCRIPT start_POSTSUBSCRIPT - 0.04 end_POSTSUBSCRIPT 1.38±0.02plus-or-minus1.380.021.38\pm 0.021.38 ± 0.02
pow Norm (2)2{}^{(2)}start_FLOATSUPERSCRIPT ( 2 ) end_FLOATSUPERSCRIPT 1.77±0.05plus-or-minus1.770.051.77\pm 0.051.77 ± 0.05 2.4±0.1plus-or-minus2.40.12.4\pm 0.12.4 ± 0.1 2.80.2+0.1subscriptsuperscript2.80.10.22.8^{+0.1}_{-0.2}2.8 start_POSTSUPERSCRIPT + 0.1 end_POSTSUPERSCRIPT start_POSTSUBSCRIPT - 0.2 end_POSTSUBSCRIPT 3.31.0+2.1subscriptsuperscript3.32.11.03.3^{+2.1}_{-1.0}3.3 start_POSTSUPERSCRIPT + 2.1 end_POSTSUPERSCRIPT start_POSTSUBSCRIPT - 1.0 end_POSTSUBSCRIPT
refl Scal (3)3{}^{(3)}start_FLOATSUPERSCRIPT ( 3 ) end_FLOATSUPERSCRIPT 0.32±0.05plus-or-minus0.320.050.32\pm 0.050.32 ± 0.05 0.36±0.06plus-or-minus0.360.060.36\pm 0.060.36 ± 0.06 0.380.05+0.06subscriptsuperscript0.380.060.050.38^{+0.06}_{-0.05}0.38 start_POSTSUPERSCRIPT + 0.06 end_POSTSUPERSCRIPT start_POSTSUBSCRIPT - 0.05 end_POSTSUBSCRIPT 0.170.05+0.04subscriptsuperscript0.170.040.050.17^{+0.04}_{-0.05}0.17 start_POSTSUPERSCRIPT + 0.04 end_POSTSUPERSCRIPT start_POSTSUBSCRIPT - 0.05 end_POSTSUBSCRIPT
hot NHsubscript𝑁HN_{\rm H}italic_N start_POSTSUBSCRIPT roman_H end_POSTSUBSCRIPT (4)4{}^{(4)}start_FLOATSUPERSCRIPT ( 4 ) end_FLOATSUPERSCRIPT 0.700.03+0.04subscriptsuperscript0.700.040.030.70^{+0.04}_{-0.03}0.70 start_POSTSUPERSCRIPT + 0.04 end_POSTSUPERSCRIPT start_POSTSUBSCRIPT - 0.03 end_POSTSUBSCRIPT 0.460.03+0.04subscriptsuperscript0.460.040.030.46^{+0.04}_{-0.03}0.46 start_POSTSUPERSCRIPT + 0.04 end_POSTSUPERSCRIPT start_POSTSUBSCRIPT - 0.03 end_POSTSUBSCRIPT - -
mbb Norm (5)5{}^{(5)}start_FLOATSUPERSCRIPT ( 5 ) end_FLOATSUPERSCRIPT 83+4×109subscriptsuperscript843superscript1098^{+4}_{-3}\times 10^{9}8 start_POSTSUPERSCRIPT + 4 end_POSTSUPERSCRIPT start_POSTSUBSCRIPT - 3 end_POSTSUBSCRIPT × 10 start_POSTSUPERSCRIPT 9 end_POSTSUPERSCRIPT 83+4×108subscriptsuperscript843superscript1088^{+4}_{-3}\times 10^{8}8 start_POSTSUPERSCRIPT + 4 end_POSTSUPERSCRIPT start_POSTSUBSCRIPT - 3 end_POSTSUBSCRIPT × 10 start_POSTSUPERSCRIPT 8 end_POSTSUPERSCRIPT <2×108absent2superscript108<2\times 10^{8}< 2 × 10 start_POSTSUPERSCRIPT 8 end_POSTSUPERSCRIPT 21+3×106subscriptsuperscript231superscript1062^{+3}_{-1}\times 10^{6}2 start_POSTSUPERSCRIPT + 3 end_POSTSUPERSCRIPT start_POSTSUBSCRIPT - 1 end_POSTSUBSCRIPT × 10 start_POSTSUPERSCRIPT 6 end_POSTSUPERSCRIPT
mbb t (6)6{}^{(6)}start_FLOATSUPERSCRIPT ( 6 ) end_FLOATSUPERSCRIPT 64±2plus-or-minus64264\pm 264 ± 2 72±3plus-or-minus72372\pm 372 ± 3 8010+20subscriptsuperscript80201080^{+20}_{-10}80 start_POSTSUPERSCRIPT + 20 end_POSTSUPERSCRIPT start_POSTSUBSCRIPT - 10 end_POSTSUBSCRIPT 14±30plus-or-minus143014\pm 3014 ± 30
pion logξ𝜉\log\xiroman_log italic_ξ (7)7{}^{(7)}start_FLOATSUPERSCRIPT ( 7 ) end_FLOATSUPERSCRIPT - 3.10.2+0.1subscriptsuperscript3.10.10.23.1^{+0.1}_{-0.2}3.1 start_POSTSUPERSCRIPT + 0.1 end_POSTSUPERSCRIPT start_POSTSUBSCRIPT - 0.2 end_POSTSUBSCRIPT 2.10.1+0.5subscriptsuperscript2.10.50.12.1^{+0.5}_{-0.1}2.1 start_POSTSUPERSCRIPT + 0.5 end_POSTSUPERSCRIPT start_POSTSUBSCRIPT - 0.1 end_POSTSUBSCRIPT 3.1±0.1plus-or-minus3.10.13.1\pm 0.13.1 ± 0.1
pion NHsubscript𝑁HN_{\rm H}italic_N start_POSTSUBSCRIPT roman_H end_POSTSUBSCRIPT (8)8{}^{(8)}start_FLOATSUPERSCRIPT ( 8 ) end_FLOATSUPERSCRIPT - 1.7±0.3plus-or-minus1.70.31.7\pm 0.31.7 ± 0.3 3.50.9+1.6subscriptsuperscript3.51.60.93.5^{+1.6}_{-0.9}3.5 start_POSTSUPERSCRIPT + 1.6 end_POSTSUPERSCRIPT start_POSTSUBSCRIPT - 0.9 end_POSTSUBSCRIPT 7040+60subscriptsuperscript70604070^{+60}_{-40}70 start_POSTSUPERSCRIPT + 60 end_POSTSUPERSCRIPT start_POSTSUBSCRIPT - 40 end_POSTSUBSCRIPT
pion voutsubscript𝑣outv_{\rm out}italic_v start_POSTSUBSCRIPT roman_out end_POSTSUBSCRIPT (9)9{}^{(9)}start_FLOATSUPERSCRIPT ( 9 ) end_FLOATSUPERSCRIPT - 21800300+200subscriptsuperscript21800200300-21800^{+200}_{-300}- 21800 start_POSTSUPERSCRIPT + 200 end_POSTSUPERSCRIPT start_POSTSUBSCRIPT - 300 end_POSTSUBSCRIPT 34700200+400subscriptsuperscript34700400200-34700^{+400}_{-200}- 34700 start_POSTSUPERSCRIPT + 400 end_POSTSUPERSCRIPT start_POSTSUBSCRIPT - 200 end_POSTSUBSCRIPT 12000300+200subscriptsuperscript12000200300-12000^{+200}_{-300}- 12000 start_POSTSUPERSCRIPT + 200 end_POSTSUPERSCRIPT start_POSTSUBSCRIPT - 300 end_POSTSUBSCRIPT
pion fcovsubscript𝑓covf_{\rm cov}italic_f start_POSTSUBSCRIPT roman_cov end_POSTSUBSCRIPT (10)10{}^{(10)}start_FLOATSUPERSCRIPT ( 10 ) end_FLOATSUPERSCRIPT - - 0.7±0.1plus-or-minus0.70.10.7\pm 0.10.7 ± 0.1 -
pion AFesubscript𝐴FeA_{\rm Fe}italic_A start_POSTSUBSCRIPT roman_Fe end_POSTSUBSCRIPT (11)11{}^{(11)}start_FLOATSUPERSCRIPT ( 11 ) end_FLOATSUPERSCRIPT - - - 41+4subscriptsuperscript4414^{+4}_{-1}4 start_POSTSUPERSCRIPT + 4 end_POSTSUPERSCRIPT start_POSTSUBSCRIPT - 1 end_POSTSUBSCRIPT
pion logξ𝜉\log\xiroman_log italic_ξ (7)7{}^{(7)}start_FLOATSUPERSCRIPT ( 7 ) end_FLOATSUPERSCRIPT - - 0.510.15+0.19subscriptsuperscript0.510.190.150.51^{+0.19}_{-0.15}0.51 start_POSTSUPERSCRIPT + 0.19 end_POSTSUPERSCRIPT start_POSTSUBSCRIPT - 0.15 end_POSTSUBSCRIPT 1.0±0.1plus-or-minus1.00.11.0\pm 0.11.0 ± 0.1
pion NHsubscript𝑁HN_{\rm H}italic_N start_POSTSUBSCRIPT roman_H end_POSTSUBSCRIPT (8)8{}^{(8)}start_FLOATSUPERSCRIPT ( 8 ) end_FLOATSUPERSCRIPT - - 0.810.06+0.08subscriptsuperscript0.810.080.060.81^{+0.08}_{-0.06}0.81 start_POSTSUPERSCRIPT + 0.08 end_POSTSUPERSCRIPT start_POSTSUBSCRIPT - 0.06 end_POSTSUBSCRIPT 1.1±0.1plus-or-minus1.10.11.1\pm 0.11.1 ± 0.1
pion voutsubscript𝑣outv_{\rm out}italic_v start_POSTSUBSCRIPT roman_out end_POSTSUBSCRIPT (9)9{}^{(9)}start_FLOATSUPERSCRIPT ( 9 ) end_FLOATSUPERSCRIPT - - 7600±200plus-or-minus7600200-7600\pm 200- 7600 ± 200 4300200+100subscriptsuperscript4300100200-4300^{+100}_{-200}- 4300 start_POSTSUPERSCRIPT + 100 end_POSTSUPERSCRIPT start_POSTSUBSCRIPT - 200 end_POSTSUBSCRIPT
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 1911.5±57.4plus-or-minus1911.557.41911.5\pm 57.41911.5 ± 57.4 1918.0±57.5plus-or-minus1918.057.51918.0\pm 57.51918.0 ± 57.5 1917.4±57.5plus-or-minus1917.457.51917.4\pm 57.51917.4 ± 57.5 1912.1±57.4plus-or-minus1912.157.41912.1\pm 57.41912.1 ± 57.4
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
ΔΔ\Deltaroman_ΔC-stat -- 68.068.0-68.0- 68.0 129.7129.7-129.7- 129.7 93.193.1-93.1- 93.1
  • 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. ΔΔ\Deltaroman_Δ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 1051photonss1keV1superscript1051photonssuperscripts1superscriptkeV110^{51}~{}{\rm photons~{}s^{-1}~{}keV^{-1}}10 start_POSTSUPERSCRIPT 51 end_POSTSUPERSCRIPT roman_photons roman_s start_POSTSUPERSCRIPT - 1 end_POSTSUPERSCRIPT roman_keV start_POSTSUPERSCRIPT - 1 end_POSTSUPERSCRIPT at 1 keV. (3) Scale parameter of the reflection component. (4) Column density NHN{{}_{H}}italic_N start_FLOATSUBSCRIPT roman_H end_FLOATSUBSCRIPT of neutral ISM gas component in the host galaxy of the AGN in 1022superscript102210^{22}10 start_POSTSUPERSCRIPT 22 end_POSTSUPERSCRIPT cm22{}^{-2}start_FLOATSUPERSCRIPT - 2 end_FLOATSUPERSCRIPT. (5) Normalization of the modified blackbody in 1034cm0.5superscript1034superscriptcm0.510^{34}~{}{\rm cm^{0.5}}10 start_POSTSUPERSCRIPT 34 end_POSTSUPERSCRIPT roman_cm start_POSTSUPERSCRIPT 0.5 end_POSTSUPERSCRIPT. (6) Temperature of the modified blackbody in eV. (7) Logarithm of the ionization parameter ξ𝜉\xiitalic_ξ of the ionized wind components in ergcms1ergcmsuperscripts1{\rm erg~{}cm~{}s^{-1}}roman_erg roman_cm roman_s start_POSTSUPERSCRIPT - 1 end_POSTSUPERSCRIPT. (8) Column density NHN{{}_{H}}italic_N start_FLOATSUBSCRIPT roman_H end_FLOATSUBSCRIPT of the ionized wind components in 1022superscript102210^{22}10 start_POSTSUPERSCRIPT 22 end_POSTSUPERSCRIPT cm22{}^{-2}start_FLOATSUPERSCRIPT - 2 end_FLOATSUPERSCRIPT. (9) Velocity of the outflowing ionized wind components in kms1kmsuperscripts1{\rm km~{}s^{-1}}roman_km roman_s start_POSTSUPERSCRIPT - 1 end_POSTSUPERSCRIPT. (10) Covering factor of the pion component. (11) Fe abundance of the pion component.

Refer to caption
Figure 3: All spectra (RGS + EPIC/pn) and model curve, wherein the black cross point, orange point, and red line represent EPIC/pn spectra, combined RGS spectra, and the M3 model, respectively. For the EPIC/pn spectra, we used the obin command for optimal binning. RGS spectra were combined by rgscombine command and the binning factor was set to 5.
Refer to caption
Figure 4: Soft X-ray spectra (RGS) and model curve, wherein the orange point, blue dotted line, magenta dash-dotted line, red line, and green line represent the combined RGS spectra, M1 model, M2 model, M3 model, and M4 model, respectively. RGS spectra were combined by rgscombine command and the binning factor was set to 5.
Refer to caption
Figure 5: Energy band spectra (EPIC/pn)of 5105105-105 - 10 keV and model curve, wherein the black cross point, blue-dotted line, magenta dash-dotted line, red line, and green line represent the EPIC/pn spectra, M1 model, M2 model, M3 model, and M4 model, respectively. For the EPIC/pn spectra, the obin command was used for optimal binning.
Refer to caption
Figure 6: SED curve of the model, wherein the black line, blue line, orange line, green line, and magenta line represent the M3 model SED, power-law model components, modified blackbody model, reflection model components, and the M4 model, respectively.
Refer to caption
Figure 7: Absorption lines of UFO, wherein the black cross point, red line, blue line, red dashed lines, and blue dashed lines represent the EPIC/pn spectra, M3 model, M4 model, UFO absorption lines (EW>2eV𝐸𝑊2eVEW>2~{}{\rm eV}italic_E italic_W > 2 roman_eV) calculated based on M3 model and UFO absorption lines (EW>2eV𝐸𝑊2eVEW>2~{}{\rm eV}italic_E italic_W > 2 roman_eV) calculated based on M3 model respectively.
Refer to caption
Figure 8: Absorption lines of WA, wherein the orange point, red line, blue line, red dashed lines, and blue dashed lines represent the combined RGS spectra, M3 model, M4 model, WA absorption lines (EW>2eV𝐸𝑊2eVEW>2~{}{\rm eV}italic_E italic_W > 2 roman_eV) calculated based on M3 model and WA absorption lines (EW>2eV𝐸𝑊2eVEW>2~{}{\rm eV}italic_E italic_W > 2 roman_eV) calculated based on M3 model, respectively. RGS spectra is combined by rgscombine command. For plotting, the binning factor of the combined RGS spectrum was set to 5.
Table 4: The UFO absorption lines from the M3 model and the M4 model (EW>2eV𝐸𝑊2eVEW>2~{}{\rm eV}italic_E italic_W > 2 roman_eV).
Absorption line Rest-frame Energy (keV) EW𝐸𝑊EWitalic_E italic_W (eV) Observed Energy (keV)
M3 model (vout=34700km/ssubscript𝑣out34700kmsv_{\rm out}=-34700~{}{\rm km/s}italic_v start_POSTSUBSCRIPT roman_out end_POSTSUBSCRIPT = - 34700 roman_km / roman_s)
Fe xix 6.4656 3.5 7.0794
Fe xviii 6.4347 3.4 7.0456
M4 model (vout=12000km/ssubscript𝑣out12000kmsv_{\rm out}=-12000~{}{\rm km/s}italic_v start_POSTSUBSCRIPT roman_out end_POSTSUBSCRIPT = - 12000 roman_km / roman_s)
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 (vout=34700km/ssubscript𝑣out34700kmsv_{\rm out}=-34700~{}{\rm km/s}italic_v start_POSTSUBSCRIPT roman_out end_POSTSUBSCRIPT = - 34700 roman_km / roman_s) and M4 model (vout=12000km/ssubscript𝑣out12000kmsv_{\rm out}=-12000~{}{\rm km/s}italic_v start_POSTSUBSCRIPT roman_out end_POSTSUBSCRIPT = - 12000 roman_km / roman_s) 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).

Table 5: Comparison of parameter of UFOs in 3C 390.3 and 4C +++74.26.
Mrk 6 3C 390.3*{}^{*}start_FLOATSUPERSCRIPT * end_FLOATSUPERSCRIPT 4C +++74.26*{}^{*}start_FLOATSUPERSCRIPT * end_FLOATSUPERSCRIPT
logξ𝜉\log\xiroman_log italic_ξ (1)1{}^{(1)}start_FLOATSUPERSCRIPT ( 1 ) end_FLOATSUPERSCRIPT 2.10.1+0.5subscriptsuperscript2.10.50.12.1^{+0.5}_{-0.1}2.1 start_POSTSUPERSCRIPT + 0.5 end_POSTSUPERSCRIPT start_POSTSUBSCRIPT - 0.1 end_POSTSUBSCRIPT 5.60.8+0.2subscriptsuperscript5.60.20.85.6^{+0.2}_{-0.8}5.6 start_POSTSUPERSCRIPT + 0.2 end_POSTSUPERSCRIPT start_POSTSUBSCRIPT - 0.8 end_POSTSUBSCRIPT 4.62±0.25plus-or-minus4.620.254.62\pm 0.254.62 ± 0.25
NHsubscript𝑁HN_{\rm H}italic_N start_POSTSUBSCRIPT roman_H end_POSTSUBSCRIPT (2)2{}^{(2)}start_FLOATSUPERSCRIPT ( 2 ) end_FLOATSUPERSCRIPT 3.50.9+1.6subscriptsuperscript3.51.60.93.5^{+1.6}_{-0.9}3.5 start_POSTSUPERSCRIPT + 1.6 end_POSTSUPERSCRIPT start_POSTSUBSCRIPT - 0.9 end_POSTSUBSCRIPT >3absent3>3> 3 >4absent4>4> 4
voutsubscript𝑣outv_{\rm out}italic_v start_POSTSUBSCRIPT roman_out end_POSTSUBSCRIPT (3)3{}^{(3)}start_FLOATSUPERSCRIPT ( 3 ) end_FLOATSUPERSCRIPT 34700200+400subscriptsuperscript34700400200-34700^{+400}_{-200}- 34700 start_POSTSUPERSCRIPT + 400 end_POSTSUPERSCRIPT start_POSTSUBSCRIPT - 200 end_POSTSUBSCRIPT 44000±1000plus-or-minus440001000-44000\pm 1000- 44000 ± 1000 13000±2000plus-or-minus130002000-13000\pm 2000- 13000 ± 2000
  • *: Tombesi et al. (2014). The xabs model was used for fitting. (1) Logarithm of the ionization parameter ξ𝜉\xiitalic_ξ of the ionized wind components in ergcms1ergcmsuperscripts1{\rm erg~{}cm~{}s^{-1}}roman_erg roman_cm roman_s start_POSTSUPERSCRIPT - 1 end_POSTSUPERSCRIPT. (2) Column density NHN{{}_{H}}italic_N start_FLOATSUBSCRIPT roman_H end_FLOATSUBSCRIPT of the ionized wind components in ×1022absentsuperscript1022\times 10^{22}~{}× 10 start_POSTSUPERSCRIPT 22 end_POSTSUPERSCRIPTcm22{}^{-2}start_FLOATSUPERSCRIPT - 2 end_FLOATSUPERSCRIPT. (3) Velocity of the outflowing ionized wind components kms1kmsuperscripts1\rm km~{}s^{-1}roman_km roman_s start_POSTSUPERSCRIPT - 1 end_POSTSUPERSCRIPT.

The ionization parameter is defined as ξ=Lion/nHr2𝜉subscript𝐿ionsubscript𝑛Hsuperscript𝑟2\xi=L_{\rm ion}/n_{\rm H}r^{2}italic_ξ = italic_L start_POSTSUBSCRIPT roman_ion end_POSTSUBSCRIPT / italic_n start_POSTSUBSCRIPT roman_H end_POSTSUBSCRIPT italic_r start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT (Tarter et al., 1969) using the hydrogen density nHsubscript𝑛Hn_{\rm H}italic_n start_POSTSUBSCRIPT roman_H end_POSTSUBSCRIPT, disk wind distance r𝑟ritalic_r, and luminosity of the light source Lionsubscript𝐿ionL_{\rm ion}italic_L start_POSTSUBSCRIPT roman_ion end_POSTSUBSCRIPT. Thus, the distance of the disk wind was calculated as

r=(LionnHξ)0.5𝑟superscriptsubscript𝐿ionsubscript𝑛H𝜉0.5\displaystyle r=\left({L_{\rm ion}\over n_{\rm H}\xi}\right)^{0.5}italic_r = ( divide start_ARG italic_L start_POSTSUBSCRIPT roman_ion end_POSTSUBSCRIPT end_ARG start_ARG italic_n start_POSTSUBSCRIPT roman_H end_POSTSUBSCRIPT italic_ξ end_ARG ) start_POSTSUPERSCRIPT 0.5 end_POSTSUPERSCRIPT (1)

Based on the M3 model result, the UFO ionization parameter was log10ξ=2.10.1+0.5ergcms1subscript10𝜉subscriptsuperscript2.10.50.1ergcmsuperscripts1\log_{10}\xi=2.1^{+0.5}_{-0.1}~{}{\rm erg~{}cm~{}s}^{-1}roman_log start_POSTSUBSCRIPT 10 end_POSTSUBSCRIPT italic_ξ = 2.1 start_POSTSUPERSCRIPT + 0.5 end_POSTSUPERSCRIPT start_POSTSUBSCRIPT - 0.1 end_POSTSUBSCRIPT roman_erg roman_cm roman_s start_POSTSUPERSCRIPT - 1 end_POSTSUPERSCRIPT and the non-filtered (before getting absorbed by the UFO pion component) ionized luminosity was calculated to Lion/ξ=3.2×1041cm1subscript𝐿ion𝜉3.2superscript1041superscriptcm1L_{\rm ion}/\xi=3.2\times 10^{41}~{}{\rm cm^{-1}}italic_L start_POSTSUBSCRIPT roman_ion end_POSTSUBSCRIPT / italic_ξ = 3.2 × 10 start_POSTSUPERSCRIPT 41 end_POSTSUPERSCRIPT roman_cm start_POSTSUPERSCRIPT - 1 end_POSTSUPERSCRIPT. Therefore, the distance of the UFO from the SMBH could be estimated as

r𝑟\displaystyle ritalic_r =\displaystyle== (3.2×1041102.1×nH)0.5superscript3.2superscript1041superscript102.1subscript𝑛H0.5\displaystyle\left({3.2\times 10^{41}\over 10^{2.1}\times n_{\rm H}}\right)^{0% .5}( divide start_ARG 3.2 × 10 start_POSTSUPERSCRIPT 41 end_POSTSUPERSCRIPT end_ARG start_ARG 10 start_POSTSUPERSCRIPT 2.1 end_POSTSUPERSCRIPT × italic_n start_POSTSUBSCRIPT roman_H end_POSTSUBSCRIPT end_ARG ) start_POSTSUPERSCRIPT 0.5 end_POSTSUPERSCRIPT (2)
=\displaystyle== 5.0×1015(nH108cm3)0.5cm5.0superscript1015superscriptsubscript𝑛Hsuperscript108superscriptcm30.5cm\displaystyle 5.0\times 10^{15}\left(\frac{n_{\rm H}}{10^{8}{\rm cm}^{-3}}% \right)^{-0.5}~{}{\rm cm}5.0 × 10 start_POSTSUPERSCRIPT 15 end_POSTSUPERSCRIPT ( divide start_ARG italic_n start_POSTSUBSCRIPT roman_H end_POSTSUBSCRIPT end_ARG start_ARG 10 start_POSTSUPERSCRIPT 8 end_POSTSUPERSCRIPT roman_cm start_POSTSUPERSCRIPT - 3 end_POSTSUPERSCRIPT end_ARG ) start_POSTSUPERSCRIPT - 0.5 end_POSTSUPERSCRIPT roman_cm

This is equivalent to 95(nH108cm3)0.5RS95superscriptsubscript𝑛Hsuperscript108superscriptcm30.5subscript𝑅S95\left(\frac{n_{\rm H}}{10^{8}{\rm cm}^{-3}}\right)^{-0.5}~{}R_{\rm S}95 ( divide start_ARG italic_n start_POSTSUBSCRIPT roman_H end_POSTSUBSCRIPT end_ARG start_ARG 10 start_POSTSUPERSCRIPT 8 end_POSTSUPERSCRIPT roman_cm start_POSTSUPERSCRIPT - 3 end_POSTSUPERSCRIPT end_ARG ) start_POSTSUPERSCRIPT - 0.5 end_POSTSUPERSCRIPT italic_R start_POSTSUBSCRIPT roman_S end_POSTSUBSCRIPT, where the Schwarzschild radius RS=2GMBH/c2=5.3×1013cmsubscript𝑅S2𝐺subscript𝑀BHsuperscript𝑐25.3superscript1013cmR_{\rm S}=2GM_{\rm BH}/c^{2}=5.3\times 10^{13}~{}{\rm cm}italic_R start_POSTSUBSCRIPT roman_S end_POSTSUBSCRIPT = 2 italic_G italic_M start_POSTSUBSCRIPT roman_BH end_POSTSUBSCRIPT / italic_c start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT = 5.3 × 10 start_POSTSUPERSCRIPT 13 end_POSTSUPERSCRIPT roman_cm. 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

rmin=2GMBHvout2subscript𝑟min2𝐺subscript𝑀BHsuperscriptsubscript𝑣out2\displaystyle r_{\rm min}={2GM_{\rm BH}\over v_{\rm out}^{2}}italic_r start_POSTSUBSCRIPT roman_min end_POSTSUBSCRIPT = divide start_ARG 2 italic_G italic_M start_POSTSUBSCRIPT roman_BH end_POSTSUBSCRIPT end_ARG start_ARG italic_v start_POSTSUBSCRIPT roman_out end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG (3)

using the gravitational constant G𝐺Gitalic_G, SMBH mass MBHsubscript𝑀BHM_{\rm BH}italic_M start_POSTSUBSCRIPT roman_BH end_POSTSUBSCRIPT, and outflow velocity voutsubscript𝑣outv_{\rm out}italic_v start_POSTSUBSCRIPT roman_out end_POSTSUBSCRIPT. Based on the analysis results, the outflow velocity of the UFO was vout=34700kms1subscript𝑣out34700kmsuperscripts1v_{\rm out}=34700~{}{\rm km~{}s^{-1}}italic_v start_POSTSUBSCRIPT roman_out end_POSTSUBSCRIPT = 34700 roman_km roman_s start_POSTSUPERSCRIPT - 1 end_POSTSUPERSCRIPT, and therefore, the lower limit of the UFO distance was calculated as rmin=4.0×1015cm(=75RS)subscript𝑟minannotated4.0superscript1015cmabsent75subscript𝑅Sr_{\rm min}=4.0\times 10^{15}~{}{\rm cm}~{}(=75~{}R_{\rm S})italic_r start_POSTSUBSCRIPT roman_min end_POSTSUBSCRIPT = 4.0 × 10 start_POSTSUPERSCRIPT 15 end_POSTSUPERSCRIPT roman_cm ( = 75 italic_R start_POSTSUBSCRIPT roman_S end_POSTSUBSCRIPT ). Assuming that the disk wind width ΔrΔ𝑟\Delta rroman_Δ italic_r did not exceed the distance r𝑟ritalic_r to the SMBH, the upper limit of the disk wind distance was computed as

rmax=LionCvξNHsubscript𝑟maxsubscript𝐿ionsubscript𝐶v𝜉subscript𝑁H\displaystyle r_{\rm max}={{L_{\rm ion}C_{\rm v}\over\xi N_{\rm H}}}italic_r start_POSTSUBSCRIPT roman_max end_POSTSUBSCRIPT = divide start_ARG italic_L start_POSTSUBSCRIPT roman_ion end_POSTSUBSCRIPT italic_C start_POSTSUBSCRIPT roman_v end_POSTSUBSCRIPT end_ARG start_ARG italic_ξ italic_N start_POSTSUBSCRIPT roman_H end_POSTSUBSCRIPT end_ARG (4)

considering the covering factor Cvsubscript𝐶vC_{\rm v}italic_C start_POSTSUBSCRIPT roman_v end_POSTSUBSCRIPT of the disk wind (Wang et al., 2022). M3 model assumed a covering factor of Cv=1subscript𝐶v1C_{\rm v}=1italic_C start_POSTSUBSCRIPT roman_v end_POSTSUBSCRIPT = 1, such that the upper limit of the UFO distance was rmax=7.2×1016cm(=1.3×103RS)subscript𝑟maxannotated7.2superscript1016cmabsent1.3superscript103subscript𝑅Sr_{\rm max}=7.2\times 10^{16}~{}{\rm cm}~{}(=1.3\times 10^{3}~{}R_{\rm S})italic_r start_POSTSUBSCRIPT roman_max end_POSTSUBSCRIPT = 7.2 × 10 start_POSTSUPERSCRIPT 16 end_POSTSUPERSCRIPT roman_cm ( = 1.3 × 10 start_POSTSUPERSCRIPT 3 end_POSTSUPERSCRIPT italic_R start_POSTSUBSCRIPT roman_S end_POSTSUBSCRIPT ), considering a column density of NH=3.5×1022cm2subscript𝑁H3.5superscript1022superscriptcm2N_{\rm H}=3.5\times 10^{22}~{}{\rm cm}^{-2}italic_N start_POSTSUBSCRIPT roman_H end_POSTSUBSCRIPT = 3.5 × 10 start_POSTSUPERSCRIPT 22 end_POSTSUPERSCRIPT roman_cm start_POSTSUPERSCRIPT - 2 end_POSTSUPERSCRIPT. Therefore, the UFOs of Mrk 6 originated at 95(nH108cm3)0.5RS95superscriptsubscript𝑛Hsuperscript108superscriptcm30.5subscript𝑅S95\left(\frac{n_{\rm H}}{10^{8}{\rm cm}^{-3}}\right)^{-0.5}~{}R_{\rm S}95 ( divide start_ARG italic_n start_POSTSUBSCRIPT roman_H end_POSTSUBSCRIPT end_ARG start_ARG 10 start_POSTSUPERSCRIPT 8 end_POSTSUPERSCRIPT roman_cm start_POSTSUPERSCRIPT - 3 end_POSTSUPERSCRIPT end_ARG ) start_POSTSUPERSCRIPT - 0.5 end_POSTSUPERSCRIPT italic_R start_POSTSUBSCRIPT roman_S end_POSTSUBSCRIPT, and the range of the UFOs is 7.5×1011.3×103RS7.5superscript1011.3superscript103subscript𝑅S7.5\times 10^{1}-1.3\times 10^{3}~{}R_{\rm S}7.5 × 10 start_POSTSUPERSCRIPT 1 end_POSTSUPERSCRIPT - 1.3 × 10 start_POSTSUPERSCRIPT 3 end_POSTSUPERSCRIPT italic_R start_POSTSUBSCRIPT roman_S end_POSTSUBSCRIPT. Based on observational studies, UFOs are considered to originate in the 102104RSsuperscript102superscript104subscript𝑅S10^{2}-10^{4}~{}R_{\rm S}10 start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT - 10 start_POSTSUPERSCRIPT 4 end_POSTSUPERSCRIPT italic_R start_POSTSUBSCRIPT roman_S end_POSTSUBSCRIPT (e.g. Gofford et al., 2015), and according to simulation research, they originate in the 30RSsimilar-toabsent30subscript𝑅S\sim 30~{}R_{\rm S}∼ 30 italic_R start_POSTSUBSCRIPT roman_S end_POSTSUBSCRIPT (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

nH=Lionξr2subscript𝑛Hsubscript𝐿ion𝜉superscript𝑟2\displaystyle n_{\rm H}={L_{\rm ion}\over\xi r^{2}}italic_n start_POSTSUBSCRIPT roman_H end_POSTSUBSCRIPT = divide start_ARG italic_L start_POSTSUBSCRIPT roman_ion end_POSTSUBSCRIPT end_ARG start_ARG italic_ξ italic_r start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG (5)

The range of the UFOs was 4.0×10157.2×1016cm4.0superscript10157.2superscript1016cm4.0\times 10^{15}-7.2\times 10^{16}~{}{\rm cm}4.0 × 10 start_POSTSUPERSCRIPT 15 end_POSTSUPERSCRIPT - 7.2 × 10 start_POSTSUPERSCRIPT 16 end_POSTSUPERSCRIPT roman_cm, and thus, the range of the hydrogen density was calculated as 4.8×1051.6×108cm34.8superscript1051.6superscript108superscriptcm34.8\times 10^{5}-1.6\times 10^{8}~{}{\rm cm}^{-3}4.8 × 10 start_POSTSUPERSCRIPT 5 end_POSTSUPERSCRIPT - 1.6 × 10 start_POSTSUPERSCRIPT 8 end_POSTSUPERSCRIPT roman_cm start_POSTSUPERSCRIPT - 3 end_POSTSUPERSCRIPT.

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).

Table 6: Comparison of parameters of WAs in the 3C 390.3 and 4C +++74.26.
Mrk 6 3C 390.3*{}^{*}start_FLOATSUPERSCRIPT * end_FLOATSUPERSCRIPT 4C +++74.26*{}^{*}start_FLOATSUPERSCRIPT * end_FLOATSUPERSCRIPT
logξ𝜉\log\xiroman_log italic_ξ (1)1{}^{(1)}start_FLOATSUPERSCRIPT ( 1 ) end_FLOATSUPERSCRIPT 0.510.15+0.19subscriptsuperscript0.510.190.150.51^{+0.19}_{-0.15}0.51 start_POSTSUPERSCRIPT + 0.19 end_POSTSUPERSCRIPT start_POSTSUBSCRIPT - 0.15 end_POSTSUBSCRIPT 1.63±0.11plus-or-minus1.630.111.63\pm 0.111.63 ± 0.11, 2.77±0.06plus-or-minus2.770.062.77\pm 0.062.77 ± 0.06 1.69±0.04plus-or-minus1.690.041.69\pm 0.041.69 ± 0.04, 2.46±0.04plus-or-minus2.460.042.46\pm 0.042.46 ± 0.04
NHsubscript𝑁HN_{\rm H}italic_N start_POSTSUBSCRIPT roman_H end_POSTSUBSCRIPT (2)2{}^{(2)}start_FLOATSUPERSCRIPT ( 2 ) end_FLOATSUPERSCRIPT 8.10.6+0.8subscriptsuperscript8.10.80.68.1^{+0.8}_{-0.6}8.1 start_POSTSUPERSCRIPT + 0.8 end_POSTSUPERSCRIPT start_POSTSUBSCRIPT - 0.6 end_POSTSUBSCRIPT 0.37±0.06plus-or-minus0.370.060.37\pm 0.060.37 ± 0.06, 1.1±0.5plus-or-minus1.10.51.1\pm 0.51.1 ± 0.5 3.6±0.6plus-or-minus3.60.63.6\pm 0.63.6 ± 0.6, 6.8±0.8plus-or-minus6.80.86.8\pm 0.86.8 ± 0.8
voutsubscript𝑣outv_{\rm out}italic_v start_POSTSUBSCRIPT roman_out end_POSTSUBSCRIPT (3)3{}^{(3)}start_FLOATSUPERSCRIPT ( 3 ) end_FLOATSUPERSCRIPT 7600±200plus-or-minus7600200-7600\pm 200- 7600 ± 200 1500±60plus-or-minus150060-1500\pm 60- 1500 ± 60, 50±100plus-or-minus5010050\pm 10050 ± 100 1490±90plus-or-minus149090-1490\pm 90- 1490 ± 90, 3000±500plus-or-minus3000500-3000\pm 500- 3000 ± 500
  • *: Mehdipour & Costantini (2019). The pion model was used for fitting. (1) Logarithm of the ionization parameter ξ𝜉\xiitalic_ξ of the ionized wind components in ergcms1ergcmsuperscripts1{\rm erg~{}cm~{}s^{-1}}roman_erg roman_cm roman_s start_POSTSUPERSCRIPT - 1 end_POSTSUPERSCRIPT. (2) Column density NHN{{}_{H}}italic_N start_FLOATSUBSCRIPT roman_H end_FLOATSUBSCRIPT of the ionized wind components in ×1022absentsuperscript1022\times 10^{22}~{}× 10 start_POSTSUPERSCRIPT 22 end_POSTSUPERSCRIPTcm22{}^{-2}start_FLOATSUPERSCRIPT - 2 end_FLOATSUPERSCRIPT. (3) Velocity of the outflowing ionized wind components kms1kmsuperscripts1\rm km~{}s^{-1}roman_km roman_s start_POSTSUPERSCRIPT - 1 end_POSTSUPERSCRIPT.

Based on the results of the analysis with the M3 model, the ionization parameter of WA was log10ξ(ergcms1)=0.510.15+0.19subscript10𝜉ergcmsuperscripts1subscriptsuperscript0.510.190.15\log_{10}\xi~{}({\rm erg~{}cm~{}s}^{-1})=0.51^{+0.19}_{-0.15}roman_log start_POSTSUBSCRIPT 10 end_POSTSUBSCRIPT italic_ξ ( roman_erg roman_cm roman_s start_POSTSUPERSCRIPT - 1 end_POSTSUPERSCRIPT ) = 0.51 start_POSTSUPERSCRIPT + 0.19 end_POSTSUPERSCRIPT start_POSTSUBSCRIPT - 0.15 end_POSTSUBSCRIPT 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 Lion/ξ=1.0×1043cm1subscript𝐿ion𝜉1.0superscript1043superscriptcm1L_{\rm ion}/\xi=1.0\times 10^{43}~{}{\rm cm^{-1}}italic_L start_POSTSUBSCRIPT roman_ion end_POSTSUBSCRIPT / italic_ξ = 1.0 × 10 start_POSTSUPERSCRIPT 43 end_POSTSUPERSCRIPT roman_cm start_POSTSUPERSCRIPT - 1 end_POSTSUPERSCRIPT. Therefore, with Equation (1), the distance of WA from SMBH was

r𝑟\displaystyle ritalic_r =\displaystyle== (1.0×1043100.5×nH)0.5superscript1.0superscript1043superscript100.5subscript𝑛H0.5\displaystyle\left({1.0\times 10^{43}\over 10^{0.5}\times n_{\rm H}}\right)^{0% .5}( divide start_ARG 1.0 × 10 start_POSTSUPERSCRIPT 43 end_POSTSUPERSCRIPT end_ARG start_ARG 10 start_POSTSUPERSCRIPT 0.5 end_POSTSUPERSCRIPT × italic_n start_POSTSUBSCRIPT roman_H end_POSTSUBSCRIPT end_ARG ) start_POSTSUPERSCRIPT 0.5 end_POSTSUPERSCRIPT (6)
=\displaystyle== 5.6×1018(nH105cm3)0.5cm5.6superscript1018superscriptsubscript𝑛Hsuperscript105superscriptcm30.5cm\displaystyle 5.6\times 10^{18}\left(\frac{n_{\rm H}}{10^{5}{\rm cm}^{-3}}% \right)^{-0.5}~{}{\rm cm}5.6 × 10 start_POSTSUPERSCRIPT 18 end_POSTSUPERSCRIPT ( divide start_ARG italic_n start_POSTSUBSCRIPT roman_H end_POSTSUBSCRIPT end_ARG start_ARG 10 start_POSTSUPERSCRIPT 5 end_POSTSUPERSCRIPT roman_cm start_POSTSUPERSCRIPT - 3 end_POSTSUPERSCRIPT end_ARG ) start_POSTSUPERSCRIPT - 0.5 end_POSTSUPERSCRIPT roman_cm

which is equivalent to 1.8(nH105cm3)0.5pc1.8superscriptsubscript𝑛Hsuperscript105superscriptcm30.5pc1.8\left(\frac{n_{\rm H}}{10^{5}{\rm cm}^{-3}}\right)^{-0.5}~{}{\rm pc}1.8 ( divide start_ARG italic_n start_POSTSUBSCRIPT roman_H end_POSTSUBSCRIPT end_ARG start_ARG 10 start_POSTSUPERSCRIPT 5 end_POSTSUPERSCRIPT roman_cm start_POSTSUPERSCRIPT - 3 end_POSTSUPERSCRIPT end_ARG ) start_POSTSUPERSCRIPT - 0.5 end_POSTSUPERSCRIPT roman_pc. The fitting results exhibited that the outflow velocity of the WA was vout=7600kms1subscript𝑣out7600kmsuperscripts1v_{\rm out}=7600~{}{\rm km~{}s^{-1}}italic_v start_POSTSUBSCRIPT roman_out end_POSTSUBSCRIPT = 7600 roman_km roman_s start_POSTSUPERSCRIPT - 1 end_POSTSUPERSCRIPT, and the lower limit of the WA distance was rmin=8.3×1016cm(=2.7×102pcr_{\rm min}=8.3\times 10^{16}~{}{\rm cm}~{}(=2.7\times 10^{-2}~{}{\rm pc}italic_r start_POSTSUBSCRIPT roman_min end_POSTSUBSCRIPT = 8.3 × 10 start_POSTSUPERSCRIPT 16 end_POSTSUPERSCRIPT roman_cm ( = 2.7 × 10 start_POSTSUPERSCRIPT - 2 end_POSTSUPERSCRIPT roman_pc) using Equation (3). As the column density of the WA was NH=8.10.6+0.8×1021cm3subscript𝑁Hsubscriptsuperscript8.10.80.6superscript1021superscriptcm3N_{\rm H}=8.1^{+0.8}_{-0.6}\times 10^{21}~{}{\rm cm}^{-3}italic_N start_POSTSUBSCRIPT roman_H end_POSTSUBSCRIPT = 8.1 start_POSTSUPERSCRIPT + 0.8 end_POSTSUPERSCRIPT start_POSTSUBSCRIPT - 0.6 end_POSTSUBSCRIPT × 10 start_POSTSUPERSCRIPT 21 end_POSTSUPERSCRIPT roman_cm start_POSTSUPERSCRIPT - 3 end_POSTSUPERSCRIPT according to the analysis and the coverage factor was assumed to be Cv=1subscript𝐶v1C_{\rm v}=1italic_C start_POSTSUBSCRIPT roman_v end_POSTSUBSCRIPT = 1, the upper limit of the WA distance was calculated to be rmax=3.9×1020cm(=1.3×102pc)subscript𝑟maxannotated3.9superscript1020cmabsent1.3superscript102pcr_{\rm max}=3.9\times 10^{20}~{}{\rm cm}~{}(=1.3\times 10^{2}~{}{\rm pc})italic_r start_POSTSUBSCRIPT roman_max end_POSTSUBSCRIPT = 3.9 × 10 start_POSTSUPERSCRIPT 20 end_POSTSUPERSCRIPT roman_cm ( = 1.3 × 10 start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT roman_pc ) using Equation (4). Therefore, the WAs of Mrk 6 originated at 1.8(nH105cm3)0.5pc1.8superscriptsubscript𝑛Hsuperscript105superscriptcm30.5pc1.8\left(\frac{n_{\rm H}}{10^{5}{\rm cm}^{-3}}\right)^{-0.5}~{}{\rm pc}1.8 ( divide start_ARG italic_n start_POSTSUBSCRIPT roman_H end_POSTSUBSCRIPT end_ARG start_ARG 10 start_POSTSUPERSCRIPT 5 end_POSTSUPERSCRIPT roman_cm start_POSTSUPERSCRIPT - 3 end_POSTSUPERSCRIPT end_ARG ) start_POSTSUPERSCRIPT - 0.5 end_POSTSUPERSCRIPT roman_pc, and the range of the WAs was 2.7×1021.3×102pc2.7superscript1021.3superscript102pc2.7\times 10^{-2}-1.3\times 10^{2}~{}{\rm pc}2.7 × 10 start_POSTSUPERSCRIPT - 2 end_POSTSUPERSCRIPT - 1.3 × 10 start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT roman_pc. WAs are considered to originate in the 0.110000.110000.1-10000.1 - 1000 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 8.3×10163.9×1020cm8.3superscript10163.9superscript1020cm8.3\times 10^{16}-3.9\times 10^{20}~{}{\rm cm}8.3 × 10 start_POSTSUPERSCRIPT 16 end_POSTSUPERSCRIPT - 3.9 × 10 start_POSTSUPERSCRIPT 20 end_POSTSUPERSCRIPT roman_cm, such that the hydrogen density of WA was calculated as 2.1×1014.6×108cm32.1superscript1014.6superscript108superscriptcm32.1\times 10^{1}-4.6\times 10^{8}~{}{\rm cm}^{-3}2.1 × 10 start_POSTSUPERSCRIPT 1 end_POSTSUPERSCRIPT - 4.6 × 10 start_POSTSUPERSCRIPT 8 end_POSTSUPERSCRIPT roman_cm start_POSTSUPERSCRIPT - 3 end_POSTSUPERSCRIPT using Equation (5). The hydrogen number density of WAs was estimated to be 1041011cm3similar-toabsentsuperscript104superscript1011superscriptcm3\sim 10^{4}-10^{11}~{}{\rm cm}^{-3}∼ 10 start_POSTSUPERSCRIPT 4 end_POSTSUPERSCRIPT - 10 start_POSTSUPERSCRIPT 11 end_POSTSUPERSCRIPT roman_cm start_POSTSUPERSCRIPT - 3 end_POSTSUPERSCRIPT (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 0.3410.00.3410.00.34-10.00.34 - 10.0 keV. For comparison, we simulated the M4 model to obtain results with Fe xxv/xxvi absorption lines at 7similar-toabsent7\sim 7∼ 7 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 NHsubscript𝑁HN_{\rm H}italic_N start_POSTSUBSCRIPT roman_H end_POSTSUBSCRIPT, which is a vital parameter for the significant detection of the disk wind, was reduced. The 7similar-toabsent7\sim 7∼ 7 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 7similar-toabsent7\sim 7∼ 7 keV absorption lines in Mrk 6 can be attributed to Fe xxv/xxvi or Fe xix/xviii.

Table 7: XRISM/Resolve simulation.
M3 M4
pow ΓΓ\Gammaroman_Γ (1)1{}^{(1)}start_FLOATSUPERSCRIPT ( 1 ) end_FLOATSUPERSCRIPT 1.56±0.02plus-or-minus1.560.021.56\pm 0.021.56 ± 0.02 1.40±0.01plus-or-minus1.400.011.40\pm 0.011.40 ± 0.01
pow Norm (2)2{}^{(2)}start_FLOATSUPERSCRIPT ( 2 ) end_FLOATSUPERSCRIPT 2.8±0.1plus-or-minus2.80.12.8\pm 0.12.8 ± 0.1 3.60.4+0.5subscriptsuperscript3.60.50.43.6^{+0.5}_{-0.4}3.6 start_POSTSUPERSCRIPT + 0.5 end_POSTSUPERSCRIPT start_POSTSUBSCRIPT - 0.4 end_POSTSUBSCRIPT
refl Scal (3)3{}^{(3)}start_FLOATSUPERSCRIPT ( 3 ) end_FLOATSUPERSCRIPT 0.37±0.02plus-or-minus0.370.020.37\pm 0.020.37 ± 0.02 0.16±0.02plus-or-minus0.160.020.16\pm 0.020.16 ± 0.02
mbb Norm (4)4{}^{(4)}start_FLOATSUPERSCRIPT ( 4 ) end_FLOATSUPERSCRIPT 1.10.4+0.6×108subscriptsuperscript1.10.60.4superscript1081.1^{+0.6}_{-0.4}\times 10^{8}1.1 start_POSTSUPERSCRIPT + 0.6 end_POSTSUPERSCRIPT start_POSTSUBSCRIPT - 0.4 end_POSTSUBSCRIPT × 10 start_POSTSUPERSCRIPT 8 end_POSTSUPERSCRIPT 3.10.7+1.0×106subscriptsuperscript3.11.00.7superscript1063.1^{+1.0}_{-0.7}\times 10^{6}3.1 start_POSTSUPERSCRIPT + 1.0 end_POSTSUPERSCRIPT start_POSTSUBSCRIPT - 0.7 end_POSTSUBSCRIPT × 10 start_POSTSUPERSCRIPT 6 end_POSTSUPERSCRIPT
mbb t (5)5{}^{(5)}start_FLOATSUPERSCRIPT ( 5 ) end_FLOATSUPERSCRIPT 826+8subscriptsuperscript828682^{+8}_{-6}82 start_POSTSUPERSCRIPT + 8 end_POSTSUPERSCRIPT start_POSTSUBSCRIPT - 6 end_POSTSUBSCRIPT 130±10plus-or-minus13010130\pm 10130 ± 10
pion logξ𝜉\log\xiroman_log italic_ξ (6)6{}^{(6)}start_FLOATSUPERSCRIPT ( 6 ) end_FLOATSUPERSCRIPT 2.090.03+0.04subscriptsuperscript2.090.040.032.09^{+0.04}_{-0.03}2.09 start_POSTSUPERSCRIPT + 0.04 end_POSTSUPERSCRIPT start_POSTSUBSCRIPT - 0.03 end_POSTSUBSCRIPT 3.10±0.04plus-or-minus3.100.043.10\pm 0.043.10 ± 0.04
pion NHsubscript𝑁HN_{\rm H}italic_N start_POSTSUBSCRIPT roman_H end_POSTSUBSCRIPT (7)7{}^{(7)}start_FLOATSUPERSCRIPT ( 7 ) end_FLOATSUPERSCRIPT 3.6±0.3plus-or-minus3.60.33.6\pm 0.33.6 ± 0.3 8010+20subscriptsuperscript80201080^{+20}_{-10}80 start_POSTSUPERSCRIPT + 20 end_POSTSUPERSCRIPT start_POSTSUBSCRIPT - 10 end_POSTSUBSCRIPT
pion voutsubscript𝑣outv_{\rm out}italic_v start_POSTSUBSCRIPT roman_out end_POSTSUBSCRIPT (8)8{}^{(8)}start_FLOATSUPERSCRIPT ( 8 ) end_FLOATSUPERSCRIPT 34650±40plus-or-minus3465040-34650\pm 40- 34650 ± 40 12040±10plus-or-minus1204010-12040\pm 10- 12040 ± 10
pion fcovsubscript𝑓covf_{\rm cov}italic_f start_POSTSUBSCRIPT roman_cov end_POSTSUBSCRIPT (9)9{}^{(9)}start_FLOATSUPERSCRIPT ( 9 ) end_FLOATSUPERSCRIPT 0.67±0.02plus-or-minus0.670.020.67\pm 0.020.67 ± 0.02 -
pion AFesubscript𝐴FeA_{\rm Fe}italic_A start_POSTSUBSCRIPT roman_Fe end_POSTSUBSCRIPT (10)10{}^{(10)}start_FLOATSUPERSCRIPT ( 10 ) end_FLOATSUPERSCRIPT - 3.30.4+0.5subscriptsuperscript3.30.50.43.3^{+0.5}_{-0.4}3.3 start_POSTSUPERSCRIPT + 0.5 end_POSTSUPERSCRIPT start_POSTSUBSCRIPT - 0.4 end_POSTSUBSCRIPT
pion logξ𝜉\log\xiroman_log italic_ξ (6)6{}^{(6)}start_FLOATSUPERSCRIPT ( 6 ) end_FLOATSUPERSCRIPT 0.460.05+0.06subscriptsuperscript0.460.060.050.46^{+0.06}_{-0.05}0.46 start_POSTSUPERSCRIPT + 0.06 end_POSTSUPERSCRIPT start_POSTSUBSCRIPT - 0.05 end_POSTSUBSCRIPT 0.97±0.03plus-or-minus0.970.030.97\pm 0.030.97 ± 0.03
pion NHsubscript𝑁HN_{\rm H}italic_N start_POSTSUBSCRIPT roman_H end_POSTSUBSCRIPT (7)7{}^{(7)}start_FLOATSUPERSCRIPT ( 7 ) end_FLOATSUPERSCRIPT 0.80±0.02plus-or-minus0.800.020.80\pm 0.020.80 ± 0.02 1.130.02+0.03subscriptsuperscript1.130.030.021.13^{+0.03}_{-0.02}1.13 start_POSTSUPERSCRIPT + 0.03 end_POSTSUPERSCRIPT start_POSTSUBSCRIPT - 0.02 end_POSTSUBSCRIPT
pion voutsubscript𝑣outv_{\rm out}italic_v start_POSTSUBSCRIPT roman_out end_POSTSUBSCRIPT (8)8{}^{(8)}start_FLOATSUPERSCRIPT ( 8 ) end_FLOATSUPERSCRIPT 7800±100plus-or-minus7800100-7800\pm 100- 7800 ± 100 4400±100plus-or-minus4400100-4400\pm 100- 4400 ± 100
Total C-stat / d.o.f. 4095.6 / 3935 3920.9 / 3935
Total Expected C-stat 3992.1±90.1plus-or-minus3992.190.13992.1\pm 90.13992.1 ± 90.1 3995.7±90.1plus-or-minus3995.790.13995.7\pm 90.13995.7 ± 90.1
  • 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 1051photonss1keV1superscript1051photonssuperscripts1superscriptkeV110^{51}~{}{\rm photons~{}s^{-1}~{}keV^{-1}}10 start_POSTSUPERSCRIPT 51 end_POSTSUPERSCRIPT roman_photons roman_s start_POSTSUPERSCRIPT - 1 end_POSTSUPERSCRIPT roman_keV start_POSTSUPERSCRIPT - 1 end_POSTSUPERSCRIPT at 1 keV. (3) Scale parameter of the reflection component. (4) Normalization of the modified blackbody in 1034cm0.5superscript1034superscriptcm0.510^{34}~{}{\rm cm^{0.5}}10 start_POSTSUPERSCRIPT 34 end_POSTSUPERSCRIPT roman_cm start_POSTSUPERSCRIPT 0.5 end_POSTSUPERSCRIPT. (5) Temperature of the modified blackbody in eV. (6) Logarithm of the ionization parameter ξ𝜉\xiitalic_ξ of the ionized wind components in ergcms1ergcmsuperscripts1{\rm erg~{}cm~{}s^{-1}}roman_erg roman_cm roman_s start_POSTSUPERSCRIPT - 1 end_POSTSUPERSCRIPT. (7) Column density NHN{{}_{H}}italic_N start_FLOATSUBSCRIPT roman_H end_FLOATSUBSCRIPT of the ionized wind components in 1022superscript102210^{22}10 start_POSTSUPERSCRIPT 22 end_POSTSUPERSCRIPT cm22{}^{-2}start_FLOATSUPERSCRIPT - 2 end_FLOATSUPERSCRIPT. (8) Velocity of the outflowing ionized wind components) in kms1kmsuperscripts1{\rm km~{}s^{-1}}roman_km roman_s start_POSTSUPERSCRIPT - 1 end_POSTSUPERSCRIPT. (9) Covering factor of the pion component. (10) Fe abundance of the pion component.

Refer to caption
Figure 9: The 6.67.36.67.36.6-7.36.6 - 7.3 keV energy band simulated spectra and model curve. The magenta point, black cross point, blue line, and red dashed line represent the XRISM simulated spectra based on the M3 model, and the XRISM simulated spectra based on the M4 model, the M3 model, and the M4 model, respectively. For plotting, the binning factor was set to 10 for the simulated XRISM spectra.
Refer to caption
Figure 10: The 0.680.940.680.940.68-0.940.68 - 0.94 keV energy band simulated spectra and model curve. The magenta point, black cross point, blue line, and red dashed line represent the XRISM simulated spectra based on the M3 model, and the XRISM simulated spectra based on the M4 model, the M3 model, and the M4 model, respectively. For plotting the simulated XRISM spectra, we used the obin command for optimal binning

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 95(n108cm3)0.5RS95superscript𝑛superscript108superscriptcm30.5subscript𝑅S95\left(\frac{n}{10^{8}{\rm cm}^{-3}}\right)^{-0.5}~{}R_{\rm S}95 ( divide start_ARG italic_n end_ARG start_ARG 10 start_POSTSUPERSCRIPT 8 end_POSTSUPERSCRIPT roman_cm start_POSTSUPERSCRIPT - 3 end_POSTSUPERSCRIPT end_ARG ) start_POSTSUPERSCRIPT - 0.5 end_POSTSUPERSCRIPT italic_R start_POSTSUBSCRIPT roman_S end_POSTSUBSCRIPT and 1.8(n105cm3)0.5pc1.8superscript𝑛superscript105superscriptcm30.5pc1.8\left(\frac{n}{10^{5}{\rm cm}^{-3}}\right)^{-0.5}~{}{\rm pc}1.8 ( divide start_ARG italic_n end_ARG start_ARG 10 start_POSTSUPERSCRIPT 5 end_POSTSUPERSCRIPT roman_cm start_POSTSUPERSCRIPT - 3 end_POSTSUPERSCRIPT end_ARG ) start_POSTSUPERSCRIPT - 0.5 end_POSTSUPERSCRIPT roman_pc, 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 7similar-toabsent7\sim 7∼ 7 keV absorption lines can be attributed to Fe xxv/xxvi absorption lines, the spectra derived herein suggested that the 7similar-toabsent7\sim 7∼ 7 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 7similar-toabsent7\sim 7∼ 7 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&ouml;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