Source author record

D. A. Dale

D. A. Dale 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

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

52 published item(s)

preprint2026arXiv

PHANGS-JWST: the largest extragalactic molecular cloud catalog traced by polycyclic aromatic hydrocarbon emission

High-resolution JWST images of nearby spiral galaxies reveal polycyclic aromatic hydrocarbon (PAH) emission structures that trace molecular gas, including CO-dark regions. We identify ISM cloud structures in PHANGS-JWST 7.7 $μ$m PAH maps for 66 galaxies, smoothed to 30 pc and at native resolution, extracting 108,466 and 146,040 clouds, respectively. Molecular properties were inferred using a linear conversion from PAH to CO. Given the tendency for clouds in galaxy centers to overlap in velocity space, we opted to flag these and omit them from the analysis in this work. The remaining clouds correspond to giant molecular clouds, such as those detected in CO(2-1) emission by ALMA, or lower surface density clouds that either fall below the ALMA detection limits of existing maps or genuinely have no molecular counterpart. Cross-matching with ALMA CO maps at 90 pc in 27 galaxies shows that 41 % of PAH clouds have CO associations. The converted molecular properties vary little across environments, but the most massive clouds are preferentially found in spiral arms. Fitting lognormal mass distributions down to $2\times10^{3} M_{\odot}$ shows that spiral arms host the highest-mass clouds, consistent with enhanced formation in arm gravitational potentials. Cloud molecular surface densities decline by a factor of $\sim 1.5-2$ toward $2 - 3 R_{e}$. However, the trend largely varies in individual galaxies, with flat, decreasing, and even no trend as a function of galactocentric radius. Factors like large-scale processes and morphologies might influence the observed trends. We publish two catalogs online, one at the common resolution of 30 pc and another at the native resolution. We expect them to have broad utility for future PAH clouds, molecular clouds, and star formation studies.

preprint2022arXiv

Fraction of Stars in Clusters for the LEGUS Dwarf Galaxies

We study the young star cluster populations in 23 dwarf and irregular galaxies observed by the HST Legacy ExtraGalactic Ultraviolet Survey (LEGUS), and examine relationships between the ensemble properties of the cluster populations and those of their host galaxies: star formation rate (SFR) density ($Σ_{\rm SFR}$). A strength of this analysis is the availability of SFRs measured from temporally resolved star formation histories which provide the means to match cluster and host-galaxy properties on several timescales (1-10, 1-100, and 10-100~Myr). Nevertheless, studies of this kind are challenging for dwarf galaxies due to the small numbers of clusters in each system. We mitigate these issues by combining the clusters across different galaxies with similar $Σ_{\rm SFR}$ properties. We find good agreement with a well-established relationship ($M_{V}^{brightest}$-SFR), but find no significant correlations between $Σ_{\rm SFR}$ and the slopes of the cluster luminosity function, mass function, nor the age distribution. We also find no significant trend between the the fraction of stars in bound clusters at different age ranges ($Γ_{1-10}$, $Γ_{10-100}$, and $Γ_{1-100}$) and $Σ_{\rm SFR}$ of the host galaxy. Our data show a decrease in $Γ$ over time (from 1-10 to 10-100~Myr) suggesting early cluster dissolution, though the presence of unbound clusters in the youngest time bin makes it difficult to quantify the degree of dissolution. While our data do not exhibit strong correlations between $Σ_{\rm SFR}$ and ensemble cluster properties, we cannot rule out that a weak trend might exist given the relatively large uncertainties due to low number statistics and the limited $Σ_{\rm SFR}$ range probed.

preprint2022arXiv

Linking stellar populations to HII regions across nearby galaxies: I. Constraining pre-supernova feedback from young clusters in NGC1672

One of the fundamental factors regulating the evolution of galaxies is stellar feedback. However, we still do not have strong observational constraints on the relative importance of the different feedback mechanisms (e.g. radiation, ionised gas pressure, stellar winds) in driving HII region evolution and molecular cloud disruption. In this letter, we constrain the relative importance of the various feedback mechanisms from young massive star populations by resolving HII regions across the disk of the nearby star-forming galaxy NGC 1672. We combine measurements of ionised gas nebular lines obtained by PHANGS-MUSE, with high-resolution imaging from the HST in both the narrow-band Hα and broad-band filters. We identify a sample of 40 isolated, compact HII regions in the HST Hα image, for which we measure the sizes that were previously unresolved in seeing-limited ground-based observations. Additionally, we identify the ionisation source(s) for each HII region from catalogues produced as part of the PHANGS-HST survey. We find that the HII regions investigated are mildly dominated by thermal or wind pressure, yet their elevation above the radiation pressure is within the expected uncertainty range. We see that radiation pressure provides a substantially higher contribution to the total pressure than previously found in the literature over similar size scales. In general, we find higher pressures within more compact HII regions, which is driven by the inherent size scaling relations of each pressure term, albeit with significant scatter introduced by the variation in the stellar population properties (e.g. luminosity, mass, age, metallicity). For nearby galaxies, here we provide a promising approach that could yield the statistics required to map out how the importance of different stellar feedback mechanisms evolve over the lifetime of an HII region.

preprint2022arXiv

Monitoring AGNs with H$β$ Asymmetry. III. Long-term Reverberation Mapping Results of 15 Palomar-Green Quasars

In this third paper of the series reporting on the reverberation mapping (RM) campaign of active galactic nuclei with asymmetric H$β$ emission-line profiles, we present results for 15 Palomar-Green (PG) quasars using spectra obtained between the end of 2016 to May 2021. This campaign combines long time spans with relatively high cadence. For 8 objects, both the time lags obtained from the entire light curves and the measurements from individual observing seasons are provided. Reverberation mapping of 9 of our targets has been attempted for the first time, while the results for 6 others can be compared with previous campaigns. We measure the H$β$ time lags over periods of years and estimate their black hole masses. The long duration of the campaign enables us to investigate their broad line region (BLR) geometry and kinematics for different years by using velocity-resolved lags, which demonstrate signatures of diverse BLR geometry and kinematics. The BLR geometry and kinematics of individual objects are discussed. In this sample, the BLR kinematics of Keplerian/virialized motion and inflow is more common than outflow.

preprint2022arXiv

Star Cluster Formation and Evolution in M101: An Investigation with the Legacy Extragalactic UV Survey

We present Hubble Space Telescope WFC3/UVIS (F275W, F336W) and ACS/WFC optical (F435W, F555W, and F814W) observations of the nearby grand-design spiral galaxy M101 as part of the Legacy Extragalactic UV Survey (LEGUS). Compact sources detected in at least four bands were classified by both human experts and the convolutional neural network StarcNet. Human experts classified the 2,351 brightest sources, retrieving $N_{c} = 965$ star clusters. StarcNet, trained on LEGUS data not including M101, classified all 4,725 sources detected in four bands, retrieving $N_{c} = 2,270$ star clusters. The combined catalog represents the most complete census to date of compact star clusters in M101. We find that for the 2,351 sources with both a visual- and ML-classification StarcNet is able to reproduce the human classifications at high levels of accuracy ($\sim 80-90\%$), which is equivalent to the level of agreement between human classifiers in LEGUS. The derived cluster age distribution implies a disruption rate of $dN/dτ\propto τ^{-0.45 \pm 0.14}$ over $10^{7} < τ< 10^{8.5}$yr for cluster masses $\geq 10^{3.55} M_{\odot}$ for the central region of M101 and $dN/dτ\propto τ^{-0.02 \pm 0.15}$ for cluster masses $\geq 10^{3.38} M_{\odot}$ in the northwest region of the galaxy. The trends we recover are weaker than those of other nearby spirals (e.g. M51) and starbursts, consistent with the M101 environment having a lower-density interstellar medium, and providing evidence in favor of environmentally-dependent cluster disruption in the central, southeast, and northwest regions of M101.

preprint2020arXiv

The Age-Dependence of Mid-Infrared Emission Around Young Star Clusters

Using the star cluster catalogs from the Hubble Space Telescope program Legacy ExtraGalactic UV survey (LEGUS) and 8 $μ$m images from the IRAC camera on the Spitzer Space Telescope for 5 galaxies within 5 Mpc, we investigate how the 8 $μ$m dust luminosity correlates with the stellar age on the 30--50 pc scale of star forming regions. We construct a sample of 97 regions centered at local peaks of 8 $μ$m emission, each containing one or more young star cluster candidates from the LEGUS catalogs. We find a tight anti-correlation with a Pearson correlation coefficient of $r=-0.84\pm0.05$ between the mass-normalized dust-only 8 $μ$m luminosity and the age of stellar clusters younger than 1 Gyr; the 8 $μ$m luminosity decreases with increasing age of the stellar population. Simple assumptions on a combination of stellar and dust emission models reproduce the observed trend. We also explore how the scatter of the observed trend depends on assumptions of stellar metallicity, PAH abundance, fraction of stellar light absorbed by dust, and instantaneous versus continuous star formation models. We find that variations in stellar metallicity have little effect on the scatter, while PAH abundance and the fraction of dust-absorbed light bracket the full range of the data. We also find that the trend is better explained by continuous star formation, rather than instantaneous burst models. We ascribe this result to the presence of multiple star clusters with different ages in many of the regions. Upper limits of the dust-only 8 $μ$m emission as a function of age are provided.

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

SOFIA/HAWC+ traces the magnetic fields in NGC 1068

We report the first detection of galactic spiral structure by means of thermal emission from magnetically aligned dust grains. Our 89 $μ$m polarimetric imaging of NGC 1068 with the High-resolution Airborne Wideband Camera/Polarimeter (HAWC+) on NASA's Stratospheric Observatory for Infrared Astronomy (SOFIA) also sheds light on magnetic field structure in the vicinity of the galaxy's inner-bar and active galactic nucleus (AGN). We find correlations between the 89 $μ$m magnetic field vectors and other tracers of spiral arms, and a symmetric polarization pattern as a function of the azimuthal angle arising from the projection and inclination of the disk field component in the plane of the sky. The observations can be fit with a logarithmic spiral model with pitch angle of $16.9^{+2.7}_{-2.8}$$^{\circ}$ and a disk inclination of $48\pm2^{\circ}$. We infer that the bulk of the interstellar medium from which the polarized dust emission originates is threaded by a magnetic field that closely follows the spiral arms. Inside the central starburst disk ($<1.6$ kpc), the degree of polarization is found to be lower than for far-infrared sources in the Milky Way, and has minima at the locations of most intense star formation near the outer ends of the inner-bar. Inside the starburst ring, the field direction deviates from the model, becoming more radial along the leading edges of the inner-bar. The polarized flux and dust temperature peak $\sim 3-6$" NE of the AGN at the location of a bow shock between the AGN outflow and the surrounding interstellar medium, but the AGN itself is weakly polarized ($< 1$%) at both 53 and 89 \um.

preprint2019arXiv

Spatial segregation of massive clusters in dwarf galaxies

The relative average minimum projected separations of star clusters in the Legacy ExtraGalactic UV Survey (LEGUS) and in tidal dwarfs around the interacting galaxy NGC 5291 are determined as a function of cluster mass to look for cluster-cluster mass segregation. Class 2 and 3 LEGUS clusters, which have a more irregular internal structure than the compact and symmetric class 1 clusters, are found to be mass segregated in low mass galaxies, which means that the more massive clusters are systematically bunched together compared to the lower mass clusters. This mass segregation is not present in high-mass galaxies nor for class 1 clusters. We consider possible causes for this segregation including differences in cluster formation and scattering in the shallow gravitational potentials of low mass galaxies.

preprint2019arXiv

Star formation histories of the LEGUS dwarf galaxies (III): the non-bursty nature of 23 star forming dwarf galaxies

We derive the recent star formation histories of 23 active dwarf galaxies using HST observations from the Legacy Extragalactic UV Survey (LEGUS). We apply a color-magnitude diagram fitting technique using two independent sets of stellar models, PARSEC-COLIBRI and MIST. Despite the non-negligible recent activity, none of the 23 star forming dwarfs show enhancements in the last 100 Myr larger than three times the 100-Myr-average. The unweighted mean of the individual SFHs in the last 100 Myr is also consistent with a rather constant activity, irrespective of the atomic gas fraction. We confirm previous results that for dwarf galaxies the CMD-based average star formation rates (SFRs) are generally higher than the FUV-based SFR. For half of the sample, the 60-Myr-average CMD-based SFR is more than two times the FUV SFR. In contrast, we find remarkable agreement between the 10-Myr-average CMD-based SFR and the H$α$-based SFR. Finally, using core helium burning stars of intermediate mass we study the pattern of star formation spatial progression over the past 60 Myr, and speculate on the possible triggers and connections of the star formation activity with the environment in which these galaxies live. Approximately half of our galaxies show spatial progression of star formation in the last 60 Myr, and/or very recent diffuse and off-center activity compared to RGB stars.

preprint2019arXiv

The WISE Extended Source Catalogue (WXSC) I: The 100 Largest Galaxies

We present mid-infrared photometry and measured global properties of the 100 largest galaxies in the sky, including the Magellanic Clouds, Local Group galaxies M31 and M33, the Fornax and Virgo Galaxy Cluster giants, and many of the most spectacular Messier objects (e.g., M51 and M83). This is the first release of a larger catalog of extended sources as imaged in the mid-infrared, called the WISE Extended Source Catalogue (WXSC). In this study we measure their global attributes, including integrated flux, surface brightness and radial distribution. The largest of the large are the LMC, SMC and the Andromeda Galaxy, which are also the brightest mid-infrared galaxies in the sky. We interrogate the large galaxies using WISE colors, which serve as proxies for four general types of galaxies: bulge-dominated spheroidals, intermediate semi-quiescent disks, star-forming spirals, and AGN-dominated. The colors reveal a tight "sequence" that spans 5 magnitudes in W2-W3 color, ranging from early to late-types, and low to high star-forming activity; we fit the functional form given by: ${\rm (W1-W2)} = [0.015 \times {\rm e}^{ \frac{\rm (W2-W3)}{1.38} }] - 0.08$. Departures from this sequence may reveal nuclear, starburst, and merging events. Physical properties and luminosity attributes are computed, notably the diameter, aggregate stellar mass and the dust-obscured star formation activity. We introduce the 'pinwheel' diagram which depicts physical properties with respect to the median value observed for WISE galaxies in the local universe. Utilized with the WXSC, this diagram will delineate between different kinds of galaxies, identifying those with similar star formation and structural properties. Finally, we present the mid-infrared photometry of the 25 brightest globular clusters in the sky, including Omega Centauri, 47 Tucanae and a number of famed night-sky targets (e.g. M 13). (Abridged)

preprint2016arXiv

Baryonic Distributions in Galaxy Dark Matter Haloes I: New Observations of Neutral and Ionized Gas Kinematics

We present a combination of new and archival neutral hydrogen (HI) observations and new ionized gas spectroscopic observations for sixteen galaxies in the statistically representative EDGES kinematic sample. HI rotation curves are derived from new and archival radio synthesis observations from the Very Large Array (VLA) as well as processed data products from the Westerbork Radio Synthesis Telescope (WSRT). The HI rotation curves are supplemented with optical spectroscopic integral field unit (IFU) observations using SparsePak on the WIYN 3.5 m telescope to constrain the central ionized gas kinematics in twelve galaxies. The full rotation curves of each galaxy are decomposed into baryonic and dark matter halo components using 3.6$μ$m images from the Spitzer Space Telescope for the stellar content, the neutral hydrogen data for the atomic gas component, and, when available, CO data from the literature for the molecular gas component. Differences in the inferred distribution of mass are illustrated under fixed stellar mass-to-light ratio (M/L) and maximum disc/bulge assumptions in the rotation curve decomposition.

preprint2016arXiv

The properties, origin and evolution of stellar clusters in galaxy simulations and observations

We investigate the properties and evolution of star particles in two simulations of isolated spiral galaxies, and two galaxies from cosmological simulations. Unlike previous numerical work, where typically each star particle represents one `cluster', for the isolated galaxies we are able to model features we term `clusters' with groups of particles. We compute the spatial distribution of stars with different ages, and cluster mass distributions, comparing our findings with observations including the recent LEGUS survey. We find that spiral structure tends to be present in older (100s Myrs) stars and clusters in the simulations compared to the observations. This likely reflects differences in the numbers of stars or clusters, the strength of spiral arms, and whether the clusters are allowed to evolve. Where we model clusters with multiple particles, we are able to study their evolution. The evolution of simulated clusters tends to follow that of their natal gas clouds. Massive, dense, long-lived clouds host massive clusters, whilst short-lived clouds host smaller clusters which readily disperse. Most clusters appear to disperse fairly quickly, in basic agreement with observational findings. We note that embedded clusters may be less inclined to disperse in simulations in a galactic environment with continuous accretion of gas onto the clouds than isolated clouds and correspondingly, massive young clusters which are no longer associated with gas tend not to occur in the simulations. Caveats of our models include that the cluster densities are lower than realistic clusters, and the simplistic implementation of stellar feedback.

preprint2015arXiv

Baryonic Distributions in the Dark Matter Halo of NGC5005

We present results from multiwavelength observations of the galaxy NGC5005. We use new neutral hydrogen (HI) observations from the Very Large Array to examine the neutral gas morphology and kinematics. We find an HI disk with a well-behaved flat rotation curve in the radial range 20\arcsec-140\arcsec. Ionized gas observations from the SparsePak integral field unit on the WIYN 3.5m telescope provide kinematics for the central 70\arcsec. We use both the SparsePak and HI velocity fields to derive a rotation curve for NGC5005. Deep 3.6μm observations from the Spitzer Space Telescope probe the faint extended stellar population of NGC5005. The images reveal a large stellar disk with a high surface brightness component that transitions to a low surface brightness component at a radius nearly 1.6 times farther than the extent of the gas disk detected in HI. The 3.6μm image is also decomposed into bulge and disk components to account for the stellar light distribution. Optical broadband B and R and narrowband Hα from the WIYN 0.9m telescope complement the 3.6μm data by providing information about the dominant stellar population and current star formation activity. The neutral and ionized gas rotation curve is used along with the stellar bulge and disk light profiles to decompose the mass distributions in NGC5005 and determine a dark matter halo model. The maximum stellar disk contribution to the total rotation curve is only about 70\%, suggesting that dark matter makes a significant contribution to the dynamics at all radii.

preprint2015arXiv

The Brightest Young Star Clusters in NGC 5253

The nearby dwarf starburst galaxy NGC5253 hosts a number of young, massive star clusters, the two youngest of which are centrally concentrated and surrounded by thermal radio emission (the `radio nebula'). To investigate the role of these clusters in the starburst energetics, we combine new and archival Hubble Space Telescope images of NGC5253 with wavelength coverage from 1500 Ang to 1.9 micron in 13 filters. These include H-alpha, P-beta, and P-alpha, and the imaging from the Hubble Treasury Program LEGUS (Legacy Extragalactic UV Survey). The extraordinarily well-sampled spectral energy distributions enable modeling with unprecedented accuracy the ages, masses, and extinctions of the 9 optically brightest clusters (M_V < -8.8) and the two young radio nebula clusters. The clusters have ages ~1-15 Myr and masses ~1x10^4 - 2.5x10^5 M_sun. The clusters' spatial location and ages indicate that star formation has become more concentrated towards the radio nebula over the last ~15 Myr. The most massive cluster is in the radio nebula; with a mass 2.5x10^5 M_sun and an age ~1 Myr, it is 2-4 times less massive and younger than previously estimated. It is within a dust cloud with A_V~50 mag, and shows a clear nearIR excess, likely from hot dust. The second radio nebula cluster is also ~1 Myr old, confirming the extreme youth of the starburst region. These two clusters account for about half of the ionizing photon rate in the radio nebula, and will eventually supply about 2/3 of the mechanical energy in present-day shocks. Additional sources are required to supply the remaining ionizing radiation, and may include very massive stars.

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

Cool dust heating and temperature mixing in nearby star-forming galaxies

Physical conditions of the interstellar medium in galaxies are closely linked to the ambient radiation field and the heating of dust grains. In order to characterize dust properties in galaxies over a wide range of physical conditions, we present here the radial surface brightness profiles of the entire sample of 61 galaxies from Key Insights into Nearby Galaxies: Far-Infrared Survey with Herschel (KINGFISH). The main goal of our work is the characterization of the grain emissivities, dust temperatures, and interstellar radiation fields responsible for heating the dust. After fitting the dust and stellar radial profiles with exponential functions, we fit the far-infrared spectral energy distribution (SED) in each annular region with single-temperature modified black bodies using both variable (MBBV) and fixed (MBBF) emissivity indices beta, as well as with physically motivated dust models. Results show that while most SED parameters decrease with radius, the emissivity index beta also decreases with radius in some galaxies, but in others is increasing, or rising in the inner regions and falling in the outer ones. Despite the fixed grain emissivity (average beta~ 2.1) of the physically-motivated models, they are well able to accommodate flat spectral slopes with beta<= 1. We find that flatter slopes (beta<= 1.5) are associated with cooler temperatures, contrary to what would be expected from the usual Tdust-beta degeneracy. This trend is related to variations in Umin since beta and Umin are very closely linked over the entire range in Umin sampled by the KINGFISH galaxies: low Umin is associated with flat beta<=1. Both these results strongly suggest that the low apparent βvalues (flat slopes) in MBBV fits are caused by temperature mixing along the line-of-sight, rather than by intrinsic variations in grain properties. Abstract truncated for arXiv.

preprint2013arXiv

Calibration of the total infrared luminosity of nearby galaxies from Spitzer and Herschel bands

We present new empirical calibrations to estimate resolved and integrated total infrared luminosities from Spitzer and Herschel bands used as monochromatic or combined tracers. We base our calibrations on resolved elements of nearby galaxies (3 to 30 Mpc) observed with Herschel. We perform a resolved SED modelling of these objects using the Draine and Li (2007) dust models and investigate the influence of the addition of SPIRE measurements in the estimation of LTIR. We find that using data up to 250 um leads to local LTIR values consistent with those obtained with a complete coverage (up to 500 um) within 10 per cent for most of our resolved elements. We then study the distribution of energy in the resolved SEDs of our galaxies. The bulk of energy (30-50 per cent) is contained in the (70-160 um) band. The (24-70 um) fraction decreases with increasing metallicity. The (160-1100 um) submillimeter band can account for up to 25 per cent of the LTIR in metal-rich galaxies. We investigate the correlation between TIR surface brightnesses/luminosities and monochromatic Spitzer and Herschel surface brightnesses/luminosities. The three PACS bands can be used as reliable monochromatic estimators of the LTIR, the 100 um band being the most reliable monochromatic tracer. There is also a strong correlation between the SPIRE 250um and LTIR, although with more scatter than for the PACS relations. We also study the ability of our monochromatic relations to reproduce integrated LTIR of nearby galaxies as well as LTIR of z=1-3 sources. Finally, we provide calibration coefficients that can be used to derive TIR surface brightnesses/luminosities from a combination of Spitzer and Herschel surface brightnesses/fluxes and analyse the associated uncertainties.

preprint2013arXiv

Mapping dust through emission and absorption in nearby galaxies

Dust has long been identified as a barrier to measuring inherent galaxy properties. However, the link between dust and attenuation is not straightforward and depends on both the amount of dust and its distribution. Herschel imaging of nearby galaxies undertaken as part of the KINGFISH project allows us to map the dust as seen in emission with unprecedented sensitivity and ~1 kpc resolution. We present here new optical integral field unit spectroscopy for eight of these galaxies that provides complementary 100-200 pc scale maps of the dust attenuation through observation of the reddening in both the Balmer decrement and the stellar continuum. The stellar continuum reddening, which is systematically less than that observed in the Balmer decrement, shows no clear correlation with the dust, suggesting that the distribution of stellar reddening acts as a poor tracer of the overall dust content. The brightest HII regions are observed to be preferentially located in dusty regions, and we do find a correlation between the Balmer line reddening and the dust mass surface density for which we provide an empirical relation. Some of the high-inclination systems in our sample exhibit high extinction, but we also find evidence that unresolved variations in the dust distribution on scales smaller than 500 pc may contribute to the scatter in this relation. We caution against the use of integrated A_V measures to infer global dust properties.

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

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

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.

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

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

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

preprint2011arXiv

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

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

preprint2011arXiv

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.

preprint2010arXiv

An Aromatic Inventory of the Local Volume

Using infrared photometry from the Spitzer Space Telescope, we perform the first inventory of aromatic feature emission (AFE, but also commonly referred to as PAH emission) for a statistically complete sample of star-forming galaxies in the local volume. The photometric methodology involved is calibrated and demonstrated to recover the aromatic fraction of the IRAC 8 micron flux with a standard deviation of 6% for a training set of 40 SINGS galaxies (ranging from stellar to dust dominated) with both suitable mid-infrared Spitzer IRS spectra and equivalent photometry. A potential factor of two improvement could be realized with suitable 5.5 and 10 micron photometry, such as what may be provided in the future by JWST. The resulting technique is then applied to mid-infrared photometry for the 258 galaxies from the Local Volume Legacy (LVL) survey, a large sample dominated in number by low-luminosity dwarf galaxies for which obtaining comparable mid-infrared spectroscopy is not feasible. We find the total LVL luminosity due to five strong aromatic features in the 8 micron complex to be 2.47E10 solar luminosities with a mean volume density of 8.8E6 solar luminosities per cubic Megaparsec. Twenty-four of the LVL galaxies, corresponding to a luminosity cut at M = -18.22 in the B band, account for 90% of the aromatic luminosity. Using oxygen abundances compiled from the literature for 129 of the 258 LVL galaxies, we find a correlation between metallicity and the aromatic to total infrared emission ratio but not the aromatic to total 8 micron dust emission ratio. A possible explanation is that metallicity plays a role in the abundance of aromatic molecules relative to the total dust content, but other factors such as star formation and/or the local radiation field affect the excitation of those molecules.

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

Nearby Galaxies with Spitzer

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

preprint2010arXiv

The Calibration of Monochromatic Far-Infrared Star Formation Rate Indicators

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

preprint2009arXiv

Radial distribution of stars, gas and dust in SINGS galaxies. I. Surface photometry and morphology

We present ultraviolet through far-infrared surface brightness profiles for the 75 galaxies in the Spitzer Infrared Nearby Galaxies Survey (SINGS). The imagery used to measure the profiles includes GALEX UV data, optical images from KPNO, CTIO and SDSS, near-IR data from 2MASS, and mid- and far-infrared images from Spitzer. Along with the radial profiles, we also provide multi-wavelength asymptotic magnitudes and several non-parametric indicators of galaxy morphology: the concentration index (C_42), the asymmetry (A), the Gini coefficient (G) and the normalized second-order moment of the brightest 20% of the galaxy's flux (M_20). Our radial profiles show a wide range of morphologies and multiple components (bulges, exponential disks, inner and outer disk truncations, etc.) that vary not only from galaxy to galaxy but also with wavelength for a given object. In the optical and near-IR, the SINGS galaxies occupy the same regions in the C_42-A-G-M_20 parameter space as other normal galaxies in previous studies. However, they appear much less centrally concentrated, more asymmetric and with larger values of G when viewed in the UV (due to star-forming clumps scattered across the disk) and in the mid-IR (due to the emission of Polycyclic Aromatic Hydrocarbons at 8.0 microns and very hot dust at 24 microns).

preprint2009arXiv

Radial distribution of stars, gas and dust in SINGS galaxies. II. Derived dust properties

We present a detailed analysis of the radial distribution of dust properties in the SINGS sample, performed on a set of UV, IR and HI surface brightness profiles, combined with published molecular gas profiles and metallicity gradients. The internal extinction, derived from the TIR-to-FUV luminosity ratio, decreases with radius, and is larger in Sb-Sbc galaxies. The TIR-to-FUV ratio correlates with the UV spectral slope beta, following a sequence shifted to redder UV colors with respect to that of starbursts. The star formation history (SFH) is identified as the main driver of this departure. We have also derived radial profiles of the total dust mass surface density, the fraction of the dust mass contributed by PAHs, the fraction of the dust mass heated by very intense starlight and the intensity of the radiation field heating the grains. The dust profiles are exponential, their radial scale-length being constant from Sb to Sd galaxies (only ~10% larger than the stellar scale-length). Many S0/a-Sab galaxies have central depressions in their dust radial distributions. The PAH abundance increases with metallicity for 12+\log(O/H)<9, and at larger metallicities the trend flattens and even reverses, with the SFH being a plausible underlying driver for this behavior. The dust-to-gas ratio is also well correlated with metallicity and therefore decreases with galactocentric radius.

preprint2009arXiv

The Spitzer Local Volume Legacy: Survey Description and Infrared Photometry

The survey description and the near-, mid-, and far-infrared flux properties are presented for the 258 galaxies in the Local Volume Legacy (LVL). LVL is a Spitzer Space Telescope legacy program that surveys the local universe out to 11 Mpc, built upon a foundation of ultraviolet, H-alpha, and HST imaging from 11HUGS (11 Mpc H-alpha and Ultraviolet Galaxy Survey) and ANGST (ACS Nearby Galaxy Survey Treasury). LVL covers an unbiased, representative, and statistically robust sample of nearby star-forming galaxies, exploiting the highest extragalactic spatial resolution achievable with Spitzer. As a result of its approximately volume-limited nature, LVL augments previous Spitzer observations of present-day galaxies with improved sampling of the low-luminosity galaxy population. The collection of LVL galaxies shows a large spread in mid-infrared colors, likely due to the conspicuous deficiency of 8um PAH emission from low-metallicity, low-luminosity galaxies. Conversely, the far-infrared emission tightly tracks the total infrared emission, with a dispersion in their flux ratio of only 0.1 dex. In terms of the relation between infrared-to-ultraviolet ratio and ultraviolet spectral slope, the LVL sample shows redder colors and/or lower infrared-to-ultraviolet ratios than starburst galaxies, suggesting that reprocessing by dust is less important in the lower mass systems that dominate the LVL sample. Comparisons with theoretical models suggest that the amplitude of deviations from the relation found for starburst galaxies correlates with the age of the stellar populations that dominate the ultraviolet/optical luminosities.

preprint2008arXiv

The relations among 8, 24, and 160 micron dust emission within nearby spiral galaxies

We investigate the relations among the stellar continuum-subtracted 8 micron polycyclic aromatic hydrocarbon (PAH 8 micron) emission, 24 micron hot dust emission, and 160 micron cold dust emission in fifteen nearby face-on spiral galaxies in the Spitzer Infrared Nearby Galaxies Survey sample. The relation between PAH 8 and 24 micron emission measured in ~2 kpc regions is found to exhibit a significant amount of scatter, and strong spatial variations are observed in the (PAH 8 micron)/24 micron surface brightness ratio. In particular, the (PAH 8 micron)/24 micron surface brightness ratio is observed to be high in the diffuse interstellar medium and low in bright star-forming regions and other locations with high 24 micron surface brightness. PAH 8 micron emission is found to be well-correlated with 160 micron emission on spatial scales of ~2 kpc, and the (PAH 8 micron)/160 micron surface brightness ratio is generally observed to increase as the 160 micron surface brightness increases. These results suggest that the PAHs are associated with the diffuse, cold dust that produces most of the 160 micron emission in these galaxies, and the variations in the (PAH 8 micron)/160 micron ratio may generally be indicative of either the intensity or the spectrum of the interstellar radiation field that is heating both the PAHs and the diffuse interstellar dust.

preprint2008arXiv

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

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

preprint2007arXiv

The Calibration of Mid-Infrared Star Formation Rate Indicators

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

preprint2007arXiv

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

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

preprint2007arXiv

Warm molecular hydrogen in the Spitzer SINGS galaxy sample

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

preprint2006arXiv

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

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

preprint2006arXiv

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

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

preprint2006arXiv

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

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

preprint2004arXiv

IRAC Observations of M81

IRAC images of M81 show three distinct morphological constituents: a smooth distribution of evolved stars with bulge, disk, and spiral arm components; a clumpy distribution of dust emission tracing the spiral arms; and a pointlike nuclear source. The bulge stellar colors are consistent with M-type giants, and the disk colors are consistent with a slightly younger population. The dust emission generally follows the blue and ultraviolet emission, but there are large areas that have dust emission without ultraviolet and smaller areas with ultraviolet but little dust emission. The former are presumably caused by extinction, and the latter may be due to cavities in the gas and dust created by supernova explosions. The nucleus appears fainter at 8 um than expected from ground-based 10 um observations made four years ago.

preprint2004arXiv

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

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

preprint2004arXiv

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

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

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.