Source author record

L. Armus

L. Armus 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

101works
5topics
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

101 published item(s)

preprint2022arXiv

Preparing for low surface brightness science with the Vera C. Rubin Observatory: characterisation of tidal features from mock images

Tidal features in the outskirts of galaxies yield unique information about their past interactions and are a key prediction of the hierarchical structure formation paradigm. The Vera C. Rubin Observatory is poised to deliver deep observations for potentially of millions of objects with visible tidal features, but the inference of galaxy interaction histories from such features is not straightforward. Utilising automated techniques and human visual classification in conjunction with realistic mock images produced using the NEWHORIZON cosmological simulation, we investigate the nature, frequency and visibility of tidal features and debris across a range of environments and stellar masses. In our simulated sample, around 80 per cent of the flux in the tidal features around Milky Way or greater mass galaxies is detected at the 10-year depth of the Legacy Survey of Space and Time (30-31 mag / sq. arcsec), falling to 60 per cent assuming a shallower final depth of 29.5 mag / sq. arcsec. The fraction of total flux found in tidal features increases towards higher masses, rising to 10 per cent for the most massive objects in our sample (M*~10^{11.5} Msun). When observed at sufficient depth, such objects frequently exhibit many distinct tidal features with complex shapes. The interpretation and characterisation of such features varies significantly with image depth and object orientation, introducing significant biases in their classification. Assuming the data reduction pipeline is properly optimised, we expect the Rubin Observatory to be capable of recovering much of the flux found in the outskirts of Milky Way mass galaxies, even at intermediate redshifts (z<0.2).

preprint2020arXiv

A Hard X-ray Test of HCN Enhancements as a Tracer of Embedded Black Hole Growth

Enhanced emission from the dense gas tracer HCN (relative to HCO$^+ $) has been proposed as a signature of active galactic nuclei (AGN). In a previous single-dish millimeter line survey we identified galaxies with HCN/HCO$ ^+ $ (1-0) intensity ratios consistent with those of many AGN but whose mid-infrared spectral diagnostics are consistent with little to no ( $\lesssim15\% $) contribution of an AGN to the bolometric luminosity. To search for putative heavily obscured AGN, we present and analyze \nustar hard X-ray (3-79 keV) observations of four such galaxies from the Great Observatories All-sky LIRG Survey. We find no X-ray evidence for AGN in three of the systems and place strong upper limits on the energetic contribution of any heavily obscured ($N_{\rm H}>10^{24}$ cm$^{-2}$) AGN to their bolometric luminosity. The X-ray flux upper limits are presently an order of magnitude below what XDR-driven chemistry model predict are necessary to drive HCN enhancements. In a fourth system we find a hard X-ray excess consistent with the presence of an AGN, but contributing only $\sim3\%$ of the bolometric luminosity. It is also unclear if the AGN is spatially associated with the HCN enhancement. We further explore the relationship between HCN/HCO$^+$ (for several $\mathrm{J}_\mathrm{upper}$ levels) and $L_\mathrm{AGN}/L_\mathrm{IR}$ for a larger sample of systems in the literature. We find no evidence for correlations between the line ratios and the AGN fraction derived from X-rays, indicating that HCN/HCO$^+$ intensity ratios are not driven by the energetic dominance of AGN, nor are they reliable indicators of whether SMBH accretion is ongoing.

preprint2020arXiv

The Molecular Gas in the NGC 6240 Merging Galaxy System at the Highest Spatial Resolution

We present the highest resolution --- 15 pc (0.03'') --- ALMA $^{12}$CO(2-1) line emission and 1.3mm continuum maps, tracers of the molecular gas and dust, respectively, in the nearby merging galaxy system NGC 6240, that hosts two supermassive black holes growing simultaneously. These observations provide an excellent spatial match to existing Hubble optical and near-infrared observations of this system. A significant molecular gas mass, $\sim$9$\times$10$^9$M$_\odot$, is located in between the two nuclei, forming a clumpy stream kinematically dominated by turbulence, rather than a smooth rotating disk as previously assumed from lower resolution data. Evidence for rotation is seen in the gas surrounding the southern nucleus, but not in the northern one. Dynamical shells can be seen, likely associated with nuclear supernovae remnants. We further detect the presence of significant high velocity outflows, some of them reaching velocities $>$500 km/s, affecting a significant fraction, $\sim$11\% of the molecular gas in the nuclear region. Inside the spheres of influence of the northern and southern supermassive black holes we find molecular masses of 7.4$\times$10$^8$M$_\odot$ and 3.3$\times$10$^9$M$_\odot$, respectively. We are thus directly imaging the reservoir of gas that can accrete onto each supermassive black hole. These new ALMA maps highlight the critical need for high resolution observations of molecular gas in order to understand the feeding of supermassive black holes and its connection to galaxy evolution in the context of a major galaxy merger.

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.

preprint2016arXiv

Comparing [CII], HI, and CO dynamics of nearby galaxies

The HI and CO components of the interstellar medium (ISM) are usually used to derive the dynamical mass M_dyn of nearby galaxies. Both components become too faint to be used as a tracer in observations of high-redshift galaxies. In those cases, the 158 $μ$m line of atomic carbon [CII] may be the only way to derive M_dyn. As the distribution and kinematics of the ISM tracer affects the determination of M_dyn, it is important to quantify the relative distributions of HI, CO and [CII]. HI and CO are well-characterised observationally, however, for [CII] only very few measurements exist. Here we compare observations of CO, HI, and [CII] emission of a sample of nearby galaxies, drawn from the HERACLES, THINGS and KINGFISH surveys. We find that within R_25, the average [CII] exponential radial profile is slightly shallower than that of the CO, but much steeper than the HI distribution. This is also reflected in the integrated spectrum ("global profile"), where the [CII] spectrum looks more like that of the CO than that of the HI. For one galaxy, a spectrally resolved comparison of integrated spectra was possible; other comparisons were limited by the intrinsic line-widths of the galaxies and the coarse velocity resolution of the [CII] data. Using high-spectral-resolution SOFIA [CII] data of a number of star forming regions in two nearby galaxies, we find that their [CII] linewidths agree better with those of the CO than the HI. As the radial extent of a given ISM tracer is a key input in deriving M_dyn from spatially unresolved data, we conclude that the relevant length-scale to use in determining M_dyn based on [CII] data, is that of the well-characterised CO distribution. This length scale is similar to that of the optical disk.

preprint2016arXiv

Morphological Classification of Local Luminous Infrared Galaxies

We present an analysis of the morphological classification of 89 luminous infrared galaxies (LIRGs) from the Great Observatories All-sky LIRG Survey (GOALS) sample using non-parametric coefficients and compare their morphology as a function of wavelength. We rely on images obtained in the optical (B- and I-band) as well as in the infrared (H-band and 5.8$μ$m). Our classification is based on the calculation of $Gini$ and the second order of light ($M_{20}$) non-parametric coefficients which we explore as a function of stellar mass ($M_\star$), infrared luminosity ($L_{IR}$) and star formation rate (SFR). We investigate the relation between $M_{20}$, the specific SFR (sSFR) and the dust temperature ($T_{dust}$) in our galaxy sample. We find that $M_{20}$ is a better morphological tracer than $Gini$, as it allows to distinguish systems formed by double systems from isolated and post-merger LIRGs. The multi-wavelength analysis allows us to identify a region in the $Gini$-$M_{20}$ parameter space where ongoing mergers reside, regardless of the band used to calculate the coefficients. In particular when measured in the H-band, this region can be used to identify ongoing mergers, with a minimal contamination from LIRGs in other stages. We also find that while the sSFR is positively correlated with $M_{20}$ when measured in the mid-infrared, i.e. star-bursting galaxies show more compact emission, it is anti-correlated with the B-band based $M_{20}$. We interpret this as the spatial decoupling between obscured and un-obscured star formation, whereby the ultraviolet/optical size of a LIRGs experience an intense dust enshrouded central starburst is larger than in the one in the mid-infrared since the contrast between the nuclear to the extended disk emission is smaller in the mid-infrared. This has important implications for high redshift surveys of dusty sources. [abridged]

preprint2016arXiv

The Far-Infrared Surveyor Mission Study: Paper I, the Genesis

This paper describes the beginning of the Far-Infrared Surveyor mission study for NASA's Astrophysics Decadal 2020. We describe the scope of the study, and the open process approach of the Science and Technology Definition Team. We are currently developing the science cases and provide some preliminary highlights here. We note key areas for technological innovation and improvements necessary to make a Far-Infrared Surveyor mission a reality.

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

Heating and cooling of the neutral ISM in the NGC4736 circumnuclear ring

The manner in which gas accretes and orbits within circumnuclear rings has direct implications for the star formation process. In particular, gas may be compressed and shocked at the inflow points, resulting in bursts of star formation at these locations. Afterwards the gas and young stars move together through the ring. In addition, star formation may occur throughout the ring, if and when the gas reaches sufficient density to collapse under gravity. These two scenarios for star formation in rings are often referred to as the `pearls on a string' and `popcorn' paradigms. In this paper, we use new Herschel PACS observations, obtained as part of the KINGFISH Open Time Key Program, along with archival Spitzer and ground-based observations from the SINGS Legacy project, to investigate the heating and cooling of the interstellar medium in the nearby star-forming ring galaxy, NGC4736. By comparing spatially resolved estimates of the stellar FUV flux available for heating, with the gas and dust cooling derived from the FIR continuum and line emission, we show that while star formation is indeed dominant at the inflow points in NGC 4736, additional star formation is needed to balance the gas heating and cooling throughout the ring. This additional component most likely arises from the general increase in gas density in the ring over its lifetime. Our data provide strong evidence, therefore, for a combination of the two paradigms for star formation in the ring in NGC4736.

preprint2015arXiv

Probing highly obscured, self-absorbed galaxy nuclei with vibrationally excited HCN

We present high resolution (0."4) IRAM PdBI and ALMA mm and submm observations of the (ultra) luminous infrared galaxies ((U)LIRGs) IRAS17208-0014, Arp220, IC860 and Zw049.057 that reveal intense line emission from vibrationally excited ($ν_2$=1) J=3-2 and 4-3 HCN. The emission is emerging from buried, compact (r<17-70 pc) nuclei that have very high implied mid-infrared surface brightness $>$$5\times 10^{13}$ L$_{\odot}$ kpc$^{-2}$. These nuclei are likely powered by accreting supermassive black holes (SMBHs) and/or hot (>200 K) extreme starbursts. Vibrational, $ν_2$=1, lines of HCN are excited by intense 14 micron mid-infrared emission and are excellent probes of the dynamics, masses, and physical conditions of (U)LIRG nuclei when H$_2$ column densities exceed $10^{24}$ cm$^{-2}$. It is clear that these lines open up a new interesting avenue to gain access to the most obscured AGNs and starbursts. Vibrationally excited HCN acts as a proxy for the absorbed mid-infrared emission from the embedded nuclei, which allows for reconstruction of the intrinsic, hotter dust SED. In contrast, we show strong evidence that the ground vibrational state ($ν$=0), J=3-2 and 4-3 rotational lines of HCN and HCO$^+$ fail to probe the highly enshrouded, compact nuclear regions owing to strong self- and continuum absorption. The HCN and HCO$^+$ line profiles are double-peaked because of the absorption and show evidence of non-circular motions - possibly in the form of in- or outflows. Detections of vibrationally excited HCN in external galaxies are so far limited to ULIRGs and early-type spiral LIRGs, and we discuss possible causes for this. We tentatively suggest that the peak of vibrationally excited HCN emission is connected to a rapid stage of nuclear growth, before the phase of strong feedback.

preprint2015arXiv

Spatially Resolved Spitzer-IRS Spectral Maps of the Superwind in M82

We have mapped the superwind/halo region of the nearby starburst galaxy M82 in the mid-infrared with $Spitzer-IRS$. The spectral regions covered include the H$_2 S(1)-S(3)$, [NeII], [NeIII] emission lines and PAH features. We estimate the total warm H$_2$ mass and the kinetic energy of the outflowing warm molecular gas to be between $M_{warm}\sim5-17\times10^6$ M$_{\odot}$ and $E_{K}\sim6-20\times10^{53}$ erg. Using the ratios of the 6.2, 7.7 and 11.3 micron PAH features in the IRS spectra, we are able to estimate the average size and ionization state of the small grains in the superwind. There are large variations in the PAH flux ratios throughout the outflow. The 11.3/7.7 and the 6.2/7.7 PAH ratios both vary by more than a factor of five across the wind region. The Northern part of the wind has a significant population of PAH's with smaller 6.2/7.7 ratios than either the starburst disk or the Southern wind, indicating that on average, PAH emitters are larger and more ionized. The warm molecular gas to PAH flux ratios (H$_2/PAH$) are enhanced in the outflow by factors of 10-100 as compared to the starburst disk. This enhancement in the H$_2/PAH$ ratio does not seem to follow the ionization of the atomic gas (as measured with the [NeIII]/[NeII] line flux ratio) in the outflow. This suggests that much of the warm H$_2$ in the outflow is excited by shocks. The observed H$_2$ line intensities can be reproduced with low velocity shocks ($v < 40$ km s$^{-1}$) driven into moderately dense molecular gas ($10^2 <n_H < 10^4$ cm$^{-3}$) entrained in the outflow.

preprint2015arXiv

The Herschel Comprehensive (U)LIRG Emission Survey (HerCULES): CO Ladders, fine structure lines, and neutral gas cooling

(Ultra) Luminous Infrared Galaxies ((U)LIRGs) are objects characterized by their extreme infrared (8-1000 $μ$m) luminosities ($L_{LIRG}>10^{11} $L$_\odot$ and $L_{ULIRG}>10^{12}$ L$_\odot$). The Herschel Comprehensive ULIRG Emission Survey (HerCULES; PI van der Werf) presents a representative flux-limited sample of 29 (U)LIRGs that spans the full luminosity range of these objects (10$^{11}\leq L_\odot \geq10^{13}$). With the \emph{Herschel Space Observatory}, we observe [CII] 157 $μ$m, [OI] 63 $μ$m, and [OI] 145 $μ$m line emission with PACS, CO J=4-3 through J=13-12, [CI] 370 $μ$m, and [CI] 609 $μ$m with SPIRE, and low-J CO transitions with ground-based telescopes. The CO ladders of the sample are separated into three classes based on their excitation level. In 13 of the galaxies, the [OI] 63 $μ$m emission line is self absorbed. Comparing the CO excitation to the IRAS 60/100 $μ$m ratio and to far infrared luminosity, we find that the CO excitation is more correlated to the far infrared colors. We present cooling budgets for the galaxies and find fine-structure line flux deficits in the [CII], [SiII], [OI], and [CI] lines in the objects with the highest far IR fluxes, but do not observe this for CO $4\leq J_{upp}\leq13$. In order to study the heating of the molecular gas, we present a combination of three diagnostic quantities to help determine the dominant heating source. Using the CO excitation, the CO J=1-0 linewidth, and the AGN contribution, we conclude that galaxies with large CO linewidths always have high-excitation CO ladders, and often low AGN contributions, suggesting that mechanical heating is important.

preprint2015arXiv

The Spitzer Archival Far-InfraRed Extragalactic Survey

We present the Spitzer Archival Far-InfraRed Extragalactic Survey (SAFIRES). This program produces refined mosaics and source lists for all far-infrared extragalactic data taken during the more than six years of the cryogenic operation of the Spitzer Space Telescope. The SAFIRES products consist of far-infrared data in two wavelength bands (70 um and 160 um) across approximately 180 square degrees of sky, with source lists containing far-infrared fluxes for almost 40,000 extragalactic point sources. Thus, SAFIRES provides a large, robust archival far-infrared data set suitable for many scientific goals.

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

Accretion-Inhibited Star Formation in the Warm Molecular Disk of the Green-valley Elliptical Galaxy NGC 3226

We present archival Spitzer photometry and spectroscopy, and Herschel photometry, of the peculiar "Green Valley" elliptical galaxy NGC~3226. The galaxy, which contains a low-luminosity AGN, forms a pair with NGC~3227, and is shown to lie in a complex web of stellar and HI filaments. Imaging at 8 and 16$μ$m reveals a curved plume structure 3 kpc in extent, embedded within the core of the galaxy, and coincident with the termination of a 30 kpc-long HI tail. In-situ star formation associated with the IR plume is identified from narrow-band HST imaging. The end of the IR-plume coincides with a warm molecular hydrogen disk and dusty ring, containing 0.7-1.1 $\times$ 10$^7$ M$_{\odot}$ detected within the central kpc. Sensitive upper limits to the detection of cold molecular gas may indicate that a large fraction of the H$_2$ is in a warm state. Photometry, derived from the UV to the far-IR, shows evidence for a low star formation rate of $\sim$0.04 M$_{\odot}$ yr$^{-1}$ averaged over the last 100 Myrs. A mid-IR component to the Spectral Energy Distribution (SED) contributes $\sim$20$\%$ of the IR luminosity of the galaxy, and is consistent with emission associated with the AGN. The current measured star formation rate is insufficient to explain NGC3226's global UV-optical "green" colors via the resurgence of star formation in a "red and dead" galaxy. This form of "cold accretion" from a tidal stream would appear to be an inefficient way to rejuvenate early-type galaxies, and may actually inhibit star formation.

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

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

Near-Infrared long-slit spectra of Seyfert Galaxies: gas excitation across the central kiloparsec

The excitation of the gas phase of the interstellar medium can be driven by various mechanisms. In galaxies with an active nucleus, such as Seyfert galaxies, both radiative and mechanical energy from the central black hole, or the stars in the disk surrounding it may play a role. We investigate the relative importance and range of influence of the active galactic nucleus for the excitation of ionized and molecular gas in the central kiloparsec of its host galaxy. We present H- and K-band long-slit spectra for a sample of 21 nearby (D <70 Mpc) Seyfert galaxies obtained with the NIRSPEC instrument on the Keck telescope. For each galaxy, we fit the nebular line emission, stellar continua, and warm molecular gas as a function of distance from the nucleus. Our analysis does not reveal a clear difference between the nucleus proper and off-nuclear environment in terms of excitation mechanisms, suggesting that the influence of an AGN reaches far into the disk of the host galaxy. The radial variations in emission line ratios indicate that, while local mechanisms do affect the gas excitation, they are often averaged out when measuring over extended regions. The fully calibrated long-slit spectra, as well as our fitting results, are made available on-line.

preprint2014arXiv

Polycyclic aromatic hydrocarbon feature deficit of starburst galaxies in the AKARI North Ecliptic Pole Deep Field

We study the behaviour of polycyclic aromatic hydrocarbon emission in galaxies at z=0.3-1.4 using 1868 samples from the revised catalogue of AKARI North Ecliptic Pole Deep survey. The continuous filter coverage at 2-24um makes it possible to measure 8um luminosity, which is dominated by polycyclic aromatic hydrocarbon emission for galaxies at up to z=2. We compare the IR8 (= LIR/L(8)) and 8um to 4.5um luminosity ratio (L(8)/L(4.5)) with the starburstiness, Rsb, defined as excess of specific star -formation rate over that of main-sequence galaxy. All AGN candidates were excluded from our sample using an SED fitting. We found L(8)/L(4.5) increases with starburstiness at log Rsb < 0.5 and stays constant at higher starburstiness. On the other hand, IR8 is constant at log Rsb < 0, while it increases with starburstiness at log Rsb > 0. This behaviour is seen in all redshift range of our study. These results indicate that starburst galaxies have deficient polycyclic aromatic hydrocarbon emission compared with main-sequence galaxies. We also find that galaxies with extremely high L(8)/L(4.5) ratio have only moderate starburstiness. These results suggest that starburst galaxies have compact star-forming regions with intense radiation, which destroys PAHs and/or have dusty HII regions resulting in a lack of ionising photons.

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

Star formation relations and CO SLEDs across the J-ladder and redshift

We present FIR-CO luminosity relations ($\log L_{\rm FIR} = α\log L'_{\rm CO} + β$) for the full CO rotational ladder from J=1-0 to J=13-12 for 62 local (z < 0.1) (Ultra) Luminous Infrared Galaxies (LIRGs) using data from Herschel SPIRE-FTS and ground-based telescopes. We extend our sample to high redshifts (z > 1) by including 35 (sub)-millimeter selected dusty star forming galaxies from the literature with robust CO observations. The addition of luminous starbursts at high redshifts enlarge the range of the FIR-CO luminosity relations towards the high-IR-luminosity end while also significantly increasing the small amount of mid-/high-J CO line data available prior to Herschel. This new data-set (both in terms of IR luminosity and J-ladder) reveals linear FIR-CO luminosity relations ($α\sim 1$) for J=1-0 up to J=5-4, with a nearly constant normalisation ($β\sim 2$). This is expected from the (also) linear FIR-(molecular line) relations found for the dense gas tracer lines (HCN and CS), as long as the dense gas mass fraction does not vary strongly within our (merger/starburst)-dominated sample. However from J=6-5 and up to J=13-12 we find an increasingly sub-linear slope and higher normalization constant with increasing J. We argue that these are caused by a warm (~100K) and dense ($>10^4{\rm cm^{-3}}$) gas component whose thermal state is unlikely to be maintained by star formation powered far-UV radiation fields (and thus is no longer directly tied to the star formation rate). We suggest that mechanical heating (e.g., supernova driven turbulence and shocks), and not cosmic rays, is the more likely source of energy for this component. The global CO spectral line energy distributions (SLEDs), which remain highly excited from J=6-5 up to J=13-12, are found to be a generic feature of the (U)LIRGs in our sample, and further support the presence of this gas component.

preprint2014arXiv

Warm Molecular Gas in Luminous Infrared Galaxies

We present our initial results on the CO rotational spectral line energy distribution (SLED) of the $J$ to $J$$-$1 transitions from $J=4$ up to $13$ from Herschel SPIRE spectroscopic observations of 65 luminous infrared galaxies (LIRGs) in the Great Observatories All-Sky LIRG Survey (GOALS). The observed SLEDs change on average from one peaking at $J \le 4$ to a broad distribution peaking around $J \sim\,$6$-$7 as the IRAS 60-to-100 um color, $C(60/100)$, increases. However, the ratios of a CO line luminosity to the total infrared luminosity, $L_{\rm IR}$, show the smallest variation for $J$ around 6 or 7. This suggests that, for most LIRGs, ongoing star formation (SF) is also responsible for a warm gas component that emits CO lines primarily in the mid-$J$ regime ($5 \lesssim J \lesssim 10$). As a result, the logarithmic ratios of the CO line luminosity summed over CO (5$-$4), (6$-$5), (7$-$6), (8$-$7) and (10$-$9) transitions to $L_{\rm IR}$, $\log R_{\rm midCO}$, remain largely independent of $C(60/100)$, and show a mean value of $-4.13$ ($\equiv \log R^{\rm SF}_{\rm midCO}$) and a sample standard deviation of only 0.10 for the SF-dominated galaxies. Including additional galaxies from the literature, we show, albeit with small number of cases, the possibility that galaxies, which bear powerful interstellar shocks unrelated to the current SF, and galaxies, in which an energetic active galactic nucleus contributes significantly to the bolometric luminosity, have their $R_{\rm midCO}$ higher and lower than $R^{\rm SF}_{\rm midCO}$, respectively.

preprint2013arXiv

Dynamical Modeling of Galaxy Mergers using Identikit

We present dynamical models of four interacting systems: NGC 5257/8, The Mice, the Antennae, and NGC 2623. The parameter space of the encounters are constrained using the Identikit model-matching and visualization tool. Identikit utilizes hybrid N-body and test particle simulations to enable rapid exploration of the parameter space of galaxy mergers. The Identikit-derived matches of these systems are reproduced with self-consistent collisionless simulations which show very similar results. The models generally reproduce the observed morphology and \HI\ kinematics of the tidal tails in these systems with reasonable properties inferred for the progenitor galaxies. The models presented here are the first to appear in the literature for NGC 5257/8 and NGC 2623, and The Mice and the Antennae are compared with previously published models. Based on the assumed mass model and our derived initial conditions, the models indicate the four systems are currently being viewed 175-260 Myr after first passage and cover a wide range of merger stages. In some instances there are mismatches between the models and the data (e.g., in the length of a tail); these are likely due to our adoption of a single mass model for all galaxies. Despite the use of a single mass model, these results demonstrate the utility of Identikit in constraining the parameter space for galaxy mergers when applied to real data.

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

Far-Infrared Properties of Type 1 Quasars

We use the Spitzer Space Telescope Enhanced Imaging Products (SEIP) and the Spitzer Archival Far-InfraRed Extragalactic Survey (SAFIRES) to study the spectral energy distributions of spectroscopically confirmed type 1 quasars selected from the Sloan Digital Sky Survey (SDSS). By combining the Spitzer and SDSS data with the 2-Micron All Sky Survey (2MASS) we are able to construct a statistically robust rest-frame 0.1-100 micron type 1 quasar template. We find the quasar population is well-described by a single power-law SED at wavelengths less than 20 microns, in good agreement with previous work. However, at longer wavelengths we find a significant excess in infrared luminosity above an extrapolated power-law, along with signifiant object-to-object dispersion in the SED. The mean excess reaches a maximum of 0.8 dex at rest-frame wavelengths near 100 microns.

preprint2013arXiv

Hubble Space Telescope ACS Imaging of the GOALS Sample: Quantitative Structural Properties of Nearby Luminous Infrared Galaxies with L_IR > 10^{11.4} L_sun

A {\it Hubble Space Telescope} ({\it HST}) / Advanced Camera for Surveys (ACS) study of the structural properties of 85 luminous and ultraluminous ($L_{\rm IR} > 10^{11.4}$ L$_\odot$) infrared galaxies (LIRGs and ULIRGs) in the Great Observatories All-sky LIRG Survey (GOALS) sample is presented. Two-dimensional GALFIT analysis has been performed on F814W "{\it I}-band" images to decompose each galaxy, as appropriate, into bulge, disk, central PSF and stellar bar components. The fraction of bulge-less disk systems is observed to be higher in LIRGs (35%) than in ULIRGs (20%), with the disk+bulge systems making up the dominant fraction of both LIRGs (55%) and ULIRGs (45%). Further, bulge+disk systems are the dominant late-stage merger galaxy type and are the dominant type for LIRGs and ULIRGs at almost every stage of galaxy-galaxy nuclear separation. The mean {\it I}-band host absolute magnitude of the GOALS galaxies is $-22.64\pm$0.62 mag (1.8$^{+1.4}_{-0.4}$ L$^*_I$), and the mean bulge absolute magnitude in GOALS galaxies is about 1.1 magnitude fainter than the mean host magnitude. Almost all ULIRGs have bulge magnitudes at the high end (-20.6 to - 23.5 mag) of the GOALS bulge magnitude range. Mass ratios in the GOALS binary systems are consistent with most of the galaxies being the result of major mergers, and an examination of the residual-to-host intensity ratios in GOALS binary systems suggests that smaller companions suffer more tidal distortion than the larger companions. We find approximately twice as many bars in GOALS disk+bulge systems (32.8%) than in pure-disk mergers (15.9%) but most of the disk+bulge systems that contain bars are disk-dominated with small bulges. The bar-to-host intensity ratio, bar half-light radius, and bar ellipticity in GOALS galaxies are similar to those found in nearby spiral galaxies.(abridged)

preprint2013arXiv

Mid-Infrared Atomic Fine-Structure Emission Line Spectra of Luminous Infrared Galaxies: Spitzer/IRS Spectra of the GOALS Sample

We present the data and our analysis of MIR fine-structure emission lines detected in Spitzer/IRS high-res spectra of 202 local LIRGs observed as part of the GOALS project. We detect emission lines of [SIV], [NeII], [NeV], [NeIII], [SIII]18.7, [OIV], [FeII], [SIII]33.5, and [SiII]. Over 75% of our galaxies are classified as starburst (SB) sources in the MIR. We compare ratios of the emission line fluxes to stellar photo- and shock-ionization models to constrain the gas properties in the SB nuclei. Comparing the [SIV]/[NeII] and [NeIII]/[NeII] ratios to the Starburst99-Mappings III models with an instantaneous burst history, the line ratios suggest that the SB in our LIRGs have ages of 1-4.5Myr, metallicities of 1-2Z_sun, and ionization parameters of 2-8e7cm/s. Based on the [SIII]/[SIII] ratios, the electron density in LIRG nuclei has a median electron density of ~300cm-3 for sources above the low density limit. We also find that strong shocks are likely present in 10 SB sources. A significant fraction of the GOALS sources have resolved neon lines and 5 show velocity differences of >200km/s in [NeIII] or [NeV] relative to [NeII]. Furthermore, 6 SB and 5 AGN LIRGs show a trend of increasing line width with ionization potential, suggesting the possibility of a compact energy source and stratified ISM in their nuclei. We confirm a strong correlation between the [NeII]+[NeIII] emission, as well as [SIII]33.5, with both the IR luminosity and the 24um warm dust emission measured from the spectra. Finally, we find no correlation between the hardness of the radiation field or the line width and the ratio of the total IR to 8um emission (IR8). This may be because the IR luminosity and the MIR fine-structure lines are sensitive to different timescales over the SB, or that IR8 is more sensitive to the geometry of the warm dust region than the radiation field producing the HII region emission.

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

Star Formation Rates in Resolved Galaxies: Calibrations with Near and Far Infrared Data for NGC5055 and NGC6946

We use the near--infrared Brγhydrogen recombination line as a reference star formation rate (SFR) indicator to test the validity and establish the calibration of the {\it Herschel} PACS 70 μm emission as a SFR tracer for sub--galactic regions in external galaxies. Brγoffers the double advantage of directly tracing ionizing photons and of being relatively insensitive to the effects of dust attenuation. For our first experiment, we use archival CFHT Brγand Ks images of two nearby galaxies: NGC\,5055 and NGC\,6946, which are also part of the {\it Herschel} program KINGFISH (Key Insights on Nearby Galaxies: a Far-Infrared Survey with Herschel). We use the extinction corrected Brγemission to derive the SFR(70) calibration for H{\sc ii} regions in these two galaxies. A comparison of the SFR(70) calibrations at different spatial scales, from 200 pc to the size of the whole galaxy, reveals that about 50% of the total 70μm emission is due to dust heated by stellar populations that are unrelated to the current star formation. We use a simple model to qualitatively relate the increase of the SFR(70) calibration coefficient with decreasing region size to the star formation timescale. We provide a calibration for an unbiased SFR indicator that combines the observed Hαwith the 70 μm emission, also for use in H{\sc ii} regions. We briefly analyze the PACS 100 and 160 μm maps and find that longer wavelengths are not as good SFR indicators as 70μm, in agreement with previous results. We find that the calibrations show about 50% difference between the two galaxies, possibly due to effects of inclination.

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.

preprint2013arXiv

Towards A Removal of Temperature Dependencies from Abundance Determinations

The metal content of a galaxy is a key property for distinguishing between viable galaxy evolutionary scenarios, and it strongly influences many of the physical processes in the interstellar medium. An absolute and robust determination of extragalactic metallicities is essential in constraining models of chemical enrichment and chemical evolution, however, current gas phase abundance determinations from optical fine-structure lines are uncertain to 0.8 dex as conversion of these optical line fluxes to abundances is strongly dependent on the electron temperature of the ionized gas. In contrast, the far-IR emission lines can be used to derive an O++ abundance that is relatively insensitive to temperature, while the ratio of the optical to far-IR lines provides a consistent temperature to be used in the derivation of an O$^+$ abundance. We present observations of the [O III] 88 \micron\ fine-structure line in NGC 628 that were obtained as part of the KINGFISH program. These data are combined with optical IFU data to derive oxygen abundances for seven HII regions. We find the abundance of these regions to all lie between the high and low values of strong line calibrations and in agreement with estimates that assume temperature fluctuations are present in the HII regions.

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

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

Evidence for CO shock excitation in NGC 6240 from Herschel SPIRE spectroscopy

We present Herschel SPIRE FTS spectroscopy of the nearby luminous infrared galaxy NGC 6240. In total 20 lines are detected, including CO J=4-3 through J=13-12, 6 H2O rotational lines, and [CI] and [NII] fine-structure lines. The CO to continuum luminosity ratio is 10 times higher in NGC 6240 than Mrk 231. Although the CO ladders of NGC 6240 and Mrk 231 are very similar, UV and/or X-ray irradiation are unlikely to be responsible for the excitation of the gas in NGC 6240. We applied both C and J shock models to the H2 v=1-0 S(1) and v=2-1 S(1) lines and the CO rotational ladder. The CO ladder is best reproduced by a model with shock velocity v_s=10 km s^-1 and a pre-shock density n_H=5 * 10^4 cm^-3. We find that the solution best fitting the H2 lines is degenerate: The shock velocities and number densities range between v_s = 17 - 47 km s^-1 and n_H=10^7 - 5 * 10^4 cm^-3, respectively. The H2 lines thus need a much more powerful shock than the CO lines. We deduce that most of the gas is currently moderately stirred up by slow (10 km s^-1) shocks while only a small fraction (< 1 percent) of the ISM is exposed to the high velocity shocks. This implies that the gas is rapidly loosing its highly turbulent motions. We argue that a high CO line-to-continuum ratio is a key diagnostic for the presence of shocks.

preprint2012arXiv

Investigation of Dual Active Nuclei, Outflows, Shock-Heated Gas, and Young Star Clusters in Markarian 266

Results of observations with the Spitzer, Hubble, GALEX, Chandra, and XMM-Newton space telescopes are presented for the Luminous Infrared Galaxy (LIRG) merger Mrk 266. The SW (Seyfert 2) and NE (LINER) nuclei reside in galaxies with Hubble types SBb (pec) and S0/a (pec), respectively. Both galaxies have L > L*, and they are inferred to each contain a ~2.5x10^8 M_sun black hole. Mrk 266 SW is likely the primary source of a bright Fe K-alpha line detected from the system, consistent with the reflection-dominated X-ray spectrum of a heavily obscured AGN. Optical knots embedded in an arc with aligned radio continuum radiation, combined with luminous H_2 line emission, provide evidence for a radiative bow shock in an AGN-driven outflow surrounding the NE nucleus. Soft X-ray emission modeled as shock-heated plasma is co-spatial with radio continuum emission between the galaxies. Mid-IR diagnostics indicate roughly equal contributions of AGN and starburst radiation powering the bolometric luminosity. Approximately 120 star clusters have been detected, with most having estimated ages < 50 Myr. Detection of 24 micron emission extending ~34 arcsec (20 kpc) north of the galaxies is interpreted as ~2x10^7 M_sun of dust entrained in an outflowing superwind; at optical wavelengths this region is resolved into a fragmented morphology indicative of Rayleigh-Taylor instabilities in an expanding shell of ionized gas. Mrk 266 demonstrates that the dust "blow-out" phase can begin in a LIRG well before the galaxies fully coalesce during a subsequent ULIRG phase, and rapid gas consumption in luminous dual AGNs (kpc scale separations) early in the merger process may explain the paucity of detected binary QSOs (sub-pc scale orbits) in large surveys. An evolutionary sequence is proposed representing a progression from dual to binary AGNs, accompanied by an increase in observed L_x/L_ir ratios by 10^4 or more.

preprint2012arXiv

Mapping the cold dust temperatures and masses of nearby Kingfish galaxies with Herschel

Taking advantage of the sensitivity and angular resolution of the Herschel Space Observatory at far-infrared and submm wavelengths, we aim to characterize the physical properties of cold dust within nearby galaxies and study the robustness of the parameters we derive using different modified blackbody models. For a pilot subsample of the KINGFISH program, we perform 2 temperature fits of the Spitzer and Herschel photometric data (24 to 500um), with a warm and a cold component, globally and in each resolution element.At global scales, we observe ranges of values for beta_c(0.8 to 2.5) and Tc(19.1 to 25.1K).We compute maps of our parameters with beta fixed or free to test the robustness of the temperature and dust surface density maps we deduce. When the emissivity is fixed, we observe temperature gradients as a function of radius.When the emissivity is fitted as a free parameter, barred galaxies tend to have uniform fitted emissivities.Gathering resolved elements in a Tc-beta_c diagram underlines an anti-correlation between the two parameters.It remains difficult to assess whether the dominant effect is the physics of dust grains, noise, or mixing along the line of sight and in the beam. We finally observe in both cases that the dust column density peaks in central regions of galaxies and bar ends (coinciding with molecular gas density enhancements usually found in these locations).We also quantify how the total dust mass varies with our assumptions about the emissivity index as well as the influence of the wavelength coverage used in the fits. We show that modified blackbody fits using a shallow emissivity (beta_c < 2.0) lead to significantly lower dust masses compared to the beta_c < 2.0 case, with dust masses lower by up to 50% if beta_c=1.5 for instance.The working resolution affects our total dust mass estimates: masses increase from global fits to spatially-resolved fits.

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

The Emission-Line Spectra of Major Mergers: Evidence for Shocked Outflows

Using a spectral decomposition technique (Soto & Martin 2012, hereafter Paper I), we investigate the physical origin of the high-velocity emission line gas in a sample of 39 gas-rich, ultraluminous infrared galaxy (ULIRG) mergers. Regions with shock-like excitation were identified in two kinematically distinct regimes, characterized by broad ($σ>$ 150 \kms) and narrow linewidths. Here we investigate the physical origin of the high-velocity (broad) emission with shock-like line ratios. Considering the large amount of extinction in these galaxies, the blueshift of the broad emission suggests an origin on the near side of the galaxy and therefore an interpretation as a galactic outflow. The large spatial extent of the broad, shocked emission component is generally inconsistent with an origin in the narrow-line region of a AGN, so we conclude that energy and momentum supplied by the starburst drives these outflows. The new data are used to examine the fraction of the supernova energy radiated by shocks and the mass loss rate in the warm-ionized phase of the wind. We show that the shocks produced by galactic outflows can be recognized in moderately high-resolution, integrated spectra of these nearby, ultraluminous starbursts. The spectral fitting technique introduced in Paper I may therefore be used to improve the accuracy of the physical properties measured for high-redshift galaxies from their (observed frame) infrared spectra.

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

A view of the narrow-line region in the infrared: active galactic nuclei with resolved fine-structure lines in the Spitzer archive

We queried the Spitzer archive for high-resolution observations with the Infrared Spectrograph of optically selected active galactic nuclei (AGN) for the purpose of identifying sources with resolved fine-structure lines that would enable studies of the narrow-line region (NLR) at mid-infrared wavelengths. By combining 298 Spitzer spectra with 6 Infrared Space Observatory spectra, we present kinematic information of the NLR for 81 z<=0.3 AGN. We used the [NeV], [OIV], [NeIII], and [SIV] lines, whose fluxes correlate well with each other, to probe gas photoionized by the AGN. We found that the widths of the lines are, on average, increasing with the ionization potential of the species that emit them. No correlation of the line width with the critical density of the corresponding transition was found. The velocity dispersion of the gas, sigma, is systematically higher than that of the stars, sigma_*, in the AGN host galaxy, and it scales with the mass of the central black hole, M_BH. Further correlations between the line widths and luminosities L, and between L and M_BH, are suggestive of a three dimensional plane connecting log(M_BH) to a linear combination of log(sigma) and log(L). Such a plane can be understood within the context of gas motions that are driven by AGN feedback mechanisms, or virialized gas motions with a power-law dependence of the NLR radius on the AGN luminosity. The M_BH estimates obtained for 35 type 2 AGN from this plane are consistent with those obtained from the M_BH-sigma_* relation.

preprint2011arXiv

Black Hole Masses and Star Formation Rates of z >1 Dust Obscured Galaxies (DOGs): Results from Keck OSIRIS Integral Field Spectroscopy

We have obtained high spatial resolution Keck OSIRIS integral field spectroscopy of four z~1.5 ultra-luminous infrared galaxies that exhibit broad H-alpha emission lines indicative of strong AGN activity. The observations were made with the Keck laser guide star adaptive optics system giving a spatial resolution of 0.1", or <1 kpc at these redshifts. These high spatial resolution observations help to spatially separate the extended narrow-line regions --- possibly powered by star formation --- from the nuclear regions, which may be powered by both star formation and AGN activity. There is no evidence for extended, rotating gas disks in these four galaxies. Assuming dust correction factors as high as A(H-alpha)=4.8 mag, the observations suggest lower limits on the black hole masses of (1 - 9) x 10^8 solar masses, and star formation rates <100 solar masses per year. The black hole masses and star formation rates of the sample galaxies appear low in comparison to other high-z galaxies with similar host luminosities. We explore possible explanations for these observations including, host galaxy fading, black hole growth, and the shut down of star formation.

preprint2011arXiv

C-GOALS: Chandra observations of a complete sample of luminous infrared galaxies from the IRAS Revised Bright Galaxy Survey

We present X-ray data for a complete sample of 44 luminous infrared galaxies (LIRGs), obtained with the Chandra X-ray Observatory. These are the X-ray observations of the high luminosity portion of the Great Observatory All-sky LIRG Survey (GOALS), which includes the most luminous infrared selected galaxies, log (Lir/Lsun) > 11.73, in the local universe, z < 0.088. X-rays were detected from 43 out of 44 objects, and their arcsec-resolution images, spectra, and radial brightness distributions are presented. With a selection by hard X-ray colour and the 6.4 keV iron line, AGN are found in 37% of the objects, with higher luminosity sources more likely to contain an AGN. These AGN also tend to be found in late-stage mergers. The AGN fraction would increase to 48% if objects with mid-IR [Ne V] detection are included. Double AGN are clearly detected only in NGC 6240 among 24 double/triple systems. Other AGN are found either in single nucleus objects or in one of the double nuclei at similar rates. Objects without conventional X-ray signatures of AGN appear to be hard X-ray quiet, relative to the X-ray to far-IR correlation for starburst galaxies, as discussed elsewhere. Most objects also show extended soft X-ray emission, which is likely related to an outflow from the nuclear region, with a metal abundance pattern suggesting enrichment by core collapse supernovae, as expected for a starburst.

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

GOODS-Herschel A population of 24micron dropout sources at z < 2

Using extremely deep PACS 100- and 160um Herschel data from the GOODS-Herschel program, we identify 21 infrared bright galaxies previously missed in the deepest 24um surveys performed by MIPS. These MIPS dropouts are predominantly found in two redshift bins, centred at z ~0.4 and ~1.3. Their S_100/S_24 flux density ratios are similar to those of local LIRGs and ULIRGs, whose silicate absorption features at 18um (at z ~ 0.4) and 9.7um (at z ~ 1.3) are shifted into the 24um MIPS band at these redshifts. The high-z sub-sample consists of 11 infrared luminous sources, accounting for ~2% of the whole GOODS-Herschel sample and putting strong upper limits on the fraction of LIRGs/ULIRGs at 1.0<z<1.7 that are missed by the 24um surveys. We find that a S_100/S_24 > 43 colour cut selects galaxies with a redshift distribution similar to that of the MIPS dropouts and when combined with a second colour cut, S_16/S_8 > 4, isolates sources at 1.0 < z < 1.7. We show that these sources have elevated specific star formation rates (sSFR) compared to main sequence galaxies at these redshifts and are likely to be compact starbursts with moderate/strong 9.7um silicate absorption features in their mid-IR spectra. Herschel data reveal that their infrared luminosities extrapolated from the 24um flux density are underestimated, on average, by a factor of ~3. These silicate break galaxies account for 16% (8%) of the ULIRG (LIRG) population in the GOODS fields, indicating a lower limit in their space density of 2.0 \times 10^(-5) Mpc^(-3). Finally, we provide estimates of the fraction of z < 2 MIPS dropout sources as a function of the 24-, 100-, 160-, 250- and 350um sensitivity limits, and conclude that previous predictions of a population of silicate break galaxies missed by the major 24um extragalactic surveys have been overestimated.

preprint2011arXiv

GOODS-Herschel: an infrared main sequence for star-forming galaxies

We present the deepest far-IR observations obtained with Herschel and examine the 3-500um SEDs of galaxies at 0<z<2.5, supplemented by a local reference sample from IRAS, ISO, Spitzer and AKARI data. We find that the ratio of total IR luminosity to rest-frame 8um luminosity, IR8 (=Lir/L8), follows a Gaussian distribution centered on IR8=4 and defines an IR main sequence (MS). A minority population (<20 %) of outliers producing a tail skewed toward higher values of IR8 consist of starbursts (SB) with compact projected star formation densities. IR8 can be used to separate galaxies with normal and extended modes of star formation from compact SBs with high-IR8, high projected IR surface brightness (>3x10^10 Lsun kpc^-2) and a high specific SFR (i.e., SBs). The rest-frame, UV-2700A size of these distant SBs is typically half that of MS galaxies, supporting the correlation between star formation density and SB activity that is measured for the local sample. Locally, (U)LIRGs are systematically in the SB mode, whereas most distant (U)LIRGs form stars in the "normal" MS mode. This confusion between two modes of star formation is the cause of the so-called "mid-IR excess" population of galaxies found at z>1.5 by previous studies. MS galaxies have strong PAH emission line features, a broad far-IR bump resulting from a combination of dust temperatures (Tdust~15-50 K), and an effective Tdust~31 K, as derived from the peak wavelength of their IR SED. Galaxies in the SB regime instead exhibit weak PAH EW and a sharper far-IR bump with an effective Tdust~40 K. Finally, we present evidence that the mid-to-far IR emission of X-ray AGNs is predominantly produced by star formation and that candidate dusty AGNs with a power-law emission in the mid-IR systematically occur in compact, dusty SBs. After correcting for the effect of SBs on IR8, we identify new candidates for extremely obscured AGNs.

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

HST Morphologies of z ~ 2 Dust-Obscured Galaxies II: Bump Sources

We present Hubble Space Telescope (HST) imaging of 22 ultra-luminous infrared galaxies (ULIRGs) at z~2 with extremely red R-[24] colors (called dust-obscured galaxies, or DOGs) which have a local maximum in their spectral energy distribution (SED) at rest-frame 1.6um associated with stellar emission. These sources, which we call "bump DOGs", have star-formation rates of 400-4000 Msun/yr and have redshifts derived from mid-IR spectra which show strong polycyclic aromatic hydrocarbon emission --- a sign of vigorous on-going star-formation. Using a uniform morphological analysis, we look for quantifiable differences between bump DOGs, power-law DOGs (Spitzer-selected ULIRGs with mid-IR SEDs dominated by a power-law and spectral features that are more typical of obscured active galactic nuclei than starbursts), sub-millimeter selected galaxies (SMGs), and other less-reddened ULIRGs from the Spitzer extragalactic First Look Survey (XFLS). Bump DOGs are larger than power-law DOGs (median Petrosian radius of 8.4 +/- 2.7 kpc vs. 5.5 +/- 2.3 kpc) and exhibit more diffuse and irregular morphologies (median M_20 of -1.08 +/- 0.05 vs. -1.48 +/- 0.05). These trends are qualitatively consistent with expectations from simulations of major mergers in which merging systems during the peak star-formation rate period evolve from M_20 = -1.0 to M_20 = -1.7. Less obscured ULIRGs (i.e., non-DOGs) tend to have more regular, centrally peaked, single-object morphologies rather than diffuse and irregular morphologies. This distinction in morphologies may imply that less obscured ULIRGs sample the merger near the end of the peak star-formation rate period. Alternatively, it may indicate that the intense star-formation in these less-obscured ULIRGs is not the result of a recent major merger.

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

Multi-wavelength GOALS Observations of Star Formation and Active Galactic Nucleus Activity in the Luminous Infrared Galaxy IC 883

New optical HST, Spitzer, GALEX, and Chandra observations of the single-nucleus, luminous infrared galaxy (LIRG) merger IC 883 are presented. The galaxy is a member of the Great Observatories All-sky LIRG Survey (GOALS), and is of particular interest for a detailed examination of a luminous late-stage merger due to the richness of the optically-visible star clusters and the extended nature of the nuclear X-ray, mid-IR, CO and radio emission. In the HST ACS images, the galaxy is shown to contain 156 optically visible star clusters distributed throughout the nuclear regions and tidal tails of the merger, with a majority of visible clusters residing in an arc ~ 3-7 kpc from the position of the mid-infrared core of the galaxy. The luminosity functions of the clusters have an alpha_F435W ~ -2.17+/-0.22 and alpha_F814W ~ -2.01+/-0.21. Further, the colors and absolute magnitudes of the majority of the clusters are consistent with instantaneous burst population synthesis model ages in the range of a few x10^7 - 10^8 yrs (for 10^5 M_sun clusters), but may be as low as few x10^6 yrs with extinction factored in. The X-ray and mid-IR spectroscopy are indicative of predominantly starburst-produced nuclear emission, and the star formation rate is ~ 80 M_sun / yr. The kinematics of the CO emission and the morphology of both the CO and radio emission are consistent with the nuclear starburst being situated in a highly inclined disk 2 kpc in diameter with an infrared surface brightness mu_IR ~ 2x10^11 L_sun kpc^-2, a factor of 10 less than that of the Orion star-forming region. Finally, the detection of the [Ne V] 14.32 um emission line is evidence that an AGN is present. The faintness of the line (i.e., [Ne V] / [Ne II] um ~ 0.01) and the small equivalent width of the 6.2 um PAH feature ($= 0.39μ$m) are both indicative of a relatively weak AGN. (abridged)

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 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 location of an active nucleus and the soft X-ray shadowing by a tidal tail in the ULIRG Mrk 273

Analysis of data from the Chandra X-ray Observatory for the double nucleus ULIRG Mrk 273 reveals an absorbed hard X-ray source coincident with the southwest nucleus, implying that this unresolved near infrared source is where an active nucleus resides while the northern nuclear region contains a powerful starburst which dominates the far infrared luminosity. There is evidence of a slight image extension in the 6-7 keV band, where a Fe K line is present, towards the northern nucleus. A large-scale, diffuse emission nebula detected in soft X-rays contains a dark lane that spatially coincides with a high surface-brightness tidal tail extending ~50 arcsec (40 kpc) to the south. The soft X-ray source is likely located behind the tidal tail which absorbs X-ray photons along the line of sight. The estimated column density of cold gas in the tidal tail responsible for shadowing the soft X-rays is nH >= 6e+21 cm-2, consistent with the tidal tail having an edge-on orientation.

preprint2011arXiv

The Star-Formation Histories of z~2 DOGs and SMGs

The Spitzer Space Telescope has identified a population of ultra-luminous infrared galaxies (ULIRGs) at z ~ 2 that may play an important role in the evolution of massive galaxies. We measure the stellar masses of two populations of Spitzer-selected ULIRGs, both of which have extremely red R-[24] colors (dust-obscured galaxies, or DOGs) and compare our results with sub-millimeter selected galaxies (SMGs). One set of 39 DOGs has a local maximum in their mid-IR spectral energy distribution (SED) at rest-frame 1.6um associated with stellar emission ("bump DOGs"), while the other set of 51 DOGs has a power-law dominated mid-IR SED with spectral features typical of obscured AGN ("power-law DOGs"). We use stellar population synthesis models applied self-consistently to broad-band photometry in the rest-frame ultra-violet, optical, and near-infrared of each of these populations and test a variety of stellar population synthesis codes, star-formation histories (SFHs), and initial mass functions (IMFs). Assuming a simple stellar population SFH and a Chabrier IMF, we find that the median and inner quartile stellar masses of SMGs, bump DOGs and power-law DOGs are given by log(M_*/M_sun) = 10.42_-0.36^+0.42, 10.62_-0.32^+0.36, and 10.71_-0.34^+0.40, respectively. Implementing more complicated SFHs with multiple age components increases these mass estimates by up to 0.5 dex. Our stellar mass estimates are consistent with physical mechanisms for the origin of z~2 ULIRGs that result in high star-formation rates for a given stellar mass. Such mechanisms are usually driven by a major merger of two gas-rich systems, rather than smooth accretion of gas and small satellites.

preprint2011arXiv

The Stellar Mass Content of Submillimeter-Selected Galaxies

We present a new study of the stellar mass in a sample of ~ 70 submillimeter-selected galaxies (SMGs) with accurate spectroscopic redshifts. We fit combinations of stellar population synthesis models and power laws to the galaxies' observed-frame optical through mid-IR spectral energy distributions to separate stellar emission from non-stellar near-IR continuum. By separating the stellar emission from the non-stellar near-IR continuum, we find that ~ 50% of our sample have non-stellar continuum contributions of less than 10% in rest-frame H-band, but ~ 10% of our sample have non-stellar contributions greater than 50%. We find that the K-band luminosity of the non-stellar continuum emission is correlated with hard X-ray luminosity, indicating an AGN origin of the emission. Upon subtracting this AGN-contributed continuum component from all of the galaxies in our sample, we determine a lower median stellar mass for SMGs than previous studies, ~ 7 x 10^10 M_sun. We use constraints of the starburst time-scale from molecular gas studies to estimate the amount of fading our sample would undergo if they passively evolve after the starburst terminates. The results suggest that typical SMGs, while among the most massive galaxies at z ~ 2, are likely to produce descendants of similar mass and luminosity to L* galaxies in the local universe.

preprint2011arXiv

X-ray observations of highly obscured 9.7 micron sources: an efficient method for selecting Compton-thick AGN ?

Spitzer/IRS has revealed many sources with very deep Si features at 9.7micron (tau>1). We set out to investigate whether a strong Si absorption feature is a good indicator for the presence of a heavily obscured AGN. We compile X-ray spectroscopic observations available in the literature on the optically-thick,tau(9.7)>1 sources from the IRAS Seyfert sample. We find that the majority of the high-tau optically confirmed Seyferts (6/9) in this sample are probably CT. Thus we provide direct evidence for a connection between mid-IR optically-thick galaxies and CT AGN, with the success rate being close to 70% in the local Universe. This is at least comparable, if not better, than other rates obtained with photometric information in the mid to far-IR, or even mid-IR to Xray. However, this technique cannot provide complete CT AGN samples,ie there are many CT AGN which do not show significant Si absorption, with the most notable example being N1068. Having assessed the validity of the high 9.7micron technique locally, we attempt to construct a sample of candidate CT AGN at higher redshifts. We compile a sample of 7 high-tau sources in the GOODS and 5 in the Spitzer FLS. All these have been selected to have no PAH features EW(6.2)<0.3 in order to maximize the probability that they are AGN. 6 out of 7 sources in the GOODS have been detected in X-rays, while for the five FLS sources only X-ray flux upper limits are available. The high X-ray luminosities of the detected GOODS sources corroborates that these are AGN. For FLS, ancillary optical spectroscopy reveals hidden nuclei in two more sources. SED fitting can support the presence of an AGN in the vast majority of sources. We cannot derive useful X-ray spectroscopy constraints on whether these are CT. However, the low LX/L6 ratios, suggest that at least 4 out of the 6 detected sources in GOODS may be associated with CT AGN.

preprint2010arXiv

Black hole accretion and star formation as drivers of gas excitation and chemistry in Mrk231

We present a full high resolution SPIRE FTS spectrum of the nearby ultraluminous infrared galaxy Mrk231. In total 25 lines are detected, including CO J=5-4 through J=13-12, 7 rotational lines of H2O, 3 of OH+ and one line each of H2O+, CH+, and HF. We find that the excitation of the CO rotational levels up to J=8 can be accounted for by UV radiation from star formation. However, the approximately flat luminosity distribution of the CO lines over the rotational ladder above J=8 requires the presence of a separate source of excitation for the highest CO lines. We explore X-ray heating by the accreting supermassive black hole in Mrk231 as a source of excitation for these lines, and find that it can reproduce the observed luminosities. We also consider a model with dense gas in a strong UV radiation field to produce the highest CO lines, but find that this model strongly overpredicts the hot dust mass in Mrk231. Our favoured model consists of a star forming disk of radius 560 pc, containing clumps of dense gas exposed to strong UV radiation, dominating the emission of CO lines up to J=8. X-rays from the accreting supermassive black hole in Mrk231 dominate the excitation and chemistry of the inner disk out to a radius of 160 pc, consistent with the X-ray power of the AGN in Mrk231. The extraordinary luminosity of the OH+ and H2O+ lines reveals the signature of X-ray driven excitation and chemistry in this region.

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

Herschel observations of water vapour in Markarian 231

The Ultra Luminous InfraRed Galaxy Mrk 231 reveals up to seven rotational lines of water (H2O) in emission, including a very high-lying (E_{upper}=640 K) line detected at a 4sigma level, within the Herschel/SPIRE wavelength range, whereas PACS observations show one H2O line at 78 microns in absorption, as found for other H2O lines previously detected by ISO. The absorption/emission dichotomy is caused by the pumping of the rotational levels by far-infrared radiation emitted by dust, and subsequent relaxation through lines at longer wavelengths, which allows us to estimate both the column density of H2O and the general characteristics of the underlying far-infrared continuum source. Radiative transfer models including excitation through both absorption of far-infrared radiation emitted by dust and collisions are used to calculate the equilibrium level populations of H2O and the corresponding line fluxes. The highest-lying H2O lines detected in emission, with levels at 300-640 K above the ground state, indicate that the source of far-infrared radiation responsible for the pumping is compact (radius=110-180 pc) and warm (T_{dust}=85-95 K), accounting for at least 45% of the bolometric luminosity. The high column density, N(H2O)~5x10^{17} cm^{-2}, found in this nuclear component, is most probably the consequence of shocks/cosmic rays, an XDR chemistry, and/or an "undepleted chemistry" where grain mantles are evaporated. A more extended region, presumably the inner region of the 1-kpc disk observed in other molecular species, could contribute to the flux observed in low-lying H2O lines through dense hot cores, and/or shocks. The H2O 78 micron line observed with PACS shows hints of a blue-shifted wing seen in absorption, possibly indicating the occurrence of H2O in the prominent outflow detected in OH (Fischer et al., this volume).

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

Mid-Infrared Properties of the Swift Burst Alert Telescope Active Galactic Nuclei Sample of the Local Universe. I. Emission-Line Diagnostics

We compare mid-infrared emission-line properties, from high-resolution Spitzer spectra of a hard X-ray (14 -- 195 keV) selected sample of nearby (z < 0.05) AGN detected by the Burst Alert Telescope (BAT) aboard Swift. The luminosity distribution for the mid-infrared emission-lines, [O IV] 25.89 micron, [Ne II] 12.81 micron, [Ne III] 15.56 micron and [Ne V] 14.32/24.32 micron, and hard X-ray continuum show no differences between Seyfert 1 and Seyfert 2 populations, however six newly discovered BAT AGNs are under-luminous in [O IV], most likely the result of dust extinction in the host galaxy. The overall tightness of the mid-infrared correlations and BAT fluxes and luminosities suggests that the emission lines primarily arise in gas ionized by the AGN. We also compare the mid-infrared emission-lines in the BAT AGNs with those from published studies of ULIRGs, PG QSOs, star-forming galaxies and LINERs. We find that the BAT AGN sample fall into a distinctive region when comparing the [Ne III]/[Ne II] and the [O IV]/[Ne III] ratios. These line ratios are lower in sources that have been previously classified in the mid-infrared/optical as AGN than those found for the BAT AGN, suggesting that, in our X-ray selected sample, the AGN represents the main contribution to the observed line emission. These ratios represent a new emission line diagnostic for distinguishing between AGN and star forming galaxies.

preprint2010arXiv

Mid-Infrared Spectral Diagnostics of Luminous Infrared Galaxies

We present a statistical analysis of the mid-infrared (MIR) spectra of 248 luminous infrared (IR) galaxies (LIRGs) which comprise the Great Observatories All-sky LIRG Survey (GOALS) observed with the Infrared Spectrograph (IRS) on-board the Spitzer Space Telescope. The GOALS sample enables a direct measurement of the relative contributions of star-formation and active galactic nuclei (AGN) to the total IR emission from a large sample of local LIRGs. The AGN contribution to the MIR emission (f-AGN) is estimated by employing several diagnostics based on the properties of the [NeV], [OIV] and [NeII] fine structure gas emission lines, the 6.2 microns PAH and the shape of the MIR continuum. We find that 18% of all LIRGs contain an AGN and that in 10% of all sources the AGN contributes more than 50% of the total IR luminosity. Summing up the total IR luminosity contributed by AGN in all our sources suggests that AGN supply ~12% of the total energy emitted by LIRGs. The average spectrum of sources with an AGN looks similar to the average spectrum of sources without an AGN, but it has lower PAH emission and a flatter MIR continuum. AGN dominated LIRGs have higher IR luminosities, warmer MIR colors and are found in interacting systems more often than pure starbursts LIRGs. However we find no linear correlations between these properties and f-AGN. We used the IRAC colors of LIRGs to confirm that finding AGN on the basis of their MIR colors may miss ~40% of AGN dominated (U)LIRGs

preprint2010arXiv

On the black hole mass - velocity dispersion relation for type-1 and type-2 AGN

We present results from infrared spectroscopic projects that aim to test the relation between the mass of a black hole, M_BH, and the velocity dispersion of the stars in its host-galaxy bulge. We demonstrate that near-infrared, high-resolution spectroscopy assisted by adaptive optics is key in populating the high-luminosity end of the relation. We show that the velocity dispersions of mid-infrared, high-ionization lines originating from gas in the narrow-line region of the active galactic nucleus follow the same relation. This result provides a way of inferring M_BH estimates for the cosmologically significant population of obscured, type-2 AGN that can be applicable to data from spectrographs on the next generation infrared telescopes.

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 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

A Spitzer high resolution mid-infrared spectral atlas of starburst galaxies

We present an atlas of Spitzer/IRS high resolution (R~600) 10-37um spectra for 24 well known starburst galaxies. The spectra are dominated by fine-structure lines, molecular hydrogen lines, and emission bands of polycyclic aromatic hydrocarbons. Six out of the eight objects with a known AGN component show emission of the high excitation [NeV] line. This line is also seen in one other object (NGC4194) with, a priori, no known AGN component. In addition to strong polycyclic aromatic hydrocarbon emission features in this wavelength range (11.3, 12.7, 16.4um), the spectra reveal other weak hydrocarbon features at 10.6, 13.5, 14.2um, and a previously unreported emission feature at 10.75um. An unidentified absorption feature at 13.7um is detected in many of the starbursts. We use the fine-structure lines to derive the abundance of neon and sulfur for 14 objects where the HI 7-6 line is detected. We further use the molecular hydrogen lines to sample the properties of the warm molecular gas. Several basic diagrams characterizing the properties of the sample are also shown. We have combined the spectra of all the pure starburst objects to create a high S/N template, which is available to the community.

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

HST Morphologies of z~2 Dust Obscured Galaxies I: Power-law Sources

We present high spatial resolution optical and near-infrared imaging obtained using the ACS, WFPC2 and NICMOS cameras aboard the Hubble Space Telescope of 31 24um--bright z~2 Dust Obscured Galaxies (DOGs) identified in the Bootes Field of the NOAO Deep Wide-Field Survey. Although this subset of DOGs have mid-IR spectral energy distributions dominated by a power-law component suggestive of an AGN, all but one of the galaxies are spatially extended and not dominated by an unresolved component at rest-frame UV or optical wavelengths. The observed V-H and I-H colors of the extended components are 0.2-3 magnitudes redder than normal star-forming galaxies. All but 1 have axial ratios >0.3, making it unlikely that DOGs are composed of an edge-on star-forming disk. We model the spatially extended component of the surface brightness distributions of the DOGs with a Sersic profile and find effective radii of 1-6 kpc. This sample of DOGs is smaller than most sub-millimeter galaxies (SMGs), but larger than quiescent high-redshift galaxies. Non-parametric measures (Gini and M20) of DOG morphologies suggest that these galaxies are more dynamically relaxed than local ULIRGs. We estimate lower limits to the stellar masses of DOGs based on the rest-frame optical photometry and find that these range from ~10^(9-11) M_sun. If major mergers are the progenitors of DOGs, then these observations suggest that DOGs may represent a post-merger evolutionary stage.

preprint2009arXiv

Local Lyman Break Galaxy Analogs: The Impact of Massive Star-forming Clumps on the Interstellar Medium and the Global Structure of Young, Forming Galaxies

We present HST UV/optical imaging, Spitzer mid-IR photometry, and optical spectroscopy of a sample of 30 low-redshift (z=0.1-0.3) galaxies chosen from SDSS/GALEX to be accurate local analogs of the high-z Lyman Break Galaxies. The Lyman Break Analogs (LBAs) are similar in mass, metallicity, dust, SFR, size and gas velocity dispersion, thus enabling a detailed investigation of processes that are important at high-z. The optical emission line properties of LBAs are also similar to those of LBGs, indicating comparable conditions in their ISM. In the UV, LBAs are characterized by complexes of massive star-forming "clumps", while in the optical they most often show evidence for (post-)mergers/interactions. In 6 cases, we find an extremely massive (>10^9 Msun) compact (R~100 pc) dominant central object (DCO). The DCOs are preferentially found in LBAs with the highest mid-IR luminosities and correspondingly high SFRs (15-100 Msun/yr). We show that the massive SF clumps (including the DCOs) have masses much larger than the nuclear super star clusters seen in normal late type galaxies. However, the DCOs have masses, sizes, and densities similar to the excess-light/central-cusps seen in typical elliptical galaxies with masses similar to the LBA galaxies. We suggest that the DCOs form in present-day examples of the dissipative mergers at high redshift that are believed to have produced the central-cusps in local ellipticals. More generally, the properties of the LBAs are consistent with the idea that instabilities in a gas-rich disk lead to very massive star-forming clumps that eventually coalesce to form a spheroid. We speculate that the DCOs are too young at present to be growing a supermassive black hole because they are still in a supernova-dominated outflow phase.

preprint2009arXiv

Spitzer-IRS Study of the Antennae Galaxies NGC 4038/39

Using the Infrared Spectrograph on the Spitzer Space Telescope, we observed the Antennae galaxies obtaining spectral maps of the entire central region and high signal-to-noise 5-38 um spectra of the two galactic nuclei and six infrared-luminous regions. The total infrared luminosity of our six IR peaks plus the two nuclei is L_IR = 3.8x10^10 L_o, with their derived star formation rates ranging between 0.2 and 2 M_o/yr, with a total of 6.6 M_o/yr. The hardest and most luminous radiation originates from two compact clusters in the southern part of the overlap region, which also have the highest dust temperatures. PAH emission and other tracers of softer radiation are spatially extended throughout and beyond the overlap region, but regions with harder and intenser radiation field show a reduced PAH strength. The strong H_2 emission is rather confined around the nucleus of NGC 4039, where shocks appear to be the dominant excitation mechanism, and the southern part of the overlap region, where it traces the most recent starburst activity. The luminosity ratio between the warm molecular gas (traced by the H_2 lines) and the total far-IR emission is ~1.6x10^-4, similar to that found in many starburst and ULIRGs. The total mass of warm H_2 in the Antennae is 2.5x10^7 M_o, with a fraction of warm to total H_2 gas mass of about 0.35%. The average warm H_2 temperature is 302+/-26 K and appears anti-correlated with the radiation field hardness, possibly due to an evolution of the PDR morphology. The previously reported tight correlation between the H_2 and PAH emission was not found but higher total PAH emission to continuum ratios were found in PDRs with warmer gas.

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

High-velocity neon line emission from the ULIRG IRAS F00183-7111: revealing the optically obscured base of a nuclear outflow

We report the first mid-IR detection of highly disturbed ionized gas in the ultraluminous infrared galaxy IRAS F00183-7111. The gas, traced by the 12.81um [NeII] and 15.56um [NeIII] lines, spans a velocity range of-3500 to 3000 km/s with respect to systemic velocity. Optical and near-IR spectroscopic studies show no evidence for similarly high velocity gas components in forbidden lines at shorter wavelengths. We interpret this as the result of strong extinction (Av=10-50) on the high-velocity gas, which identifies the base of the outflow traced in 5007A [OIII] as a plausible origin. Unusual excitation conditions are implied by a comparison of the mid-infrared low-excitation neon line emission and the PAH emission for a sample of 56 ULIRGs. For IRAS F00183, the neon/PAH ratio is 8 times higher than the average ratio. Similar mid-infrared kinematic and excitation characteristics are found for only 2 other ULIRGs in our sample: IRAS 12127-1412NE and IRAS 13451+1232. Both sources have an elevated neon/PAH ratio and exhibit pronounced blue wings in their 15.56um [NeIII] line profiles. IRAS 13451 even shows a strongly blue shifted and broad 14.32um [NeV] line. While for IRAS 13451 the observed [NeIII]/[NeII] and [NeV]/[NeII] line ratios indicate exposure of the blue shifted gas to direct radiation from the AGN, for IRAS F00183 and 12127 the neon line ratios are consistent with an origin in fast shocks in a high-density environment, and with an evolutionary scenario in which strongly blue shifted [Ne II] and [Ne III] emission trace the (partial) disruption of the obscuring medium around a buried AGN. The detection of strongly blue shifted [Ne V] emission in IRAS 13451 would then indicate this process to be much further advanced in this ULIRG than in IRAS F00183 and 12127, where this line is undetected.

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

A Spitzer Spectrum of the Exoplanet HD 189733b

We report on the measurement of the 7.5-14.7 micron spectrum for the transiting extrasolar giant planet HD 189733b using the Infrared Spectrograph on the Spitzer Space Telescope. Though the observations comprise only 12 hours of telescope time, the continuum is well measured and has a flux ranging from 0.6 mJy to 1.8 mJy over the wavelength range, or 0.49 +/- 0.02% of the flux of the parent star. The variation in the measured fractional flux is very nearly flat over the entire wavelength range and shows no indication of significant absorption by water or methane, in contrast with the predictions of most atmospheric models. Models with strong day/night differences appear to be disfavored by the data, suggesting that heat redistribution to the night side of the planet is highly efficient.

preprint2007arXiv

Dust and Atomic Gas in Dwarf Irregular Galaxies of the M81 Group: The SINGS and THINGS view

[abridged] We present observations of the dust and atomic gas phase in seven dwarf irregular galaxies of the M81 group from the SINGS and THINGS surveys. The Spitzer observations provide a first glimpse of the nature of the non-atomic ISM in these metal-poor (Z~0.1 Z_sun), quiescent (SFR~0.001-0.1 M_sun/yr) dwarf galaxies. Dust emission is detected in five out of the seven targets. Most detected dust emission is restricted to HI column densities >1x10^21 cm^-2. Spitzer spectroscopy of two regions in the brightest galaxies (IC 2574 and Holmberg II) show distinctly different spectral shapes. The spectrum of IC 2574 shows aromatic features that are less luminous (relative to the FIR luminosity) compared to an average SINGS spiral galaxy by a factor of \~7 . The aromatic features in Holmberg~II (which has only a slightly lower gas-phase metallicity) are fainter than in IC 2574 by an order of magnitude. This result emphazises that the strength of the aromatic features is not a simple linear function of metallicity. We estimate dust masses of ~10^4-10^6 M_sun for the M81 dwarf galaxies, resulting in an average dust--to--gas ratio (M_dust/M_HI) of ~3x10^-4 (1.5x10^-3 if only the HI that is associated with dust emission is considered); this is an order of magnitude lower than the typical value derived for the SINGS spirals. The dwarf galaxies are underluminous per unit star formation rate at 70um as compared to the more massive galaxies in SINGS by a factor of ~2. However, the average 70um/160um ratio in the sample dwarf galaxies is higher than what is found in the other galaxies of the SINGS sample. This can be explained by a combination of a lower dust content in conjunction with a higher dust temperature in the dwarfs.

preprint2007arXiv

Infrared Molecular Starburst Fingerprints in Deeply Obscured (U)LIRG Nuclei

High resolution spectra of the Spitzer Space Telescope show vibration-rotation absorption bands of gaseous C2H2, HCN, and CO2 molecules toward a sample of deeply obscured (U)LIRG nuclei. The observed bands reveal the presence of dense (n>~ 10^7 cm^-3), warm (T_ex = 200-700 K) molecular gas with high column densities of these molecules ranging from a few 10^15 - 10^17 cm^-2. Abundances relative to H2, inferred from the silicate optical depth, range from ~10^-7 to 10^-6 and show no correlation with temperature. Theoretical studies show that the high abundances of both C2H2 and HCN exclude a X-ray dominated region (XDR) associated with the toroid surrounding an AGN as the origin of this dense warm molecular gas. Galactic massive protostars in the so-called Hot Core phase have similar physical characteristics with comparable high abundances of C2H2, HCN, and CO2 in the hot phase. However, the abundances of C2H2 and HCN and the C2H2/CO2 and HCN/CO2 ratios are much higher toward the (U)LIRGs in the cooler (T_ex <= 400 K) phase. We suggest that the warm dense molecular gas revealed by the mid-IR absorption lines is associated with a phase of deeply embedded star formation where the extreme pressures and densities of the nuclear starburst environment have inhibited the expansion of HII regions and the global disruption of the star forming molecular cloud cores, and `trapped' the star formation process in an `extended' Hot Core phase.

preprint2007arXiv

Measuring PAH Emission in Ultradeep Spitzer IRS Spectroscopy of High Redshift IR Luminous Galaxies

The study of the dominant population of high redshift IR-luminous galaxies (10^11 - 10^12 Lsun at 1<z<3), requires observation of sources at the ~0.1 mJy level in the mid-IR. We present the deepest spectra taken to date with the Infrared Spectrograph (IRS) on the Spitzer Space Telescope. We targeted two faint (f24~0.15 mJy) sources in the Southern GOODS field at z=1.09 and z=2.69. Spectra of the lower redshift target were taken in the observed-frame 8--21 micron range, while the spectrum of the higher redshift target covered 21--37 microns. We also present the spectra of two secondary sources within the slit. We detect strong PAH emission in all four targets, and compare the spectra to those of local galaxies observed by the IRS. The z=1.09 source appears to be a typical, star-formation dominated IR-luminous galaxy, while the z=2.69 source is a composite source with strong star formation and a prominent AGN. The IRAC colors of this source show no evidence of rest-frame near-infrared stellar photospheric emission. We demonstrate that an AGN which contributes only a small (~10%) fraction of the bolometric luminosity can produce enough hot dust emission to overwhelm the near-IR photospheric emission from stars. Such sources would be excluded from photometric surveys which rely on the near-IR bump to identify starbursts, leading to an underestimate of the star formation rate density.

preprint2007arXiv

Optical Line Diagnostics of z~2 Optically Faint ULIRGs in the Spitzer Bootes Survey

We present near-infrared spectroscopic observations for a sample of ten optically faint luminous infrared galaxies (R-[24]> 14) using Keck NIRSPEC and Gemini NIRI. The sample is selected from a 24 micron Spitzer MIPS imaging survey of the NDWFS Bootes field. We measure accurate redshifts in the range 1.3<z<3.4. Based on either emission line widths or line diagnostics, we find that all ten galaxies harbor luminous AGN. Seven sources are type I AGN, exhibiting broad (>1900 km/s) Halpha or Hbeta emission lines; the remaining three are type II AGN. Given their large mid-IR luminosities and faint optical magnitudes, we might expect these sources to be heavily extincted quasars, and therefore only visible as type II AGN. The visibility of broad lines in 70% of the sources suggests that it is unlikely that these AGN are being viewed through the mid-plane of a dusty torus. For four of the sources we constrain the Halpha/Hbeta Balmer decrement and estimate the extinction to the emission line region to be large for both type I and type II AGN, with A_Halpha > 2.4-5 mag. Since the narrow-line region is also extincted and the UV continuum emission from the host galaxies is extremely faint, this suggests that much of the obscuration is contributed by dust on large (~kpc) scales within the host galaxies. These sources may be examples of "host-obscured" AGN which could have space densities comparable or greater to that of optically luminous type I AGN with similar bolometric luminosities.

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

The Role of Galaxy Interactions and Mergers in Star Formation at z<1.3: Mid-Infrared Properties in the Spitzer First Look Survey

By combining the 0.12 square degree F814W Hubble Space Telescope (HST) and Spitzer MIPS 24 micron imaging in the First Look Survey (FLS), we investigate the properties of interacting and merging Mid-Infrared bright and faint sources at 0.2< z <1.3. We find a marginally significant increase in the pair fraction for MIPS 24 micron detected, optically selected close pairs, pair fraction=0.25+/-0.10 at z~1, in contrast to 0.11+/-0.08 at z~0.4, while galaxies below our 24 micron MIPS detection limit show a pair fraction consistent with zero at all redshifts. Additionally, 24 micron detected galaxies with fluxes >0.1mJy are on average five times more likely to be in a close galaxy pair between 0.2< z <1.3 than galaxies below this flux limit. Using the 24 micron flux to derive the total Far-IR luminosity we find that paired galaxies (early stage mergers) are responsible for 27% +/-9% of the IR luminosity density resulting from star formation at z~1 while morphologically classified (late stage) mergers make up 34%+/-11%. This implies that 61%+/-14% of the infrared luminosity density and in turn ~40% of the star formation rate density at z~1 can be attributed to galaxies at some stage of a major merger or interaction. We argue that, close pairs/mergers in a LIRG/ULIRG phase become increasingly important contributers to the IR luminosity and star formation rate density of the Universe at z>0.7.

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

Observations of Ultraluminous Infrared Galaxies with the Infrared Spectrograph on the Spitzer Space Telescope II: The IRAS Bright Galaxy Sample

We present spectra taken with the Infrared Spectrograph on Spitzer covering the 5-38 micron region of the ten Ultraluminous Infrared Galaxies (ULIRGs) found in the IRAS Bright Galaxy Sample. Among the BGS ULIRGs, we find a factor of 50 spread in the rest-frame mid to far-infrared spectral slope. The 9.7 micron silicate optical depths range from less than 0.4 more than 4.2, implying line of sight extinctions of A(V) ~ 8 - 78 mag. There is evidence for water ice and hydrocarbon absorption and C2H2 and HCN absorption features in four and possibly six of the 10 BGS ULIRGs, indicating shielded molecular clouds and a warm, dense ISM. We have detected [NeV] emission in three of the ten BGS ULIRGs, at flux levels of 5-18E-14 erg/cm^2/sec and [NeV] 14.3/[NeII] 12.8 line flux ratios of 0.12-0.85. The remaining BGS ULIRGs have limits on their [NeV]/[NeII] line flux ratios which range from less than 0.15 to less than 0.01. Among the BGS ULIRGs, the AGN fractions implied by either the [NeV]/[NeII] or [OIV]/[NeII] line flux ratios (or their upper limits) are significantly lower than implied by the MIR slope or EQW of the 6.2 micron PAH feature. Fitting the SEDs, we see evidence for hot (T > 300K) dust in five of the BGS ULIRGs, with the fraction of hot dust to total dust luminosity ranging from ~1-23%, before correcting for extinction. When integrated over the IRAC-8, IRS blue peakup, and MIPS-24 filter bandpasses, the IRS spectra imply very blue colors for some ULIRGs at z ~ 1.3. This is most extreme for sources with significant amounts of warm dust and deep silicate absorption.

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.

preprint2006arXiv

The Mid-IR Properties of Starburst Galaxies from Spitzer-IRS Spectroscopy

We present 5-38um mid-infrared spectra at a spectral resolution of R~65-130 of a large sample of 22 starburst nuclei taken with the Infrared Spectrograph IRS on board the Spitzer Space Telescope. The spectra show a vast range in starburst SEDs. The silicate absorption ranges from essentially no absorption to heavily obscured systems with an optical depth of tau(9.8um)~5. The spectral slopes can be used to discriminate between starburst and AGN powered sources. The monochromatic continuum fluxes at 15um and 30um enable a remarkably accurate estimate of the total infrared luminosity of the starburst. We find that the PAH equivalent width is independent of the total starburst luminosity L_IR as both continuum and PAH feature scale proportionally. However, the luminosity of the 6.2um feature scales with L_IR and can be used to approximate the total infrared luminosity of the starburst. Although our starburst sample covers about a factor of ten difference in the [NeIII]/[NeII] ratio, we found no systematic correlation between the radiation field hardness and the PAH equivalent width or the 7.7um/11.3um PAH ratio. These results are based on spatially integrated diagnostics over an entire starburst region, and local variations may be ``averaged out''. It is presumably due to this effect that unresolved starburst nuclei with significantly different global properties appear spectrally as rather similar members of one class of objects.

preprint2005arXiv

Identifying silicate-absorbed ULIRGs at z~1-2 in the Bootes Field using Spitzer/IRS

Using the 16$μ$m peakup imager on the Infrared Spectrograph (IRS) on Spitzer, we present a serendipitous survey of 0.0392 deg$^{2}$ within the area of the NOAO Deep Wide Field Survey in Bootes. Combining our results with the available Multiband Imaging Photometer for Spitzer (MIPS) 24$μ$m survey of this area, we produce a catalog of 150 16$μ$m sources brighter than 0.18 mJy (3$σ$) for which we derive measures or limits on the 16/24$μ$m colors. Such colors are especially useful in determining redshifts for sources whose mid infrared spectra contain strong emission or absorption features that characterize these colors as a function of redshift. We find that the 9.7$μ$m silicate absorption feature in Ultraluminous Infrared Galaxies (ULIRGs) results in sources brighter at 16$μ$m than at 24$μ$m at z $\sim$ 1--1.8 by at least 20%. With a threshold flux ratio of 1.2, restricting our analysis to $>5σ$ detections at 16$μ$m, and using a $3σ$ limit on 24$μ$m non-detections, the number of silicate-absorbed ULIRG candidates is 36. This defines a strong upper limit of $\sim$920 sources deg$^{-2}$, on the population of silicate-absorbed ULIRGs at z $\sim$ 1--1.8. This source count is about half of the total number of sources predicted at z $\sim$ 1--2 by various phenomenological models. We note that the high 16/24$μ$m colors measured cannot be reproduced by any of the mid-IR spectral energy distributions assumed by these models, which points to the strong limitations currently affecting our phenomenological and theoretical understanding of infrared galaxy evolution.

preprint2005arXiv

Spitzer 70 and 160-micron Observations of the Extragalactic First Look Survey

We present Spitzer 70um and 160um observations of the Spitzer extragalactic First Look Survey (xFLS). The data reduction techniques and the methods for producing co-added mosaics and source catalogs are discussed. Currently, 26% of the 70um sample and 49% of the 160um-selected sources have redshifts. The majority of sources with redshifts are star-forming galaxies at z<0.5, while about 5% have infrared colors consistent with AGN. The observed infrared colors agree with the spectral energy distribution (SEDs) of local galaxies previously determined from IRAS and ISO data. The average 160um/70um color temperature for the dust is Td~= 30+/-5 K, and the average 70um/24um spectral index is alpha~= 2.4+/-0.4. The observed infrared to radio correlation varies with redshift as expected out to z~1 based on the SEDs of local galaxies. The xFLS number counts at 70um and 160um are consistent within uncertainties with the models of galaxy evolution, but there are indications that the current models may require slight modifications. Deeper 70um observations are needed to constrain the models, and redshifts for the faint sources are required to measure the evolution of the infrared luminosity function.

preprint2005arXiv

The Detection of Crystalline Silicates in Ultra-Luminous Infrared Galaxies

Silicates are an important component of interstellar dust and the structure of these grains -- amorphous versus crystalline -- is sensitive to the local physical conditions. We have studied the infrared spectra of a sample of ultra-luminous infrared galaxies. Here, we report the discovery of weak, narrow absorption features at 11, 16, 19, 23, and 28 microns, characteristic of crystalline silicates, superimposed on the broad absorption bands at 10 and 18 microns due to amorphous silicates in a subset of this sample. These features betray the presence of forsterite (Mg_2SiO_4), the magnesium-rich end member of the olivines. Previously, crystalline silicates have only been observed in circumstellar environments. The derived fraction of forsterite to amorphous silicates is typically 0.1 in these ULIRGs. This is much larger than the upper limit for this ratio in the interstellar medium of the Milky Way, 0.01. These results suggest that the timescale for injection of crystalline silicates into the ISM is short in a merger-driven starburst environment (e.g., as compared to the total time to dissipate the gas), pointing towards massive stars as a prominent source of crystalline silicates. Furthermore, amorphization due to cosmic rays, which is thought to be of prime importance for the local ISM, lags in vigorous starburst environments.

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.

preprint2003arXiv

SINGS: The SIRTF Nearby Galaxies Survey

The SIRTF Nearby Galaxy Survey is a comprehensive infrared imaging and spectroscopic survey of 75 nearby galaxies. Its primary goal is to characterize the infrared emission of galaxies and their principal infrared-emitting components, across a broad range of galaxy properties and star formation environments. SINGS will provide new insights into the physical processes connecting star formation to the interstellar medium properties of galaxies, and provide a vital foundation for understanding infrared observations of the distant universe and ultraluminous and active galaxies. The galaxy sample and observing strategy have been designed to maximize the scientific and archival value of the data set for the SIRTF user community at large. The SIRTF images and spectra will be supplemented by a comprehensive multi-wavelength library of ancillary and complementary observations, including radio continuum, HI, CO, submillimeter, BVRIJHK, H-alpha, Paschen-alpha, ultraviolet, and X-ray data. This paper describes the main astrophysical issues to be addressed by SINGS, the galaxy sample and the observing strategy, and the SIRTF and other ancillary data products.

preprint2001arXiv

Are Starburst Galaxies the Hosts of Gamma-Ray Bursts ?

We present deep 2.2 micron imaging of twelve gamma-ray burst host galaxies. Template spectral energy distributions are fit to the multiband photometry between visible and near-infrared wavelengths to derive a better constraint on the stellar mass of these galaxies. The internal extinction in the host galaxies is estimated using the rest-frame ultraviolet (UV) slope. We find that the extinction corrected star-formation rates (SFRs) of the galaxies are significantly larger than rates derived from emission lines in the visible or the UV continuum. The ratio between the extinction corrected SFRs and stellar mass for 7 of the host galaxies is high compared to local starbursts and 3 of the hosts have derived far-infrared luminosities comparable to infrared luminous galaxies. In addition, existing observational data reveal that at least 6 of the 11 putative hosts seem to be disturbed or have companion galaxies within a projected angular separation of about 2.5 arcsec. If we assume that the host and the companion are at similar redshifts, this corresponds to a physical separation of less than 20 kpc, providing some evidence for an ongoing/recent tidal encounter. We conclude that tidally-induced starbursts such as those found in infrared luminous galaxies might be popular birthplaces for gamma-ray bursts. The age of the stellar population in 4 out of 6 galaxies is rather young, of order 10 Myr. This favors models where gamma-ray bursts are due to the core-collapse of isolated, massive stars and explosion of the resultant black hole-accretion disk system.

preprint2000arXiv

K band Spectroscopy of Ultraluminous Infrared Galaxies: The 2 Jy Sample

We present near-infrared spectroscopy for a complete sample of 33 ultraluminous infrared galaxies at a resolution of R\approx 1000. Most of the wavelength range from 1.80-2.20 microns in the rest frame is covered, including the Pa-alpha and Br-gamma hydrogen recombination lines, and the molecular hydrogen vibration-rotation 1-0 S(1) and S(3) lines. Other species, such as He I, [Fe II], and [Si VI] appear in the spectra as well, in addition to a number of weaker molecular hydrogen lines. Nuclear extractions for each of the individual galaxies are presented here, along with spectra of secondary nuclei, where available. The Pa-alpha emission is seen to be highly concentrated on the nuclei, typically with very little emision extending beyond a radius of 1 kpc. Signatures of active nuclei are rare in the present sample, occurring in only two of the 33 galaxies. It is found that visual extinctions to the nuclei via the Pa-alpha/Br-gamma line ratio in excess of 10 magnitudes are relatively common among ULIRGs, and that visual extinctions greater than 25 mag are necessary to conceal a QSO emitting half the total bolometric luminosity. The vibration-rotation lines of molecular hydrogen appear to be predominantly thermal in origin, with effective temperatures generally around 2200 K. The relative nuclear velocities between double nucleus ULIRGs are investigated, through which it is inferred that the maximum deprojected velocity difference is about 200 km/s. This figure is lower than the velocities predicted by physical models of strong interactions/mergers of large, gas-rich galaxies.

preprint1999arXiv

Near-Infrared Spectra of Ultraluminous Infrared Galaxies

Near infrared spectra with resolution R ~ 1100 in the rest wavelength range 1.8-2.2 microns have been obtained for a complete sample of 33 ultraluminous infrared galaxies. Of the 33 objects observed, 2 show evidence of a central AGN through either a broad Paschen-alpha line or emission in the 1.963 micron fine structure of [Si VI]. In the median spectrum of the remaining 31 objects, the lines present are recombination lines of Hydrogen, neutral Helium, vibration-rotation lines of H_2, and [Fe II]. There is no indication of AGN activity in the median spectrum, either through broad atomic recombination lines or high ionization lines. No trends in luminosity are apparent when subsets of the 31 non-AGN ULIRGs are binned by luminosity and median combined. When secondary nuclei exist in ULIRGs, they typically have spectra very much like those seen in the primary nuclei.

preprint1997arXiv

The Unusual Near-Infrared Morphology of the Radio Loud Quasar 4C+09.17

Near-infrared images of the luminous, high redshift (z=2.1108) radio loud quasar 4C+09.17 reveal a complex structure. The quasar (K=15.76 mag) is surrounded by three "companion" objects having 17.9< K < 20.2 mag at radii of 1.7" < r < 2.9", as well as bright, diffuse emission. The brightest companion has a redshift of z=0.8384 (Lehnert & Becker 1997) and its optical-infrared colors (Lehnert et al. 1997) are consistent with a late-type spiral galaxy at this redshift with a luminosity of about 2L*. This object is likely the galaxy responsible for the strongest MgII absorption line system seen in the spectrum of 4C+09.17 by Barthel et al. (1990). Redshifts are not available for the remaining two companions. The red colors of the second brightest companion appear most consistent with a high redshift star-forming galaxy at z > 1.5. If this object is at the redshift of 4C +09.17 it has a luminosity of about 7L*. The faintest companion has colors which are unlike those expected from either a spiral or an E/S0 galaxy at any redshift associated with the 4C+09.17 system. Since this object lies along the same direction as the radio jet/lobe of 4C+09.17, as well as the extended Ly-alpha emission mapped by Heckman et al. (1991) we suggest that this component can be explained as a combination of strong line emission and scattered QSO light. The resolved, diffuse emission surrounding 4C+09.17 is bright, K~17.0 mag, and about one magnitude redder in J-K than the quasar. If this emission is starlight, a very luminous elliptical host galaxy is implied for 4C+09.17. Scattered and reddened AGN light, emission line gas, and flux from absorbing galaxies along the line of sight may all contribute to this emission.