Arxiv:1009.2968V2 [Astro-Ph.CO] 25 Oct 2010 01 Hne L 03.Tewi,Woetmeaueand Al

arXiv:1009.2968v2 [astro-ph.CO] 25 Oct 2010 01 hne l 03.TeWI,woetmeaueand al. temperature et Davé whose about 2006; WHIM, be to The Ostriker thought is & 2003). density al. Cen et 2006; Chen Fang 2001; & Ostr Cen & Cen 1999; (e.g. cos- formation of structure simulations large-scale dynamics gas mological and N-body several on (WHIM) based medium intergalactic inter warm-hot the as called of medium, filaments tic as structure large-scale the along ist Introduction 1. ftegscmae ihtehg ee fteteGlci for Galactic the the of level high the densi with low compared extremely gas gas, the the WHIM of of because the not from could others emission the report the while (2009) al. of al. detection et et Galeazzi significant Galeazzi et and the al. 2008; (2008) Takei al. al. et 2007; et et Werner al. Kaastra Werner 2009). et 2007; 2002; Bregman Galeazzi al. 2007a; 2003; al. et search al. et (Sołtan to Finoguenov emission us 2003; WHIM enabled the spectroscopy for imaging resolution high ob- universe. the if of structure WHIM, evolution the and the about mation information of rich us properties bring chemical would served, cam and CCD Thermal such detectors present-day eras. the with detect to hard is

PASJ: c 08 srnmclSceyo Japan. of Society Astronomical 2018. ag ato ayn ntelcluies stogtt ex- to thought is universe local the in baryons of part large A eetCadaadXMNwo bevtoswt the with observations XMM-Newton and Chandra Recent ul srn o.Jpn,1– , Japan Soc. Astron. Publ. background cm aatcadCBcmoet.W osrie h intensitie filament the the constrained for We fits components. spectral CXB the and clus and the Galactic of radii, radii virial virial memb the the the of toward of outside Scu radii region the virial the near in the taken emission to were determined which were profiles data Galactic dance background the the both used included also which We the background in the taken as spectra used The were emission. redshifte Galactic the the the of by intensity fected measured the rad because virial estimation the beyond regions the including Nov.–Dec., 2005 uaucridotfu bevtosi h uecutr n supercluster: the in observations four out carried Suzaku × 10 − e words: Key esuidtehg eprtr lsai h ieto fth of direction the in plasma temperature high the studied We − 3 8 ( Z htn cm photons / 0 . 4 1 eateto hsc,TkoMtooia nvriy 1-1 University, Metropolitan Tokyo Physics, of Department Z ⊙ ) 1 oioY Y Y. Noriko 3 − aais lses niiulA81 20,adA81 – A2811) and A2804, individual(A2801, clusters: galaxies: al nttt o srpyisadSaeRsac,Massa Research, Space and Astrophysics for Institute Kavli 1 nttt fSaeadAtoatclSine(SS,Japan (ISAS), Science Astronautical and Space of Institute td fteItalse n neglci Medium Intergalactic and Intracluster the of Study / − 2 2 ( s L − ?? / T 1 oueS Kosuke 5Mpc) 25 5 — , arcmin = hsc eatet nvriyo im,CrlGbe,F 3 FL Gables, Coral Miami, of University Department, Physics nteSupo uecutrwt Suzaku with Supercluster Sculptor the in 10 2 amasaki NASA -- ohndi hok,Sgmhr,Kngw 252-5210 Kanagawa Sagamihara, Chuo-ku, Yoshinodai, 3-1-1 5 − 7 − − n 10 and K 2 Rcie 00Jl 4 cetd21 etme 15) September 2010 accepted 14; July 2010 (Received ato 1 epciey s2 as respectively, , 7MsahstsAe,Cmrde A019 USA 02139, MA Cambridge, Ave., Massachusetts 77 / / odr pc lgtCne,Genet D271 USA 20771, MD Greenbelt, Center, Flight Space Goddard 2 , hc orsod oa vrdensity, over an to corresponds which , 3 , 1 aaaO Takaya ihr .K L. Richard − 6 –10 − hashi 4 galac- cm k elley iker σ [email protected] for- Abstract , − 4 ed e- ty 3 pe iis h nest fO of intensity The limits. upper - naiG Anjali , , , es sarsl,tetmeaueo h lsesdecreased clusters the of temperature the result, a As ters. 2 ie,epcal hs rmoye,wr togyaf- strongly were oxygen, from those especially lines, d io hs lses h td eddpeiebackground precise needed study The clusters. these of ii 92 tardhf of redshift a at 1992) l 09 ohn ta.20)bte hntepeiu satel previous the than better 2009) al. et Yoshino 2009; Suzaku al. that shown e have can studies constrai Past XIS th in emission. of important WHIM estimation the most Reliable ing is 2008). emission inten al. Galactic the et on foreground Tamura constraint regio tight also outer a (see set the sity al. they from Suzaku, et de- lines with Reiprich (2007b) A2218 O of al. al 2009; redshifted et et Takei significant al. Although no (Sato et tected 2010). radii Bautz al. virial et 2007b; Hoshino the al. 2009; near et to clusters Sato emission of (ICM) 2007a; region medium outer intra-cluster 2007) the the al. observe et (Koyama can keV 1 we below sensitivity higher and level a et et Zappacosta Takei 2010; 2007; al. al. et Fang et 2010). 2009; Rasmussen al. et 2007; Kaast Buote 2005; al. 2007a; al. et al. et Fang Nicastro 2006; 2004; al. al. et (Fang et significantly Fujimoto more 2002; density al. possi- WHIM et WHIM the the restricted of and study ble absorption-line other the spectrom- made the grating also On with eters observations XMM background. and intrinsic Chandra hand, the and emission ground ml,A81 21 o A2811 A2811, amely, o T Yoh ein usd ftevra ai ftemme clusters member the of radii virial the of outside regions einddntrqieetacmoetohrta the than other component extra require not did region h cltrsprlse Shekr&Ot19;Seitter 1991; Ott & (Schuecker supercluster Sculptor The eas uauXSi hrceie ylwrbackground lower by characterized is XIS Suzaku Because fO of s rcutr,adte esace o h xgnline oxygen the for searched we then and clusters, er n omcXryBcgon CB components. (CXB) Background X-ray Cosmic and upta akei cltrsprlse at supercluster Sculptor e iaiOaa ahoi oy 192-0397 Tokyo Hachioji, Minami-Osawa, po uecutr eprtr n ea abun- metal and Temperature supercluster. lptor , VII 5 ff , 3 n asmlaoG Massimiliano and hstsIsiueo Technology, of Institute chusetts cieycntanteGlci msin(ut et (Gupta emission Galactic the constrain ectively eopc xlrto Agency, Exploration Aerospace aauiT Takayuki n O and < δ neglci eim–Xry:di X-rays: – medium intergalactic VIII 0( 60 z 14 USA 3124, = ie ob esta . n 5.1 and 8.1 than less be to lines amura Z ff 0 / . e,A84 20 ein in regions A2801 A2804, set, 1i n ftercetsupercluster, richest the of one is 11 0 VII . 1 , Z indicates 3 aleazzi ⊙ z ) − = 1 / 0 2 . ( 0 ihSuzaku. with 108 5 L n / H 5Mpc) 25 < 1 . 6 × − ff 10 1 / use 2 − . 5 lites. n- ra l. n e - . , 2 K. Sato et al. [Vol. ,

Table 1. Basic properties of the clusters.

redshift virial radius∗ Average temperature NH 2 (Mpc/arcmin) k T (keV) (cm− ) h i 20 A2811 0.1079 1.96/16′5 5.0 1.54 10 × 20 A2804 0.108 0.88/7′3 1.0 1.78 10 A2801 0.108 1.25/10 5 2.0 1.70×1020 ′ × ∗ At a redshift of z = 0.108, 1′ corresponds to 119 kpc, and the virial radii, r180 = 1.95 h 1 √k T /10 keV Mpc (Markevitch et al. 1998) 100− h i

Table 2. Suzaku Observation logs of the Sculptor supercluster.

Region Seq.No. Obs.date (RA,Dec)∗ Exp. Afterscreening J2000 ksec (BI/FI) ksec h m s A2811 800005010 2005-11-28T04:01:40 (00 41 52.9, 28◦33′19′′) 30.6 29.3/29.2 h m s − A2811 oﬀset 800006010 2005-11-28T21:10:16 (00 40 42.7, 28◦44′00′′) 33.9 32.4/32.3 h m s − A2804 800007010 2005-12-02T08:51:05 (00 39 35.2, 28◦54′35′′) 31.6 30.0/29.8 h m s − A2801 800008010 2005-11-27T08:30:47 (00 38 31.3, 29◦02 48 ) 34.8 33.4/29.5 − ′ ′′ ∗ Average pointing direction of the XIS, written in the RA NOM and DEC NOM keywords of the event FITS ﬁles.

-28°20′

-28°30′

-28°40′

-28°50′

-29°00′

-29°10′

-29°20′ 0h42m 0h40m 0h38m

Fig. 1. Left: Combined XIS image of the four observations of the Sculptor supercluster in the 0.5–5.0 keV energy range. The observed XIS0, 1, 2, and 3 images were added on the sky coordinate, and smoothed with σ = 16 pixel 17′′ Gaussian. The instrumental background (NXB) were subtracted, and the exposure was corrected, though vignetting was not corrected.≃ The green circles show the extracted region for the spectral analysis as shown in subsection 3.4. The subtracted point source regions are indicated by the white dashed circles. Right: Galaxy distribution associated with the Sculptor supercluster from 2dFGRS catalogue. The galaxies are selected between 0.103 < z < 0.113, and the nearer (z < 0.108) and farther (z > 0.108) galaxies than the sculptor redshift correspond to the solid circles and crosses, respectively. The blue boxes show the observation regions with Suzaku, and the green circle and dashed circles are same as the left panel. composed of more than 20 clusters of galaxies (Einasto et al. lapping with our Suzaku field. On the other hand, Zappacosta 1997) including 5 Abell systems. It is located near the south et al. (2005) reported extended soft X-ray emission in 0.1 –0.3 galactic pole. Since WHIM is thought to distribute along the keV band with ROSAT PSPC, and the region corresponds to filamentary structures over several tens of Mpc scale connect- 2–4 degrees east of our observations. They suggested that ing clusters of galaxies, it is likely that superclusters contain the∼ diffuse emission was tracing the WHIM. significant amount of WHIM. The Sculptor supercluster has Fang et al. (2010) reported a 3σ detection of an oxygen ab- been observed with ROSAT and ASCA. Spiekermann (1996) sorption line from the WHIM in the “Sculptor Wall”at z 0.03 and Obayashi et al. (1998) found no indications of a largely ex- using the Chandra and XMM grating observations of a∼ blazer tended X-ray emission in their observed regions, which is over- behind this structure. The Sculptor Wall is located much No.] SuzakuobservationsoftheSculptorSupercluster 3 closer than the Sculptor supercluster in redshift and 10 de- rors. In order to subtract the non X-ray background (NXB), we grees westward from our observations. The estimated∼ hydro- employed the dark Earth database by the “xisnxbgen” Ftools 6 3 gen density is n = 6 10− cm− under the assumption of the task. uniform absorber distribution× along the line of sight between For the spectral fits of the ICM emission, we generated 0.028 < z < 0.032 using the O VII column density. Zappacosta Ancillary Response Files (ARFs) for A2801/A2811 assumed et al. (2010) also reported results from the WHIM search in the β-model surface brightness profile, β = 0.8/0.82 and rc = “Pisces–Cetus Supercluster” at z = 0.062 and “Farther Sculptor 0.′88/1.′6, respectively, reported in Obayashi et al. (1998), by Wall” at z = 0.128 in the same line of sight as “Sculptor Wall”. the “xissimarfgen” Ftools task (Ishisaki et al. 2007). We also In contrast to the result of “Sculptor Wall”, significant absorp- generated the ARFs for each observation assumed an uni- tions (σ> 3) were not detected from any of the strong H-like form sky of 20′ 20′ for the Galactic and Cosmic X-ray and He-like transitions of C, O, and Ne at the redshift. Background emission.× Note that we used the uniform ARFs 1 1 We use H0 = 70kms− Mpc− , ΩΛ = 1 ΩM = 0.73in thispa- for the spectral fits of A2804. We included the effect of the per. The redshifts of individual clusters− are 0.1079 for A2811 contaminations on the optical blocking filter of the XISs in the and the offset observations, 0.108 for A2801 and A2804, re- ARFs. Since the energy resolution also slowly degraded after spectively, taken from NED as shown table 1. At a redshift the launch, due to radiation damage, this effect was included of z = 0.108, 1′ corresponds to 119 kpc, and the virial radii, in the Redistribution Matrix File (RMF) by the “xisrmfgen” r = 1.95 h 1 √k T /10 keV Mpc (Markevitch et al. 1998), is Ftools task. 180 100− h i 1.96 Mpc (16.5′) for an average temperature of k T = 5.0 keV h i of A2811, 0.88 Mpc (7.3′) for k T = 1.0 keV of A2804, and 3. Results and Discussion 1.25 Mpc (10.5 ) for k T = 2.0 keVh i of A2801. Because of this ′ 3.1. Strategy of the spectral fits redshift, we can separateh i the WHIM O line from the supercluster relatively easily from the Galactic emission. Throughout We analyzed the spectra of the Sculptor supercluster using this paper we adopt the Galactic hydrogen column density of the following way. N = 1.54,1.78,1.70 1020 cm 2 (Dickey & Lockman 1990) in H − 1. Extraction of the point-like sources for each observation. the direction of A2811,× A2801, A2804, respectively. Because 2. Background estimation for the Galactic and Cosmic X- the system is located near the south galactic pole, the Galactic ray background (CXB). absorption and emission are both low and have no complex feature. Unless noted otherwise, the solar abundance table is (a) Spectral fits of the outer region of the virial radii given by Anders & Grevesse (1989), and the errors are in the of the member clusters in A2811 offset region as 90% confidence region for a single interesting parameter. shown in figure 1. (b) Estimation of the background models using the 2. Observations and Data Reduction Suzaku observations near the Sculptor supercluster as shown in figure 3. 2.1. Observations (c) Comparing the above two results, we determine the Suzaku observed the Sculptor supercluster in November and background model. December 2005. The observation logs are given in table 2, and 3. Spectral analysis of the ICM emission of A2801, A2804, the XIS image in 0.5–5.0keV is shownin figure1 left. Figure 1 A2811 to the virial radius. right shows a galaxy density map in 0.103 < z < 0.113 using the 4. Constraining the O emissions from WHIM in the fila- 2dF Galaxy Redshift Survey (2dFGRS) 1. The XIS was oper- ment region between the clusters. ated in the Normal clocking mode (8 s exposure per frame), with the standard 5 5 or 3 3 editing mode. During these 3.2. Fore- and Background Estimation observations, the significant× e×ffect of the Solar Wind Charge eXchange (SWCX) was not confirmed in ACE data 2, while Using ”wavdetect” tool in CIAO3, we searched for point-like it was known that the SWCX affected some Suzaku spectra in sources in the Suzaku images in the energy range 0.5-2.0 keV. the lower energy range as reported in Fujimoto et al. (2007), As shown by the white dashed circles in figure 1, we detected Gupta et al. (2009), Koutroumpa et al. (2009) and Yoshino et and subtracted 7 point-like sources using 2 arcmin radius which al. (2009). corresponded to the HPD of the Suzaku XRT. The flux level of the sources were > 1.5 10 12 erg cm 2 s 1 in 0.5-2.0 keV, and 2.2. Data Reduction − − − contribution from unresolved× sources was taken into account We used version 2.1 processing data, and the analysis was when we subtract the CXB spectrum. After subtracting the performed with HEAsoft version 6.6.3 and XSPEC 12.5.0ac. point sources, we fitted the spectra integrated over the A2811 In the analysis of XIS data, we selected ELEVATION > 15◦ offset region excluding the virial region (r < r180) of A2811, of the data set to remove stray-light from the day Earth limb. where a few galaxies existed as shown in figure 1 right, with Event screening with cut-off rigidity (COR) was not performed the following models including the local hot bubble (LHB) and in our data. Note that, even if we select the data with COR>8 milky way halo (MWH) as the Galactic foregrounds and the GV, the resultant fit parameters all stay within the statistical er- CXB as a background: apecLHB +phabs (apecMWH +powCXB). It was suggested that SWCX made a significant× contribution to 1 http://www.mso.anu.edu.au/2dFGRS/ 2 http://www.srl.caltech.edu/ACE/ASC/ 3 http://cxc.harvard.edu/ciao/ 4 K. Sato et al. [Vol. ,

case 1 case 2 case 3 case 4 −1 −1 −1 −1 0.1 0.1 0.1 0.1 keV keV keV keV −1 −1 −1 −1

0.01 0.01 0.01 0.01

normalized counts s 10−3 normalized counts s 10−3 normalized counts s 10−3 normalized counts s 10−3 4 4 4 4 2 2 2 2

χ 0 χ 0 χ 0 χ 0 −2 −2 −2 −2 −4 −4 −4 −4 0.51 2 0.51 2 0.51 2 0.51 2 Energy (keV) Energy (keV) Energy (keV) Energy (keV)

case 5 −1 0.1 keV −1 Fig. 2. The panels show the resultant ﬁts of the outer region of the virial radii of the member clusters with 0.01 the models as shown in table 3. The data are plotted by black and red crosses for BI and FI, and the black and red lines show the best-ﬁt model for the BI and FI spectra, respectively. The CXB components are normalized counts s 10−3 shown in a magenta line, and the Galactic emissions are indicated by cyan and blue lines, respectively. 4 2

χ 0 −2 −4 0.51 2 Energy (keV)

Table 3. Summary of the parameters of the ﬁts of the outer region of the virial radii of the clusters. We assumed zero redshift and a solar metallicity in the ﬁts.

case1 case2 case3 case4 case5 LHB +0.01 kT (keV) 0.10 0.01 0.10(fixed) (fixed) (fixed) (fixed) 3 +−1.57 +0.37 ← +0.36 ← +0.31 ← +0.32 Norm∗ ( 10− ) 2.62 1.00 2.61 0.37 2.36 0.36 2.32 0.31 2.32 0.32 MWH × − − − − − +0.03 +0.03 +0.03 kT (keV) 0.32 0.03 0.32 0.02 0.31 0.03 0.30(fixed) (fixed) 4 +−0.71 +−0.67 +−0.63 +0.50 ← +0.51 Norm∗ ( 10− ) 4.18 0.39 4.18 0.67 3.60 0.63 3.69 0.50 3.54 0.50 CXB × − − − − − +0.08 +0.08 Γ 1.22 0.06 1.22 0.08 1.4(fixed) (fixed) (fixed) Norm , ( 10 4) 7.93+−0.64 7.93+−0.65 8.99+0.48 ←9.01+0.48 9.6← (fixed) ∗ † × − 0.62 0.62 0.47 0.47 χ2/ν 218/−167 218/−168 231/−169 231/−170 235/171

Table 4. List of χ2/dof in case 4 in table 3 for the ﬁts of the nominal and considering the systematic errors such as the NXB level and contaminant of OBF. For details, see text.

uncertainties nominal NXB contaminant +10% -10% +10% -10% χ2/dof ...... 231/170 247/170 227/170 231/170 231/170 No.] SuzakuobservationsoftheSculptorSupercluster 5

Table 5. Suzaku Observation logs of NGC 253 oﬀset and Q0109-3518.

Region Seq.No. Obs.date (RA,Dec)∗ Exp. Afterscreening Distance† J2000 ksec (BI/FI)ksec degree h m s NGC 253 oﬀset 803004010 2008-12-29 (00 48 22.2, 25◦02′51′′) 56.7 53.7/53.7 3.75 h m s − Q0109-3518 703037010 2008-05-20 (01 11 45.0, 35◦03′10′′) 30.0 26.0/25.9 9.01 − ∗ Average pointing direction of the XIS, written in the RA NOM and DEC NOM keywords of the event FITS ﬁles.

† Distance from A2811

Fig. 3. Extracted regions of NGC 253 oﬀset (Left) and Q0109-3518 (Right). Green boxes corresponds to the extracted region for the spectral ﬁts, and the elliptical regions are eliminated from the spectral region.

NGC253 offset region; case 2 NGC253 offset region; case 5 Q0109−3518 region; case 2 Q0109−3518 region; case 5 −1 −1 −1 −1 0.1 0.1 0.1 0.1 keV keV keV keV −1 −1 −1 −1

0.01 0.01 0.01 0.01

normalized counts s 10−3 normalized counts s 10−3 normalized counts s 10−3 normalized counts s 10−3 4 4 4 4 2 2 2 2

χ 0 χ 0 χ 0 χ 0 −2 −2 −2 −2 −4 −4 −4 −4 0.51 2 0.51 2 0.51 2 0.51 2 Energy (keV) Energy (keV) Energy (keV) Energy (keV)

Fig. 4. The spectra of NGC 253 offset and Q0109-3518. The notations in the figure are same as figure 2.

Table 6. Summary of the parameters of the ﬁts of the NGC 253 oﬀset region and Q0109-3518 region.

Region NGC253offset Q0109-3518 2 20 20 NH (cm− ) 1.42 10 1.93 10 case2× case5 case2× case5 LHB kT (keV) 0.10(fixed) (fixed) (fixed) (fixed) 3 +0.74 ← +0.65 ← +0.76 ← +0.72 Norm∗ ( 10− ) 2.45 0.78 2.56 0.66 2.93 0.86 2.57 0.72 MWH × − − − − +0.08 +0.07 kT (keV) 0.32 0.06 0.3(fixed) 0.34 0.06 0.3 (fixed) 4 +−1.14 +0.72 +−1.28 +0.84 Norm∗ ( 10− ) 2.59 1.02 3.00 0.71 2.98 1.17 2.79 0.83 CXB × − − − − +0.09 +0.13 Γ 1.38 0.09 1.4(fixed) 1.30 0.12 1.4 (fixed) Norm , ( 10 4) 10.26−+0.94 9.6(fixed) 8.49+−1.04 9.6 (fixed) ∗ † × − 0.89 0.98 χ2/ν 180/155− 185/158 204/−155 207/158

u ∗ Normalization of the apec component divided by the solid angle, Ω , assumed in the uniform-sky ARF calcula- 2 2 u 14 5 2 tion (20′ radius), Norm = R nenHdV /(4π(1 + z) DA)/Ω 10− cm− 400π arcmin− , where DA is the angular distance to the source. × 4 † The referred CXB normalization in Kushino et al. (2002) is 9.6 10− for Γ = 1.4 in the above unit. × 6 K. Sato et al. [Vol. , unabsorbed emission-line component (LHB) even in the case We also estimated the fraction of photons entering from out- of no flare in the solar-wind flux (e.g. Yoshino et al. 2009). side of the extracted regions using a simulator of the Suzaku However, since SWCX and thermal plasma is hardly distin- XRT/XIS system “xissim”. The photon contamination can guished with a limited spectral resolution of CCDs (Gupta et cause a significant effect in the case of cluster observations, al. 2009), we model the sum of SWCX and LHB as a single since much bright center is outside of the field and the point- unabsorbed thermal plasma. Here, we assumed zero redshift spread function of the Suzaku XRT shows an extended tail (see and a solar metallicity. We examined the spectral fits with the e.g. Sato et al. 2007a). If we assume a beta-model of A2811 model, by changing the parameters to be free or fixed. We emission as described in subsection 2.2 extending beyond the also assumed that there was no intracluster and/or intergalac- virial radius, a half of the simulated photons in the extracted tic emission outside of the virial radii. We discuss possible region, which is outside of the virial radius, is originated from intergalactic emission in subsection 3.5. The resultant param- outside of the extracted region. However, this emission has a eters of the fits are shown in table 3 and figure 2. As a result, flux of 1/50 of the estimated CXB level below 2 keV. Thus, the spectra were well-fitted by the sum of Galactic and CXB the effect∼ of the stray light is negligible in the spectral fits. components. We also tested the inferred photon index, Γ, and 3.3. Observations close to the Sculptor supercluster normalization of the CXB model in Kushino et al. (2002) after the point-source subtraction. Our results gave slightly lower in- In order to estimate the Galactic and CXB spectra pre- tensity than the Kushino et al. (2002) level, with poor statistics cisely, we looked into other observations which were carried above 2 keV. Thus, we employed the parameters of the case 5 out near the Sculptor supercluster. There are two observations: in table 3, derived with fixed spectral shape for the Galactic NGC 253 offset region, and Q0109-3518. Their observation emission and fixed shape and intensity for the CXB compo- logs are shown in table 5. We analyzed these data with the nent, respectively. We performed spectral fit for the clusters same criteria used for the fore- and background spectral esti- in the Sculptor supercluster with this fixed background in the mation as described in subsection 3.2. As mentioned in sub- following analysis. section 2.2, the contamination on the XIS optical blocking fil- We also examined whether or not an additional apec com- ter has increased after the launch, and the effective area in the ponent (with redshift either zero or supercluster value) gave a energy range below 1 keV were smaller in these offset data as significant improvement in χ2 against a change of degree of shown in figure 4. However, because this effect was properly freedom ∆ν. In the case of zero-redshift apec model, the im- included in the ARF, the observationswere useful in improving provement was ∆χ2 = 1 with ∆ν = 3 over the value of case 4 in the background estimation. table 3. This is not a significant change considering the F-test. After subtracting point-like sources as shown in figure 3, On the other hand, in the case of adding a cluster-redshift apec we carried out spectral fits with the same model, consisting of model (z = 0.108), there was no improvement (∆χ2 < 1) over LHB, MWH and CXB, by changing the parameters to be either the case 4 value, for a temperature range between 0.1–0.8 keV. free or fixed. The resultant parameters are shown in table 6 and Since the fit is not formally acceptable (χ2/ν = 231/170 in in figure 4. As a result, the spectra of these near-by regions are case 4), we investigated the χ2 behavior in 2 energy bands, also well-fitted with the backgroundmodel as shown in subsec- 0.4–1.0 and 1.0–5.0 keV, separately. As a result, the χ2-value tion 3.2. The intensities of the Galactic and CXB components in 1.0–5.0 keV (the reduced χ2 1.5) was much larger than are consistent with those for the region outside of the virial ra- that in 0.4–1.0 keV (the reduced∼χ2 1.0). Changing the pa- dius of the clusters in the Sculptor supercluster. rameters of the CXB component did∼ not improve the fit. The 3.4. Radial Profile of ICM Properties for each cluster spectral data in this high-energy range seem to fluctuate from bin to bin by more than the statistical limit. There are some fea- We extracted spectra from the regions 0 < r < 0.1 r180,0.1 < tures around the edges of Si-K (1.8 keV) and Au-L (3.4 keV), r < 0.3 r180, 0.3 < r < 0.5 r180, and 0.5 < r < 1 r180, centered h m s suggesting a problem in the XIS response function. However, on (RA, Dec) = (0 42 09.3, 28◦32 33 ) for A2811, 0 < r < − ′ ′′ the model reproduces average feature of the observed spectrum 0.4 r180, 0.4 < r < 0.6 r180, and 0.6 < r < 1 r180, centered on h m s quite well, yielding the CXB spectrum consistent with the pre- (RA, Dec) = (0 38 32.3, 29◦04′44′′) for A2801, and 0< r < − h m s vious Kushino et al. (2002) result, and the fit is good in the en- 0.5 r180,0.5 < r < 1.0 r180 centered on (RA, Dec) = (0 39 38.7, ergy range of our interest (< 1 keV). Therefore, we proceeded 28◦53′21′′) for A2804, respectively. The spectra with BI and with the spectral fit using these data. −FI for all regions were fitted simultaneously. Energy ranges We also investigated the effect of possible incorrect calibra- used for the fit are 0.4–7.1 keV for A2811 and A2801, and tion, i.e., NXB level and OBF contaminations, by artificially 0.4–5.0 keV for A2804, respectively. In the simultaneous fit, changing the values by 10% and comparing the resultant χ2 the Galactic and CXB components based on the case 4 values value. While the temperatures± of the LHB and MWH did not as described in subsection 3.2 were commonly included for all change within < 1% compared with the values for the case 1 in regions. We excluded the narrow energy band around Si K- table 3, the normalization and Γ of the CXB models changed edge (1.825–1.840 keV) because its response was not modeled by 10%. We also examined the other uncertainty in the OBF correctly. The energy range below 0.4 keV was also excluded contaminant∼ by changing the absorber thickness by 10%. The becausethe C edge(0.284keV) seen in the BI spectra could not temperatures of LHB and MWH and Γ and normalization± of be reproduced well in our data. These are the known problems the CXB did not change within < 1%, the normalizations of in the XIS response. the LHB and MWH changed by 10%. A list of χ2/dof by The BI and FI spectra of the ICM in A2801 and A2811were changing the systematic errors in±∼ case 4 is presented in table 4. well-represented by a single temperature model as shown in No.] SuzakuobservationsoftheSculptorSupercluster 7

Abell 2801 Abell 2801 Abell 2801

0.1 0 < r < 0.4 r 0.1 0.4 < r < 0.6 r 0.1 0.6 < r < 1 r

−1 180 −1 180 −1 180 keV keV keV −1 −1 −1

0.01 0.01 0.01

10−3 10−3 10−3 normalized counts s normalized counts s normalized counts s 4 4 4 2 2 2

χ 0 χ 0 χ 0 −2 −2 −2 −4 −4 −4 0.51 2 5 0.51 2 5 0.51 2 5 Energy (keV) Energy (keV) Energy (keV)

Abell 2804 (ICM 2T model + zbremss) Abell 2804 Abell 2811

0.1 0 < r < 0.5 r 0.1 0.5 < r < 1 r 0 < r < 0.1 r180 −1 −1 −1 180 180 keV keV keV 0.1 −1 −1 −1

0.01 0.01

0.01

10−3 10−3 normalized counts s normalized counts s normalized counts s 10−3 4 4 4 2 2 2 χ χ χ 0 0 0 −2 −2 −2 −4 −4 −4 0.51 2 5 0.51 2 5 0.51 2 5 Energy (keV) Energy (keV) Energy (keV)

Abell 2811 Abell 2811 Abell 2811

0.1 < r < 0.3 r180 0.3 < r < 0.5 r180 0.5 < r < 1 r180 −1 −1 −1

keV 0.1 keV 0.1 keV 0.1 −1 −1 −1

0.01 0.01 0.01 normalized counts s normalized counts s normalized counts s 10−3 10−3 10−3 4 4 4 2 2 2

χ 0 χ 0 χ 0 −2 −2 −2 −4 −4 −4 0.51 2 5 0.51 2 5 0.51 2 5 Energy (keV) Energy (keV) Energy (keV)

Fig. 5. The panels show the observed spectra for the regions of A2801, A2804, and A2811 which are denoted in the panels, and the data are plotted by black and red crosses for BI and FI, respectively. The estimated NXB components are subtracted, and the black and red lines show the best-fit model for the BI and FI spectra, respectively. The BI spectra of the ICM component are shown in an orange line. The CXB components are shown in a magenta line, and the Galactic emissions are indicated by cyan and blue lines, respectively. For A2804, the additional ICM and LMXB components are shown in orange and magenta dashed lines, respectively. The energy range around the Si K-edge (1.825–1.840 keV) is ignored in the spectral fits. The lower panels show the fit residuals in units of σ.

figure 5. The spectra except for the r < 0.3 r180 region of drop to about 60% of the peak temperature in the three clus- A2811 were all fitted with apec model, while these inner re- ters. This tendency is consistent with the previous Suzaku region spectra of A2811 needed fitting with vapec model. In sults (Bautz et al. 2009; Reiprich et al. 2009; Hoshino et al. this fit, the metal abundances were linked in the following way, 2010). While the metal abundances of α elements in the cen- Mg=Al, S=Ar=Ca, and Fe=Ni. As for A2804, the spectrum in tral region of A2811 are about the solar value, the iron abun- 0.5 < r < 1 r was represented by a single temperature model, dance is 0.4 times the solar value (see table 8). The obtained 180 ∼ but the data in the inner r < 0.5 r180 region needed an addi- temperatures are consistent with the ASCA results (Obayashi tional low mass X-ray binary (LMXB) component represented et al. 1998), but the present results give better accuracy and ra- by zbremss with kT = 10 keV as well as the apec ICM compo- dial distribution of temperature. The metal abundance in the nent. The A2804 emission is significantly affected by point- central region of A2801 is also 0.4 solar, which is mostly de- like sources due to their lower ICM temperature than those termined by iron. The metal abundance in the outer region of in other clusters. Results of the temperature and abundance A2801 and the wholeregionof A2804had to befixed at 0.3 so- from the spectral fit for individual regions are summarized in lar because of the poor photon statistics. The iron abundances table 7–8 and figure 6. The typical systematic error on the tem- of these clusters in the r < 0.1 r180 region may be slightly lower perature of A2811 due to the OBF contamination and change than those in other clusters (e.g. Tamura et al. 2008; Hoshino of NXB level by 10% are indicated. et al. 2010), while the alpha elements in A2811 show some The ICM temperatures± of the clusters clearly decrease to- enhanced abundance compared with the iron. wards the virial radii. The temperatures in 0.5–1 r180 region Kawaharada et al. (2010) reported a higher temperature and 8 K. Sato et al. [Vol. ,

Table 7. Summary of the parameters of the ﬁts of A2801 and A2804. Errors are 90% conﬁdence range of statistical errors, and do not include systematic errors. The solar abundance ratio of angr was assumed.

A2801 Region 0< r < 0.4 r180 0.4< r < 0.6 r180 0.6< r < 1 r180 +0.15 +0.24 +0.44 kT (keV) 2.33 0.15 1.43 0.11 1.08 0.25 +−0.12 − − Z (solar) 0.43 0.10 0.3(fix) 0.3(fix) −+1.2 +1.44 1 +5.10 2 Norm∗ 15.9 1.1 7.01 1.30 10− 6.51 4.32 10− χ2/do f 211/188− 219− /189× 244− /189× A2804 Region 0< r < 0.5 r180 0.5< r < 1 r180 +0.31 +0.23 +0.06 kT (keV) 1.15 0.13/0.63 0.28 0.91 0.07 Z (solar)− 0.3(fix)− 0.3(fix)− 3 +1.02 +2.19 +0.27 Norm† (10− ) 5.37 1.87/1.35 0.79 2.15 0.26 χ2/do f −176/168− 169/−171

∗ Normalization of the vapec component scaled with a factor of the selected region comparing to the assumed 2 2 20 5 2 image in “xissimarfgen”, Norm = f actor R nenHdV / [4π (1 + z) DA] 10− cm− arcmin− , where DA is the angular distance to the source. × u † Normalization of the apec component divided by the solid angle, Ω , assumed in the uniform-sky ARF calcula- 2 2 u 14 5 2 tion (20′ radius), Norm = R nenHdV /(4π(1 + z) DA)/Ω 10− cm− 400π arcmin− , where DA is the angular distance to the source. ×

Table 8. Summary of the parameters of the ﬁts of A2811. Errors are 90% conﬁdence range of statistical errors, and do not include systematic errors. The solar abundance ratio of angr was assumed.

A2811 Region 0< r < 0.1 r180 0.1< r < 0.3 r180 0.3< r < 0.5 r180 0.5< r < 1 r180 +0.22 +0.12 +0.39 +0.82 kT (keV) 5.38 0.16 5.17 0.12 4.97 0.37 2.63 0.57 +−0.60 +−0.40 +−0.13 − O(solar) 0.70 0.56 0.33 0.33 0.28 0.13 0.3 (ﬁx) +−0.66 +−0.47 − Ne(solar) 0.99 0.62 1.08 0.45 −+0.67 −+0.48 ↓ ↓ Mg,Al(solar) 0.90 0.64 0.80 0.46 −+0.47 −+0.34 ↓ ↓ Si(solar) 0.85 0.46 0.60 0.33 −+0.46 −+0.33 ↓ ↓ S,Ar,Ca(solar) 0.46 0.46 0.53 0.33 −+0.05 −+0.04 ↓ ↓ Fe,Ni(solar) 0.34 0.05 0.35 0.04 −+30 −+2.7 ↓+0.38 ↓+0.169 Norm∗ 648 30 81.5 2.7 8.08 0.38 0.968 0.167 χ2/do f 194/183− 235/183− 240/−188 215/−189

∗ Same as the unit of the normalization of A2801 in figure 6 entropy region in the northeastern edge of A1689 which was nents without an additional model required as shown in subsec- connected to an overdense filamentary structure of galaxies tion 3.2, however we investigate the upper limit of the intensi- outside of the cluster. As shown in figure 1 right, the galaxy ties of the additional O VII and O VIII lines at the cluster redshift distribution around A2811 shows an extension to the northwest from the WHIM origin. We fitted the data with the following direction rather than the southwest direction. We investigated model: apec +phabs (apec +pow +gaussian + LHB × MWH CXB OVII the difference of temperature between these directions in the gaussianOVIII). Here, the parameters of the Galactic and CXB radius range of 0.3 < r < 0.5 r180 and 0.5 < r < 1 r180 from components followed our previous fit in case 5. The tempera- A2811. Metal abundance was fixed to 0.3 solar because of the tures of the Galactic components and the photon index of the poor statistics. As a result, we found no temperature variation CXB component were fixed, and the normalizations were free between the two directions within the statistical error. in the fit. The redshifted central energies of the O VII and O VIII gaussian lines were fixed at 518 and 589 eV, respectively, with 3.5. Upper limit of Oxygen emission lines no intrinsic width of the lines assumed. We employed the in- Finally, we fitted the spectra taken outside of the virial radii crement of ∆χ2 = 4 as the measure for the 2σ upper limits of of the clusters. The region is the A2811 offset region, supposed the line intensities. This value also corresponds to the 95% to be the filament connecting the clusters. This region is the limit for an additional single parameter (single line intensity) same as the one we used to estimate the fore- and background in the F test. As a result, we determined the upper limit of the 8 2 1 emissions as shown in subsection 3.2. The spectrum in this O VII and O VIII lines to be 8.1 and 5.1 10− photons cm− s− 2 × region is well-represented by the Galactic and CXB compo- arcmin− as the 2σ confidence limits as shown in table 9 and No.] SuzakuobservationsoftheSculptorSupercluster 9

10 Table 9. Resultant intensities of O VII and O VIII emission lines with 2σ conﬁdence range. Abell 2811

5 O VII Center(keV) 0.518(fixed) Abell 2801 Sigma 0 (fixed) 8 2 1 2 +3.9 2 I( 10− photons cm− s− arcmin− ) 4.2 × 3.9 O VIII − Center(keV) 0.589(fixed) 1 Temperature (keV) Abell 2804 Sigma 0 (fixed) 8 2 1 2 +5.1 I( 10− photons cm− s− arcmin− ) 0.0 0.0 × − 0.5 0 0.2 0.4 0.6 0.8 1

Projected radius (r/r180) and electron-to-hydrogen number density ratio n /n = 1.2 for Fig. 6. Radial temperature proﬁles of A2811, A2801, and e H A2804 against the projected radius in unit of the virial ra- ionized gas give the following condition: dius. Black and Light gray lines show typical systematic 5 3 1/2 1/2 nH < 1.6 10− cm− (Z/0.1 Z )− (L/25 Mpc)− . (1) changes of the best-ﬁt values by varying the thickness of the × ⊙ OBF contaminant and NXB levels by 10%, respectively. The corresponding overdensity δ n /n¯ is ± ≡ H H 1/2 1/2 δ< 60 (Z/0.1 Z )− (L/25 Mpc)− . (2) ⊙

−1 Compared with the previous Suzaku observations (Takei et 0.1 al. 2007b; Tamura et al. 2008), our upper limits for the O VII keV −1 and O VIII lines imposed slightly tighter constraints on the WHIM emission. The present result is consistent with the ex- 0.01 pected values from numerical simulations, indicating a baryon over-density ρ/ρcrit 10 100 by Dolag et al. (2006) for exam- ple. We note that the∼ dense− part of WHIM is found only in the

normalized counts s 10−3 superclusters and the gas is likely to be inhomogeneous with 4 a wide range of temperature and density. In the case of the 2

χ 0 bridge region between A222 and A223, Werner et al. (2008) −2 −4 showed the temperature to be 0.9 keV. Such a high temper- 0.51 2 Energy (keV) ature (also with high density) yields lower emissivity of O VII and O VIII lines. Thus, estimation of temperature and density is Fig. 7. The panels show the resultant fits for the constraint of the inten- not straightforward. sities of the O VII and O VIII lines in the filament region between clusters As mentioned in section 1, Fang et al. (2010) reported a of the Sculptor supercluster. The O VII emission line is shown in a green line, and the notations of the other lines are same as shown in figure 2. significant detection of O VII absorption line in the direction of the “Sculptor Wall” at z 0.03. Zappacosta et al. (2010) also searched for the C, O,∼ and Ne absorption lines from the figure 7. Even if we allowed the redshift of the O VII and O VIII “Pisces–Cetus Supercluster” and the “ Farther Sculptor Wall” lines to vary to give the peak energies in 490–540 and 560– in the same line of sight as the “Sculptor Wall”. Even though 620 eV, respectively, the upper limit did not change within the the Sculptor supercluster and the above studied filaments are statistical errors. Our resultant values are almost consistent or not physically connected with each other, we can compare the give slightly stronger constraint compared with the results in estimated hydrogen density. Our upper limit on the density is Takei et al. (2007b). This is due to the combination of less twice larger than the “Sculptor Wall” value reported by Fang coupling between the redshifted and the Galactic O lines, to et al. (2010). Also, the overdensities derived for the “Pisces– slightly lower Galactic foreground intensity and to lower inter- Cetus Supercluster” and “Farther Sculptor Wall” are 1/4 stellar absorption. times of our result. This indicates a higher sensitivity of t∼he ab- In order to constrain the WHIM density based on the O line sorption line measurement in the detection of WHIM. In fact, upper limits, we show a wedge diagram of galaxy distribution Suzaku cannot probably detect the O line emission from the in figure 8. This is expressed in R.A. vs. redshift projected for “Sculptor wall” because the redshift is too low to separate from the declination range 32◦ < Dec.< 26◦ in the direction of the Galactic emission lines. Thus, future quantitative studies the Sculptor supercluster.− The galaxy concentration− around the of the WHIM will need X-ray optics with larger effective area Sculptor supercluster around z 0.11is obviousand we assume and instruments with much higher energy resolution such as 25 Mpc for the line-of-sight depth∼ of the WHIM distribution, microcalorimeters, which would play a key role in future in- which∼ corresponds to 0.105 < z < 0.111. As mentioned in Takei vestigation of the cosmic plasma diagnostics. et al. (2007b), we constrain the density of the WHIM cloud assuming the temperature to be T = 2 106 K. The line intensity, I < 8.1 10 8 photons cm 2 s 1 ×arcmin 2 at z=0.1079, × − − − − 10 K. Sato et al. [Vol. , References

Anders, E., & Grevesse, N. 1989, Geochim. Cosmochim. Acta, 53, 197 Bautz, M. W., et al. 2009, PASJ, 61, 1117 Bregman, J. N. 2007, ARA&A, 45, 221 Buote, D. A., Zappacosta, L., Fang, T., Humphrey, P. J., Gastaldello, F., & Tagliaferri, G. 2009, ApJ, 695, 1351 Cen, R., & Ostriker, J. P. 1999, ApJ, 514, 1 Cen, R., & Fang, T. 2006, ApJ, 650, 573 Cen, R., & Ostriker, J. P. 2006, ApJ, 650, 560 Chen, X., Weinberg, D. H., Katz, N., & Davé, R. 2003, ApJ, 594, 42 Davé, R., et al. 2001, ApJ, 552, 473 Dickey, J. M., & Lockman, F. J. 1990, ARA&A, 28, 215 Dolag, K., Meneghetti, M., Moscardini, L., Rasia, E., & Bonaldi, A. 2006, MNRAS, 370, 656 Einasto, M., Tago, E., Jaaniste, J., Einasto, J., & Andernach, H. 1997, A&AS, 123, 119 Fang, T., Davis, D. S., Lee, J. C., Marshall, H. L., Bryan, G. L., & Canizares, C. R. 2002, ApJ, 565, 86 Fang, T., Canizares, C. R., & Yao, Y. 2007, ApJ, 670, 992 Fang, T., Buote, D. A., Humphrey, P. J., Canizares, C. R., Zappacosta, L., Maiolino, R., Tagliaferri, G., & Gastaldello, F. 2010, ApJ, 714, Fig. 8. Wedge diagram expressed in R.A. vs. redshift projected on 1715 32◦ < Dec. < 26◦. The data are taken from 2dFGRS as mentioned in Finoguenov, A., Briel, U. G., & Henry, J. P. 2003, A&A, 410, 777 −subsection 2.1 and− figure 1. About 44,000 galaxies exist in the direction Fujimoto, R., et al. 2004, PASJ, 56, L29 of the map and the redshift, z < 0.2. The red and green dashed line show the direction of our observations and “the Sculptor Wall” at z 0.03 Fujimoto, R., et al. 2007, PASJ, 59, 133 in Fang et al. (2010), respectively. The red and green arrows∼ also Galeazzi, M., Gupta, A., Covey, K., & Ursino, E. 2007, ApJ, 658, show the Sculptor supercluster and “the Sculptor Wall”, respectively. 1081 Galeazzi, M., Gupta, A., & Ursino, E. 2009, ApJ, 695, 1127 Gupta, A., Galeazzi, M., Koutroumpa, D., Smith, R., & Lallement, R. 4. Summary 2009, ApJ, 707, 644 Hoshino, A., et al. 2010, PASJ, 62, 371 Suzaku carried out observations of the Sculptor superclus- Ishisaki, Y., et al. 2007, PASJ, 59, 113 ter in four fields: A2811, A2811 offset, A2804, and A2801 Kaastra, J. S., Lieu, R., Tamura, T., Paerels, F. B. S., & den Herder, regions. After estimation of the foreground and background J. W. 2003, A&A, 397, 445 Kaastra, J. S., Werner, N., Herder, J. W. A. d., Paerels, F. B. S., de emission, we examined the annular cluster spectra to the virial Plaa, J., Rasmussen, A. P., & de Vries, C. P. 2006, ApJ, 652, 189 radius. The resultant temperature of the clusters showed sig- Kawaharada, M., et al. 2010, ApJ, 714, 423 nificant decrease from the central region toward the virial ra- Koutroumpa, D., Lallement, R., Raymond, J. C., & Kharchenko, V. dius, consistent with the previous report for other clusters. Iron 2009, ApJ, 696, 1517 abundance suggests a nearly uniform distribution. We also Koyama, K., et al. 2007, PASJ, 59, 23 searched for the oxygen line emission outside of the virial Kushino, A., Ishisaki, Y., Morita, U., Yamasaki, N. Y., Ishida, M., radii of the clusters, correspondingto the filament region in the Ohashi, T., & Ueda, Y. 2002, PASJ, 54, 327 Sculptor supercluster. The spectral fits on the filament region Markevitch, M., et al. 1998, ApJ, 503, 77 data did not require an additional component except for the Nicastro, F., et al. 2005, Nature, 433, 495 Galactic and CXB components. We constrained the intensities Obayashi, H., Makishima, K., & Tamura, T. 1998, PASJ, 50, 537 8 Rasmussen, A. P., Kahn, S. M., Paerels, F., Herder, J. W. d., Kaastra, of redshifted O VII and O VIII lines to be less than 8.1/5.1 10− 2 1 2 × J., & de Vries, C. 2007, ApJ, 656, 129 photons cm− s− arcmin− , respectively, as the 2σ upper lim- 5 3 Reiprich, T. H., et al. 2009, A&A, 501, 899 its. The intensity of O VII line indicates nH < 1.6 10− cm− Sato, K., et al. 2007a, PASJ, 59, 299 1/2 1/2 × (Z/0.1 Z )− (L/25 Mpc)− , and the over density, δ < 60 Sato, K., Tokoi, K., Matsushita, K., Ishisaki, Y., Yamasaki, N. Y., ⊙ 1/2 1/2 (Z/0.1 Z )− (L/25 Mpc)− . Ishida, M., & Ohashi, T. 2007b, ApJL, 667, L41 ⊙ Schuecker, P., & Ott, H.-A. 1991, ApJL, 378, L1 Authors are grateful to Noel Bermann for supporting the Seitter, W. C. 1992, Digitised Optical Sky Surveys, 174, 367 data analysis. We also thank the referee for providing valu- Spiekermann, G. 1996, Roentgenstrahlung from the Universe, 615 able comments. K.S is supported by a JSPS Postdoctral fel- Sołtan, A. M., Freyberg, M. J., & Hasinger, G. 2002, A&A, 395, 475 lowship for research abroad. Part of this work was finan- Takei, Y., Henry, J. P., Finoguenov, A., Mitsuda, K., Tamura, T., cially supported by the Ministry of Education, Culture, Sports, Fujimoto, R., & Briel, U. G. 2007a, ApJ, 655, 831 Takei, Y., et al. 2007b, PASJ, 59, 339 Science and Technology, Grant-in-Aid for Scientific Research Tamura, T., et al. 2008, PASJ, 60, 695 Nos. 20340068, 20340041, 20840051, 21224003, 21740134. Tawa, N., 2008, Ph.D thesis, Osaka university Yoshino, T., et al. 2009, PASJ, 61, 805 No.] SuzakuobservationsoftheSculptorSupercluster 11

Werner, N., Finoguenov, A., Kaastra, J. S., Simionescu, A., Dietrich, J. P., Vink, J., Böhringer, H. 2008, A&A, 482, L29 Zappacosta, L., Maiolino, R., Mannucci, F., Gilli, R., & Schuecker, P. 2005, MNRAS, 357, 929 Zappacosta, L., Nicastro, F., Maiolino, R., Tagliaferri, G., Buote, D. A., Fang, T., Humphrey, P. J., & Gastaldello, F. 2010, ApJ, 717, 74