NGC 4388 - SPECTRAL STUDIES OF THE FIRST SEYFERT 2 SEEN BY INTEGRAL

V. Beckmann'a2, N. Gehrels', P. Fa~re~,~,T. J.-L. Cour~oisier~~~,R. Walter3>4, J. Malzac5i6, and P.-0. Petrucci7 'NASA Goddard Space Flight Center, Code 661, Greenbelt, MD 20771, USA Joint Center for Astrophysics, Department of Physics, University of Maryland, Baltimore County, MD 21250 31NTEGRAL Science Data Centre, Chemin d'Ecogia 16, 1290 Versoix, Switzerland 40bservatoire de Gen&ve,51 Ch. des Maillettes, 1290 Sauverny, Switzerland 'Centre d'Etude Spatiale des Rayonnements, 31028 Toulouse, France 'Institute of Astronomy, University of Cambridge, Madingley Road, Cambridge CB3 OHA, United Kingdom 7Laboratoire d'Astrophysique de Grenoble, BP 53X, 38041 Grenoble Cedex, France

ABSTRACT is believed to be due to a geometrical effect. The emission of the super massive black hole at the AGN core is probably absorbed by a surrounding torus. We present first INTEGRAL and XMM-Newton ob- The torus not only hides the broad line region, caus- servations of a Seyfert , the type 2 AGN NGC ing the characteristic Seyfert 2 optical spectra with 4388. Several INTEGRAL observations performed narrow emission lines, but also absorbes of the in 2003 allow us to study the spectrum in the 20 low X-ray emission. As the soft X-ray region is - 300 keV range. In addition two XMM-Newton very sensitive to absorption in the line of sight, the observations give detailed insight into the 0.2 - 10 study of the hard X-ray emission above - 3 keV is keV emission. Comparison with previous observa- a powerful tool to investigate the central engine of tions by BeppoSAX, SIGMA and CGROIOSSE show Seyfert 2 (Ghisellini, Haardt & Matt 1994). In addi- that the overall spectrum for soft X-rays up to the tion to this the variation of the Fe Ko line seems -rays can be described by a highly absorbed to vary differently with respect to the continuum (NH 'v 2.7 x 1023cm-2) and variable non-thermal in several cases (Iwasawa et al. 2003; Fabian 1977; component in addition to constant non-absorbed Malzac & Petrucci 2002), which might be caused by thermal emission (T E 0.8keV) of low abundance the interaction of the core emission with the absorb- (2 - 5%Za), plus a constant Fe Ko line. The ing and reflecting material. hard X-ray component is well described by a sim- ple power law with a mean photon index of K' = 1.7. At higher energies absorption becomes iess During the INTEGRAL observations the flux at 100 important revealing the emission from the keV increased by a factor of 1.5. The analysis of close environment of the central black hole. XMM-Newton data implies that the emission below Comptonized hard X-ray radiation from a hot 3 keV is decoupled from the AGN and probably due (- lOOkeV) plasma presumably forming an ac- to extended emission as seen in Chandra observa- cretion disk corona (Haardt et al. 1994), and tions. The constant iron line emission is apparently Compton reflection of these hard X-rays on also decoupled from the direct emission of the central a cool accretion disk (George & Fabian 1991; engine and likely to be generated in the obscuring Magdziarz & Zdziarski 1995). In Seyfert 2 , material, e.g. in the molecular torus or even further reflection of the central radiationon the molecular away. torus often leads to an additional, possibly dom- inant, reflection component (Ghisellini, Haardt & Matt 1994). Key words: galaxies: active; galaxies: individual (NGC 4388); gamma rays: observations; X-rays: NGC 4388 (z = 0.00842; Phillips & Malin 1982) is galaxies. a good probe to test these scenarios, as it is one of the brightest Seyfert 2 galaxies at hard X-rays. This barred spiral of morphological class SB(s)b 1. INTRODUCTION pec is located at RA = 12h25'46.7'' and DEC = 12'39'44'' (J2000.0) with an apparent optical magni- tude of B = 11.9mag and is a member of the Seyfert 2 galaxies are among the faintest Active . The energy range up to - 1OkeV Galactic Nuclei (AGN) seen in the Universe. This has been studied with X-ray missions like ROSAT (Antonelli et al. 1997), ASCA (Iwasawa et al. 1997; obs. date rev. SPI exD... Doint. observation Forster et al. 1999), BeppoSAX (Risaliti 2002), and start time [sec] mode lately with Chandra (Iwasawa et al. 2003). Spectral 05/01/03 28 80520 39 5x5 studies of NGC 4388 of the hardest X-ray energies 11/01/03 30 10558 2 staring up to several hundreds of keV have been performed 17/01/03 32 182643 43 5x5 by SIGMA (Paul et al. 1991), by the OSSE experi- 16/06/03 82 20000 1 staring ment (Johnson et al. 1993) on the Compton Gamma 06/07/03 89 152950 47 5x5 Ray Observatory (CGRO; Gehrels et al. 1993), and 09/07/03 90 59989 17 5x5 up to - 150keV by the Phoswich Detector Sys- 18/07/03 93 5283 2 staring tem (PDS; Frontera et al. 1997) on-board Bep- poSAX (Butler & Scarsi 1990; Boella et al. 1997a). Table 1. INTEGRAL observations The two BeppoSAX observations in 1999 and 2000 showed a high energy spectrum with a power law of r = 1.6 and 1.5, respectively. Due to the en- ergy range covered by the PDS it was only possi- ble to give a lower limit for a possible high energy cut-off (Ec > 109 keV). The INTEGRAL mission (Winkler et al. 2003), launched in October 2002, of- fers the unique opportunity to study the entire spec- trum from 3 keV up to several MeV simultaneously. NGC 4388 was the first Seyfert 2 galaxy detected by INTEGRAL in January 2003. In this paper we present data from several INTEGRAL observa- tions and compare the results with previous SIGMA, OSSE, and BeppoSAX measurements. The data also allow us to study spectral variations between the 2003 observations. In addition XMM-Newton data are presented to study the connection with soft X- rays.

Figure 1. INTEGRAL/ISGRI significance contour map an the 15 - 40 keV energy band, based on the JUNO3 data set. The coordinates are equatorial 2. SIMULTANEOUS X-RAY AND (J2000.0). GAMMA-RAY OBSERVATIONS

2.1. IBIS NGC 4388 was detected by INTEGRAL in the course of the 3C273 observation (Courvoisier et al. 2003a). Version 3.0 of ISDC's (Courvoisier et al. 2003b) Of- Though NGC 4388 lies 10.4 degrees north of the fline Science Analysis (OSA) software has been 3C273, valuable data for the Seyfert 2 galaxy partly used for analysing the data. Due to known were obtained due to the large field of view (FOV) problems in this software version, an alternative ap- of IBIS (19" x 19", partially coded FOV) and SPI proach has been used when extracting spectra from (35" x 35", partially coded FOV), and because of the the ISGRI data. The main difference is that we com- dithering observation strategy which is optimised for bined the data of different pointings first before ex- the spectrometer in order to allow a proper back- tracting fluxes in several energy bands. ground determination. Observations were performed during seven INTEGRAL revolutions in January, NGC 4388 was too faint to be seen by IBIS/PICsIT. June and July 2003 (Tab. 1). The total amount of ex- Summing all available IBIS/ISGRI data together posure time is 512 ksec with 274 ksec and 238 ksec in shows that NGC 4388 is detected up to 200 keV, the January and June/ July campaign, respectively. though the data above 150 keV are consistent with Most of the observations were carried out in a dither- zero flux on a la level. A fit to a single power ing mode. 36 ksec of observation were carried out in law gives a photon index of I? = 1.70+::$ and a staring mode, and are therefore only partly useful for flux of 7.8 f 0.5mCrab in the 20 - 40 keV energy SPI analysis, but they allow a more precise flux ex- band. Analysing the data according to the JAN03 traction from the ISGRI data. Due to the smaller and JUNO3 campaign shows that the source was in FOV of the two X-ray monitors (JEM-X) and of a lower flux state in January 2003 with a similar the optical monitor (OMC) on-board INTEGRAL, spectral slo e while the June spectrum is steeper NGC 4388 has not been observed by those instru- (I? = 2.06-,,,,)+82d but has a higher flux. Note however ments. We will discuss the INTEGRAL data mainly that the spectral shape is still consistent with the in terms of the two observation campaigns in Jan- = 1.7 model on a 2a level. For details see Table 2. uary 2003 (JAN03) and June/July 2003 (JUN03). Figure 1 shows the ISGRI significance contour map a-

for the JUN03 observation, based on the 15 - 40 keV Table 3. XMM-Newton observations energy band. obs. date MOS-1 MOS-2 PN start exp. [sec] exp. [sec] exp. [sec] 07/07/02 9850 9937 4540 12/12/02 11667 11667 8291 2.2. SPI 3. XMM-NE WTON OBSERVATIONS In order to reduce the noise in the data reduction process of the SPI data, only those pointings were Two observations with XMM-Newton considered where NGC 4388 was less than 15 de- (Jansen et al. 2001) were performed in July grees off axis. This decreases the total amount of and December 2002 as described in Table 3. useable pointings to 122 with an exposure time of Concerning the connection between the soft and 178 ksec in JAN03 and 134 ksec in JUNO3. OSA hard X-ray emission, the most interesting data 3.0 has been applied, except for some updates. In for our study come from the two MOS cameras the binning of the data SPIHIST 3.1.2 was used (Turner et al. 2001) and from the PN detector in order to avoid a known problem with the ear- (Struder et al. 2001), as they cover the energy lier version. In the case of the image reconstruc- range 0.2 - 1OkeV. The data have been reduced tion program SPIROS (Skinner & Connell 2003), we using the XMM-Newton Science Analysis Software used the more recent version 6.0, applying the mean version 5.4.1 and the model fitting using XSPEC count modulation background model and the max- 11.3 was done simultaneously for the PN and MOS imum likelihood method in order to find the best data. Both observations are well represented by a fitting model. Only single detector events were used Raymond-Smith model (Raymond & Smith 1977) for the analysis and the dead time correction has at energies below - 2.5keV with low abundance been applied for each detector individually, resulting (4% ) and a temperature of about 0.8keV, in a mean dead time fraction per detector of 12.1%. applying only Galactic absorption (2.7 x lo2' cm-2). We used an updated response function which takes At harder X-rays, the spectrum is dominated by into account the results from the in-flight calibra- an absorbed power law. The column density of the tion on the Crab. For details of the SPI analysis warm absorber is 2.8 x 1023cm-2. The power-law procedures see Diehl et al. (2003) and Beckmann index is not well constrained by the XMM-Newton (2002). The SPI spectrum for the combined data of data, so the photon index of r = 1.7 measured JAN03 and JUNO3 is well described b a single power at hard X-rays, has been applied. In addition, a law with photon index r = 1.68-0.35+O"' and a flux of gaussian line at 6.39f::i: keV is apparent in both 9 f 2.6 mCrab in the 20 - 40 keV energy band, which observations. Details of the fit are listed in Table 4. is consistent with the ISGRI results. Analysing the This is consistent with a recent measurement by JAN03 and JUN03 subset it turned out that the Chandra, which showed an Fe line of centroid in source strength is not sufficient for SPI to achieve 6.36 :",:E keV (Iwasawa et al. 2003). The variation a reasonable signal-to-noise for a spectrum. The re- between the July and December data occurs mainly sults are shown in Table 2. The trend in the spectral in the high energy part. This is clearly seen when evolution is opposite to the one detected in the IS- comparing the December data with the model for &L 1..1-. -1-4.- I :- n\ GFU data, but, the large error bars OI?the SPI resu!t_c UC: JUIY U~W(FI~. L). '&Lk iiieit: are iittie or no should be taken into account. changes below - 2.5keV7 the spectrum varied sig- nificantly at higher energies. The iron fluorescence line at 6.39keV was not affected by the variations. In both spectra there was no sign of a further line at 6.9keV or at 7.1 keV, as reported from Chandra 2.3. RXTE data by Iwasawa et al. (2003).

The All Sky Monitor (ASM) on board the Rossi X- 4. X-RAY TO GAMMA-RAY SPECTRUM ray Timing Explorer (RXTE)scans about 80% of the sky every orbit. This offers an unique way to mon- itor the emission of bright X-ray sources like NGC The INTEGRAL data alone do not allow us to re- 4388 in the 1.5 - 12 keV energy range. We extracted construct the complete X-ray spectrum of NGC 4388, fluxes from the RXTE/ASM data base. The fluxes mainly because the statistics of the ISGRI and SPI have been averaged over the same time periods as the data do not constrain existing models. JAN03 and JUN03 INTEGRAL observations in or- der to have comparable results. The weighted mean Comparison with former observations by SIGMA of the JAN03 data was 3.5f 1.2 mCrab, while during (Lebrun et al. 1992), OSSE (Johnson et al. 1994), the JUNO3 campaign the flux was 2.2 f 1.9mCrab. and BeppoSAXIPDS (Risaliti 2002) show a similar Due to low statistics the analysis in the three dif- spectrum at hardest X-rays as observed with IN- ferent ASM bands did not give further information TEGRAL. Measurements by BATSE gave a flux of about the flux variability. f20-100kev = (2.6 f 1.5) x 10-3photonscm-2s-', Table 2. Fit results for a single-power law model Observation Photon F2 G F2 x; (dof1 Index 20-40 keV 40-100 keV 100-200 keV Jan.03 / ISGRI 1.69 z:::: 4.2 f 0.9 7.2 f 3.1 7.0f5.9 0.63 (5) Jan.03 / SPI 1.56 'E::",: 9.7 f 3.6 17.0 f 7.8 19.9 f 6.2 0.78 (2) Jun.03 / ISGRI 2.06 :;::: 11.2 f 1.4 14.1 f 3.1 10.1 f 8.8 1.00 (7) Jun.03 / SPI 1.36 Ti:$ 4.9 f 1.3 22.6 f 8.7 30.9 f 11.0 0.52 (5) summed / ISGRI 1.70 :::$ 7.0 f 0.4 11.7 f 1.4 11.3 f 3.3 1.69 (14) summed / SPI 1.68 Ti::: 7.7 f 2.1 12.2 f 4.8 12.6 f 7.7 1.30 (4) a un-absorbed flux in ergcm-2 s-'

Table 4. Fit results for the XMM-Newton EPIC data Observation T[kT] abundance " fFe - line EW u x; (dof) ~m-~] ph cmF2s-l] [eV] [eV] July 2002 0.83+:::: 0.05:;:;; 2.49 :::;I 6.7 f 1.2 552 67f16 1.04 (559) Dec. 2002 0.817:::: 0.04 E;::: 2.84 z;::: 7.8 f 1.4 205 67f12 1.31 (893)

NGC 4388 XUM-Newton spectrum December 2002 do10 VC~SYSJuly 2002 mooel

~ ~~ L ...... I 05 I 2 5 .0 channel energy (keV)

Figure 2. Comparison of XMM-Newton MOS data from December 2002 with respect to the model for the July 2002 EPIC/MOS data. compared with an INTEGRAL value of (2.7 f0.3) x photons cm-2 s-l. This indicates that the soft gamma-ray spectrum does not vary dramatically as shown in the lightcurves for the 60 keV and 100 keV flux in Fig. 3. The flux measurements are consistent on a 20 level and the spectral shape is conserved during the 13 years of various observations. time [year] In order to extend the spectrum into the energy range below 20 keV we added the two XMM-Newton Figure 3. Lightcurve for NGC 4388 at 60 keV (up- EPIC data sets, described in the previous section. per solid line) and at 100 keV (lower dashed line). Data taken from SIGMA (circle), OSSE (pentagon), The data were fit simultaneously in XSPEC 11.3 BeppoSAX/PDS (hexagon), SPI (square), and IS- applying the same model as for the XMM data GRI (triangle). alone: a Raymond-Smith model for a hot plasma with Galactic absorption, an absorbed power law, plus a gaussian emission line. Figure 4 shows the combined photon spectrum of NGC 4388. The fit, which allowed for flux variations between the various non-simultaneous observations (see Sect. 5), gave a x; = 1.41 for 1431 degrees of freedom. The nor- malisation factors are in the range 2.5 to 4.2 relative to the XMM-Newton July 2002 data. The power I

L

SPI * I-

t I. 20 50 1W 2w t" 1w Figure 6. NGC 4388 combined high energy spectrum, Figure 5. NGC 4388 combined spectrum in E2fE us. including BeppoSAX PDS, CGRO OSSE, SIGMA, E. The data cannot cnfinn or disprove the existence and INTEGRAL ISGRI/SPI. The spectrum is shown of a high energy cut-off. in instrument dependent counts s-l keV-' .

law, dominating the emission above - 2.5 keV, has The fact that the data of the different missions give a photon index of r = l.SS'::,"; and an absorption an acceptable fit using the same normalisation is of NH = 2.73f:::; x cm-2. Figure 5 shows the based on the fact that the error bars in the high en- combined spectrum in E2 versus E. ergy domain are still rather large. But the good cor- fE respondence between the ISGRI and SPI data leads Using a more complicate model with a cut-off power to the assumption that the INTEGRAL calibration law plus reflection from cold material (the so-called is reasonable and that normalisation problems as re- PEXRAV model; Magdziarz & Zdziarski 1995) in- ported e.g. in Courvoisier et al. (2003) are solved, stead of the simple power law, does not improve the allthough the large error bars on the high energy flux fit significantly. Applying a single-power law plus measurements (Fig. 3) still leave open the possibility an unresolved set of lines, represented by a broad of some errors in the calibration. Gaussian line around 1keV, instead of the Raymond- Smith component does not give acceptable fit results The spectral data of NGC 4388 do not drop drastically above 200 keV, but are not of high 2.0). (xt >> enough significance to distinguish whether there is Using only the high energy (> 20 keV) data gives a a cut-off, as expected from Comptonisation models (Petrucci et al. 2001). The XMM-Newton data show single power law with l? = 1.72::::; with = 1.5 x$ variability only in the energy region above - 2.5 keV. (Fig. 6). Also in this energy region alone the more This variability cannot be explained by a change in complex PEXR4V model does not give a signifi- the absorption column density. But using the same cantly better fit. Statistics do not allow us to dis- model parameters of the July observation (Tab. 4) tinguish betwen a single power law a~dB PEXR4V for the December 2002 data and allowing the normal- model. The differences are smaller than the 10 error ization of the single power-law to vary, shows that a bars of the spectra as seen in Fig. 6 and Fig. 5. flux increase by a factor of 3.4 is sufficient to model the December data (x: = 1.1; Fig. 2). The spectral shape of the power law, which describes the hard 5. DISCUSSION X-rays, the Fe Kcu line strength, and the soft X-ray component does not seem to vary.

The combination of the INTEGRAL data with pre- The observed spectrum can be compared to the one vious observations by XMM-Newton, BeppoSAX, derived by Ghisellini, Haardt, & Matt (1994) based CGRO, and SIGMA shows that the spectrum from on Monte-Carlo simulations for a 0.2 keV up to several hundred keV is well repre- with inclination i = 60" - 63", column density of sented by a thermal Raymond-Smith model with 1024cm-2 and a torus with half-opening angle of only Galactic absorption applied, plus a highly 0 = 30". Based on these assumptions, an equiva- absorbed (NH 2: 2.8 x 1023cm-2) power law lent width of the iron fluorescence line of the order like emission (I' = 1.7), and a Gaussian line of EW - lOOeV would be expected, which would to model the iron fluorescence line at 6.39keV. also mean a parallel evolution of the line flux with The results are consistent with previous stud- respect to the continuum. But the Fe Ka does ies of NGC 4388 using subsets of these data not show variability, while the underlying contin- (Risaliti 2002; Lebrun et al. 1992). Also observa- uum varies by a factor of - 4 in the line region tions with Chandra (Iwasawa et al. 2003) and ASCA (see Fig. 7). The line flux of (6.7 f 1.2) phcm-2 s-l (Iwasawa et al. 1997; Forster et al. 1999) show a and (7.8 f 1.4)phcm-2s-' for the July and De- similar spectral behavior. cember measurement is also consistent with the NGC 4388 combined spectrum XMM-Newton EPIC MOS/PN, BeppoSAX/PDS, INTEGRAL SPI/ISGRI, OSSE, and SIGMA I I I '1 I

C

1 10 100 channel energy (keV)

Figure 4. NGC 4388 combined photon spectrum, including XMM-Newton EPIC MOS/PN, BeppoSAX PDS, CGRO OSSE, SIGMA, and INTEGRAL ISGRI/SPI. The high energies are dominated by the highly absorbed component (sufficiently modeled by a single power law) plus the iron fluorescence line, while the lower energies show emission characteristic for a hot plasma with low abundance. XMM data from July 2002 are not shown in this plot for better visibility. I

sented here alone, to give a preference to a thermal or photoionized plasma, the latter one has the advan- tage of avoiding the low abundance in the vicinity (within several kpc) of the AGN core. The Chan- dra data also allowed a space resolved study at soft X-rays, which was not possible with XMM. We there- fore refer to the more simple model, applied to the data presented here, which gives sufficient good fit results.

The INTEGRALIISGRI measurements indicate that a lower flux in the soft gamma-ray region is accom- panied by a harder spectrum, but the errors on the spectral slope measurement are too large to be sta- tistically significant. Since the SPI data at the same 1999 2000 2001 2002 2003 time show the contrary spectral evolution, it is likely time [year] that the spectral slope did not change and that there is only a change in flux (20 - 200 keV) by a factor of Figure 7. Flux in the 2-10 keV band (triangles) 2. and of the 6.4 keV Fe Ka line (circles) measured by - BeppoSAX (1999 and 2000), Chandra (2001), and This variation is not seen in the RXTE data but the XMM-Newton (2002). The underlying continuum is results are still consistent with the ASM values, as highly variable, while the line flux shows no sign$- the error bars are large. cant variation.

one measured by Chandra ((9.3 f 1.9) ph cmd2 s-l) 6. CONCLUSION (Iwasawa et al. 2003). The equivalent width there- fore varied strongly (EW = 190 - 720eV). This im- The complex X-ray spectrum of NGC 4388 is com- plies that the line emitting region is separated from posed out of three major components: the continuum source of the object, probably being even more remote than the molecular torus. How- ever, the iron line can still include contribution from 0 Hot plasma component: Emission below - the torus, but the variability expected to be seen 2.5keV is dominated by emission of a hot (lo7 in this line component is outshined by the stronger, OK), optically thin plasma with low abundance constant component. (about 4 % solar). This component does not seem to be variable. In ROSAT/HFU and Chan- Another difference of the spectrum observed here to dra observations this emission appeared to be an the simulated one by Ghisellini, Haardt, & Matt extended X-ray out to several kpc from (1994) is the fact that the spectrum of NGC 4388 the AGN. Chandra observations imply that this does not seem to have a cut-off in the hard X-rays. component might be due to photoionized gas. This could be caused by a higher temperature of the corona (>> 100kpV) in the cpp of NC,C 4-388.-A- pp- Nnn-thermal emission: Emissicn sh,c-;c -- flection component is not detectable in the X-rays, 2.5 keV follows a highly absorbed (NH - 2.8 x which might be a hint for a non isotropic radiation 1023cm-2) power law with photon index I? = or, generally speaking, a more complex geometry (see 1.7. While the flux of this component varied by e.g. George & Fabian 1991). Apparently there is no a factor -" 4 over the past ten years, the spectral scattering of the hard X-rays away from the line of shape has been constant within the error bars. sight. This indicates that in the case of NGC 4388 we indeed observe the unscattered emission of the 0 Iron line emission: The Fe Kcr fluores- central engine at hard X-rays. Also no decline due cence line at 6.39 keV, consistent with the red- to Klein-Nishina processes appears above 50 keV. shift of the AGN, has a constant flux of - 7.5 ph cm-2 s-l. No line broadening is seen. The constant soft X-ray emission is most likely not linked to the AGN. This is also supported by the fact, The INTEGRAL data do not allow to constrain the that this emission is extended, as has been reported model on the high energy emission, i.e. there is no from ROSAT/HRI measurements (Matt et al. 1994) sign of a Compton reflection component or a cutoff and lately from Chandra (Iwasawa et al. 2003). The at higher energies. low abundance in the extended emission has been seen already in the ASCA observations though with Following the observations reported here, NGC 4388 much larger uncertainties (2 = 0.052,;,52@; Iwa- was studied in a dedicated INTEGRAL observation, sawa et al. 1997). Iwasawa et al. (2003) argue that which was performed in staring mode in revolution the extended emission at low engergies is more likely 94 in July 2003. Though staring mode observations to origin from photoionized plasma. Even though it usually do not produce useful SPI spectra for sources is not possible, based on the XMM-Newton data pre- as faint as NGC 4388, the combination of ISGRI and JEM-X data might show in more detail the connec- Johnson, W. N., Kinzer, R. L., Kurfess, J. D., et al. tion between the hard and soft X-rays, as seen in 1993, ApJS, 86, 693 this work from simultaneous INTEGRAL and RXTE Johnson, W.N., et al. 1994, in The 2nd Compton data. Another 500 ksec observation of 3C273 in AO- Observatory Symposium, C.E. Fichtel et al. (eds.) 2 might give further insights into the cut-off of the New York: AIP, p. 515 NGC 4388 spectrum at high energies, e.g. when com- bined with the already existing data. Lebrun, F., Ballet, J., Paul, J., et al. 1992, A&A 264, L22 Magdziarz, P., & Zdziarski, A. A. 1995, MNRAS, ACKNOWLEDGMENTS 273, 837 Malzac, J., & Petrucci, P.-0. 2002, MNRAS, 336, 1209 This research has made use of the NASA/IPAC Ex- tragalactic Database (NED) which is operated by the Matt, G., Piro, L., Antonelli, L. A., et al. 1994, Jet Propulsion Laboratory. A&A, 292, L13 Parmar, A.N., Martin, D. D. E., Bavdaz, M., et al. 1997, A&AS, 122, 309 REFERENCES Paul, J. et al. 1991, Adv.Sp.Res. 11, 8, 289 Petrucci, P. O., Haardt, F., Maraschi, L., et al. Antonelli, L. A., Matt, G., & Piro, L. 1997, A&A, 2001, ApJ, 556, 716 317,686 Phillips, M. M., & Malin, D. F. 1982, MNRAS, 199, Beckmann, V. 2002, proc. of the XXII Moriond As- 905 trophysics Meeting "The Gamma-Ray Universe", eds. A. Goldwurm, D. Neumann, and J. Tran Raymond, J. C., & Smith, B. W. 1977, ApJS, 35, Thanh Van, p. 417 419 Boella, G., Butler, R. C., Perola, G.C., et al. 1997, Risaliti, G. 2002, A&A, 386, 379 A&AS, 122, 299 Skinner, G., & Connell, P. 2003, A&A, 411, L123 Boella, G., Chiapetti, L., Conti, G., et al. 199713, Struder, L., Briel, U., Dennerl, K., et al. 2001, A&AS, 122, 327 .4&A, 365, L18 Butler, C., Scarsi, L., 1990, SPIE, 1344, 46 Turner, M.J.L., Abbey, A., Arnaud, M., et al. 2001, Courvoisier, T.J.-L., Beckmann, V., Bourban, G., A&A, 365, L27 et al. 2003a, A&A, 411, L343 Winkler, C., Courvoisier, T. J.-L., Di Cocco, G., et Courvoisier, T.J.-L., Walter, R., Beckmann, V., et al. 2003, A&A, 411, L1 al. 2003b, A&A, 411, L53 Diehl, R., Baby, N., Beckmann, V., et al. 2003, A&A, 411, L117 Fabian, A. C. 1977, Nature, 269, 672 Forster, K., Leighly, K. M., & Kay, L. E. 1999, ApJ, 523, 521 Frontera, F., Costa, E., & Dal Fiume, D. 1997, A&AS, 122, 357 Gehrels, N., Chipman, E., & Kniffen, D. A. 1993, A&AS, 97, 5 George, I. M., & Fabian, A. C. 1991, MNRAS, 249, 352 Ghisellini, G., Haardt, F., & Matt, G. 1994, MN- RAS, 267, 743 Haardt, F., Maraschi, L., & Ghisellini, G. 1994, ApJ, 432, L95 Hanson, C. G., Skinner, G. K., Eyles, C. J., & Will- more, A. P. 1990, MNRAS, 242, 262 Iwasawa, K., Fabian, A. C., Ueno, S., et d. 1997, MNRAS, 285, 683 Iwasawa, K., Wilson, A. S., Fabian, A. C., Young, A. J. 2003, MNRAS, 345, 369 Jansen, F., Lumb, D., Altieri, B., et al. 2001, A&A. 365, L1