Mon. Not. R. Astron. Soc. 000, 000–000 (0000) Printed 17 November 2007 (MN LATEX style file v2.2) OH megamasers, starburst and AGN activity in Markarian 231 A. M. S. Richards1⋆, J. H. Knapen2, J. A. Yates3, R. J. Cohen1, J. L. Collett2, M. M. Wright4, M. D. Gray5 and D. Field6 1 Jodrell Bank Observatory, The University of Manchester, Jodrell Bank, Cheshire, SK11 9DL, UK. 2Centre for Astrophysical Research, University of Hertfordshire, Hatfield, AL10 9AB, UK. 3 Department of Physics and Astronomy, University College London, Gower Street, London WC1E 6BT, UK. 4 Department of Physics, University of Bristol, Tyndall Avenue, Bristol, BS8 1TL, UK. 5 Astrophysics Group, The University of Manchester, Sackville Street, P.O. Box 88, Manchester M60 1QD, UK 6 Intsitute of Physics and Astronomy, University of Aarhus, DK-8000, Aarhus C, Denmark. Accepted ... Received ... in original form ... ABSTRACT We present MERLIN observations of OH maser and radio continuum emission within a few hundred pc of the core of the Ultra Luminous IR Galaxy Markarian 231. This is the only known OH megamaser galaxy classed as a Seyfert 1. Maser emission is identified with the 1665- and 1667-MHz transitions over a velocity extent of 720 km s−1. Both lines show a similar position-velocity structure including a gradient of 1.7 km s−1 pc−1 from NW to SE along the 420-pc major axis. The (unresolved) inner few tens of pc possess a much steeper velocity gradient. The maser distribution is modelled as a torus rotating about an axis inclined at ∼ 45◦ to the plane of the −3 sky. We estimate the enclosed mass density to be 320 ± 90 M⊙ pc in a flattened < 6 distribution. This includes a central unresolved mass of ∼ 8 × 10 M⊙. All the maser emission is projected against a region with a radio continuum brightness temperature ≥ 105 K, giving a maser gain of ≤ 2.2. The 1667:1665-MHz line ratio is close to 1.8 (the value predicted for thermal emission) consistent with radiatively pumped, unsaturated masers. This behaviour and the kinematics of the torus suggest that the . arXiv:astro-ph/0509327v1 13 Sep 2005 size of individual masing regions is in the range 0 25 − 4 pc with a covering factor close to unity. There are no very bright compact masers, in contrast to galaxies such as the Seyfert 2 Markarian 273 where the masing torus is viewed nearer edge-on. The comparatively modest maser amplification seen from Markarian 231 is consistent with its classification in the unification scheme for Seyfert galaxies. Most of the radio continuum emission on 50-500 pc scales is probably of starburst origin but the compact peak is 0.4 per cent polarized by a magnetic field running north-south, similar to the jet direction on these scales. There is no close correlation between maser and continuum intensity, suggesting that much of the radio continuum must originate in the foreground and indeed the relative continuum brightness is slightly greater in the direction of the approaching jet. Comparisons with other data show that the jet changes direction close the nucleus and suggest that the sub-kpc disc hosting the masers and starburst activity is severely warped. Key words: masers — galaxies: active — galaxies: kinematics and dynamics — radio continuum: galaxies — radio lines: galaxies — galaxies: individual: Markarian 231 1 INTRODUCTION Ultraluminous infrared galaxies (ULIRGs) were first discov- ered using IRAS (see review by Sanders & Mirabel 1996). ⋆ 12 email: [email protected] They are defined as having luminosities LIR ≥ 10 L⊙,. 2 A. M. S. Richards et al. They are thought to be powered mainly by bursts of extreme flux density of 48 mJy. Kl¨ockner, Baan & Garrett (2003) ob- star formation, with contributions from active nuclei in at tained similar results using the Westerbork array at a reso- least some cases (e.g. Genzel et al. 1998). Seemingly without lution of 14 arcsec, but the European VLBI Network (EVN) exception, and in contrast to lower-luminosity starbursts, at 39 milli-arcsec (mas) resolution only detected just over ULIRGs show evidence for dramatic intergalactic collisions half the OH emission within a region of maximum angular or mergers when observed in the optical or NIR (Sanders & size 150 mas. Mirabel 1996; Knapen 2004). At such extreme luminosities, The Multi Element Radio Linked Interferometer Net- ULIRGS would be detectable out to redshifts ≫ 1 by new work (MERLIN) is sensitive to 1.4–1.6 GHz emission on IR telescopes such as Spitzer. The galaxy interaction rate is scales up to 2 arcsec (1.6 kpc) with a resolution of 100 − 200 expected to be greater at larger look-back times (Schweizer mas. We describe our observations using MERLIN in Sec- 1998). The study of these rare objects in the local universe tion 2. We present the continuum and line results in Sec- is crucial to our understanding of galaxy interactions and tion 3. In Section 4 we explain the radio morphology and the triggering of star-formation at all epochs and, by impli- polarization and analyse the maser kinematics in compari- cation, the development of structure in the early universe son with other observations of the galaxy on various scales, (Kneib et al. 2004). which provides limits on the size of masing clouds and the Sanders et al. (1988) proposed that ULIRGs are part enclosed mass. In Section 5 we analyse our findings and in- of a sequence in which spiral galaxies merge, evolve through formation from the literature in order to estimate the maser a starburst phase and become optically visible QSOs. They optical depths and other characteristics, further constrain- postulate that gas is funnelled into the merging nuclei caus- ing the size of masing regions and their number density. We ing nuclear starbursts and the formation of a self-gravitating make some comparisons with other megamaser galaxies in nuclear disc on a scale of 1 kpc. The disc is thought to fuel Section 6 and summarise our conclusions in Section 7. the activity of the nucleus. Furthermore, starburst activity produces enough fast-evolving stars to give rise to a high rate of mass loss back into the interstellar medium (ISM) leading to sustained fuelling of the central black hole (Nor- man & Scoville 1988). The merger activity produces a warm (∼ 50 K) dusty environment in the nuclear region, providing 2 OBSERVATIONS good conditions for OH megamaser activity (Randell et al. 2.1 Continuum observations 1995). Markarian (Mrk) 231 is one of the most luminous and We observed Mrk 231 at 1658 MHz on 1997 June 20 using 6 best studied IRAS galaxies, with log(LFIR/L⊙) = 12.04. It antennas of MERLIN for 10 hr with a 16-MHz bandwidth. has a redshift of 0.042, yielding an estimated distance of We also retrieved MERLIN archive observations at 1420 and 1 170 Mpc . Optical spectroscopy shows broad lines character- 1658 MHz made on 1993 May 23-25 for ≈ 7 hr at each istic of a Seyfert 1 nucleus (Boksenberg et al. 1977). A radio frequency, using 8 antennas including the Lovell and War- jet is apparent in VLBI and VLA images on scales of sub-pc dle telescopes, giving improved sensitivity on angular scales to hundreds of pc, with a weak counter-jet (Ulvestad, Wrobel up to 6 arcsec. J1302+578 was used as the phase-reference & Carilli 1999 and references therein), providing unusually source, observed alternately with Mrk 231 in a 10-12-min strong evidence for the presence of an AGN in an OH mega- cycle. Its position is known to better than 1 mas (Ma et al. maser region. The IR luminosity has been partly attributed 1998) and it is only 1◦. 3 from Mrk 231. Details of the sources to starburst activity (Downes & Solomon 1998) and partly used for calibration of the bandpass and the flux scale (by to the AGN (Soifer et al. 2000), discussed further in Sec- comparision with with 3C286, Baars et al. 1977) are given < tion 4.1. On scales ∼ 100 pc the radio-continuum emission is in Table 1. The data were reduced using standard MERLIN AGN dominated (Lonsdale et al. 2003). On greater than kpc procedures (Diamond et al. 2002) and the aips package. For scales it is starburst-dominated although the contribution each epoch, we self-calibrated J1302+578 and interpolated from the dissipated radio jet may be significant (Ulvestad the phase and amplitude solutions over the Mrk 231 data. et al. 1999a). The spectral indices of the intermediate-scale We also corrected for polarization leakage and polarization radio continuum suggest a diffuse component powered by angle offset (with respect to 3C286). The pointing position a disc or halo of star formation as well as fragmented jet for Mrk 231 (referred to hereafter as the reference position) emission to the south (Taylor et al. 1999). was 12h 56m14s.2383 +56◦ 52′ 25′′. 210 (J2000). The absolute Mrk 231 shows powerful OH megamaser emission first position uncertainty of our images is about 10 mas, taking detected by Baan (1985). Baan, Haschick & Henckel (1992) into account the errors in the positions of the antennas and reported emission in all 4 ground-state transitions; the 1667- J1302+578 and in transferring the phase corrections to Mrk MHz line was stronger than the 1665-, 1720- and 1612-MHz 231. lines by factors of approximately 2.3, 13 and 30, respec- Table 1 gives the properties of images made using the tively. Using the Lovell telescope with a beam size of 10 same restoring beam at all epochs/frequencies. This shows arcmin, Stavely-Smith et al. (1987) measured emission over that the peak flux densities are the same within the flux scale 6.9 MHz and estimated the total velocity width of the 1667- uncertainy, but the total flux detected in 1997 is only 86 per −1 MHz emission to be 760 ± 100 km s with a peak maser cent of that detected in 1993.