Source author record

Daniel A. Dale

Daniel 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

55works
7topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

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

Building this map preview

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

Published work

55 published item(s)

preprint2025arXiv

Evolved Supergiants in PHANGS I: Red Supergiants in 19 Galaxies between 5-20 Mpc with HST and JWST

Red supergiants (RSGs) are important for our understanding of supernova progenitors, stellar populations, stellar evolution, mass loss and dust production. Extragalactic surveys of RSGs have a long history in the Local Group, but few studies exist beyond that due to the limited resolution and sensitivity of ground-based and previous space-based infrared observatories. Here we demonstrate the combined power of HST and JWST to push systematic searches of RSGs out to $\sim$20 Mpc. We introduce a catalog of 97057 RSGs -- the largest single-survey release of RSGs -- with masses $\gtrsim$10 M$_{\odot}$ in 19 galaxies from the PHANGS HST+JWST Treasury program. We use HST F814W and JWST F200W photometry to select stars as RSGs based on predicted colors and magnitudes from PARSEC isochrones. The spatial distribution of our recovered RSGs follow the familiar pattern of mostly being concentrated in active star-forming regions such as spiral arms and central starburst rings. The RSG number density on kpc-scales is strongly correlated ($r_s$$\sim$0.82) with local star-formation rate density ($Σ_{SFR}$) traced by extinction-corrected far-ultraviolet (FUV) from GALEX+WISE, and weakly correlated ($r_s$$\sim$0.57) with the total stellar mass density ($Σ_*$), traced by near-infrared emission from WISE+Spitzer. The number of RSGs per mass of stellar populations with ages 6-30 Myr (the likely age range of RSGs $>$10 M$_{\odot}$) is $\sim$1 per 10$^{3.77\pm0.27}$ M$_{\odot}$, assuming constant star-formation rates from FUV+W4. Our sample will be a useful resource for tracking progenitors and feedback sites of future supernovae in PHANGS, age-dating stellar populations, and more.

preprint2024arXiv

ALMA-LEGUS I: The Influence of Galaxy Morphology on Molecular Cloud Properties

We present a comparative study of the molecular gas in two galaxies from the LEGUS sample: barred spiral NGC 1313 and flocculent spiral NGC 7793. These two galaxies have similar masses, metallicities, and star formation rates, but NGC 1313 is forming significantly more massive star clusters than NGC 7793, especially young massive clusters (<10 Myr, >10^4 Msol). Using ALMA CO(2-1) observations of the two galaxies with the same sensitivities and resolutions of 13 pc, we directly compare the molecular gas in these two similar galaxies to determine the physical conditions responsible for their large disparity in cluster formation. By fitting size-linewidth relations for the clouds in each galaxy, we find that NGC 1313 has a higher intercept than NGC 7793, implying that its clouds have higher kinetic energies at a given size scale. NGC 1313 also has more clouds near virial equilibrium than NGC 7793, which may be connected to its higher rate of massive cluster formation. However, these virially bound clouds do not show a stronger correlation with young clusters than that of the general cloud population. We find surprisingly small differences between the distributions of molecular cloud populations in the two galaxies, though the largest of those differences are that NGC 1313 has higher surface densities and lower free-fall times.

preprint2024arXiv

ALMA-LEGUS II: The Influence of Sub-Galactic Environment on Molecular Cloud Properties

We compare the molecular cloud properties in sub-galactic regions of two galaxies, barred spiral NGC 1313, which is forming many massive clusters, and flocculent spiral NGC 7793, which is forming significantly fewer massive clusters despite having a similar star formation rate to NGC 1313. We find that there are larger variations in cloud properties between different regions within each galaxy than there are between the galaxies on a global scale, especially for NGC 1313. There are higher masses, linewidths, pressures, and virial parameters in the arms of NGC 1313 and center of NGC 7793 than in the interarm and outer regions of the galaxies. The massive cluster formation of NGC 1313 may be driven by its greater variation in environments, allowing more clouds with the necessary conditions to arise, although no one parameter seems primarily responsible for the difference in star formation. Meanwhile NGC 7793 has clouds that are as massive and have as much kinetic energy as clouds in the arms of NGC 1313, but have densities and pressures more similar to the interarm regions and so are less inclined to collapse and form stars. The cloud properties in NGC 1313 and NGC 7793 suggest that spiral arms, bars, interarm regions, and flocculent spirals each represent distinct environments with regard to molecular cloud populations. We see surprisingly little difference in surface densities between the regions, suggesting that the differences in surface densities frequently seen between arm and interarm regions of lower-resolution studies are indicative of the sparsity of molecular clouds, rather than differences in their true surface density.

preprint2022arXiv

A tale of two DIGs: The relative role of HII regions and low-mass hot evolved stars in powering the diffuse ionised gas (DIG) in PHANGS-MUSE galaxies

We use integral field spectroscopy from the PHANGS-MUSE survey, which resolves the ionised interstellar medium at ${\sim}50$ pc resolution in 19 nearby spiral galaxies, to study the origin of the diffuse ionised gas (DIG). We examine the physical conditions of the diffuse gas by first removing morphologically-defined HII regions and then binning the low-surface-brightness areas to achieve significant detections of the key nebular lines. A simple model for the leakage and propagation of ionising radiation from HII regions is able to reproduce the observed distribution of H$α$ in the DIG. Leaking radiation from HII regions also explains the observed decrease in line ratios of low-ionisation species ([SII]/H$α$, [NII]/H$α$ and [OI]/H$α$) with increasing H$α$ surface brightness ($Σ_{Hα}$). Emission from hot low-mass evolved stars, however, is required to explain: (1) the enhanced low-ionisation line ratios observed in the central regions of some galaxies; (2) the observed trends of a flat or decreasing [OIII]/H$β$ with $Σ_{Hα}$; and (3) the offset of some DIG regions from the locus of HII regions in the Baldwin-Phillips-Terlevich (BPT) diagram, extending into the area of low-ionisation (nuclear) emission-line regions (LI[N]ERs). Hot low-mass evolved stars make a small contribution to the energy budget of the DIG (2% of the galaxy-integrated H$α$ emission), but their harder spectra make them fundamental contributors to [OIII] emission. The DIG might result from a superposition of two components, an energetically dominant contribution from young stars and a more diffuse background of harder ionising photons from old stars. This unified framework bridges observations of the Milky Way DIG with LI(N)ER-like emission observed in nearby galaxy bulges.

preprint2022arXiv

After The Fall: Resolving the Molecular Gas in Post-Starburst Galaxies

Post-starburst (PSB), or 'E+A', galaxies represent a rapid transitional phase between major, gas-rich mergers and gas-poor, quiescent early-type galaxies. Surprisingly, many PSBs have been shown to host a significant interstellar medium (ISM), despite theoretical predictions that the majority of star-forming gas should be expelled in AGN- or starburst-driven outflows. To-date, the resolved properties of this surviving ISM have remained unknown. We present high resolution ALMA continuum and CO(2$-$1) observations in six gas- and dust-rich PSBs, revealing for the first time the spatial and kinematic structure of their ISM on sub-kpc scales. We find extremely compact molecular reservoirs, with dust and gas surface densities rivaling those found in (ultra-)luminous infrared galaxies. We observe spatial and kinematic disturbances in all sources, with some also displaying disk-like kinematics. Estimates of the internal turbulent pressure in the gas exceed those of normal star-forming disks by at least 2 orders of magnitude, and rival the turbulent gas found in local interacting galaxies, such as the Antennae. Though the source of this high turbulent pressure remains uncertain, we suggest that the high incidence of tidal disruption events (TDEs) in PSBs could play a role. The star formation in these PSBs' turbulent central molecular reservoirs is suppressed, forming stars only 10% as efficiently as starburst galaxies with similar gas surface densities. "The fall" of star formation in these galaxies was not precipitated by complete gas expulsion or redistribution. Rather, this high-resolution view of PSBs' ISM indicates that star formation in their remaining compact gas reservoirs is suppressed by significant turbulent heating.

preprint2022arXiv

Environmental dependence of the molecular cloud lifecycle in 54 main sequence galaxies

The processes of star formation and feedback, regulating the cycle of matter between gas and stars on the scales of giant molecular clouds (GMCs; $\sim$100pc), play a major role in governing galaxy evolution. Measuring the time-scales of GMC evolution is important to identify and characterise the specific physical mechanisms that drive this transition. By applying a robust statistical method to high-resolution CO and narrow-band H$α$ imaging from the PHANGS survey, we systematically measure the evolutionary timeline from molecular clouds to exposed young stellar regions on GMC scales, across the discs of an unprecedented sample of 54 star-forming main-sequence galaxies (excluding their unresolved centres). We find that clouds live for about $1{-}3$ GMC turbulence crossing times ($5{-}30$Myr) and are efficiently dispersed by stellar feedback within $1{-}5$Myr once the star-forming region becomes partially exposed, resulting in integrated star formation efficiencies of $1{-}8$%. These ranges reflect physical galaxy-to-galaxy variation. In order to evaluate whether galactic environment influences GMC evolution, we correlate our measurements with average properties of the GMCs and their local galactic environment. We find several strong correlations that can be physically understood, revealing a quantitative link between galactic-scale environmental properties and the small-scale GMC evolution. Notably, the measured CO-visible cloud lifetimes become shorter with decreasing galaxy mass, mostly due to the increasing presence of CO-dark molecular gas in such environment. Our results represent a first step towards a comprehensive picture of cloud assembly and dispersal, which requires further extension and refinement with tracers of the atomic gas, dust, and deeply-embedded stars.

preprint2022arXiv

Extragalactic magnetism with SOFIA (SALSA Legacy Program) -- III: First data release and on-the-fly polarization mapping characterization

We describe the data processing of the Survey on extragALactic magnetiSm with SOFIA (SALSA Legacy Program). This first data release presents 33% (51.34h out of 155.7h, including overheads) of the total awarded time taken from January 2020 to December 2021. Our observations were performed using the newly implemented on-the-fly mapping (OTFMAP) technique in the polarimetric mode. We present the pipeline steps to obtain homogeneously reduced high-level data products of polarimetric maps of galaxies for use in scientific analysis. Our approach has a general design and can be applied to sources smaller than the field-of-view of the HAWC+ array in any given band. We estimate that the OTFMAP polarimetric mode offers a reduction of observing overheads by a factor 2.34, and an improvement in sensitivity by a factor 1.80 when compared to previously obtained polarimetric observations using the chopping and nodding mode. The OTFMAP is a significant optimization of the polarimetric mode of HAWC+ as it ultimately reduces the cost of operations of SOFIA/HAWC+ by increasing the science collected per hour of observation up to an overall factor of 2.49. The OTFMAP polarimetric mode is the standard observing strategy of SALSA. The results and quantitative analysis of this first data release are presented in Papers IV and V of the series.

preprint2022arXiv

Extragalactic magnetism with SOFIA (SALSA Legacy Program) -- IV: Program overview and first results on the polarization fraction

We present the first data release of the Survey on extragALactic magnetiSm with SOFIA (SALSA Legacy Program) with a set of 14 nearby ($<20$ Mpc) galaxies with resolved imaging polarimetric observations using HAWC+ from $53$ to $214$ $μ$m at a resolution of $5-18$" ($90$ pc $-$ $1$ kpc). We introduce the definitions and background on extragalactic magnetism, and present the scientific motivation and sample selection of the program. Here, we focus on the general trends in the emissive polarization fraction. Far-infrared polarimetric observations trace the thermal polarized emission of magnetically aligned dust grains across the galaxy disks with polarization fractions of $P=0-15$% in the cold, $T_{\rm d} = [19,48]$ K, and dense, $\log_{10}(N_{\rm HI+H_{2}}) = [19.96,22.91]$, interstellar medium. The spiral galaxies show a median $\langle P_{154μm} \rangle = 3.3\pm0.9 $% across the disks. We report the first polarized spectrum of starburst galaxies showing a minimum within $89-154$ $μ$m. The falling $53-154$ $μ$m polarized spectrum may be due to a decrease in the dust grain alignment efficiency produced by variations in dust temperatures along the line-of-sight in the galactic outflow. We find that the starburst galaxies and the star-forming regions within normal galaxies have the lowest polarization fractions. We find that 50% (7 out of 14) of the galaxies require a broken power-law in the $P-N_{HI+H_{2}}$ and $P-T_{d}$ relations with three different trends. Group 1 has a relative increase of anisotropic random B-fields produced by compression or shear of B-fields in the galactic outflows, starburst rings, and inner-bar of galaxies; and Groups 2 and 3 have a relative increase of isotropic random B-fields driven by star-forming regions in the spiral arms, and/or an increase of dust grain alignment efficiency caused by shock-driven regions or evolutionary stages of a galaxy.

preprint2022arXiv

Extragalactic magnetism with SOFIA (SALSA Legacy Program). VI. The magnetic fields in the multi-phase interstellar medium of the Antennae galaxies

Mergers are thought to be a fundamental channel for galaxy growth, perturbing the gas dynamics and the magnetic fields (B-fields) in the interstellar medium (ISM). However, the mechanisms that amplify and dissipate B-fields during a merger remain unclear. We characterize the morphology of the ordered B-fields in the multi-phase ISM of the closest merger of two spiral galaxies, the Antennae galaxies. We compare the inferred B-fields using $154~μ$m thermal dust and $11$ cm radio synchrotron emission polarimetric observations. We find that the $154~μ$m B-fields are more ordered across the Antennae galaxies than the $11$ cm B-fields. The turbulent-to-ordered $154~μ$m B-field increases at the galaxy cores and star-forming regions. The relic spiral arm has an ordered spiral $154~μ$m B-field, while the $11$ cm B-field is radial. The $154~μ$m B-field may be dominated by turbulent dynamos with high $^{12}$CO(1-0) velocity dispersion driven by star-forming regions, while the $11$ cm B-field is cospatial with high HI velocity dispersion driven by galaxy interaction. This result shows the dissociation between the warm gas mainly disturbed by the merger, and the dense gas still following the dynamics of the relic spiral arm. We find a $\sim8.9$ kpc scale ordered B-field connecting the two galaxies. The base of the tidal tail is cospatial with the HI and $^{12}$CO(1-0) emission and has compressed and/or sheared $154~μ$m and $11$ cm B-fields driven by the merger. We suggest that amplify B-fields, with respect to the rest of the system and other spiral galaxies, may be supporting the gas flow between both galaxies and the tidal tail.

preprint2022arXiv

The Gas-Star Formation Cycle in Nearby Star-forming Galaxies II. Resolved Distributions of CO and H$α$ Emission for 49 PHANGS Galaxies

The relative distribution of molecular gas and star formation in galaxies gives insight into the physical processes and timescales of the cycle between gas and stars. In this work, we track the relative spatial configuration of CO and H$α$ emission at high resolution in each of our galaxy targets, and use these measurements to quantify the distributions of regions in different evolutionary stages of star formation: from molecular gas without star formation traced by H$α$ to star-forming gas, and to HII regions. The large sample, drawn from the Physics at High Angular resolution in Nearby GalaxieS ALMA and narrowband H$α$ (PHANGS-ALMA and PHANGS-H$α$) surveys, spans a wide range of stellar mass and morphological types, allowing us to investigate the dependencies of the gas-star formation cycle on global galaxy properties. At a resolution of 150 pc, the incidence of regions in different stages shows a dependence on stellar mass and Hubble type of galaxies over the radial range probed. Massive and/or earlier-type galaxies exhibit a significant reservoir of molecular gas without star formation traced by H$α$, while lower-mass galaxies harbor substantial HII regions that may have dispersed their birth clouds or formed from low-mass, more isolated clouds. Galactic structures add a further layer of complexity to relative distribution of CO and H$α$ emission. Trends between galaxy properties and distributions of gas traced by CO and H$α$ are visible only when the observed spatial scale is $\ll$ 500 pc, reflecting the critical resolution requirement to distinguish stages of star formation process.

preprint2022arXiv

The PHANGS-MUSE survey -- Probing the chemo-dynamical evolution of disc galaxies

We present the PHANGS-MUSE survey, a programme using the MUSE IFS at the ESO VLT to map 19 massive $(9.4 < \log(M_{*}/M_\odot) < 11.0)$ nearby (D < 20 Mpc) star-forming disc galaxies. The survey consists of 168 MUSE pointings (1'x1' each), a total of nearly 15 Million spectra, covering ~1.5 Million independent spectra. PHANGS-MUSE provides the first IFS view of star formation across different local environments (including galaxy centres, bars, spiral arms) in external galaxies at a median resolution of 50~pc, better than the mean inter-cloud distance in the ionised interstellar medium. This `cloud-scale' resolution allows detailed demographics and characterisations of HII regions and other ionised nebulae. PHANGS-MUSE further delivers a unique view on the associated gas and stellar kinematics, and provides constraints on the star formation history. The PHANGS-MUSE survey is complemented by dedicated ALMA CO(2-1) and multi-band HST observations, therefore allowing us to probe the key stages of the star formation process from molecular clouds to HII regions and star clusters. This paper describes the scientific motivation, sample selection, observational strategy, data reduction and analysis process of the PHANGS-MUSE survey. We present our bespoke automated data-reduction framework, which is built on the reduction recipes provided by ESO, but additionally allows for mosaicking and homogenisation of the point spread function. We further present a detailed quality assessment and a brief illustration of the potential scientific applications of the large set of PHANGS-MUSE data products generated by our data analysis framework. The data cubes and analysis data products described in this paper represent the basis for the first PHANGS-MUSE public data release and are available in the ESO archive and via the Canadian Astronomy Data Centre.

preprint2021arXiv

Giant Molecular Cloud Catalogues for PHANGS-ALMA: Methods and Initial Results

We present improved methods for segmenting CO emission from galaxies into individual molecular clouds, providing an update to the CPROPS algorithms presented by Rosolowsky & Leroy (2006; arXiv:astro-ph/0601706 ). The new code enables both homogenization of the noise and spatial resolution among data, which allows for rigorous comparative analysis. The code also models the completeness of the data via false source injection and includes an updated segmentation approach to better deal with blended emission. These improved algorithms are implemented in a publicly available python package, PYCPROPS. We apply these methods to ten of the nearest galaxies in the PHANGS-ALMA survey, cataloguing CO emission at a common 90 pc resolution and a matched noise level. We measure the properties of 4986 individual clouds identified in these targets. We investigate the scaling relations among cloud properties and the cloud mass distributions in each galaxy. The physical properties of clouds vary among galaxies, both as a function of galactocentric radius and as a function of dynamical environment. Overall, the clouds in our target galaxies are well-described by approximate energy equipartition, although clouds in stellar bars and galaxy centres show elevated line widths and virial parameters. The mass distribution of clouds in spiral arms has a typical mass scale that is 2.5x larger than interarm clouds and spiral arms clouds show slightly lower median virial parameters compared to interarm clouds (1.2 versus 1.4).

preprint2021arXiv

Low-J CO Line Ratios From Single Dish CO Mapping Surveys and PHANGS-ALMA

We measure the low-J CO line ratio R21=CO(2-1)/CO(1-0), R32=CO(3-2)/CO(2-1), and R31 = CO(3-2)/CO(1-0) using whole-disk CO maps of nearby galaxies. We draw CO(2-1) from PHANGS--ALMA, HERACLES, and follow-up IRAM surveys; CO(1-0) from COMING and the Nobeyama CO Atlas of Nearby Spiral Galaxies; and CO(3-2) from the JCMT NGLS and APEX LASMA mapping. Altogether this yields 76, 47, and 29 maps of R21, R32, and R31 at 20" \sim 1.3 kpc resolution, covering 43, 34, and 20 galaxies. Disk galaxies with high stellar mass, log10 M_* [Msun]=10.25-11 and star formation rate, SFR=1-5 Msun/yr, dominate the sample. We find galaxy-integrated mean values and 16%-84% range of R21 = 0.65 (0.50-0.83), R32=0.50 (0.23-0.59), and R31=0.31 (0.20-0.42). We identify weak trends relating galaxy-integrated line ratios to properties expected to correlate with excitation, including SFR/M_* and SFR/L_CO. Within galaxies, we measure central enhancements with respect to the galaxy-averaged value of \sim 0.18^{+0.09}_{-0.14} dex for R21, 0.27^{+0.13}_{-0.15} dex for R31, and 0.08^{+0.11}_{-0.09} dex for R32. All three line ratios anti-correlate with galactocentric radius and positively correlate with the local star formation rate surface density and specific star formation rate, and we provide approximate fits to these relations. The observed ratios can be reasonably reproduced by models with low temperature, moderate opacity, and moderate densities, in good agreement with expectations for the cold ISM. Because the line ratios are expected to anti-correlate with the CO(1-0)-to-H_2 conversion factor, alphaCO^(1-0), these results have general implications for the interpretation of CO emission from galaxies.

preprint2021arXiv

PHANGS-HST: Star Cluster Spectral Energy Distribution Fitting with CIGALE

The sensitivity and angular resolution of photometric surveys executed by the Hubble Space Telescope (HST) enable studies of individual star clusters in galaxies out to a few tens of megaparsecs. The fitting of spectral energy distributions (SEDs) of star clusters is essential for measuring their physical properties and studying their evolution. We report on the use of the publicly available Code Investigating GALaxy Emission (CIGALE) SED fitting package to derive ages, stellar masses, and reddenings for star clusters identified in the Physics at High Angular resolution in Nearby GalaxieS-HST (PHANGS-HST) survey. Using samples of star clusters in the galaxy NGC 3351, we present results of benchmark analyses performed to validate the code and a comparison to SED fitting results from the Legacy ExtraGalactic Ultraviolet Survey (LEGUS). We consider procedures for the PHANGS-HST SED fitting pipeline, e.g., the choice of single stellar population models, the treatment of nebular emission and dust, and the use of fluxes versus magnitudes for the SED fitting. We report on the properties of clusters in NGC 3351 and find, on average, the clusters residing in the inner star-forming ring of NGC 3351 are young ($< 10$ Myr) and massive ($10^{5} M_{\odot}$) while clusters in the stellar bulge are significantly older. Cluster mass function fits yield $β$ values around -2, consistent with prior results with a tendency to be shallower at the youngest ages. Finally, we explore a Bayesian analysis with additional physically-motivated priors for the distribution of ages and masses and analyze the resulting cluster distributions.

preprint2021arXiv

PHANGS-MUSE: the HII region luminosity function of local star-forming galaxies

We use an unprecedented sample of about 23,000 HII regions detected at an average physical resolution of 67pc in the PHANGS-MUSE sample to study the extragalactic HII region Ha luminosity function (LF). Our observations probe the star-forming disk of 19 nearby spiral galaxies with low inclination and located close to the star formation main sequence at z=0. The mean LF slope $α$ in our sample is =1.73 with a $σ$ of 0.15. We find that $α$ decreases with the galaxy's star formation rate surface density and argue that this is driven by an enhanced clustering of young stars at high gas surface densities. Looking at the HII regions within single galaxies we find that no significant variations occur between the LF of the inner and outer part of the star-forming disk, whereas the LF in the spiral arm areas is shallower than in the inter-arm areas for six out of the 13 galaxies with clearly visible spiral arms. We attribute these variations to the spiral arms increasing the molecular clouds arm--inter-arm mass contrast and find suggestive evidence that they are more evident for galaxies with stronger spiral arms. Furthermore, we find systematic variations in $α$ between samples of HII regions with high and low ionization parameter q and argue that they are driven by the aging of HII regions.

preprint2021arXiv

The Case for Thermalization as a Contributor to the [CII] Deficit

The [CII] deficit, which describes the observed decrease in the ratio of [CII] 158 micron emission to continuum infrared emission in galaxies with high star formation surface densities, places a significant challenge to the interpretation of [CII] detections from across the observable universe. In an attempt to further decode the cause of the [CII] deficit, the [CII] and dust continuum emission from 18 Local Volume galaxies has been split based on conditions within the interstellar medium where it originated. This is completed using the Key Insights in Nearby Galaxies: a Far-Infrared Survey with Herschel (KINGFISH) and Beyond the Peak (BtP) surveys and the wide-range of wavelength information, from UV to far-infrared emission lines, available for a selection of star-forming regions within these samples. By comparing these subdivided [CII] emissions to isolated infrared emission and other properties, we find that the thermalization (collisional de-excitation) of the [CII] line in HII regions plays a significant role in the deficit observed in our sample.

preprint2021arXiv

The PHANGS-HST Survey: Physics at High Angular resolution in Nearby GalaxieS with the Hubble Space Telescope

The PHANGS program is building the first dataset to enable the multi-phase, multi-scale study of star formation across the nearby spiral galaxy population. This effort is enabled by large survey programs with ALMA, VLT/MUSE, and HST, with which we have obtained CO(2-1) imaging, optical spectroscopic mapping, and high resolution UV-optical imaging, respectively. Here, we present PHANGS-HST, which has obtained five band NUV-U-B-V-I imaging of the disks of 38 spiral galaxies at distances of 4-23 Mpc, and parallel V and I band imaging of their halos, to provide a census of tens of thousands of compact star clusters, and multi-scale stellar associations. The combination of HST, ALMA, and VLT/MUSE observations will yield an unprecedented joint catalog of the observed and physical properties of ~100,000 star clusters, associations, HII regions, and molecular clouds. With these basic units of star formation, PHANGS will systematically chart the evolutionary cycling between gas and stars, across a diversity of galactic environments found in nearby galaxies. We discuss the design of the PHANGS-HST survey, and provide an overview of the HST data processing pipeline and first results, highlighting new methods for selecting star cluster candidates, morphological classification of candidates with convolutional neural networks, and identification of stellar associations over a range of physical scales with a watershed algorithm. We describe the cross-observatory imaging, catalogs, and software products to be released. These high-level science products will seed a broad range of investigations, in particular, the study of embedded stellar populations and dust with JWST, for which a PHANGS Cycle 1 Treasury program to obtain eight band 2-21 $μ$m imaging has been approved.

preprint2020arXiv

[CI](1-0) and [CI](2-1) in resolved local galaxies

We present resolved [CI] line intensities of 18 nearby galaxies observed with the SPIRE FTS spectrometer on the Herschel Space Observatory. We use these data along with resolved CO line intensities from $J_\mathrm{up} = 1$ to 7 to interpret what phase of the interstellar medium the [CI] lines trace within typical local galaxies. A tight, linear relation is found between the intensities of the CO(4-3) and [CI](2-1) lines; we hypothesize this is due to the similar upper level temperature of these two lines. We modeled the [CI] and CO line emission using large velocity gradient models combined with an empirical template. According to this modeling, the [CI](1-0) line is clearly dominated by the low-excitation component. We determine [CI] to molecular mass conversion factors for both the [CI](1-0) and [CI](2-1) lines, with mean values of $α_{\mathrm{[CI](1-0)}} = 7.3$ M$_{\mathrm{sun}}$ K$^{-1}$ km$^{-1}$ s pc$^{-2}$ and $α_{\mathrm{[CI](2-1)}} = 34 $ M$_{\mathrm{sun}}$ K$^{-1}$ km$^{-1}$ s pc$^{-2}$ with logarithmic root-mean-square spreads of 0.20 and 0.32 dex, respectively. The similar spread of $α_{\mathrm{[CI](1-0)}}$ to $α_{\mathrm{CO}}$ (derived using the CO(2-1) line) suggests that [CI](1-0) may be just as good a tracer of cold molecular gas as CO(2-1) in galaxies of this type. On the other hand, the wider spread of $α_{\mathrm{[CI](2-1)}}$ and the tight relation found between [CI](2-1) and CO(4-3) suggest that much of the [CI](2-1) emission may originate in warmer molecular gas.

preprint2020arXiv

Candidate LBV stars in galaxy NGC 7793 found via HST photometry + MUSE spectroscopy

Only about 19 Galactic and 25 extra-galactic bona-fide Luminous Blue Variables (LBVs) are known to date. This incomplete census prevents our understanding of this crucial phase of massive star evolution which leads to the formation of heavy binary black holes via the classical channel. With large samples of LBVs one could better determine the duration and maximum stellar luminosity which characterize this phase. We search for candidate LBVs (cLBVs) in a new galaxy, NGC 7793. For this purpose, we combine high spatial resolution images from two Hubble Space Telescope (HST) programs with optical spectroscopy from the Multi Unit Spectroscopic Explorer (MUSE). By combining PSF-fitting photometry measured on F547M, F657N, and F814W images, with restrictions on point-like appearance (at HST resolution) and Hα luminosity, we find 100 potential cLBVs, 36 of which fall in the MUSE fields. Five of the latter 36 sources are promising cLBVs which have M(V) less than or equal to -7 and a combination of: Hα with a P-Cygni profile; no [O I] 6300 emission; weak or no [O III] 5007 emission; large [N II]/Hα relative to H II regions; and [S II] 6716 / [S II] 6731 similar to 1. It is not clear if these five cLBVs are isolated from O-type stars, which would favor the binary formation scenario of LBVs. Our study, which approximately covers one fourth of the optical disc of NGC 7793, demonstrates how by combining the above HST surveys with multi-object spectroscopy from 8-m class telescopes, one can efficiently find large samples of cLBVs in nearby galaxies.

preprint2020arXiv

Deep Transfer Learning for Star Cluster Classification: I. Application to the PHANGS-HST Survey

We present the results of a proof-of-concept experiment which demonstrates that deep learning can successfully be used for production-scale classification of compact star clusters detected in HST UV-optical imaging of nearby spiral galaxies (D < 20 Mpc) in the PHANGS-HST survey. Given the relatively small nature of existing, human-labelled star cluster samples, we transfer the knowledge of state-of-the-art neural network models for real-object recognition to classify star clusters candidates into four morphological classes. We perform a series of experiments to determine the dependence of classification performance on: neural network architecture (ResNet18 and VGG19-BN); training data sets curated by either a single expert or three astronomers; and the size of the images used for training. We find that the overall classification accuracies are not significantly affected by these choices. The networks are used to classify star cluster candidates in the PHANGS-HST galaxy NGC 1559, which was not included in the training samples. The resulting prediction accuracies are 70%, 40%, 40-50%, 50-70% for class 1, 2, 3 star clusters, and class 4 non-clusters respectively. This performance is competitive with consistency achieved in previously published human and automated quantitative classification of star cluster candidate samples (70-80%, 40-50%, 40-50%, and 60-70%). The methods introduced herein lay the foundations to automate classification for star clusters at scale, and exhibit the need to prepare a standardized dataset of human-labelled star cluster classifications, agreed upon by a full range of experts in the field, to further improve the performance of the networks introduced in this study.

preprint2020arXiv

Radial Star Formation Histories in 32 Nearby Galaxies

The spatially resolved star formation histories are studied for 32 normal star-forming galaxies drawn from the the Spitzer Extended Disk Galaxy Exploration Science survey. At surface brightness sensitivities fainter than 28 mag arcsec$^{-2}$, the new optical photometry is deep enough to complement archival ultraviolet and infrared imaging and to explore the properties of the emission well beyond the traditional optical extents of these nearby galaxies. Fits to the spectral energy distributions using a delayed star formation history model indicate a subtle but interesting average radial trend for the spiral galaxies: the inner stellar systems decrease in age with increasing radius, consistent with inside-out disk formation, but the trend reverses in the outermost regions with the stellar age nearly as old as the innermost stars. These results suggest an old stellar outer disk population formed through radial migration and/or the cumulative history of minor mergers and accretions of satellite dwarf galaxies. The subset of S0 galaxies studied here show the opposite trend compared to what is inferred for spirals: characteristic stellar ages that are increasingly older with radius for the inner portions of the galaxies, and increasingly younger stellar ages for the outer portions. This result suggests that either S0 galaxies are not well modeled by a delayed-$τ$ model, and/or that S0 galaxies have a more complicated formation history than spiral galaxies.

preprint2019arXiv

Mapping electron temperature variations across a spiral arm in NGC 1672

We report one of the first extragalactic observations of electron temperature variations across a spiral arm. Using MUSE mosaic observations of the nearby galaxy NGC 1672, we measure the [N II]5755 auroral line in a sample of 80 HII regions in the eastern spiral arm of NGC1672. We discover systematic temperature variations as a function of distance perpendicular to the spiral arm. The electron temperature is lowest on the spiral arm itself and highest on the downstream side. Photoionization models of different metallicity, pressure, and age of the ionizing source are explored to understand what properties of the interstellar medium drive the observed temperature variations. An azimuthally varying metallicity appears to be the most likely cause of the temperature variations. The electron temperature measurements solidify recent discoveries of azimuthal variations of oxygen abundance based on strong lines, and rule out the possibility that the abundance variations are artefacts of the strong-line calibrations.

preprint2019arXiv

The headlight cloud in NGC 628: An extreme giant molecular cloud in a typical galaxy disk

Cloud-scale surveys of molecular gas reveal the link between molecular clouds properties and star formation (SF) across a range of galactic environments. Cloud populations in galaxy disks are considered to be representative of the `normal' SF. At high resolution, however, clouds with exceptional gas properties and SF activity may also be observed in normal disk environments. In this paper, we study the brightest cloud traced in CO emission in the disk of NGC628. The cloud is spatially coincident with an extremely bright HII region. We characterize its molecular gas properties and investigate how feedback and large-scale processes influence the properties of the molecular gas. High resolution CO ALMA observations are used to characterize its mass and dynamical state, which are compared to other clouds in NGC628. A LVG analysis is used to constrain the beam-diluted density and temperature of the molecular gas. We analyze the MUSE spectrum using Starburst99 to characterize the young stellar population associated with the HII region. The cloud is massive ($1-2\times10^7$M$_{\odot}$), with a beam-diluted density of $n_{\rm H_2}=5\times10^4$ cm$^{-3}$. It has a low virial parameter, suggesting that its CO emission may be overluminous due to heating by the HII region. A young ($2-4$ Myr), massive $3\times10^{5}$ M$_{\odot}$ stellar population is associated. We argue that the cloud is currently being destroyed by feedback from young massive stars. Due to the cloud's large mass, this phase of the cloud's evolution is long enough for the impact of feedback on the excitation of the gas to be observed. Its high mass may be related to its location at a spiral co-rotation radius, where gas experiences reduced galactic shear compared to other regions of the disk, and receives a sustained inflow of gas that can promote the cloud's mass growth.

preprint2019arXiv

Using [CII] 158 micron Emission From Isolated ISM Phases as a Star-Formation Rate Indicator

The brightest observed emission line in many star-forming galaxies is the [CII] 158 micron line, making it detectable up to z~7. In order to better understand and quantify the [CII] emission as a tracer of star-formation, the theoretical ratio between the [NII] 205 micron emission and the [CII] 158 micron emission has been employed to empirically determine the fraction of [CII] emission that originates from the ionized and neutral phases of the ISM. Sub-kiloparsec measurements of the [CII] 158 micron and [NII] 205 micron line in nearby galaxies have recently become available as part of the Key Insights in Nearby Galaxies: a Far Infrared Survey with Herschel (KINGFISH) and Beyond the Peak (BtP) programs. With the information from these two far-infrared lines along with the multi-wavelength suite of KINGFISH data, a calibration of the [CII] emission line as a star formation rate indicator and a better understanding of the [CII] deficit are pursued. [CII] emission is also compared to PAH emission in these regions to compare photoelectric heating from PAH molecules to cooling by [CII] in the neutral and ionized phases of the ISM. We find that the [CII] emission originating in the neutral phase of the ISM does not exhibit a deficit with respect to the infrared luminosity and is therefore preferred over the [CII] emission originating in the ionized phase of the ISM as a star formation rate indicator for the normal star-forming galaxies included in this sample.

preprint2016arXiv

A comprehensive search for stellar bowshock nebulae in the Milky Way: a catalog of 709 mid-infrared selected candidates

We identify 709 arc-shaped mid-infrared nebula in 24 micron Spitzer Space Telescope or 22 micron Wide Field Infrared Explorer surveys of the Galactic Plane as probable dusty interstellar bowshocks powered by early-type stars. About 20% are visible at 8 microns or shorter mid-infrared wavelengths as well. The vast majority (660) have no previous identification in the literature. These extended infrared sources are strongly concentrated near Galactic mid-Plane with an angular scale height of ~0.6 degrees. All host a symmetrically placed star implicated as the source of a stellar wind sweeping up interstellar material. These are candidate "runaway" stars potentially having high velocities in the reference frame of the local medium. Among the 286 objects with measured proper motions, we find an unambiguous excess having velocity vectors aligned with the infrared morphology --- kinematic evidence that many of these are "runaway" stars with large peculiar motions responsible for the bowshock signature. We discuss a population of "in-situ" bowshocks (103 objects) that face giant HII regions where the relative motions between the star and ISM may be caused by bulk outflows from an overpressured bubble. We also identify 58 objects that face 8 micron bright-rimmed clouds and apparently constitute a sub-class of in-situ bowshocks where the stellar wind interacts with a photo-evaporative flow from an eroding molecular cloud interface (i.e., "PEF bowshocks"). Orientations of the arcuate nebulae exhibit a correlation over small angular scales, indicating that external influences such as HII regions are responsible for producing some bowshock nebulae. However, the vast majority of this sample appear to be isolated (499 objects) from obvious external influences.

preprint2016arXiv

The Connection Between Galaxy Environment and the Luminosity Function Slopes of Star-Forming Regions

We present the first study of GALEX far ultra-violet (FUV) luminosity functions of individual star-forming regions within a sample of 258 nearby galaxies spanning a large range in total stellar mass and star formation properties. We identify ~65,000 star-forming regions (i.e., FUV sources), measure each galaxy's luminosity function, and characterize the relationships between the luminosity function slope (alpha) and several global galaxy properties. A final sample of 82 galaxies with reliable luminosity functions are used to define these relationships and represent the largest sample of galaxies with the largest range of galaxy properties used to study the connection between luminosity function properties and galaxy environment. We find that alpha correlates with global star formation properties, where galaxies with higher star formation rates and star formation rate densities (Sigma_SFR) tend to have flatter luminosity function slopes. In addition, we find that neither stochastic sampling of the luminosity function in galaxies with low-number statistics nor the effects of blending due to distance can fully account for these trends. We hypothesize that the flatter slopes in high Sigma_SFR galaxies is due to higher gas densities and higher star formation efficiencies which result in proportionally greater numbers of bright star-forming regions. Finally, we create a composite luminosity function composed of star-forming regions from many galaxies and find a break in the luminosity function at brighter luminosities. However, we find that this break is an artifact of varying detection limits for galaxies at different distances.

preprint2016arXiv

The Spatially Resolved [CII] Cooling Line Deficit in Galaxies

We present [CII] 158um measurements from over 15,000 resolved regions within 54 nearby galaxies of the KINGFISH program to investigate the so-called [CII] "line cooling deficit" long known to occur in galaxies with different luminosities. The [CII]/TIR ratio ranges from above 1% to below 0.1% in the sample, with a mean value of 0.48+-0.21%. We find that the surface density of 24um emission dominates this trend, with [CII]/TIR dropping as nuInu{24um} increases. Deviations from this overall decline are correlated with changes in the gas phase metal abundance, with higher metallicity associated with deeper deficits at a fixed surface brightness. We supplement the local sample with resolved [CII] measurements from nearby luminous infrared galaxies and high redshift sources from z=1.8-6.4, and find that star formation rate density drives a continuous trend of deepening [CII] deficit across six orders of magnitude in SFRD. The tightness of this correlation suggests that an approximate star formation rate density can be estimated directly from global measurements of [CII]/TIR, and a relation is provided to do so. Several low-luminosity AGN hosts in the sample show additional and significant central suppression of [CII]/TIR, but these deficit enhancements occur not in those AGN with the highest X-ray luminosities, but instead those with the highest central starlight intensities. Taken together, these results demonstrate that the [CII] cooling line deficit in galaxies likely arises from local physical phenomena in interstellar gas.

preprint2015arXiv

The Stellar Halo and Tidal Streams of Messier 63

We present new near-infrared (NIR) observations of M63 from the Extended Disk Galaxy Exploration Science (EDGES) Survey. The extremely deep 3.6 $μ$m mosaic reaches 29 AB mag arcsec$^{-2}$ at the outer reaches of the azimuthally-averaged surface brightness profile. At this depth the consequences of galactic accretion are found within a nearby tidal stream and an up-bending break in the slope of the surface brightness profile. This break occurs at a semi-major axis length of $\sim$8', and is evidence of either an enhanced outer disc or an inner stellar halo. Simulations of galaxy evolution, along with our observations, support an inner halo as the explanation for the up-bending break. The mass of this halo component is the largest found in an individual galaxy thus far. Additionally, our observations detect a nearby tidal stream. The mass of the stream suggests that a handful of such accretion events are necessary to populate the inner stellar halo. We also find that the accretion rate of the galaxy from the stream alone underestimates the accretion rate required to build M63's inner stellar halo.

preprint2014arXiv

A Two-Parameter Model for the Infrared/Submillimeter/Radio Spectral Energy Distributions of Galaxies and AGN

A two-parameter semi-empirical model is presented for the spectral energy distributions of galaxies with contributions to their infrared-submillimeter-radio emission from both star formation and accretion disk-powered activity. This model builds upon a previous one-parameter family of models for star-forming galaxies, and includes an update to the mid-infrared emission using an average template obtained from Spitzer Space Telescope observations of normal galaxies. Star-forming/AGN diagnostics based on PAH equivalent widths and broadband infrared colors are presented, and example mid-infrared AGN fractional contributions are estimated from model fits to the GOALS sample of nearby U/LIRGS and the 5MUSES sample of 24um-selected sources at redshifts 0 < z < 2.

preprint2014arXiv

Empirical ugri-UBVRc Transformations for Galaxies

We present empirical color transformations between Sloan Digital Sky Survey ugri and Johnson-Cousins UBVRc photometry for nearby galaxies (D < 11 Mpc). We use the Local Volume Legacy (LVL) galaxy sample where there are 90 galaxies with overlapping observational coverage for these two filter sets. The LVL galaxy sample consists of normal, non-starbursting galaxies. We also examine how well the LVL galaxy colors are described by previous transformations derived from standard calibration stars and model-based galaxy templates. We find significant galaxy color scatter around most of the previous transformation relationships. In addition, the previous transformations show systematic offsets between transformed and observed galaxy colors which are visible in observed color-color trends. The LVL-based $galaxy$ transformations show no systematic color offsets and reproduce the observed color-color galaxy trends.

preprint2014arXiv

Spitzer Local Volume Legacy (LVL) SEDs and Physical Properties

We present the panchromatic spectral energy distributions (SEDs) of the Local Volume Legacy (LVL) survey which consists of 258 nearby galaxies ($D<$11 Mpc). The wavelength coverage spans the ultraviolet to the infrared (1500 $\textrmÅ$ to 24 $μ$m) which is utilized to derive global physical properties (i.e., star formation rate, stellar mass, internal extinction due to dust.). With these data, we find color-color relationships and correlated trends between observed and physical properties (i.e., optical magnitudes and dust properties, optical color and specific star formation rate, and ultraviolet-infrared color and metallicity). The SEDs are binned by different galaxy properties to reveal how each property affects the observed shape of these SEDs. In addition, due to the volume-limited nature of LVL, we utilize the dwarf-dominated galaxy sample to test star formation relationships established with higher-mass galaxy samples. We find good agreement with the star-forming "main-sequence" relationship, but find a systematic deviation in the infrared "main-sequence" at low luminosities. This deviation is attributed to suppressed polycyclic aromatic hydrocarbon (PAH) formation in low metallicity environments and/or the destruction of PAHs in more intense radiation fields occurring near a suggested threshold in sSFR at a value of log($sSFR$) $\sim$ $-$10.2.

preprint2014arXiv

The Relationship between Stellar Mass, Gas Metallicity, and Star Formation Rate for Halpha-selected Galaxies at z~0.8 from the NewHalpha Survey

Using a sample of 299 Ha-selected galaxies at z~0.8, we study the relationship between galaxy stellar mass, gas-phase metallicity, and star formation rate (SFR), and compare to previous results. We use deep optical spectra obtained with the IMACS spectrograph at the Magellan telescope to measure strong oxygen lines. We combine these spectra and metallicities with (1) rest-frame UV-to-optical imaging, which allows us to determine stellar masses and dust attenuation corrections, and (2) Ha narrowband imaging, which provides a robust measure of the instantaneous SFR. Our sample spans stellar masses of 10^9 to 6*10^11 solar masses, SFRs of 0.4 to 270 solar masses per year, and metal abundances of 12+log(O/H)~8.3-9.1 (~0.4-2.6 solar metallicity). The correlations that we find between the Ha-based SFR and stellar mass (i.e., the star-forming "main sequence"), and between the stellar mass and metallicity, are both consistent with previous z~1 studies of star-forming galaxies. We then study the relationship between the three properties using various plane-fitting techniques (Lara-Lopez et al.) and a curve-fitting projection (Mannucci et al.). In all cases, we exclude strong dependence of the M-Z relation on SFR, but are unable to distinguish between moderate and no dependence. Our results are consistent with previous mass-metallicity-SFR studies. We check whether dataset limitations may obscure a strong dependence on the SFR by using mock samples drawn from the SDSS. These experiments reveal that the adopted signal-to-noise cuts may have a significant effect on the measured dependence. Further work is needed to investigate these results, and to test whether a "fundamental metallicity relation" or a "fundamental plane" describes star-forming galaxies across cosmic time.

preprint2014arXiv

The Spitzer Local Volume Legacy (LVL) Global Optical Photometry

We present the global optical photometry of 246 galaxies in the Local Volume Legacy (LVL) survey. The full volume-limited sample consists of 258 nearby (D < 11 Mpc) galaxies whose absolute B-band magnitude span a range of -9.6 < M_B < -20.7 mag. A composite optical (UBVR) data set is constructed from observed UBVR and SDSS ugriz imaging, where the ugriz magnitudes are transformed into UBVR. We present photometry within three galaxy apertures defined at UV, optical, and IR wavelengths. Flux comparisons between these apertures reveal that the traditional optical R25 galaxy apertures do not fully encompass extended sources. Using the larger IR apertures we find color-color relationships where later-type spiral and irregular galaxies tend to be bluer than earlier-type galaxies. These data provide the missing optical emission from which future LVL studies can construct the full panchromatic (UV-optical-IR) spectral energy distributions.

preprint2014arXiv

Toward Complete Statistics of Massive Binary Stars: Penultimate Results from the Cygnus OB2 Radial Velocity Survey

We analyze orbital solutions for 48 massive multiple-star systems in the Cygnus OB2 Association, 23 of which are newly presented here, to find that the observed distribution of orbital periods is approximately uniform in log P for P<45 d, but it is not scale-free. Inflections in the cumulative distribution near 6 d, 14, d, and 45 d, suggest key physical scales of ~0.2, ~0.4, and ~1 A.U. where yet-to-be-identified phenomena create distinct features. No single power law provides a statistically compelling prescription, but if features are ignored, a power law with exponent beta = -0.22 provides a crude approximation over P=1.4 -- 2000 d, as does a piece-wise linear function with a break near 45 d. The cumulative period distribution flattens at P > 45 d, even after correction for completeness, indicating either a lower binary fraction or a shift toward low-mass companions. A high degree of similarity (91% likelihood) between the Cyg OB2 period distribution and that of other surveys suggests that the binary properties at P<25 d are determined by local physics of disk/clump fragmentation and are relatively insensitive to environmental and evolutionary factors. Fully 30% of the unbiased parent sample is a binary with period P < 45 d. Completeness corrections imply a binary fraction near 55% for P < 5000 d. The observed distribution of mass ratios 0.2<q<1 is consistent with uniform, while the observed distribution of eccentricities 0.1<e<0.6 is consistent with uniform plus an excess of e ~ 0 systems. We identify six stars, all supergiants, that exhibit aperiodic velocity variations of ~30 km/s attributed to atmospheric fluctuations.

preprint2013arXiv

Measuring Galaxy Star Formation Rates From Integrated Photometry: Insights from Color-Magnitude Diagrams of Resolved Stars

We use empirical star formation histories (SFHs), measured from HST-based resolved star color-magnitude diagrams, as input into population synthesis codes to model the broadband spectral energy distributions (SEDs) of ~50 nearby dwarf galaxies (6.5 < log M/M_* < 8.5, with metallicities ~10% solar). In the presence of realistic SFHs, we compare the modeled and observed SEDs from the ultraviolet (UV) through near-infrared (NIR) and assess the reliability of widely used UV-based star formation rate (SFR) indicators. In the FUV through i bands, we find that the observed and modeled SEDs are in excellent agreement. In the Spitzer 3.6micron and 4.5micron bands, we find that modeled SEDs systematically over-predict observed luminosities by up to ~0.2 dex, depending on treatment of the TP-AGB stars in the synthesis models. We assess the reliability of UV luminosity as a SFR indicator, in light of independently constrained SFHs. We find that fluctuations in the SFHs alone can cause factor of ~2 variations in the UV luminosities relative to the assumption of a constant SFH over the past 100 Myr. These variations are not strongly correlated with UV-optical colors, implying that correcting UV-based SFRs for the effects of realistic SFHs is difficult using only the broadband SED. Additionally, for this diverse sample of galaxies, we find that stars older than 100 Myr can contribute from <5% to100% of the present day UV luminosity, highlighting the challenges in defining a characteristic star formation timescale associated with UV emission. We do find a relationship between UV emission timescale and broadband UV-optical color, though it is different than predictions based on exponentially declining SFH models. Our findings have significant implications for the comparison of UV-based SFRs across low-metallicity populations with diverse SFHs.

preprint2013arXiv

Mid-infrared Spectral Properties of Post-Starburst Quasars

We present Spitzer InfraRed Spectrograph (IRS) low-resolution spectra of 16 spectroscopically selected post-starburst quasars (PSQs) at z ~ 0.3. The optical spectra of these broad-lined active galactic nuclei (AGNs) simultaneously show spectral signatures of massive intermediate-aged stellar populations making them good candidates for studying the connections between AGNs and their hosts. The resulting spectra show relatively strong polycyclic aromatic hydrocarbon (PAH) emission features at 6.2 and 11.3\micron and a very weak silicate feature, indicative of ongoing star formation and low dust obscuration levels for the AGNs. We find that the mid-infrared composite spectrum of PSQs has spectral properties between ULIRGs and QSOs, suggesting that PSQs are hybrid AGN and starburst systems as also seen in their optical spectra. We also find that PSQs in early-type host galaxies tend to have relatively strong AGN activities, while those in spiral hosts have stronger PAH emission, indicating more star formation.

preprint2013arXiv

The Nature of the Second Parameter in the IRX-beta Relation for Local Galaxies

We present an analysis of 98 galaxies of low-dust content, selected from the Spitzer Local Volume Legacy survey, aimed at examining the relation between the ultraviolet (UV) color, beta, and dust attenuation in normal star-forming galaxies. The IRX-beta diagram relates the total dust attenuation in a galaxy, traced by the far-IR (FIR) to UV ratio, to the observed UV color, indicated by beta. Previous research has indicated that while starburst galaxies exhibit a relatively tight IRX-beta relation, normal star-forming galaxies do not, and have a much larger spread in the total-IR to far-UV (FUV) luminosity for a fixed UV color. We examine the role that the age of the stellar population plays as the "second parameter" responsible for the observed deviation and spread of star-forming galaxies from the starburst relation. We model the FUV to FIR spectral energy distribution (SED) of each galaxy according to two broad bins of star formation history (SFH): constant and instantaneous burst. We find clear trends between stellar population mean age estimators (extinction-corrected FUV/NIR, U-B, and EW(H-alpha)) and the UV color, beta; the trends are mostly driven by the galaxies best-described by instantaneous burst populations. We also find a significant correlation between beta and the mean age directly determined from the best-fit instantaneous models. As already indicated by other authors, the UV attenuation in star-forming galaxies may not be recovered with the UV color alone and is highly influenced by the stellar population's mean age and SFH. Overall, the scatter in the IRX-beta diagram is better correlated with beta than with the perpendicular distance.

preprint2012arXiv

Direct Oxygen Abundances for Low Luminosity LVL Galaxies

We present MMT spectroscopic observations of HII regions in 42 low luminosity galaxies in the LVL. For 31 galaxies, we measured the temperature sensitive [O III] line at a strength of 4 sigma or greater, and thus determine direct oxygen abundances. Our results provide the first direct estimates of oxygen abundance for 19 galaxies. Oxygen abundances were compared to B-band and 4.5 micron luminosities and stellar masses in order to characterize the luminosity-metallicity (L-Z) and mass-metallicity (M-Z) relationships at low-luminosity. We present and analyze a "Combined Select" sample composed of 38 objects (drawn from our parent sample and the literature) with direct oxygen abundances and reliable distance determinations (TRGB or Ceph). Consistent with previous studies, the B-band and 4.5 micron L-Z relationships were found to be 12+log(O/H)=(6.27+/-0.21)+(-0.11+/-0.01)M_B and 12+log(O/H)=(6.10+/-0.21)+(-0.10+/-0.01)M_[4.5] (sigma=0.15 and 0.14). For this sample, we derive a M-Z relationship of 12+log(O/H)=(5.61+/-0.24)+(0.29+/-0.03)log(M*), which agrees with previous studies; however, the dispersion (sigma=0.15) is not significantly lower than that of the L-Z relationships. Because of the low dispersions in these relationships, if an accurate distance is available, the luminosity of a low-luminosity galaxy is often a better indicator of metallicity than that derived using certain strong-line methods, so significant departures from the L-Z relationships may indicate that caution is prudent in such cases. We also revisit the 70/160 micron color metallicity relationship. Additionally, we examine N/O abundance trends with respect to oxygen abundance and B-V color. We find a positive correlation between N/O ratio and B-V color for 0.05\lesssimB-V\lesssim0.75: log(N/O)=(1.18+/-0.9)x(B-V)+(-1.92+/-0.08), with a dispersion of sigma=0.14, that is in agreement with previous studies.

preprint2012arXiv

Nebular Attenuation in Hα-selected Star-forming Galaxies at z=0.8 from the NewHα Survey

We present measurements of the dust attenuation of Hα-selected emission-line galaxies at z=0.8 from the NewHα narrowband survey. The analysis is based on deep follow-up spectroscopy with Magellan/IMACS, which captures the strong rest-frame optical emission lines from [OII] λ3727 to [OIII] λ5007. The spectroscopic sample used in this analysis consists of 341 confirmed Hα emitters. We place constraints on the AGN fraction using diagnostics which can be applied at intermediate redshift. We find that at least 5% of the objects in our spectroscopic sample can be classified as AGN and 2% are composite, i.e. powered by a combination of star-formation and AGN activity. We measure the dust attenuation for individual objects from the ratios of the higher order Balmer lines. The Hβ and Hγ pair of lines is detected with S/N>5 in 55 individual objects and the Hβ and Hδ pair is detected in 50 individual objects. We also create stacked spectra to probe the attenuation in objects without individual detections. The median attenuation at Hα based on the objects with individually detected lines is A(Hα)=0.9+-1.0 magnitudes, in good agreement with the attenuation found in local samples of star-forming galaxies. We find that the z=0.8 galaxies occupy a similar locus of attenuation as a function of magnitude, mass and SFR as a comparison sample drawn from the SDSS DR4. Both the results from the individual z=0.8 galaxies and from the stacked spectra show consistency with the mass -- attenuation and SFR -- attenuation relations found in the local Universe, indicating that these relations are also applicable at intermediate redshift.

preprint2012arXiv

The ACS Nearby Galaxy Survey Treasury. X. Quantifying the Star Cluster Formation Efficiency of Nearby Dwarf Galaxies

We study the relationship between the field star formation and cluster formation properties in a large sample of nearby dwarf galaxies. We use optical data from the Hubble Space Telescope and from ground-based telescopes to derive the ages and masses of the young (t_age < 100Myr) cluster sample. Our data provides the first constraints on two proposed relationships between the star formation rate of galaxies and the properties of their cluster systems in the low star formation rate regime. The data show broad agreement with these relationships, but significant galaxy-to-galaxy scatter exists. In part, this scatter can be accounted for by simulating the small number of clusters detected from stochastically sampling the cluster mass function. However, this stochasticity does not fully account for the observed scatter in our data suggesting there may be true variations in the fraction of stars formed in clusters in dwarf galaxies. Comparison of the cluster formation and the brightest cluster in our sample galaxies also provide constraints on cluster destruction models.

preprint2011arXiv

Average Metallicity and Star Formation Rate of Lya Emitters Probed by a Triple Narrow-Band Survey

We present the average metallicity and star-formation rate of Lya emitters (LAEs) measured from our large-area survey with three narrow-band (NB) filters covering the Lya, [OII]3727, and Ha+[NII] lines of LAEs at z=2.2. We select 919 z=2.2 LAEs from Subaru/Suprime-Cam NB data in conjunction with Magellan/IMACS spectroscopy. Of these LAEs, 561 and 105 are observed with KPNO/NEWFIRM near-infrared NB filters whose central wavelengths are matched to redshifted [OII] and Ha nebular lines, respectively. By stacking the near-infrared images of the LAEs, we successfully obtain average nebular-line fluxes of LAEs, the majority of which are too faint to be identified individually by narrow-band imaging or deep spectroscopy. The stacked object has an Ha luminosity of 1.7x10^{42} erg s^{-1} corresponding to a star formation rate (SFR) of 14 M_{sun} yr^{-1}. We place, for the first time, a firm lower limit to the average metallicity of LAEs of Z>~0.09 Z_{sun} (2sigma) based on the [OII]/(Ha+[NII]) index together with photo-ionization models and empirical relations. This lower limit of metallicity rules out the hypothesis that LAEs, so far observed at z~2, are extremely metal poor (Z<2x10^{-2} Z_{sun}) young galaxies at the 4sigma level. This limit is higher than a simple extrapolation of the observed mass-metallicity relation of z~2 UV-selected galaxies toward lower masses (5x10^{8} M_{sun}), but roughly consistent with a recently proposed fundamental mass-metallicity relation when the LAEs' relatively low SFR is taken into account. The Ha and Lya luminosities of our NB-selected LAEs indicate that the escape fraction of Lya photons is ~12-30 %, much higher than the values derived for other galaxy populations at z~2.

preprint2011arXiv

Dust-Corrected Star Formation Rates of Galaxies. II. Combinations of Ultraviolet and Infrared Tracers

We present new calibrations of far-ultraviolet (FUV) attenuation as derived from the total infrared to FUV luminosity ratio (IRX) and the FUV-NUV color. We find that the IRX-corrected FUV luminosities are tightly and linearly correlated with the attenuation-corrected Hα luminosities (as measured from the Balmer decrement), with a rms scatter of $\pm 0.09$ dex. The ratios of these attenuation-corrected FUV to Hα luminosities are consistent with evolutionary synthesis model predictions, assuming a constant star formation rate over 100 Myr, solar metallicity and either a Salpeter or a Kroupa IMF with lower and upper mass limits of 0.1 and 100\msun. The IRX-corrected FUV to Balmer-corrected Hα luminosity ratios do not show any trend with other galactic properties over the ranges covered by our sample objects. In contrast, FUV attenuation derived from the FUV-NUV color (UV spectral slope) show much larger random and systematic uncertainties. When compared to either Balmer-corrected Hα luminosities or IRX-corrected FUV luminosities the color-corrected FUV luminosities show $\sim 2.5$ times larger rms scatter, and systematic nonlinear deviations as functions of luminosity and other parameters. Linear combinations of 25um and 1.4GHz radio continuum luminosities with the observed FUV luminosities are also well correlated with the Balmer-corrected Hα luminosities. These results provide useful prescriptions for deriving attenuation-corrected star formation rates of galaxies based on linear combinations of UV and IR or radio luminosities, which are presented in convenient tabular form. Comparisons of our calibrations with attenuation corrections in the literature and with dust attenuation laws are also made.

preprint2011arXiv

Modeling the Effects of Star Formation Histories on Halpha and Ultra-Violet Fluxes in Nearby Dwarf Galaxies

We consider the effects of non-constant star formation histories (SFHs) on Halpha and GALEX far ultra-violet (FUV) star formation rate (SFR) indicators. Under the assumption of a fully populated Chabrier IMF, we compare the distribution of Halpha-to-FUV flux ratios from ~ 1500 simple, periodic model SFHs with observations of 185 galaxies from the Spitzer Local Volume Legacy survey. We find a set of SFH models that are well matched to the data, such that more massive galaxies are best characterized by nearly constant SFHs, while low mass systems experience bursts amplitudes of ~ 30 (i.e., an increase in the SFR by a factor of 30 over the SFR during the inter-burst period), burst durations of tens of Myr, and periods of ~ 250 Myr; these SFHs are broadly consistent with the increased stochastic star formation expected in systems with lower SFRs. We analyze the predicted temporal evolution of galaxy stellar mass, R-band surface brightness, Halpha-derived SFR, and blue luminosity, and find that they provide a reasonable match to observed flux distributions. We find that our model SFHs are generally able to reproduce both the observed systematic decline and increased scatter in Halpha-to-FUV ratios toward low mass systems, without invoking other physical mechanisms. We also compare our predictions with those from the Integrated Galactic IMF theory with a constant SFR. We find that while both predict a systematic decline in the observed ratios, only the time variable SFH models are capable of producing the observed population of low mass galaxies ($M_{*}$ < 10$^{7}$ Msun) with normal Halpha-to-FUV ratios. These results demonstrate that a variable IMF alone has difficulty explaining the observed scatter in the Halpha-to-FUV ratios. We conclude by considering the limitations of the model SFHs, and discuss the use of additional empirical constraints to improve future SFH modeling efforts.

preprint2011arXiv

The Evolution of Stellar Populations in the Outer Disks of Spiral Galaxies

We investigate recent star formation in the extended ultraviolet (XUV) disks of five nearby galaxies (NGC 0628, NGC 2090, NGC 2841, NGC 3621, and NGC 5055) using a long wavelength baseline comprised of ultraviolet and mid-infrared imaging from the Galaxy Evolution Explorer and the Spitzer Infrared Array Camera. We identify 229 unresolved stellar complexes across targeted portions of their XUV disks and utilize spectral energy distribution fitting to measure their stellar ages and masses through comparison with Starburst99 population synthesis models of instantaneous burst populations. We find that the median age of outer disk associations in our sample is ~100 Myr with a large dispersion that spans the entire range of our models (1 Myr-1 Gyr). This relatively evolved state for most associations addresses the observed dearth of Halpha emission in some outer disks, as Halpha can only be observed in star forming regions younger than ~10 Myr. The large age dispersion is robust against variations in extinction (in the range E(B-V)=0-0.3 mag) and variations in the upper end of the stellar Initial Mass Function (IMF). In particular, we demonstrate that the age dispersion is insensitive to steepening of the IMF, up to extreme slopes.

preprint2011arXiv

The H-alpha Luminosity Function and Star Formation Rate Volume Density at z=0.8 from the NEWFIRM H-alpha Survey

[Abridged] We present new measurements of the H-alpha luminosity function (LF) and SFR volume density for galaxies at z~0.8. Our analysis is based on 1.18$μ$m narrowband data from the NEWFIRM H-alpha Survey, a comprehensive program designed to capture deep samples of intermediate redshift emission-line galaxies using narrowband imaging in the near-infrared. The combination of depth ($\approx1.9\times10^{-17}$ erg s$^{-1}$ cm$^{-2}$ in H-alpha at 3$σ$) and areal coverage (0.82 deg$^2$) complements other recent H-alpha studies at similar redshifts, and enables us to minimize the impact of cosmic variance and place robust constraints on the shape of the LF. The present sample contains 818 NB118 excess objects, 394 of which are selected as H-alpha emitters. Optical spectroscopy has been obtained for 62% of the NB118 excess objects. Empirical optical broadband color classification is used to sort the remainder of the sample. A comparison of the LFs constructed for the four individual fields reveals significant cosmic variance, emphasizing that multiple, widely separated observations are required. The dust-corrected LF is well-described by a Schechter function with L*=10^{43.00\pm0.52} ergs s^{-1}, ϕ*=10^{-3.20\pm0.54} Mpc^{-3}, and α=-1.6\pm0.19. We compare our H-alpha LF and SFR density to those at z<1, and find a rise in the SFR density \propto(1+z)^{3.4}, which we attribute to significant L* evolution. Our H-alpha SFR density of 10^{-1.00\pm0.18} M_sun yr^{-1} Mpc^{-3} is consistent with UV and [O II] measurements at z~1. We discuss how these results compare to other H-alpha surveys at z~0.8, and find that the different methods used to determine survey completeness can lead to inconsistent results. This suggests that future surveys probing fainter luminosities are needed, and more rigorous methods of estimating the completeness should be adopted as standard procedure.

preprint2011arXiv

The mid-IR Luminosity Function at z<0.3 from 5MUSES: Understanding the Star-formation/AGN Balance from a Spectroscopic View

We present rest-frame 15 and 24 um luminosity functions and the corresponding star-forming luminosity functions at z<0.3 derived from the 5MUSES sample. Spectroscopic redshifts have been obtained for ~98% of the objects and the median redshift is ~0.12. The 5-35 um IRS spectra allow us to estimate accurately the luminosities and build the luminosity functions. Using a combination of starburst and quasar templates, we quantify the star-formation and AGN contributions in the mid-IR SED. We then compute the star-formation luminosity functions at 15 um and 24 um, and compare with the total 15 um and 24 um luminosity functions. When we remove the contribution of AGN, the bright end of the luminosity function exhibits a strong decline, consistent with the exponential cutoff of a Schechter function. Integrating the differential luminosity function, we find that the fractional contribution by star formation to the energy density is 58% at 15 um and 78% at 24 um, while it goes up to ~86% when we extrapolate our mid-IR results to the total IR luminosity density. We confirm that the active galactic nuclei play more important roles energetically at high luminosities. Finally, we compare our results with work at z~0.7 and confirm that evolution on both luminosity and density is required to explain the difference in the LFs at different redshifts.

preprint2010arXiv

Optical Spectroscopy and Nebular Oxygen Abundances of the Spitzer/SINGS Galaxies

We present intermediate-resolution optical spectrophotometry of 65 galaxies obtained in support of the Spitzer Infrared Nearby Galaxies Survey (SINGS). For each galaxy we obtain a nuclear, circumnuclear, and semi-integrated optical spectrum designed to coincide spatially with mid- and far-infrared spectroscopy from the Spitzer Space Telescope. We make the reduced, spectrophotometrically calibrated one-dimensional spectra, as well as measurements of the fluxes and equivalent widths of the strong nebular emission lines, publically available. We use optical emission-line ratios measured on all three spatial scales to classify the sample into star-forming, active galactic nuclei (AGN), and galaxies with a mixture of star formation and nuclear activity. We find that the relative fraction of the sample classified as star-forming versus AGN is a strong function of the integrated light enclosed by the spectroscopic aperture. We supplement our observations with a large database of nebular emission-line measurements of individual HII regions in the SINGS galaxies culled from the literature. We use these ancillary data to conduct a detailed analysis of the radial abundance gradients and average HII-region abundances of a large fraction of the sample. We combine these results with our new integrated spectra to estimate the central and characteristic (globally-averaged) gas-phase oxygen abundances of all 75 SINGS galaxies. We conclude with an in-depth discussion of the absolute uncertainty in the nebular oxygen abundance scale.

preprint2010arXiv

The Wyoming Survey for H-alpha. III. A Multi-wavelength Look at Attenuation by Dust in Galaxies out to z~0.4

We report results from the Wyoming Survey for H-alpha (WySH), a comprehensive four-square degree survey to probe the evolution of star-forming galaxies over the latter half of the age of the Universe. We have supplemented the H-alpha data from WySH with infrared data from the Spitzer Wide-area Infrared Extragalactic (SWIRE) Survey and ultraviolet data from the Galaxy Evolution Explorer (GALEX) Deep Imaging Survey. This dataset provides a multi-wavelength look at the evolution of the attenuation by dust, and here we compare a traditional measure of dust attenuation (L(TIR)/L(FUV)) to a diagnostic based on a recently-developed robust star formation rate (SFR) indicator, [H-alpha_obs+24-micron]/H-alpha_obs. With such data over multiple epochs, the evolution in the attenuation by dust with redshift can be assessed. We present results from the ELAIS-N1 and Lockman Hole regions at z~0.16, 0.24, 0.32 and 0.40. While the ensemble averages of both diagnostics are relatively constant from epoch to epoch, each epoch individually exhibits a larger attenuation by dust for higher star formation rates. Hence, an epoch to epoch comparison at a fixed star formation rate suggests a mild decrease in dust attenuation with redshift.

preprint2010arXiv

The Wyoming Survey for H-alpha. III. H-alpha Luminosity Functions at z ~ 0.16, 0.24, 0.32, and 0.40

The Wyoming Survey for H-alpha, or WySH, is a large-area, ground-based imaging survey for H-alpha-emitting galaxies at redshifts of z ~ 0.16, 0.24, 0.32, and 0.40. The survey spans up to four square degrees in a set of fields of low Galactic cirrus emission, using twin narrowband filters at each epoch for improved stellar continuum subtraction. H-alpha luminosity functions are presented for each Delta(z) ~ 0.02 epoch based on a total of nearly 1200 galaxies. These data clearly show an evolution with lookback time in the volume-averaged cosmic star formation rate. Integrals of Schechter fits to the incompleteness- and extinction-corrected H-alpha luminosity functions indicate star formation rates per co-moving volume of 0.010, 0.013, 0.020, 0.022 h_70 M_sun yr^{-1} Mpc^{-3} at z ~ 0.16, 0.24, 0.32, and 0.40, respectively. Statistical and systematic measurement uncertainties combined are on the order of 25% while the effects of cosmic variance are at the 20% level. The bulk of this evolution is driven by changes in the characteristic luminosity L_* of the H-alpha luminosity functions, with L_* for the earlier two epochs being a factor of two larger than L_* at the latter two epochs; it is more difficult with this data set to decipher systematic evolutionary differences in the luminosity function amplitude and faint-end slope. Coupling these results with a comprehensive compilation of results from the literature on emission line surveys, the evolution in the cosmic star formation rate density over 0 < z < 1.5 is measured to be rho_dot_SFR(z) = rho_dot_SFR(0) (1+z)^{3.4+/-0.4}.

preprint2009arXiv

Dust-Corrected Star Formation Rates of Galaxies. I. Combinations of H-alpha and Infrared Tracers

We combine H-alpha emission-line and infrared continuum measurements of two samples of nearby galaxies to derive dust attenuation-corrected star formation rates (SFRs). We use a simple energy balance based method that has been applied previously to HII regions in the Spitzer Infrared Nearby Galaxies Survey (SINGS), and extend the methodology to integrated measurements of galaxies. We find that our composite H-alpha+IR based SFRs are in excellent agreement with attenuation-corrected SFRs derived from integrated spectrophotometry, over the full range of SFRs (0.01 -- 80 solar mass per year) and attenuations (0 -- 2.5 mag) studied. We find that the combination of H-alpha and total infrared luminosities provides the most robust SFR measurements, but combinations of H-alpha measurements with monochromatic luminosities at 24 micron and 8 micron perform nearly as well. The calibrations differ significantly from those obtained for HII regions (Calzetti et al. 2007), with the difference attributable to a more evolved population of stars heating the dust. Our results are consistent with a significant component of diffuse dust (the `IR cirrus' component) that is heated by a non-star-forming population. The same methodology can be applied to [OII]$λ$3727 emission-line measurements, and the radio continuum fluxes of galaxies can be applied in place of IR fluxes when the latter are not available. We assess the precision and systematic reliability of all of these composite methods.

preprint2007arXiv

Star Formation in NGC 5194 (M51a). II. The Spatially-Resolved Star Formation Law

We have studied the relationship between the star formation rate (SFR) surface density and gas surface density in the spiral galaxy M51a (NGC 5194), using multi-wavelength data obtained as part of the Spitzer Infrared Nearby Galaxies Survey (SINGS). We introduce a new SFR index based on a linear combination of H-alpha emission-line and 24 micron continuum luminosities, that provides reliable extinction-corrected ionizing fluxes and SFR densities over a wide range of dust attenuations. The combination of these extinction-corrected SFR densities with aperture synthesis HI and CO maps has allowed us to probe the form of the spatially-resolved star formation law on scales of 0.5 to 2 kpc. We find that the resolved SFR vs gas surface density relation is well represented by a Schmidt power law, which is similar in form and dispersion to the disk-averaged Schmidt law. We observe a comparably strong correlation of the SFR surface density with the molecular gas surface density, but no significant correlation with the surface density of atomic gas. The best-fitting slope of the Schmidt law varies from N = 1.37 to 1.56, with zeropoint and slope that change systematically with the spatial sampling scale. We tentatively attribute these variations to the effects of areal sampling and averaging of a nonlinear intrinsic star formation law. Our data can also be fitted by an alternative parametrization of the SFR surface density in terms of the ratio of gas surface density to local dynamical time, but with a considerable dispersion.

preprint2006arXiv

The Spectral Energy Distribution of Dust Emission in the Edge-on Spiral Galaxy NGC 4631 as seen with Spitzer and JCMT

We explore the nature of variations in dust emission within an individual galaxy using 3.6 - 160 micron Spitzer Space Telescope observations and 450 and 850 micron James Clerk Maxwell Telescope observations of the edge-on Sd spiral galaxy NGC 4631 with the goals of understanding the relation between polycyclic aromatic hydrocarbons (PAH) and dust emission, studying the variations in the colors of the dust emission, and searching for possible excess submillimeter emission compared to what is anticipated based on dust models applied to the mid- and far-infrared data. PAH emission at 8 micron is found to correlate best with hot dust emission at 24 micron on kiloparsec scales, although the relation breaks down on scales equal to hundreds of parsecs, possibly because of differences in the mean free paths between the photons that excite the PAHs and heat the dust and possibly because the PAHs are destroyed by the hard radiation fields in the centers of some star formation regions. The ratio of 8 micron PAH emission to 160 micron cool dust emission appears to be a function of radius. The 70/160 and 160/450 micron flux density ratios are remarkably constant even though the surface brightness varies by a factor of 25 in each wave band, which suggests that the emission is from dust heated by a nearly-uniform radiation field. Globally, we find an excess of 850 micron emission relative to what would be predicted by dust models. The 850 micron excess is highest in regions with low 160 micron surface brightness, although the strength and statistical significance of this result depends on the model fit to the data. We rule out variable emissivity functions or ~4 K dust as the possible origins of this 850 micron emission, but we do discuss the other possible mechanisms that could produce the emission.

preprint2006arXiv

Warm Dust and Spatially Variable PAH Emission in the Dwarf Starburst Galaxy NGC 1705

We present Spitzer observations of the dwarf starburst galaxy NGC 1705 obtained as part of SINGS. The galaxy morphology is very different shortward and longward of ~5 microns: short-wavelength imaging shows an underlying red stellar population, with the central super star cluster (SSC) dominating the luminosity; longer-wavelength data reveals warm dust emission arising from two off-nuclear regions offset by ~250 pc from the SSC. These regions show little extinction at optical wavelengths. The galaxy has a relatively low global dust mass (~2E5 solar masses, implying a global dust-to-gas mass ratio ~2--4 times lower than the Milky Way average). The off-nuclear dust emission appears to be powered by photons from the same stellar population responsible for the excitation of the observed H Alpha emission; these photons are unassociated with the SSC (though a contribution from embedded sources to the IR luminosity of the off-nuclear regions cannot be ruled out). Low-resolution IRS spectroscopy shows moderate-strength PAH emission in the 11.3 micron band in the eastern peak; no PAH emission is detected in the SSC or the western dust emission complex. There is significant diffuse 8 micron emission after scaling and subtracting shorter wavelength data; the spatially variable PAH emission strengths revealed by the IRS data suggest caution in the interpretation of diffuse 8 micron emission as arising from PAH carriers alone. The metallicity of NGC 1705 falls at the transition level of 35% solar found by Engelbracht and collaborators; the fact that a system at this metallicity shows spatially variable PAH emission demonstrates the complexity of interpreting diffuse 8 micron emission. A radio continuum non-detection, NGC 1705 deviates significantly from the canonical far-IR vs. radio correlation. (Abridged)

preprint2005arXiv

Spitzer Observations Of The Supergiant Shell Region In IC 2574

We present spatially resolved Spitzer imaging of the supergiant shell region of the M81 group dwarf galaxy IC 2574 obtained as part of the Spitzer Infrared Nearby Galaxies Survey. This region harbors one of the best nearby examples of a kinematically distinct HI shell, with an associated remnant stellar cluster; the shell is initiating sequential star formation as it interacts with the surrounding interstellar medium. This region dominates the infrared luminosity of IC 2574 and is spatially resolved in all Spitzer imaging bands. We study the differences in dust temperature as a function of local environment and compare local star formation rates as inferred from H Alpha and total infrared luminosities. We find that the strong H Alpha sources are associated with regions of warm dust; however, the most luminous infrared and H Alpha sources are not necessarilyco-spatial. The coolest dust is found in the regions farthest from the rim of the shell; these regions show the best agreement between star formation rates derived from H Alpha and from total infrared luminosities (although discrepancies at the factor of 3-4 level still exist). There is considerable variation in the radio-far infrared correlation in different regions surrounding the shell. The low dust content of the region may influence the scatter seen in these relations; these data demonstrate that the expanding shell is dramatically affecting its surroundings by triggering star formation and altering the dust temperature.

preprint2003arXiv

Near-Infrared Integral Field Spectroscopy of Star-Forming Galaxies

The Palomar Integral Field Spectrograph was used to probe a variety of environments in nine nearby galaxies that span a range of morphological types, luminosities, metallicities, and infrared-to-blue ratios. For the first time, near-infrared spectroscopy was obtained for nuclear or bright HII regions in star-forming galaxies over two spatial dimensions (5.7"x10.0") in the [FeII] (1.257um), [FeII] (1.644um), Paschen beta (1.282um), H2 (2.122um), and Brackett gamma (2.166um) transition lines. These data yield constraints on various characteristics of the star-forming episodes in these regions, including their strength, maturity, spatial variability, and extinction. The HII regions stand out from the nuclei. Unlike observations of nuclear regions, HII region near-infrared observations do not show a spatial coincidence of the line and continuum emission; the continuum and line maps of HII regions usually show distinct and sometimes spatially-separated morphologies. Gauging from Paschen beta and Brackett gamma equivalent widths and luminosities, the HII regions have younger episodes of star formation than the nuclei and more intense radiation fields. Near-infrared line ratio diagnostics suggest that HII regions have "purer" starbursting properties. The correlation between ionizing photon density and mid-infrared color is consistent with the star formation activity level being higher for HII regions than for nuclei. And though the interpretation is complicated, on a purely empirical basis the HII regions show lower Fe1+ abundances than nuclei by an order of magnitude.