arXiv:1302.4436v2 [astro-ph.CO] 12 Jun 2013 .F S. abig B H,UK 0HA, CB3 Cambridge odnS72Z UK 2AZ, SW7 London od uhmD13E UK 3LE, DH1 Durham Road, ne rso S T,UK 1TL, BS8 Bristol enue, lcfr il dnug H H,UK 3HJ, EH9 Edinburgh Hill, Blackford oni,RylOsraoy lcfr il dnug EH9 Edinburgh Hill, Blackford Observatory, Royal Council, ugl7,511Bn,Germany Bonn, 53121 71, Huegel e esy ictwy J08854-8019 NJ Piscataway, Jersey, New isUiest,30 .CalsSre,Blioe D2121 MD Baltimore, Street, Charles N. 3400 University, kins iyPr,Ntiga G R,UK 2RD, NG7 Nottingham Park, sity taa Y14853-6801 NY Ithaca, vneCalsAdé -96 an ei aa ee,Fran Cedex, Laval Genis Saint F-69561 André, Charles Avenue rde A02138 MA bridge, 80 an-atndHrs France d’Hères, Saint-Martin 38406 92697 .P K J.-P. rpittpstuigL using typeset Preprint S .J I J. R. HERSCHEL .K M. .G S G. J. BITDTO UBMITTED 21 20 4 19 3 18 2 12 1 17 5 16 14 15 13 6 11 10 9 8 7 LEUREN eto srnm,Uiest fMrln,CleePr,MD Park, College Maryland, of University Astronomy, of Dept nttt o opttoa omlg,Dra University Durham Cosmology, Computational for Institute nttt o srnm,Uiest fEibrh oa Ob Royal Edinburgh, of University Astronomy, for Institute KAtooyTcnlg ete cec n ehooyFac Technology and Science Centre, Technology Astronomy UK avr-mtsna etrfrAtohsc,6 adnSt Garden 60 Astrophysics, for Center Harvard-Smithsonian eto hsc srnm,Uiest fClfri,Irvi California, of University Astronomy, & Physics of Dept twr bevtr,Uiest fAioa usn Z857 AZ Tucson, Arizona, of University Observatory, Steward eto srnm,Uiest fMneoa inaoi,M Minneapolis, Minnesota, of University Astronomy, of Dept ntttd aiatooi ilmtiu,30red aPi la de rue 300 Millimétrique, Radioastronomie de Institut nttt fAtooy nvriyo abig,Madingle Cambridge, of University Astronomy, of Institute nvriä in ntttfrAtohsk -10We,A Wien, A-1180 Astrophysik, für Institut Wien, Universität hsc et meilCleeLno,SuhKnigo Ca Kensington South London, College Imperial Dept, Physics .H il hsc aoaoy nvriyo rso,Tynd Bristol, of University Laboratory, Physics Wills H. H. eto srnm,SaeSineBidn,CrelUnivers Cornell Building, Science Space Astronomy, of Dept readrIsiuefrAtooy onUiest,Auf University, Bonn Astronomy, for Argelander-Institute eto hsc n srnm,Rtes h tt Universit State The Rutgers, Astronomy, and Physics of Dept h er .RwadDp fPyis&Atooy on Hop- Johns Astronomy, & Physics of Dept Rowland A. Henry The eto hsclSine,TeOe nvriy itnKeyn Milton University, Open The Sciences, Physical of Dept colo hsc srnm,Uiest fNtiga,Un Nottingham, of University Astronomy, & Physics of School eted ehrh srpyiu eLo,Uiest Lyon Université Lyon, de Astrophysique Recherche de Centre NF sevtroAtooiod aoa -52 Padova, I-35122 Padova, di Astronomico Osservatorio INAF, .L C L. D. RIPS VISON NEIB otnu n CO and continuum itdt edt h otitnesabrt,wihwl the will which starbursts, intense most the veloci to rotational lead large to dicted to partly and galaxies constituent e aeldt h icvr fa es oritisclyl intrinsically four least at of discovery the to fro led HyLIRG have candidate vey best the of observations Panchromatic × aegsadsa-omto ufc este.Tetobrig two The densities. surface star-formation and sub-arc Via gas proto-ellipticals. rate starbursting of cluster a – fcece r ifiutt eocl ihasml volumet simple a with g reconcile large to comprises difficult supposedly are which efficiencies sequence, main the above iuain ugs hssse stepoeio faB(ina a of progenitor the is system this suggest simulations ofie oternce,a ls oteEdntnlmt The limit. Eddington the to close are at galaxies nuclei, these their radii: to their confined with commensurate scales on Keywords: TAGUHN 10 7 26 .N D. , 32 11 1,2 LEMENTS .F A. , .M S. , ALS IAYHLR IPITN LSE FSTARBURSTI OF CLUSTER A PINPOINTING HYLIRG BINARY A -ATLAS: M T .M S M. A. , HE 14 ⊙ A ARAYANAN .A T A. T. , T n a rcin of fractions gas and , E RANCESCHINI A tl mltajv 12/16/11 v. emulateapj style X ai otnu:glxe ai ie:galaxies lines: radio — galaxies continuum: radio ADDOX aais ihrdhf aais trus submillim — starburst galaxies: — high-redshift galaxies: STROPHYSICAL 19 .D H. , WINBANK ARGETT 23 J 9 .O A. , =4 .N M. , ANNERBAUER − 27 ,addul hti CO in that double and 3, 3 2 .E G E. J. , J .A S. , I , MONT OURNAL .V I. EGRELLO AN ALTCHANOV S MBER ∼ 33 MAIL EACH .S D. , 0.Vlct iprin r oets h ik r unstab are disks the so modest are dispersions Velocity 40%. 20 10 .D A. , 15 3 .N R. , 28 umte oTeAtohsclJournal Astrophysical The to Submitted .J B J. A. , .I H I. A. , COTT .D G D. R. , H,UK 3HJ, 8 35 ARIUSH ce ERI servatory, .V J. , ustria 34 55414 N 21 Road, y South , Cam- , l Av- all Italy .S S. , s UK es, eCA ne AKER 20742 ARRIS mpus, 7 scine, ilities .P J. , ABSTRACT iver- of y dem ,9 1, EORGE IEIRA ity, 21 J ERJEANT .D G. , =1 16 4 EÑARRUBIA .S B S. R. , eoditreoerciaigw aemaue accu- measured have we imaging interferometric second .B N. , 36 − isi w one-oaiggsdss–aseai pre- scenario a – disks gas counter-rotating two in ties 2 .TebodC iei u atyt h uttd of multitude the to partly due is line CO broad The 0. , .C H C. J. , uminous h Rlmnst smaue across measured is luminosity IR the aoaoiswt hc ocnrn h otrcn hydro- recent most the confront to excellent provide which They with extreme. laboratories most and rarest the amongst (HyLIRGs, galaxies (IR) infrared nous E AND eoecm npis h ik have disks The pairs. in come refore ts aaisspan galaxies htest ry)-type Z C03,90DiledDie lcsug A24061 VA Blacksburg, Drive, Drillfield 910 0435, MC ns h aae5 adf F43A UK 3AA, CF24 Cardiff 5, Parade The ings, edA1 A,UK 9AB, AL10 field a 80 hithrh e Zealand New Christchurch, 4800, Bag Italy ieEdRa,Lno 14S UK 4NS, E1 London Road, End Mile rfsc,Vcñ akna46,Csla36 atao22, Santiago 306, Casilla 4860, Mackenna Vicuña trofísica, Switzerland edn h Netherlands The Leiden, aaea A91125 CA Pasadena, 89 ilnead aCñd,Spain Cañada, la de Villanueva 28691 ovrVT11 Canada 1Z1, V6T couver OURNE ftekondnzn fteglx o,hyperlumi- zoo, galaxy the of denizons known the Of 15 negigeteebrt fsa omto,not formation, star of bursts extreme undergoing i trfrainlaw. star-formation ric h widest the m OTTI rseicsa-omto ae lc them place rates star-formation specific ir 26 25 29 24 28 23 22 27 32 31 30 34 33 37 36 35 USSMANN .W .S L. W. M. , P. sdsslk hs.Terhg star-formation high Their these. like disks as colo hsc srnm,QenMr nvriyo Lond of University Mary Queen Astronomy, & Physics of School eto hsc,Vrii oyehi nttt tt Uni State & Institute Polytechnic Virginia Physics, of Dept nvria aóiad hl,Dpraet eAstronomía de Departamento Chile, de Católica Universidad colo hsc srnm,CrifUiest,Queen’s University, Cardiff Astronomy, & Physics of School etefrAtohsc eerh nvriyo Hertfords of University Research, Astrophysics for Centre eto hsc n srnm,Uiest fCnebr,Pr Canterbury, of University Astronomy, and Physics of Italy Trieste, Dept I-34136 265, Bonomea via SISSA, iatmnod srnma nvria iPdv,312P 35122 Padova, di Universita’ Astronomia, di Dipartimento PLS PPLSR,Osraor avry H19 Versoi CH-1290 Sauverny, Observatoire LASTRO, IPEP SB EPFL eto hsc,Uiest fteWsenCp,SuhAfric South Cape, Western the of University Physics, of Dept srpyis eto hsc,KbeRa,Ofr X R,U 3RH, OX1 Oxford Road, Keble Physics, of Dept Astrophysics, eto hsc srnm,Uiest fBiihColumbia British of University Astronomy, France & Paris, Physics F75014 IAP, of UMR7095 Dept CNRS & 06 Paris Univ. UPMC ednOsraoy ednUiest,PO o 53 NL-2 9513, Box P.O. University, Leiden Observatory, Leiden aionaIsiueo ehooy 20Es aionaB California East 1200 Technology, of Institute California uoenSaeAtooyCentre, Astronomy Space European ELLY 2 A DER VAN .R J. , 22,10 13 z ≈ 2 = .B F. , .D L. , 3 10 .H R. , ICHARD tr aais—ifae:glxe — galaxies infrared: — galaxies eter: 5 . 1glxe cosa across galaxies 41 .C A. , 14 MITH ERTOLDI W . Herschel 6 UNNE -M OPWOOD ERF 1. ∼ 11 OORAY 24 ⊙ p WMi submm/radio in FWHM kpc 3 .A R A. D. , 37 .A T A. M. , INTRODUCTION 23 cluster. 17 .A E A. S. , .B M. , xrglci mgn sur- imaging extragalactic 19,15 6 N .C P. , bd n akmatter dark and -body GPROTO-ELLIPTICALS NG IECHERS .I E. , HOMPSON REMER ALES OX ≈ BAR Herschel 7 0-p region 100-kpc M .F H. , 24 λ e potentially le, 12 18 dyn .F D. , L 1,29 rest .R K. , .A C A. J. , IR 28 fseveral of .J J J. M. , U .V E. , 8–1000 = ≥ cec ete ESA, Centre, Science 6 ARRAH .K A. , 10 OWLANDS ∼ > ALIANTE 13 ALANOG 5 Chile ARVIS × L OVÁCS 25 ⊙ µ ie Hat- hire, , where , oulevard, )are m) 0 RA 300 a versity, 30,31 adova, Build- 13 Van- , As- y 24 6 8 ivate , , , K , on, x, , 2 IVISONETAL. dynamic simulations of isolated and merging galaxies (e.g. 2. OBSERVATIONS AND DATA REDUCTION Hayward et al.2011). The LIR of a HyLIRG implies a stagger- −1 Herschel imaging, undertaken as part of the wide-field H- ing star-formation rate, SFR ∼> 1000M⊙ yr (Kennicutt 1998, ATLAS imaging survey (Eales et al. 2010), led to the se- for a Chabrier 2003 initial mass function, IMF), unless the lection of HATLASJ084933 as a potentially distant, lensed IMF is top-heavy, or there is a substantial contribution to LIR starburst: distant, because it is an example of a so-called from a deeply obscured AGN, in which case we are seeing ‘350-µm peaker’, with its thermal dust peak between the intense star formation accompanied by the rapid growth of a 250- and 500-µm bands; lensed, because models of the far- massive black hole (Alexander et al. 2008). Either way, we IR/submm source counts suggest that the majority of objects are witnessing galaxy formation at its most extreme. with S500 > 100mJy are expected to be either local (z < 0.1) How best can we identify the most luminous star-forming or lensed (e.g. Negrello et al. 2010), along with the occasional HyLIRGs? A promising method involves searching amongst flat-spectrum radio quasar. the brightest sources detected in the widest surveys with Her- On the basis that it represented the best chanceof finding an schel. These comprise blazars, low-redshift spirals and some intrinsically luminous system rather than a lensed galaxy, fol- intrinsically fainter sources that have been strongly lensed lowing the arguments laid out in §1, HATLASJ084933 was (e.g. Negrello et al. 2010). However, there is also the in- selected for observations with the Jansky Very Large Array triguing possibility that such samples may contain unlensed (JVLA). Here we detail the JVLA observations and those ob- and thus intrinsically luminous galaxies – the rarest HyLIRGs − − tained with other facilities that followed as a consequence of are predicted to have a space density, ≈ 10 7 Mpc 3 and our initial findings, presenting these datasets in order of de- may mark the sites of today’s clusters (Negrello et al. 2005; creasing wavelength. Lapi et al. 2011). Of the facilities providing redshifts for bright Herschel 2.1. JVLA CO J =1−0 and 5-GHz continuum imaging sources, the 100-m Robert C. Byrd Green Bank Telescope Whilst the National Radio Astronomy Observatory’s (GBT), together with the ultrawide-bandwidth Zpectrome- 39 ter cross-correlation spectrometer and Ka-band receiver, have (NRAO’s) JVLA was in its DnC, C, B, BnA and A configu- been amongst the most effective (e.g. Swinbank et al. 2010; rations, between 2012 January and 2013 January, we acquired Frayer et al. 2011). Having added considerably to the sam- ≈30hr of Ka-band data, scheduled dynamically to ensure ex- ple of such sources with known redshifts, Harris et al. (2012) cellent atmospheric phase and pointing stability. We recorded ′ 12 16 − 38 two sets of eight contiguous baseband pairs, 1024 × 2-MHz identified a flat trend of LCO with C O J =1 0 line width. This contrasted with the steep power-law relation between dual-polarization channels in total. We tuned the first set of 13 16 12 18 − 12 16 L′ and CO line width foundbyBothwell et al. (2013) forun- basebands to cover C O and C O J =1 0. The C O CO − lensed SMGs, as expected if the most gas-richgalaxies tend to J =1 0 transition (νrest = 115.271203GHz, Morton & Noreau live in the most massive gravitational potentials. Harris et al. 1994) was placed in the second set of baseband pairs, offset- argued that intrinsically fainter sources require a higher lens- ing down by 64MHz to 33.740GHz to avoid baseband edges. ing magnification, µ, to rise above the observational detection Around 4hr of A-configuration C-band data were also ob- threshold for CO J =1−0, which is approximately constant tained, recording 1024 × 2-MHz dual-polarization channels in flux. The flat trend seen in the GBT sample is then best every 1sec across 4.2–6.5GHz, with a small gap. Typically, interpreted as lensing magnification acting on a population 0.3GHz was lost to severe radio-frequency interference near with intrinsically steep number counts. Those GBT sources 6.1GHz, yielding a band center of 5.1GHz (5.9cm). closest to the power-law fit seen for SMGs – typically those PKSJ0825+0309 and PKSJ0839+0319 were observed ev- with the widest lines – are likely to be suffering the least lens- ery few minutes to determine accurate complex gain solutions ing magnification. This is where we might expect to find any and bandpass corrections in the Ka and C bands, respectively. HyLIRGs that are lurking amongst the lensed starbursts. 3C286 was observed to set the absolute flux density scale, Our approach here, therefore, is to search for HyLIRGs and the pointing accuracy was checked locally every hour. amongst those bright Herschel lens candidates with the The data were reduced using AIPS. The basebands were broadest CO lines, starting with the best example in knitted together using the NOIFS and VBGLU tasks, yielding the sample observed with GBT by Harrisetal., HAT- uv datasets with two intermediate frequencies, each compris- LASJ084933.4+021443 (hereafter HATLASJ084933, with ing 512 × 2-MHz channels. 350-µm flux density, S = 293mJy), which displays a line For the spectral-line data, these were then imaged in groups 350 −1 consistent with CO J =1−0 at zLSR =2.410 ± 0.003, with a of four channels (71kms ), with natural weighting. A vari- − full width at half maximum (FWHM) of 1180 ± 320kms 1 ety of Gaussian tapers were used to weight the data, with dis- and SCO1−0 =0.83 ± 0.19mJy. The amplification predicted by tances to the 30 per cent point of the Gaussian ranging from 100 to 1,000kλ, to form a number of cubes with FWHM spa- Harris et al. for HATLASJ084933 is consistent with unity, ′′ ′′ ′′ ′′ ± tial resolutions ranging from 0.53 × 0.51 to 2.0 × 1.8 . µ =2 1. ′′ ′′ In the next section we describe an extensive set of ob- The line-integrated map shown in Fig. 1 has a 1.1 × 1.0 servations. We present, analyse, interpret and discuss our synthesized beam. reduced images, spectra and cubes in §3, finishing with The naturally-weighted 5.1-GHz pseudo-continuum data our conclusions in §4. We adopt a cosmology with H0 = were imaged, averaging all of the channels, loosely following −1 −1 Ω Ω ′′ the techniques outlined by Owen & Morrison (2008). The re- 71kms Mpc , m =0.27 and Λ =0.73, so 1 equates to ′′ ′′ ◦ 8.25kpc at z =2.41. sulting map has a spatial resolution of 0.53 ×0.42 (PA, 26 ) and an r.m.s. noise level of 2.9 µJybeam−1 (Fig. 1).

39 This work is based on observations carried out with the JVLA. The NRAO is a facility of the NSF operated under cooperative agreement by As- 38 Hereafter, CO refers to 12C16O unless stated otherwise. sociated Universities, Inc. HYLIRGS AMIDST DISTANT STARBURSTING CLUSTER 3

2.2. CO J =3−2 imaging from CARMA 2.5. 350-µm continuum imaging with CSO/SHARC-2 40 We used CARMA to observe the CO J =3−2 line (νrest = A total of 80min of data were collected at 350 µm during 345.795991GHz, redshifted to νobs=101.406GHz). Obser- periods of excellent weather between 2012 January 12–13 at vations were carried out with the 3-mm receivers and the the 10.4-m Caltech Submillimeter Observatory (CSO43), us- CARMA spectral-line correlator, with an effective band- ing the SHARC-2 camera (Dowell et al. 2003). We took ad- width of 3.7GHz per sideband and a spectral resolution of vantage of the open-loop, actuated surface of the telescope 15kms−1. Four tracks were obtained between 2012 Febru- (Leonget al. 2006) and we checked and updated the focus ary 13–17 in the C configuration (15–352-m baselines), with settings several times each night, resulting in near-Gaussian, 14 usable antennas and a total on-source observing time of diffraction-limited (8.5′′ FWHM) beam profiles. Data were 13.5hr. The nearby quasar, PKSJ0757+0956, was used for reduced using CRUSH (Kovács 2008) and the resulting maps complex gain calibration. The bandpass shape and abso- have a noise level, σ = 8mJybeam−1 (Fig. 1). lute flux calibration (the latter good to ∼15%) were derived from observations of Mars, PKSJ0854+2006 and 3C84. Us- 2.6. Continuum imaging from Herschel ing natural weighting, we obtained a FWHM beam size of The acquisition and reduction of Herschel (Pilbratt et al. − 2.1′′ ×1.8′′ (PA, 126◦), with a noise level, σ =0.3mJybeam 1 2010) parallel-mode data from SPIRE (Griffin et al. 2010) (Fig. 1). and PACS (Poglitsch et al. 2010) for the 9-hr Science Demon- stration Phase field of H-ATLAS (Eales et al. 2010) are de- − 2.3. CO J =4 3 imaging from IRAM PdBI scribed in detail by Ibar et al. (2010a), Pascale et al. (2011) During 2012 February we obtained data using six 15-m an- and Rigby et al. (2011). Here, we use the original SPIRE tennas in the most extended configuration of the Institut de imaging and have obtained much deeper PACS data from the Radioastronomie Millimétrique’s Plateau de Bure Interferom- OT1 programme, OT1_RIVISON_1. We recorded data simul- eter (IRAM41 PdBI), with baselines ranging up to 800m. taneously at 100 and 160 µm whilst tracking across the target The observing frequency was set to 135.203GHz, the at 20′′ s−1 in ten 3-arcmin strips, each offset orthogonally by ′′ redshifted frequency of the CO J = 4 − 3 line (νrest = 4 . A total of 180s on source was acquired for each of two 461.04077GHz). The weather conditions were excep- near-orthogonal scans. These data were tackled with a vari- tionally good with a precipitable water vapour of around ant of the pipeline developed by Ibar et al. (2010b), reaching 2mm and Tsys ∼ 100K on average. PKSJ0825+0309 σ ≈ 4 and 7mJy at 100 and 160 µm, respectively. The 160- and PKSJ0909+0121 were used to calibrate the complex µm PACS image is shown in Fig. 1. gains, while 3C84 and MWC349 were chosen as band- pass and absolute flux calibrators. The r.m.s. noise level 2.7. Infrared continuum imaging from Spitzer, VISTA and is 0.97mJybeam−1 in 4-MHz-wide spectral channels, for a the Hubble Space Telescope synthesized beam measuring 1.0′′ × 0.5′′ FWHM (PA, 24◦). 3.6- and 4.5-µm images were acquired using the Infrared Using only the line-free channels, the r.m.s. noise level is Array Camera (IRAC – Fazio et al. 2004) aboard Spitzer44 47 µJybeam−1. (Werneretal. 2004) on 2012 June 24 as part of program We also report continuum photometry at 1.36mm from ob- 80156. The imaging involved a 36-position dither pattern, servations obtained throughout 2012 as part of a search for with a total exposure time of just over 1ks, reaching r.m.s. water that will be reported elsewhere. depths of 0.3 and 0.4 µJy at 3.6and 4.5 µm, respectively. Cor- rected basic calibrated data, pre-processed by the Spitzer Sci- 2.4. 870-µm continuum imaging from the SMA ence Center, were spatially aligned and combined into mo- saics (shown in Fig. 1) with a re-sampled pixel size of 0.6′′ Approximately 1.4, 4.9 and 1.7hr of integration were ob- ′′ tained using the Submillimeter Array (SMA42) in its com- and angular resolution of 2–2.5 , using version 18.5.0 of pact, extended and very-extended array configurations during MOPEX (Makovoz & Marleau 2005). 2011 December 8, 2012 January 31 and 2012 March 31, re- Images were obtained in Z,Y,J,H and Ks as part of spectively (2011B-S044). The receivers were tuned such that VIKING, a public survey with the 4-m Visible and Infrared Survey Telescope for Astronomy (VISTA) at Paranal, Chile the upper sideband was centered at 345GHz (870 µm), with ′′ 8GHz of single-polarization bandwidth in total. 3C84 was (Emerson & Sutherland 2010), with an image quality of 0.9 used as the bandpass calibrator and Titan was used for abso- FWHM, as shown in Fig. 1. lute flux calibration. PKSJ0825+0309 and PKSJ0909+0121 A SNAPshot observation was obtained with the Hub- were again used to track and check the phase and gain. ble Space Telescope45 (HST) on 2012 February 08, as part ′′ A ROBUST = 2 weighting scheme resulted in a 0.83 × of Cycle-19 proposal 12488, using Wide-Field Camera 3 0.67′′ beam FWHM (PA, 86◦) and a noise level, σ = (WFC3) with its wide J filter, F110W. The total exposure 0.77mJybeam−1 (Fig. 1). time was 252s. Data were reduced using the IRAF MultiDriz- zle package. Individual frames were corrected for distortion, 40 CARMA construction was derived from the states of Maryland, Califor- nia and Illinois, the James S. McDonnell Foundation, the Gordon and Betty 43 This material is based upon work at the CSO, which is operated by Moore Foundation, the Kenneth T. and Eileen L. Norris Foundation, the Uni- the California Institute of Technology under cooperative agreement with the versity of Chicago, the Associates of the California Institute of Technology, National Science Foundation (AST-0838261). and the NSF. Ongoing CARMA development and operations are supported 44 This work is based in part on observations made with the Spitzer Space by the NSF under a cooperative agreement, and by the CARMA partner uni- Telescope, which is operated by the Jet Propulsion Laboratory, California versities. Institute of Technology under a contract with NASA. 41 IRAM is supported by INSU/CNRS (France), MPG (Germany) and 45 Based on observations made with the NASA/ESA Hubble Space Tele- IGN (Spain). scope, obtained at the Space Telescope Science Institute, which is operated 42 The SMA is a joint project between the Smithsonian Astrophysical Ob- by the Association of Universities for Research in Astronomy, Inc., under servatory and the Academia Sinica Institute of Astronomy and Astrophysics NASA contract NAS 5-26555. These observations are associated with pro- and is funded by the Smithsonian Institution and the Academia Sinica. gram 12488. 4 IVISONETAL. cleaned by cosmic rays and other artifacts and median com- bined. The resulting ∼ 2′ ×2′ image was re-sampled to a finer Table 1 pixel scale of 0.064′′ and is shown in Fig. 1. Continuum flux densities WT 2.8. Optical spectroscopy from Keck 0.88 µm 8.16 1.37 5.80 0.98 µJy; VISTA Z We observed several optical/IR sources in the vicinity of 1.02 µm 10.5± 1.8 8.65 ± 1.45 µJy; VISTA Y 1.1 µm 13.6 ± 2.3 11.2± 1.9 µJy; F110W HATLASJ084933 using the DEIMOS spectrograph on the ± ± 46 1.25 µm 16.8 2.8 13.0 2.2 µJy; VISTA J 10-m Keck-II telescope during 2012 March 02. DEIMOS 1.65 µm 23.5 ± 4.0 19.2 ± 3.2 µJy; VISTA H −1 ± ± was used with its 600ZD grating (ruled with 600 linesmm ), 2.15 µm 38.1 6.4 26.1 4.4 µJy; VISTA Ks the GG455 filter and a 1′′ slit at a central wavelength of 3.6 µm 67.1 ±18.6 47.7 ±17.4 µJy; IRAC 4.5 µm 76.6 ± 20.1 45.8 ± 15.2 µJy; IRAC 720nm. Two positions were targeted, each with a different 100 µm 23± 7a 19± 7a mJy; PACS PA, for 40 and 75min. 160 µm 91 ±12a 40 ±12a mJy; PACS 250 µm± 242 18b ± mJy; SPIRE ± b 3. RESULTS, ANALYSIS AND DISCUSSION 350 µm 293 22 mJy; SPIRE 500 µm 231 ± 19b mJy; SPIRE 3.1. Basic morphological description 870 µm 25 2± 19 2 mJy;SMA 1.36 mm 8.3 ±0.5 7.5 ±0.5 mJy; PdBI At the spatial resolution of SPIRE, HATLAS J08493 is 2.21 mm 1.1 ± 0.1 0.8 ± 0.1 mJy; PdBI essentially unresolved. However, with only modestly im- 3.0 mm 0.33 ± 0.14 3σ <±0.42 mJy; CARMA proved spatial resolution, via deep imaging with PACS at 8.8 mm 0.043 ± 0.010 3σ < 0.030 mJy;JVLA 5.9 cm 0.177 ± 0.015 0.065 0.015 mJy;JVLA 100–160µm and CSO/SHARC-2 at 350 µm, the system is ± ± resolved into two sources, referred to hereafter as W and T C M (Fig. 1). 0.88 µm 5.96 1.00 – µJy; VISTA Z 1.02 µm 5.46 ± 0.92 – µJy; VISTA Y Moving to higher spatial resolution, via interferometric 1.1 µm 7.66 ± 1.29 – µJy; F110W imaging with SMA, IRAM PdBI and CARMA at rest-frame 1.25 µm 5.72 ± 0.96 – µJy; VISTA J 1.65 µm 5.09 ± 0.86 – µJy; VISTA H 255–880µm, we see a consistent morphological picture for ± the emission from cool dust in the HATLASJ084933 sys- 2.15 µm 11.1 1.9 – µJy; VISTA Ks 3.6 µm 19.5 ± 5.9 – µJy; IRAC tem (Fig. 1). The dust continuum emission is dominated 4.5 µm 17.7 ± 5.3 – µJy; IRAC ′′ by W and T, which are separated by 10.5 , or ∼85kpc in 870 µm 4.6 ±1.9 6.9 1.6 mJy;SMA the plane of the sky, with W the brightest of the two from 1.36 mm 1.8 ± 0.2 1.5 ± 0.2 mJy; PdBI 2.21 mm 0.25 ± 0.06 0.31 ± 0.08 mJy; PdBI the optical all the way through to the radio regime. W is 5.9 cm 0.022 ± 0.006 0.020 ± 0.005 mJy;JVLA amongst the brightest unlensed SMGs ever observed inter- ± ± a ferometrically at submm wavelengths, with T not far behind Errors determined by placing many apertures across the map, with 3- and 5-% calibration uncertainties added in quadrature at 100 and 160 µm, respec- (Gear et al. 2000; Lutz et al. 2001; Younger et al. 2007, 2008, tively. 2009; Wang et al. 2007; Aravena et al. 2010; Smolciˇ c´ et al. b Errors include the contribution due to confusion and a 7-% calibration un- 2012a,b). Amongst the unlensed SMG population, its flux certainty has been added in quadrature (Valiante et al., in preparation). density has been equalled only by SMMJ123711+622212 (also known as GN20, S870 = 22.9 ± 2.8 – Pope et al. 2005; Iono et al. 2006). Continuum emission from dust is also seen east of W, from Aided by our interferometric astrometry we can see that at a galaxy we label M in Fig. 1. M is detected securely in 870- optical/IR wavelengths the westernmost clump, T, is barely µm, 1.36-mm, 2.21-mm and 5.9-cm continuum (Table 1). visible as a red extension to a lenticular galaxy (Fig. 1). W Still further to the east, emission can be seen from a very coincideswith the extremely red componentof a red/blue pair, red, morphologically distorted galaxy – labelled C in Fig. 1 separated by ≈2′′. Based on experience with SMGs such as – which bears a strong resemblence to the Antennae (e.g. SMM J02399−0136 (e.g. Ferkinhoff et al., in preparation), we Whitmore& Schweizer 1995; Klaas et al. 2010). C is de- expected to find evidence of an AGN towards the compact tected weakly in 870-µm, 1.36-mm, 2.21-mm and 5.9-cm blue source, but found instead that it is a star – our optical continuum (Table 1). spectroscopy (§2.8) reveals absorption due to Mg as well as As we shall see in more detail later (§3.5), C and M are also the Na D doublet. detected in our CO J =1−0 and J =4−3 spectral-line imaging, There is no obvious trend with wavelength for the relative i.e. they also lie at z =2.41, alongside W and T. IR/submm flux density contributions from W and T (listed Both M and C are coincident with IRAC 3.6–4.5-µm emis- in Table 1), SW/ST =1.22 ± 0.33, 2.22 ± 0.42, 1.27 ± 0.19, sion, suggesting that these are not tidal tails or gaseous 1.32 ± 0.17, 1.13 ± 0.10 and 1.23 ± 0.25 at 100 µm, 160 µm, streams, but rather significant concentrations of stars, gas and 350 µm, 870 µm, 1.36mm and 2mm, respectively. The dust – luminous SMGs in their own right. error-weighted mean for the IR/submm wavelength regime is Our deep rest-frame 133–289-nm Keck spectroscopy 1.22±0.07. Moving longward to rest-frame1.7cm, where we (§2.8), with slits covering W, T and C, reveals no significant expect the emission to be dominated by synchrotron radiation line emission from T and C; the only line visible from W, al- related to supernova (SN) remnants, HATLASJ084933 is also beit very faint, is C II at 232.6nm. − dominated by component W (Fig. 1), with SW/ST =2.7±0.7. Our CO J =3 2 image from CARMA is shown in Fig. 1. It is interesting to note that the low-level southern extension to 46 Data presented herein were obtained at the W. M. Keck Observatory, the J =3−2 emission from T is mirrored at 3.6–4.5-µm, which which is operated as a scientific partnership among the California Institute of suggests that both are real. Technology, the University of California, and NASA. The Observatory was made possible by the generous financial support of the W. M. Keck Founda- tion. 3.2. Lensing model HYLIRGS AMIDST DISTANT STARBURSTING CLUSTER 5

Figure 1. Morphology of HATLAS J084933 at moderate-to-high spatial resolution. Top left: 160-, 350- and 870-µm continuum emission from dust, as detected by Herschel/PACS, CSO/SHARC-2 and the SMA, superimposed on a three-color representation of the VISTA z,J and Ks data. We see two bright components: W is centered on a red galaxy, near a blue star; T lies at the end of a lenticular galaxy at z = 0.3478, coincident with its red southern tip. Dust emission can also be seen coincident with two further galaxies, labelled C and M, to the east of W. Both C and M can also seen in continuum at 5.9 cm (see lower right), at 1.36 mm, 2.21 mm, and in CO J =1−0 and J =4−3. Top right: CO J =3−2 emission superimposed on a three-color representation of the Spitzer IRAC 3.6- and 4.5-µm imaging and a heavily smoothed VISTA J + H + Ks image as the blue channel. CO emission is coincident with the two submm-bright clumps, W and T. An apparent extension south of T is mirrored in the IRAC image, suggesting that both are real. Lower left: JVLA CO J =1−0 imaging superimposed on our HST WFC3 F110W snapshot, wherein we see faint diffraction spikes from the stars. Distorted, low-surface-brightness emission is seen coincident with galaxy C, characteristic of an interaction. Lower right: 1.25′′ 1.25′′ stamps centered on M, C, W and T, showing our JVLA 5.8-cm continuum imaging as blue contours superimposed on the HST F110W imaging.≈ At the× high spatial resolution available here via interferometric imaging, we see a consistent morphological picture for the emission from cool gas and dust in the HATLAS J084933 system. The emission is dominated by W and T, which are separated from one another by 85 kpc in the plane of the sky, with significant emission also visible from two fainter galaxies, C and M. Individually, W and T are well resolved spatially. ∼ In all panels, contours are plotted at −3,3 σ, with √2-spaced increments thereafter, where σ is the local noise level. Beam FWHM are shown as solid ellipses. × ◦ ◦ N is up; E to the left; offsets from α2000 = 132.3889 ,δ2000 = 2.2457 are marked in arcseconds. Known stars are labelled. 6 IVISONETAL.

Figure 2. Radio-through-optical SEDs of the dominant components of HAT- LAS J084933, W (red points, solid line) and T (blue points, dotted line). Our measurements cover the SED peak, constraining LIR and Tdust well. The SEDs can be fitted adequately using a combination of thermal dust and syn- chrotron emission components (see §3.3) plus the stellar population mod- elling described in §3.4. The key resulting characteristics, LIR, Tdust, Mdust, Mstars and qIR are listed in Table 2.

The lenticular galaxy that lies along the line of sight to component T lies at z =0.3478, as revealed by our optical spectroscopic observations (§2.8), where the slit was placed along the major axis of the lenticular. Correcting for the in- strumental resolution, the lines due to Mg and Ca H and K suggest the foreground galaxy has a velocity dispersion of −1 11 190 ± 80kms . Its stellar mass is (1.1 ± 0.4) × 10 M⊙, adopting the approach outlined later in §3.4. We have mod- elled this galaxy with an elliptical (e =0.3) mass profile, trun- cated smoothly at a radiusof 25kpc (theexact choicehas little effect since the lensed image is well within this radius). The line of sight to the z =2.41 starburst lies close to the semi- major axis of this foreground structure and suffers a magnifi- cation of µ =1.5 ± 0.2, where the uncertainty was calculated by varying the model parameters – orientation, ellipticity and velocity dispersion – by 10%. Figure 3. CO spectra of HATLAS J084933. Upper two panels: CO J =1−0 Galaxy W suffers no significant lensing magnification, as spectra of galaxies W and T superimposed on solid representations of the CO J =4−3 spectra (the former scaled by 16 to put them on the same Tb scale). far as it is possible to discern from our extensive, wide-field, The two galaxies both have extraordinarily× broad CO lines but the profiles panchromatic imaging data (§2.7), which are not easily rec- are noticably different, so they are not lensed images of the same galaxy. The onciled with a massive foreground galaxy group or cluster. profiles of the CO transitions are indistinguishable for each individual com- ponent; their Tb ratios are also similar (Table 2). Middle two panels: CO J =4−3 spectra of galaxies M and C, with line centers significantly redwards 3.3. Spectral energy distributions of W and T of W and T. Lower panel: JVLA CO J = 1−0 spectra of W and T along- side the significantly different profile from GBT/Zpectrometer (Harris et al. Fig. 2 shows the SED of HATLASJ084933, concentrating 2012). Emission at velocities redward of W and T – from M and C and sev- on the brightest components, W and T. Our CSO/SHARC- eral other faint clumps seen in the JVLA data, spanning a 100-kpc region – 2 350-µm imaging provides evidence that – at the peak of can account for the differences between the GBT and JVLA∼ line profiles. their SEDs – the relative SPIRE contributions of components W and T follow the average ratio seen in the resolved imag- as is LIR), and to Mdust (for a characteristic photon cross- 2 −1 ing, 1.22±0.07, as described in §3.1. Perhaps unsurprisingly, section to mass ratio, κ850 = 0.077m kg – Dunneet al. then, we then find they contribute in roughly the same ratio, 2000; Dunne & Eales 2001 – where we fixed the frequency 1.55 ± 0.44, to an immense overall bolometric luminosity of dependence of the dust emissivity, β, to be 2.0, which is 13 HATLAS J084933, LIR =(4.8±1.0)×10 L⊙ (Table 2), hav- physically plausible and consistent with the best-fit value, ing corrected the contribution from T for gravitational ampli- 2.08 ± 0.15). We fixed the synchrotron power-law index, α = α fication (§3.2). −0.75, where Sν ∝ ν , and we also fixed γ =7.2 – the power- The IR-through-radio photometry can be described ade- law index of a dust temperature distribution appropriate for −γ quately with a model comprising synchrotron and thermal local starbursts, dMdust/dTdust ∝ Tdust, which is designed to dust emission, following Kovács et al. (2010). Table 2 sum- offer a physically motivated treatment of the Wien side of the marises a simultaneous fit to the dominant cold dust tem- thermal emission spectrum. Both W and T have Tdust ∼ 36– perature, Tdust, to LIR and qIR (as defined by Helou et al. 40 K, commensurate with similarly luminous, dusty starbursts 1985, but where SIR is measured across λrest = 8–1000 µm, at z ∼ 2–3, when calculated using a power-law temperature HYLIRGS AMIDST DISTANT STARBURSTING CLUSTER 7

Table 2 Properties of HATLAS J084933.4+021443

Property W T M C R.A. (J2000) 08:49:33.59 08:49:32.96 08:49:33.80 08:49:33.91 Dec. (J2000) +02:14:44.6 +02:14:39.7 +02:14:45.6 +02:14:45.0 Mean FWHM at 5.9 cm, 870 µm and in CO J =4−3 2.9 0.4 kpc 3.8 0.6 kpca b b FWHM in CO J =1−0 /kpc 7.0 ± 2.1 kpc 6.2 ± 1.4 kpca b b ± ± a log LIR /L⊙ 13.52 0.04 13.16 0.05 12.9 0.2 12.8 0.2 −1 ± ± a ± ± SFR /M⊙ yr (Kennicutt 1998, Chabrier 2003 IMF) 3400 1500 800 640 b b Tdust /K 39.8 1.0 36.1 1.1 ± ± b b qIR 2.30 0.08 2.53 0.13 ± ± a b b log Mdust /M⊙ 9.32 0.05 9.10 0.05 −1 ± ± CO J =1−0 ICO /Jy km s 0.49 0.06 0.56 0.07 0.057 0.013 0.079 0.014 − ′ 9 −1 2 ± ± ± ± CO J =1 0 LCO /10 Kkms pc 138 17 157 20 16.0 3.6 22.2 3.9 6 ± ± ± ± CO J =1−0 LCO /10 L⊙ 6.72 0.82 7.65 0.96 0.779 0.178 0.108 0.019 CO J =1−0 FWHM /kms−1 825 ± 115 610± 55 320± 70 250 ± 100 ± ± ± ± CO J =1−0 zLSR 2.4066 0.0006 2.4090 0.0003 2.4176 0.0004 2.4138 0.0003 13 16 12 18 −1 ± c ± c ± c ± c C O and C O J =1−0 ICO /Jy km s 3σ < 0.24 3σ < 0.24 3σ < 0.24 3σ < 0.24 −1 CO J =3−2 ICO /Jy km s 4.08 0.92 5.73 1.13 0.66 0.27 1.16 0.37 CO J =3−2 FWHM /kms−1 830 ± 155 555± 75 260± 95 620 ± 160 ± ± ± ± d CO J =3−2 zLSR 2.4077 0.0009 2.4096 0.0004 2.4173 0.0006 2.414 −1 ± ± ± CO J =4−3 ICO /Jy km s 3.76 0.27 5.10 0.37 0.66 0.10 0.95 0.15 CO J =4−3 FWHM /kms−1 985 ± 115 545± 30 320± 50 450± 90 ± ± ± ± CO J =4−3 zLSR 2.4068 0.0002 2.4090 0.0002 2.4178 0.0003 2.4149 0.0004 ± ± ± ± Tb r − − 0.93 0.24 1.14 0.27 1.29 0.60 1.63 0.59 3 2/1 0 ± ± ± ± Tb r4−3/1−0 0.48 0.07 0.57 0.08 0.72 0.20 0.75 0.18 −1 ± ± ± ± vmax /km s 350 10 220 15 20 10 80 20 R /kpc 7.8 ±1.0 6.8 ±0.8a± b± b ± ± a log Mdyn /M⊙ 11.51 0.05 11.13 0.05 11.11 0.15 11.11 0.20 −1 2 −1 ± ± a ± ± log MH2+He /M⊙, for αCO = 0.8 M⊙ (K km s pc ) 11.04 0.05 10.92 0.06 10.11 0.09 10.25 0.07 −1 ± ± a ± ± SFE /L⊙ M⊙ 370 210 760 440 a log Mstars /M⊙ 11.38 0.12 11.01 0.12 10 10.36 0.18 ± ± ≈ ± a Corrected for µ = 1.5 ± 0.2. b Insufficient S/N to make a useful measurement. c 12 16 Adopting Stotal/Speak as seen for C O. d Fixed. distribution (Magnelli et al. 2012). SHARC-2 recovers less flux than SPIRE at 350 µm, a possi- We are sensitive to AGN activity via our optical imaging ble manifestation of the proto-cluster. and spectroscopy (§2.8) and our radio/mid-IR imaging (§2.1, §2.7). In relation to the far-IR/radio flux ratio, W and T have qIR =2.30±0.08and 2.53±0.13, respectively, spanning (and consistent with) the tight correlation for star-forming galaxies 3.4. Stellar masses (qIR ≈ 2.4, e.g. Yun et al. 2001). Summarizing, we have no The rest-frame ultraviolet-through-near-IR SEDs of W and 47 evidence that an AGN contributes significantly to the LIR T are illustrated in Fig. 2, alongside the far-IR–submm–radio of HATLASJ084933, though of course we have no evidence data. Flux densities are listed in in Table 1. Our high-quality of absence since even a powerful AGN can be hidden very HST, VIKING and IRAC images at 0.88–4.5µm are compro- effectively by the quantities of dust-rich gas in this system. mised to a certain extent by the unfortunatesuper-position ofa If contributions to LIR from AGN are small, then W and T star and a lenticular galaxy near W and T, respectively, but the −1 are generating stars at a rate of ≈3400 and ≈1500M⊙ yr , considerable uncertainties associated with estimating stellar respectively, for a Chabrier (2003) IMF (Kennicutt 1998, mass dominate over the photometric uncertainties (see, e.g., though see Dwek et al. 2011). The estimates of LIR for C and Hainline et al. 2011; Michalowski et al. 2012). M given in Table 2 are scaled from their submm photometry For components W, T and C we find best-fitting monochro- relative to W and T, assuming the same SED shape. This sug- matic rest-frame H-band luminosities, logLH, of 11.82, gests that both are close to being HyLIRGs in their own right, 11.46 and 11.12L⊙, respectively. Stellar masses were de- 12 ≈ 7 × 10 L⊙. termined using MAGPHYS with the most recent version of Before closing our SED discussion, we add a word of cau- the stellar population models of Bruzual & Charlot (2003), tion regarding our estimated IR luminosities. If W and T lie a Chabrier (2003) IMF and the SED prior distributions in the core of a proto-cluster, as we shall see in what follows, described by da Cunha, Charlot, & Elbaz (2008), following then lower-luminosity star-forming cluster galaxies may con- Rowlandsetal. (2012). The models encompass a wide tribute to the flux densities measured by the relatively large range of exponentially declining star-formation histories SPIRE beam at 250–500µm (see Negrello et al. 2005). Al- (SFHs), with bursts superimposed (see da Cunha et al. 2008, though the calibration cannot compare with that of Herschel, for details). The implied stellar masses are logMstars = +0.11 +0.12 +0.18 11.38−0.13,11.01−0.11 and 10.36−0.17 M⊙ for W, T and C, 47 Observations of higher-J CO lines would be useful in terms of AGN where the errors are derived from the marginalised probabil- diagnostics. ity density function, which incorporates the uncertainties in 8 IVISONETAL. the SFH48 and photometry. Systematic shifts in stellar mass 3.6. Resolved imaging of the gas and dust of ≈ 0.2dex result from adopting a Salpeter IMF or ignoring Our interferometric images from JVLA, IRAM PdBI and thermally pulsating stars on the asymptotic giant branch, and the SMA have a spatial resolution of ∼ 0.5′′ FWHM, or the magnitude of photometric contamination by buried AGN ∼ 4kpc. They allow us to compare the sizes of the regions could well be similar, so we estimate that our stellar mass es- responsible for emission from the total molecular gas reser- timates are accurate to ≈ 0.3dex. voir (via CO J =1−0 emission), from the star-forming gas (via Emission from M can be seen at 3.6–4.5 µm but the photo- CO J =4−3), from cool dust produced by SNe in regions of metric uncertainties due to blending with C and W are large recent star formation (via submm continuum emission) and so we state only that its stellar mass must be considerable, 10 from relativistic electrons spiralling in or near recent SN rem- ≈ 10 M⊙. nants (via radio continuum emission). There is little sign of the power-law SEDs that can betray We find that emission from W and T is resolved in all these the presence of AGN, supporting the lack of AGN spectral wavebands, with total flux densities significantly higher than features reported in §3.1, though of course this does not rule the peak flux densities. In Table 2 we quote their deconvolved out the presence of a deeply buried, highly obscured AGN. FWHM50 sizes at 5.9cm, 870 µm and in CO J =4−3 and J = 1−0, as determined with 2-D Gaussian fits. 3.5. Spectral-line characteristics In the absence of significant AGN-powered activity, ra- Fig. 3 shows the JVLA CO J =1−0 spectra of W and T dio continuum observations trace regions where massive stars alongside our IRAM CO J =4−3 spectra. These spectral-line have recently been formed, with no obscuration. Of the trac- data bring us velocity information for the first time, enabling ers available to us, we might anticipate the radio emission us to see that although W and T do lie at approximately the to be the most compact, as indeed is measured. At 870 µm same redshift and share similar line profiles and Tb ratios, they we see no evidence of a more extended emission component, have significantly different line widths (Table 2). Each one is heated by older stars – the deconvolved sizes of W and T are thus a distinct, extraordinarily luminous starburst. While the consistent, as measured in 5.9-cm and 870-µm continuum, as CO line profile of T is typical of SMGs (e.g. Bothwell et al. one might expect given the far-IR/radio correlation (Table 2). − − 2013), that of W is broader than most, with ≈1000kms 1 The CO J =4 3 FWHM are slightly larger, as seen for local FWZI. Both profiles are flat-topped or perhaps double-peaked; ULIRGs (e.g. Wilson et al. 2008), but are consistent with the they can not be described well by a single Gaussian. radio and submm values to within the uncertainties. The CO J =1−0 fluxes of W and T measured in our JVLA The error-weighted average of the 5.9-cm and 870-µm con- −1 tinuum and CO J =4−3 FWHM measurements is (3.4 ± data sum to ICO1−0 = 1.05 ± 0.09Jykms . At first sight ◦ this is consistent with the Harris et al. GBT measurement, 0.4)kpc × (2.4 ± 0.4)kpc (PA, 9 ± 15 ) for W; for T, before 1.04 ± 0.37Jykms−1. However, the CO J =1−0 line profiles correcting for the amplification, we measure (5.6±0.7)kpc× ◦ − observed by JVLA and GBT appear noticeably different in (3.5 ± 0.7)kpc (PA, 48 ± 16 ). In CO J =1 0, however, W Fig. 3: emission in the GBT spectrum continues significantly and T have considerably larger deconvolved sizes than those redward (by ≈500kms−1) of the emission that can attributed measuredin the submm or radio continuum, as seen for SMGs to W and T. generally (e.g. Ivison et al. 2011; Riechers et al. 2011). Lookingat ourIRAM PdBIimaging–see Fig. 4– aftersub- Summarising, we have powerful, direct evidence of star- tracting the significant continuum emission (several submm bursts covering ≈ 3–4 kpc FWHM, with still larger reservoirs continuum flux densities are listed in Table 1), we see that of gas available to fuel future star formation, on scales of ≈ 6– both C and M are detected securely in CO J =4−3 (at 9.0 and 7kpc FWHM. 7.2 σ, respectively; we also see 3.6–4.5-µm, submm and ra- The large sizes measured in CO have been mentioned dio continuum emission – recall Fig. 1). These galaxies are alongside the relatively high L[CII]/LIR ratios measured to- therefore unquestionably real and they share approximately wards a number of SMGs (e.g. Stacey et al. 2010). Both have the same redshift as components W and T. Their positions are been taken as evidence that SMGs form stars across larger listed in Table 2, along with those of the two brighter compo- spatial scales than the compact, nuclear events seen in local − ULIRGs. nents. Their CO emission is centered ≈500kms 1 redward of W and T. 3.7. Dynamical characteristics As well as C and M, a number of other faint emission fea- tures can seen in the integrated CO J =1−0 image (Fig. 1). The most remarkable characteristics of HATLASJ084933 Summing C, M and these faint clumps leads to tolerable are revealed by the superlative spatial resolution (§3.6), ve- locity resolution and sensitivity of our JVLA CO J =1−0 and agreement between the GBT and JVLA line profiles. Al- − though we must add 0.87 ± 0.19Jykms−1 to IJVLA, we must IRAM PdBI CO J =4 3 data, as shown in Fig. 4. CO To measure the 2-D velocity structure, we fit the CO emis- GBT ≈ ∼ also increase ICO by 15% to account for the 8-per-cent sion lines, pixel by pixel, in each of the cubes. Initially, we lower flux density assumed for 3C286 by Harris et al. (2012), attempt to identify an emission line in 0.2′′ ×0.2′′ regions, in- and for losses due to attenuation of the faint clumps by the 49 creasing this slowly until we reach S/N > 5, at which pointwe GBT primary beam. The agreement between the total GBT fit the CO profile, allowing the intensity, centroid and width ≈ and JVLA line intensities remains consistent to within 1σ to vary. Uncertainties are calculated by perturbing each pa- ≈ ′′ and the profile measured by the FWHM 22 GBT primary rameter independently, allowing the remaining parameters to beam is thereby reconciled with the JVLA data. find their optimum values, until ∆χ2 = 1 is reached. 48 Stellar mass estimates are generally robust to changes in SFH (e.g. 50 We use half-light radii to calculate the area appropriate for gas and star- Ilbert et al. 2010, 2013; Pforr et al. 2012). formation surface densities, alongside a 0.5 correction to the SFR or gas 49 Note that the low-level GBT spectrometer baseline (at negative veloci- mass. We use FWHM simply as a convenient× way to compare sizes in several GBT ties relative to HATLAS J084933) may have caused an underestimate of ICO . wavebands. HYLIRGS AMIDST DISTANT STARBURSTING CLUSTER 9

Figure 4. CO J = 1−0 and J = 4−3 imaging of HATLAS J084933. Left: Velocity-integrated continuum-subtracted CO emission from our JVLA and PdBI observations. In both panels, we mark the positions of the brightest four galaxies, C, M, W and T. Right: Dynamical maps of the galaxies in CO J = 1−0 and J = 4−3. To construct these velocity fields, we fit the CO emission in each cube at each pixel, recording the intensity, velocity and line width where S/N > 5. North is up; East is left. Inset: One-dimensional rotation curves for W and T, extracted along the major kinematic axis identified in the modeling. Both W and T have rotation curves which resemble disks; indeed, W and T are best described as counter-rotating disks, which Mihos & Hernquist (1996) predicted would lead to the most intense starbursts.

Figure 5. CO J =4−3 imaging of galaxies W and T, showing left to right their intensities, resolved velocity profiles, line-of-sight velocity dispersions (σ), our best-fit kinematic model (see text) and the residuals. The color scale for the velocity fields is shown in each panel. We note that there is evidence for lensing shear in the velocity field and line-of-sight velocity dispersion for T. 10 IVISONETAL.

We model the two-dimensional velocity fields of W and unity in stable situations (in models with and without feed- T by constructing and fitting two-dimensional disk mod- back – Hopkinset al. 2012), we find values of ≈ 0.35 for −1 els using an arctan function, v(r)=2π vasym arctan(r/rt ), both W and T, with a plausible range that falls just short of to describe the kinematics (where vasym is the asymptotic unity. This suggests the disks are unstable and that the gas inclination-corrected rotational velocity and rt is the effective therein is prone to condense and form stars on a timescale radius at which the rotation curve turns over – see Courteau shorter than the rotational period. 1997). We fit six free parameters to the disk velocity field: vasym, rt, x,y centroid, PA and inclination to the sky. The x,y 3.9. Gas properties and αCO centroid and r must be within the field of view and the max- t The brightness temperature (T ) ratios measured for both imum circular velocity must be 2× the maximum velocity b < W and T are higher than expected (Table 2), where we usu- seen in the data. Quoted uncertainties reflect the range of ac- ally see r − − ∼ 0.52 ± 0.09 and r − − ∼ 0.41 ± 0.07 for ceptable models from all of the attempted model fits. 3 2/1 0 4 3/1 0 SMGs (Harris et al. 2010; Ivison et al. 2011; Bothwell et al. Our velocity cubes reveal that the wide CO line that led 2013), which suggests that the thermalised CO J =3−2 fill- to this system being selected as a candidate HyLIRG is due ing factors may be close to unity and that we may have found partly to the multitude of constituent galaxies and partly similar sizes in CO J=1−0 and J=3−2 if our spatial resolution to large regular, systematic rotational velocities in the two in the latter had been 2× higher. brightest constituent galaxies. Both W and T are resolved For a CO to H2+He conversion factor, αCO = 0.8M⊙ into counter-rotating gas disks – the scenario that many stud- −1 2 −1 ies have predicted should lead to the most intense starbursts (Kkms pc ) , as commonly adopted for IR-luminous star- (e.g. Mihos & Hernquist 1994, 1996; Taniguchi & Shioya bursts where the gas is not in virialized individual clouds (Bolatto et al. 2013), the gas masses we determine for W, T, 1998; Borne et al. 2000; Bekki 2001; Di Matteo et al. 2007; 9 Salomé et al. 2012). C and M are 110,83.7,12.8 and 17.8×10 M⊙ (corrected for As shown in Fig. 5, we find acceptablefits forW and Twith lensing). ◦ ◦ If we were instead to adopt the all-sample average disk inclinations of 56 ± 10 and 49 ± 10 , at PAs of 26 and + ◦ −1 60 145 , with vasym = 415 and 370kms , respectively. Compo- gas/dust mass ratio for nearby galaxies (Mgas/Mdust = 91−36 – nent T is sheared slightly by the foreground lens, as evident Sandstrom et al. 2013), we would determinegas masses forW 9 from the velocity and σ maps. Both galaxies show small-scale andTof19.0and11.5×10 M⊙, yielding gas/dust-motivated −1 −1 2 −1 (30–60-kms , r.m.s.) deviations from the best-fit model. The αCO values of 1.4 and 1.1M⊙ (Kkms pc ) , respectively. 51 dynamical masses of W and T are listed in Table 2: several Due to the uncertainty in κ850, etc., we regard these as ac- 11 2 2 curate to a factor 2× at best. For an environment with sub- ×10 M⊙, where Mdyn = Rvmax/sin (i). Their velocity dis- persions, σ, corrected for the instrumental velocity resolution, Solar metallicity, for example, αCO would rise in line with are ≈ 65kms−1. Mgas/Mdust. −1 A dynamically-motivated estimate of αCO can be made by M has no strong velocity gradient. Its vmax is 20kms − − determining the difference between the dynamical and stellar (40kms 1, peak to peak), σ = 140 ± 15kms 1 for both CO masses in each system, since this represents the total plausi- J =1−0 and J =4−3, and M ≈ 1 3 × 1011 M . For C, in CO dyn . ⊙ ble mass of gas54. Via this method, we find α = 0.6 and − −1 −1 CO J =1 0 we find vmax = 80kms (160kms , peak to peak), 0.3M (Kkms−1 pc2)−1 for W and T, respectively, where the −1 11 − ⊙ σ = 130±20kms and Mdyn ≈ 1.3×10 M⊙ (in CO J =4 3 uncertainty, ±0.3dex, is due almost entirely to the notoriously we find slightly higher values for σ and Mdyn). difficult estimate of stellar mass (§3.4). Finally, we can estimate αCO using its relationship with 3.8. Star formation in disks off the main sequence gas-phase metallicity and CO surface brightness, as devel- For main-sequence galaxies, star formation is not confined oped by Narayanan et al. (2012b), exploiting the resolved to a compact starbursting nucleus (e.g. Elbaz et al. 2011); at imaging in CO J =1−0 acquired here. For Solar metallicity z ∼ 2.4, the specific SFR is approximately 2.8Gyr−1. W (Chapman et al. 2005) and the sizes that have been measured −1 accurately in our JVLA imaging (§3.6; Table 2) we find αCO and T have sSFRs of 14.2 and 14.7Gyr , so their starbursti- − − 52 1 2 1 ness, sSFR/sSFRMS ∼> 5, placing them well above the main = 0.8 and 0.7M⊙ (Kkms pc ) for W and T, respectively. −1 2 −1 sequence. As we saw in §3.6 and §3.7 – their star forma- Our mean is αCO =0.8 ± 0.4M⊙ (Kkms pc ) , in line tion is not centrally concentrated as would be the case in lo- with the value usually adopted for starburst galaxies, with a cal starbursts. Instead, W and T are undergoing widespread range similar to that reported by Downes & Solomon (1998), −1 2 −1 star formation in disks that are supported partly by rota- 0.3–1.3M⊙ (Kkms pc ) . To reconcile our dynamical −1 2 −1 tion, as envisaged for main-sequence galaxies, and share measurement of αCO with 3.2M⊙ (Kkms pc ) , which many of the properties of massive star-forming disks (e.g. Genzel et al. (2010) argue is appropriate for all disks, regard- Förster Schreiber et al. 2009). Some other factor – presum- less of redshift, would require a buried AGN to dominate the ably their interaction – has taken them into the realm of sSFR near-IR light, causing us to over-estimate Mstars. reserved for starbursts. Adopting our best estimate for αCO, we find that W and For the Q parameter53 (Toomre 1964; T are both approaching the accepted criterion for maximal Goldreich&Lynden-Bell 1965), which describes the starbursts, with star-formation efficiencies, LIR/Mgas, of 370 stability of disks and which has been found to tend towards −1 and 210L⊙ M⊙ . C and M are higher still, though with con- siderable uncertainties. The Eddington limit is believed to 51 − R was taken to be twice the deconvolved half-light radius in CO J =1 0. be around 500L M−1; above this level, the radiation pres- 52 We cannot determine the IR8 starburst indicator (Elbaz et al. 2011) for ⊙ ⊙ HATLAS J084933 since we lack sufficiently deep rest-frame 8-µm imaging. 53 Usually ascribed to a seminal paper on the stability of stellar disks, 54 Bearing in mind that any AGN contribution to the rest-frame near-IR Toomre (1964), which contains no mention of Q. luminosity will have led to an over-estimate of stellar mass. HYLIRGS AMIDST DISTANT STARBURSTING CLUSTER 11 sure should quickly expel the gas via radiation pressure from young, massive stars on dust (Scoville 2004). Feedback due to the aforementioned radiation pressure will lengthen the total duration of star-formation episodes, or lead to a series of short bursts. Ignoring this, the gas- consumption timescale implied for the galaxies comprising HATLASJ084933 is τgas ∼ Mgas/SFR ≈ 15–45Myr, consis- tent with an interaction-driven period of enhanced star forma- tion (Scoville, Sanders, & Clemens 1986), as already hinted by the morphologically tortured state of C. The stellar masses determined in §3.4 for W, T and C suggest gas fractions, fgas = Mgas/(Mgas + Mstars = 36– 45%. These gas fractions are the top end of the range that Narayanan, Bothwell, & Davé (2012a) showed to be difficult to reconcile with cosmological simulations.

3.10. Relevance for star-formation laws Data of the quality gathered here have the potential to test whether the rules and prescriptions developed to describe star formation in the local Universe (e.g. Schmidt 1959; Kennicutt 1989) can be applied to to the interstellar medium of distant, gas-rich starbursts like HATLASJ084933. Are the data con- sistent with a simple, volumetric star-formation law in which the SFR is ∼1% of the molecular gas mass per local free- fall time, as argued by Krumholz, Dekel, & McKee (2012)? Are their SFRs more closely related to the orbital periods of their entire galactic disks (e.g. Genzel et al. 2010; Daddi et al. 2010), or to the quantity of gas above some density threshold (e.g. Lada et al. 2010; Heiderman et al. 2010)? We begin this discussion in the upper panel of Fig. 6: a ′ luminosity-luminosity plot of the observables, LIR and LCO, for low- and high-redshift star-forming galaxies (SFGs), e.g. z = 0 disks, LIRGs and ULIRGs and z = 1–2.5 BzK galaxies, taken from Genzel et al. (2010), and for SMGs. For this plot we need make no assumption about the appropriate value of αCO. The quality and depth of our continuum and CO data mean that the errorsbars for W and T can barely be discerned. ′ − Using LCO measurements made only in CO J =1 0, thereby avoiding the need to adopt a Tb ratio, alongside self-consistent measurements of LIR, we find no compelling evidence that different relations apply to SMGs and other SFGs, in line with thefindingsof Ivison et al. (2011). Here, we use the sampleof SFGs from Genzel et al. (2010) and find essentially the same result. It is noteworthy that while componentT contributes most to ′ ICO or LCO in all three observed CO transitions, component W is significantly brighter in the continuum at IR-through-radio wavelengths. We find a similar variation between W and T in relation to the far-IR/radio flux ratio, with qIR = 2.30 ± 0.08 and 2.53 ± 0.13, respectively, spanning (and consistent with) the tight correlation for star-forming galaxies (qIR ≈ 2.4, e.g. Yun et al. 2001). One might think their different line-to- continuum ratios would lead to a relatively large difference Figure 6. Plots relating the rate or surface density of star formation with the ′ in their positions on the plot of LIR versus LCO (Fig. 6, upper available gas. Components W and T are shown as red squares, with error bars, panel). In fact, W and T both lie comfortably amongst the in each plot. Circled red symbols denote high-redshift measurements made directly in CO J =1−0. Low- and high-redshift star-forming galaxies (SFGs), scatter of SMGs, suggesting that this scatter may be intrinsic e.g. z = 0 disks, LIRGs and ULIRGs and z = 1–2.5 BzK galaxies, are taken rather than due to measurement error. ′ from Genzel et al. (2010) for all three plots. Top: LIR versus LCO, where we Moving to the standard surface-density Schmidt-Kennicutt find no compelling evidence that SMGs follow a different relation than the ∝ N other star-forming galaxies. Middle: Surface-density S-K relation, where (S-K) relation, PSFR Pgas, shown as the middle panel of W and T lie amongst hyperluminous SMGs such as SMM J02399−0136 Fig. 6, it is now necessary to adopt value(s) for αCO appro- and GN 20, offset from the relation for other star-forming galaxies. Lower: priate for SFGs and SMGs. Genzel et al. (2010) argue that Elmegreen-Silk relation between star-formation surface density and the ra- −1 2 −1 tio of gas surface density and dynamical timescale. W and T exhibit extreme 3.2M⊙ (Kkms pc ) is appropriate for disks, regardless of star-formation efficiencies, even relative to SMGs. The fraction of their avail- redshift; our dynamical αCO constraints preclude the adoption able gas converted into stars per dynamical timescale, ǫτ , is > 0.1. 12 IVISONETAL.

−1 of this value for W and T. We use 3.2 and 0.8M⊙ (Kkms pc2)−1 for SFGs and SMGs, respectively, acknowledging that µ = 3

55 ] this is overly simplistic (e.g. Narayanan et al. 2012b). We 2 µ = 20

find W and T nestled amongst the other SMGs again (Fig. 6), pc

56 1 − T but this time the SMGs are offset significantly from the re- W lation for other SFGs. The lower panel of Fig. 6 shows the Elmegreen-Silk (E- K km s

S) relation between star-formation surface density and the ra- 10 tio of gas surface density and dynamical timescale57. Taken at face value, W and T exhibit extreme star-formation ef- ficiencies. Several other SMGs are equally extreme, con-

verting > 10% of their available gas into stars per dynami- [10 L'(CO)

cal timescale (ǫτ > 0.1). Krumholz et al. (2012) have argued 1 2 5 10 20 50 in favour of a simple volumetric star-formation law, linking Galactic clouds, SMGs, and all scales between the two, with 200 500 1000 − projection effects leading to the observed scatter. Here we FWHM linewidth [km s 1] have discovered two disks exhibiting star formation of such ′ − ferocity that they are capable of converting > 10% of their Figure 7. LCO versus CO J = 1 0 FWHM line widths for SMGs with red- molecular gas into stars with each rotation. We require a num- shifts, fluxes and magnifications from the literature (red squares – for de- tails, see Harris et al. 2012, from whence this figure has been adapted). The ber of parameters to conspire together if we are to escape the power-law fit to these data represents an estimate of the intrinsic line lumi- conclusion that the gas in these disks does not obey the simple nosity versus line width relation for typical SMGs. Blue circles represent the volumetric star-formation law proposed by Krumholz et al. lensed sample from Harris et al. (2012). The point for HATLAS J084933 has (2012). been replaced with those appropriate for components W and T. W has moved well within the scatter of the intrinsic relation, with a predicted amplification Before leaving the S-K and E-S laws, we note that al- consistent with unity. The predicted amplification for T is consistent with the though LIR is an effective tracer of PSFR in Eddington-limited lens model presented in §3.2, µ = 1.5 0.2. star-forming disks, Ballantyne, Armour, & Indergaard (2013) ± have argued that velocity-integrated CO line intensity is a poor proxy for Pgas. Resolved observations of the mid-IR ro- sion that the system is bound. tational lines of H2 may provide the ultimate probe of the star- Taking another approach, we searched the (500/h)3-Mpc3 formation relations in high-redshift disks, requiring a far-IR Millenium Simulation at z = 2 (Springel et al. 2005). In each space interferometer such as FIRI (Helmich & Ivison 2008). friends-of-friends halo we identified all of the sub-halos, as well as the most massive sub-halo, which we take to be 3.11. Future evolution of the system akin to W. We then computed the velocity and spatial off- What does the HATLASJ084933 system look like in the set between the most massive sub-halo and all the other sub- present day? halos and tracked the ≈ 200 cases where there are three sub- The extreme starbursts that brought this galaxy to our at- halos within a radius of 300kpc (co-moving) with velocities tention via Herschel are likely to have been triggered by > 90, > 600 and > 800kms−1 in projection (i.e. taking into one or more interactions (Engel et al. 2010). Certainly, C is account likely line-of-sight distances and transverse veloci- 13 distorted, morphologically, and the distances between these ties). Those halos, typically ≈ 1.6 × 10 M⊙ at z = 2, end galaxies are consistent with periods of intense star formation +5.5 14 up with a median mass of ≈ 3.8−1.9 × 10 M⊙ in the present during merger simulations (e.g. Springel et al. 2005). 14 day (90% lower limit, 1.6 × 10 M⊙) – structures that are However, without knowledge of the dark-matter halo they −7 −3 inhabit and the transverse velocities of W, T, C and M, it is thought to have a space density of a few ×10 Mpc (e.g. impossible to be sure whether these galaxies are gravitation- Reiprich & Böhringer 2002; Vikhlinin et al. 2009). ally bound, let alone whether they will merge. We find that The scale over which we see the galaxies comprising HAT- − LAS J084933, ≈100kpc, is similar to the core radius of rich the total energy, 1/2v2 +Φ (kms 1)2, of each of these systems clustersof galaxies,and the separationof W and T (in terms of is negative – indicating that the galaxies are bound to one an- −1 other – if their dynamical masses represent as little as 50% of kpcandkms ) is typical of binary galaxies(e.g. Turner 1976; Schweizer 1987) and within a factor 2× of the projected sep- the total. This is confirmed by N-body simulations, starting 58 from the observed velocities and positions on the sky, and as- aration of the two dominant galaxies in the Coma cluster. suming that the transverse velocities do not exceed those in We conclude that such systems can survive for a long time the line of sight. If these galaxies inhabit an expected dark- and are not necessarily in the process of merging. The four 13 starbursts in HATLASJ084933 may instead be signposting a matter halo of mass ≈ 10 M⊙ (e.g. Amblard et al. 2011; Hildebrandt et al. 2013), it is difficult to escape the conclu- proto-cluster, where the majority of members are likely well below the detection thresholds of our panchromatic imaging. 55 Changing αCO by a factor two (in opposite directions) for the SFGs and It is ironic that deliberate attempts to select clusters at submm SMGs would bring them in line with one another. Because of the potential wavelengths (targeting radio galaxies and quasars, for exam- for a large systematic error in measurements of stellar mass, amongst other ple – Stevens et al. 2003; De Breuck et al. 2004; Greve et al. things, we view such a change as plausible. 2007; Priddey et al. 2008) have met with only limited success, 56 Any change in the area used to calculate the gas or SFR surface density results in a point shifting parallel to the power-law relation seen in this plot. while serendipitous discoveries are becoming relatively com- 57 Although we have estimated the orbital timescale, in this calculation, we use the dynamical timescale of the disk, calculated using the half-light radius 58 Coma (Abell 1656) is a B(inary)-type cluster, with two bright, super- and the velocity dispersionm since Krumholz et al. (2012) show that use of giant galaxies, NGC 4874 and NGC 4889, close to each other in the cluster the orbital time can artificially increase the apparent star-formation efficiency. core (Rood & Sastry 1971). HYLIRGS AMIDST DISTANT STARBURSTING CLUSTER 13 monplace (e.g. Riechers et al. 2010). arcsecond spectral and spatial resolution, which have led to the discovery of at least four starbursting galax- ≈ 3.12. On the selection of HyLIRGs, and their rarity ies across a 100-kpc region at z =2.41, each with a significant mass of stars. Having resolved HATLAS J084933 into its constituent parts, if we place W and T individually on the Harris et al. 2. The two brightest galaxies, W and T, are separated by ′ plot of LCO versus FWHM line width (Fig. 7) we find that W ∼85kpc on the sky. W suffers no gravitational amplifi- has moved further within the scatter of the intrinsic relation, cation. T is marginally lensed by a foreground lenticu- with a predicted amplification close to unity. The predicted lar. amplification for T, µ ≈ 2, is consistent with the lens model presented in §3.2. 3. The panchromaticSEDsof W and T revealthat bothare We can estimate the space density of objects as rare as associated with high-luminosity events – HyLIRGs. If HATLASJ084933 by utilising the semi-empirical model of contributions to LIR from AGN are small, and to date Hopkinset al. (2010) and Hayward et al. (2013). We start we have found no reason to suspect otherwise, then W with halo mass functions and merger rates over the red- and T have logLIR ∼ 13.52 and 13.16, Tdust = 36–40K, −1 shift range of interest derived from the Millenium Simula- and are generating ≈3400 and ≈1500M⊙ yr of stars, tion (Springel et al. 2005), assigning galaxies to each halo respectively. following the abundance-matching techniques described in Conroy & Wechsler (2009) to yield the galaxy-galaxy merger 4. Two other gas-rich galaxies, M and C, detected rate as a function of mass, merger-mass ratio and red- nearby, make HATLASJ084933 reminiscent of the first SMG found and the best-studied HyLIRG to date, shift. Utilising the stellar mass functions of Marchesini et al. − (2009), with extrapolations to higher redshifts following SMM J02399 0136,which comprises a complex galac- Fontana et al. (2006), we calculate the merger rate for all tic nursery with three massive stellar components, each galaxies with a stellar mass ratio equal to or greater than that seen during a different evolutionary stage (Ferkinhoff of W and T, with one component at least as massive as W. et al., in preparation). The four galaxies known to com- Integrating the merger rates between z = 2–6 (t = 2.4Gyr), prise HATLASJ084933 are most likely bound. N-body where we are sensitive to such events via H-ATLAS, results simulations and comparison with the Millenium Simu- in a space density of approximately 10−7 Mpc−3, similar to lation suggest that in the present day the system may 14.6 the space density predicted for extreme HyLIRGs (Lapi et al. resemble a B(inary)-type cluster of mass, ≈ 10 M⊙. 2011, recall also §3.11). For the typical gas-consumption 5. Sub-arcsecond interferometric imaging reveals that the timescale found in §3.9, τ ≈15–45Myr, we then expect SMG two brightest galaxies span ∼3 kpc FWHM in CO J = ≈ ≈ 2 to find 3–9 such galaxies in the 100deg explored so far 4−3, and in submm and radio continuum, and roughly in H-ATLAS. The vast majority of bright Herschel galaxies double that in CO J =1−0. Alongside the detection of have been shown to be lensed (Negrello et al. 2010), but the [C II], this supports a scenario involving a widespread low end of this range is plausible. This suggests that the burst of intense star formation. most luminous phase of these galaxies has a short duration, ∼< 15Myr, and/or that such systems require some other rare 6. Exquisite 3-D imaging from JVLA and IRAM PdBI property – counter-rotating disks being one possibility. Stel- reveal counter-rotating gas disks in W and T, a sce- lar mass functions at high redshift have significant uncertain- nario that has long been predicted to lead to the most ties associated with them; nonetheless, this calculation em- intense bursts of star formation. This suggests that sim- phasises the rarity of major mergers such as that between W ilarly luminous galaxies will often be found in pairs. and T and illustrates the efficiency of the method employed Their modest velocity dispersions mean that the disks here to find extreme HyLIRGs. The observed surface den- are prone to instabilities on many scales, even that of sity of these systems is such that even with short durations for the disks themselves; they are undergoing coherent, ex- their most luminous phases, they can end up as a subset of the treme starbursts, at close to the Eddington limit, which massive current-day clusters discussed in §3.11. are not confined to their nuclei. Despite their disk- HATLASJ084933was selected on the basis of its broad CO like morphologies and dynamics, their specific star- line profile as a low-magnification, intrinsically luminous sys- formation rates place them ∼ 5× above the main se- tem, so it is ironic that its wide line profile is in part due to the quence. multiple nature of the source and that future searches for sin- gle maximal-luminosity galaxies may require the selection of 7. W and T have CO line intensities and widths typical less extreme profiles. However, fundamentally we conclude of the brightest SMGs, with slightly higher Tb ratios. that the method explored here for the selection of the most Three independent estimates of the CO-to-H2 conver- extreme starburst events in the Universe shows promise. sion factor – exploiting the dynamical mass, the mass of dust and a recent metallicity-dependent relation – 4. CONCLUSIONS −1 2 −1 suggest αCO =0.8±0.4M⊙ (Kkms pc ) and hence 11 1. Exploiting the relationship between CO luminosity and gas masses in W and T of ≈ 10 M⊙. Their gas frac- line width determined for unlensed starbursts, we have tions, ∼ 40%, are difficult to reconcile with cosmolog- identified and removed gravitationally lensed sources ical simulations. from the brightest galaxies found in the widest ex- tragalactic Herschel survey to yield a sample of in- 8. We have determined that the single-dish CO line profile trinsically luminous galaxies. Here, we report deep of HATLASJ084933 is exceptionally broad because W panchromatic follow-up observations of the best can- and T have relatively wide lines, broadened further in didate HyLIRG system, HATLASJ084933, with sub- the Harris et al. (2012) integrated GBT spectrum by the 14 IVISONETAL.

modest velocity offset between them, and by contribu- Borne, K. D., Bushouse, H., Lucas, R. A., & Colina, L. 2000, ApJ, 529, L77 tions from a number of faint CO clumps slightly red- Bothwell, M. S., Smail, I., Chapman, S. C., et al. 2013, MNRAS, 429, 3047 Bruzual, G., & Charlot, S. 2003, MNRAS, 344, 1000 ward of the two dominant components. We can thus Chabrier, G. 2003, PASP, 115, 763 reconcile the GBT CO line profile with that of the Chapman, S. C., Blain, A., Smail, I., & Ivison, R. 2005, ApJ, 622, 772 Conroy, C., & Wechsler, R. H. 2009, ApJ, 696, 620 various galaxies found to comprise HATLASJ084933. Courteau, S. 1997, AJ, 114, 2402 These situations likely happen only rarely, even in the da Cunha, E., Charlot, S., & Elbaz, D. 2008, MNRAS, 388, 1595 widest extragalactic Herschel surveys. We estimate that Daddi, E., Elbaz, D., Walter, F., et al. 2010, ApJ, 714, L118 De Breuck, C., Bertoldi, F., Carilli, C., et al. 2004, A&A, 424, 1 only a few systems of this nature will be found in every Di Matteo, P., Combes, F., Melchior, A.-L., & Semelin, B. 2007, A&A, 468, 100deg2 surveyed. 61 Dowell, C. D., Allen, C. A., Babu, R. S., et al. 2003, in Society of 9. Placing W and T individually on the Harris et al. plot Photo-Optical Instrumentation Engineers (SPIE) Conference Series, Vol. ′ 4855, Society of Photo-Optical Instrumentation Engineers (SPIE) of LCO versus FWHM CO line width, W and T move Conference Series, ed. T. G. Phillips & J. Zmuidzinas, 73–87 further within the scatter of the intrinsic relation, with a Downes, D., & Solomon, P. M. 1998, ApJ, 507, 615 Dunne, L., Eales, S., Edmunds, M., et al. 2000, MNRAS, 315, 115 predicted amplification close to unity. We conclude that Dunne, L., & Eales, S. A. 2001, MNRAS, 327, 697 the method outlined here for the selection of extreme Dwek, E., Staguhn, J. G., Arendt, R. G., et al. 2011, ApJ, 738, 36 Eales, S., Dunne, L., Clements, D., et al. 2010, PASP, 122, 499 starbursts can indeed separate HyLIRGs from the more Elbaz, D., Dickinson, M., Hwang, H. S., et al. 2011, A&A, 533, A119 numerous, less luminous, strongly lensed population of Emerson, J. P., & Sutherland, W. J. 2010, in SPIE Conference Series, Vol. IR-bright galaxies, and may also pinpoint some distant 7733 Engel, H., Tacconi, L. J., Davies, R. I., et al. 2010, ApJ, 724, 233 clusters of starbursting proto-ellipticals. Fazio, G. G., Hora, J. L., Allen, L. E., et al. 2004, ApJS, 154, 10 ′ Fontana, A., Salimbeni, S., Grazian, A., et al. 2006, A&A, 459, 745 10. Plotting LIR versus LCO, WandTlie amongsttheSMGs Förster Schreiber, N. M., Genzel, R., Bouché, N., et al. 2009, ApJ, 706, 1364 on a single relation that is consistent with all star- Frayer, D. T., Harris, A. I., Baker, A. J., et al. 2011, ApJ, 726, L22 Gear, W. K., Lilly, S. J., Stevens, J. A., et al. 2000, MNRAS, 316, L51 forming galaxies. They also lie amongst the SMGs Genzel, R., Tacconi, L. J., Gracia-Carpio, J., et al. 2010, MNRAS, 407, 2091 on the classic Schmidt-Kennicutt and Elmegreen-Silk Goldreich, P., & Lynden-Bell, D. 1965, MNRAS, 130, 97 Greve, T. R., Stern, D., Ivison, R. J., et al. 2007, MNRAS, 382, 48 plots, but the SMGs now follow a relation that is dis- Griffin, M. J., Abergel, A., Abreu, A., et al. 2010, A&A, 518, L3+ tinct from those of other star-forming galaxies. Their Hainline, L. J., Blain, A. W., Smail, I., et al. 2011, ApJ, 740, 96 star-formation efficiency is significantly higher, and is Harris, A. I., Baker, A. J., Zonak, S. G., et al. 2010, ApJ, 723, 1139 Harris, A. I., Baker, A. J., Frayer, D. T., et al. 2012, ApJ, 752, 152 difficult to reconcile with a simple volumetric star- Hayward, C. C., Kereš, D., Jonsson, P., et al. 2011, ApJ, 743, 159 formation law in which ∼1% of the available gas is Hayward, C. C., Narayanan, D., Kereš, D., et al. 2013, MNRAS, 428, 2529 converted into stars during each local free-fall time. Heiderman, A., Evans, II, N. J., Allen, L. E., Huard, T., & Heyer, M. 2010, ApJ, 723, 1019 Helmich, F. P., & Ivison, R. J. 2008, Experimental Astronomy, 15 Understanding systems as luminous as HATLASJ084933 Helou, G., Soifer, B. T., & Rowan-Robinson, M. 1985, ApJ, 298, L7 Hildebrandt, H., van Waerbeke, L., Scott, D., et al. 2013, MNRAS, 488 represents a challenge. However, the observations presented Hopkins, P. F., Quataert, E., & Murray, N. 2012, MNRAS, 421, 3488 here have shed considerable light on the characteristics that Hopkins, P. F., Younger, J. D., Hayward, C. C., Narayanan, D., & Hernquist, likely lead to the formation – however briefly – of the most L. 2010, MNRAS, 402, 1693 Ibar, E., Ivison, R. J., Best, P. N., et al. 2010a, MNRAS, 401, L53 luminous star-forming galaxies in the Universe. Observa- Ibar, E., Ivison, R. J., Cava, A., et al. 2010b, MNRAS, 409, 38 tions scheduled for Cycle 1 with the Atacama Large Millime- Ilbert, O., Salvato, M., Le Floc’h, E., et al. 2010, ApJ, 709, 644 ter/Submillimeter Array promise to reveal more about this Ilbert, O., McCracken, H. J., Le Fevre, O., et al. 2013, ArXiv e-prints, arXiv:1301.3157 dramatic and complex proto-galactic environment, probing Iono, D., Peck, A. B., Pope, A., et al. 2006, ApJ, 640, L1 the structure of the disks and sensitive to the presence and Ivison, R. J., Papadopoulos, P. P., Smail, I., et al. 2011, MNRAS, 412, 1913 Kennicutt, Jr., R. C. 1989, ApJ, 344, 685 influence of buried AGN. —. 1998, ApJ, 498, 541 Klaas, U., Nielbock, M., Haas, M., Krause, O., & Schreiber, J. 2010, A&A, ACKNOWLEDGEMENTS 518, L44 Kovács, A. 2008, in SPIE Conference Series, Vol. 7020 RJI acknowledges support from the European Research Kovács, A., Omont, A., Beelen, A., et al. 2010, ApJ, 717, 29 COSMICISM Krumholz, M. R., Dekel, A., & McKee, C. F. 2012, ApJ, 745, 69 Council (ERC) in the form of Advanced Grant, . Lada, C. J., Lombardi, M., & Alves, J. F. 2010, ApJ, 724, 687 The Herschel-ATLAS is a project with Herschel, which is Lapi, A., González-Nuevo, J., Fan, L., et al. 2011, ApJ, 742, 24 an ESA space observatory with science instruments provided Leong, M., Peng, R., Houde, M., et al. 2006, in Society of Photo-Optical Instrumentation Engineers (SPIE) Conference Series, Vol. 6275, Society by European-led Principal Investigator consortia and with of Photo-Optical Instrumentation Engineers (SPIE) Conference Series important participation from NASA. The H-ATLAS website Lutz, D., Dunlop, J. S., Almaini, O., et al. 2001, A&A, 378, 70 is www.h-atlas.org. US participants in H-ATLAS ac- Magnelli, B., Lutz, D., Santini, P., et al. 2012, A&A, 539, A155 Makovoz, D., & Marleau, F. R. 2005, PASP, 117, 1113 knowledge support from NASA through a contract from JPL. Marchesini, D., van Dokkum, P. G., Förster Schreiber, N. M., et al. 2009, IRS acknowledges support from STFC and ERC. MN and ApJ, 701, 1765 Michalowski, M. J., Dunlop, J. S., Cirasuolo, M., et al. 2012, A&A, 541, GdZ acknowledge financial support from ASI/INAF agree- A85 ment I/072/09/0. Mihos, J. C., & Hernquist, L. 1994, ApJ, 431, L9 Facilities: CARMA, CSO, Herschel, Hubble Space Tele- —. 1996, ApJ, 464, 641 Morton, D. C., & Noreau, L. 1994, ApJS, 95, 301 scope, IRAM PdBI, JVLA, Keck, SMA, Spitzer, VISTA. Narayanan, D., Bothwell, M., & Davé, R. 2012a, MNRAS, 426, 1178 Narayanan, D., Krumholz, M. R., Ostriker, E. C., & Hernquist, L. 2012b, REFERENCES MNRAS, 421, 3127 Negrello, M., González-Nuevo, J., Magliocchetti, M., et al. 2005, MNRAS, 358, 869 Alexander, D. M., Brandt, W. N., Smail, I., et al. 2008, AJ, 135, 1968 Negrello, M., Hopwood, R., De Zotti, G., et al. 2010, Science, 330, 800 Amblard, A., Cooray, A., Serra, P., et al. 2011, Nature, 470, 510 Owen, F. N., & Morrison, G. E. 2008, AJ, 136, 1889 Aravena, M., Younger, J. D., Fazio, G. G., et al. 2010, ApJ, 719, L15 Pascale, E., Auld, R., Dariush, A., et al. 2011, MNRAS, 415, 911 Ballantyne, D. R., Armour, J. N., & Indergaard, J. 2013, ApJ, 765, 138 Pforr, J., Maraston, C., & Tonini, C. 2012, MNRAS, 422, 3285 Bekki, K. 2001, ApJ, 546, 189 Pilbratt, G. L., Riedinger, J. R., Passvogel, T., et al. 2010, A&A, 518, L1+ Bolatto, A. D., Wolfire, M., & Leroy, A. K. 2013, ARA&A, Poglitsch, A., Waelkens, C., Geis, N., et al. 2010, A&A, 518, L2+ arXiv:1301.3498 HYLIRGS AMIDST DISTANT STARBURSTING CLUSTER 15

Pope, A., Borys, C., Scott, D., et al. 2005, MNRAS, 358, 149 Smolciˇ c,´ V., Navarrete, F., Aravena, M., et al. 2012b, ApJS, 200, 10 Priddey, R. S., Ivison, R. J., & Isaak, K. G. 2008, MNRAS, 383, 289 Springel, V., White, S. D. M., Jenkins, A., et al. 2005, Nature, 435, 629 Reiprich, T. H., & Böhringer, H. 2002, ApJ, 567, 716 Stacey, G. J., Hailey-Dunsheath, S., Ferkinhoff, C., et al. 2010, ApJ, 724, Riechers, D. A., Hodge, J., Walter, F., Carilli, C. L., & Bertoldi, F. 2011, 957 ApJ, 739, L31 Stevens, J. A., Ivison, R. J., Dunlop, J. S., et al. 2003, Nature, 425, 264 Riechers, D. A., Capak, P. L., Carilli, C. L., et al. 2010, ApJ, 720, L131 Swinbank, A. M., Smail, I., Longmore, S., et al. 2010, Nature, 464, 733 Rigby, E. E., Maddox, S. J., Dunne, L., et al. 2011, MNRAS, 415, 2336 Taniguchi, Y., & Shioya, Y. 1998, ApJ, 501, L167 Rood, H. J., & Sastry, G. N. 1971, PASP, 83, 313 Toomre, A. 1964, ApJ, 139, 1217 Rowlands, K., Dunne, L., Maddox, S., et al. 2012, MNRAS, 419, 2545 Turner, E. L. 1976, ApJ, 208, 20 Salomé, P., Guélin, M., Downes, D., et al. 2012, A&A, 545, A57 Vikhlinin, A., Kravtsov, A. V., Burenin, R. A., et al. 2009, ApJ, 692, 1060 Sandstrom, K. M., Leroy, A. K., Walter, F., et al. 2013, ApJ, Wang, W.-H., Cowie, L. L., van Saders, J., Barger, A. J., & Williams, J. P. arXiv:1212.1208 2007, ApJ, 670, L89 Schmidt, M. 1959, ApJ, 129, 243 Werner, M. W., Roellig, T. L., Low, F. J., et al. 2004, ApJS, 154, 1 Schweizer, L. Y. 1987, ApJS, 64, 411 Whitmore, B. C., & Schweizer, F. 1995, AJ, 109, 960 Scoville, N. 2004, in Astronomical Society of the Pacific Conference Series, Wilson, C. D., Petitpas, G. R., Iono, D., et al. 2008, ApJS, 178, 189 Vol. 320, The Neutral ISM in Starburst Galaxies, ed. S. Aalto, Younger, J. D., Fazio, G. G., Huang, J.-S., et al. 2007, ApJ, 671, 1531 S. Huttemeister, & A. Pedlar, 253 Younger, J. D., Fazio, G. G., Wilner, D. J., et al. 2008, ApJ, 688, 59 Scoville, N. Z., Sanders, D. B., & Clemens, D. P. 1986, ApJ, 310, L77 Younger, J. D., Fazio, G. G., Huang, J.-S., et al. 2009, ApJ, 704, 803 Smolciˇ c,´ V., Aravena, M., Navarrete, F., et al. 2012a, A&A, 548, A4 Yun, M. S., Reddy, N. A., & Condon, J. J. 2001, ApJ, 554, 803