Source author record

E. J. Murphy

E. J. Murphy appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

82works
7topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

82 published item(s)

preprint2020arXiv

The e-MERLIN Galaxy Evolution Survey (e-MERGE): Overview and Survey Description

We present an overview and description of the eMERLIN Galaxy Evolution survey (eMERGE) Data Release 1 (DR1), a large program of high-resolution 1.5 GHz radio observations of the GOODS-N field comprising $\sim140$ hours of observations with eMERLIN and $\sim40$ hours with the Very Large Array (VLA). We combine the long baselines of eMERLIN (providing high angular resolution) with the relatively closely-packed antennas of the VLA (providing excellent surface brightness sensitivity) to produce a deep 1.5 GHz radio survey with the sensitivity ($\sim 1.5μ$Jy beam$^{-1}$), angular resolution ($0.2"$--$0.7"$) and field-of-view ($\sim15' \times 15'$) to detect and spatially resolve star-forming galaxies and AGN at $z\gtrsim 1$. The goal of eMERGE is to provide new constraints on the deep, sub-arcsecond radio sky which will be surveyed by SKA1-mid. In this initial publication, we discuss our data analysis techniques, including steps taken to model in-beam source variability over a $\sim20$ year baseline and the development of new point spread function/primary beam models to seamlessly merge eMERLIN and VLA data in the $uv$ plane. We present early science results, including measurements of the luminosities and/or linear sizes of $\sim500$ galaxes selected at 1.5 GHz. In combination with deep Hubble Space Telescope observations, we measure a mean radio-to-optical size ratio of $r_{\rm eMERGE}/r_{\rm HST}\sim1.02\pm0.03$, suggesting that in most high-redshift galaxies, the $\sim$GHz continuum emission traces the stellar light seen in optical imaging. This is the first in a series of papers which will explore the $\sim$kpc-scale radio properties of star-forming galaxies and AGN in the GOODS-N field observed by eMERGE DR1.

preprint2020arXiv

The Star Formation in Radio Survey: 3 - 33 GHz Imaging of Nearby Galaxy Nuclei and Extranuclear Star-forming Regions

We present 3, 15, and 33 GHz imaging towards galaxy nuclei and extranuclear star-forming regions using the Karl G. Jansky Very Large Array as part of the Star Formation in Radio Survey. With $3-33$ GHz radio spectra, we measured the spectral indices and corresponding thermal (free-free) emission fractions for a sample of 335 discrete regions having significant detections in at least two radio bands. After removing 14 likely background galaxies, we find that the median thermal fraction at 33 GHz is $92 \pm 0.8\%$ with a median absolute deviation of $11\%$, when a two-component power-law model is adopted to fit the radio spectrum. Limiting the sample to 238 sources that are confidently identified as star-forming regions, and not affected by potential AGN contamination (i.e., having galactocentric radii $r_{\rm G} \geq 250$ pc), results in a median thermal fraction of $93 \pm 0.8 \%$ with a median absolute deviation of $10\%$. We further measure the thermal fraction at 33 GHz for 163 regions identified at 7" resolution to be $94 \pm 0.8 \%$ with a median absolute deviation of $8\%$. Together, these results confirm that free-free emission dominates the radio spectra of star-forming regions on scales up to $\sim$500 pc in normal star-forming galaxies. We additionally find a factor of $\sim$1.6 increase in the scatter of the measured spectral index and thermal fraction distributions as a function of decreasing galactocentric radius. This trend is likely reflective of the continuous star-formation activity occurring in the galaxy centers, resulting a larger contribution of diffuse nonthermal emission relative to star-forming regions in the disk.

preprint2019arXiv

Modeling Dust and Starlight in Galaxies Observed by Spitzer and Herschel: The KINGFISH Sample

Dust and starlight are modeled for the KINGFISH project galaxies. With data from 3.6 micron to 500 micron, models are strongly constrained. For each pixel in each galaxy we estimate (1) dust surface density; (2) q_PAH, the dust mass fraction in PAHs; (3) distribution of starlight intensities heating the dust; (4) luminosity emitted by the dust; and (5) dust luminosity from regions with high starlight intensity. The models successfully reproduce both global and resolved spectral energy distributions. We provide well-resolved maps for the dust properties. As in previous studies, we find q_PAH to be an increasing function of metallicity, above a threshold Z/Z_sol approx 0.15. Dust masses are obtained by summing the dust mass over the map pixels; these "resolved" dust masses are consistent with the masses inferred from model fits to the global photometry. The global dust-to-gas ratios obtained from this study correlate with galaxy metallicities. Systems with Z/Z_sol > 0.5 have most of their refractory elements locked up in dust, whereas when Z/Z_sol < 0.3 most of these elements tend to remain in the gas phase. Within galaxies, we find that q_PAH is suppressed in regions with unusually warm dust with nu L_nu(70 um) > 0.4L_dust. With knowledge of one long-wavelength flux density ratio (e.g., f_{160}/f_{500}), the minimum starlight intensity heating the dust (U_min) can be estimated to within ~50%. For the adopted dust model, dust masses can be estimated to within ~0.07 dex accuracy using the 500 micron luminosity nu L_nu(500) alone. There are additional systematic errors arising from the choice of dust model, but these are hard to estimate. These calibrated prescriptions may be useful for studies of high-redshift galaxies.

preprint2019arXiv

Star-forming Clumps in Local Luminous Infrared Galaxies

We present $HST$ narrow-band near-infrared imaging of Pa$α$ and Pa$β$ emission of 48 local Luminous Infrared Galaxies (LIRGs) from the Great Observatories All-Sky LIRG Survey (GOALS). These data allow us to measure the properties of 810 spatially resolved star-forming regions (59 nuclei and 751 extra-nuclear clumps), and directly compare their properties to those found in both local and high-redshift star-forming galaxies. We find that in LIRGs, the star-forming clumps have radii ranging from $\sim90-900$ pc and star formation rates (SFRs) of $\sim1\times10^{-3}$ to 10 M$_\odot$yr$^{-1}$, with median values for extra-nuclear clumps of 170 pc and 0.03 M$_\odot$yr$^{-1}$. The detected star-forming clumps are young, with a median stellar age of $8.7$ Myrs, and a median stellar mass of $5\times10^{5}$ M$_\odot$. The SFRs span the range of those found in normal local star-forming galaxies to those found in high-redshift star-forming galaxies at $\rm{z}=1-3$. The luminosity function of the LIRG clumps has a flatter slope than found in lower-luminosity, star-forming galaxies, indicating a relative excess of luminous star-forming clumps. In order to predict the possible range of star-forming histories and gas fractions, we compare the star-forming clumps to those measured in the MassiveFIRE high-resolution cosmological simulation. The star-forming clumps in MassiveFIRE cover the same range of SFRs and sizes found in the local LIRGs and have total gas fractions that extend from 10 to 90%. If local LIRGs are similar to these simulated galaxies, we expect future observations with ALMA will find a large range of gas fractions, and corresponding star formation efficiencies, among the star-forming clumps in LIRGs.

preprint2019arXiv

The Karl G. Jansky Very Large Array Sky Survey (VLASS). Science case and survey design

The Very Large Array Sky Survey (VLASS) is a synoptic, all-sky radio sky survey with a unique combination of high angular resolution ($\approx$2.5"), sensitivity (a 1$σ$ goal of 70 $μ$Jy/beam in the coadded data), full linear Stokes polarimetry, time domain coverage, and wide bandwidth (2-4 GHz). The first observations began in September 2017, and observing for the survey will finish in 2024. VLASS will use approximately 5500 hours of time on the Karl G. Jansky Very Large Array (VLA) to cover the whole sky visible to the VLA (Declination $>-40^{\circ}$), a total of 33,885 deg$^2$. The data will be taken in three epochs to allow the discovery of variable and transient radio sources. The survey is designed to engage radio astronomy experts, multi-wavelength astronomers, and citizen scientists alike. By utilizing an "on the fly" interferometry mode, the observing overheads are much reduced compared to a conventional pointed survey. In this paper, we present the science case and observational strategy for the survey, and also results from early survey observations.

preprint2016arXiv

A survey of the cold molecular gas in gravitationally lensed star-forming galaxies at z>2

Using the Australia Telescope Compact Array (ATCA), we conducted a survey of CO J=1-0 and J=2-1 line emission towards strongly lensed high-redshift dusty star forming galaxies (DSFGs) previously discovered with the South Pole Telescope (SPT). Our sample comprises 17 sources that had CO-based spectroscopic redshifts obtained with the Atacama Large Millimeter/submillimeter Array (ALMA) and the Atacama Pathfinder Experiment (APEX). We detect all sources with known redshifts in either CO J=1-0 or J=2-1. Twelve sources are detected in the 7-mm continuum. The derived CO luminosities imply gas masses in the range (0.5-11)x10^{10} M_sun and gas depletion timescales <200 Myr, using a CO to gas mass conversion factor alpha_CO=0.8 M_sun (K km/s pc^2)^{-1}. Combining the CO luminosities and dust masses, along with a fixed gas-to-dust ratio, we derive alpha_CO factors in the range 0.4-1.8, similar to what is found in other starbursting systems. We find small scatter in alpha_CO values within the sample, even though inherent variations in the spatial distribution of dust and gas in individual cases could bias the dust-based alpha_CO estimates. We find that lensing magnification factors based on the CO linewidth to luminosity relation (mu_CO) are highly unreliable, but particularly when mu<5. Finally, comparison of the gas and dynamical masses suggest that the average molecular gas fraction stays relatively constant at z=2-5 in the SPT DSFG sample.

preprint2016arXiv

An ALMA view of the interstellar medium of the z=4.77 lensed starburst SPT-S J213242-5802.9

We present ALMA detections of the [NII] 205$μ$m and CO(12-11) emission lines, and the tentative detection of [CI] $^3$P$_1$ - $^3$P$_0$ for the strongly lensed ($μ$=5.7$\pm$0.5) dusty, star-forming galaxy SPT-S J213242-5802.9 (hereafter SPT2132-58) at z=4.77. The [NII] and CO(12-11) lines are detected at 11.5 and 8.5 $σ$ level, respectively, by our band-6 observations. The [CI] line is detected at 3.2 $σ$ after a re-analysis of existing band-3 data. The [CI] luminosity implies a gas mass of 3.8$\pm$1.2$\times$10$^{10}$ M$_\odot$, and consequently a very short depletion timescale of 34$\pm$13 Myr and a CO-luminosity-to-gas-mass conversion factor $α_{\rm CO}$ of 1.0$\pm$0.3 M$_\odot$ (K km s$^{-1}$ pc$^{2}$)$^{-1}$. SPT2132-58 is an extreme starburst with an intrinsic star formation rate of 1100$\pm$200 M$_\odot$/yr. We find a [CII]/[NII] ratio of 26$\pm$7, which is the highest reported at z$>$4. This suggests that SPT2132-58 hosts an evolved interstellar medium (0.5 Z$_\odot$<Z<1.5 Z$_\odot$), which may be dominated by photodissociation regions. The CO(2-1) and CO(5-4) transitions have lower CO-to-far-infrared ratios than local and high-redshift samples, while CO(12-11) is similar to these samples, suggesting the presence of an additional very excited component or an AGN.

preprint2016arXiv

The redshift distribution of dusty star forming galaxies from the SPT survey

We use the Atacama Large Millimeter/submillimeter Array (ALMA) in Cycle 1 to determine spectroscopic redshifts of high-redshift dusty star-forming galaxies (DSFGs) selected by their 1.4mm continuum emission in the South Pole Telescope (SPT) survey. We present ALMA 3mm spectral scans between 84-114GHz for 15 galaxies and targeted ALMA 1mm observations for an additional eight sources. Our observations yield 30 new line detections from CO, [CI] , [NII] , H_2O and NH_3. We further present APEX [CII] and CO mid-J observations for seven sources for which only a single line was detected in spectral-scan data from ALMA Cycle 0 or Cycle 1. We combine the new observations with previously published and new mm/submm line and photometric data of the SPT-selected DSFGs to study their redshift distribution. The combined data yield 39 spectroscopic redshifts from molecular lines, a success rate of >85%. Our sample represents the largest data set of its kind today and has the highest spectroscopic completeness among all redshift surveys of high-z DSFGs. The median of the redshift distribution is z=3.9+/-0.4, and the highest-redshift source in our sample is at z=5.8. We discuss how the selection of our sources affects the redshift distribution, focusing on source brightness, selection wavelength, and strong gravitational lensing. We correct for the effect of gravitational lensing and find the redshift distribution for 1.4mm-selected sources with a median redshift of z=3.1+/-0.3. Comparing to redshift distributions selected at shorter wavelengths from the literature, we show that selection wavelength affects the shape of the redshift distribution.

preprint2015arXiv

A Cryogenic Space Telescope for Far-Infrared Astrophysics: A Vision for NASA in the 2020 Decade

Many of the transformative processes in the Universe have taken place in regions obscured by dust, and are best studied with far-IR spectroscopy. We present the Cryogenic-Aperture Large Infrared-Submillimeter Telescope Observatory (CALISTO), a 5-meter class, space-borne telescope actively cooled to 4 K, emphasizing moderate-resolution spectroscopy in the crucial 35 to 600 micron band. CALISTO will enable NASA and the world to study the rise of heavy elements in the Universe's first billion years, chart star formation and black hole growth in dust-obscured galaxies through cosmic time, and conduct a census of forming planetary systems in our region of the Galaxy. CALISTO will capitalize on rapid progress in both format and sensitivity of far-IR detectors. Arrays with a total count of a few 100,000 detector pixels will form the heart of a suite of imaging spectrometers in which each detector reaches the photon background limit. This document contains a large overview paper on CALISTO, as well as six 2-3 page scientific white papers, all prepared in response to NASA's Cosmic Origins Program Analysis Group (COPAG's) request for input on future mission concepts. The Far-IR Science Interest Group will meet from 3-5 June 2015 with the intention of reaching consensus on the architecture for the Far-IR Surveyor mission. This white paper describes one of the architectures to be considered by the community. One or more companion papers will describe alternative architectures.

preprint2015arXiv

CHANG-ES VI: Probing Supernova Energy Deposition in Spiral Galaxies Through Multi-Wavelength Relationships

How a galaxy regulates its SNe energy into different interstellar/circumgalactic medium components strongly affects galaxy evolution. Based on the JVLA D-configuration C- (6 GHz) and L-band (1.6 GHz) continuum observations, we perform statistical analysis comparing multi-wavelength properties of the CHANG-ES galaxies. The high-quality JVLA data and edge-on orientation enable us for the first time to include the halo into the energy budget for a complete radio-flux-limited sample. We find tight correlations of $L_{\rm radio}$ with the mid-IR-based SFR. The normalization of our $I_{\rm 1.6GHz}/{\rm W~Hz^{-1}}-{\rm SFR}$ relation is $\sim$2-3 times of those obtained for face-on galaxies, probably a result of enhanced IR extinction at high inclination. We also find tight correlations between $L_{\rm radio}$ and the SNe energy injection rate $\dot{E}_{\rm SN(Ia+CC)}$, indicating the energy loss via synchrotron radio continuum accounts for $\sim0.1\%$ of $\dot{E}_{\rm SN}$, comparable to the energy contained in CR electrons. The integrated C-to-L-band spectral index is $α\sim0.5-1.1$ for non-AGN galaxies, indicating a dominance by the diffuse synchrotron component. The low-scatter $L_{\rm radio}-{\rm SFR}$/$L_{\rm radio}-\dot{E}_{\rm SN (Ia+CC)}$ relationships have super-linear logarithmic slopes at $\sim2~σ$ in L-band ($1.132\pm0.067$/$1.175\pm0.102$) while consistent with linear in C-band ($1.057\pm0.075$/$1.100\pm0.123$). The super-linearity could be naturally reproduced with non-calorimeter models for galaxy disks. Using Chandra halo X-ray measurements, we find sub-linear $L_{\rm X}-L_{\rm radio}$ relations. These results indicate that the observed radio halo of a starburst galaxy is close to electron calorimeter, and a galaxy with higher SFR tends to distribute an increased fraction of SNe energy into radio emission (than X-ray).

preprint2015arXiv

Probing star formation in the dense environments of z~1 lensing halos aligned with dusty star-forming galaxies detected with the South Pole Telescope

We probe star formation in the environments of massive $\sim10^{13}\,M_{\odot}$ dark matter halos at redshifts of $z$$\sim$$1$. This star formation is linked to a sub-millimetre clustering signal which we detect in maps of the Planck High Frequency Instrument that are stacked at the positions of a sample of high-redshift ($z$$>$$2$) strongly-lensed dusty star-forming galaxies (DSFGs) selected from the South Pole Telescope (SPT) 2500 deg$^2$ survey. The clustering signal has sub-millimetre colours which are consistent with the mean redshift of the foreground lensing halos ($z$$\sim$$1$). We report a mean excess of star formation rate (SFR) compared to the field, of $(2700\pm700)\,M_{\odot}\,{yr}^{-1}$ from all galaxies contributing to this clustering signal within a radius of 3.5' from the SPT DSFGs. The magnitude of the Planck excess is in broad agreement with predictions of a current model of the cosmic infrared background. The model predicts that 80$\%$ of the excess emission measured by Planck originates from galaxies lying in the neighbouring halos of the lensing halo. Using Herschel maps of the same fields, we find a clear excess, relative to the field, of individual sources which contribute to the Planck excess. The mean excess SFR compared to the field is measured to be ($370\pm40)$$\,M_{\odot}\,{yr}^{-1}$ per resolved, clustered source. Our findings suggest that the environments around these massive $z$$\sim$$1$ lensing halos host intense star formation out to about $2\,$Mpc. The flux enhancement due to clustering should also be considered when measuring flux densities of galaxies in Planck data.

preprint2015arXiv

SKA studies of in-situ synchrotron radiation from molecular clouds

Observations of the properties of dense molecular clouds are critical in understanding the process of star-formation. One of the most important, but least understood, is the role of the magnetic fields. We discuss the possibility of using high-resolution, high-sensitivity radio observations with the SKA to measure for the first time the in-situ synchrotron radiation from these molecular clouds. If the cosmic-ray (CR) particles penetrate clouds as expected, then we can measure the B-field strength directly using radio data. So far, this signature has never been detected from the collapsing clouds themselves and would be a unique probe of the magnetic field. Dense cores are typically ~0.05 pc in size, corresponding to ~arcsec at ~kpc distances, and flux density estimates are ~mJy at 1 GHz. The SKA should be able to readily detect directly, for the first time, along lines-of-sight that are not contaminated by thermal emission or complex foreground/background synchrotron emission. Polarised synchrotron may also be detectable providing additional information about the regular/turbulent fields.

preprint2014arXiv

[CII] 158 $μ$m Emission as a Star Formation Tracer

The [CII] 157.74 $μ$m transition is the dominant coolant of the neutral interstellar gas, and has great potential as a star formation rate (SFR) tracer. Using the Herschel KINGFISH sample of 46 nearby galaxies, we investigate the relation of [CII] surface brightness and luminosity with SFR. We conclude that [CII] can be used for measurements of SFR on both global and kiloparsec scales in normal star-forming galaxies in the absence of strong active galactic nuclei (AGN). The uncertainty of the $Σ_{\rm [CII]}-Σ_{\rm SFR}$ calibration is $\pm$0.21 dex. The main source of scatter in the correlation is associated with regions that exhibit warm IR colors, and we provide an adjustment based on IR color that reduces the scatter. We show that the color-adjusted $Σ_{\rm[CII]}-Σ_{\rm SFR}$ correlation is valid over almost 5 orders of magnitude in $Σ_{\rm SFR}$, holding for both normal star-forming galaxies and non-AGN luminous infrared galaxies. Using [CII] luminosity instead of surface brightness to estimate SFR suffers from worse systematics, frequently underpredicting SFR in luminous infrared galaxies even after IR color adjustment (although this depends on the SFR measure employed). We suspect that surface brightness relations are better behaved than the luminosity relations because the former are more closely related to the local far-UV field strength, most likely the main parameter controlling the efficiency of the conversion of far-UV radiation into gas heating. A simple model based on Starburst99 population-synthesis code to connect SFR to [CII] finds that heating efficiencies are $1\%-3\%$ in normal galaxies.

preprint2014arXiv

A far-IR view of the starburst driven superwind in NGC 2146

NGC 2146, a nearby luminous infrared galaxy (LIRG), presents evidence for outflows along the disk minor axis in all gas phases (ionized, neutral atomic and molecular). We present an analysis of the multi-phase starburst driven superwind in the central 5 kpc as traced in spatially resolved spectral line observations, using far-IR Herschel PACS spectroscopy, to probe the effects on the atomic and ionized gas, and optical integral field spectroscopy to examine the ionized gas through diagnostic line ratios. We observe an increased ~250 km/s velocity dispersion in the [OI] 63 micron, [OIII] 88 micron, [NII] 122 micron and [CII] 158 micron fine-structure lines that is spatially coincident with high excitation gas above and below the disk. We model this with a slow ~200 km/s shock and trace the superwind to the edge of our field of view 2.5 kpc above the disk. We present new SOFIA 37 micron observations to explore the warm dust distribution, and detect no clear dust entrainment in the outflow. The stellar kinematics appear decoupled from the regular disk rotation seen in all gas phases, consistent with a recent merger event disrupting the system. We consider the role of the superwind in the evolution of NGC 2146 and speculate on the evolutionary future of the system. Our observations of NGC 2146 in the far-IR allow an unobscured view of the wind, crucial for tracing the superwind to the launching region at the disk center, and provide a local analog for future ALMA observations of outflows in high redshift systems.

preprint2014arXiv

ALMA Observations of Warm Dense Gas in NGC 1614 --- Breaking of Star Formation Law in the Central kpc

We present ALMA Cycle-0 observations of the CO (6-5) line emission and of the 435um dust continuum emission in the central kpc of NGC 1614, a local luminous infrared galaxy (LIRG) at a distance of 67.8 Mpc (1 arcsec = 329 pc). The CO emission is well resolved by the ALMA beam (0".26 x 0".20) into a circum-nuclear ring, with an integrated flux of f_{CO(6-5)} = 898 (+-153) Jy km/s, which is 63(+-12)% of the total CO(6-5) flux measured by Herschel. The molecular ring, located between 100pc < r < 350pc from the nucleus, looks clumpy and includes seven unresolved (or marginally resolved) knots with median velocity dispersion of 40 km/s. These knots are associated with strong star formation regions with Σ_{SFR} 100 M_\sun/yr/kpc^{2} and Σ_{Gas} 1.0E4 M_\sun/pc^{2}. The non-detections of the nucleus in both the CO (6-5) line emission and the 435um continuum rule out, with relatively high confidence, a Compton-thick AGN in NGC 1614. Comparisons with radio continuum emission show a strong deviation from an expected local correlation between Σ_{Gas} and Σ_{SFR}, indicating a breakdown of the Kennicutt-Schmidt law on the linear scale of 100 pc.

preprint2014arXiv

ALMA Observations of Warm Molecular Gas and Cold Dust in NGC 34

We present ALMA Cycle-0 observations of the CO (6-5) line emission (rest-frame frequency = 691.473 GHz) and of the 435$μm$ dust continuum emission in the nuclear region of NGC 34, a local luminous infrared galaxy (LIRG) at a distance of 84 Mpc (1" = 407 pc) which contains a Seyfert 2 active galactic nucleus (AGN) and a nuclear starburst. The CO emission is well resolved by the ALMA beam ($\rm 0."26\times 0."23$), with an integrated flux of $\rm f_{CO~(6-5)} = 1004\; (\pm 151) \; Jy\; km\; s^{-1}$. Both the morphology and kinematics of the CO (6-5) emission are rather regular, consistent with a compact rotating disk with a size of 200 pc. A significant emission feature is detected on the red-shifted wing of the line profile at the frequency of the $\rm H^{13}CN\; (8-7)$ line, with an integrated flux of $\rm 17.7 \pm 2.1 (random) \pm 2.7 (sysmatic)\; Jy\;km\; s^{-1}$. However, it cannot be ruled out that the feature is due to an outflow of warm dense gas with a mean velocity of $\rm 400\; km\; s^{-1}$. The continuum is resolved into an elongated configuration, and the observed flux corresponds to a dust mass of $\rm M_{dust} = 10^{6.97\pm 0.13}\; M_{sun}$. An unresolved central core ($\rm radius \simeq 50\; pc$) contributes $28\%$ of the continuum flux and $19\%$ of the CO (6-5) flux, consistent with insignificant contributions of the AGN to both emissions. Both the CO (6-5) and continuum spatial distributions suggest a very high gas column density ($\rm >= 10^4\; M_{sun}\; pc^{-2}$) in the nuclear region at $\rm radius <= 100\; pc$.

preprint2014arXiv

An ALMA survey of submillimetre galaxies in the Extended Chandra Deep Field South: radio properties and the far-infrared/radio correlation

We present a study of the radio properties of 870$μ$m-selected submillimetre galaxies (SMGs), observed at high resolution with ALMA in the Extended Chandra Deep Field South. From our initial sample of 76 ALMA SMGs, we detect 52 SMGs at $>3σ$ significance in VLA 1400MHz imaging, of which 35 are also detected at $>3σ$ in new 610MHz GMRT imaging. Within this sample of radio-detected SMGs, we measure a median radio spectral index $α_{610}^{1400} = -0.79 \pm 0.06$, (with inter-quartile range $α=[-1.16,-0.56]$) and investigate the far-infrared/radio correlation via the parameter $q_{\rm IR}$, the logarithmic ratio of the rest-frame 8-1000$μ$m flux and monochromatic radio flux. Our median $q_{\rm IR} = 2.56 \pm 0.05$ (inter-quartile range $q_{\rm IR}=[2.42,2.78]$) is higher than that typically seen in single-dish 870$μ$m-selected sources ($q_{\rm IR} \sim 2.4$), which may reflect the fact that our ALMA-based study is not biased to radio-bright counterparts, as previous samples were. Finally, we search for evidence that $q_{\rm IR}$ and $α$ evolve with age in a co-dependent manner, as predicted by starburst models: the data populate the predicted region of parameter space, with the stellar mass tending to increase along tracks of $q_{\rm IR}$ versus $α$ in the direction expected, providing the first observational evidence in support of these models.

preprint2014arXiv

Extended [CII] Emission in Local Luminous Infrared Galaxies

We present Herschel/PACS observations of extended [CII]157.7μm line emission detected on ~ 1 - 10 kpc scales in 60 local luminous infrared galaxies (LIRGs) from the Great Observatories All-sky LIRG Survey (GOALS). We find that most of the extra-nuclear emission show [CII]/FIR ratios >~ 4 x 10^-3, larger than the mean ratio seen in the nuclei, and similar to those found in the extended disks of normal star-forming galaxies and the diffuse inter-stellar medium (ISM) of our Galaxy. The [CII] "deficits" found in the most luminous local LIRGs are therefore restricted to their nuclei. There is a trend for LIRGs with warmer nuclei to show larger differences between their nuclear and extra-nuclear [CII]/FIR ratios. We find an anti-correlation between [CII]/FIR and the luminosity surface density, Σ_IR, for the extended emission in the spatially-resolved galaxies. However, there is an offset between this trend and that found for the LIRG nuclei. We use this offset to derive a beam filling-factor for the star-forming regions within the LIRG disks of ~ 6 % relative to their nuclei. We confront the observed trend to photo-dissociation region (PDR) models and find that the slope of the correlation is much shallower than the model predictions. Finally, we compare the correlation found between [CII]/FIR and Σ_IR with measurements of high-redshift starbursting IR-luminous galaxies.

preprint2014arXiv

Radio continuum properties of luminous infrared galaxies. Identifying the presence of an AGN in the radio

Luminous infrared galaxies are systems enshrouded in dust, which absorbs most of their optical/UV emission and re-radiates it in the mid- and far-infrared. Radio observations are largely unaffected by dust obscuration, enabling us to study the central regions of LIRGs in an unbiased manner. The main goal of this project is to examine how the radio properties of local LIRGs relate to their infrared spectral characteristics. Here we present an analysis of the radio continuum properties of a subset of the Great Observatories All-sky LIRG Survey (GOALS), which consists of 202 nearby systems (z<0.088). Our radio sample consists of 35 systems, or 46 individual galaxies, that were observed at both 1.49 and 8.44 GHz with the VLA with a resolution of about 1 arcsec (FWHM). The aim of the project is to use the radio imagery to probe the central kpc of these LIRGs in search of active galactic nuclei. We used the archival data at 1.49 and 8.44 GHz to create radio-spectral-index maps using the standard relation between flux density Sv and frequency v, S~v^-a, where a is the radio spectral index. By studying the spatial variations in a, we classified the objects as radio-AGN, radio-SB, and AGN/SB (a mixture). We identified the presence of an active nucleus using the radio morphology, deviations from the radio/infrared correlation, and spatially resolved spectral index maps, and then correlated this to the usual mid-infrared ([NeV]/[NeII] and [OIV]/[NeII] line ratios and EQW of the 6.2 um PAH feature) and optical (BPT diagram) AGN diagnostics. We find that 21 out of the 46 objects in our sample are radio-AGN, 9 are classified as starbursts (SB), and 16 are AGN/SB. After comparing to other AGN diagnostics we find 3 objects out of the 46 that are identified as AGN based on the radio analysis, but are not classified as such based on the mid-infrared and optical AGN diagnostics presented in this study.

preprint2014arXiv

SKA studies of nearby galaxies: star-formation, accretion processes and molecular gas across all environments

The SKA will be a transformational instrument in the study of our local Universe. In particular, by virtue of its high sensitivity (both to point sources and diffuse low surface brightness emission), angular resolution and the frequency ranges covered, the SKA will undertake a very wide range of astrophysical research in the field of nearby galaxies. By surveying vast numbers of nearby galaxies of all types with $μ$Jy sensitivity and sub-arcsecond angular resolutions at radio wavelengths, the SKA will provide the cornerstone of our understanding of star-formation and accretion activity in the local Universe. In this chapter we outline the key continuum and molecular line science areas where the SKA, both during phase-1 and when it becomes the full SKA, will have a significant scientific impact.

preprint2014arXiv

Synchrotron radiation from molecular clouds

Observations of the properties of dense molecular clouds are critical in understanding the process of star-formation. One of the most important, but least understood, is the role of the magnetic fields. We discuss the possibility of using high-resolution, high-sensitivity radio observations to measure the in-situ synchrotron radiation from these molecular clouds. If the cosmic-ray particles penetrate clouds as expected, then we can measure the B-field strength directly using radio data. So far, this signature has never been detected from the collapsing clouds themselves and would be a unique probe of the magnetic field. Dense cores are typically ~0.05 pc in size, corresponding to ~arcsec at ~kpc distances, and flux density estimates are ~mJy at 1 GHz. They should be detectable, for example with the Square Kilometre Array.

preprint2014arXiv

The Rest-Frame Submillimeter Spectrum of High-Redshift, Dusty, Star-Forming Galaxies

We present the average rest-frame spectrum of high-redshift dusty, star-forming galaxies from 250-770GHz. This spectrum was constructed by stacking ALMA 3mm spectra of 22 such sources discovered by the South Pole Telescope and spanning z=2.0-5.7. In addition to multiple bright spectral features of 12CO, [CI], and H2O, we also detect several faint transitions of 13CO, HCN, HNC, HCO+, and CN, and use the observed line strengths to characterize the typical properties of the interstellar medium of these high-redshift starburst galaxies. We find that the 13CO brightness in these objects is comparable to that of the only other z>2 star-forming galaxy in which 13CO has been observed. We show that the emission from the high-critical density molecules HCN, HNC, HCO+, and CN is consistent with a warm, dense medium with T_kin ~ 55K and n_H2 >~ 10^5.5 cm^-3. High molecular hydrogen densities are required to reproduce the observed line ratios, and we demonstrate that alternatives to purely collisional excitation are unlikely to be significant for the bulk of these systems. We quantify the average emission from several species with no individually detected transitions, and find emission from the hydride CH and the linear molecule CCH for the first time at high redshift, indicating that these molecules may be powerful probes of interstellar chemistry in high-redshift systems. These observations represent the first constraints on many molecular species with rest-frame transitions from 0.4-1.2mm in star-forming systems at high redshift, and will be invaluable in making effective use of ALMA in full science operations.

preprint2013arXiv

A Detailed Study of the Radio--FIR Correlation in NGC6946 with Herschel-PACS/SPIRE from KINGFISH

We derive the distribution of the synchrotron spectral index across NGC6946 and investigate the correlation between the radio continuum (synchrotron) and far-infrared (FIR) emission using the KINGFISH Herschel PACS and SPIRE data. The radio--FIR correlation is studied as a function of star formation rate, magnetic field strength, radiation field strength, and the total gas surface brightness. The synchrotron emission follows both star-forming regions and the so-called magnetic arms present in the inter-arm regions. The synchrotron spectral index is steepest along the magnetic arms ($α_n \sim 1$), while it is flat in places of giant H{\sc ii} regions and in the center of the galaxy ($α_n \sim 0.6-0.7$). The map of $α_n$ provides an observational evidence for aging and energy loss of cosmic ray electrons propagating in the disk of the galaxy. Variations in the synchrotron--FIR correlation across the galaxy are shown to be a function of both star formation and magnetic fields. We find that the synchrotron emission correlates better with cold rather than with warm dust emission, when the interstellar radiation field is the main heating source of dust. The synchrotron--FIR correlation suggests a coupling between the magnetic field and the gas density. NGC6946 shows a power-law behavior between the total (turbulent) magnetic field strength B and the star formation rate surface density $Σ_{\rm SFR}$ with an index of 0.14\,(0.16)$\pm$0.01. This indicates an efficient production of the turbulent magnetic field with the increasing gas turbulence expected in actively star forming regions. The scale-by-scale analysis of the synchrotron--FIR correlation indicates that the ISM affects the propagation of old/diffused cosmic ray electrons, resulting in a diffusion coefficient of $D_0=4.6\times 10^{28}$\,cm$^2$\,s$^{-1}$ for 2.2\,GeV CREs.

preprint2013arXiv

A Dust-Obscured Massive Maximum-Starburst Galaxy at a Redshift of 6.34

Massive present-day early-type (elliptical and lenticular) galaxies probably gained the bulk of their stellar mass and heavy elements through intense, dust-enshrouded starbursts - that is, increased rates of star formation - in the most massive dark matter halos at early epochs. However, it remains unknown how soon after the Big Bang such massive starburst progenitors exist. The measured redshift distribution of dusty, massive starbursts has long been suspected to be biased low in redshift owing to selection effects, as confirmed by recent findings of systems out to redshift z~5. Here we report the identification of a massive starburst galaxy at redshift 6.34 through a submillimeter color-selection technique. We unambiguously determined the redshift from a suite of molecular and atomic fine structure cooling lines. These measurements reveal a hundred billion solar masses of highly excited, chemically evolved interstellar medium in this galaxy, which constitutes at least 40% of the baryonic mass. A "maximum starburst" converts the gas into stars at a rate more than 2,000 times that of the Milky Way, a rate among the highest observed at any epoch. Despite the overall downturn of cosmic star formation towards the highest redshifts, it seems that environments mature enough to form the most massive, intense starbursts existed at least as early as 880 million years after the Big Bang.

preprint2013arXiv

ALMA Observations of SPT-Discovered, Strongly Lensed, Dusty, Star-Forming Galaxies

We present Atacama Large Millimeter/submillimeter Array (ALMA) 860 micrometer imaging of four high-redshift (z=2.8-5.7) dusty sources that were detected using the South Pole Telescope (SPT) at 1.4 mm and are not seen in existing radio to far-infrared catalogs. At 1.5 arcsec resolution, the ALMA data reveal multiple images of each submillimeter source, separated by 1-3 arcsec, consistent with strong lensing by intervening galaxies visible in near-IR imaging of these sources. We describe a gravitational lens modeling procedure that operates on the measured visibilities and incorporates self-calibration-like antenna phase corrections as part of the model optimization, which we use to interpret the source structure. Lens models indicate that SPT0346-52, located at z=5.7, is one of the most luminous and intensely star-forming sources in the universe with a lensing corrected FIR luminosity of 3.7 X 10^13 L_sun and star formation surface density of 4200 M_sun yr^-1 kpc^-2. We find magnification factors of 5 to 22, with lens Einstein radii of 1.1-2.0 arcsec and Einstein enclosed masses of 1.6-7.2x10^11 M_sun. These observations confirm the lensing origin of these objects, allow us to measure the their intrinsic sizes and luminosities, and demonstrate the important role that ALMA will play in the interpretation of lensed submillimeter sources.

preprint2013arXiv

ALMA redshifts of millimeter-selected galaxies from the SPT survey: The redshift distribution of dusty star-forming galaxies

Using the Atacama Large Millimeter/submillimeter Array (ALMA), we have conducted a blind redshift survey in the 3 mm atmospheric transmission window for 26 strongly lensd dusty star-forming galaxies (DSFGs) selected with the South Pole Telescope (SPT). The sources were selected to have S_1.4mm>20 mJy and a dust-like spectrum and, to remove low-z sources, not have bright radio (S_843MHz<6mJy) or far-infrared counterparts (S_100um<1 Jy, S_60um<200mJy). We robustly detect 44 line features in our survey, which we identify as redshifted emission lines of 12CO, 13CO, [CI], H2O, and H2O+. We find one or more spectral features in 23 sources yielding a ~90% detection rate for this survey; in 12 of these sources we detect multiple lines, while in 11 sources we detect only a single line. For the sources with only one detected line, we break the redshift degeneracy with additional spectroscopic observations if available, or infer the most likely line identification based on photometric data. This yields secure redshifts for ~70% of the sample. The three sources with no lines detected are tentatively placed in the redshift desert between 1.7<z<2.0. The resulting mean redshift of our sample is <z>=3.5. This finding is in contrast to the redshift distribution of radio-identified DSFGs, which have a significantly lower mean redshift of <z>=2.3 and for which only 10-15% of the population is expected to be at z>3. We discuss the effect of gravitational lensing on the redshift distribution and compare our measured redshift distribution to that of models in the literature.

preprint2013arXiv

Explaining the [CII]158um Deficit in Luminous Infrared Galaxies - First Results from a Herschel/PACS Study of the GOALS Sample

We present the first results of a survey of the [CII]158um emission line in 241 luminous infrared galaxies (LIRGs) comprising the Great Observatories All-sky Survey (GOALS) sample, obtained with the PACS instrument on board Herschel. The [CII] luminosities of the LIRGs in GOALS range from ~10^7 to 2x10^9 Lsun. We find that LIRGs show a tight correlation of [CII]/FIR with far-IR flux density ratios, with a strong negative trend spanning from ~10^-2 to 10^-4, as the average temperature of dust increases. We find correlations between the [CII]/FIR ratio and the strength of the 9.7um silicate absorption feature as well as with the luminosity surface density of the mid-IR emitting region (Sigma_MIR), suggesting that warmer, more compact starbursts have substantially smaller [CII]/FIR ratios. Pure star-forming (SF) LIRGs have a mean [CII]/FIR ~ 4x10^-3, while galaxies with low 6.2um PAH equivalent widths (EWs), indicative of the presence of active galactic nuclei (AGN), span the full range in [CII]/FIR. However, we show that even when only pure SF galaxies are considered, the [CII]/FIR ratio drops by an order of magnitude, from 10^-2 to 10^-3, with Sigma_MIR and Sigma_IR, implying that the [CII] luminosity is not a good indicator of the star formation rate (SFR) for most LIRGs, for it does not scale linearly with the warm dust emission. Moreover, even in LIRGs in which we detect an AGN in the mid-IR, the majority (2/3) of galaxies show [CII]/FIR >= 10^-3 typical of high 6.2um PAH EW sources, suggesting that most AGNs do not contribute significantly to the far-IR emission. We provide an empirical relation between the [CII]/FIR and the specific SFR (SSFR) for SF LIRGs. Finally, we present predictions for the starburst size based on the observed [CII] and far-IR luminosities which should be useful for comparing with results from future surveys of high-redshift galaxies with ALMA and CCAT.

preprint2013arXiv

High Spectral Resolution Measurement of the Sunyaev-Zel'dovich Effect Null with Z-Spec

The Sunyaev-Zel'dovich (SZ) effect spectrum crosses through a null where dT_CMB = 0 near nu_0 = 217 GHz. In a cluster of galaxies, nu_0 can be shifted from the canonical thermal SZ effect value by corrections to the SZ effect scattering due to the properties of the inter-cluster medium. We have measured the SZ effect in the hot galaxy cluster RX J1347.5-1145 with Z-Spec, an R ~ 300 grating spectrometer sensitive between 185 and 305 GHz. These data comprise a high spectral resolution measurement around the null of the SZ effect and clearly exhibit the transition from negative to positive dT_CMB over the Z-Spec band. The SZ null position is measured to be nu_0 = 225.8 \pm 2.5 (stat.) \pm 1.2 (sys.) GHz, which differs from the canonical null frequency by 3.0 sigma and is evidence for modifications to the canonical thermal SZ effect shape. Assuming the measured shift in nu_0 is due only to relativistic corrections to the SZ spectrum, we place the limit T_e = 17.1 \pm 5.3 keV from the zero-point measurement alone. By simulating the response of the instrument to the sky, we are able to generate likelihood functions in {y_0, T_e, v_pec} space. For v_pec = 0 km/s, we measure the best fitting SZ model to be y_0 = 4.6 (+0.6, -0.9) x 10^-4, T_e,0 = 15.2 (+12,-7.4) keV. When v_pec is allowed to vary, a most probable value of v_pec = +450 \pm 810 km/s is found.

preprint2013arXiv

Large gas reservoirs and free-free emission in two lensed star-forming galaxies at z=2.7

We report the detection of CO(1-0) line emission in the bright, lensed star-forming galaxies SPT-S 233227-5358.5 (z=2.73) and SPT-S 053816-5030.8 (z=2.78), using the Australia Telescope Compact Array (ATCA). Both galaxies were discovered in a large-area millimeter survey with the South Pole Telescope (SPT) and found to be gravitationally lensed by intervening structures. The measured CO intensities imply galaxies with molecular gas masses of (3.2 \pm 0.5)x10^10 (mu/15)^{-1}(X_CO/0.8) M_sun and (1.7 \pm 0.3)x10^10 (mu/20)^{-1}(X_CO/0.8) M_sun, and gas depletion timescales of 4.9x10^7 (X_CO/0.8) yr and 2.6x10^7 (X_CO/0.8) yr, respectively, where mu corresponds to the lens magnification and X_CO is the CO luminosity to gas mass conversion factor. In the case of SPT-S 053816-5030.8, we also obtained significant detections of the rest-frame 115.7 and 132.4 GHz radio continuum. Based on the radio to infrared spectral energy distribution and an assumed synchrotron spectral index, we find that 42 \pm 10 % and 55 \pm 13 % of the flux at rest-frame 115.7 and 132.4 GHz arise from free-free emission. We find a radio-derived intrinsic star formation rate (SFR) of 470 \pm 170 M_sun/yr, consistent within the uncertainties with the infrared estimate. Based on the morphology of this object in the source plane, the derived gas mass and the possible flattening of the radio spectral index towards low frequencies, we argue that SPT-S 053816-5030.8 exhibits properties compatible with a scaled-up local ultra-luminous infrared galaxy.

preprint2013arXiv

Mid-Infrared Properties of Nearby Luminous Infrared Galaxies I: Spitzer IRS Spectra for the GOALS Sample

The Great Observatories All-Sky LIRG Survey (GOALS) is a multiwavelength study of luminous infrared galaxies (LIRGs) in the local universe. Here we present low resolution Spitzer spectra covering 5-38um and provide a basic analysis of the mid-IR spectral properties for nearby LIRGs. In a companion paper, we discuss detailed fits to the spectra. The GOALS sample of 244 nuclei in 180 luminous and 22 ultraluminous IR galaxies represents a complete subset of the IRAS RBGS and covers a range of merger stages, morphologies and spectral types. The majority (>60%) of GOALS LIRGs have high 6.2um PAH equivalent widths (EQW > 0.4um) and low levels of silicate absorption (s_9.7um >-1.0). There is a general trend among the U/LIRGs for silicate depth and MIR slope to increase with LIR. U/LIRGs in the late stages of a merger also have on average steeper MIR slopes and higher levels of dust obscuration. Together these trends suggest that as gas & dust is funneled towards the center of a coalescing merger, the nuclei become more compact and obscured. The sources that depart from these correlations have very low PAH EQW (EQW < 0.1um) consistent with their MIR emission being dominated by an AGN. The most heavily dust obscured sources are the most compact in their MIR emission, suggesting that the obscuring (cool) dust is associated with the outer regions of the starburst. As the merger progresses a marked decline is seen for the fraction of high EQW (star formation dominated) sources while the fraction of composite sources increases but the fraction of AGN-dominated sources remains low. When compared to the MIR spectra of submillimeter galaxies (SMGs) at z~2, the average GOALS LIRG is more absorbed at 9.7um and has more PAH emission. However, when the AGN contributions to both the local LIRGs and the high-z SMGs are removed, the average local starbursting LIRG closely resembles the starbursting SMGs.

preprint2013arXiv

Radio and Mid-Infrared Properties of Compact Starbursts: Distancing Themselves from the Main Sequence

We investigate the relationship between 8.44\,GHz brightness temperatures and 1.4 to 8.44\,GHz radio spectral indices with 6.2\,$μ$m polycyclic aromatic hydrocarbon (PAH) emission and 9.7\,$μ$m silicate absorption features for a sample of 36 local luminous and ultra-luminous infrared galaxies. We find that galaxies having small 6.2\,$μ$m PAH equivalent widths (EQWs), which signal the presence of weak PAH emission and/or an excess of very hot dust, also have flat spectral indices. The three active galactic nuclei (AGN) identified through their excessively large 8.44\,GHz brightness temperatures are also identified as AGN via their small 6.2\,$μ$m PAH EQWs. We also find that the flattening of the radio spectrum increases with increasing silicate optical depth, 8.44\,GHz brightness temperature, and decreasing size of the radio source even after removing potential AGN, supporting the idea that compact starbursts show spectral flattening as the result of increased free-free absorption. These correlations additionally suggest that the dust obscuration in these galaxies must largely be coming from the vicinity of the compact starburst itself, and is not distributed throughout the (foreground) disk of the galaxy. Finally, we investigate the location of these infrared-bright systems relative to the main sequence (star formation rate vs. stellar mass) of star-forming galaxies in the local universe. We find that the radio spectral indices of galaxies flattens with increasing distance above the main sequence, or in other words, with increasing specific star formation rate. This indicates that galaxies located above the main sequence, having high specific star formation rates, are typically compact starbursts hosting deeply embedded star formation that becomes more optically thick in the radio and infrared with increased distance above the main sequence.

preprint2013arXiv

Shock Excited Molecules in NGC 1266: ULIRG conditions at the center of a Bulge Dominated Galaxy

We investigate the far infrared spectrum of NGC 1266, a S0 galaxy that contains a massive reservoir of highly excited molecular gas. Using the SPIRE-FTS, we detect the $^{12}$CO ladder up to J=(13-12), [C I] and [N II] lines, and also strong water lines more characteristic of UltraLuminous IR Galaxies (ULIRGs). The 12CO line emission is modeled with a combination of a low-velocity C-shock and a PDR. Shocks are required to produce the H2O and most of the high-J 12CO emission. Despite having an infrared luminosity thirty times less than a typical ULIRG, the spectral characteristics and physical conditions of the ISM of NGC 1266 closely resemble those of ULIRGs, which often harbor strong shocks and large-scale outflows.

preprint2013arXiv

Similarity Analysis of Macroecology

We perform a full similarity analysis of an idealized ecosystem using Buckingham's $Π$ theorem to obtain dimensionless similarity parameters given that some (non- unique) method exists that can differentiate different functional groups of individuals within an ecosystem. We then obtain the relationship between the similarity parameters under the assumptions of (i) that the ecosystem is in a dynamically balanced steady state and (ii) that these functional groups are connected to each other by the flow of resource. The expression that we obtain relates the level of complexity that the ecosystem can support to intrinsic macroscopic variables such as density, diversity and characteristic length scales for foraging or dispersal, and extrinsic macroscopic variables such as habitat size and the rate of supply of resource. This expression relates these macroscopic variables to each other, generating commonly observed macroecological patterns; these broad trends simply reflect the similarity property of ecosystems. We thus find that details of the ecosystem function are not required to obtain these broad macroecological patterns this may explain why they are ubiquitous. Departures from our relationship may indicate that the ecosystem is in a state of rapid change, i.e., abundance or diversity explosion or collapse. Our result provides normalised variables that can be used to isolate the trend in one ecosystem variable from another, providing a new method for isolating macroecological patterns in data. A dimensionless control parameter for ecosystem complexity emerges from our analysis and this will be a control parameter in dynamical models for ecosystems based on energy flow and conservation and will order the emergent behaviour of these models.

preprint2013arXiv

SPT 0538-50: Physical conditions in the ISM of a strongly lensed dusty star-forming galaxy at z=2.8

We present observations of SPT-S J053816-5030.8, a gravitationally-lensed dusty star forming galaxy (DSFG) at z = 2.7817, first discovered at millimeter wavelengths by the South Pole Telescope. SPT 0538-50 is typical of the brightest sources found by wide-field millimeter-wavelength surveys, being lensed by an intervening galaxy at moderate redshift (in this instance, at z = 0.441). We present a wide array of multi-wavelength spectroscopic and photometric data on SPT 0538-50, including data from ALMA, Herschel PACS and SPIRE, Hubble, Spitzer, VLT, ATCA, APEX, and the SMA. We use high resolution imaging from HST to de-blend SPT 0538-50, separating DSFG emission from that of the foreground lens. Combined with a source model derived from ALMA imaging (which suggests a magnification factor of 21 +/- 4), we derive the intrinsic properties of SPT 0538-50, including the stellar mass, far-IR luminosity, star formation rate, molecular gas mass, and - using molecular line fluxes - the excitation conditions within the ISM. The derived physical properties argue that we are witnessing compact, merger-driven star formation in SPT 0538-50, similar to local starburst galaxies, and unlike that seen in some other DSFGs at this epoch.

preprint2013arXiv

The Build-Up of Nuclear Stellar Cusps in Extreme Starburst Galaxies and Major Mergers

Nuclear stellar cusps are defined as central excess light component in the stellar light profiles of galaxies and are suggested to be stellar relics of intense compact starbursts in the central ~100-500 pc region of gas-rich major mergers. Here we probe the build-up of nuclear cusps during the actual starburst phase for a complete sample of Luminous Infrared Galaxy systems (85 LIRGs, with 11.4<log[LIR/L_sun]<12.5) in the GOALS sample. Cusp properties are derived via 2-dimensional fitting of the nuclear stellar light imaged in the near-infrared by the Hubble Space Telescope and have been combined with mid-IR diagnostics for AGN/starburst characterization. We find that nuclear stellar cusps are resolved in 76% of LIRGs (merger and non-interacting galaxies). The cusp strength and luminosity increases with far-IR luminosity (excluding AGN) and merger stage, confirming theoretical models that starburst activity is associated with the build-up of nuclear stellar cusps. Evidence for ultra compact nuclear starbursts is found in ~13% of LIRGs, which have a strong unresolved central NIR light component but no significant contribution of an AGN. The nuclear near-IR surface density (measured within 1 kpc radius) increases by a factor of ~5 towards late merger stages. A careful comparison to local early-type galaxies with comparable masses reveals (a) that local (U)LIRGs have a significantly larger cusp fraction and (b) that the majority of the cusp LIRGs have host galaxy luminosities (H-band) similar to core ellipticals which is roughly one order in magnitude larger than for cusp ellipticals.

preprint2012arXiv

A Study of Heating and Cooling of the ISM in NGC 1097 with Herschel-PACS and Spitzer-IRS

NGC 1097 is a nearby Seyfert 1 galaxy with a bright circumnuclear starburst ring, a strong large-scale bar and an active nucleus. We present a detailed study of the spatial variation of the far infrared (FIR) [CII]158um and [OI]63um lines and mid-infrared H2 emission lines as tracers of gas cooling, and of the polycyclic aromatic hydrocarbon (PAH) bands as tracers of the photoelectric heating, using Herschel-PACS, and Spitzer-IRS infrared spectral maps. We focus on the nucleus and the ring, and two star forming regions (Enuc N and Enuc S). We estimated a photoelectric gas heating efficiency ([CII]158um+[OI]63um)/PAH in the ring about 50% lower than in Enuc N and S. The average 11.3/7.7um PAH ratio is also lower in the ring, which may suggest a larger fraction of ionized PAHs, but no clear correlation with [CII]158μm/PAH(5.5 - 14um) is found. PAHs in the ring are responsible for a factor of two more [CII]158um and [OI]63um emission per unit mass than PAHs in the Enuc S. SED modeling indicates that at most 25% of the FIR power in the ring and Enuc S can come from high intensity photodissociation regions (PDRs), in which case G0 ~ 10^2.3 and nH ~ 10^3.5 cm^-3 in the ring. For these values of G0 and nH PDR models cannot reproduce the observed H2 emission. Much of the the H2 emission in the starburst ring could come from warm regions in the diffuse ISM that are heated by turbulent dissipation or shocks.

preprint2012arXiv

Characterizing Cosmic-Ray Propagation in Massive Star-Forming Regions: The Case of 30 Doradus and the Large Magellanic Cloud

Using infrared, radio, and gamma-ray data,we investigate the propagation characteristics of cosmic-ray (CR) electrons and nuclei in the 30 Doradus (30\,Dor) star-forming region in the Large Magellanic Cloud (LMC) using a phenomenological model based on the radio-far-infrared correlation within galaxies. Employing a correlation analysis, we derive an average propagation length of \sim 100-140 pc for \sim 3 GeV CR electrons resident in 30 Dor from consideration of the radio and infrared data. Assuming that the observed gamma-ray emission towards 30 Dor is associated with the star-forming region, and applying the same methodology to the infrared and gamma-ray data, we estimate a \sim 20 GeV propagation length of 200-320 pc for the CR nuclei. This is approximately twice as large as for\sim 3 GeV CR electrons, corresponding to a spatial diffusion coefficient that is \sim 4 times higher, scaling as (R/GV)^δ with δ\approx 0.7-0.8 depending on the smearing kernel used in the correlation analysis. This value is in agreement with the results found by extending the correlation analysis to include \sim 70 GeV CR nuclei traced by the 3-10 GeV gamma-ray data (δ\approx 0.66+/-0.23). Using the mean age of the stellar populations in 30 Dor and the results from our correlation analysis, we estimate a diffusion coefficient D_{R} \approx 0.9-1.0 \times10^27 (R/GV)^{0.7} cm^2 s^-1. We compare the values of the CR electron propagation length and surface brightness for 30 Dor and the LMC as a whole with those of entire disk galaxies. We find that the trend of decreasing average CR propagation distance with increasing disk-averaged star formation activity holds for the LMC, and extends down to single star-forming regions, at least for the case of 30 Dor.

preprint2012arXiv

Continuum Halos in Nearby Galaxies -- an EVLA Survey (CHANG-ES) -- I: Introduction to the Survey

We introduce a new survey to map the radio continuum halos of a sample of 35 edge-on spiral galaxies at 1.5 GHz and 6 GHz in all polarization products. The survey is exploiting the new wide bandwidth capabilities of the Karl G. Jansky Very Large Array (i.e. the Expanded Very Large Array, or EVLA) in a variety of array configurations (B, C, and D) in order to compile the most comprehensive data set yet obtained for the study of radio halo properties. This is the first survey of radio halos to include all polarization products. In this first paper, we outline the scientific motivation of the survey, the specific science goals, and the expected improvements in noise levels and spatial coverage from the survey. Our goals include investigating the physical conditions and origin of halos, characterizing cosmic ray transport and wind speed, measuring Faraday rotation and mapping the magnetic field, probing the in-disk and extraplanar far-infrared - radio continuum relation, and reconciling non-thermal radio emission with high-energy gamma-ray models. The sample size allows us to search for correlations between radio halos and other properties, including environment, star formation rate, and the presence of AGNs. In a companion paper (Paper II) we outline the data reduction steps and present the first results of the survey for the galaxy, NGC 4631.

preprint2012arXiv

Continuum Halos in Nearby Galaxies -- an EVLA Survey (CHANG-ES) -- II: First Results on NGC 4631

We present the first results from the CHANG-ES survey, a new survey of 35 edge-on galaxies to search for both in-disk as well as extra-planar radio continuum emission. The motivation and science case for the survey are presented in a companion paper (Paper I). In this paper (Paper II), we outline the observations and data reduction steps required for wide-band calibration and mapping of EVLA data, including polarization, based on C-array test observations of NGC 4631. With modest on-source observing times (30 minutes at 1.5 GHz and 75 minutes at 6 GHz for the test data) we have achieved best rms noise levels of 22 and 3.5 $μ$Jy beam$^{-1}$ at 1.5 GHz and 6 GHz, respectively. New disk-halo features have been detected, among them two at 1.5 GHz that appear as loops in projection. We present the first 1.5 GHz spectral index map of NGC 4631 to be formed from a single wide-band observation in a single array configuration. This map represents tangent slopes to the intensities within the band centered at 1.5 GHz, rather than fits across widely separated frequencies as has been done in the past and is also the highest spatial resolution spectral index map yet presented for this galaxy. The average spectral index in the disk is $\barα_{1.5 GHz}\,=\,-0.84\,\pm\,0.05$ indicating that the emission is largely non-thermal, but a small global thermal contribution is sufficient to explain a positive curvature term in the spectral index over the band. Two specific star forming regions have spectral indices that are consistent with thermal emission. Polarization results (uncorrected for internal Faraday rotation) are consistent with previous observations and also reveal some new features. On broad scales, we find strong support for the notion that magnetic fields constrain the X-ray emitting hot gas.

preprint2012arXiv

Cool Dust in the Outer Ring of NGC 1291

We examine Herschel Space Observatory images of one nearby prototypical outer ring galaxy, NGC 1291, and show that the ring becomes more prominent at wavelengths longer than 160um. The mass of cool dust in the ring dominates the total dust mass of the galaxy, accounting for at least 70% of it. The temperature of the emitting dust in the ring (T=19.5+/-0.3K) is cooler than that of the inner galaxy (T=25.7+/-0.7K). We discuss several explanations for the difference in dust temperature, including age and density differences in the stellar populations of the ring versus the bulge.

preprint2012arXiv

GOODS-Herschel: Ultra-deep XMM-Newton observations reveal AGN/star-formation connection

Models of galaxy evolution assume some connection between the AGN and star formation activity in galaxies. We use the multi-wavelength information of the CDFS to assess this issue. We select the AGNs from the 3Ms XMM-Newton survey and measure the star-formation rates of their hosts using data that probe rest-frame wavelengths longward of 20 um. Star-formation rates are obtained from spectral energy distribution fits, identifying and subtracting an AGN component. We divide the star-formation rates by the stellar masses of the hosts to derive specific star-formation rates (sSFR) and find evidence for a positive correlation between the AGN activity (proxied by the X-ray luminosity) and the sSFR for the most active systems with X-ray luminosities exceeding Lx=10^43 erg/s and redshifts z~1. We do not find evidence for such a correlation for lower luminosity systems or those at lower redshifts. We do not find any correlation between the SFR (or the sSFR) and the X-ray absorption derived from high-quality XMM-Newton spectra either, showing that the absorption is likely to be linked to the nuclear region rather than the host, while the star-formation is not nuclear. Comparing the sSFR of the hosts to the characteristic sSFR of star-forming galaxies at the same redshift we find that the AGNs reside mostly in main-sequence and starburst hosts, reflecting the AGN - sSFR connection. Limiting our analysis to the highest X-ray luminosity AGNs (X-ray QSOs with Lx>10^44 erg/s), we find that the highest-redshift QSOs (with z>2) reside predominantly in starburst hosts, with an average sSFR more than double that of the "main sequence", and we find a few cases of QSOs at z~1.5 with specific star-formation rates compatible with the main-sequence, or even in the "quiescent" region. (abridged)

preprint2012arXiv

GOODS-Herschel:dust attenuation properties of UV selected high redshift galaxies

We study dust attenuation at UV wavelengths at high redshift, where the UV is redshifted to the observed visible. In particular, we search for a UV bump and related implications for dust attenuation determinations. We use data in the CDFS, obtained in intermediate and broad band filters by the MUSYC project, to sample the UV rest-frame of 751 galaxies with 0.95<z<2.2. When available, Herschel/PACS data (GOODS-Herschel project), and Spitzer/MIPS measurements, are used to estimate the dust emission. The SED of each source is fit using the CIGALE code. The amount of dust attenuation and the dust attenuation curve are obtained as outputs of the SED fitting process, together with other parameters linked to the SFH. The global amount of dust attenuation at UV wavelengths is found to increase with stellar mass and to decrease as UV luminosity increases. A UV bump at 2175A is securely detected in 20% of the galaxies, and the mean amplitude of the bump for the sample is similar to that observed in the LMC supershell region. This amplitude is found to be lower in galaxies with very high SSFRs, and 90% of the galaxies exhibiting a secure bump are at z<1.5. The attenuation curve is confirmed to be steeper than that of local starburst galaxies for 20$% of the galaxies. The large dispersion found for these two parameters describing the attenuation law is likely to reflect a wide diversity of attenuation laws among galaxies. The relations between dust attenuation, IR-to-UV flux ratio, and the slope of the UV continuum are derived for the mean attenuation curve found for our sample. Deviations from the average trends are found to correlate with the age of the young stellar population and the shape of the attenuation curve.(abriged)

preprint2012arXiv

Measurements of CO redshifts with Z-Spec for lensed submillimeter galaxies discovered in the H-ATLAS survey

We present new observations from Z-Spec, a broadband 185-305 GHz spectrometer, of five sub-millimeter bright lensed sources selected from the Herschel Astrophysical Terahertz Large Area Survey (H-ATLAS) science demonstration phase (SDP) catalog. We construct a redshift finding algorithm using combinations of the signal-to-noise of all the lines falling in the Z-Spec bandpass to determine redshifts with high confidence, even in cases where the signal-to-noise in individual lines is low. We measure the dust continuum in all sources and secure CO redshifts for four out of five (z~1.5-3). In one source, SDP.17, we tentatively identify two independent redshifts and a water line, confirmed at z=2.308. Our sources have properties characteristic of dusty starburst galaxies, with magnification-corrected star formation rates of 10^(2-3) M_sun/yr. Lower limits for the dust masses (~a few 10^8 M_sun) and spatial extents (~1 kpc equivalent radius) are derived from the continuum spectral energy distributions, corresponding to dust temperatures between 54 and 69K. In the LTE approximation, we derive relatively low CO excitation temperatures (< 100 K) and optical depths (tau<1). Performing a non-LTE excitation analysis using RADEX, we find that the CO lines measured by Z-Spec (from J=4->3 to 10->9, depending on the galaxy) localize the best solutions to either a high-temperature / low-density region, or a low-temperature / high-density region near the LTE solution, with the optical depth varying accordingly. Observations of additional CO lines, CO(1-0) in particular, are needed to constrain the non-LTE models.

preprint2012arXiv

Modeling Dust and Starlight in Galaxies Observed by Spitzer and Herschel: NGC 628 and NGC 6946

We characterize the dust in NGC628 and NGC6946, two nearby spiral galaxies in the KINGFISH sample. With data from 3.6um to 500um, dust models are strongly constrained. Using the Draine & Li (2007) dust model, (amorphous silicate and carbonaceous grains), for each pixel in each galaxy we estimate (1) dust mass surface density, (2) dust mass fraction contributed by polycyclic aromatic hydrocarbons (PAH)s, (3) distribution of starlight intensities heating the dust, (4) total infrared (IR) luminosity emitted by the dust, and (5) IR luminosity originating in regions with high starlight intensity. We obtain maps for the dust properties, which trace the spiral structure of the galaxies. The dust models successfully reproduce the observed global and resolved spectral energy distributions (SEDs). The overall dust/H mass ratio is estimated to be 0.0082+/-0.0017 for NGC628, and 0.0063+/-0.0009 for NGC6946, consistent with what is expected for galaxies of near-solar metallicity. Our derived dust masses are larger (by up to a factor 3) than estimates based on single-temperature modified blackbody fits. We show that the SED fits are significantly improved if the starlight intensity distribution includes a (single intensity) "delta function" component. We find no evidence for significant masses of cold dust T<12K. Discrepancies between PACS and MIPS photometry in both low and high surface brightness areas result in large uncertainties when the modeling is done at PACS resolutions, in which case SPIRE, MIPS70 and MIPS160 data cannot be used. We recommend against attempting to model dust at the angular resolution of PACS.

preprint2012arXiv

Resolving the Far-IR Line Deficit: Photoelectric Heating and Far-IR Line Cooling in NGC 1097 and NGC 4559

The physical state of interstellar gas and dust is dependent on the processes which heat and cool this medium. To probe heating and cooling of the ISM over a large range of infrared surface brightness, on sub-kiloparsec scales, we employ line maps of [C \ii] 158 $μ$m, [O \one] 63 $μ$m, and [N \ii] 122 $μ$m in NGC 1097 and NGC 4559, obtained with the PACS spectrometer onboard {\it Herschel}. We matched new observations to existing Spitzer-IRS data that trace the total emission of polycyclic aromatic hydrocarbons (PAHs). We confirm at small scales in these galaxies that the canonical measure of photoelectric heating efficiency, ([C \ii] + [O \one])/TIR, decreases as the far-infrared color, $νf_ν$(70 $μ$m)/$νf_ν$(100 $μ$m), increases. In contrast, the ratio of far-infrared (far-IR) cooling to total PAH emission, ([C \ii] + [O \one])/PAH, is a near constant $\sim$6% over a wide range of far-infrared color, 0.5 \textless\ $νf_ν$(70 $μ$m)/$νf_ν$(100 $μ$m) $\lesssim$ 0.95. In the warmest regions, where $νf_ν$(70 $μ$m)/$νf_ν$(100 $μ$m) $\gtrsim$ 0.95, the ratio ([C \ii] + [O \one])/PAH drops rapidly to 4%. We derived representative values of the local UV radiation density, $G_0$, and the gas density, $n_H$, by comparing our observations to models of photodissociation regions. The ratio $G_0/n_H$, derived from fine-structure lines, is found to correlate with the mean dust-weighted starlight intensity, $<U>$ derived from models of the IR SED. Emission from regions that exhibit a line deficit is characterized by an intense radiation field, indicating that small grains are susceptible to ionization effects. We note that there is a shift in the 7.7 / 11.3 $μ$m PAH ratio in regions that exhibit a deficit in ([C \ii] + [O \one])/PAH, suggesting that small grains are ionized in these environments.

preprint2012arXiv

Submillimeter Observations of Millimeter Bright Galaxies Discovered by the South Pole Telescope

We present APEX SABOCA 350micron and LABOCA 870micron observations of 11 representative examples of the rare, extremely bright (S_1.4mm > 15mJy), dust-dominated millimeter-selected galaxies recently discovered by the South Pole Telescope (SPT). All 11 sources are robustly detected with LABOCA with 40 < S_870micron < 130mJy, approximately an order of magnitude higher than the canonical submillimeter galaxy (SMG) population. Six of the sources are also detected by SABOCA at >3sigma, with the detections or upper limits providing a key constraint on the shape of the spectral energy distribution (SED) near its peak. We model the SEDs of these galaxies using a simple modified blackbody and perform the same analysis on samples of SMGs of known redshift from the literature. These calibration samples inform the distribution of dust temperature for similar SMG populations, and this dust temperature prior allows us to derive photometric redshift estimates and far infrared luminosities for the sources. We find a median redshift of <z> = 3.0, higher than the <z> = 2.2 inferred for the normal SMG population. We also derive the apparent size of the sources from the temperature and apparent luminosity, finding them to appear larger than our unlensed calibration sample, which supports the idea that these sources are gravitationally magnified by massive structures along the line of sight.

preprint2012arXiv

The Star Formation in Radio Survey: GBT 33 GHz Observations of Nearby Galaxy Nuclei and Extranuclear Star-Forming Regions

We present 33\,GHz photometry of 103 galaxy nuclei and extranuclear star-forming complexes taken with the Green Bank Telescope (GBT) as part of the Star Formation in Radio Survey (SFRS). Among the sources without evidence for an AGN, and also having lower frequency radio data, we find a median thermal fraction at 33GHz of ~76% with a dispersion of ~24%. For all sources resolved on scales <0.5kpc, the thermal fraction is even larger, being >90%. This suggests that the rest-frame 33GHz emission provides a sensitive measure of the ionizing photon rate from young star-forming regions, thus making it a robust star formation rate indicator. Taking the 33GHz star formation rates as a reference, we investigate other empirical calibrations relying on different combinations of warm 24μm dust, total infrared (IR; 8-1000μm), Hα line, and far-UV continuum emission. The recipes derived here generally agree with others found in the literature, albeit with a large dispersion that most likely stems from a combination of effects. Comparing the 33GHz to total IR flux ratios as a function of the radio spectral index, measured between 1.7 and 33GHz, we find that the ratio increases as the radio spectral index flattens which does not appear to be a distance effect. Consequently, the ratio of non-thermal to total IR emission appears relatively constant, suggesting only moderate variations in the cosmic-ray electron injection spectrum and ratio of synchrotron to total cooling processes among star-forming complexes. Assuming that this trend solely arises from an increase in the thermal fraction sets a maximum on the scatter of the non-thermal spectral indices among the star-forming regions of σ_α^{NT} < 0.13.

preprint2011arXiv

An Accounting of the Dust-Obscured Star Formation and Accretion Histories Over the Last ~11~Billion Years

(Abridged) We report on an accounting of the star formation and accretion driven energetics of 24um detected sources in GOODS North. For sources having infrared (IR; 8-1000um) luminosities >3x10^12 L_sun when derived by fitting local SEDs to 24um photometry alone, we find these IR luminosity estimates to be a factor of ~4 times larger than those estimated when the SED fitting includes additional 16 and 70um data (and in some cases mid-infrared spectroscopy and 850um data). This discrepancy arises from the fact that high luminosity sources at z>>0 appear to have far- to mid-infrared ratios, as well as aromatic feature equivalent widths, typical of lower luminosity galaxies in the local Universe. Using our improved estimates for IR luminosity and AGN contributions, we investigate the evolution of the IR luminosity density versus redshift arising from star formation and AGN processes alone. We find that, within the uncertainties, the total star formation driven IR luminosity density is constant between 1.15 < z < 2.35, although our results suggest a slightly larger value at z>2. AGN appear to account for <18% of the total IR luminosity density integrated between 0< z < 2.35, contributing <25% at each epoch. LIRG appear to dominate the star formation rate (SFR) density along with normal star-forming galaxies (L_IR < 10^11 L_sun) between 0.6 < z < 1.15. Once beyond z >2, the contribution from ultraluminous infrared galaxies ULIRGs becomes comparable with that of LIRGs. Using our improved IR luminosity estimates, we find existing calibrations for UV extinction corrections based on measurements of the UV spectral slope typically overcorrect UV luminosities by a factor of ~2, on average, for our sample of 24um-selected sources; accordingly we have derived a new UV extinction correction more appropriate for our sample.

preprint2011arXiv

Calibrating Extinction-Free Star Formation Rate Diagnostics with 33GHz Free-Free Emission in NGC6946

Abridged: Using free-free emission measured in the Ka-band (26-40GHz) for 10 star-forming regions in the nearby galaxy NGC6946, including its starbursting nucleus, we compare a number of SFR diagnostics that are typically considered to be unaffected by interstellar extinction: i.e., non-thermal radio (i.e., 1.4GHz), total infrared (IR; 8-1000um), and warm dust (i.e., 24um) emission, along with the hybrid (obscured + unobscured) indicators of Hα+24um and UV+IR. The 33GHz free-free emission is assumed to provide the most accurate measure of the current SFR. Among the extranuclear star-forming regions, the 24um, Hα+24um and UV+IR SFR calibrations are in good agreement with the 33GHz free-free SFRs. However, each of the SFR calibrations relying on some form of dust emission overestimate the nuclear SFR by a factor of ~2. This is more likely the result of excess dust heating through an accumulation of non-ionizing stars associated with an extended episode of star formation in the nucleus rather than increased competition for ionizing photons by dust. SFR calibrations using the non-thermal radio continuum yield values which only agree with the free-free SFRs for the nucleus, and underestimate the SFRs from the extranuclear star-forming regions by a factor of ~2. This result likely arises from the CR electrons decaying within the starburst region with negligible escape compared to the young extranuclear star-forming regions. Finally, we find that the SFRs estimated using the total 33GHz emission agree well with the free-free SFRs due to the large thermal fractions present at these frequencies even when local diffuse backgrounds are not removed. Thus, rest-frame 33GHz observations may act as a reliable method to measure the SFRs of galaxies at increasingly high redshift without the need of ancillary radio data to account for the non-thermal emission.

preprint2011arXiv

Detection of an ultra-bright submillimeter galaxy in the Subaru/XMM-Newton Deep Field using AzTEC/ASTE

We report the detection of an extremely bright ($\sim$37 mJy at 1100 $μ$m and $\sim$91 mJy at 880 $μ$m) submillimeter galaxy (SMG), AzTEC-ASTE-SXDF1100.001 (hereafter referred to as SXDF1100.001 or Orochi), discovered in 1100 $μ$m observations of the Subaru/XMM-Newton Deep Field using AzTEC on ASTE. Subsequent CARMA 1300 $μ$m and SMA 880 $μ$m observations successfully pinpoint the location of Orochi and suggest that it has two components, one extended (FWHM of $\sim$ 4$^{\prime\prime}$) and one compact (unresolved). Z-Spec on CSO has also been used to obtain a wide band spectrum from 190 to 308 GHz, although no significant emission/absorption lines are found. The derived upper limit to the line-to-continuum flux ratio is 0.1--0.3 (2 $σ$) across the Z-Spec band. Based on the analysis of the derived spectral energy distribution from optical to radio wavelengths of possible counterparts near the SMA/CARMA peak position, we suggest that Orochi is a lensed, optically dark SMG lying at $z \sim 3.4$ behind a foreground, optically visible (but red) galaxy at $z \sim 1.4$. The deduced apparent (i.e., no correction for magnification) infrared luminosity ($L_{\rm IR}$) and star formation rate (SFR) are $6 \times 10^{13}$ $L_{\odot}$ and 11000 $M_{\odot}$ yr$^{-1}$, respectively, assuming that the $L_{\rm IR}$ is dominated by star formation. These values suggest that Orochi will consume its gas reservoir within a short time scale ($3 \times 10^{7}$ yr), which is indeed comparable to those in extreme starbursts like the centres of local ULIRGs.

preprint2011arXiv

Herschel Far-Infrared and Sub-millimeter Photometry for the KINGFISH Sample of Nearby Galaxies

New far-infrared and sub-millimeter photometry from the Herschel Space Observatory is presented for 61 nearby galaxies from the Key Insights on Nearby Galaxies: A Far-Infrared Survey with Herschel (KINGFISH) sample. The spatially-integrated fluxes are largely consistent with expectations based on Spitzer far-infrared photometry and extrapolations to longer wavelengths using popular dust emission models. Dwarf irregular galaxies are notable exceptions, as already noted by other authors, as their 500um emission shows evidence for a sub-millimeter excess. In addition, the fraction of dust heating attributed to intense radiation fields associated with photo-dissociation regions is found to be (21+/-4)% larger when Herschel data are included in the analysis. Dust masses obtained from the dust emission models of Draine & Li are found to be on average nearly a factor of two higher than those based on single-temperature modified blackbodies, as single blackbody curves do not capture the full range of dust temperatures inherent to any galaxy. The discrepancy is largest for galaxies exhibiting the coolest far-infrared colors.

preprint2011arXiv

KINGFISH -- Key Insights on Nearby Galaxies: A Far-Infrared Survey with Herschel: Survey Description and Image Atlas

The KINGFISH project (Key Insights on Nearby Galaxies: a Far-Infrared Survey with Herschel) is an imaging and spectroscopic survey of 61 nearby (d < 30 Mpc) galaxies, chosen to cover a wide range of galaxy properties and local interstellar medium (ISM) environments found in the nearby Universe. Its broad goals are to characterize the ISM of present-day galaxies, the heating and cooling of their gaseous and dust components, and to better understand the physical processes linking star formation and the ISM. KINGFISH is a direct descendant of the Spitzer Infrared Nearby Galaxies Survey (SINGS), which produced complete Spitzer imaging and spectroscopic mapping and a comprehensive set of multi-wavelength ancillary observations for the sample. The Herschel imaging consists of complete maps for the galaxies at 70, 100, 160, 250, 350, and 500 microns. The spectal line imaging of the principal atomic ISM cooling lines ([OI]63um, [OIII]88um, [NII]122,205um, and [CII]158um) covers the subregions in the centers and disks that already have been mapped in the mid-infrared with Spitzer. The KINGFISH and SINGS multi-wavelength datasets combined provide panchromatic mapping of the galaxies sufficient to resolve individual star-forming regions, and tracing the important heating and cooling channels of the ISM, across a wide range of local extragalactic ISM environments. This paper summarizes the scientific strategy for KINGFISH, the properties of the galaxy sample, the observing strategy, and data processing and products. It also presents a combined Spitzer and Herschel image atlas for the KINGFISH galaxies, covering the wavelength range 3.6 -- 500 microns. All imaging and spectroscopy data products will be released to the Herschel user generated product archives.

preprint2011arXiv

Redshift Determination and CO Line Excitation Modeling for the Multiply-Lensed Galaxy HLSW-01

We report on the redshift measurement and CO line excitation of HERMES J105751.1+573027 (HLSW-01), a strongly lensed submillimeter galaxy discovered in Herschel/SPIRE observations as part of the Herschel Multi-tiered Extragalactic Survey (HerMES). HLSW-01 is an ultra-luminous galaxy with an intrinsic far-infrared luminosity of 1.4x10^(13) solar luminosities, and is lensed by a massive group of galaxies into at least four images with a total magnification of 10.9+/-0.7. With the 100 GHz instantaneous bandwidth of the Z-Spec instrument on the Caltech Submillimeter Observatory, we robustly identify a redshift of z=2.958+/-0.007 for this source, using the simultaneous detection of four CO emission lines (J = 7-6, J = 8-7, J = 9-8, and J = 10-9). Combining the measured line fluxes for these high-J transitions with the J = 1-0, J = 3-2 and J = 5-4 line fluxes measured with the Green Bank Telescope, the Combined Array for Research in Millimeter Astronomy, and the Plateau de Bure Interferometer, respectively, we model the physical properties of the molecular gas in this galaxy. We find that the full CO spectral line energy distribution is well described by warm, moderate-density gas with Tkin = 86-235 K and n(H2) = (1.1-3.5)x10^3 cm^(-3). However, it is possible that the highest-J transitions are tracing a small fraction of very dense gas in molecular cloud cores, and two-component models that include a warm/dense molecular gas phase with Tkin ~ 200 K, n(H2) ~ 10^5 cm^(-3) are also consistent with these data. Higher signal-to-noise measurements of the J(upper) > 7 transitions with high spectral resolution, combined with high spatial resolution CO maps, are needed to improve our understanding of the gas excitation, morphology, and dynamics of this interesting high-redshift galaxy.

preprint2011arXiv

Spitzer IRS Spectral Mapping of the Toomre Sequence: Spatial Variations of PAH, Gas, and Dust Properties in Nearby Major Mergers

We have mapped the key mid-IR diagnostics in eight major merger systems of the Toomre Sequence (NGC4676, NGC7592, NGC6621, NGC2623, NGC6240, NGC520, NGC3921, and NGC7252) using the Spitzer Infrared Spectrograph (IRS). With these maps, we explore the variation of the ionized-gas, PAH, and warm-gas (H_2) properties across the sequence and within the galaxies. While the global PAH interband strength and ionized gas flux ratios ([Ne III]/[Ne II]) are similar to those of normal star forming galaxies, the distribution of the spatially resolved PAH and fine structure line flux ratios is significant different from one system to the other. Rather than a constant H_2/PAH flux ratio, we find that the relation between the H_2 and PAH fluxes is characterized by a power law with a roughly constant exponent (0.61+/-0.05) over all merger components and spatial scales. While following the same power law on local scales, three galaxies have a factor of ten larger integrated (i.e. global) H_2/PAH flux ratio than the rest of the sample, even larger than what it is in most nearby AGNs. These findings suggest a common dominant excitation mechanism for H_2 emission over a large range of global H_2/PAH flux ratios in major mergers. Early merger systems show a different distribution between the cold (CO J=1-0) and warm (H_2) molecular gas component, which is likely due to the merger interaction. Strong evidence for buried star formation in the overlap region of the merging galaxies is found in two merger systems (NGC6621 and NGC7592) as seen in the PAH, [Ne II], [Ne III], and warm gas line emission, but with no apparent corresponding CO (J=1-0) emission. Our findings also demonstrate that the variations of the physical conditions within a merger are much larger than any systematic trends along the Toomre Sequence.

preprint2011arXiv

The Dense Molecular Gas in the Circumnuclear Disk of NGC 1068

We present a 190-307 GHz broadband spectrum obtained with Z-Spec of NGC 1068 with new measurements of molecular rotational transitions. After combining our measurements with those previously published and considering the specific geometry of this Seyfert 2 galaxy, we conduct a multi-species Bayesian likelihood analysis of the density, temperature, and relative molecular abundances of HCN, HNC, CS, and HCO+. We find that these molecules trace warm (T > 100 K) gas of H2 number densities 10^4.2 - 10^4.9 cm^-3. Our models also place strong constraints on the column densities and relative abundances of these molecules, as well as on the total mass in the circumnuclear disk. Using the uniform calibration afforded by the broad Z-Spec bandpass, we compare our line ratios to X-ray dominated region (XDR) and photon-dominated region models. The majority of our line ratios are consistent with the XDR models at the densities indicated by the likelihood analysis, lending substantial support to the emerging interpretation that the energetics in the circumnuclear disk of NGC 1068 are dominated by accretion onto an active galactic nucleus.

preprint2011arXiv

The Displaced Dusty ISM of NGC 3077: Tidal Stripping in the M 81 Triplet

We present the detection of extended (~30 kpc^2) dust emission in the tidal \hi\ arm near NGC 3077 (member of the M\,81 triplet) using SPIRE on board Herschel. Dust emission in the tidal arm is typically detected where the \hi\ column densities are >10^21 cm^-2. The SPIRE band ratios show that the dust in the tidal arm is significantly colder (~13 K) than in NGC 3077 itself (~31 K), consistent with the lower radiation field in the tidal arm. The total dust mass in the tidal arm is ~1.8 x 10^6 M_sun (assuming beta=2), i.e. substantially larger than the dust mass associated with NGC 3077 (~2 x 10^5 M_sun). Where dust is detected, the dust-to-gas ratio is 6+/-3 x 10^-3, consistent within the uncertainties with what is found in NGC 3077 and nearby spiral galaxies with Galactic metallicities. The faint HII regions in the tidal arm can not be responsible for the detected enriched material and are not the main source of the dust heating in the tidal arm. We conclude that the interstellar medium (atomic HI, molecules and dust) in this tidal feature was pre-enriched and stripped off NGC 3077 during its recent interaction (~3 x 10^8 yr ago) with M 82 and M 81. This implies that interaction can efficiently remove heavy elements and enriched material (dust, molecular gas) from galaxies. As interactions were more frequent at large lookback times, it is conceivable that they could substantially contribute (along with galactic outflows) to the enrichment of the intergalactic medium.

preprint2011arXiv

The Water Vapor Spectrum of APM 08279+5255: X-Ray Heating and Infrared Pumping over Hundreds of Parsecs

We present the rest-frame 200--320 \mm\ spectrum of the z=3.91 quasar \apm, obtained with Z-Spec at the Caltech Submillimeter Observatory. In addition to the \jeight\ to \jthirteen\ CO rotational transitions which dominate the CO cooling, we find six transitions of water originating at energy levels ranging up to 643 K. Most are first detections at high redshift, and we have confirmed one transition with CARMA. The CO cooling is well-described by our XDR model, assuming L$_{\rm 1-100\,keV}\sim1\times10^{46}\rm\,erg\,s^{-1}$, and that the gas is distributed over a 550-pc sizescale, per the now-favored $μ$=4 lensing model. The total observed cooling in water corresponds to 6.5$\times10^{9}$ \ls, comparable to that of CO. We compare the water spectrum with that of Mrk 231, finding that the intensity ratios among the high-lying lines are similar, but with a total luminosity scaled up by a factor of $\sim$50. Using this scaling, we estimate an average water abundance relative to \hh\ of 1.4$\times10^{-7}$, a good match to the prediction of the chemical network in the XDR model. As with Mrk 231, the high-lying water transitions are excited radiatively via absorption in the rest-frame far-infrared, and we show that the powerful dust continuum in \apm\ is more than sufficient to pump this massive reservoir of warm water vapor.

preprint2010arXiv

Balancing the Energy Budget: Star-Formation versus AGN in High Redshift Infrared Luminous Galaxies

We present deep {\it Spitzer} mid-infrared spectroscopy, along with 16, 24, 70, and 850\,$\micron$\ photometry, for 22 galaxies located in the Great Observatories Origins Deep Survey-North (GOODS-N) field. The sample spans a redshift range of $0.6\la z \la 2.6$, 24~$μ$m flux densities between $\sim$0.2$-$1.2 mJy, and consists of submillimeter galaxies (SMGs), X-ray or optically selected active galactic nuclei (AGN), and optically faint ($z_{AB}>25$\,mag) sources. We find that infrared (IR; $8-1000~\micron$) luminosities derived by fitting local spectral energy distributions (SEDs) with 24~$\micron$ photometry alone are well matched to those when additional mid-infrared spectroscopic and longer wavelength photometric data is used for galaxies having $z\la1.4$ and 24~$\micron$-derived IR luminosities typically $\la 3\times 10^{12}~L_{\sun}$. However, for galaxies in the redshift range between $1.4\la z \la 2.6$, typically having 24~$\micron$-derived IR luminosities $\ga 3\times 10^{12}~L_{\sun}$, IR luminosities are overestimated by an average factor of $\sim$5 when SED fitting with 24~$\micron$ photometry alone. This result arises partly due to the fact that high redshift galaxies exhibit aromatic feature equivalent widths that are large compared to local galaxies of similar luminosities. Through a spectral decomposition of mid-infrared spectroscopic data, we are able to isolate the fraction of IR luminosity arising from an AGN as opposed to star formation activity. This fraction is only able to account for $\sim$30\% of the total IR luminosity among the entire sample.

preprint2010arXiv

Enhanced dust heating in the bulges of early-type spiral galaxies

Stellar density and bar strength should affect the temperatures of the cool (T ~ 20-30 K) dust component in the inner regions of galaxies, which implies that the ratio of temperatures in the circumnuclear regions to the disk should depend on Hubble type. We investigate the differences between cool dust temperatures in the central 3 kpc and disk of 13 nearby galaxies by fitting models to measurements between 70 and 500 microns. We attempt to quantify temperature trends in nearby disk galaxies, with archival data from Spitzer/MIPS and new observations with Herschel/SPIRE, which were acquired during the first phases of the Herschel observations for the KINGFISH (key insights in nearby galaxies: a far-infrared survey with Herschel) sample. We fit single-temperature modified blackbodies to far-infrared and submillimeter measurements of the central and disk regions of galaxies to determine the temperature of the component(s) emitting at those wavelengths. We present the ratio of central-region-to-disk-temperatures of the cool dust component of 13 nearby galaxies as a function of morphological type. We find a significant temperature gradient in the cool dust component in all galaxies, with a mean center-to-disk temperature ratio of 1.15 +/- 0.03. The cool dust temperatures in the central ~3 kpc of nearby galaxies are 23(+/-3)% hotter for morphological types earlier than Sc, and only 9(+/-3)% hotter for later types. The temperature ratio is also correlated with bar strength, with only strongly barred galaxies having a ratio over 1.2. The strong radiation field in the high stellar density of a galactic bulge tends to heat the cool dust component to higher temperatures, at least in early-type spirals with relatively large bulges, especially when paired with a strong bar.

preprint2010arXiv

Far-Infrared Line Imaging of the Starburst Ring in NGC 1097 with the Herschel/PACS Spectrometer

NGC 1097 is a nearby SBb galaxy with a Seyfert nucleus and a bright starburst ring. We study the physical properties of the interstellar medium (ISM) in the ring using spatially resolved far-infrared spectral maps of the circumnuclear starburst ring of NGC 1097, obtained with the PACS spectrometer on board the Herschel Space Telescope. In particular, we map the important ISM cooling and diagnostic emission lines of [OI] 63 $μ$m, [OIII] 88 $μ$m, [NII] 122 $μ$m, [CII] 158 $μ$m and [NII] 205 $μ$m. We observe that in the [OI] 63 $μ$m, [OIII] 88 $μ$m, and [NII] 122 $μ$m line maps, the emission is enhanced in clumps along the NE part of the ring. We observe evidence of rapid rotation in the circumnuclear ring, with a rotation velocity of ~220$ km s$^{-1}$ (inclination uncorrected) measured in all lines. The [OI] 63 $μ$m/[CII] 158 $μ$m ratio varies smoothly throughout the central region, and is enhanced on the northeastern part of the ring, which may indicate a stronger radiation field. This enhancement coincides with peaks in the [OI] 63 $μ$m and [OIII] 88 $μ$m maps. Variations of the [NII] 122 $μ$m/[NII] 205 $μ$m ratio correspond to a range in the ionized gas density between 150 and 400 cm$^{-3}$.

preprint2010arXiv

Global cosmic-ray related luminosity and energy budget of the Milky Way

We use the GALPROP code for cosmic-ray (CR) propagation to calculate the broad-band luminosity spectrum of the Milky Way related to CR propagation and interactions in the interstellar medium. This includes gamma-ray emission from the production and subsequent decay of neutral pions, bremsstrahlung, and inverse Compton scattering, and synchrotron radiation. The Galaxy is found to be nearly a CR electron calorimeter, but {\it only} if gamma ray emitting processes are taken into account. Synchrotron radiation alone accounts for only one third of the total electron energy losses with ~10-20% of the total synchrotron emission from secondary CR electrons and positrons. The relationship between far-infrared and radio luminosity that we find from our models is consistent with that found for galaxies in general. The results will be useful for understanding the connection between diffuse emissions from radio through gamma rays in ``normal'' (non-AGN dominated) galaxies, as well as for estimating the broad-band extragalactic diffuse background from these kinds of galaxies.

preprint2010arXiv

Mapping far-IR emission from the central kiloparsec of NGC 1097

Using photometry of NGC 1097 from the Herschel PACS (Photodetector Array Camera and Spectrometer) instrument, we study the resolved properties of thermal dust continuum emission from a circumnuclear starburst ring with a radius ~ 900 pc. These observations are the first to resolve the structure of a circumnuclear ring at wavelengths that probe the peak (i.e. lambda ~ 100 micron) of the dust spectral energy distribution. The ring dominates the far-infrared (far-IR) emission from the galaxy - the high angular resolution of PACS allows us to isolate the ring's contribution and we find it is responsible for 75, 60 and 55% of the total flux of NGC 1097 at 70, 100 and 160 micron, respectively. We compare the far-IR structure of the ring to what is seen at other wavelengths and identify a sequence of far-IR bright knots that correspond to those seen in radio and mid-IR images. The mid- and far-IR band ratios in the ring vary by less than +/- 20% azimuthally, indicating modest variation in the radiation field heating the dust on ~ 600 pc scales. We explore various explanations for the azimuthal uniformity in the far-IR colors of the ring including a lack of well-defined age gradients in the young stellar cluster population, a dominant contribution to the far-IR emission from dust heated by older (> 10 Myr) stars and/or a quick smoothing of local enhancements in dust temperature due to the short orbital period of the ring. Finally, we improve previous limits on the far-IR flux from the inner ~ 600 pc of NGC 1097 by an order of magnitude, providing a better estimate of the total bolometric emission arising from the active galactic nucleus and its associated central starburst.

preprint2010arXiv

Nearby Galaxies with Spitzer

We review the main advances brought by the Spitzer Space Telescope in the field of nearby galaxies studies, concentrating on a few subject areas, including: (1) the physics of the Polycyclic Aromatic Hydrocarbons that generate the mid-infrared features between ~3.5 micron and ~20 micron; (2) the use of the mid- and far-infrared emission from galaxies as star formation rate indicators; and (3) the improvement of mid-infrared diagnostics to discriminate between thermal (star-formation) and non-thermal (AGN) emission in galaxies and galaxy centers.

preprint2010arXiv

The Buried Starburst in the Interacting Galaxy II Zw 096 as Revealed by the Spitzer Space Telescope

An analysis of data from the Spitzer Space Telescope, Hubble Space Telescope, Chandra X-ray Observatory, and AKARI Infrared Astronomy Satellite is presented for the z=0.036 merging galaxy system II Zw 096 (CGCG 448-020). Because II Zw 096 has an infrared luminosity of log(L_IR/L_sun) = 11.94, it is classified as a Luminous Infrared Galaxy (LIRG), and was observed as part of the Great Observatories All-sky LIRG Survey (GOALS). The Spitzer data suggest that 80% of the total infrared luminosity comes from an extremely compact, red source not associated with the nuclei of the merging galaxies. The Spitzer mid-infrared spectra indicate no high-ionization lines from a buried active galactic nucleus in this source. The strong detection of the 3.3 micron and 6.2 micron PAH emission features in the AKARI and Spitzer spectra also implies that the energy source of II Zw 096 is a starburst. Based on Spitzer infrared imaging and AKARI near-infrared spectroscopy, the star formation rate is estimated to be 120 M_sun/yr and > 45 M_sun/yr, respectively. Finally, the high-resolution B, I, and H-band images show many star clusters in the interacting system. The colors of these clusters suggest at least two populations - one with an age of 1-5 Myr and one with an age of 20-500 Myr, reddened by 0-2 magnitudes of visual extinction. The masses of these clusters span a range between 10^6-10^8 M_sun. This starburst source is reminiscent of the extra-nuclear starburst seen in NGC 4038/9 (the Antennae Galaxies) and Arp 299 but approximately an order of magnitude more luminous than the Antennae. The source is remarkable in that the off-nuclear infrared luminosity dominates the enitre system.

preprint2010arXiv

The Calibration of Monochromatic Far-Infrared Star Formation Rate Indicators

(Abridged) Spitzer data at 24, 70, and 160 micron and ground-based H-alpha images are analyzed for a sample of 189 nearby star-forming and starburst galaxies to investigate whether reliable star formation rate (SFR) indicators can be defined using the monochromatic infrared dust emission centered at 70 and 160 micron. We compare recently published recipes for SFR measures using combinations of the 24 micron and observed H-alpha luminosities with those using 24 micron luminosity alone. From these comparisons, we derive a reference SFR indicator for use in our analysis. Linear correlations between SFR and the 70 and 160 micron luminosity are found for L(70)>=1.4x10^{42} erg/s and L(160)>=2x10^{42} erg/s, corresponding to SFR>=0.1-0.3 M_sun/yr. Below those two luminosity limits, the relation between SFR and 70 micron (160 micron) luminosity is non-linear and SFR calibrations become problematic. The dispersion of the data around the mean trend increases for increasing wavelength, becoming about 25% (factor ~2) larger at 70 (160) micron than at 24 micron. The increasing dispersion is likely an effect of the increasing contribution to the infrared emission of dust heated by stellar populations not associated with the current star formation. The non-linear relation between SFR and the 70 and 160 micron emission at faint galaxy luminosities suggests that the increasing transparency of the interstellar medium, decreasing effective dust temperature, and decreasing filling factor of star forming regions across the galaxy become important factors for decreasing luminosity. The SFR calibrations are provided for galaxies with oxygen abundance 12+Log(O/H)>8.1. At lower metallicity the infrared luminosity no longer reliably traces the SFR because galaxies are less dusty and more transparent.

preprint2010arXiv

The Detection of a Population of Submillimeter-Bright, Strongly-Lensed Galaxies

Gravitational lensing is a powerful astrophysical and cosmological probe and is particularly valuable at submillimeter wavelengths for the study of the statistical and individual properties of dusty starforming galaxies. However the identification of gravitational lenses is often time-intensive, involving the sifting of large volumes of imaging or spectroscopic data to find few candidates. We used early data from the Herschel Astrophysical Terahertz Large Area Survey to demonstrate that wide-area submillimeter surveys can simply and easily detect strong gravitational lensing events, with close to 100% efficiency.

preprint2010arXiv

The Nuclear Structure in Nearby Luminous Infrared Galaxies: HST NICMOS Imaging of the GOALS Sample

We present results of Hubble Space Telescope NICMOS H-band imaging of 73 of most luminous (i.e., log[L_IR/L_0]>11.4) Infrared Galaxies (LIRGs) in the Great Observatories All-sky LIRG Survey (GOALS). This dataset combines multi-wavelength imaging and spectroscopic data from space (Spitzer, HST, GALEX, and Chandra) and ground-based telescopes. In this paper we use the high-resolution near-infrared data to recover nuclear structure that is obscured by dust at optical wavelengths and measure the evolution in this structure along the merger sequence. A large fraction of all galaxies in our sample possess double nuclei (~63%) or show evidence for triple nuclei (~6%). Half of these double nuclei are not visible in the HST B-band images due to dust obscuration. The majority of interacting LIRGs have remaining merger timescales of 0.3 to 1.3 Gyrs, based on the projected nuclear separations and the mass ratio of nuclei. We find that the bulge luminosity surface density increases significantly along the merger sequence (primarily due to a decrease of the bulge radius), while the bulge luminosity shows a small increase towards late merger stages. No significant increase of the bulge Sersic index is found. LIRGs that show no interaction features have on average a significantly larger bulge luminosity, suggesting that non merging LIRGs have larger bulge masses than merging LIRGs. This may be related to the flux limited nature of the sample and the fact that mergers can significantly boost the IR luminosity of otherwise low luminosity galaxies. We find that the projected nuclear separation is significantly smaller for ULIRGs (median value of 1.2 kpc) than for LIRGs (mean value of 6.7 kpc), suggesting that the LIRG phase appears earlier in mergers than the ULIRG phase.

preprint2010arXiv

The Spatial Extent of (U)LIRGs in the mid-Infrared I: The Continuum Emission

We present an analysis of the extended mid-infrared (MIR) emission of the Great Observatories All-Sky LIRG Survey (GOALS) sample based on 5-15um low resolution spectra obtained with the IRS on Spitzer. We calculate the fraction of extended emission as a function of wavelength for the galaxies in the sample, FEE_lambda. We can identify 3 general types of FEE_lambda: one where it is constant, one where features due to emission lines and PAHs appear more extended than the continuum, and a third which is characteristic of sources with deep silicate absorption at 9.7um. More than 30% of the galaxies have a median FEE_lambda larger than 0.5 implying that at least half of their MIR emission is extended. Luminous Infrared Galaxies (LIRGs) display a wide range of FEE in their warm dust continuum (0<=FEE_13.2um<=0.85). The large values of FEE_13.2um that we find in many LIRGs suggest that their extended MIR continuum emission originates in scales up to 10kpc. The mean size of the LIRG cores at 13.2um is 2.6kpc. However, once the LIR of the systems reaches the threshold of ~10^11.8Lsun, all sources become clearly more compact, with FEE_13.2um<=0.2, and their cores are unresolved. Our estimated upper limit for the core size of ULIRGs is less than 1.5kpc. The analysis indicates that the compactness of systems with LIR>~10^11.25Lsun strongly increases in those classified as mergers in their final stage of interaction. The FEE_13.2um is also related to the contribution of an active galactic nucleus (AGN) to the MIR. Galaxies which are more AGN-dominated are less extended, independently of their LIR. We finally find that the extent of the MIR continuum emission is correlated with the far-IR IRAS log(f_60um/f_100um) color. This enables us to place a lower limit to the area in a galaxy from where the cold dust emission may originate, a prediction which can be tested soon with the Herschel Space Telescope.

preprint2009arXiv

Balancing the Energy Budget between Star-Formation and AGN in High Redshift Infrared Luminous Galaxies

(abridged) We present deep Spitzer mid-infrared spectroscopy, along with 16, 24, 70, and 850um photometry, for 22 galaxies located in GOODS-N. The sample spans a redshift range of 0.6 < z < 2.6, 24um flux densities between ~0.2-1.2 mJy, and consists of SMGs, AGN, and optically faint (z_AB > 25) sources. We find that IR luminosities derived by fitting local SEDs with 24um photometry alone are well matched to those when additional mid-infrared spectroscopic and longer wavelength photometric data is used for galaxies having z < 1.4 and 24um-derived IR luminosities typically > 3x10^12 L_sun. However, for galaxies in the redshift range between 1.4 < z < 2.6, typically having 24um-derived IR luminosities > 3x10^12 L_sun, IR luminosities are overestimated by an average factor of ~5 when SED fitting with 24um photometry alone. This result arises partly due to the fact that high redshift galaxies exhibit PAH EQWs that are large compared to local galaxies of similar luminosities. Using improved estimates for the IR luminosities of these sources, we investigate whether their IR emission is found to be in excess relative to that expected based on extinction corrected UV SFRs, possibly suggesting the presence of an obscured AGN. Through a spectral decomposition of IRS spectroscopic data, we are able to isolate the fraction of IR luminosity arising from an AGN as opposed to star formation activity. This fraction is only ~30% of the total IR luminosity among the entire sample, on average. Of the sources identified as having mid-infrared excesses, half are accounted for by using proper bolometric corrections while half show the presence of obscured AGN. We do not find evidence for evolution in the FIR-radio correlation over this redshift range, although the SMGs have IR/radio ratios which are, on average, ~3 times lower than the nominal value.

preprint2009arXiv

GOALS: The Great Observatories All-Sky LIRG Survey

The Great Observatories All-sky LIRG Survey (GOALS) combines data from NASA's Spitzer, Chandra, Hubble and GALEX observatories, together with ground-based data into a comprehensive imaging and spectroscopic survey of over 200 low redshift Luminous Infrared Galaxies (LIRGs). The LIRGs are a complete subset of the IRAS Revised Bright Galaxy Sample (RBGS). The LIRGs targeted in GOALS span the full range of nuclear spectral types defined via traditional optical line-ratio diagrams as well as interaction stages. They provide an unbiased picture of the processes responsible for enhanced infrared emission in galaxies in the local Universe. As an example of the analytic power of the multi-wavelength GOALS dataset, we present data for the interacting system VV 340 (IRAS F14547+2449). Between 80-95% of the total far-infrared emission (or about 5E11 solar luminosities) originates in VV 340 North. While the IRAC colors of VV 340 North and South are consistent with star-forming galaxies, both the Spitzer IRS and Chandra ACIS data indicate the presence of a buried AGN in VV 340 North. The GALEX far and near-UV fluxes imply a extremely large infrared "excess" (IRX) for the system (IR/FUV = 81) which is well above the correlation seen in starburst galaxies. Most of this excess is driven by VV 340 N, which alone has an IR excess of nearly 400. The VV 340 system seems to be comprised of two very different galaxies - an infrared luminous edge-on galaxy (VV 340 North) that dominates the long-wavelength emission from the system and which hosts a buried AGN, and a face-on starburst (VV 340 South) that dominates the short-wavelength emission.

preprint2009arXiv

The Detection of Anomalous Dust Emission in the Nearby Galaxy NGC 6946

We report on the Ka-band (26-40 GHz) emission properties for 10 star-forming regions in the nearby galaxy NGC 6946. From a radio spectral decomposition, we find that the 33 GHz flux densities are typically dominated by thermal (free-free) radiation. However, we also detect excess Ka-band emission for an outer-disk star-forming region relative to what is expected given existing radio, submillimeter, and infrared data. Among the 10 targeted regions, measurable excess emission at 33 GHz is detected for half of them, but in only one region is the excess found to be statistically significant ($\approx7σ$). We interpret this as the first likely detection of so called `anomalous' dust emission outside of the Milky Way. We find that models explaining this feature as the result of dipole emission from rapidly rotating ultrasmall grains are able to reproduce the observations for reasonable interstellar medium conditions. While these results suggest that the use of Ka-band data as a measure of star formation activity in external galaxies may be complicated by the presence of anomalous dust, it is unclear how significant a factor this will be for globally integrated measurements as the excess emission accounts for $\la$10% of the total Ka-band flux density from all 10 regions.

preprint2008arXiv

Environmental Effects in Clusters: Modified Far-Infrared--Radio Relations within Virgo Cluster Galaxies

(abridged) We present a study on the effects of the intracluster medium (ICM) on the interstellar medium (ISM) of 10 Virgo cluster spiral galaxies using {\it Spitzer} far-infrared (FIR) and VLA radio continuum imaging. Relying on the FIR-radio correlation within normal galaxies, we use our infrared data to create model radio maps which we compare to the observed radio images. For 6 of our sample galaxies we find regions along their outer edges that are highly deficient in the radio compared with our models. We believe these observations are the signatures of ICM ram pressure. For NGC 4522 we find the radio deficit region to lie just exterior to a region of high radio polarization and flat radio spectral index, although the total 20 cm radio continuum in this region does not appear strongly enhanced. These characteristics seem consistent for other galaxies with radio polarization data in the literature. The strength of the radio deficit is inversely correlated with the time since peak pressure as inferred from stellar population studies and gas stripping simulations, suggesting the strength of the radio deficit is good indicator of the strength of the current ram pressure. We also find that galaxies having {\it local} radio {\it deficits} appear to have {\it enhanced global} radio fluxes. Our preferred physical picture is that the observed radio deficit regions arise from the ICM wind sweeping away cosmic-ray (CR) electrons and the associated magnetic field, thereby creating synchrotron tails as observed for some of our galaxies. We propose that CR particles are also re-accelerated by ICM-driven shocklets behind the observed radio deficit regions which in turn enhances the remaining radio disk brightness.

preprint2008arXiv

The Spitzer Infrared Nearby Galaxies Survey: A High-Resolution Spectroscopy Anthology

High resolution mid-infrared spectra are presented for 155 nuclear and extranuclear regions from the Spitzer Infrared Nearby Galaxies Survey (SINGS). The fluxes for nine atomic forbidden and three molecular hydrogen mid-infrared emission lines are also provided, along with upper limits in key lines for infrared-faint targets. The SINGS sample shows a wide range in the ratio of [SIII]18.71um/[SIII]33.48um, but the average ratio of the ensemble indicates a typical interstellar electron density of 300-400 cm^{-3} on ~23"x15" scales and 500-600 cm^{-3} using ~11"x9" apertures, independent of whether the region probed is a star-forming nuclear, a star-forming extranuclear, or an AGN environment. Evidence is provided that variations in gas-phase metallicity play an important role in driving variations in radiation field hardness, as indicated by [NeIII]15.56um/[NeII]12.81um, for regions powered by star formation. Conversely, the radiation hardness for galaxy nuclei powered by accretion around a massive black hole is independent of metal abundance. Furthermore, for metal-rich environments AGN are distinguishable from star-forming regions by significantly larger [NeIII]15.56um/[NeII]12.81um ratios. Finally, [FeII]25.99um/[NeII]12.81um versus [SiII]34.82um/[SIII]33.48um also provides an empirical method for discerning AGN from normal star-forming sources. However, similar to [NeIII]15.56um/[NeII]12.81um, these mid-infrared line ratios lose their AGN/star-formation diagnostic powers for very low metallicity star-forming systems with hard radiation fields.

preprint2007arXiv

The Calibration of Mid-Infrared Star Formation Rate Indicators

With the goal of investigating the degree to which the mid-infrared emission traces the star formation rate (SFR), we analyze Spitzer 8 um and 24 um data of star-forming regions in a sample of 33 nearby galaxies with available HST/NICMOS images in the Paschen-alpha (1.8756 um) emission line. The galaxies are drawn from the Spitzer Infrared Nearby Galaxies Survey (SINGS) sample, and cover a range of morphologies and a factor ~10 in oxygen abundance. Published data on local low-metallicity starburst galaxies and Luminous Infrared Galaxies are also included in the analysis. Both the stellar-continuum-subtracted 8 um emission and the 24 um emission correlate with the extinction-corrected Pa-alpha line emission, although neither relationship is linear. Simple models of stellar populations and dust extinction and emission are able to reproduce the observed non-linear trend of the 24 um emission versus number of ionizing photons, including the modest deficiency of 24 um emission in the low metallicity regions, which results from a combination of decreasing dust opacity and dust temperature at low luminosities. Conversely, the trend of the 8 um emission as a function of the number of ionizing photons is not well reproduced by the same models. The 8 um emission is contributed, in larger measure than the 24 um emission, by dust heated by non-ionizing stellar populations, in agreement with previous findings. Two SFR calibrations, one using the 24 um emission and the other using a combination of the 24 um and H-alpha luminosities (Kennicutt et al. 2007), are presented. No calibration is presented for the 8 um emission, because of its significant dependence on both metallicity and environment. The calibrations presented here should be directly applicable to systems dominated by on-going star formation.

preprint2007arXiv

Variations in 24 micron morphologies among galaxies in the Spitzer Infrared Nearby Galaxies Survey: New insights into the Hubble sequence

To study the distribution of star formation and dust emission within nearby galaxies, we measured five morphological parameters in the 3.6 and 24 micron wave bands for 65 galaxies in the Spitzer Infrared Nearby Galaxies Survey (SINGS) and 8 galaxies that were serendipitously observed by SINGS. The morphological parameters demonstrate strong variations along the Hubble sequence, including statistically significant differences between S0/a-Sab and Sc-Sd galaxies. Early-type galaxies are generally found to be compact, centralized, symmetric sources in the 24 micron band, while late-type galaxies are generally found to be extended, asymmetric sources. These results suggest that the processes that increase the real or apparent sizes of galaxies' bulges also lead to more centralized 24 micron dust emission. Several phenomena, such as strong nuclear star formation, Seyfert activity, or outer ring structures, may cause galaxies to deviate from the general morphological trends observed at 24 microns. We also note that the 24 micron morphologies of Sdm-Im galaxies are quite varied, with some objects appearing very compact and symmetric while others appear diffuse and asymmetric. These variations reflect the wide variation in star formation in irregular galaxies as observed at other wavelengths. The variations in the 24 micron morphological parameters across the Hubble sequence mirror many of the morphological trends seen in other tracers of the ISM and in stellar emission. However, the 24 micron morphological parameters for the galaxies in this sample do not match the morphological parameters measured in the stellar wave bands. This implies that the distribution of dust emission is related to but not equivalent to the distribution of stellar emission.

preprint2007arXiv

Warm molecular hydrogen in the Spitzer SINGS galaxy sample

(simplified) Results on the properties of warm H2 in 57 normal galaxies are derived from H2 rotational transitions, obtained as part of SINGS. This study extends previous extragalactic surveys of H2, the most abundant constituent of the molecular ISM, to more common systems (L_FIR = e7 to 6e10 L_sun) of all morphological and nuclear types. The S(1) transition is securely detected in the nuclear regions of 86% of SINGS galaxies with stellar masses above 10^9.5 M_sun. The derived column densities of warm H2 (T > ~100 K), even though averaged over kiloparsec-scale areas, are commensurate with those of resolved PDRs; the median of the sample is 3e20 cm-2. They amount to between 1% and >30% of the total H2. The power emitted in the sum of the S(0) to S(2) transitions is on average 30% of the [SiII] line power, and ~4e-4 of the total infrared power (TIR) within the same area for star-forming galaxies, which is consistent with excitation in PDRs. The fact that H2 emission scales tightly with PAH emission, even though the average radiation field intensity varies by a factor ten, can also be understood if both tracers originate predominantly in PDRs, either dense or diffuse. A large fraction of the 25 LINER/Sy targets, however, strongly depart from the rest of the sample, in having warmer H2 in the excited states, and an excess of H2 emission with respect to PAHs, TIR and [SiII]. We propose a threshold in H2 to PAH power ratios, allowing the identification of low-luminosity AGNs by an excess H2 excitation. A dominant contribution from shock heating is favored in these objects. Finally, we detect, in nearly half the star-forming targets, non-equilibrium ortho to para ratios, consistent with FUV pumping combined with incomplete ortho-para thermalization by collisions, or possibly non-equilibrium PDR fronts advancing into cold gas.

preprint2006arXiv

Spitzer and JCMT Observations of the Active Galactic Nucleus in the Sombrero Galaxy (NGC 4594)

We present Spitzer 3.6-160 micron images, Spitzer mid-infrared spectra, and JCMT SCUBA 850 micron images of the Sombrero Galaxy (NGC 4594), an Sa galaxy with a 10^9 M_solar low luminosity active galactic nucleus (AGN). The brightest infrared sources in the galaxy are the nucleus and the dust ring. The spectral energy distribution of the AGN demonstrates that, while the environment around the AGN is a prominent source of mid-infrared emission, it is a relatively weak source of far-infrared emission, as had been inferred for AGN in previous research. The weak nuclear 160 micron emission and the negligible polycyclic aromatic hydrocarbon emission from the nucleus also implies that the nucleus is a site of only weak star formation activity and the nucleus contains relatively little cool interstellar gas needed to fuel such activity. We propose that this galaxy may be representative of a subset of low ionization nuclear emission region galaxies that are in a quiescent AGN phase because of the lack of gas needed to fuel circumnuclear star formation and Seyfert-like AGN activity. Surprisingly, the AGN is the predominant source of 850 micron emission. We examine the possible emission mechanisms that could give rise to the 850 micron emission and find that neither thermal dust emission, CO line emission, bremsstrahlung emission, nor the synchrotron emission observed at radio wavelengths can adequately explain the measured 850 micron flux density by themselves. The remaining possibilities for the source of the 850 micron emission include a combination of known emission mechanisms, synchrotron emission that is self-absorbed at wavelengths longer than 850 microns, or unidentified spectral lines in the 850 micron band.

preprint2006arXiv

The Effect of Star Formation on the Far-Infrared--Radio Correlation within Galaxies

Using data obtained for twelve galaxies as part of the {\it Spitzer} Infrared Nearby Galaxies Survey (SINGS) and the Westerbork Synthesis Radio Telescope (WSRT)-SINGS radio continuum survey, we study how star formation activity affects the far-infrared (FIR)--radio correlation {\it within} galaxies by testing a phenomenological model, which describes the radio image as a smeared version of the FIR image. The physical basis of this description is that cosmic-ray (CR) electrons will diffuse measurably farther than the mean free path of dust-heating photons before decaying by synchrotron radiation. This description works well in general. Galaxies with higher infrared surface brightnesses have best-fit smoothing scale-lengths of a few hundred parsecs, substantially shorter than those for lower surface brightness galaxies. We interpret this result to suggest that galaxies with higher disk averaged star formation rates have had a recent episode of enhanced star formation and are characterized by a higher fraction of young CR electrons that have traveled only a few hundred parsecs from their acceleration sites in supernova remnants compared to galaxies with lower star formation activity.

preprint2006arXiv

The Mid-Infrared Spectrum of Star-Forming Galaxies: Global Properties of PAH Emission

We present a sample of low-resolution 5-38um Spitzer IRS spectra of the inner few square kiloparsecs of 59 nearby galaxies spanning a large range of star formation properties. A robust method for decomposing mid-infrared galaxy spectra is described, and used to explore the behavior of PAH emission and the prevalence of silicate dust extinction. Evidence for silicate extinction is found in ~1/8 of the sample, at strengths which indicate most normal galaxies undergo A_V < ~3 magnitudes averaged over their centers. The contribution of PAH emission to the total infrared power is found to peak near 10% and extend up to ~20%, and is suppressed at metallicities Z < ~Z_sun/4, as well as in low-luminosity AGN environments. Strong inter-band PAH feature strength variations (2-5x) are observed, with the presence of a weak AGN and, to a lesser degree, increasing metallicity shifting power to the longer wavelength bands. A peculiar PAH emission spectrum with markedly diminished 5-8um features arises among the sample solely in systems with relatively hard radiation fields harboring low-luminosity AGN. The AGN may modify the emitting grain distribution and provide the direct excitation source of the unusual PAH emission, which cautions against using absolute PAH strength to estimate star formation rates in systems harboring active nuclei. Alternatively, the low star formation intensity often associated with weak AGN may affect the spectrum. The effect of variations in the mid-infrared spectrum on broadband infrared surveys is modeled, and points to more than a factor of two uncertainty in results which assume a fixed PAH emission spectrum, for redshifts z=0-2.5.

preprint2005arXiv

The Discovery of the Optical and Near-IR Afterglows of the First Swift Gamma-Ray Bursts

We present optical and near-infrared searches for afterglow emission from the first four Swift bursts with accurate positions from the X-ray Telescope (XRT). Using telescopes at Las Campanas, Keck, and Palomar observatories we rapidly identified and followed up afterglows for three of the four bursts. The burst positions were also observed with the Very Large Array, but no radio afterglow emission was detected. The optical/NIR afterglows are fainter than about 75% of all afterglows detected to date, with GRB 050126 being the faintest, and were identified thanks to accurate and rapid positions from the XRT and rapid response with >1-m telescopes. This suggests that the fraction of dust-obscured bursts is small, <10% when combined with afterglows localized by the HETE-2 Soft X-ray Camera. The X-ray fluxes are typical of the known population, with the exception of GRB 050126 which has the faintest X-ray afterglow to date (normalized to t=10 hr), and was detected thanks to a response time of only 130 s after the burst. Finally, we find that all three optical/NIR afterglows are located <2 arcsec away from the nominal XRT positions, suggesting that the XRT is capable of delivering highly accurate positions, which will revolutionize afterglow studies.

preprint2004arXiv

Mid-Infrared IRS Spectroscopy of NGC 7331: A First Look at the SINGS Legacy

The nearby spiral galaxy NGC 7331 was spectrally mapped from 5-38um using all modules of Spitzer's IRS spectrograph. A strong new dust emission feature, presumed due to PAHs, was discovered at 17.1um. The feature's intensity is nearly half that of the ubiquitous 11.3um band. The 7-14um spectral maps revealed significant variation in the 7.7 and 11.3um PAH features between the stellar ring and nucleus. Weak [OIV] 25.9um line emission was found to be centrally concentrated in the nucleus, with an observed strength over 10% of the combined neon line flux, indicating an AGN or unusually active massive star photo-ionization. Two [SIII] lines fix the characteristic electron density in the HII regions at n_e < ~200 cm^-3. Three detected H_2 rotational lines, tracing warm molecular gas, together with the observed IR continuum, are difficult to match with standard PDR models. Either additional PDR heating or shocks are required to simultaneously match lines and continuum.

preprint2004arXiv

Spatially Resolved Ultraviolet, H-alpha, Infrared, and Radio Star Formation in M81

We present Multiband Imaging Photometer for Spitzer (MIPS) observations of M81 at 24, 70, and 160 micron. The grand design nature of M81 is clearly seen, showing two well resolved spiral arms containing numerous bright star forming regions. The MIPS images reveal a significant amount of cold dust associated with the spiral arms. We investigate the variation of the ultraviolet (UV), H-alpha, and infrared (IR) luminosities and star formation rate (SFR) indicators across the face of M81 using the MIPS images and archival UV and H-alpha images. For regions in M81, we find that UV and H-alpha SFRs (uncorrected for dust attenuation) are always lower than the IR SFR. The cause of this behavior is dust attenuation and/or using SFR calibrations appropriate for entire galaxies, not regions in galaxies. The characteristics of the dust attenuation for the regions indicate the dust grains and/or geometry are different from those in starburst galaxies. The behavior of the infrared-radio correlation in M81 is seen to vary from the global average, with variations correlated with the morphology of M81.