EDP Sciences
Press Release
Free Access
Volume 534, October 2011
Article Number A36
Number of page(s) 10
Section Extragalactic astronomy
DOI https://doi.org/10.1051/0004-6361/201116869
Published online 29 September 2011

© ESO, 2011

1. Introduction

Active galactic nuclei (AGN) are compact sources with very high luminosities, located at the centres of galaxies. Accretion onto the super-massive black holes (SMBH) at their centres is generally believed to be the driving process for the activity. Thanks to their brightness, they form one of the richest laboratories for studying astrophysical processes. In this paper we present one of the deepest multiwavelength campaigns of an AGN, the Seyfert 1 galaxy Mrk 509. It is among the best for these studies because it is unique in combining X-ray brightness, outflow features, and significant but moderate variability. Below we introduce the most important astrophysical processes that are addressed by our study.

After our introduction of the relevant astrophysics, we briefly provide an overview in Sect. 3 of the target of our campaign, Mrk 509, followed by an overview of the observations in Sect. 4, and we present the spectral energy distribution (SED) and light curve in the subsequent sections. Details about various aspects of our study are deferred to subsequent papers of this series.

2. Astrophysics of AGN in the context of our campaign

2.1. Outflows from AGN

In recent years the potential importance of AGN outflows for the growth of SMBH (Silk & Rees 1998; Blandford & Begelman 1999; King 2003; Blandford & Begelman 2004), the enrichment of the intergalactic medium (Furlanetto & Loeb 2001; Cavaliere et al. 2002), the evolution of the host galaxy (Scannapieco & Oh 2004), cluster cooling flows (Wu et al. 2000; Bower et al. 2001; Ciotti & Ostriker 2001; Borgani et al. 2002; Platania et al. 2002), the magnetisation of cluster and galactic gas (Daly & Loeb 1990; Furlanetto & Loeb 2001; Kronberg et al. 2001), and the luminosity function of AGN (Wyithe & Loeb 2003) has been widely recognised. However, for the lack of a better alternative, theoretical studies use the physical properties of the outflow (metallicity, mass, and kinetic energy flux) as free parameters because there are few observational constraints. To assess the importance of AGN outflows on the processes mentioned above, it is essential to establish the real mass flux, metallicity, and kinetic luminosity (Ėk) of AGN winds. Determining Ėk requires answers to some fundamental questions. What is the physical state of the outflowing gas and what is its total column density? What is the distance of the gas from the central source?

The Ėk of a shell-like, non-accelerating outflow is given by (1)where Ω is the solid angle occupied by the outflow, R the distance from the central source, NH the total hydrogen column density, mp the proton mass, fm = ρ/nHmp (1.43 for a plasma with proto-solar abundances) with ρ the mass density and nH the hydrogen density, and v the outflow velocity. Spectral observations straightforwardly determine v, and Ω is expected to be  ~π since 50% of all Seyfert 1s show outflow signatures (Crenshaw et al. 1999). Our campaign focuses on determining the two most uncertain quantities: NH and R. Determining these is tightly connected to the physical state and location of the outflow, as explained below.

In the best studied AGN at least two to three ionisation components are needed to model the rich X-ray absorption spectrum. In NGC 3783 there are three discrete components, which are most likely in pressure equilibrium, representing different phases of gas at the same distance from the nucleus (Krongold et al. 2003; Netzer et al. 2003). However, in NGC 5548 at least five ionisation components are needed if the X-ray absorber is modelled by a finite number of discrete absorption systems, and these absorbers cannot be in pressure equilibrium (Kaastra et al. 2002). Instead, a continuous, power-law distribution of NH versus ionisation parameter ξ gives a better description of the data (Steenbrugge et al. 2003, 2005). In yet another case, Mrk 279, the distribution is continuous but more complicated than a simple power-law (Costantini et al. 2007).

Photo-ionisation modelling of the outflow yields the ionisation parameter ξ ≡ L/nR2, where L is the 1–1000 Ryd ionising luminosity and n the hydrogen density. Independent measurements of n then determine R. The use of density sensitive X-ray lines is difficult and has not (yet) delivered robust results (Kaastra et al. 2004). Alternatively, when L changes, the outflow has to adjust to the new situation. How fast this happens depends on the recombination time scale, which scales as  ~n-1. This method has been applied to a 280 ks XMM-Newton observation of NGC 3783 with wildly different results. From the RGS data R > 10 pc was inferred based on the lack of change in the deep O vii/O viii absorption edges and the Fe-M UTA absorption complex (Behar et al. 2003). On the other hand, variability in the Fe xxv resonance line at 6.7 keV detected in the EPIC data (Reeves et al. 2004) implies R < 0.2 pc. From a 100 ks XMM-Newton observation of NGC 4051, R ≃ 0.001 pc was deduced for the highest ionisation gas (Krongold et al. 2007). However, the large amplitude variations (max/min = 12) and the short time scales (down to 100 s) make these results rather model-dependent; an analysis of Chandra LETGS data yielded distances in the range of 0.02−1 pc (Steenbrugge et al. 2009).

Density-sensitive lines have given more robust results in the UV, and in at least one case (NGC 3783) they yield a distance comparable to the one determined from densities based on recombination time scales. Using metastable transitions in [C iii] λ1176, Gabel et al. (2005) find a distance of  ~25 pc for the absorbjing gas in one of the components in NGC 3783. They find a comparable distance based on the recombination time for Si iv. Again using the [C iii] metastable transitions, gas in one of the components in NGC 4151 lies at a distance of  <0.1 pc (Kraemer et al. 2006). Using metastable levels in [Si ii and [Fe ii], Moe et al. (2009) find a distance for a low-ionisation outflow in the quasar SDSS J0838+2955 of  ~3 kpc. Similarly, large distances for outflows have been found using similar density diagnostics in luminous quasars (Hamann et al. 2001; Dunn et al. 2010; Aoki et al. 2011; Arav et al. 2011).

In order to obtain a reliable density, we need both accurate ionic column densities, which yield the ionisation structure and total outflowing column density, and time variability on a suitable time scale, from which the location of the outflow can be constrained through monitoring. With the current generation of X-ray telescopes, only a handful of AGN with outflows have high enough fluxes to yield suitable data for such an ambitious programme. Of these targets, Mrk 509 is the most promising due to the following three attributes. 1) It is among the brightest X-ray AGNs in the sky. 2) It has an excellent X-ray line structure (Pounds et al. 2001; Smith et al. 2007; Detmers et al. 2010), not as deep and blended as NGC 3783, but not as shallow and sparse as Mrk 279. Mrk 509 is a representative example of a moderate outflow. 3) Mrk 509 varies on a time scale of a few days making it ideal for a spectral timing campaign. Other targets like MCG −6-30-15, Mrk 766, or NGC 4051 with long XMM-Newton exposure times vary too rapidly to do this: within the characteristic time scales of 100−1000 s, it is simply not possible to get a high-quality grating spectrum.

2.2. Kinematics of the outflow as traced in the UV

Our campaign on Mrk 509 includes HST high-resolution ultraviolet spectral observations that are simultaneously with Chandra’s grating X-ray observations. The ultraviolet spectra provide a more detailed view of the kinematics of the outflowing gas. Prior high-resolution spectral observations with the Far Ultraviolet Spectroscopic Explorer (FUSE) (Kriss et al. 2000) and the Space Telescope Imaging Spectrograph (STIS) on the Hubble Space Telescope (HST) (Kraemer et al. 2003) provide a good baseline for our observations, but these were both obtained nearly a decade before to our current campaign. The new UV spectra obtained with the Cosmic Origins Spectrograph (COS) on HST provide an updated view of the outflow, and enable us to examine long-term characteristics of the variability.

Although the UV absorbers are not necessarily a perfect tracer of the X-ray absorbing gas (Kriss et al. 2000; Kraemer et al. 2003), they provide insight into both lower ionisation absorbing components and portions of the high-ionisation outflow. Our COS observations yield accurate measurements of (or limits on) the column densities of C ii, C iv, N v, O i, Si ii – Si iv, S ii, S iii, Fe ii, Fe iii, and H i Lyα. Our limits or detections on the low-ionisation ions of Fe, Si, and S nicely complement the high-ionisation range covered by the RGS and the LETGS. Velocity-resolved measurements of the Li-like doublets of C iv and N v can be used to determine column densities and covering fractions of the outflowing gas (e.g., Arav et al. 2007).

2.3. Abundances in AGN

Studying the abundances, especially those of C, N, O, and Fe, of the gas in galaxies informs us about the ongoing enrichment processes through AGB star winds and supernova type Ia and II explosions. Therefore there has long been an interest in studying abundances with redshift to determine the star formation rate through history. Obvious sources for the study of the abundances are AGN, by taking advantage of their luminosity. Generally, the abundances in AGN have been determined from broad emission lines, however several assumptions make this method subject to systematic errors. Thanks to the improvement in atomic data, especially in the dielectronic recombination rates for iron, the excellent statistics in the combined RGS spectrum, and the simultaneous optical and hard X-ray fluxes to constrain the SED, we can determine reliable and accurate relative abundances from the narrow absorption lines of the outflow observed in Mrk 509.

Since absolute abundances (i.e., metal abundances relative to hydrogen) cannot be measured directly from X-ray spectra, our UV observations of the Lyα line provide a constraint on the total hydrogen column density. The first absolute abundance estimate for an outflow was derived by our group from the UV spectra of Mrk 279, where absorption lines of H, C, N, and O were carefully modelled (Arav et al. 2007) to find abundances relative to solar for carbon (2.2    ±    0.7), nitrogen, (3.5    ±    1.1) and oxygen (1.6    ±    0.8), which fully agree with the relative abundances of C, N, and O derived from the simultaneous X-ray spectra (Costantini et al. 2007). The apparently enhanced N/H and N/O ratios in the Mrk 279 outflow may hint at strong contributions from stellar winds from massive stars and AGB stars.

2.4. Broad emission lines

Broad lines, especially at the energy of the O vii triplet and the O viii Lyα line, have been recently detected in high-resolution X-ray spectra (Kaastra et al. 2002; Ogle et al. 2004; Steenbrugge et al. 2005; Costantini et al. 2007; Longinotti et al. 2010). The width of the lines (about 1 Å, 10 000 km s-1) is comparable to the width of the broad lines detected in the UV band and known to be produced in the so-called broad line region. In at least one case, the soft X-ray emission lines could be related, via a physical model (the locally optimally emitting clouds model, Baldwin et al. 1995), to the UV lines (Costantini et al. 2007). The contribution of the broad line region gas to the iron Kα line at 6.4 keV is not completely understood. For classical Seyfert 1 galaxies this contribution seems modest (e.g.,  <20% in Mrk 279, Costantini et al. 2010), while it may account for most of the line emission in specific objects (e.g., NGC 7213, Bianchi et al. 2008). The quality of the data of Mrk 509 allows us to significantly detect broad emission features around the O vii resonant line, the O viii, N vii Lyα lines and the Ne ix triplet. The simultaneous observation of HST/COS and Chandra/LETGS allow us to connect the broad lines detected in the UV with the soft energy lines and even the iron Kα line.

2.5. Iron K complex

Mrk 509 also shows a rich variety of emission and absorption components in the Fe-K band. In particular, EPIC data from previous XMM-Newton observations shows evidence of a complex Fe-K emission line, with a narrow and neutral component possibly produced far from the source, plus a broad, ionised, and variable component possibly originating in the accretion disk (Pounds et al. 2001; Page et al. 2003; Ponti et al. 2009).

In addition, strong absorption features were found in the same data set at rest-frame energies 8−8.5 keV and 9.7 keV. These were interpreted as being produced by H-like iron K and K-shell absorptions associated with an outflow with mildly relativistic velocity of 0.14−0.2c. The lines were found to be variable in energy and marginal in intensity, implying that variations in either the column density, geometry, and/or ionisation structure of the outflow were maybe common in this source (Dadina et al. 2005; Cappi et al. 2009).

The above properties, combined with the source brightness (F2−10   keV ~ 2−5 × 10-14 W m-2) make Mrk 509 unique in attempting time-resolved spectral studies in the Fe-K band to follow the time evolution of both emission and absorption Fe-K features, and in disentangling the different physical components present in this source.

2.6. Continuum emission of Seyfert galaxies

The physical process at the origin of the X-ray emission of Seyfert galaxies is generally believed to be thermal Comptonisation. In this process, the soft UV photons coming from the accretion disk are up-scattered into the X-ray domain by the hot thermal electrons of the corona. A thermal distribution is generally preferred to a non-thermal one since the discovery, more than 15 years ago, of a high-energy cut-off near 100 keV in NGC 4151 by OSSE and Sigma (Jourdain et al. 1992; Maisack et al. 1993). Such a cut-off was not expected by the non-thermal models developed at that time to explain the power-law-like X-ray spectrum (e.g. Zdziarski et al. 1990). The presence of a high-energy cut-off has now been observed in a large number of objects (Perola et al. 2002; Beckmann et al. 2005; Dadina 2008).

The cold (disk) and hot (corona) phases are expected to be radiatively linked: part of the cold emission, which gives birth to the UV bump, is produced by the reprocessing of part of the hot emission. Inversely the hot emission, at the origin of the broad-band X-ray continuum, is believed to be produced by Compton up-scattering of the soft photons, emitted by the cold phase, on the coronal energetic electrons. The system must then satisfy equilibrium energy balance equations, which depend on geometry and on the ratio of direct heating of the disk to that of the corona (Haardt & Maraschi 1991). In the limiting case of a “passive” disk, the amplification of the Comptonisation process is only fixed by geometry. Therefore, if the corona is in energy balance, its temperature Te and optical depth τ must satisfy a relation that can be computed for different geometries of the disk  +  corona configuration (e.g., Svensson 1996).

Realistic thermal Comptonisation spectra have been computed for more than two decades, and important effects like the anisotropy of the soft photon field have been precisely taken into account (Haardt 1993; Stern et al. 1995; Poutanen & Svensson 1996). These effects appear far from negligible. Noticeably, the spectral shape is significantly different for different disk-corona geometries but also for different viewing angles. Moreover, these effects underline the differences between realistic thermal Comptonisation spectra and the cut-off power-law approximation generally used to mimic them (Petrucci et al. 2001). This may have important consequences in the study of superimposed spectral components like the soft X-ray excess and/or the the iron line and the reflection hump, which require a precise determination of the underlying continuum.

A well known characteristic of thermal Comptonisation spectra, however, is that they are strongly degenerate, i.e., significantly different combinations of temperature and optical depth of the corona give the same power law slope in the 2–10 keV band. To break this degeneracy requires broad-band observations from UV to hard X-rays, the X-ray/γ-ray shape constraining the high-energy cut-off, which is directly linked to the corona temperature. The ratio of the UV to soft X-ray flux constrains the optical depth better. Stronger and less ambiguous constraints on the nature of the coronal plasma can also be obtained from multiwavelength (from UV to hard X-rays/γ-rays) variability studies, as they provide direct insight into the way the emitting particles are heated and cooled. Indeed, variations in the X-ray spectral shape may be produced by intrinsic changes in the hot corona properties (e.g., changes of the heating process efficiency) and/or by variations in the external environment, such as changes in the soft photon flux (and consequently in the coronal cooling) produced by the cold phase (Malzac & Jourdain 2000; Petrucci et al. 2000). For example, thermal Comptonisation models, where hot and cold phases are in radiative equilibrium, predict that the X-ray spectrum of the sources should harden when the energy of the high-energy cut-off increases (e.g., Haardt et al. 1997). Indeed a correlation between Γ and Te has been observed in different objects (Petrucci et al. 2000; Zdziarski & Grandi 2001), thus giving strong support to a thermal nature of the corona particle distribution. The analysis of the one month, simultaneous IUE/RXTE monitoring campaign on NGC 7469, performed in 1996, is also in agreement with thermal Comptonisation emission (Nandra et al. 2000; Petrucci et al. 2004).

The XMM-Newton/INTEGRAL monitoring of Mrk 509 presented here provides an ideal data set to test Comptonisation models and to derive further constraints on the physical parameters and geometry of the source, as well as on the precise shape of spectral components such as the soft-X-ray excess, the outflow, and the reflection component. It should also be noted that a good knowledge of the SED up to 100 keV is very important input for the photo-ionisation modelling of the absorption/emission features produced by these outflows.

3. Mrk 509

Mrk 509 was detected as a Seyfert 1 galaxy with a photographic magnitude mpg = 13 and a size of 10″ (Markarian & Lipovetskij 1973). It has a redshift of 0.034397 (Huchra et al. 1993) and broad (FWHM 4500 km s-1) and narrow optical emission lines (Markarian 1973). Its optical luminosity is high for a Seyfert 1 galaxy, and it puts Mrk 509 close to the limit between Seyfert galaxies and QSOs (Kopylov et al. 1974).

The host galaxy extends approximately east-west (Magnitskaia & Saakian 1976) and has an axial ratio b/a ~ 0.85 (Dahari & De Robertis 1988; Kirhakos & Steiner 1990). Mrk 509 presently shows no nearby neighbours (Magnitskaia & Saakian 1976), the closest spiral galaxy is at a projected distance of 0.3 Mpc (Fuentes-Williams & Stocke 1988), but the asymmetric halo to the south (MacKenty 1990) is direct evidence for strong interactions at some point in the recent past.

In almost any wavelength band that opened its window, Mrk 509 has been one of the first AGN to be studied, thanks to its brightness, and several well-known astrophysicists today have investigated this object in earlier stages of their careers.

In the radio band, Mrk 509 was discovered at 3.9 cm in 1976 (Mingaliev et al. 1978). Spatially resolved measurements at 6 and 20 cm show an asymmetric structure of 1.4″ × 0.5″ in PA 124° (Ulvestad & Wilson 1984; Unger et al. 1987), corresponding to a linear size of about 1 kpc. This core is surrounded by more diffuse emission, visible at 20 cm, with a size comparable to the optical size (10 × 8″, or 7 × 5 kpc) (Singh & Westergaard 1992). The radio emission seems to be preferentially aligned at PA −40° to −60°, coinciding with the preferred direction of the optical polarisation (Singh & Westergaard 1992; Martin et al. 1983).

In the infrared the first observations date back to 1975 (Stein & Weedman 1976; Allen 1976). The spectrum in the 1−10 μm band is approximately a power law (Rieke 1978; Glass et al. 1982; McAlary et al. 1983). Contrary to several other Seyferts, Mrk 509 shows no evidence of strong dust features in infrared spectra (Roche et al. 1984; Moorwood 1986). The variability in the mid-IR is much weaker than at shorter wavelengths, with a possible delay of two to three months (Glass 2004). This suggests that the bulk of the IR emission originates far from the nucleus.

Quantitative measurements of the optical broad and narrow emission lines were first presented by Osterbrock (1977) and de Bruyn & Sargent (1978). Since then, dozens of papers have appeared with refined measurements and models for the line profiles and intensity ratios. The size of the broad emission line region in the Hβ line has been estimated as 80 light days (0.07 pc) based on reverberation studies (Peterson et al. 1998).

Perhaps the most detailed study of interest for our project has been the mapping of the [O iii] λ5007 line by Phillips et al. (1983). They show that the line has two components. There is a rotating, low-ionisation gas disk, coinciding with hot stellar components, which is photo-ionised by the UV radiation of young hot stars. The rotation axis has approximately the same direction as the radio emission and the preferred polarisation angle of the optical emission mentioned before. The second component is expanding, high-ionisation gas, which is photo-ionised by the nucleus. The line profile of that component is consistent with outflowing gas up to velocities of −800 km s-1. It extends out to about 5 kpc from the nucleus. Outflowing gas is seen at opposite sides of the nucleus, indicating that the outflow is approximately face-on. A similar spatial extent is seen in a few other narrow emission lines (Winge et al. 2000). Approximately the same region also shows broad Hβ line emission, most likely caused by scattered nuclear light in the narrow line region (Mediavilla et al. 1998). X-rays from the nucleus travel through these media, and the lowewst ionised X-ray absorbers may show the imprint of this material on the spectrum.

The first UV spectra were presented in 1980, with measurements of the Lyα line flux using the IUE satellite (Wu et al. 1980). Significantly better spectra were subsequently obtained with HST (Crenshaw et al. 1995, 1999; Kraemer et al. 2003) and FUSE (Kriss et al. 2000).

In X-rays, the first detection was made with Ariel V in 1974−1976 in the 2−10 keV band (Cooke et al. 1978). Soon the first (power-law) spectra were fitted (Mushotzky et al. 1980), and the source was detected at higher energies, beyond 40 keV (Dil et al. 1981) and up to about 200 keV (Pietsch et al. 1981; Rothschild et al. 1983). It was one of the first sources where a soft X-ray excess was discovered with HEAO1-A2 (Singh et al. 1985), the first iron line detection dated to 1987 with EXOSAT (Morini et al. 1987), and the reflection component was revealed by Ginga (Pounds et al. 1994).

Variability was discovered first in the optical band (Magnitskaia & Saakian 1976). They found variations of about a magnitude on a time scale of several months, but a comparison with observations 21 years before showed a similar flux level. This is one of the desired properties for the purpose of our campaign: significant variability on suitable time scales, with a limited range of flux levels. Since then, variability has been found in all energy bands from radio (Mingaliev et al. 1978), infrared continuum (Glass 2004), optical polarisation (Martin et al. 1983), broad emission lines (Peterson et al. 1982, 1984; Rosenblatt et al. 1992), UV (Chapman et al. 1985), and X-rays (Dower et al. 1980).

There have been many papers that derive the mass of the central SMBH using a variety of methods. Derived numbers range from 1.4 × 109 M in one of the oldest determinations, based on the correlation of the Hβ intensity with the full-width at zero intensity (Liu 1983) to 1.43 ± 0.12 × 108 M for a modern estimate based on reverberation mapping of optical broad lines (Peterson et al. 2004).

4. Observations

Table 1

Observation log of our Mrk 509 campaign.

We have obtained data from seven different observatories during our campaign. At the core of our programme are ten observations of approximately 60 ks each, with XMM-Newton, spaced by four days. We used the data of all instruments on XMM-Newton: Reflection Grating Spectrometer (RGS), EPIC (pn and MOS), and Optical Monitor (OM). For this last instrument, all filters except for the white light filter were used, and we also obtained a spectrum with the optical grism for each observation. Observations with the UV grism were not allowed for operational reasons.

Simultaneously with our ten XMM-Newton observations, we obtained data with INTEGRAL to observe the hard X-rays. Three weeks after the end of the XMM-Newton monitoring, we did simultaneous observations with Chandra (Low Energy Transmission Grating (LETGS) with the HRC-S camera), together with the HST COS. Due to operational constraints, these observations could not coincide with the XMM-Newton/INTEGRAL monitoring.

Before our XMM-Newton/INTEGRAL observations, and in between the XMM-Newton and Chandra/HST observations, we monitored Mrk 509 with Swift, using both its X-ray telescope (XRT) and the UltraViolet and Optical Telescope (UVOT). This was in order to obtain a continuous monitoring and to allow us to study the flux before the start of our campaign. This is important because absorption components may respond with a delay to continuum variations.

thumbnail Fig. 1

Timeline of our monitoring campaign of Mrk 509. The first observation, with Swift, started on September 4, 2009, and the last observation, with Chandra, ended on December 13, 2009.

Open with DEXTER

Finally, throughout our high-energy campaign, we obtained optical spectroscopy and photometry. We have five observations with the 4.2 m William Herschel Telescope (WHT) at La Palma, photometry with the ACAM camera using the Sloan g, r, i and Z filters, and low-resolution spectroscopy with the VPH disperser. Furthermore, we obtained nine observations with the 1.3 m PAIRITEL (Peters Automated IR Imaging TELescope), with photometry in the J, H, and K bands. Table 1 gives some details on our observations, and in Fig. 1 we show a graphical overview of the timeline of our campaign.

5. Light curve

The light curves in the X-ray band (0.3–1.0 keV, 2–10 keV and 20–60 keV) and in UV (2310 Å) are shown in Fig. 2. More details and the interpretation of these light curves are given in subsequent papers of our series (e.g. Mehdipour et al. 2011; Petrucci et al., in prep.). The flux was at a typical flux level, compared with archival observations. Right in the middle of our

thumbnail Fig. 2

Top panel: UV light curve of Mrk 509 as obtained with Swift/UVOT, XMM-Newton/OM, and HST/COS. The data have been corrected for extinction. Bottom panel: X-ray light curve as obtained by Swift/XRT, XMM-Newton/pn, Chandra/LETGS, and INTEGRAL/ISGRI. These data have also been corrected for Galactic absorption. For the pn data, the squares correspond to the 0.3–1.0 keV band, the stars to the 2–10 keV band. Note that 1012 JyHz equals 10-14 W m-2 (10-11 in c.g.s. units).

Open with DEXTER

observing campaign, the source showed a significant brightening by  ~50%, in particular in the soft X-ray band, and with a smaller amplitude in the UV and hard X-ray bands. This allowed us to study the response of the absorption components to this brightening, and thanks to our broad-band measurements we can constrain the emission mechanisms for the continuum components.

6. Spectral energy distribution

The broad-band SED is essential for obtaining the ionisation balance needed for the photo-ionisation modelling of the outflow. The multiwavelength nature of our campaign gave us almost simultaneous coverage of the total spectrum, which allowed us to constrain the exact shape of the SED to high accuracy.

We need two different SEDs for our campaign. The first describes the average SED during the XMM-Newton part of the campaign, which is needed for the analysis of the stacked RGS spectrum. The second one is needed for the Chandra LETGS and HST COS part of the campaign, when the source had a lower flux level. Table 2 shows the flux points of every instrument used to create the SED shown in Fig. 3. The statistical uncertainties on the flux are much lower than the reported values (except in the case of INTEGRAL), therefore we list the systematic uncertainties due to calibration uncertainties.

Table 2

Continuum fluxes corrected for Galactic and intrinsic absorption during the XMM-Newton and Chandra parts of the campaign.

thumbnail Fig. 3

SED used of Mrk 509 for the time-averaged XMM-Newton and Chandra observations. Data points from Table 2 are indicated with triangles and squares.

Open with DEXTER

6.1. X-ray SED

We used the best fit RGS model (model 2 of Detmers et al. 2011) for the spectrum between 7−38 Å. We correct for both the intrinsic absorption (the ionised outflow) and the Galactic ISM absorption. For the flux between 1.2 Å (10 keV) and 7 Å (1.77 keV), we used the EPIC-pn data. Above 10 keV we used the INTEGRAL data, which were obtained simultaneous with the XMM-Newton observations. This gives us the continuum flux up to  ~200 keV. Beyond that, the SED is extrapolated using the model described in Petrucci et al. (in prep.).

We compared the Chandra LETGS flux in two different energy bands with the non-contemporaneous RGS flux to obtain the flux variations in the soft (20−35 Å) and hard (7−10 Å) X-ray bands. The RGS flux in these bands is 66% and 30% higher than the LETGS flux, respectively. The soft X-ray flux for the RGS observation for λ > 37 Å is obtained by scaling the LETGS continuum by a factor of 1.66. As we have no information on the X-ray flux above 10 keV during the LETGS observation, and the uncertainty on the LETGS spectrum for λ < 7 Å increases, we estimate the flux for the Chandra observation for λ < 7 Å by dividing our model for XMM-Newton/INTEGRAL by a factor of 1.30.

6.2. EUV extrapolation

The EUV spectrum produces most of the ionising flux; however, it is also the most uncertain part of the SED, as there are no data between 50 and 912 Å. LETGS formally measures up to 175 Å, but due to the Galactic absorption, the flux is low and also the modelling of the higher spectral orders becomes more uncertain at longer wavelengths. We therefore have to interpolate our SED between the soft X-rays and the UV. There are several options for doing this. Our baseline model is a power-law interpolation between 50−943 Å (the last FUSE data point). For the XMM-Newton epoch this power law has a photon index of 2.39, for the Chandra epoch 2.53. Alternatively, we consider some kind of big blue bump, by extrapolating the RGS spectrum with a power law, using its slope near 30 Å (photon index 2.58) up to the point (193 Å) where it matches the extrapolation from the UV with a photon index of −2. The different extrapolations between the X-ray and UV data are shown in Fig. 3.

6.3. UV SED

The UV part of the SED was derived from archival FUSE data and from the HST COS data of our campaign. As FUSE and COS share a common wavelength band, the FUSE archival fluxes were scaled to the COS flux level, assuming that the spectral shape remained the same.

The full description of the COS data reduction can be found in Kriss et al. (2011). In short, the COS data were flat-field corrected and additional wavelength calibration performed. Time-dependent sensitivity corrections were applied to the COS data, resulting in an absolute flux accuracy of 5%. As the COS data were taken during the LETGS observations, we needed to adjust them to the higher flux levels during the XMM-Newton observations. The OM data at somewhat longer wavelengths (see below) show on average 6% higher flux during the XMM-Newton observations compared to the LETGS observations. Using the near-simultaneous FUSE/HST/optical spectrum of Mrk 509 from Shang et al. (2005) normalised to the OM flux, we find that increasing the COS flux by 10% is appropriate if assuming a constant shape for the UV spectrum.

6.4. Optical SED

The optical part of the SED was obtained from the OM data. The data were corrected for interstellar absorption and de-reddened using the reddening curve of Cardelli et al. (1989), including the near-UV update by O’Donnell (1994). The total colour excess E(B − V) is 0.057 mag, and RV = A(V)/E(B − V) was fixed at 3.1. The total neutral hydrogen column density in the direction of Mrk 509 is 4.44 × 1024 m-2, as given by Murphy et al. (1995). The host galaxy correction is based on the results of Bentz et al. (2009) and Kinney et al. (1996). For Mrk 509 the host-galaxy contribution at 2231 Å is 7 × 10-20 W m-2  Å-1, which is negligible compared to the AGN flux of 8.45 × 10-17 W m-2 Å-1. At longer wavelengths the host galaxy contribution increases, to 3.2 × 10-18 W m-2 Å-1 at 5500 Å, where the AGN flux is 1.7 × 10-17 W m-2 Å-1. For a more thorough description of the optical data reduction, see Mehdipour et al. (2011). The optical fluxes (from Swift UVOT) at the time of the LETGS observations have to be adjusted for the lower flux level of the source. As XMM-Newton observation 2 has a similar flux to that of Swift UVOT observation 18, which is close in time to the LETGS observations, we assume that the optical flux during the LETGS observations is the same as that of XMM-Newton OM observation 2, which is 6% smaller than the average optical flux during the XMM observations.

6.5. Infrared and radio SED

At IR wavelengths shortward of 1 μm, the spectrum shows a small upturn due to emission from the torus (see, e.g., the spectra by Landt et al. 2011). The infrared fluxes that we use here are based on photometry with host-galaxy subtraction in the J, H, and K bands. We use the average from the observations of Danese et al. (1992) and Kotilainen et al. (1992), which are in good agreement with each other. For 100 μm, we use the IRAS flux (Moshir et al. 1990). Between the K-band and the 100 μm flux points, the SED shows an almost constant value in terms of νFν. In the 6−35 μm band Spitzer data exist (Wu et al. 2009), with  ~30% higher flux compared to our adopted SED. However, these data, taken with a relatively large aperture of 10″, are not corrected for the stellar contribution. As our results are not very sensitive to the details of the IR spectrum, we can safely use our simplified SED in the mid-IR band as shown in Fig. 3.

Beyond 100 μm, the flux drops rapidly; following usual practice in this region, we extrapolate the flux from 100 μm to lower frequencies with an energy index of −2.5 down to 1 cm. This gives a flux of only about a factor of two above the observed value at 14.9 GHz (Barvainis et al. 1996). Finally, we also consider an SED without IR emission, because it is unclear how much of the IR emission is seen by the outflow (see next section).

6.6. Effect of the SED on the ionisation balance

Because a different SED can have a strong effect on the ionisation balance of the outflow, we investigated the effects of different SEDs, focusing on the uncertainties in the EUV and infrared parts of the spectrum. In the EUV band we have no measurements, so we have to rely on interpolation. Depending on the location of the outflow relative to the dusty torus expected to surround the AGN, it may receive either all or only a small fraction of the IR flux emitted by the torus. To investigate the effects of different assumptions about this, we compare four SEDs (Fig. 3):

  • 1.

    Our standard SED (Table 2,Col. 3 and Fig. 2, solid line), whichuses a simple power-law interpolation between the soft X-ray andUV bands;

  • 2.

    A SED with a stronger EUV flux, using the broken power-law approximation with a break at 193 Å mentioned in Sect. 6.2, to mimic a stronger EUV flux (EUV bump in Fig. 3);

  • 3.

    Same as model 1, but with the infrared flux essentially set to zero, for an absorber that does not receive emission from the torus (No IR in Fig. 3);

  • 4.

    the SED for the Chandra observation (Table 2, Col. 4), to see the effects of time variability.

The ionisation balance calculations were performed using Cloudy1 (Ferland et al. 1998) version C08.00, with Lodders et al. (2009) abundances. The results are shown in Fig. 4. As can be seen, the differences between the four cooling curves are small. The case without infrared flux (model 3) starts deviating from the default case above 106 K, because of enhanced cooling by inverse Compton scattering of the infrared photons. The case with a stronger EUV flux (model 2) has a similar shape to model 1, but has on average log Ξ higher by 0.08 due to the enhanced ionising EUV flux. This results in a small shift towards the right of the figure. Finally, the lower flux during the Chandra observations leads to a very similar cooling curve, because the peak of the UV spectrum is very similar.

We are therefore confident that our derived SED is an accurate description of the true source continuum of Mrk 509 and that the photo-ionisation results obtained from the analysis of Detmers et al. (2011) are not seriously affected by the uncertainties in the assumed SED.

thumbnail Fig. 4

The four different cooling curves for the different assumed SEDs.

Open with DEXTER

thumbnail Fig. 5

The two different cooling curves for solar (solid line) and 3 times solar metallicity (dashed line).

Open with DEXTER

As an additional test we checked how metallicity influences the cooling curve. We did this because we do not know a priori the metallicity of the outflow, which can have a significant effect on the shape of the cooling curve (see e.g., Chakravorty et al. 2009). We adopted 3 times the solar abundances, cf. typical abundances found in Mrk 279 (Arav et al. 2007). Figure 5 shows the results compared to the solar metallicity run. There are slight differences in the shape of the two cooling curves, with the higher metallicity curve having bigger unstable branches (where the slope is negative) due to enhanced line cooling. However, the resulting photo-ionisation model will give no significant differences in derived parameters, except for the hydrogen column densities, which are a factor of three lower owing to the three times higher metal abundances. This again shows that the ionisation structure of the outflow in Mrk 509 is very stable to the uncertainties that are present in either the assumed SED or metallicity.

7. Summary

In this paper we have introduced our multiwavelength study of the Seyfert 1 galaxy Mrk 509, which used XMM-Newton, INTEGRAL, Chandra, HST, Swift, WHT and PAIRITEL. Our observations spanned 100 days of monitoring from September 2009 to December, covering band-passes from 2 μm to 200 keV. The core of our programme consisted of ten simultaneous observations with XMM-Newton and INTEGRAL, followed by a long Chandra LETGS spectrum obtained simultaneously with an HST/COS far-UV spectrum. The high-resolution spectroscopy combined with the time variability monitoring enabled us to disentangle the different absorption and emission components in Mrk 509.

Using our data we produced a continuous light curve in the X-ray and UV bands and a comprehensive SED. Thanks to our broad-band coverage, we deduced that uncertainties on the SED have little effect on the photo-ionisation equilibrium that applies to our subsequent models of the ionised outflow. Our light curve shows a 60% flux increase in the soft X-ray band correlated with an enhancement of the UV flux.

A series of subsequent papers will elaborate on these results and others from our campaign. The stacked XMM-Newton RGS spectrum is presented in Kaastra et al. (2011, Paper II). The time-averaged ionisation and velocity structure of the outflow deduced from this spectrum is presented by Detmers et al. (2011, Paper III), the broad-band continuum variability by Mehdipour et al. (2011, Paper IV). Ebrero et al. (2011, Paper V) describes the Chandra LETGS data, and Kriss et al. (2011, Paper VI) presents the analysis of the HST/COS data. More papers are in preparation.


This work is based on observations obtained with XMM-Newton, an ESA science mission with instruments and contributions directly funded by ESA Member States and the USA (NASA). It is also based on observations with INTEGRAL, an ESA project with instrument and science data centre funded by ESA member states (especially the PI countries: Denmark, France, Germany, Italy, Switzerland, Spain), Czech Republic, and Poland and with the participation of Russia and the USA. This work made use of data supplied by the UK Swift Science Data Centre at the University if Leicester. SRON is supported financially by NWO, the Netherlands Organization for Scientific Research. J.S. Kaastra thanks the PI of Swift, Neil Gehrels, for approving the TOO observations, and the duty scientists at the William Herschel Telescope for performing the service observations. P.-O. Petrucci acknowledges financial support from CNES and the French GDR PCHE. M. Cappi, M. Dadina, S. Bianchi, and G. Ponti acknowledge financial support from contract ASI-INAF n. I/088/06/0. N. Arav and G. Kriss gratefully acknowledge support from NASA/XMM-Newton Guest Investigator grant NNX09AR01G. Support for HST Program number 12022 was provided by NASA through grants from the Space Telescope Science Institute, which is operated by the Association of Universities for Research in Astronomy, Inc., under NASA contract NAS5-26555. E. Behar was supported by a grant from the ISF. A. Blustin acknowledges the support of a STFC Postdoctoral Fellowship. P. Lubiński has been supported by the Polish MNiSW grants NN203065933 and 362/1/N-INTEGRAL/2008/09/0. M. Mehdipour acknowledges the support of a PhD studentship awarded by the UK Science & Technology Facilities Council (STFC). G. Ponti acknowledges support via an EU Marie Curie Intra-European Fellowship under contract No. FP7-PEOPLE-2009-IEF-254279. K. Steenbrugge acknowledges the support of Comité Mixto ESO – Gobierno de Chile.


All Tables

Table 1

Observation log of our Mrk 509 campaign.

Table 2

Continuum fluxes corrected for Galactic and intrinsic absorption during the XMM-Newton and Chandra parts of the campaign.

All Figures

thumbnail Fig. 1

Timeline of our monitoring campaign of Mrk 509. The first observation, with Swift, started on September 4, 2009, and the last observation, with Chandra, ended on December 13, 2009.

Open with DEXTER
In the text
thumbnail Fig. 2

Top panel: UV light curve of Mrk 509 as obtained with Swift/UVOT, XMM-Newton/OM, and HST/COS. The data have been corrected for extinction. Bottom panel: X-ray light curve as obtained by Swift/XRT, XMM-Newton/pn, Chandra/LETGS, and INTEGRAL/ISGRI. These data have also been corrected for Galactic absorption. For the pn data, the squares correspond to the 0.3–1.0 keV band, the stars to the 2–10 keV band. Note that 1012 JyHz equals 10-14 W m-2 (10-11 in c.g.s. units).

Open with DEXTER
In the text
thumbnail Fig. 3

SED used of Mrk 509 for the time-averaged XMM-Newton and Chandra observations. Data points from Table 2 are indicated with triangles and squares.

Open with DEXTER
In the text
thumbnail Fig. 4

The four different cooling curves for the different assumed SEDs.

Open with DEXTER
In the text
thumbnail Fig. 5

The two different cooling curves for solar (solid line) and 3 times solar metallicity (dashed line).

Open with DEXTER
In the text

Current usage metrics show cumulative count of Article Views (full-text article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.

Data correspond to usage on the plateform after 2015. The current usage metrics is available 48-96 hours after online publication and is updated daily on week days.

Initial download of the metrics may take a while.