Draft version June 20, 2018 Preprint typeset using LATEX style emulateapj v. 5/2/11 GLOBAL STAR FORMATION RATES AND DUST EMISSION OVER THE GALAXY INTERACTION SEQUENCE Lauranne Lanz1, Andreas Zezas1,2,3, Nicola Brassington4, Howard A. Smith1, Matthew L. N. Ashby1, Elisabete da Cunha5, Giovanni G. Fazio1, Christopher C. Hayward6, Lars Hernquist1, Patrik Jonsson7 Draft version June 20, 2018 ABSTRACT We measured and modeled the spectral energy distributions (SEDs) in 28 bands from the ultraviolet to the far-infrared (FIR) for 31 interacting galaxies in 14 systems. The sample is drawn from the Spitzer Interacting Galaxy Survey, which probes a range of galaxy interaction parameters at multiple wavelengths with an emphasis on the infrared bands. The subset presented in this paper consists of all galaxies for which FIR Herschel SPIRE observations are publicly available. Our SEDs combine the Herschel photometry with multi-wavelength data from Spitzer, GALEX, Swift UVOT, and 2MASS. While the shapes of the SEDs are broadly similar across our sample, strongly interacting galaxies typically have more mid-infrared emission relative to their near-infrared and FIR emission than weakly or moderately interacting galaxies. We modeled the full SEDs to derive host galaxy star formation rates (SFR), specific star formation rates (sSFR), stellar masses, dust temperatures, dust luminosities, and dust masses. We find increases in the dust luminosity and mass, SFR, and cold (15-25 K) dust temperature as the interaction progresses from moderately to strongly interacting and between non- interacting and strongly interacting galaxies. We also find increases in the SFR between weakly and strongly interacting galaxies. In contrast, the sSFR remains unchanged across all the interaction stages. The ultraviolet photometry is crucial for constraining the age of the stellar population and the SFR, while dust mass is primarily determined by SPIRE photometry. The SFR derived from the SED modeling agrees well with rates estimated by proportionality relations that depend on infrared emission. Subject headings: infrared: galaxies, galaxies: interactions, galaxies: photometry, galaxies: star for- mation, ultraviolet: galaxies 1. INTRODUCTION vatory8, are thought to be predominantly mergers (e.g., Galaxy evolution is believed to be heavily influenced Blain et al. 1999), although the relative contribution of by interactions between galaxies, both for local systems mergers of different stages to their numbers is still an and for distant objects at earlier cosmological times. In open question (e.g., Hayward et al. 2012a, 2012b). the canonical view, interactions between galaxies have The second effect is that interactions significantly af- three primary observable effects. In the most dramatic fect the subsequent evolution of galaxies, which may lead cases, interactions stimulate star formation in a burst of to significant changes in their morphology. Disturbed activity that is presumed to power the high infrared (IR) galaxies have long been associated with mergers (e.g., luminosities typically seen in such systems. Many local Toomre & Toomre 1972). Numerical simulations of in- ultra-luminous IR galaxies (L ≥ 1012 L ; ULIRGs) and teractions (e.g., Hopkins et al. 2006; Hopkins 2012; Mi- hos & Hernquist 1994, 1996; Barnes 1992; Barnes & luminous IR galaxies (1011 L ≤ L ≤ 1012 L ; LIRGs) Hernquist 1996; Sanders 1999) show a variety of mor- show evidence of galaxy interactions (e.g., Veilleux et phological distortions as well as variable amounts of star al. 2002). Similarly, their high-redshift counterparts, formation. The simulations also demonstrate the com- sub-millimeter galaxies, first detected by SCUBA and plexity of the problem: the degree of induced activity now studied extensively by the Spectral and Photometric and distortion varies greatly with the parameters of the Imaging Receiver (SPIRE) on the Herschel Space Obser- arXiv:1302.5011v1 [astro-ph.CO] 20 Feb 2013 encounter, the phase of the interactions, the molecular gas content (\wetness"), and the mass of the progenitor galaxies among many other properties. Third, the canonical picture, as seen in many simula- 1 Harvard-Smithsonian Center for Astrophysics, 60 Garden St., Cambridge, MA 02138, USA; [email protected] tions (e.g., Di Matteo et al. 2005, Springel et al. 2005), 2 University of Crete, Physics Department & Institute of Theo- involves merger-driven gas inflow to the central regions, retical & Computational Physics, 71003 Heraklion, Crete, Greece resulting in heightened activity of the central supermas- 3 Foundation for Research and Technology-Hellas, 71110 Her- sive black hole as well as starburst activity due to the in- aklion, Crete, Greece 4 School of Physics, Astronomy and Mathematics, University creased central gas density and possibly turbulence. The of Hertfordshire, College Lane, Hatfield, AL10 9AB, UK process in principle converts a low-luminosity nucleus 5 Max Planck Institute for Astronomy (MPIA), K¨onigstuhl17, into an active galactic nucleus (AGN) but one whose lu- 69117, Heidelberg, Germany 6 Heidelberger Institut f¨ur Theoretische Studien, Schloss- Wolfsbrunnenweg 35, 69118, Heidelberg, Germany 8 Herschel is an ESA space observatory with science instruments 7 Space Exploration Technologies, 1 Rocket Road, Hawthorne, provided by European-led Principal Investigator consortia and with CA, 90250 important participation from NASA. 2 Lanz et al. minosity might range widely depending on the stage of (not least) previous interaction histories (e.g., Di Mat- the interaction. Indeed, observations of merging galaxies teo et al. 2007). Systems are ordered into an evolution- over the years have tended to provide evidence support- ary sequence using intuition provided by simulations and ing the conclusion that, at least on a statistical level, physical models, which are themselves based on observa- interactions trigger an enhancement in the formation of tions of particular systems. Selection criteria, however, stars as well as nuclear activity. However, the recent lit- can introduce a bias for more luminous, morphologically erature includes works that argue both for and against a disturbed systems and, hence, towards the most active strong connection between nuclear activity and mergers phases of interactions. Therefore, a selection criterion (e.g., Li et al. 2008, Kocevski et al. 2012, Ellison et al. not associated with either morphological disturbance or 2011, Scudder et al. 2012, Silverman et al. 2011). degree of activity is critical for obtaining a sample con- Therefore, despite many previous studies (e.g., Dahari taining systems throughout the full interaction sequence. 1985, Sanders et al. 1988, Kewley et al. 2001, Lambas et We have undertaken a program to take advantage of al. 2003), both observational and through simulations, all these developments: full, multi-wavelength datasets our understanding of the evolution of the physical ac- of an interacting galaxy sample selected with few biases; tivity during the course of a galaxy-galaxy interaction hydrodynamic simulations; and radiative transfer model- remains incomplete. In the past decade, two new devel- ing, in a systematic effort to better understand systems opments have dramatically changed our understanding across a range of interaction stages and to iterate re- of star formation and accretion activity around galactic finements to the various modeling and radiative transfer nuclei, which are the two dominant processes at work in codes. We have chosen a representative sample of ob- controlling the observed emission. The first is the suc- jects spanning the interaction sequence, obtained their cess of space missions, in particular, the Spitzer Space full spectral energy distributions (SEDs), and are com- Telescope9 (Werner et al. 2004) and the Herschel Space paring the results against a variety of models - based on Observatory (Pilbratt et al. 2010) in the IR, as well as both templates/stellar evolution/radiative transfer and the Galaxy Evolution Explorer 10 (GALEX ; Martin et on diagnostic features. al. 2005) and Swift (Gehrels et al. 2004) in the ultraviolet This first paper of the project presents results and con- (UV), providing photometry across the complete spectral clusions for a sample of thirty-one interacting galaxies range from UV to millimeter wavelengths. Most impor- in fourteen systems for which there are currently com- tantly, detailed imaging and high sensitivity photometry plete multi-wavelength data that can be used to study now available at the critical far-infrared (FIR) emission the variations in their star formation and dust heating. peak resulting from warm, luminous dust heated by star- This paper is organized as follows. We describe the full bursts provides crucial information regarding dust heat- Spitzer Interacting Galaxy Survey (SIGS) sample in Sec- ing and embedded star formation. The combination of tion 2 and the classification of each of the sources in the UV and IR observations is essential to obtain a complete interaction sequence. Section 3 describes the wide range census of recent and ongoing star formation by capturing of observational photometry used to construct the SEDs. both the unobscured and obscured emission from young It also describes the issues associated with obtaining reli- stars. able photometry from the diverse datasets. In Section 4, The second development has been the success of com- we model the SEDs of these objects. Section 5 discusses putational codes. We have new tools for the derivation