arXiv:astro-ph/9911045v1 3 Nov 1999 ihRslto i-nrrdIaigo lrlmnu Inf Ultraluminous of Imaging Mid-Infrared Resolution High aoa bevtr,Clfri nttt fTcnlg,320-47 Technology, of Institute California Observatory, Palomar IT cec etr aionaIsiueo ehooy 314-6, Technology, of Institute California Center, Science SIRTF eateto hsc n srnm,Uiest fClfri Los California of University Astronomy, and Physics of Department iio fPyis ahmtc n srnm,Clfri Institut California Astronomy, and Mathematics Physics, of Divison [email protected],gnmpcleheu [email protected] [email protected], [email protected], e rplinLb 6-0,40 a rv r,Psdn,C 91 CA Pasadena, Dr., Grove Oak 4800 169-506, Lab, Propulsion Jet [email protected],[email protected] .T Soifer T. B. [email protected],[email protected] [email protected], [email protected],[email protected] .E eki,A .Weinberger J. A. Becklin, E. E. 2 .Nueae,K ates .Egami E. Matthews, K. Neugebauer, G. , 0-4 aaea A91125 CA Pasadena, 105-24, .S Evans S. A. [email protected] .Rslr .W Werner W. M. Ressler, M. [email protected] o nee,C 90095 CA Angeles, Los Galaxies ..Surace J.A. 3 and ..Scoville N.Z. , 1 aaea A91125 CA Pasadena, , aaea A91125 CA Pasadena, fTechnology, of e nee,156205 Angeles, tech.edu, 109 rared –2–

J.J. Condon

National Radio Astronomy Observatory, 520 Edgemont Road, Charlottesville, VA 22903

[email protected]

Received ; accepted

1Based in part on observations obtained at the W. M. Keck Observatory which is operated as a scientific partnership among the California Institute of Technology, the University of California and the National Aeronautics and Space Administration. 2also at SIRTF Science Center, California Institute of Technology, 314-6, Pasadena, CA 91125 3current address: Department of Physics and Astronomy, SUNY Stony Brook, Stony Brook NY 11794-3800 –3–

ABSTRACT

Observations of ultraluminous infrared galaxies (ULIRGs) with an achieved resolution approaching the diffraction limit in the mid-infrared from 8 - 25 µm using the Keck Telescopes are reported. We find extremely compact structures, with spatial scales of < 0.3′′ (diameter) in six of the seven ULIRGs observed. These compact sources emit between 30% and 100% of the mid-infrared energy from these galaxies. We have utilized the compact mid-infrared structures as a diagnostic of whether an AGN or a compact (100 – 300 pc) starburst is the primary power source in these ULIRGs. In Markarian 231, the upper limit on the diameter of the 12.5 µm source , 0.13′′, shows that the size of the infrared source must increase with increasing wavelength, consistent with AGN models. In IRAS 05189-2524 and IRAS 08572+3915 there is strong evidence that the source size increases with increasing wavelength. This suggests heating by a central source rather than an extended luminosity source, consistent with the optical classification as an AGN. The compact mid-infrared sources seen in the other galaxies cannot be used to distinguish the ultimate luminosity source. If these ULIRGs are powered by compact starbursts, the star formation rates seen in the central few hundred parsecs far exceed the global rates seen in nearby starburst galaxies, and approach the surface brightness of individual clusters in nearby starburst galaxies.

Subject headings: luminous infrared galaxies, infrared, galaxies individual: Arp 220; Markarian 231; Markarian 273; UGC 5101; IRAS 05189-2524; IRAS 08572+3915; IRAS 17208-0014 –4–

1. Introduction

Ultraluminous Infrared Galaxies (ULIRGs; defined as those systems with luminosity

12 from 8 - 1000 µm L − 10 L⊙) were discovered to be a significant class of objects 8 1000µm ≥ in the local Universe in the IRAS all sky survey (Soifer et al. 1987, Sanders et al. 1988). These systems have recently been suggested as a major component in the Universe at z> 2 through deep submillimeter surveys ( c.f. Lilly et al. 1999, Blain et al. 1999, Barger et al. 1998).

Whether ULIRGs are powered by dust enshrouded AGN or starbursts has been debated since their discovery(e.g. Sanders 1999, Joseph 1999). This question is made all the more important by the realization that ULIRGs could produce a major fraction of the radiant energy in the Universe (Lilly et al. 1999, Heckmann 1998).

A wide variety of observational techniques have been employed to probe the nature of the central luminosity source in ULIRGS, ranging from optical spectroscopy (Sanders et al. 1988) to cm radio continuum imaging (Condon et al. 1991, hereafter C91), to VLBI observations (Smith, Lonsdale and Lonsdale 1998) and hard X-ray observations (Rieke 1988, Nakagawa et al. 1998, Isawawa 1999). The nuclear environments of these galaxies have been effectively imaged with high spatial resolution in the optical and near infrared (Surace et al. 1998, Surace and Sanders 1999, Scoville et al. 1999, hereafter S99).

Observations in the thermal infrared have the potential to address the question of the origin of the luminosity by establishing the size of the regions that are radiating the infrared luminosity. The IRAS observations showed that the mid and far infrared luminosity in these systems emerging predominantly between 50 and 200 µm was due to thermal reradiation by dust. However, the spatial resolution of those observations, 1′ 2′, was inadequate ∼ − to place meaningful size constraints on the luminosity sources. Because the mid-infrared wavelengths (8–25µm) carry a significant fraction of the total bolometric luminosity in –5–

ULIRGs, ranging from 10% to > 30% of the 8–1000µm luminosity, observations at these ∼ wavelengths with a resolution 1′′ have the potential to constrain the sizes of the emitting ≤ regions and thereby address the nature of the underlying sources that heat the radiating dust. Observations in the 10µm atmospheric window, with a resolution 1′′, have been ∼ reported for several ULIRGs by Miles et al. (1996), and Keto et al. (1992). In most cases they found the sources to be unresolved implying that the sources were either very compact starbursts or dust enshrouded AGN.

With a diffraction limit of 0.24′′ FWHM at 10µm, the Keck Telescope provides a substantial improvement in spatial resolution over previous mid-infrared observations, probing the distribution of the thermal emission at the 100-300 pc scale in nearby ULIRGs. With the introduction of mid-infrared imaging on the Keck Telescope, we have started a program to address this problem. In the first paper from this program we reported imaging of the nucleus of the closest ULIRG - Arp 220 (Soifer et al. 1999). The observations were able to spatially resolve the emission at 24.5µm, and demonstrated that the central thermal source is optically thick at this wavelength. These observations placed severe constraints on the surface brightness in the infrared luminous nuclei.

In this paper we report imaging observations from 8 - 24.5µm of a sample of the closest ULIRGs at spatial resolutions of 0.3-0.6′′. These data provide the highest spatial resolution yet achieved on the thermal emission from the closest ULIRGs and trace the

spatial distribution of the emergent luminosity in these systems. We adopt Ho =75 km s−1Mpc−1. –6–

2. The Sample

The objects observed were taken from the IRAS Bright Galaxy Sample (Soifer et al. 1987) with the addition of IRAS 17208-0014 (Sanders et al. 1995), which is at a low Galactic latitude. The basic information for the seven objects in the sample is given in Table 1. The objects selected for observation were chosen on the basis of 1) meeting the

12 luminosity definition of the ULIRGs, i.e. L − 10 L⊙, 2) having a 25µm flux 8 1000µm ≥ density > 1 Jy based on IRAS measurements, 3) being at the lowest possible redshift for the best linear resolution, and 4) being available at the time of the observations. This sample contains the five closest ULIRGs in the IRAS Bright Galaxy sample of Soifer et al., and seven of the ten ULIRGs known with δ 30 deg and z 0.06. The projected linear scales ≥ − ≤ of these objects range from 360 pc arcsec−1 to 1200 pc arcsec−1, so that the point spread function/angular resolution of the Keck observations, 0.3′′ at 12µm and 0.6′′ at 24.5µm provide linear resolutions of 100 pc - 700 pc, depending on the object and the wavelength observed. Three of the objects in the sample have been classified as Seyfert nuclei, three as LINERS, and one as an HII region, based on optical spectroscopy.

3. Observations and Data Reduction

The observations were made primarily using the MIRLIN mid-infrared camera (Ressler et al. 1994) at the f/40 bent Cassegrain visitor port of the Keck II Telescope. The camera uses a 128 128 Si:As array with a plate scale of 0.138′′/pixel for a total field of view of × 17′′ 17′′. Table 2 reports the central wavelengths for the filters through which each object × was observed. At each wavelength the observing procedure was the same. A secondary with a square wave chop of amplitude 6′′ in the north-south or east-west direction at 4 Hz was employed for fast beam switching. The frames sampling each chop position were coadded separately in hardware, resulting in two images. After an interval of approximately –7– a minute, the telescope was nodded perpendicular to the chop direction (east-west or north-south) by 6′′ and a second pair of images was obtained in order to cancel residuals in the sky and to subtract telescope emission. This procedure was repeated a number of times at each wavelength. The data were reduced by differencing the two images obtained within the chop pairs at each nod location, and then coadding the resulting positive images, with the positions appropriately adjusted to a common location, to yield a positive image centered in a field approximately 6′′ 6′′. Because of the chopper and telescope nod spacings × employed for the observations, the data are not capable of measuring flux outside a 6′′ diameter region.

Observations with MIRLIN of all the ULIRGs except IRAS 05189-2524 were obtained in two nights in March 1998. The observations of the targets were interleaved with observations of nearby bright stars that were used to establish the point spread function (PSF) for the observations. The PSF stars were observed only three times per night in each of the bands at 12.5, 17.9 and 24.5 µm, but an attempt was made to accompany the critical observations with measurements of a nearby PSF star. At 24.5 µm the PSF was always equal to the diffraction limit of the telescope to be 0.62′′ FWHM. At the shorter wavelengths, the size of the PSF image was set by diffraction and the pixel sampling about half of the time; at other times atmospheric seeing affected the image size.

IRAS 05189-2524 was observed in October 1998. During this run the PSF stars were observed more frequently, but for wavelengths shorter than λ 17.9 µm significant ∼ variation in the measured size of the PSF star was seen.

The variation in the measured size of the PSF calibration star is the largest uncertainty in the measurement of the source sizes; this issue is addressed when observations of a particular object are described.

The MIRLIN observations were made under photometric conditions. The photometry –8– was calibrated based on observations of four bright stars, α Tau, α Boo, β Peg, and α Cet whose magnitudes, in turn, were based on IRAS photometry. The uncertainties in the photometry, based on the internal consistency of the observations, is estimated to be 5% at λ 17.9 µm and 10% at 24.5 µm. The flux density corresponding to 0.0 mag (Vega-based) ≤ was taken to follow the prescription given in the Explanatory Supplement to the IRAS Catalogs and Atlases (Joint IRAS Science Team, 1989), and is given in Soifer et al. (1999).

In addition to the MIRLIN observations, UGC 5101, Markarian 231 and IRAS 08572+3915 were observed using the Long Wavelength Spectrograph on the Keck I Telescope on the night of 01 April 1999. The night was not photometric and so only measurements of the sizes of the compact sources were made; all the observations were carried out at 12.5 µm. The chopper was set to an amplitude of 5′′ chopping north-south at a frequency 5 Hz. Observations were made in a fashion similar to the MIRLIN ∼ observations, except that the telescope nodding was in the same direction as the chopping and the nodding amplitude was the same as the chopping amplitude. The data were reduced in a manner similar to the MIRLIN data. The pixel scale is 0.08′′ pixel−1.

Observations of Markarian 273 were made at 3.45 µm using the near-infrared camera at the Cassegrain focus of the 200-inch Hale Telescope at Palomar Observatory in May 1999. The camera was read out in a 64 64 pixel subarray because of the high read rate × required to avoid saturation. The scale is 0.125′′ pixel−1. The chopper throw was 15′′ at a frequency of 1 Hz. Observations of stars given by Elias et al. (1982) were used to establish ∼ the photometric scale. –9–

4. Results

The mid-infrared observations of the ULIRGs reported here show a variety of characteristics, generally consistent with their properties seen at high spatial resolution at radio wavelengths (c.f. C91) and in the near infrared using NICMOS on HST (Scoville et al. 1998, S99). The mid-infrared data on all the ULIRGs are presented in Figures 1, 2 and 3. In Figure 1 we present the 12.5 µm images of the ULIRGs and, where the radio and near infrared data are available, compare these maps to the equivalent maps at 8.4 GHz from the data of C91 and to the NICMOS 2.2 µm data from S99 and Scoville et al (1998). In the case of three galaxies, IRAS 05189-2524, IRAS 08572+3915 and Markarian 231, where there is a single compact nucleus that emits a substantial fraction of the 12.5µm flux density, we have presented the corresponding PSF star at 12.5µm for comparison. In Figure 2 we compare, for the five single nucleus sources, the integrated flux at 12.5 µm as a function of beam radius with the same quantity for a corresponding point source to establish the spatial distribution of the mid-infrared radiation. Figure 3 displays the spectral energy distributions (SEDs) obtained from the mid-infrared data and compares these data with the corresponding IRAS data and data at other wavelengths. The flux densities measured for the ULIRGs in a 4′′ diameter beam are presented in Table 3.

The mid-infrared observations did not establish accurate absolute positions for the sources. In some cases, described below, we use arguments regarding the similarities of the structures seen at vastly different wavelengths to argue for the physical coincidence of the sources seen at different wavelengths.

In the following sections we give results for the individual objects. – 10 –

4.1. IRAS 05189 2524 −

MIRLIN observations of IRAS 05189 2524 were made at 12.5 and 24.5 µm. The − shorter wavelength data have the higher signal-to-noise ratio, and so are presented in Figure 1. Images at 12.5µm (MIRLIN), 2.2µm (NICMOS, S99) and 8.4GHz (VLA, C91) are shown in Figure 1a. The dominant central source is clear in all three maps. At 8.4GHz, the single source is concentrated with a measured size of 0.2′′ 0.17′′(C91). At 2.2 µm, the unresolved × source is < 0.2′′ and contains 75% of the total flux measured in a 5′′ diameter beam, and ∼ 70% of the flux contained in an 11′′ diameter beam (S99).

The source appears virtually unresolved in the MIRLIN image. As can be seen from Figure 2, there is excellent agreement between the curve of growth for IRAS 05189-2524 and the PSF star observation. The apparent size of the PSF star varied during the time of the observations, so the limit on the size of IRAS 05189-2524 is not as stringent as would otherwise be achieved. At 12.5µm the size is <0.2′′ (Table 4). The source size at 24.5 µm is an upper limit as well, and is larger and thus consistent with that established at 12.5 µm.

As seen in Figure 3 and Table 3, the flux densities measured by MIRLIN at 12.5 µm and 24.5 µm agree well with the IRAS observations of this source; virtually all of the luminosity of this source at 12.5 µm and 24.5 µm is contained within the unresolved source measured here.

The MIRLIN observations of IRAS 05189 2524 lead to a limit on the source diameter − of <170 pc in the mid-infrared (Table 5). The combination of flux density and upper limit on the source sizes(Table 4), imply the brightness temperatures Tb >123 K and >85 K at 12.5 µm and 24.5 µm respectively. The flux ratio between 12.5 µm and 24.5 µm corresponds to a color temperature Tc =169 K, while the color temperature derived from the 25µm and – 11 –

60µm IRAS measurements is 83 K1.

A central source with the bolometric luminosity of IRAS 05189 2524 would heat − silicate or graphite grains to the apparent color temperature between 12 and 24.5µm, T c ∼ 170 K, at a radius of 30-60 pc using the grain parameters of Draine and Lee (1984). This is much smaller than the observational limit. Alternatively, the size of an optically thin source which produces the 25 60 µm flux ratio with a color temperature of 83 K is significantly − larger, 170-250 pc in radius, again using the Draine and Lee grain parameters.

4.2. IRAS 08572+3915

MIRLIN observations of IRAS 08572+3915 were obtained at five wavelengths (Table 2). In addition, LWS observations of IRAS 08572+3915 were obtained at 12.5 µm. Optical and near infrared imaging (Sanders et al. 1988, Carico et al. 1990, S99) show two galaxies separated by about 5′′ in a northwest-southeast orientation. The southeast source is not detected in the radio (C91) or in the present data.

The 12.5µm LWS image in Figure 1b shows a bright central source (the northwest source in the near infrared) that is unresolved. The size limit at 17.9 µm is consistent with that found at 12.5 µm. Figure 1b and Figure 2 show that the central source is dominant in the mid-infrared and radio, containing 80% of the total flux measured in a 4′′ diameter ∼ beam. At 2.2 µm, the central unresolved source contains 60% of the total flux in an 11.4′′ ∼ diameter beam (S99).

1The brightness temperature, T , is defined by f = Ω B (T ) where B (T ) is the b ν × ν b ν b blackbody emission at frequency ν for temperature Tb and Ω is the solid angle subtended by the source. The color temperature of a source, Tc, measured between frequencies ν1 and

fν1 Bν1 (Tc) ν2 is found from = . fν2 Bν2 (Tc) – 12 –

The most stringent observations constraining the size and spatial distribution of the source are the LWS observations at 12.5µm. Five observations of the PSF star (BS3275) were interleaved with eight observations of the object. The average FWHM of the object (0.39 0.04′′) was marginally greater than the average FWHM of the PSF (0.34 0.01′′). ± ± Formally this leads to a size of the 12.5µm source of 0.19′′, but this is only a 1σ result. A 3σ limit on the size of the compact source is 0.3′′. As illustrated in Figure 2, there appears to ∼ be a measureable flux, approximately 20% of the total at 12.5 µm, that extends beyond the central core, and is contained within 4′′. The figure shows clearly that the flux measured in IRAS 08572+3915 continues to increase beyond the beam where the flux in the PSF star has stopped increasing.

Figure 3 and Table 3 compare the Keck photometry in a 4′′ diameter beam for IRAS 08572+3915 with the IRAS measurements. The MIRLIN photometry accounts for 70% of the IRAS flux density at 12.5µm. The MIRLIN observations were taken with comparatively narrow filters, i.e., ∆λ/λ < 0.1, from 7.9µm to 17.9µm, while the IRAS ∼ 12µm filter bandpass extends from 8-15µm. Given the uncertainties in the corrections for the different bandpasses, the MIRLIN data are consistent with all the flux detected by IRAS being confined to a 4′′ diameter beam surrounding the detected source although it is possible that some of the flux could be accounted for by extended emission associated with the northwest source or from the southeastern source.

Neither the MIRLIN nor the LWS observations detected the galaxy approximately 5′′ southeast of the object centered in Figure 1b. The MIRLIN data limit the contribution of the southeastern source to <30% of the total 12.5µm flux; it is likely that the southeast galaxy does not contribute significantly to the total infrared luminosity in this system. This is consistent with the optical and near infrared morphology and colors of this galaxy, which are indicative of quiescent galaxy nucleus (Surace et al. 1998, Surace and Sanders 1999). – 13 –

4.3. UGC 5101 = IRAS 09321+6134

MIRLIN observations of UGC 5101 were obtained at 12.5 and 17.9 µm. Contour plots of this source in Figure 1c show that at all wavelengths the structure is basically the same: a bright central core centered on an apparently elliptical structure having a major axis at a position angle of 80◦. While the position angle of the elliptical structures agree well at ∼ 2.2µm and 8.4 Ghz, there is a suggestion that the position angle of the ellipse at 12.5 µm is 15◦ smaller. At 8.4 Ghz the central source contains 80% of the radio flux(C91), while at ∼ 2.2µm, the unresolved peak contains 35% of the flux in a 5′′ diameter beam (S99).

As shown in Figure 2, at 12.5 µm approximately half of the total flux measured in a 4′′ diameter beam is contained in the central compact source. The inner central compact source is indistinguishable from the PSF (Figure 2). From the dispersion in the sizes determined for the PSF star (BS3775) the compact source has a size (3σ) less that 0.22′′. The size measured at 17.9µm is consistent with that measured at 12.5 µm.

The data of Figure 3 and Table 3 show that at 12.5 µm the flux measured by MIRLIN is 60% of that measured by IRAS. This suggests that there is low surface brightness ∼ emission extended on a larger scale than 4′′. The 12.5µm flux density contained in the compact source is 30% of the total flux density at 12µm measured by IRAS. ∼

4.4. Markarian 231 = UGC 8058 = IRAS 12541+5708

Markarian 231 has long been known both as an infrared luminous system (Rieke and Low, 1972) and a Seyfert 1 system (Boksenberg et al. 1977). The IRAS survey found that

12 it is the most luminous object within 300 Mpc, with a bolometric luminosity of 3 10 L⊙. × MIRLIN images of Markarian 231 were obtained at seven wavelengths from 7.9 to 17.9µm. In addition, LWS imaging was obtained at 12.5µm. Figure 1d shows contour plots – 14 – of Markarian 231 at 12.5 µm from the LWS data and data at 8.4 Ghz(C91). The emission at all the mid-infrared wavelengths is dominated by the very bright central source, and is effectively unresolved. The radio continuum is entirely from the bright nucleus.

LWS imaging was used to establish the size of the central source at 12.5µm. The LWS observations consisted of eleven measurements of the size of Markarian 231, interleaved with nine measurements of a near-by PSF star(BS 4905) over a period of two hours. The average FWHM of the PSF star and Markarian 231 were indistinguishable and were slightly greater than the diffraction limit of the telescope and the rms (population) scatter in the FWHM was 0.04′′. The central source thus is unresolved in these observations. We have chosen to set a 3σ limit on its FWHM of 0.13′′ on the basis that a difference three times the average standard deviation of the FWHM would have been easily detected.

The SED of Markarian 231 is plotted in Figure 3, and shows clearly that all the flux measured by IRAS is detected in the 4′′ diameter beam. The SED shows a dip at 10µm apparently due to foreground silicate absorption of the central source.

4.5. Markarian 273 = UGC 8696 = IRAS 13428+5608

The MIRLIN observations of Markarian 273 consist of imaging at seven wavelengths from 7.9µm to 17.9µm. Figure 1e shows a comparision of the MIRLIN 12.5µm data with the 8.4 GHz map (C91) and the NICMOS 2.2µm image (S99), while a montage of images of Markarian 273 is shown in Figure 4. At 12.5µm the source structure is basically double, with the northern source brighter than the southwestern source at most wavelengths. A similar morphology is seen at 2.2µm (S99), but not in the radio (C91). In comparing the 12.5µm, NICMOS and 8.4 Ghz data we have assumed the spatial coincidence between the bright northern source in each image. The overall ageement between the 12.5µm and NICMOS – 15 – images is excellent, while the discrepancy between the positions of the southwestern source in the infared and radio is real. The position angle between the northern and southwestern sources is the same (to within the uncertainties) in the 12.5µm, 2.2µm and 8.4 GHz images, while the source separations are 1.0′′, 1.0′′ and 1.2′′ at these same wavelengths. The apparent separation between the infrared and radio sources, 0.2′′, corresponds to 160pc at the distance of Markarian 273. Note also that the flux ratio between the northern and southwestern source is considerably different between the infrared and radio.

As shown in Figure 4, the structure seen in Markarian 273 at all infrared wavelengths from 3.4µm to 17.9µm is basically the same. There is no evidence for a mid-infrared counterpart to the southeastern source seen at 8.4 Ghz by C91 (Figure 1e) and at 5 Ghz by Knapen et al. (1997). The bright source present to the southeast in the radio map is suggested to be a background radio object by Knapen et al., although the statistical likelihood of a background source falling within 1′′ of the location of the radio core of Markarian 273 is < 10−5. There is a faint compact object found in the NICMOS images (S99) at the location of southeast radio source, although there is no evidence for any mid-infrared emission at this location. The upper limit on the 12.5µm flux density at the location of the southeastern radio source of 15 mJy is 1/10 that of the northern source. This limit, coupled with the 8.4GHz flux density of this source of 3.2 mJy (C91), leads to fν (12.5µm) < a limit on the ratio of flux densities fν (8.4GHz) 4.7, consistent with the observed ratio fν (12.5µm) fν (8.4GHz) =4.5 for the northern source. By comparison, if the radio ( 0.2mJy at 8.4 GHz, C91, Knapen et al.) and infrared (58 mJy at 12.5µm) flux densities of the southwestern

fν (12.5µm) source are assumed to be associated with the same physical source, the ratio fν (8.4GHz) = 290, a factor of 65 times greater than that seen in the northern source, making this nucleus far more infrared bright relative to the radio than the northern source. If the southwestern radio and infrared sources are not related, the southwestern infrared source has an even

fν (12.5µm) larger fν (8.4GHz) ratio. – 16 –

As can be seen in Figures 1e and 4, the relative brightnesses of the nuclei change substantially as a function of wavelength. In order to get an estimate of the relative flux emitted by the two nuclei, the flux in two identically sized square beams was measured at the same locations on the images at all wavelengths presented; the beams, 1.6′′ on a side, are shown in Figure 4. The sky was taken northwest and southeast of each of the squares. The efficiency of the square beams was estimated by placing the PSF of the appropriate wavelength at the locations of the nuclei within the square beam. The photometry of the NICMOS image was arbitrarily normalized so that the sum of the fluxes from the two squares equals the sum of the fluxes give by Knapen et al. (1997) and the flux density in a 4′′ diameter beam equals that of Knapen et al. The results of these photometric measurments are shown in Figure 5.

At 3.4µm and 10.3µm the southwestern source is slightly (10%) brighter than the northern source. At the other mid-infrared wavelengths the northern source is between a factor of 1.5 and 3 times brighter than the southwestern source. The different SED’s of the components implied by their changing relative brightnesses is shown directly in Figure 5, where photometry of the two nuclei is shown. The northern source clearly shows a signficantly deeper absorption at 10µm than the southwestern source. This is consistent with the distribution of CO emission which strongly peaks on the northern source (Downes and Solomon, 1998), but is inconsistent with the near infrared colors found by S99, where the northern source has a smaller J-K color than does the southwest source.

The SED of Figure 5, particularly the relatively strong emission at 3.5 and 11.7 µm, suggests aromatic hydrocarbon emission (PAH) emission associated with the southwest nucleus. Such emission is seen in the spectra of many ULIRGs (Genzel et al. 1998). The emission from the southwest source could account for all of the 11.3 µm emission seen in the spectrum of the combined nuclei by Dudley(1999).The northern source does not show a – 17 – similar peak at 11.7µm, but clearly shows the effects of substantial overlaying cold silicate absorption.

The MIRLIN photometry shown in Figure 3 was obtained using a 4′′ diameter beam centered between the two nuclei. The comparison of the MIRLIN and IRAS data (Figure 3) demonstrates that the mid-infrared emission in this source emerges virtually entirely from the central region. Because of the proximity of the two nuclei and the fairly low signal-to-noise ratio of the nuclei in some of the images, it is difficult to cleanly determine the sizes of the two sources. The data suggest that both nuclei have unresolved mid-infrared cores with sizes <0.4′′. In the northern nucleus approximately 2/3 of the 12.5µm flux is unresolved, while in the southwestern nucleus, about half the flux is unresolved.

4.6. Arp 220 = UGC 9913 = IRAS 15328+2340

The data for this source were discussed in detail in Soifer et al. (1999) and are included in Figures 1f and 3 for completeness.

4.7. IRAS 17208-0014

MIRLIN imaging was obtained at 7.9 and 12.5µm. The 12.5µm image is shown in Figure 1g. The 2.2µm NICMOS image of IRAS 17208-0014 (Figure 1g, S99) shows it to be an extended source of roughly 2′′ diameter with a number of luminous starburst centers to the northeast of the nucleus. The mid-infrared image is apparently also similar, but the limitation due to the chop-nod spacing and the relatively low infrared surface brightness results in an image with low signal-to-noise ratio. The size of the 12.5 µm source of 2′′ ∼ or 1.6 kpc FWHM (Figure 2) is quite similar as well to that found in CO by Downs and Solomon (1998) and Jogee et al.(2000). – 18 –

The best illustration of the spatial distribution of the mid-infrared emission is shown in Figure 2. The source flux smoothly increases up to 5′′ diameter, with a FWHM of 2′′; the apparent leveling off of the enclosed flux at a diamteter of 5′′ is caused by the artificial imposition of the sky level between diameters of 5 and 6′′.

Evidence that the source is not extended much beyond a diameter of 5′′ is provided by Figure 3 which shows that, within the uncertainties, the enclosed flux within a 4′′ diameter beam is the same as within the IRAS beam. The lack of a dominant compact component in the mid-infrared make this object very different from the other ULIRGs in this sample.

The radio emission from IRAS 17208-0014 is much more compact than the infrared emission. At least 90% of the 1.4GHz flux comes from a region smaller than 2′′ 2′′ (FWHM, × Condon et al. 1996). Martin et al. (1989) measure a radio source with a size of 0.26′′ 0.32′′ × at 1.6 GHz that accounts for 90% of the 1.4 GHz flux density, suggesting a compact radio ∼ core. This core is clearly not dominant in the mid-infrared imaging. The near infrared NICMOS imaging of this galaxy (S99) suggests that the nucleus could be obscured at 2.2µm, but the lack of detection of a strong compact core at 12.5µm is quite peculiar, since the radio source is so compact.

5. Discussion

5.1. Compact Sources in ULIRGs

The data presented here, with angular resolution of <0.3′′, represent the highest spatial resolution achieved to date in probing the central regions of ULIRGs in the thermal infrared. Table 4 summarizes the key observations, the sizes of the compact sources found in the ULIRGs at 12.5 µm, 17.9 µm and 24.5 µm, as well as the fraction of the total flux at these wavelengths emitted from these cores. All of the ULIRGs except IRAS 17208 0014 have − – 19 – a nuclear source in the mid-infrared of size 200 pc that produces a substantial fraction, ≤ if not all, of the galaxy’s mid-infrared luminosity. In three sources, IRAS 05189 2524, − IRAS 08572+3915 and Markarian 231, a single unresolved nucleus produces between 80% and 100% of the total mid-infrared luminosity. In Markarian 273 and Arp 220, two nuclei produce nearly all the mid-infrared emission in the galaxies. In the case of Markarian 273 the nuclei are substantially unresolved, while in Arp 220 at least one nucleus is resolved. Because of its low redshift the linear size of the resolved nucleus in Arp 220 is smaller than the limit on all the unresolved nuclei except Markarian 231(see Table 5). In only one compact system, UGC 5101, does the compact nucleus produce significantly less than half the total mid-infrared emission from the system.

As noted previously, IRAS 17208 0014 presents an interesting puzzle. There is no − mid-infrared evidence for a compact nuclear source, but the radio data (C91, Martin et al. 1989) strongly suggest a power source with a size <0.3′′. In the following discussion, we will adopt the morphology implied by the 12.5µm image, but we note that extremely high 12.5µm extinction could be leading to an erroneous interpretation of the nuclear structure in IRAS 17208 0014. − Table 5 shows that the compact nuclear sources have linear sizes typically <200 pc in diameter. In Arp 220 the diameter of the western source is 140 pc at 24.5µm (Soifer et al. 1999), while in Markarian 231, the diameter of the 12.5µm source is <100 pc.

5.2. AGNs as the Luminosity Source in ULIRGs

One of the goals of the observations reported here is to distinguish the ultimate source of the infrared luminosity in these systems by means of the spatial extent of the infrared emission. If the dust is heated by a bright compact object such as an AGN, the dust – 20 – emission will be compact, and the size of the infrared source will increase with increasing wavelength roughly as r λ2.5, since the dust temperature will decrease with radius as ∼ T r−0.4 and the wavelength of maximum emission λ T −1. If, on the other hand, ∼ max ∼ the heating is due to many hot stars in a starburst region, then the dust emission will be distributed on the same scale as the young, luminous stars. Since the dust heating is due to the proximity of the dust to the luminous young stars there will be little or no temperature gradient seen. Spatially extended 10 µm emission therefore implies starburst activity, while spatially compact 10 µm emission could be due to either a compact starburst or an AGN.

To quantify the compact mid-infrared emission we utilize the brightness temperatures derived from our observations and the color temperatures determined from the IRAS 25

and 60µm observations of these objects. We assume that Tc measures the temperature of the grains heated by the luminosity source; T should be T if the optical depth τ < 1. b ≤ c ν Only in the limit of large optical depth should the brightness temperature approach the color temperature.

Based on the MIRLIN observations of Arp 220, Soifer et al. (1999) showed that the dominant luminosity source, Arp 220W, is optically thick with a brightness temperature at 25 µm equal to the color temperature between 25 and 60 µm. The observations reported here permit analogous constraints to be placed on some of the nuclear sources in other ULIRGs in this sample. In Markarian 231, IRAS 05189 2524 and IRAS 08572+3915, the − limits on the sizes of the central infrared sources are sufficiently small that it is worthwhile considering how these sizes constrain the far infrared sources.

The size limit on Markarian 231 of 0.13′′ at 12.5µm places the most severe limits on models of the infrared emitting region. In this case the observed flux density and size limit lead to a lower limit on the brightness temperature of T 141K at 12.5µm. If the b ≥ infrared emitting source has the same size at 17.9 µm, 25 µm and 60 µm, the corresponding – 21 – brightness temperatures are 124 K, 128 K and 241 K. The brightness temperatures at 12.5µm, 17.9µm and 25 µm are consistent with each other and with a nearly optically thick source over this wavelength range, but the 60 µm brightness temperature is much too large to be consistent with dust models with reasonable grain properties. This implies that the 60µm source size is much larger than the 12.5µm size.

The color temperature measured between 25 and 60 µm, determined from the IRAS measurements, is 85 K. The size of a blackbody of 85 K required to produce the observed 60µm flux density is 0.4′′, significantly larger than the limit on the source size at 12.5µm.

The source size determined if Tc = Tb represents a lower limit on the size of the 60µm source. For example, if the size of the 60µm source were equal to the 1′′ size found in ∼ CO by Downs and Solomon (1998), the 60µm brightness temperature would be 53K, much less than that inferred by assuming all the 25 and 60 µm emission is emerging from the same region. The size limits are based on the assumption of large optical depth (τ >>1); smaller optical depth would imply an even larger source. Since significant emission is seen

at λ < 60µm, either τ60µm <1 or we are viewing a tilted optically thick disk, as suggested for Arp 220 (Scoville et al. 1998, Downs and Solomon, 1998, Sakamoto et al. 1999, Soifer et al. 1999). In either case, the size constraints show that in Markarian 231 the infrared source increases in size with increasing wavelength. This is consistent with the picture of

12 an AGN with luminosity of 3 10 L⊙ heating the surrounding dust and producing a radial × temperature gradient.

In this sample there are two additional ULIRGs where significantly more than 50% of the mid-infrared luminosity emerges from a compact central source, IRAS 05189-2524 and IRAS 08572+3915. In IRAS 05189-2524, if the size of the 60µm emitting source is the same

as the limit on the 12.5µm source measured here, Tb(60µm)=101 K, which is significantly

greater than Tc =85 K from the IRAS 25 and 60µm data. For the same assumptions, in – 22 –

IRAS 08572+3915 Tb(60µm)=77 K is slightly less than Tc =83 K determined from the IRAS data at 25 and 60µm. For both sources, using dust properties described by Draine and Lee (1984) the sizes calculated for silicate and graphite dust grains heated by central sources of the observed bolometric luminosity, and having dust temperatures equal to the observed color temperatures are larger by a factor of 1.5 to 3 than the observed sizes of these sources at 12.5µm. This suggests that the sizes of these infrared sources also increase with increasing mid and far-infrared wavelength.

High apparent surface brightnesses and a radial temperature gradient are natural consequences of centrally heated dust emission surrounding a high luminosity source, i.e., an AGN. In Markarian 231 there is direct evidence that the size of the infrared source increases with increasing wavelength in the mid and far infrared. In IRAS 05189-2524 and IRAS 08572+3915, there is strong evidence for such radial temperature gradients based on the compact sources observed in the mid-infrared. The observed luminosities are adequate to quantitatively account for the inferred variation of dust temperature with radius if the dust is heated by a central source in these objects.

5.3. Starbursts as the Luminosity Source in ULIRGs and Energy Production Rates

The compact nature of the sources in the objects observed leads to very high luminosities produced in small volumes. By comparing these energy production rates with those found in local examples of starburst environments, we can determine whether known starburst environments can generate the apparent energy production in the ULIRGs.

5 In the core of Orion the luminosity of 2 10 L⊙ emerges from a 0.3 pc size region × (Werner et al. 1976), while in the nuclear starburst in M82, the nearest starburst galaxy, the – 23 –

10 ′′ ′′ luminosity of 3 10 L⊙ is emitted over a projected area of 450 75pc (30 5 ) (Rieke et × × × al. 1980). These observations lead to apparent surface brightnesses of emergent luminosity

12 2 of 1 3 10 L⊙/kpc . If, as seems likely, M82 is being viewed edge-on, then the mean − × 11 2 surface brightness of this system, if viewed face on, is closer to 2 10 L⊙/kpc . × A survey of the maximum surface brightness in nearby starburst galaxies in the UV, infrared and radio by Meurer et al. (1997) leads to peak global averaged surface brightnesses

11 2 in starburst galaxies of 2 10 L⊙/kpc . Meurer et al. find that the peak surface ∼ × 13 2 brightnesses of clusters in starburst galaxies reaches 5 10 L⊙/kpc over regions 10 ∼ × ∼ pc in diameter. Meurer et al. suggest that the global peak surface brightnesses is due to a mechanism whereby the pressure from stellar winds and supernovae is able to regulate the star formation rates in galaxy disks, while the formation of massive stars in giant HII regions is regulated by the “Kennicutt Law”, i.e. star formation occurs when the surface gas density exceeds a critical value that depends on the rotational velocity and shear (Kennicutt, 1989).

The surface brightnesses observed in the ULIRGs can be estimated if we assume that the size (or limit) of the region measured in the mid-infrared is representative of the infrared emitting region as a whole. With this assumption, we have presented the surface brightnesses for the observed ULIRGs in Table 5. We include two values of the surface brightness, a “maximum” which assumes that the fraction of the total infrared luminosity emerging from the compact core is the same as the fraction of the mid-infrared flux density emerging from the core, and a “minimum” that assumes that only the observed mid-infrared luminosity emerges from the compact core. In the case of Markarian 231, the former number is unphysically large (section 5.2). In the other systems, the two values represent likely upper and very conservative lower limits on the apparent surface brightness of the infrared emission. – 24 –

Only in the case of IRAS 17208 0014, where the near and mid-infrared observations − suggest the luminosity is due to star formation, is the measured infrared surface brightness

11 2 comparable to the peak global value 2 10 L⊙/kpc found by Meurer et al.(1997). Other × than this case, the inferred surface brightnesses are much greater than found globally in disk starburst systems. The highest apparent surface brightnesses in the ULIRGs are found in the western nucleus of Arp 220, which is most likely powered by a nuclear starburst (Smith et al. 1998b) and the cores of Markarian 231 and IRAS 05189 2524, which are probably − powered by AGN. In these cores the apparent surface brightness approaches or surpasses the surface brightness of clusters in starburst galaxies, while in the other systems (excluding

13 2 IRAS 17208 0014) the surface brightness most probably exceeds 10 L⊙/kpc . In all these − cases the surface brightness over regions of hundreds of parsecs across are comparable to the brightness found in 10 pc sized clusters within starburst systems by Meurer et al. ∼ For the ULIRGs that are powered by star formation, the large surface brightnesses and sizes suggest that the global starburst intensity regulating mechanism of Meurer et al. (1997) does not apply. Gravitational instability allows for higher gas densities and hence higher star-formation rates in the solid-body (no shear) rotation nuclear regions. Since the nuclei of most ULIRGs are mergers, far from equilibrium, the gas densities may be very high, leading to high star formation rates and correspondingly high surface brightnesses for short times. ULIRGs are clearly extraordinary events in the evolution of galaxies. Several studies (e.g., Sanders et al. 1988, Scoville, Yun and Bryant 1997, Downes and Solomon 1998) suggest that the nuclear disks are unstable in some ULIRGs. The large gas masses in these environments might then lead to nuclear star formation rates much greater than those found in less extreme systems. – 25 –

5.4. A Panchromatic View of ULIRGs

Table 6 summarizes the characterization of the ULIRGs applying many different diagnositics of the underlying energy sources. Here we have summarized the structure in the mid-ir from the data presented here, in the near infrared from S99 and Surace and Sanders (1999), and in the radio from the VLBI data of Smith et al. (1998). The optical spectral classifications are from Sanders et al. (1988) and Veilleux et al (1995). The near-infrared spectral classifications are from Veilleux et al. (1999) and Murphy (private communication). The ISO spectral classifications are based on the strength of the PAH features and are taken from Genzel et al. (1998) and Lutz, Veilleux and Genzel (1999). The X-ray classifications are from ASCA results and are taken from Iwasawa (1999) and Nakagawa et al. (1999).

In none of the objects in this sample, excepting perhaps Markarian 273 and IRAS17208- 0014, is there a completely consistent set of “classifications” based on the wide range of wavelengths over which these objects are observed. In the object that is the clearest case for a central AGN luminosity source, Markarian 231, the hard X-ray emission is not consistent with the simplest view of such power sources (Nakagawa et al. 1999). In the predominantly Seyfert nucleus IRAS 05189-2524 there is not a dominant compact radio core. In Markarian 273 the two nuclei clearly seen in high resolution imaging complicates the classification, since none of the spectral classifications resolved the separate nuclei. In UGC 5101, another source where the optical, infrared and X-ray indications point towards a starburst power source, a compact radio core characteristic of a central AGN exists. Compact mid-infrared sources occur in six of seven of these sources, where other evidence points strongly to both AGN and starbursts as the dominant underlying power source. Table 6 shows that most ULIRGs do not present a picture over all wavelengths consistent exclusively with either an AGN or starburst power source. In most objects observed here, it is likely that both contribute to some extent to their total bolometric luminosity, though Markarian 231 and – 26 –

IRAS 17208-0014 appear to represent the extremes of dominant AGN and starburst sources.

6. Summary and Conclusions

Diffraction limited mid-infrared observations of a sample of seven ULIRGs obtained on the Keck 10m telescopes have shown:

1 - Extremely compact structures, corresponding to emission on spatial scales of <100 200 pc, are seen in six of the seven systems observed. The mid-infrared emission in − these galaxies generally emerges from a single region of diameter < 1 kpc, and in most cases from regions <200 pc in diameter.

2 - The upper limit on the diameter of the 12.5 µm source in Markarian 231, 0.13′′, is too small to be consistent with the observed 60 µm flux density emerging from a region this size, and demonstrates that in this AGN the source size increases with increasing wavelength (decreasing dust temperature).

3 - In IRAS 05189-2524 and IRAS 08572+3915 there is strong circumstantial evidence that the size of the infrared source increases with wavelength between 12.5 and 60 µm, suggesting that these too are predominantly AGN powered ULIRGs.

4 - In the other objects observed, the angular size limits for the compact mid-infrared sources are not yet adequate to distinguish between AGN and starburst luminosity sources, although the small sizes are consistent with a central luminosity source heating the surrounding dust as in AGN models in all the sources except IRAS 17208-0014.

5 - If star formation powers the luminosity in the very compact infrared sources seen in these systems, the star formation rates averaged over a few hundred parsecs is comparable to that seen in the brightest 10 pc sized clusters in nearby starbursts, and exceeds by ∼ – 27 – factors of up to 100, the global star formation rates in nearby starburst galaxies.

7. Acknowledgments

We thank J. Aycock and R. Campbell for assistance with the observations, R. Goodrich, R. Moskitis and the entire Keck staff for their help establishing the visitor port and to Barbara Jones, Rick Puetter and the Keck team that brought the LWS into service, enabling these observations. The W. M. Keck Observatory is operated as a scientific partnership between the California Institute of Technology, the University of California and the National Aeronautics and Space Administration. It was made possible by the generous financial support of the W. M. Keck Foundation. This research has made use of the NASA/IPAC Extragalactic Database which is operated by the Jet Propulsion Laboratory, Caltech under contract with NASA.

B.T.S, G.N., K.M. and E.E. are supported by grants from the NSF and NASA. J.A.S. and B.T.S. are supported by the SIRTF Science Center at Caltech. SIRTF is carried out at the Jet Propulsion Laboratory. E.E.B., A.J.W, N.Z.S. and A.S.E. were supported by NASA grant NAG5-3042. This work was carried out in part (M.W.W., M.R.) at the Jet Propulsion Laboratory, operated by the California Institute of Technology, under an agreement with NASA. The development of MIRLIN was suported by NASA’s Office of Space Science. The National Radio Astronomy Observatory (J.J.C.) is a facility of the National Science Foundation operated under cooperative agreement by Associated Universities, Inc. – 28 –

Fig. 1.—

Figure 1a. Images of IRAS 05189-2524 obtained at 2.2µm (S99, upper left), 8.4 GHz (C91, lower left), and 12.5µm (this work, upper right). The lower right panel is the PSF star image obtained at 12.5µm contemporaneously with those of IRAS 05189-2524 and used to determine the 12.5µm size of IRAS 05189-2524. The peak of each image is at the location 0,0 in each panel. The contours in each panel are set so that the peak contour is 90% of the maximum surface brightness. The contours are spaced down by a factor of 1.34 so that the third contour down from the peak is at 50% of the peak surface brightness. The ninth coutour down is 8.5% of the peak. Below this level the contours are spaced by a factor of two. The faintest contour is set at about 3 times the noise in the image. The mid-infrared data are from the MIRLIN observations. The hatched circle show the FWHM of the PSF for the 2.2µm observations.

Figure 1b. Images of IRAS 08572+3915. The format, source locations and contour levels are the same as in Figure 1a. The mid-infrared data are from the LWS observations. The 2.2µm and 8.4 GHz data are from S99 and C91 respectively. The hatched circles show the FWHM of the PSF for the appropriate observations.

Figure 1c. Images of UGC 5101. The format is 2.2µm image(S99) - left, 12.5µmimage - center, 8.4 GHz image(C91) - right. The source locations and contours are set as in Figure 1a. The mid-infrared data are from the MIRLIN observations. The 2.2µm and 8.4 GHz data are from S99 and C91 respectively. The hatched circles show the FWHM of the PSF for the observations.

Figure 1d. Images of Markarian 231. The format is 12.5µm image - left, 12.5µm PSF star - center, 8.4 GHz image - right. The source locations and contours are set as in Figure 1a.The mid-infrared data are from the LWS observations. The 8.4 GHz data are from C91. The hatched circle in the 8.4 GHz panel shows the FWHM of the PSF for this observation. – 29 –

Figure 1e. Images of Markarian 273. The format, source locations and contour levels are the same as in Figure 1c. The mid-infrared data are from the MIRLIN observations. The 2.2µm and 8.4 GHz data are from S99 and C91 respectively. The + symbols are at the same coordinates in each panel, and are centered on the peaks of the 2.2µm image. There is excellent agreement between the locations of the 2.2µm and 12.5µm peaks. The discrepancy between the positions of the infrared peaks and the radio peak is real (see text). The hatched circles show the FWHM of the PSF for the observations.

Figure 1f. Images of Arp 220. The format, source locations and contour levels are the same as in Figure 1c. The mid-infrared data are from the MIRLIN observations and are from Soifer et al. (1999). The 2.2µm and 8.4 GHz data are from S99 and C91 respectively. The hatched circle in the 12.5µm panel shows the FWHM of the PSF for the observation.

Figure 1g. Images of IRAS 17208-0014. The left image is 2.2µm, the right image is 12.5µm. The source locations and contour levels are the same as in Figure 1a. The mid-infrared data are from the MIRLIN observations. The 2.2µm data are from S99. The hatched circles show the FWHM of the PSF for the observations.

Fig. 2.—

The curves of growth for the objects (solid) and PSF’s (dashed) for the five ULIRGS which are single nuclei. In all five the PSF’s have been normalized so the enclosed flux of the PSF within a radius of 2.5′′ equals that of the object. For the objects, the sky level was set to be the value in an annulus between 2.5 and 3′′. In the plots for IRAS 08572+3915 and UGC 5101 two normalizations are used for the PSF; in the lower curve the PSF is set to match the central profile of the galaxy, while the upper curve matches the enclosed flux from the PSF and the galaxy at the outer radius of the measurement beams, as for the other galaxies. The data of IRAS0857+39 and Markarian 231 were obtained with LWS as outlined in the text; the other three data sets are from MIRLIN observations. – 30 –

Fig. 3.—

The spectral energy distributions of the seven ULIRGS studied are shown as a function of the rest wavelength from 1 µm to 1 mm. Filled circles represent the MIRLIN data, ∼ while open circles have been taken from Carico et al. (1988) (near infrared data, 5′′ diameter beam) and Soifer et al. (1989) (IRAS data 1′ diameter beam). The MIRLIN photometry ∼ has been calculated for a 4′′ diameter beam centered on the object or between the two nuclei if multiple nuclei are present. As discussed in the text, the photometry is based on IRAS standards.

Fig. 4.—

Montage of images of Markarian 273. The images at wavelengths from 7.9µm through 17.9µm obtained with MIRLIN. The image at 3.4µm was obtained on the 200-inch Hale Telescope using the Cassegrain infrared camera. The images are centered so that the northern peak is at the location 0,0 in each panel. The locations of the boxes in the 12.5µm image show where the photometry of the peaks was obtained in all the images.

Fig. 5.—

Spectral energy distributions of the separate peaks in Markarian 273 as determined from the images presented in Figure 4. The total flux (also presented in Figure 3b) was measured in a 4′′ diameter beam centered between the peaks and was distributed between the two sources by doing photometry in non-overlapping boxes (displayed in the 12.5µm panel of Figure 4) centered on the two peaks. The uncertainties shown are statistical. The systematic uncertainty depends on the wavelength, but is probably in the range 15 - 20 %. – 31 –

REFERENCES

Barger, A.J., Cowie, L.L., Sanders, D.B., Fulton, E., Taniguchi, Y., Sato, Y., Kawara, K. and Okuda, H. 1998, Nature, 394, 248

Blain, A.W., Smail, I., Ivision, R.J. and Kneib, J.-P. 1999, MNRAS, 302, 632

Boksenberg, A., Carswell, R.F., Allen, D.A., Fosbury, R.A.E., Penston, M.V. and Sargent, W.L.W. 1977, MNRAS, 178, 451

Carico, D.P., Sanders, D.B., Soifer, B.T., Elias, J.H., Matthews, K. and Neugebauer, G. 1988, AJ, 95, 356

Carico, D.P., Graham, J.R., Matthews, K., Wilson, T.D., Soifer, B.T., Neugebauer, G., and Sanders, D.B. 1990, ApJ, 349, L39

Condon, J. J., Huang, Z. P., Yin, Q. F., Thuan, T. X., 1991, ApJ, 378, 65(C91)

Condon, J. J., Helou, G., Sanders, D.B., and Soifer, B.T. 1996, ApJS, 103, 81

Downes, D. and Solomon, P. M., 1998, ApJ, 507, 615

Draine, B.T. and Lee, H.M., 1984, ApJ,285 , 89

Dudley, C.C., 1999, MNRAS, in press, astro-ph9903250

Elias,J.H., Frogel, J.A., Mattthews, K. et al. 1982, AJ, 87, 1029

Genzel, R., Lutz, D., Sturm, E., Egami, E., Kunze, D., Moorwood, A. F. M., Rigopoulou, D., Spoon, H. W. W., Sternberg, A., Tacconi-Garman, L. E., Tacconi, L. and Thatte, N., 1998, ApJ, 498, 579

Heckman, T. 1998, to appear in “Ultraluminous Galaxies: Monsters or Babies? ”, astro-ph9903041 – 32 –

Iwasawa, K. 1999, MNRAS, 302, 96

Joint IRAS Science Team, 1989, IRAS Point Source Catalog, Version 2 (Washington, DC, US Government Printing Office)

Jogee, S., Bell, G.R. and Scoville, N.Z. 2000, in preparation

Joseph, R.D. 1999, to appear in “Ultraluminous Galaxies: Monsters or Babies? ”, preprint

Keto, E., Ball, R., Arens, J., Jernigan, G., and Meixner, M., 1992, ApJ, 389, 223

Kennicutt, R.C. 1989, ApJ, 344, 685

Kim, D.-C., Veilleux, S. and Sanders, D.B. 1998, ApJ, 508, 627

Knapen, J.H., Laine, S., Yates, J.A., Robinson, A., Richards, A.M.S., Doyon, R. and Nadeau, D. 1997, ApJ, 490, L29

Lilly, S.J., Eales, S.A., Gear, W.K., Webb, T.M., Bond, J.R., & Dunne, L., 1999, to appear in STScI workshop “Formation of Galactic Bulges”, astro-ph9903157

Lutz, D., Veillleux, S. and Genzel, R., 1999, ApJ, 517, L13

Martin, J.M., Bottinelli, L., Dennefeld, M., Gouguenheim, L. and Le Squeren, A.M. 1989, A&A, 208, 39

Meurer, G.R., Heckman, T.M., Lehnert, M.D., Leitherer, C. and Lowenthal, J. 1997, AJ, 114, 54

Miles, J. W., Houck, J. R., Hayward, T. L., Ashby, M. L. N., 1996, ApJ, 465, 191

Nakagawa, T., Kii, T.,Fujimoto, R., Miyazaki, T., Inoue, H., Ogasaka, Y., Arnaud, K., and Kawabe, R. 1999, to appear in proceedings of IAU Symposium 186, eds. Barnes, J.E. and Sanders, D.B. – 33 –

Ressler, M. E.,Werner, M. W., Van Cleve, J. and Choa, H., 1994, Experimental Astronomy, 3, 277

Rieke, G.H. 1988, ApJ, 331, L5

Rieke, G.H. and Low, F.J. 1972, ApJ,

Rieke, G.H., Lebofsky, M.J., Thompson, R.I., Low, F.J. and Tokanaga, A.T., 1980 ApJ, 238, 24

Sakamoto, K., Scoville, N. Z., Yun, M. S., Crosas, M., Genzel, R., and Tacconi, L. J., 1999, ApJ, 514, 68

Sanders, D. B. 1999, to appear in “Ultraluminous Galaxies: Monsters or Babies? ”, astro-ph9908297

Sanders, D. B., Soifer, B. T., Elias, J., Madore, B., Matthews, K., Neugebauer, G., and Scoville, N., 1988, ApJ, 325, 74

Sanders, D.B., Egami, E., Lipari, S., Mirabel, I.F., & Soifer, B.T., 1995, AJ, 110, 1993

Scoville, N. Z., Evans, A. S., Dinshaw, N., Thompson, R., Rieke, M., Schneider, G., Low, F. J., Hines, D., Stobie, B., Becklin, E.E., and Epps, H., 1998, ApJ, 492, L107

Scoville, N. Z., Evans, A. S., Thompson, R., Rieke, M., Hines, D.C., Low, F.J., Dinshaw, N., Surace, J.A. and Armus, L.. S99, 1999, AJ, submitted

Scoville, N. Z., Yun, M. S., and Bryant, P. M., 1997, ApJ,484,702

Smith, H. E., Lonsdale, C. J. and Lonsdale, C. J., 1998, ApJ, 492, 137

Smith, H. E., Lonsdale, C. J. Lonsdale, C. J. and Diamond, 1998b, ApJ, 493, L17

Soifer, B. T., Boehmer, L., Neugebauer, G. and Sanders, D.B. 1989, AJ, 98, 766 – 34 –

Soifer, B.T., Neugebauer, G., Matthews, K., Becklin, E.E., Ressler, M., Werner, M.W., Weinberger, A.J., and Egami, E. 1999, ApJ, 513, 207

Soifer, B. T., Sanders, D. B., Madore, B. F., Neugebauer, G., Danielson, G. E., Elias, J. H., Lonsdale, C. J., and Rice, W. L., 1987, ApJ, 320, 238

Surace, J.R. & Sanders, D.B. 1999, ApJ, 512, 162

Surace, J.R., Sanders, D.B., Vaca, W.D. & Veilleux, S., 1998, ApJ, 492, 116

Veilleux, S., Kim, D.-C., Sanders, D.B., Mazzarella, J.M. and Soifer, B.T. 1995, ApJS, 98, 171

Veilleux, S., Sanders, D.B., and Kim, D.-C. 1999, ApJ, 522, 139

Werner, M.W., Gatley, I., Becklin, E.E., Harper, D.A., Loewenstein, R.F., Telesco, C.M. and Thronson, H.A. 1976, ApJ, 204, 402

This manuscript was prepared with the AAS LATEX macros v5.0. – 35 –

Table 1

Table 1: Basic Properties of ULIRGS Observed

Name z log L Spectrum linear scale

′′ Lbol[L⊙] Kpc/ IRAS05189-2524 0.043 12.08 Sey1 830 IRAS08572+3915 0.058 12.10 LINER 1200 UGC5101 0.040 12.00 LINER 820 Markarian 231 0.040 12.52 Sey 1 800 Markarian 273 0.037 12.13 Sey 800 Arp220 0.018 12.17 LINER 360 IRAS17208-0014 0.040 12.46 HII 800 – 36 –

Table 2

Table 2: Log of Keck II mid-infrared Observations of ULIRGS

Object 7.9µm 8.8µm 9.7µm 10.3µm 11.7µm 12.5µm 17.9µm 24.5µm IRAS05189 √ √ IRAS08572 √√√ √∗ √ UGC 5101 √∗ √ Mrk 231 √√√√ √ √∗ √ Mrk 273 √√√√ √ √ √ Arp 220 √√√√ √ √ √ √ IRAS17208 √ √

∗size determination made with the Long Wavelength Spectrograph in March 99 – 37 –

Table 3

Table 3: MIRLIN and IRAS Flux Densities for ULIRGS

Object 7.9µm 8.8µm 9.7µm 10.3µm 11.7µm 12.5µm 12.5µm∗ 17.9µm 24.5µm 25µm∗ 60µm∗ [mJy] [Jy] IRAS05189 720 740 3.2±0.2 3.5 13.9 IRAS08572 615±18 434±10 221±10 218±10 340 552±46 1.8 7.5 UGC 5101 185±10 260 256±42 1.1 12.1 Mrk 231 1370±43 1488±18 1081±30 1075±21 1355±15 1836±20 1930 3249±49 8.8 33.6 Mrk 273 267±15 141±10 92±10 38±6 150±5 210 ±10 230 442±31 2.3 23.7 Arp 220 521±15 155±6 124±14 59±10 175±8 404 ±12 500 1170±50 9.8±0.2 8.0 104.0 IRAS17208 315±25 202±13 200 1.7 34.1

The MIRLIN photometry refers to a 4′′ diameter beam. The uncertainties listed in the table are statistical only. Photometric uncertainties in the MIRLIN data are 5% for λ< 20µm, and 10% for λ> 20µm. The ± ± uncertainties in the IRAS data are 5-10%.

∗Flux densities from IRAS Catalog and are not color corrected – 38 –

Table 4

Table 4: Properties of Compact Sources in ULIRGS

Object 12.5µm 17.9µm 24.5µm 8.4 Ghz 2.2 µm % flux+ size %flux+ size %flux+ size size size in core ′′ in core ′′ in core ′′ ′′ ′′

IRAS05189 100 < 0.20 – – 100 <0.26 0.20 0.17 <0.2 × IRAS08572 80 <0.22∗ 100 <0.25 – – 0.09 0.07 <0.2 ∼ ∼ × UGC 5101 50 <0.25∗ 60 <0.25 – – 0.14 0.11 <0.2 × Mrk 231 100 <0.13∗ 85 <0.25 – – 0.07 0.06 <0.2 ≤ × Mrk 273N 67 <0.3 – – – – 0.32 0.18 0.3 × ∼ Mrk 273SW 50 <0.3 – – – – 0.18 0.11 <0.2 × Arp 220 E – – – – 25 – 0.32 0.19 <0.2 × Arp 220 W – – – – 75 0.40 0.21 0.14 0.4 × IRAS17208 100 2.0 – – – – 0.26 0.32 x 2 ×

+percentage of flux measured in 4′′ diameter beam

∗size measurement from LWS observations xsize measurement at 1.6 GHz from Martin et al. 1989 – 39 –

Table 5

Table 5: Infrared Sizes of ULRIGs

−2 Object size of mid-ir core L/A[L⊙Kpc ] ′′ pc max min

IRAS05189 <0.2 <166 5.6 1013 2.1 1013 × × IRAS08572 <0.22 <260 > 2.8 1013 > 6.0 1012 × × UGC 5101 <0.25 <205 > 1.1 1013 > 1.0 1012 × × Mrk 231 <0.13 <104 3.9 1014 1.0 1014 × × Mrk 273N <0.3 <240 > 2.2 1013 > 1.0 1012 × × Arp 220 0.40 140 6.0 1013 1.2 1013 × × IRAS17208 2.0 1600 1.2 1012 3.0 1010 × × – 40 –

Table 6

Table 6: Summary of Properties of ULIRGS

Object Compact core? Spectral Class MIR/NIR/Radio Optical Near IR ISO ASCA

IRAS05189 Y/Y/N Sey2 Sey – AGN IRAS08572 Y/Y/N LINER HII – SB UGC 5101 Y/Y/Y LINER – SB SB Mrk 231 Y/Y/Y Sey1 Sey1 AGN SB? Mrk 273 Y/Y/Y Sey2 Sey? AGN AGN Arp 220 Y/N/Y LINER – SB SB? IRAS17208 N/N/? HII HII SB –