Using the star cluster catalogs from the Hubble Space Telescope program Legacy ExtraGalactic UV survey (LEGUS) and 8 $\mu$m images from the IRAC camera on the Spitzer Space Telescope for 5 galaxies within 5 Mpc, we investigate how the 8 $\mu$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 $\mu$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 $\mu$m luminosity and the age of stellar clusters younger than 1 Gyr; the 8 $\mu$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 $\mu$m emission as a function of age are provided.

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-$\tau$ model, and/or that S0 galaxies have a more complicated formation history than spiral galaxies.

We compare the observed turbulent pressure in molecular gas, $P_\mathrm{turb}$, to the required pressure for the interstellar gas to stay in equilibrium in the gravitational potential of a galaxy, $P_\mathrm{DE}$. To do this, we combine arcsecond resolution CO data from PHANGS-ALMA with multi-wavelength data that traces the atomic gas, stellar structure, and star formation rate (SFR) for 28 nearby star-forming galaxies. We find that $P_\mathrm{turb}$ correlates with, but almost always exceeds the estimated $P_\mathrm{DE}$ on kiloparsec scales. This indicates that the molecular gas is over-pressurized relative to the large-scale environment. We show that this over-pressurization can be explained by the clumpy nature of molecular gas; a revised estimate of $P_\mathrm{DE}$ on cloud scales, which accounts for molecular gas self-gravity, external gravity, and ambient pressure, agrees well with the observed $P_\mathrm{turb}$ in galaxy disks. We also find that molecular gas with cloud-scale ${P_\mathrm{turb}}\approx{P_\mathrm{DE}}\gtrsim{10^5\,k_\mathrm{B}\,\mathrm{K\,cm^{-3}}}$ in our sample is more likely to be self-gravitating, whereas gas at lower pressure appears more influenced by ambient pressure and/or external gravity. Furthermore, we show that the ratio between $P_\mathrm{turb}$ and the observed SFR surface density, $\Sigma_\mathrm{SFR}$, is compatible with stellar feedback-driven momentum injection in most cases, while a subset of the regions may show evidence of turbulence driven by additional sources. The correlation between $\Sigma_\mathrm{SFR}$ and kpc-scale $P_\mathrm{DE}$ in galaxy disks is consistent with the expectation from self-regulated star formation models. Finally, we confirm the empirical correlation between molecular-to-atomic gas ratio and kpc-scale $P_\mathrm{DE}$ reported in previous works.

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{\alpha} 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{\alpha} with a P-Cygni profile; no [O I] 6300 emission; weak or no [O III] 5007 emission; large [N II]/H{\alpha} 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.

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 $\alpha_{\mathrm{[CI](1-0)}} = 7.3$ M$_{\mathrm{sun}}$ K$^{-1}$ km$^{-1}$ s pc$^{-2}$ and $\alpha_{\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 $\alpha_{\mathrm{[CI](1-0)}}$ to $\alpha_{\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 $\alpha_{\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.

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.

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.

We report the first detection of galactic spiral structure by means of thermal emission from magnetically aligned dust grains. Our 89 $\mu$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 $\mu$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.

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${\alpha}$-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.

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.

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)

The processes regulating star formation in galaxies are thought to act across a hierarchy of spatial scales. To connect extragalactic star formation relations from global and kpc-scale measurements to recent cloud-scale resolution studies, we have developed a simple, robust method that quantifies the scale dependence of the relative spatial distributions of molecular gas and recent star formation. In this paper, we apply this method to eight galaxies with roughly 1 arcsec resolution molecular gas imaging from the PHANGS-ALMA and PAWS surveys that have matched resolution, high quality narrowband Halpha imaging. At a common scale of 140pc, our massive (log(Mstar/Msun)=9.3-10.7), normally star-forming (SFR/Msun/yr=0.3-5.9) galaxies exhibit a significant reservoir of quiescent molecular gas not associated with star formation as traced by Halpha emission. Galactic structures act as backbones for both molecular and HII region distributions. As we degrade the spatial resolution, the quiescent molecular gas disappears, with the most rapid changes occurring for resolutions up to about 0.5kpc. As the resolution becomes poorer, the morphological features become indistinct for spatial scales larger than about 1kpc. The method is a promising tool to search for relationships between the quiescent or star-forming molecular reservoir and galaxy properties, but requires a larger sample size to identify robust correlations between the star-forming molecular gas fraction and global galaxy parameters.

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.

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.

The morphology of HII regions around young star clusters provides insight into the timescales and physical processes that clear a cluster's natal gas. We study ~700 young clusters (<10Myr) in three nearby spiral galaxies (NGC 7793, NGC 4395, and NGC 1313) using Hubble Space Telescope (HST) imaging from LEGUS (Legacy ExtraGalactic Ultraviolet Survey). Clusters are classified by their H$\alpha$ morphology (concentrated, partially exposed, no-emission) and whether they have neighboring clusters (which could affect the clearing timescales). Through visual inspection of the HST images, and analysis of ages, reddenings, and stellar masses from spectral energy distributions fitting, together with the (U-B), (V-I) colors, we find: 1) the median ages indicate a progression from concentrated (~3 Myr), to partially exposed (~4 Myr), to no H$\alpha$ emission (>5Myr), consistent with the expected temporal evolution of HII regions and previous results. However, 2) similarities in the age distributions for clusters with concentrated and partially exposed H$\alpha$ morphologies imply a short timescale for gas clearing (<1Myr). 3) our cluster sample's median mass is ~1000 M, and a significant fraction (~20%) contain one or more bright red sources (presumably supergiants), which can mimic reddening effects. Finally, 4) the median E(B-V) values for clusters with concentrated H$\alpha$ and those without H$\alpha$ emission appear to be more similar than expected (~0.18 vs. ~0.14, respectively), but when accounting for stochastic effects, clusters without H$\alpha$ emission are less reddened. To mitigate stochastic effects, we experiment with synthesizing more massive clusters by stacking fluxes of clusters within each H$\alpha$ morphological class. Composite isolated clusters also reveal a color and age progression for H$\alpha$ morphological classes, consistent with analysis of the individual clusters.

We present results from a joint ALMA-HST study of the nearby spiral galaxy NGC 628. We combine the HST LEGUS database of over 1000 stellar clusters in NGC 628 with ALMA Cycle 4 millimeter/submillimeter observations of the cold dust continuum that span ~15 square kpc including the nuclear region and western portions of the galaxy's disk. The resolution -- 1.1" or approximately 50 pc at the distance of NGC 628 -- allows us to constrain the spatial variations in the slope of the millimeter dust continuum as a function of the ages and masses of the nearby stellar clusters. Our results indicate an excess of dust emission in the millimeter assuming a typical cold dust model for a normal star-forming galaxy, but little correlation of the dust continuum slope with stellar cluster age or mass. For the depth and spatial coverage of these observations, we cannot substantiate the millimeter/submillimeter excess arising from the processing of dust grains by the local interstellar radiation field. We detect a bright unknown source in NGC 628 in ALMA bands 4 and 7 with no counterparts at other wavelengths from ancillary data. We speculate this is possibly a dust obscured supernova.

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 in the PHANGS-HST survey. Given the relatively small and unbalanced nature of existing, human-labelled star cluster datasets, we transfer the knowledge of neural network models for real-object recognition to classify star clusters candidates into four morphological classes. We show that human classification is at the 66%:37%:40%:61% agreement level for the four classes considered. Our findings indicate that deep learning algorithms achieve 76%:63%:59%:70% for a star cluster sample within 4Mpc < D <10Mpc. We tested the robustness of our deep learning algorithms to generalize to different cluster images using the first data obtained by PHANGS-HST of NGC1559, which is more distant at D = 19Mpc, and found that deep learning produces classification accuracies 73%:42%:52%:67%. We furnish evidence for the robustness of these analyses by using two different neural network models for image classification, trained multiple times from the ground up to assess the variance and stability of our results. We quantified the importance of the NUV, U, B, V and I images for morphological classification with our deep learning models, and find that the V-band is the key contributor as human classifications are based on images taken in that filter. This work lays the foundations to automate classification for these objects at scale, and the creation of a standardized dataset.

We report on polarimetric maps made with HAWC+/SOFIA toward Rho Oph A, the densest portion of the Rho Ophiuchi molecular complex. We employed HAWC+ bands C (89 $\mu$m) and D (154 $\mu$m). The slope of the polarization spectrum was investigated by defining the quantity R_DC = p_D/p_C, where p_C and p_D represent polarization degrees in bands C and D, respectively. We find a clear correlation between R_DC and the molecular hydrogen column density across the cloud. A positive slope (R_DC > 1) dominates the lower density and well illuminated portions of the cloud, that are heated by the high mass star Oph S1, whereas a transition to a negative slope (R_DC < 1) is observed toward the denser and less evenly illuminated cloud core. We interpret the trends as due to a combination of: (1) Warm grains at the cloud outskirts, which are efficiently aligned by the abundant exposure to radiation from Oph S1, as proposed in the radiative torques theory; and (2) Cold grains deep in the cloud core, which are poorly aligned due to shielding from external radiation. To assess this interpretation, we developed a very simple toy model using a spherically symmetric cloud core based on Herschel data, and verified that the predicted variation of R_DC is consistent with the observations. This result introduces a new method that can be used to probe the grain alignment efficiency in molecular clouds, based on the analysis of trends in the far-infrared polarization spectrum.

We present the Census of the Local Universe (CLU) narrow-band survey to search for emission-line (\ha) galaxies. CLU-\ha~has imaged $\approx$3$\pi$ of the sky (26,470~deg$^2$) with 4 narrow-band filters that probe a distance out to 200~Mpc. We have obtained spectroscopic follow-up for galaxy candidates in 14 preliminary fields (101.6~deg$^2$) to characterize the limits and completeness of the survey. In these preliminary fields, CLU can identify emission lines down to an \ha~flux limit of $10^{-14}$~$\rm{erg~s^{-1}~cm^{-2}}$ at 90\% completeness, and recovers 83\% (67\%) of the \ha~flux from catalogued galaxies in our search volume at the $\Sigma$=2.5 ($\Sigma$=5) color excess levels. The contamination from galaxies with no emission lines is 61\% (12\%) for $\Sigma$=2.5 ($\Sigma$=5). Also, in the regions of overlap between our preliminary fields and previous emission-line surveys, we recover the majority of the galaxies found in previous surveys and identify an additional $\approx$300 galaxies. In total, we find 90 galaxies with no previous distance information, several of which are interesting objects: 7 blue compact dwarfs, 1 green pea, and a Seyfert galaxy; we also identified a known planetary nebula. These objects show that the CLU-\ha~survey can be a discovery machine for objects in our own Galaxy and extreme galaxies out to intermediate redshifts. However, the majority of the CLU-\ha~galaxies identified in this work show properties consistent with normal star-forming galaxies. CLU-\ha~galaxies with new redshifts will be added to existing galaxy catalogs to focus the search for the electromagnetic counterpart to gravitational wave events.

We present a detailed study of the flocculent spiral galaxy NGC 7793, part of the Sculptor group. By analyzing the resolved stellar populations of the galaxy, located at a distance of ~3.7 Mpc, we infer for the first time its radial star formation history (SFH) from Hubble Space Telescope photometry, thanks to both archival and new data from the Legacy ExtraGalactic UV Survey. We determine an average star formation rate (SFR) for the galaxy portion covered by our F555W and F814W data of 0.23 +- 0.02 Msun/yr over the whole Hubble time, corresponding to a total stellar mass of 3.09 +- 0.33 x 10^9 Msun in agreement with previous determinations. Thanks to the new data extending to the F336W band, we are able to analyze the youngest stellar populations with a higher time resolution. Most importantly, we recover the resolved SFH in different radial regions of the galaxy; this shows an indication of a growing trend of the present-to-past SFR ratio, increasing from internal to more external regions, supporting previous findings of the inside-out growth of the galaxy.

The survey description and near-infrared properties for 92 galaxies are presented for the Extended Disc Galaxy Exploration Science (EDGES) Survey, along with an investigation into the properties of the stellar halos of these galaxies. EDGES is a Spitzer Space Telescope Warm Mission program designed to reach the surface brightness limit ($\sim$0.5 kJy sr$^{-1}$ or 29 AB mag arcsec$^{-2}$) of the Infrared Array Camera (IRAC) 3.6 and 4.5 $\mu$m bands for a wide range of galaxy types found within the local volume. The surface brightness profiles exhibit a large range in disc scalelength, with breaks more frequently seen than in previous studies, owing in large part to the extremely deep near-infrared imaging. A number of these surface brightness profile breaks may be due to stellar halos, up to 7 galaxies out of the full sample of 92 galaxies, and we explore these implications in relation to current cosmological models. We also report the discovery of a new tidal stream near NGC 3953.

How do stars manage to form within low-density, HI-dominated gas? Such environments provide a laboratory for studying star formation with physical conditions distinct from starbursts and the metal-rich disks of spiral galaxies where most effort has been invested. Here we outline fundamental open questions about the nature of star formation at low-density. We describe the wide-field, high-resolution UV-optical-IR-radio observations of stars, star clusters and gas clouds in nearby galaxies needed in the 2020's to provide definitive answers, essential for development of a complete theory of star formation.

We present the star cluster catalogs for 17 dwarf and irregular galaxies in the $HST$ Treasury Program "Legacy ExtraGalactic UV Survey" (LEGUS). Cluster identification and photometry in this subsample are similar to that of the entire LEGUS sample, but special methods were developed to provide robust catalogs with accurate fluxes due to low cluster statistics. The colors and ages are largely consistent for two widely used aperture corrections, but a significant fraction of the clusters are more compact than the average training cluster. However, the ensemble luminosity, mass, and age distributions are consistent suggesting that the systematics between the two methods are less than the random errors. When compared with the clusters from previous dwarf galaxy samples, we find that the LEGUS catalogs are more complete and provide more accurate total fluxes. Combining all clusters into a composite dwarf galaxy, we find that the luminosity and mass functions can be described by a power law with the canonical index of $-2$ independent of age and global SFR binning. The age distribution declines as a power law, with an index of $\approx-0.80\pm0.15$, independent of cluster mass and global SFR binning. This decline of clusters is dominated by cluster disruption since the combined star formation histories and integrated-light SFRs are both approximately constant over the last few hundred Myr. Finally, we find little evidence for an upper-mass cutoff ($<2\sigma$) in the composite cluster mass function, and can rule out a truncation mass below $\approx10^{4.5}$M$_{\odot}$ but cannot rule out the existence of a truncation at higher masses.

We report new polarimetric and photometric maps of the massive star-forming region OMC-1 using the HAWC+ instrument on the Stratospheric Observatory for Infrared Astronomy (SOFIA). We present continuum polarimetric and photometric measurements of this region at 53, 89, 154, and 214 microns at angular resolutions of 5.1, 7.9, 14.0, and 18.7 arcseconds for the four bands, respectively. The photometric maps enable the computation of improved SEDs for the region. We find that at the longer wavelengths, the inferred magnetic field configuration matches the `hourglass' configuration seen in previous studies, indicating magnetically-regulated star formation. The field morphology differs at the shorter wavelengths. The magnetic field inferred at these wavelengths traces the bipolar structure of the explosive Becklin-Neugebauer (BN)/Kleinman-Low (KL) outflow emerging from OMC-1 behind the Orion Nebula. Using statistical methods to estimate the field strength in the region, we find that the explosion dominates the magnetic field near the center of the feature. Farther out, the magnetic field is close to energetic equilibrium with the ejecta and may be providing confinement to the explosion. The correlation between polarization fraction and the local polarization angle dispersion indicates that the depolarization as a function of unpolarized intensity is a result of intrinsic field geometry as opposed to decreases in grain alignment efficiency in denser regions.

We present a study correlating the spatial locations of young star clusters with those of molecular clouds in NGC~5194, in order to investigate the timescale over which clusters separate from their birth clouds. The star cluster catalogues are from the Legacy ExtraGalactic UV Survey (LEGUS) and the molecular clouds from the Plateau de Bure Interefrometer Arcsecond Whirpool Survey (PAWS). We find that younger star clusters are spatially closer to molecular clouds than older star clusters. The median ages for clusters associated with clouds is 4~Myr whereas it is 50~Myr for clusters that are sufficiently separated from a molecular cloud to be considered unassociated. After $\sim$6~Myr, the majority of the star clusters lose association with their molecular gas. Younger star clusters are also preferentially located in stellar spiral arms where they are hierarchically distributed in kpc-size regions for 50-100~Myr before dispersing. The youngest star clusters are more strongly clustered, yielding a two-point correlation function with $\alpha=-0.28\pm0.04$, than the GMCs ($\alpha=-0.09\pm0.03$) within the same PAWS field. However, the clustering strength of the most massive GMCs, supposedly the progenitors of the young clusters for a star formation efficiency of a few percent, is comparable ($\alpha=-0.35\pm0.05$) to that of the clusters. We find a galactocentric-dependence for the coherence of star formation, in which clusters located in the inner region of the galaxy reside in smaller star-forming complexes and display more homogeneous distributions than clusters further from the centre. This result suggests a correlation between the survival of a cluster complex and its environment.

We present new multi-broadband transit photometry of HD 189733b obtained at the Wyoming Infrared Observatory. Using an ensemble of five-band Sloan filter observations across multiple transits we have created an "ultra-low" resolution transmission spectrum to discern the nature of the exoplanet atmosphere. The observations were taken over three transit events and total 108 u', 120 g', 120 r', 110 i', and 116 z' images with an average exposure cadence of seven minutes for an entire series. The analysis was performed with a Markov-Chain Monte-Carlo method assisted by a Gaussian processes regression model. We find the apparent planet radius to increase from 0.154 +0.000920-0.00096 R* at z'-band to 0.157 +0.00074-0.00078 R* at u'-band. Whether this apparent radius implies an enhanced Rayleigh scattering or clear or grey planet atmosphere is highly dependent on stellar spot modeling assumptions, but our results are consistent with the literature for HD 189733b. This set of observations demonstrates the ability of our 2.3-m ground-based observatory to measure atmospheres of large exoplanets.

We have fit the far-ultraviolet (FUV) to sub-millimeter (850 micron) spectral energy distributions (SEDs) of the 61 galaxies from the "Key Insights on Nearby Galaxies: A Far-Infrared Survey with Herschel" (KINGFISH). The fitting has been performed using three models: the Code for Investigating GALaxy Evolution (CIGALE), the GRAphite-SILicate approach (GRASIL), and the Multi-wavelength Analysis of Galaxy PHYSical properties (MAGPHYS). We have analyzed the results of the three codes in terms of the SED shapes, and by comparing the derived quantities with simple "recipes" for stellar mass (Mstar), star-formation rate (SFR), dust mass (Mdust), and monochromatic luminosities. Although the algorithms rely on different assumptions for star-formation history, dust attenuation and dust reprocessing, they all well approximate the observed SEDs and are in generally good agreement for the associated quantities. However, the three codes show very different behavior in the mid-infrared regime, in particular between 25 and 70 micron where there are no observational constraints for the KINGFISH sample. We find that different algorithms give discordant SFR estimates for galaxies with low specific SFR, and that the standard "recipes" for calculating FUV absorption overestimate the extinction compared to the SED-fitting results. Results also suggest that assuming a "standard" constant stellar mass-to-light ratio overestimates Mstar relative to the SED fitting, and we provide new SED-based formulations for estimating Mstar from WISE W1 (3.4 micron) luminosities and colors. From a Principal Component Analysis of Mstar, SFR, Mdust, and O/H, we reproduce previous scaling relations among Mstar, SFR, and O/H, and find that Mdust can be predicted to within roughly 0.3 dex using only Mstar and SFR.

The science case and associated science requirements for a next-generation Very Large Array (ngVLA) are described, highlighting the five key science goals developed out of a community-driven vision of the highest scientific priorities in the next decade. Building on the superb cm observing conditions and existing infrastructure of the VLA site in the U.S. Southwest, the ngVLA is envisaged to be an interferometric array with more than 10 times the sensitivity and spatial resolution of the current VLA and ALMA, operating at frequencies spanning $\sim1.2 - 116$\,GHz with extended baselines reaching across North America. The ngVLA will be optimized for observations at wavelengths between the exquisite performance of ALMA at submm wavelengths, and the future SKA-1 at decimeter to meter wavelengths, thus lending itself to be highly complementary with these facilities. The ngVLA will be the only facility in the world that can tackle a broad range of outstanding scientific questions in modern astronomy by simultaneously delivering the capability to: (1) unveil the formation of Solar System analogues; (2) probe the initial conditions for planetary systems and life with astrochemistry; (3) characterize the assembly, structure, and evolution of galaxies from the first billion years to the present; (4) use pulsars in the Galactic center as fundamental tests of gravity; and (5) understand the formation and evolution of stellar and supermassive blackholes in the era of multi-messenger astronomy.

We present an investigation of the relationship between giant molecular cloud (GMC) properties and the associated stellar clusters in the nearby flocculent galaxy NGC 7793. We combine the star cluster catalog from the HST LEGUS (Legacy ExtraGalactic UV Survey) program with the 15 parsec resolution ALMA CO(2-1) observations. We find a strong spatial correlation between young star clusters and GMCs such that all clusters still associated with a GMC are younger than 11 Myr and display a median age of 2 Myr. The age distribution increases gradually as the cluster-GMC distance increases, with star clusters that are spatially unassociated with molecular gas exhibiting a median age of 7 Myr. Thus, star clusters are able to emerge from their natal clouds long before the timescale required for clouds to disperse. To investigate if the hierarchy observed in the stellar components is inherited from the GMCs, we quantify the amount of clustering in the spatial distributions of the components and find that the star clusters have a fractal dimension slope of $-0.35 \pm 0.03$, significantly more clustered than the molecular cloud hierarchy with slope of $-0.18 \pm 0.04$ over the range 40-800 pc. We find, however, that the spatial clustering becomes comparable in strength for GMCs and star clusters with slopes of $-0.44\pm0.03$ and $-0.45\pm0.06$ respectively, when we compare massive ($>$10$^5$ M$_{\odot}$) GMCs to massive and young star clusters. This shows that massive star clusters trace the same hierarchy as their parent GMCs, under the assumption that the star formation efficiency is a few percent.

We have explored the role environmental factors play in determining characteristics of young stellar objects in nearby dwarf irregular and Blue Compact Dwarf galaxies. Star clusters are characterized by concentrations, masses, and formation rates, OB associations by mass and mass surface density, O stars by their numbers and near-ultraviolet absolute magnitudes, and HII regions by Halpha surface brightnesses. These characteristics are compared to surrounding galactic pressure, stellar mass density, HI surface density, and star formation rate surface density. We find no trend of cluster characteristics with environmental properties, implying that larger scale effects are more important in determining cluster characteristics or that rapid dynamical evolution erases memory of the initial conditions. On the other hand, the most massive OB associations are found at higher pressure and HI surface density, and there is a trend of higher HII region Halpha surface brightness with higher pressure, suggesting that a higher concentration of massive stars and gas are found preferentially in regions of higher pressure. At low pressures we find massive stars but not bound clusters and OB associations. We do not find evidence for an increase of cluster formation efficiency as a function of star formation rate density. However, there is an increase in the ratio of the number of clusters to number of O stars with pressure, perhaps reflecting an increase in clustering properties with star formation rate.

The traditional picture of post-starburst galaxies as dust- and gas-poor merger remnants, rapidly transitioning to quiescence, has been recently challenged. Unexpected detections of a significant ISM in many post-starbursts raise important questions. Are they truly quiescent and, if so, what mechanisms inhibit further star formation? What processes dominate their ISM energetics? We present an infrared spectroscopic and photometric survey of 33 SDSS-selected E+A post-starbursts, aimed at resolving these questions. We find compact, warm dust reservoirs with high PAH abundances, and total gas and dust masses significantly higher than expected from stellar recycling alone. Both PAH/TIR and dust-to-burst stellar mass ratios are seen to decrease with post-burst age, indicative of the accumulating effects of dust destruction and an incipient transition to hot, early-type ISM properties. Their infrared spectral properties are unique, with dominant PAH emission, very weak nebular lines, unusually strong H$_{2}$ rotational emission, and deep ${\rm [C\, II]}$ deficits. There is substantial scatter among SFR indicators, and both PAH and TIR luminosities provide overestimates. Even as potential upper limits, all tracers show that the SFR has typically experienced a more than two order-of-magnitude decline since the starburst, and that the SFR is considerably lower than expected given both their stellar masses and molecular gas densities. These results paint a coherent picture of systems in which star formation was, indeed, rapidly truncated, but in which the ISM was $\textit{not}$ completely expelled, and is instead supported against collapse by latent or continued injection of turbulent or mechanical heating. The resulting aging burst populations provide a "high-soft" radiation field which seemingly dominates the E+As' unusual ISM energetics.

We present accurate resolved $WISE$ photometry of galaxies in the combined SINGS and KINGFISH sample. The luminosities in the W3 12$\mu$m and W4 23$\mu$m bands are calibrated to star formation rates (SFRs) derived using the total infrared luminosity, avoiding UV/optical uncertainties due to dust extinction corrections. The W3 relation has a 1-$\sigma$ scatter of 0.15 dex over nearly 5 orders of magnitude in SFR and 12$\mu$m luminosity, and a range in host stellar mass from dwarf (10$^7$ M$_\odot$) to $\sim3\times$M$_\star$ (10$^{11.5}$ M$_\odot$) galaxies. In the absence of deep silicate absorption features and powerful active galactic nuclei, we expect this to be a reliable SFR indicator chiefly due to the broad nature of the W3 band. By contrast the W4 SFR relation shows more scatter (1-$\sigma =$ 0.18 dex). Both relations show reasonable agreement with radio continuum-derived SFRs and excellent accordance with so-called "hybrid" H$\alpha + 24 \mu$m and FUV$+$24$\mu$m indicators. Moreover, the $WISE$ SFR relations appear to be insensitive to the metallicity range in the sample. We also compare our results with IRAS-selected luminous infrared galaxies, showing that the $WISE$ relations maintain concordance, but systematically deviate for the most extreme galaxies. Given the all-sky coverage of $WISE$ and the performance of the W3 band as a SFR indicator, the $L_{12\mu \rm m}$ SFR relation could be of great use to studies of nearby galaxies and forthcoming large area surveys at optical and radio wavelengths.

One of the main theories for explaining the formation of spiral arms in galaxies is the stationary density wave theory. This theory predicts the existence of an age gradient across the arms. We use the stellar cluster catalogues of the galaxies NGC 1566, M51a, and NGC 628 from the Legacy Extragalactic UV Survey (LEGUS) program. In order to test for the possible existence of an age sequence across the spiral arms, we quantified the azimuthal offset between star clusters of different ages in our target galaxies. We found that NGC 1566, a grand-design spiral galaxy with bisymmetric arms and a strong bar, shows a significant age gradient across the spiral arms that appears to be consistent with the prediction of the stationary density wave theory. In contrast, M51a with its two well-defined spiral arms and a weaker bar does not show an age gradient across the arms. In addition, a comparison with non LEGUS star cluster catalogues for M51a yields similar results. We believe that the spiral structure of M51a is not the result of a stationary density wave with a fixed pattern speed. Instead, tidal interactions could be the dominant mechanism for the formation of spiral arms. We also found no offset in the azimuthal distribution of star clusters with different ages across the weak spiral arms of NGC 628.

It has recently been established that the properties of young star clusters (YSCs) can vary as a function of the galactic environment in which they are found. We use the cluster catalogue produced by the Legacy Extragalactic UV Survey (LEGUS) collaboration to investigate cluster properties in the spiral galaxy M51. We analyse the cluster population as a function of galactocentric distance and in arm and inter-arm regions. The cluster mass function exhibits a similar shape at all radial bins, described by a power law with a slope close to $-2$ and an exponential truncation around $10^5\ \rm{M}_{\odot}$ . While the mass functions of the YSCs in the spiral arm and inter-arm regions have similar truncation masses, the inter-arm region mass function has a significantly steeper slope than the one in the arm region; a trend that is also observed in the giant molecular cloud mass function and predicted by simulations. The age distribution of clusters is dependent on the region considered, and is consistent with rapid disruption only in dense regions, while little disruption is observed at large galactocentric distances and in the inter-arm region. The fraction of stars forming in clusters does not show radial variations, despite the drop in the $H_2$ surface density measured as function of galactocentric distance. We suggest that the higher disruption rate observed in the inner part of the galaxy is likely at the origin of the observed flat cluster formation efficiency radial profile.

We present a detailed study of the Magellanic irregular galaxy NGC 4449 based on both archival and new photometric data from the Legacy Extragalactic UV Survey, obtained with the Hubble Space Telescope Advanced Camera for Surveys and Wide Field Camera 3. Thanks to its proximity ($D=3.82\pm 0.27$ Mpc) we reach stars 3 magnitudes fainter than the tip of the red giant branch in the F814W filter. The recovered star formation history spans the whole Hubble time, but due to the age-metallicity degeneracy of the red giant branch stars, it is robust only over the lookback time reached by our photometry, i.e. $\sim 3$ Gyr. The most recent peak of star formation is around 10 Myr ago. The average surface density star formation rate over the whole galaxy lifetime is $0.01$ M$_{\odot}$ yr$^{-1}$ kpc$^{-2}$. From our study it emerges that NGC 4449 has experienced a fairly continuous star formation regime in the last 1 Gyr with peaks and dips whose star formation rates differ only by a factor of a few. The very complex and disturbed morphology of NGC 4449 makes it an interesting galaxy for studies of the relationship between interactions and starbursts, and our detailed and spatially resolved analysis of its star formation history does indeed provide some hints on the connection between these two phenomena in this peculiar dwarf galaxy.

We present a study of the dust-to-gas ratios in five nearby galaxies NGC 628 (M74), NGC 6503, NGC 7793, UGC 5139 (Holmberg I), and UGC 4305 (Holmberg II). Using Hubble Space Telescope broad band WFC3/UVIS UV and optical images from the Treasury program LEGUS (Legacy ExtraGalactic UV Survey) combined with archival HST/ACS data, we correct thousands of individual stars for extinction across these five galaxies using an isochrone-matching (reddening-free Q) method. We generate extinction maps for each galaxy from the individual stellar extinctions using both adaptive and fixed resolution techniques, and correlate these maps with neutral HI and CO gas maps from literature, including The HI Nearby Galaxy Survey (THINGS) and the HERA CO-Line Extragalactic Survey (HERACLES). We calculate dust-to-gas ratios and investigate variations in the dust-to-gas ratio with galaxy metallicity. We find a power law relationship between dust-to-gas ratio and metallicity, consistent with other studies of dust-to-gas ratio compared to metallicity. We find a change in the relation when H$_2$ is not included. This implies that underestimation of $N_{H_2}$ in low-metallicity dwarfs from a too-low CO-to-H$_2$ conversion factor $X_{CO}$ could have produced too low a slope in the derived relationship between dust-to-gas ratio and metallicity. We also compare our extinctions to those derived from fitting the spectral energy distribution (SED) using the Bayesian Extinction and Stellar Tool (BEAST) for NGC 7793 and find systematically lower extinctions from SED-fitting as compared to isochrone matching.

We use HST observations from the Legacy Extragalactic UV Survey to reconstruct the recent star formation histories (SFHs) of three actively star-forming dwarf galaxies, NGC4449, Holmberg II and NGC1705, from their UV color-magnitude diagrams (CMDs). We apply a CMD fitting technique using two independent sets of stellar isochrones, PARSEC-COLIBRI and MIST, to assess the uncertainties related to stellar evolution modelling. Irrespective of the adopted stellar models, all the three dwarfs are found to have had almost constant star formation rates (SFRs) in the last 100-200 Myr, with modest enhancements (a factor of $\sim$2) above the 100 Myr-averaged-SFR. Significant differences among the three dwarfs are found in the overall SFR, the timing of the most recent peak and the SFR$/$area. The Initial Mass Function (IMF) of NGC1705 and Holmberg II is consistent with a Salpeter slope down to $\approx$ 5 M$_{\odot}$, whereas it is slightly flatter, s$=-2.0$, in NGC4449. The SFHs derived with the two different sets of stellar models are consistent with each other, except for some quantitative details, attributable to their input assumptions. They also share the drawback that all synthetic diagrams predict a clear separation in color between upper main sequence and helium burning stars, which is not apparent in the data. Since differential reddening, significant in NGC4449, or unresolved binaries don't appear to be sufficient to fill the gap, we suggest this calls for a revision of both sets of stellar evolutionary tracks.

Two dwarf irregular galaxies DDO 187 and NGC 3738 exhibit a striking pattern of star formation: intense star formation is taking place in a large region occupying roughly half of the inner part of the optical galaxy. We use data on the HI distribution and kinematics and stellar images and colors to examine the properties of the environment in the high star formation rate (HSF) halves of the galaxies in comparison with the low star formation rate (LSF) halves. We find that the pressure and gas density are higher on the HSF sides by 30-70%. In addition we find in both galaxies that the HI velocity fields exhibit significant deviations from ordered rotation and there are large regions of high velocity dispersion and multiple velocity components in the gas beyond the inner regions of the galaxies. The conditions in the HSF regions are likely the result of large-scale external processes affecting the internal environment of the galaxies and enabling the current star formation there.

The Legacy ExtraGalactic UV Survey (LEGUS) is a multiwavelength Cycle 21 Treasury program on the Hubble Space Telescope. It studied 50 nearby star-forming galaxies in five bands from the near UV to the I-band, combining new Wide Field Camera 3 observations with archival Advanced Camera for Surveys data. LEGUS was designed to investigate how star formation occurs and develops on both small and large scales, and how it relates to the galactic environments. In this paper we present the photometric catalogs for all the apparently single stars identified in the 50 LEGUS galaxies. Photometric catalogs and mosaicked images for all filters are available for download. We present optical and near UV color-magnitude diagrams for all the galaxies. For each galaxy we derived the distance from the tip of the red giant branch. We then used the NUV color-magnitude diagrams to identify stars more massive than 14 Mo, and compared their number with the number of massive stars expected from the GALEX FUV luminosity. Our analysis shows that the fraction of massive stars forming in star clusters and stellar associations is about constant with the star formation rate. This lack of a relation suggests that the time scale for evaporation of unbound structures is comparable or longer than 10 Myr. At low star formation rates this translates to an excess of mass in clustered environments as compared to model predictions of cluster evolution, suggesting that a significant fraction of stars form in unbound systems.

Recently acquired WFC3 UV (F275W and F336W) imaging mosaics under the Legacy Extragalactic UV Survey (LEGUS) combined with archival ACS data of M51 are used to study the young star cluster (YSC) population of this interacting system. Our newly extracted source catalogue contains 2834 cluster candidates, morphologically classified to be compact and uniform in colour, for which ages, masses and extinction are derived. In this first work we study the main properties of the YSC population of the whole galaxy, considering a mass-limited sample. Both luminosity and mass functions follow a power law shape with slope -2, but at high luminosities and masses a dearth of sources is observed. The analysis of the mass function suggests that it is best fitted by a Schechter function with slope -2 and a truncation mass at $1.00\pm0.12\times10^5\ M_\odot$. Through Monte Carlo simulations we confirm this result and link the shape of the luminosity function to the presence of a truncation in the mass function. A mass limited age function analysis, between 10 and 200 Myr, suggests that the cluster population is undergoing only moderate disruption. We observe little variation in the shape of the mass function at masses above $1\times10^4\ M_\odot$, over this age range. The fraction of star formation happening in the form of bound clusters in M51 is $\sim20\%$ in the age range 10 to 100 Myr and little variation is observed over the whole range from 1 to 200 Myr.

The [CII] 158um fine-structure line is the brightest emission line observed in local star-forming galaxies. As a major coolant of the gas-phase interstellar medium, [CII] balances the heating, including that due to far-ultraviolet photons, which heat the gas via the photoelectric effect. However, the origin of [CII] emission remains unclear, because C+ can be found in multiple phases of the interstellar medium. Here we measure the fractions of [CII] emission originating in the ionized and neutral gas phases of a sample of nearby galaxies. We use the [NII] 205um fine-structure line to trace the ionized medium, thereby eliminating the strong density dependence that exists in the ratio of [CII]/[NII] 122um. Using the FIR [CII] and [NII] emission detected by the KINGFISH and Beyond the Peak Herschel programs, we show that 60-80% of [CII] emission originates from neutral gas. We find that the fraction of [CII] originating in the neutral medium has a weak dependence on dust temperature and the surface density of star formation, and a stronger dependence on the gas-phase metallicity. In metal-rich environments, the relatively cooler ionized gas makes substantially larger contributions to total [CII] emission than at low abundance, contrary to prior expectations. Approximate calibrations of this metallicity trend are provided.

At $z=1-3$, the formation of new stars is dominated by dusty galaxies whose far-IR emission indicates they contain colder dust than local galaxies of a similar luminosity. We explore the reasons for the evolving IR emission of similar galaxies over cosmic time using: 1) Local galaxies from GOALS $(L_{\rm IR}=10^{11}-10^{12}\,L_\odot)$; 2) Galaxies at $z\sim0.1-0.5$ from the 5MUSES ($L_{\rm IR}=10^{10}-10^{12}\,L_\odot$); 3) IR luminous galaxies spanning $z=0.5-3$ from GOODS and Spitzer xFLS ($L_{\rm IR}>10^{11}\,L_\odot$). All samples have Spitzer mid-IR spectra, and Herschel and ground-based submillimeter imaging covering the full IR spectral energy distribution, allowing us to robustly measure $L_{\rm IR}^{\rm\scriptscriptstyle SF}$, $T_{\rm dust}$, and $M_{\rm dust}$ for every galaxy. Despite similar infrared luminosities, $z>0.5$ dusty star forming galaxies have a factor of 5 higher dust masses and 5K colder temperatures. The increase in dust mass is linked with an increase in the gas fractions with redshift, and we do not observe a similar increase in stellar mass or star formation efficiency. $L_{160}^{\rm\scriptscriptstyle SF}/L_{70}^{\rm\scriptscriptstyle SF}$, a proxy for $T_{\rm dust}$, is strongly correlated with $L_{\rm IR}^{\rm\scriptscriptstyle SF}/M_{\rm dust}$ independently of redshift. We measure merger classification and galaxy size for a subsample, and there is no obvious correlation between these parameters and $L_{\rm IR}^{\rm \scriptscriptstyle SF}/M_{\rm dust}$ or $L_{160}^{\rm\scriptscriptstyle SF}/L_{70}^{\rm\scriptscriptstyle SF}$. In dusty star forming galaxies, the change in $L_{\rm IR}^{\rm\scriptscriptstyle SF}/M_{\rm dust}$ can fully account for the observed colder dust temperatures, suggesting that any change in the spatial extent of the interstellar medium is a second order effect.

We present an analysis of the positions and ages of young star clusters in eight local galaxies to investigate the connection between the age difference and separation of cluster pairs. We find that star clusters do not form uniformly but instead are distributed such that the age difference increases with the cluster pair separation to the 0.25-0.6 power, and that the maximum size over which star formation is physically correlated ranges from ~200 pc to ~1 kpc. The observed trends between age difference and separation suggest that cluster formation is hierarchical both in space and time: clusters that are close to each other are more similar in age than clusters born further apart. The temporal correlations between stellar aggregates have slopes that are consistent with turbulence acting as the primary driver of star formation. The velocity associated with the maximum size is proportional to the galaxy's shear, suggesting that the galactic environment influences the maximum size of the star-forming structures.

We present the implementation of a Bayesian formalism within the Stochastically Lighting Up Galaxies (SLUG) stellar population synthesis code, which is designed to investigate variations in the initial mass function (IMF) of star clusters. By comparing observed cluster photometry to large libraries of clusters simulated with a continuously varying IMF, our formalism yields the posterior probability distribution function (PDF) of the cluster mass, age, and extinction, jointly with the parameters describing the IMF. We apply this formalism to a sample of star clusters from the nearby galaxy NGC 628, for which broad-band photometry in five filters is available as part of the Legacy ExtraGalactic UV Survey (LEGUS). After allowing the upper-end slope of the IMF ($\alpha_3$) to vary, we recover PDFs for the mass, age, and extinction that are broadly consistent with what is found when assuming an invariant Kroupa IMF. However, the posterior PDF for $\alpha_3$ is very broad due to a strong degeneracy with the cluster mass, and it is found to be sensitive to the choice of priors, particularly on the cluster mass. We find only a modest improvement in the constraining power of $\alpha_3$ when adding H$\alpha$ photometry from the companion H$\alpha$-LEGUS survey. Conversely, H$\alpha$ photometry significantly improves the age determination, reducing the frequency of multi-modal PDFs. With the aid of mock clusters we quantify the degeneracy between physical parameters, showing how constraints on the cluster mass that are independent of photometry can be used to pin down the IMF properties of star clusters.

We present a study of the effective (half-light) radii and other structural properties of a systematically selected sample of young, massive star clusters (YMCs, $\geq$$5\times10^3$ M$_{\odot}$ and $\leq$200 Myr) in two nearby spiral galaxies, NGC 628 and NGC 1313. We use Hubble Space Telescope WFC3/UVIS and archival ACS/WFC data obtained by the Legacy Extragalactic UV Survey (LEGUS), an HST Treasury Program. We measure effective radii with GALFIT, a two-dimensional image-fitting package, and with a new technique to estimate effective radii from the concentration index (CI) of observed clusters. The distribution of effective radii from both techniques spans $\sim$0.5-10 pc and peaks at 2-3 pc for both galaxies. We find slight positive correlations between effective radius and cluster age in both galaxies, but no significant relationship between effective radius and galactocentric distance. Clusters in NGC 1313 display a mild increase in effective radius with cluster mass, but the trend disappears when the sample is divided into age bins. We show that the vast majority of the clusters in both galaxies are much older than their dynamical times, suggesting they are gravitationally bound objects. We find that about half of the clusters in NGC 628 are underfilling their Roche lobes, based on their Jacobi radii. Our results suggest that the young, massive clusters in NGC 628 and NGC 1313 are expanding due to stellar mass loss or two-body relaxation and are not significantly influenced by the tidal fields of their host galaxies.

We report the large effort which is producing comprehensive high-level young star cluster (YSC) catalogues for a significant fraction of galaxies observed with the Legacy ExtraGalactic UV Survey (LEGUS) Hubble treasury program. We present the methodology developed to extract cluster positions, verify their genuine nature, produce multiband photometry (from NUV to NIR), and derive their physical properties via spectral energy distribution fitting analyses. We use the nearby spiral galaxy NGC628 as a test case for demonstrating the impact that LEGUS will have on our understanding of the formation and evolution of YSCs and compact stellar associations within their host galaxy. Our analysis of the cluster luminosity function from the UV to the NIR finds a steepening at the bright end and at all wavelengths suggesting a dearth of luminous clusters. The cluster mass function of NGC628 is consistent with a power-law distribution of slopes $\sim -2$ and a truncation of a few times $10^5$ M$_\odot$. After their formation YSCs and compact associations follow different evolutionary paths. YSCs survive for a longer timeframe, confirming their being potentially bound systems. Associations disappear on time scales comparable to hierarchically organized star-forming regions, suggesting that they are expanding systems. We find mass-independent cluster disruption in the inner region of NGC628, while in the outer part of the galaxy there is little or no disruption. We observe faster disruption rates for low mass ($\leq$ $10^4$ M$_\odot$) clusters suggesting that a mass-dependent component is necessary to fully describe the YSC disruption process in NGC628.

We present a study of the hierarchical clustering of the young stellar clusters in six local (3--15 Mpc) star-forming galaxies using Hubble Space Telescope broad band WFC3/UVIS UV and optical images from the Treasury Program LEGUS (Legacy ExtraGalactic UV Survey). We have identified 3685 likely clusters and associations, each visually classified by their morphology, and we use the angular two-point correlation function to study the clustering of these stellar systems. We find that the spatial distribution of the young clusters and associations are clustered with respect to each other, forming large, unbound hierarchical star-forming complexes that are in general very young. The strength of the clustering decreases with increasing age of the star clusters and stellar associations, becoming more homogeneously distributed after ~40--60 Myr and on scales larger than a few hundred parsecs. In all galaxies, the associations exhibit a global behavior that is distinct and more strongly correlated from compact clusters. Thus, populations of clusters are more evolved than associations in terms of their spatial distribution, traveling significantly from their birth site within a few tens of Myr whereas associations show evidence of disruption occurring very quickly after their formation. The clustering of the stellar systems resembles that of a turbulent interstellar medium that drives the star formation process, correlating the components in unbound star-forming complexes in a hierarchical manner, dispersing shortly after formation, suggestive of a single, continuous mode of star formation across all galaxies.

We investigate how star formation is spatially organized in the grand-design spiral NGC 1566 from deep HST photometry with the Legacy ExtraGalactic UV Survey (LEGUS). Our contour-based clustering analysis reveals 890 distinct stellar conglomerations at various levels of significance. These star-forming complexes are organized in a hierarchical fashion with the larger congregations consisting of smaller structures, which themselves fragment into even smaller and more compact stellar groupings. Their size distribution, covering a wide range in length-scales, shows a power-law as expected from scale-free processes. We explain this shape with a simple "fragmentation and enrichment" model. The hierarchical morphology of the complexes is confirmed by their mass--size relation which can be represented by a power-law with a fractional exponent, analogous to that determined for fractal molecular clouds. The surface stellar density distribution of the complexes shows a log-normal shape similar to that for supersonic non-gravitating turbulent gas. Between 50 and 65 per cent of the recently-formed stars, as well as about 90 per cent of the young star clusters, are found inside the stellar complexes, located along the spiral arms. We find an age-difference between young stars inside the complexes and those in their direct vicinity in the arms of at least 10 Myr. This timescale may relate to the minimum time for stellar evaporation, although we cannot exclude the in situ formation of stars. As expected, star formation preferentially occurs in spiral arms. Our findings reveal turbulent-driven hierarchical star formation along the arms of a grand-design galaxy.

We present an update to the ultraviolet-to-radio database of global broadband photometry for the 79 nearby galaxies that comprise the union of the KINGFISH (Key Insights on Nearby Galaxies: A Far-Infrared Survey with Herschel) and SINGS (Spitzer Infrared Nearby Galaxies Survey) samples. The 34-band dataset presented here includes contributions from observational work carried out with a variety of facilities including GALEX, SDSS, PS, NOAO, 2MASS, WISE, Spitzer, Herschel, Planck, JCMT, and the VLA. Improvements of note include recalibrations of previously-published SINGS BVRcIc and KINGFISH far-infrared/submillimeter photometry. Similar to previous results in the literature, an excess of submillimeter emission above model predictions is seen primarily for low-metallicity dwarf/irregular galaxies. This 34-band photometric dataset for the combined KINGFISH$+$SINGS sample serves as an important multi-wavelength reference for the variety of galaxies observed at low redshift. A thorough analysis of the observed spectral energy distributions is carried out in a companion paper.

Dynamic and thermal processes regulate the structure of the multi-phase interstellar medium (ISM), and ultimately establish how galaxies evolve through star formation. Thus, to constrain ISM models and better understand the interplay of these processes, it is of great interest to measure the thermal pressure ($P_{\rm th}$) of the diffuse, neutral gas. By combining [C II] 158 $\mu$m, HI, and CO data from 31 galaxies selected from the Herschel KINGFISH sample, we have measured thermal pressures in 534 predominantly atomic regions with typical sizes of $\sim$1 kiloparsec. We find a distribution of thermal pressures in the $P_{\rm th}/k\sim10^3-10^5$ K cm$^{-3}$ range. For a sub-sample of regions with conditions similar to those of the diffuse, neutral gas in the Galactic plane, we find thermal pressures that follow a log-normal distribution with a median value of $P_{\rm th}/k\approx3600$ K cm$^{-3}$. These results are consistent with thermal pressure measurements using other observational methods. We find that $P_{\rm th}$ increases with radiation field strength and star formation activity, as expected from the close link between the heating of the gas and the star formation rate. Our thermal pressure measurements fall in the regime where a two-phase ISM with cold and warm neutral medium could exist in pressure equilibrium. Finally, we find that the midplane thermal pressure of the diffuse gas is about $\sim30$% of the vertical weight of the overlying ISM, consistent with results from hydrodynamical simulations of self-regulated star formation in galactic disks.

We study the spectral energy distribution (SED) of the radio continuum emission from the KINGFISH sample of nearby galaxies to understand the energetics and origin of this emission. Effelsberg multi-wavelength observations at 1.4GHz, 4.8GHz, 8.5GHz, and 10.5GHz combined with archive data allow us, for the first time, to determine the mid-radio continuum (1-10 GHz, MRC) bolometric luminosities and further present calibration relations vs. the monochromatic radio luminosities. The 1-10 GHz radio SED is fitted using a Bayesian Markov Chain Monte Carlo (MCMC) technique leading to measurements for the nonthermal spectral index and the thermal fraction f_th with mean values of alpha_nt=0.97+-0.16 (0.79+-0.15 for the total spectral index) and f_th= 10% +- 9% at 1.4 GHz. The MRC luminosity changes over ~3 orders of magnitude in the sample. The thermal emission is responsible for ~23% of the MRC on average. We also compare the extinction-corrected diagnostics of star formation rate with the thermal and nonthermal radio tracers and derive the first star formation calibration relations using the MRC radio luminosity. The nonthermal spectral index flattens with increasing star formation rate surface density, indicating the effect of the star formation feedback on the cosmic ray electron population in galaxies. Comparing the radio and IR SEDs, we find that the FIR-to-MRC ratio could decrease with star formation rate, due to the amplification of the magnetic fields in star forming regions. This particularly implies a decrease in the ratio at high redshifts, where mostly luminous/star forming galaxies are detected.

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.

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.

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.

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.

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$\mu$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.

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

To compute the SFR of galaxies from the rest-frame UV it is essential to take into account the obscuration by dust. To do so, one of the most popular methods consists in combining the UV with the emission from the dust itself in the IR. Yet, different studies have derived different estimators, showing that no such hybrid estimator is truly universal. In this paper we aim at understanding and quantifying what physical processes drive the variations between different hybrid estimators. Doing so, we aim at deriving new universal UV+IR hybrid estimators to correct the UV for dust attenuation, taking into account the intrinsic physical properties of galaxies. We use the CIGALE code to model the spatially-resolved FUV to FIR SED of eight nearby star-forming galaxies drawn from the KINGFISH sample. This allows us to determine their local physical properties, and in particular their UV attenuation, average SFR, average specific SFR (sSFR), and their stellar mass. We then examine how hybrid estimators depend on said properties. We find that hybrid UV+IR estimators strongly depend on the stellar mass surface density (in particular at 70 and 100 micron) and on the sSFR (in particular at 24 micron and the TIR). Consequently, the IR scaling coefficients for UV obscuration can vary by almost an order of magnitude. This result contrasts with other groups who found relatively constant coefficients with small deviations. We exploit these variations to construct a new class of hybrid estimators based on observed UV to near-IR colours and near-IR luminosity densities per unit area. We find that they can reliably be extended to entire galaxies. The new estimators provide better estimates of attenuation-corrected UV emission than classical hybrid estimators. Naturally taking into account the variable impact of dust heated by old stellar populations, they constitute a step towards universal estimators.

New deep optical and near-infrared imaging is combined with archival ultraviolet and infrared data for fifteen nearby galaxies mapped in the Spitzer Extended Disk Galaxy Exploration Science survey. These images are particularly deep and thus excellent for studying the low surface brightness outskirts of these disk-dominated galaxies with stellar masses ranging between 10^8 and 10^11 Msun. The spectral energy distributions derived from this dataset are modeled to investigate the radial variations in the galaxy colors and star formation histories. Taken as a whole, the sample shows bluer and younger stars for larger radii until reversing near the optical radius, whereafter the trend is for redder and older stars for larger galacto-centric distances. These results are consistent with an inside-out disk formation scenario coupled with an old stellar outer disk population formed through radial migration and/or the cumulative history of minor mergers and accretions of satellite dwarf galaxies. However, these trends are quite modest and the variation from galaxy to galaxy is substantial. Additional data for a larger sample of galaxies are needed to confirm or dismiss these modest sample-wide trends.

We present a study of the spatial distribution of the stellar cluster populations in the star forming galaxy NGC 628. Using Hubble Space Telescope broad band WFC3/UVIS UV and optical images from the Treasury Program LEGUS (Legacy ExtraGalactic UV Survey), we have identified 1392 potential young (<100 Myr) stellar clusters within the galaxy, identified from a combination of visual inspection and automatic selection. We investigate the clustering of these young stellar clusters and quantify the strength and change of clustering strength with scale using the two-point correlation function. We also investigate how image boundary conditions and dust lanes affect the observed clustering. The distribution of the clusters is well fit by a broken power law with negative exponent $\alpha$. We recover a weighted mean index of $\alpha$ ~ -0.8 for all spatial scales below the break at 3".3 (158 pc at a distance of 9.9 Mpc) and an index of $\alpha$ ~ -0.18 above 158 pc for the accumulation of all cluster types. The strength of the clustering increases with decreasing age and clusters older than 40 Myr lose their clustered structure very rapidly and tend to be randomly distributed in this galaxy whereas the mass of the star cluster has little effect on the clustering strength. This is consistent with results from other studies that the morphological hierarchy in stellar clustering resembles the same hierarchy as the turbulent interstellar medium.

We present new near-infrared (NIR) observations of M63 from the Extended Disk Galaxy Exploration Science (EDGES) Survey. The extremely deep 3.6 $\mu$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.

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.

We present a detailed clustering analysis of the young stellar population across the star-forming ring galaxy NGC 6503, based on the deep HST photometry obtained with the Legacy ExtraGalactic UV Survey (LEGUS). We apply a contour-based map analysis technique and identify in the stellar surface density map 244 distinct star-forming structures at various levels of significance. These stellar complexes are found to be organized in a hierarchical fashion with 95% being members of three dominant super-structures located along the star-forming ring. The size distribution of the identified structures and the correlation between their radii and numbers of stellar members show power-law behaviors, as expected from scale-free processes. The self-similar distribution of young stars is further quantified from their autocorrelation function, with a fractal dimension of ~1.7 for length-scales between ~20 pc and 2.5 kpc. The young stellar radial distribution sets the extent of the star-forming ring at radial distances between 1 and 2.5 kpc. About 60% of the young stars belong to the detected stellar structures, while the remaining stars are distributed among the complexes, still inside the ring of the galaxy. The analysis of the time-dependent clustering of young populations shows a significant change from a more clustered to a more distributed behavior in a time-scale of ~60 Myr. The observed hierarchy in stellar clustering is consistent with star formation being regulated by turbulence across the ring. The rotational velocity difference between the edges of the ring suggests shear as the driving mechanism for this process. Our findings reveal the interesting case of an inner ring forming stars in a hierarchical fashion.

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

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{\mu}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{\mu}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{\alpha} from the WIYN 0.9m telescope complement the 3.6{\mu}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.

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

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.

We use a sample of 36 galaxies from the KINGFISH (Herschel IR), HERACLES (IRAM CO), and THINGS (VLA HI) surveys to study empirical relations between Herschel infrared (IR) luminosities and the total mass of the interstellar gas (H2+HI). Such a comparison provides a simple empirical relationship without introducing the uncertainty of dust model fitting. We find tight correlations, and provide fits to these relations, between Herschel luminosities and the total gas mass integrated over entire galaxies, with the tightest, almost linear, correlation found for the longest wavelength data (SPIRE500). However, we find that accounting for the gas-phase metallicity (affecting the dust-to-gas ratio) is crucial when applying these relations to low-mass, and presumably high-redshift, galaxies. The molecular (H2) gas mass is found to be better correlated with the peak of the IR emission (e.g. PACS160), driven mostly by the correlation of stellar mass and mean dust temperature. When examining these relations as a function of galactocentric radius we find the same correlations, albeit with a larger scatter, up to a radius of 0.7 r_25 (within which most of the galaxy's baryonic mass resides). However, beyond this radius the same correlations no longer hold, with the gas mass (predominantly HI) increasing relative to the infrared emission. The tight relations found for the bulk of the galaxy's baryonic content suggest that the total gas masses of disk-like (non-merging) galaxies can be inferred from far-infrared continuum measurements in situations where only the latter are available, e.g. in ALMA continuum observations of high-redshift galaxies.

The [CII] 157.74 $\mu$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 $\Sigma_{\rm [CII]}-\Sigma_{\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 $\Sigma_{\rm[CII]}-\Sigma_{\rm SFR}$ correlation is valid over almost 5 orders of magnitude in $\Sigma_{\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.

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 \beta 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.

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{\AA}$ to 24 $\mu$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.

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.

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.

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.

We investigate the far-infrared (IR) dust emission for 20 local star forming galaxies from the Key Insights on Nearby Galaxies: A Far-IR Survey with Herschel (KINGFISH) sample. We model the far-IR/submillimeter spectral energy distribution (SED) using images from Spitzer Space Telescope and Herschel Space Observatory. We calculate the cold dust temperature (T(cold)) and emissivity (beta) on a pixel by pixel basis (where each pixel ranges from 0.1-3 kpc^2) using a two temperature modified blackbody fitting routine. Our fitting method allows us to investigate the resolved nature of temperature and emissivity variations by modeling from the galaxy centers to the outskirts (physical scales of ~15-50 kpc, depending on the size of the galaxy). We fit each SED in two ways: (1) fit T(cold) and beta simultaneously, (2) hold beta constant and fit T(cold). We compare T(cold) and beta with star formation rates (calculated from L(Halpha) and L(24)), the luminosity of the old stellar population (traced through L(3.6), and the dust mass surface density (traced by 500 micron luminosity, L(500)). We find a significant trend between SFR/L(500) and T(cold), implying that the flux of hard UV photons relative to the amount of dust is significantly contributing to the heating of the cold, or diffuse, dust component. We also see a trend between L(3.6)/L(500) and beta, indicating that the old stellar population contributes to the heating at far-IR/submillimeter wavelengths. Finally, we find that when beta is held constant, T(cold) exhibits a strongly decreasing radial trend, illustrating that the shape of the far-IR SED is changing radially through a galaxy, thus confirming on a sample almost double in size the trends observed in Galametz et al. (2012).

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.

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.

We model the infrared to submillimeter spectral energy distribution of 11 nearby galaxies of the KINGFISH sample using Spitzer and Herschel data and compare model extrapolations at 870um (using different fitting techniques) with LABOCA 870um observations. We investigate how the differences between predictions and observations vary with model assumptions or environment. At global scales, we find that modified blackbody models using realistic cold emissivity indices (beta_c=2 or 1.5) are able to reproduce the 870um observed emission within the uncertainties for most of the sample. Low values (beta_c<1.3) would be required in NGC0337, NGC1512 and NGC7793. At local scales, we observe a systematic 870um excess when using beta_=2.0. The beta_c=1.5 or the Draine and Li (2007) models can reconcile predictions with observations in part of the disks. Some of the remaining excesses occur towards the centres and can be partly or fully accounted for by non-dust contributions such as CO(3-2) or, to a lesser extent, free-free or synchrotron emission. In three non-barred galaxies, the remaining excesses rather occur in the disk outskirts. This could be a sign of a flattening of the submm slope (and decrease of the effective emissivity index) with radius in these objects.

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.

Submillimeter excess emission has been reported at 500 microns in a handful of local galaxies, and previous studies suggest that it could be correlated with metal abundance. We investigate the presence of an excess submillimeter emission at 500 microns for a sample of 20 galaxies from the Key Insights on Nearby Galaxies: a Far Infrared Survey with Herschel (KINGFISH) that span a range of morphologies and metallicities (12+log(O/H)=7.8-8.7). We probe the far-infrared (IR) emission using images from the Spitzer Space Telescope and Herschel Space Observatory in the wavelength range 24-500 microns. We model the far-IR peak of the dust emission with a two-temperature modified blackbody and measure excess of the 500 micron photometry relative to that predicted by our model. We compare the submillimeter excess, where present, with global galaxy metallicity and, where available, resolved metallicity measurements. We do not find any correlation between the 500 micron excess and metallicity. A few individual sources do show excess (10-20%) at 500 microns; conversely, for other sources, the model overpredicts the measured 500 micron flux density by as much as 20%, creating a 500 micron "deficit". None of our sources has an excess larger than the calculated 1-sigma uncertainty, leading us to conclude that there is no substantial excess at submillimeter wavelengths at or shorter than 500 microns in our sample. Our results differ from previous studies detecting 500 micron excess in KINGFISH galaxies largely due to new, improved photometry used in this study.

We study the mid- to far-IR properties of a 24um-selected flux-limited sample (S24 > 5mJy) of 154 intermediate redshift (<z>~0.15), infrared luminous galaxies, drawn from the 5MUSES survey. By combining existing mid-IR spectroscopy and new Herschel SPIRE submm photometry from the HerMES program, we derived robust total infrared luminosity (LIR) and dust mass (Md) estimates and infered the relative contribution of the AGN to the infrared energy budget of the sources. We found that the total infrared emission of galaxies with weak 6.2um PAH emission (EW<0.2um) is dominated by AGN activity, while for galaxies with EW>0.2um more than 50% of the LIR arises from star formation. We also found that for galaxies detected in the 250-500um Herschel bands an AGN has a statistically insignificant effect on the temperature of the cold dust and the far-IR colours of the host galaxy, which are primarily shaped by star formation activity. For star-forming galaxies we reveal an anti-correlation between the LIR-to-rest-frame 8um luminosity ratio, IR8 = LIR\L8, and the strength of PAH features. We found that this anti-correlation is primarily driven by variations in the PAHs emission, and not by variations in the 5-15um mid-IR continuum emission. Using the [NeIII]/[NeII] line ratio as a tracer of the hardness of the radiation field, we confirm that galaxies with harder radiation fields tend to exhibit weaker PAH features, and found that they have higher IR8 values and higher dust-mass-weighted luminosities (LIR/Md), the latter being a proxy for the dust temperature (Td). We argue that these trends originate either from variations in the environment of the star-forming regions or are caused by variations in the age of the starburst. Finally, we provide scaling relations that will allow estimating LIR, based on single-band observations with the mid-infrared instrument, on board the upcoming JWST.

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.

We present kiloparsec (kpc) spatial resolution maps of the CO-to-H2 conversion factor (alpha_co) and dust-to-gas ratio (DGR) in 26 nearby, star-forming galaxies. We have simultaneously solved for alpha_co and DGR by assuming that the DGR is approximately constant on kpc scales. With this assumption, we can combine maps of dust mass surface density, CO integrated intensity and HI column density to solve for both alpha_co and DGR with no assumptions about their value or dependence on metallicity or other parameters. Such a study has just become possible with the availability of high resolution far-IR maps from the Herschel key program KINGFISH, 12CO J=(2-1) maps from the IRAM 30m large program HERACLES and HI 21-cm line maps from THINGS. We use a fixed ratio between the (2-1) and (1-0) lines to present our alpha_co results on the more typically used 12CO J=(1-0) scale and show using literature measurements that variations in the line ratio do not effect our results. In total, we derive 782 individual solutions for alpha_co and DGR. On average, alpha_co = 3.1 Msun pc^-2 (K km s^-1)^-1 for our sample with a standard deviation of 0.3 dex. Within galaxies we observe a generally flat profile of alpha_co as a function of galactocentric radius. However, most galaxies exhibit a lower alpha_co in the central kpc---a factor of ~2 below the galaxy mean, on average. In some cases, the central alpha_co value can be factors of 5 to 10 below the standard Milky Way (MW) value of alpha_co,MW =4.4 Msun pc^-2 (K km s^-1)^-1. While for alpha_co we find only weak correlations with metallicity, DGR is well-correlated with metallicity, with an approximately linear slope. Finally, we present several recommendations for choosing an appropriate alpha_co for studies of nearby 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.

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.

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.

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.

We use the near--infrared Br\gamma 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 \mu m emission as a SFR tracer for sub--galactic regions in external galaxies. Br\gamma 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\gamma 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\gamma 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\mu 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\alpha with the 70 \mu m emission, also for use in H{\sc ii} regions. We briefly analyze the PACS 100 and 160 \mu m maps and find that longer wavelengths are not as good SFR indicators as 70\mu 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.

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.

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

We present an extragalactic population model of the cosmic background light to interpret the rich high-quality survey data in the X-ray and IR bands. The model incorporates star-formation and supermassive black hole (SMBH) accretion in a co-evolution scenario to fit simultaneously 617 data points of number counts, redshift distributions and local luminosity functions (LFs) with 19 free parameters. The model has four main components, the total IR LF, the SMBH accretion energy fraction in the IR band, the star-formation SED and the unobscured SMBH SED extinguished with a HI column density distribution. As a result of the observational uncertainties about the star-formation and SMBH SEDs, we present several variants of the model. The best-fit reduced chi^2 reaches as small as 2.7-2.9 of which a significant amount (>0.8) is contributed by cosmic variances or caveats associated with data. Compared to previous models, the unique result of this model is to constrain the SMBH energy fraction in the IR band that is found to increase with the IR luminosity but decrease with redshift up to z ~ 1.5; this result is separately verified using aromatic feature equivalent width data. The joint modelling of X-ray and mid-IR data allows for improved constraints on the obscured AGN, especially the Compton-thick AGN population. All variants of the model require that Compton-thick AGN fractions decrease with the SMBH luminosity but increase with redshift while the type-1 AGN fraction has the reverse trend.

We present measurements of the dust attenuation of H\alpha-selected emission-line galaxies at z=0.8 from the NewH\alpha\ 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] \lambda 3727 to [OIII] \lambda 5007. The spectroscopic sample used in this analysis consists of 341 confirmed H\alpha\ 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\beta\ and H\gamma\ pair of lines is detected with S/N>5 in 55 individual objects and the H\beta\ and H\delta\ 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\alpha\ based on the objects with individually detected lines is A(H\alpha)=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.

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.

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.

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.

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.

We present first results from a narrowband imaging program for intermediate redshift emission-line galaxies using the newly commissioned FourStar infrared camera at the 6.5m Magellan telescope. To enable prompt identification of H\alpha\ emitters, a pair of custom 1% filters, which sample low-airglow atmospheric windows at 1.19 \mu m and 2.10 \mu m, is used to detect both H\alpha\ and [OII]\lambda 3727 emission from the same redshift volume at z=2.2. Initial observations are taken over a 130 arcmin^2 area in the CANDELS-COSMOS field. The exquisite image quality resulting from the combination of the instrument, telescope, and standard site conditions (~0.55" FWHM) allows the 1.19 \mu m and 2.10 \mu m data to probe 3\sigma\ emission-line depths down to 1.0e-17 erg/s/cm^2 and 1.2e-17 erg/s/cm^2 respectively, in less than 10 hours of integration time in each narrowband. For H\alpha\ at z=0.8 and z=2.2, these fluxes correspond to observed star formation rates of ~0.3 and ~4 Msun/yr respectively. We find 122 sources with a 1.19 \mu m excess, and 136 with a 2.10 \mu m excess, 41 of which show an excess in both bands. The dual narrowband technique, as implemented here, is estimated to identify about >80% of z=2.2 H\alpha\ emitters in the narrowband excess population. With the most secure such sample obtained to-date, we compute constraints on the faint-end slope of the z=2.2 H\alpha\ luminosity function. These "narrow-deep" FourStar observations have been obtained as part of the larger NewH\alpha\ Survey, which will combine the data with "wide-shallow" imaging through a similar narrowband filter pair with NEWFIRM at the KPNO/CTIO 4m telescopes, to enable study of both luminous (but rare) and faint emission-line galaxies in the intermediate redshift universe. [Abridged]

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{\mu}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.

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.

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 $\mu$m, [O \one] 63 $\mu$m, and [N \ii] 122 $\mu$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, $\nu f_\nu$(70 $\mu$m)/$\nu f_\nu$(100 $\mu$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\ $\nu f_\nu$(70 $\mu$m)/$\nu f_\nu$(100 $\mu$m) $\lesssim$ 0.95. In the warmest regions, where $\nu f_\nu$(70 $\mu$m)/$\nu f_\nu$(100 $\mu$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 $\mu$m PAH ratio in regions that exhibit a deficit in ([C \ii] + [O \one])/PAH, suggesting that small grains are ionized in these environments.

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.

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.

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.

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.

Using new far-infrared imaging from the Herschel Space Observatory with ancillary data from ultraviolet to submillimeter wavelengths, we estimate the total emission from dust and stars of 62 nearby galaxies in the KINGFISH survey in a way that is as empirical and model-independent as possible. We collect and exploit these data in order to measure from the spectral energy distributions (SEDs) precisely how much stellar radiation is intercepted and re-radiated by dust, and how this quantity varies with galaxy properties. By including SPIRE data, we are more sensitive to emission from cold dust grains than previous analyses at shorter wavelengths, allowing for more accurate estimates of dust temperatures and masses. The dust/stellar flux ratio, which we measure by integrating the SEDs, has a range of nearly three decades. The inclusion of SPIRE data shows that estimates based on data not reaching these far-IR wavelengths are biased low. We find that the dust/stellar flux ratio varies with morphology and total IR luminosity. We also find that dust/stellar flux ratios are related to gas-phase metallicity, while the dust/stellar mass ratios are less so. The substantial scatter between dust/stellar flux and dust/stellar mass indicates that the former is a poor proxy of the latter. Comparing the dust/stellar flux ratios and dust temperatures, we show that early-types tend to have slightly warmer temperatures than spiral galaxies, which may be due to more intense interstellar radiation fields, or to different dust grain compositions. Finally, we show that early-types and early-type spirals have a strong correlation between the dust/stellar flux ratio and specific star formation rate, which suggests that the relatively bright far-IR emission of some of these galaxies is due to ongoing star formation and the radiation field from older stars.

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\alpha\ 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\alpha\ 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\alpha\ 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\alpha\ 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\alpha\ 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.

We have produced the next generation of quasar spectral energy distributions (SEDs), essentially updating the work of Elvis et al. (1994) by using high-quality data obtained with several space and ground-based telescopes, including NASA's Great Observatories. We present an atlas of SEDs of 85 optically bright, non-blazar quasars over the electromagnetic spectrum from radio to X-rays. The heterogeneous sample includes 27 radio-quiet and 58 radio-loud quasars. Most objects have quasi-simultaneous ultraviolet-optical spectroscopic data, supplemented with some far-ultraviolet spectra, and more than half also have Spitzer mid-infrared IRS spectra. The X-ray spectral parameters are collected from the literature where available. The radio, far-infrared, and near-infrared photometric data are also obtained from either the literature or new observations. We construct composite spectral energy distributions for radio-loud and radio-quiet objects and compare these to those of Elvis et al., finding that ours have similar overall shapes, but our improved spectral resolution reveals more detailed features, especially in the mid and near-infrared.

[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$\mu$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$\sigma$) 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}, \phi*=10^{-3.20\pm0.54} Mpc^{-3}, and \alpha=-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.

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\alpha+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\alpha+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.

We present total infrared (IR) and ultraviolet (UV) luminosity functions derived from large representative samples of galaxies at z ~ 0, selected at IR and UV wavelengths from the IRAS IIFSCz catalogue, and the GALEX AIS respectively. We augment these with deep Spitzer and GALEX imaging of galaxies in the 11 Mpc Local Volume Legacy Survey (LVL), allowing us to extend these luminosity functions to lower luminosities (~10^6 L_sun), and providing good constraints on the slope of the luminosity function at the extreme faint end for the first time. Using conventional star formation prescriptions, we generate from our data the SFR distribution function for the local Universe. We find that it has a Schechter form, that the faint-end slope has a constant value (to the limits of our data) of {\alpha} = -1.51 \pm 0.08, and the 'characteristic' SFR is 9.2 M_sun/yr. We also show the distribution function of the SFR volume density; we then use this to calculate a value for the total SFR volume density at z ~ 0 of 0.025 \pm 0.0016 M_sun/yr/Mpc^-3, of which ~ 20% is occurring in starbursts. Decomposing the total star formation by infrared luminosity, it can be seen that 9 \pm 1% is due to LIRGs, and 0.7 \pm 0.2% is occuring in ULIRGs. By comparing UV and IR emission for galaxies in our sample, we also calculate the fraction of star formation occurring in dust obscured environments, and examine the distribution of dusty star formation: we find a very shallow slope at the highly extincted end, which may be attributable to line of sight orientation effects as well as conventional internal extinction.

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.

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.

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.

We use Spitzer 24 $\mu$m, 70 $\mu$m and ground based H$\alpha$ data for a sample of 40 SINGS galaxies to establish a star formation rate (SFR) indicator using 70 $\mu$m emission for sub--galactic ($\sim0.05-2\ \rm{kpc}$) line-emitting regions and to investigate limits in application. A linear correlation between 70 $\mu$m and SFR is found and a star formation indicator SFR(70) is proposed for line-emitting sub-galactic regions as $\rm \Sigma(SFR)\ ({M_{\odot}\cdot yr^{-1}\cdot kpc^{-2}})=9.4\times10^{-44}\ \Sigma(70)\ \rm{(ergs\cdot s^{-1}\cdot kpc^{-2})}$, for regions with $12+\rm{log(O/H)}\gtrsim8.4$ and $\rm \Sigma(SFR)\gtrsim10^{-3}\ (M_{\odot}\cdot yr^{-1}\cdot kpc^{-2})$, with a 1-$\sigma$ dispersion around the calibration of $\sim0.16$ dex. We also discuss the influence of metallicity on the scatter of the data. Comparing with the SFR indicator at 70 $\mu$m for integrated light from galaxies, we find that there is $\sim40%$ excess 70 $\mu$m emission in galaxies, which can be attributed to stellar populations not involved in the current star formation activity.

We study a 24\,$\mu$m selected sample of 330 galaxies observed with the Infrared Spectrograph for the 5\,mJy Unbiased Spitzer Extragalactic Survey. We estimate accurate total infrared luminosities by combining mid-IR spectroscopy and mid-to-far infrared photometry, and by utilizing new empirical spectral templates from {\em Spitzer} data. The infrared luminosities of this sample range mostly from 10$^9$L$_\odot$ to $10^{13.5}$L$_\odot$, with 83% in the range 10$^{10}$L$_\odot$$<$L$_{\rm IR}$$<10^{12}$L$_\odot$. The redshifts range from 0.008 to 4.27, with a median of 0.144. The equivalent widths of the 6.2\,$\mu$m aromatic feature have a bimodal distribution. We use the 6.2\,$\mu$m PAH EW to classify our objects as SB-dominated (44%), SB-AGN composite (22%), and AGN-dominated (34%). The high EW objects (SB-dominated) tend to have steeper mid-IR to far-IR spectral slopes and lower L$_{\rm IR}$ and redshifts. The low EW objects (AGN-dominated) tend to have less steep spectral slopes and higher L$_{\rm IR}$ and redshifts. This dichotomy leads to a gross correlation between EW and slope, which does not hold within either group. AGN dominated sources tend to have lower log(L$_{\rm PAH 7.7\mu m}$/L$_{\rm PAH 11.3\mu m}$) ratios than star-forming galaxies, possibly due to preferential destruction of the smaller aromatics by the AGN. The log(L$_{\rm PAH 7.7\mu m}$/L$_{\rm PAH 11.3\mu m}$) ratios for star-forming galaxies are lower in our sample than the ratios measured from the nuclear spectra of nearby normal galaxies, most probably indicating a difference in the ionization state or grain size distribution between the nuclear regions and the entire galaxy. Finally, we provide a calibration relating the monochromatic 5.8, 8, 14 and 24um continuum or Aromatic Feature luminosity to L$_{\rm IR}$ for different types of objects.

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.

Fitting the spectral energy distributions (SEDs) of galaxies is an almost universally used technique that has matured significantly in the last decade. Model predictions and fitting procedures have improved significantly over this time, attempting to keep up with the vastly increased volume and quality of available data. We review here the field of SED fitting, describing the modelling of ultraviolet to infrared galaxy SEDs, the creation of multiwavelength data sets, and the methods used to fit model SEDs to observed galaxy data sets. We touch upon the achievements and challenges in the major ingredients of SED fitting, with a special emphasis on describing the interplay between the quality of the available data, the quality of the available models, and the best fitting technique to use in order to obtain a realistic measurement as well as realistic uncertainties. We conclude that SED fitting can be used effectively to derive a range of physical properties of galaxies, such as redshift, stellar masses, star formation rates, dust masses, and metallicities, with care taken not to over-interpret the available data. Yet there still exist many issues such as estimating the age of the oldest stars in a galaxy, finer details ofdust properties and dust-star geometry, and the influences of poorly understood, luminous stellar types and phases. The challenge for the coming years will be to improve both the models and the observational data sets to resolve these uncertainties. The present review will be made available on an interactive, moderated web page (sedfitting.org), where the community can access and change the text. The intention is to expand the text and keep it up to date over the coming years.

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.

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.

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

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.

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 $\mu$m, [OIII] 88 $\mu$m, [NII] 122 $\mu$m, [CII] 158 $\mu$m and [NII] 205 $\mu$m. We observe that in the [OI] 63 $\mu$m, [OIII] 88 $\mu$m, and [NII] 122 $\mu$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 $\mu$m/[CII] 158 $\mu$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 $\mu$m and [OIII] 88 $\mu$m maps. Variations of the [NII] 122 $\mu$m/[NII] 205 $\mu$m ratio correspond to a range in the ionized gas density between 150 and 400 cm$^{-3}$.

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.

The total infrared (TIR) luminosity from galaxies can be used to examine both star formation and dust physics. We provide here new relations to estimate the TIR luminosity from various Spitzer bands, in particular from the 8 micron and 24 micron bands. To do so, we use 45" subregions within a subsample of nearby face-on spiral galaxies from the Spitzer Infrared Nearby Galaxies Survey (SINGS) that have known oxygen abundances as well as integrated galaxy data from the SINGS, the Local Volume Legacy Survey (LVL) and Engelbracht et al. (2008) samples. Taking into account the oxygen abundances of the subregions, the star formation rate intensity, and the relative emission of the polycyclic aromatic hydrocarbons at 8 micron, the warm dust at 24 micron and the cold dust at 70 micron and 160 micron we derive new relations to estimate the TIR luminosity from just one or two of the Spitzer bands. We also show that the metallicity and the star formation intensity must be taken into account when estimating the TIR luminosity from two wave bands, especially when data longward of 24 micron are not available.

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

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

We carried out an unbiased, spectroscopic survey using the low-resolution module of the infrared spectrograph (IRS) on board Spitzer targeting two 2.6 square arcminute regions in the GOODS-North field. IRS was used in spectral mapping mode with 5 hours of effective integration time per pixel. One region was covered between 14 and 21 microns and the other between 20 and 35 microns. We extracted spectra for 45 sources. About 84% of the sources have reported detections by GOODS at 24 microns, with a median F_nu(24um) ~ 100 uJy. All but one source are detected in all four IRAC bands, 3.6 to 8 microns. We use a new cross-correlation technique to measure redshifts and estimate IRS spectral types; this was successful for ~60% of the spectra. Fourteen sources show significant PAH emission, four mostly SiO absorption, eight present mixed spectral signatures (low PAH and/or SiO) and two show a single line in emission. For the remaining 17, no spectral features were detected. Redshifts range from z ~ 0.2 to z ~ 2.2, with a median of 1. IR Luminosities are roughly estimated from 24 microns flux densities, and have median values of 2.2 x 10^{11} L_{\odot} and 7.5 x 10^{11} L_{\odot} at z ~ 1 and z ~ 2 respectively. This sample has fewer AGN than previous faint samples observed with IRS, which we attribute to the fainter luminosities reached here.

Compared to starburst galaxies, normal star forming galaxies have been shown to display a much larger dispersion of the dust attenuation at fixed reddening through studies of the IRX-beta diagram (the IR/UV ratio "IRX" versus the UV color "beta"). To investigate the causes of this larger dispersion and attempt to isolate second parameters, we have used GALEX UV, ground-based optical, and Spitzer infrared imaging of 8 nearby galaxies, and examined the properties of individual UV and 24 micron selected star forming regions. We concentrated on star-forming regions, in order to isolate simpler star formation histories than those that characterize whole galaxies. We find that 1) the dispersion is not correlated with the mean age of the stellar populations, 2) a range of dust geometries and dust extinction curves are the most likely causes for the observed dispersion in the IRX-beta diagram 3) together with some potential dilution of the most recent star-forming population by older unrelated bursts, at least in the case of star-forming regions within galaxies, 4) we also recover some general characteristics of the regions, including a tight positive correlation between the amount of dust attenuation and the metal content. Although generalizing our results to whole galaxies may not be immediate, the possibility of a range of dust extinction laws and geometries should be accounted for in the latter systems as well.

(abridged) Using a complete sample of ~300 star-forming galaxies within 11 Mpc, we evaluate the consistency between star formation rates (SFRs) inferred from the far ultraviolet (FUV) non-ionizing continuum and H-alpha nebular emission, assuming standard conversion recipes in which the SFR scales linearly with luminosity at a given wavelength. Our analysis probes SFRs over 5 orders of magnitude, down to ultra-low activities on the order of ~0.0001 M_sun/yr. The data are drawn from the 11 Mpc H-alpha and Ultraviolet Galaxy Survey (11HUGS), which has obtained H-alpha fluxes from ground-based narrowband imaging, and UV fluxes from imaging with GALEX. For normal spiral galaxies (SFR~1 M_sun/yr), our results are consistent with previous work which has shown that FUV SFRs tend to be lower than H-alpha SFRs before accounting for internal dust attenuation, but that there is relative consistency between the two tracers after proper corrections are applied. However, a puzzle is encountered at the faint end of the luminosity function. As lower luminosity dwarf galaxies, roughly less active than the Small Magellanic Cloud, are examined, H-alpha tends to increasingly under-predict the SFR relative to the FUV. Although past studies have suggested similar trends, this is the first time this effect is probed with a statistical sample for galaxies with SFR~<0.1 M_sun/yr. A range of standard explanations does not appear to be able to account for the magnitude of the systematic. Some recent work has argued for an IMF which is deficient in high mass stars in dwarf and low surface brightness galaxies, and we also consider this scenario.

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

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.

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.

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]$\lambda$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.

When, and how, did the first galaxies and supermassive black holes (SMBH) form, and how did they reionization the Universe? First galaxy formation and cosmic reionization are among the last frontiers in studies of cosmic structure formation. We delineate the detailed astrophysical probes of early galaxy and SMBH formation afforded by observations at centimeter through submillimeter wavelengths. These observations include studies of the molecular gas (= the fuel for star formation in galaxies), atomic fine structure lines (= the dominant ISM gas coolant), thermal dust continuum emission (= an ideal star formation rate estimator), and radio continuum emission from star formation and relativistic jets. High resolution spectroscopic imaging can be used to study galaxy dynamics and star formation on sub-kpc scales. These cm and mm observations are the necessary compliment to near-IR observations, which probe the stars and ionized gas, and X-ray observations, which reveal the AGN. Together, a suite of revolutionary observatories planned for the next decade from centimeter to X-ray wavelengths will provide the requisite panchromatic view of the complex processes involved in the formation of the first generation of galaxies and SMBHs, and cosmic reionization.

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.

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.

Far-infrared line and continuum fluxes are presented for a sample of 227 galaxies observed with the Long Wavelength Spectrometer on the Infrared Space Observatory. The galaxy sample includes normal star-forming systems, starbursts, and active galactic nuclei covering a wide range of colors and morphologies. The dataset spans some 1300 line fluxes, 600 line upper limits, and 800 continuum fluxes. Several fine structure emission lines are detected that arise in either photodissociation or HII regions: [OIII]52um, [NIII]57um, [OI]63um, [OIII]88um, [NII]122um, [OI]145um, and [CII]158um. Molecular lines such as OH at 53um, 79um, 84um, 119um, and 163um, and H2O at 58um, 66um, 75um, 101um, and 108um are also detected in some galaxies. In addition to those lines emitted by the target galaxies, serendipitous detections of Milky Way [CII]158um and an unidentified line near 74um in NGC1068 are also reported. Finally, continuum fluxes at 52um, 57um, 63um, 88um, 122um, 145um, 158um, and 170um are derived for a subset of galaxies in which the far-infrared emission is contained within the ~75" ISO LWS beam. The statistics of this large database of continuum and line fluxes, including trends in line ratios with the far-infrared color and infrared-to-optical ratio, are explored.

The Wyoming Survey for H-alpha, or WySH, is a large-area, ground-based, narrowband imaging survey for H-alpha-emitting galaxies over the latter half of the age of the Universe. The survey spans several square degrees in a set of fields of low Galactic cirrus emission. The observing program focuses on multiple dz~0.02 epochs from z~0.16 to z~0.81 down to a uniform (continuum+line) luminosity at each epoch of ~10^33 W uncorrected for extinction (3sigma for a 3" diameter aperture). First results are presented here for 98+208 galaxies observed over approximately 2 square degrees at redshifts z~0.16 and 0.24, including preliminary luminosity functions at these two epochs. These data clearly show an evolution with lookback time in the volume-averaged cosmic star formation rate. Integrals of Schechter fits to the extinction-corrected H-alpha luminosity functions indicate star formation rates per co-moving volume of 0.009 and 0.014 h_70 M_sun/yr/Mpc^3 at z~0.16 and 0.24, respectively. The formal uncertainties in the Schechter fits, based on this initial subset of the survey, correspond to uncertainties in the cosmic star formation rate density at the >~40% level; the tentative uncertainty due to cosmic variance is 25%, estimated from separately carrying out the analysis on data from the first two fields with substantial datasets.

We have mapped the warm molecular gas traced by the H_2 S(0) - H_2 S(5) pure rotational mid-infrared emission lines over a radial strip across the nucleus and disk of M51 (NGC 5194) using the Infrared Spectrograph (IRS) on the Spitzer Space Telescope. The six H_2 lines have markedly different emission distributions. We obtained the H_2 temperature and surface density distributions by assuming a two temperature model: a warm (T = 100 - 300 K) phase traced by the low J (S(0) - S(2)) lines and a hot phase (T = 400 - 1000 K) traced by the high J (S(2) - S(5)) lines. The lowest molecular gas temperatures are found within the spiral arms (T ~ 155 K), while the highest temperatures are found in the inter-arm regions (T > 700 K). The warm gas surface density reaches a maximum of 11 M_sun/pc^2 in the northwestern spiral arm, whereas the hot gas surface density peaks at 0.24 M_sun/pc^2 at the nucleus. The spatial offset between the peaks in the warm and hot phases and the differences in the distributions of the H_2 line emission suggest that the warm phase is mostly produced by UV photons in star forming regions while the hot phase is mostly produced by shocks or X-rays associated with nuclear activity. The warm H_2 is found in the dust lanes of M51, spatially offset from the brightest HII regions. The warm H_2 is generally spatially coincident with the cold molecular gas traced by CO (J = 1 - 0) emission, consistent with excitation of the warm phase in dense photodissociation regions (PDRs). In contrast, the hot H_2 is most prominent in the nuclear region. Here, over a 0.5 kpc radius around the nucleus of M51, the hot H_2 coincides with [O IV](25.89 micron) and X-ray emission indicating that shocks and/or X-rays are responsible for exciting this phase.

Three dimensional spectroscopy of extended sources is typically performed with dedicated integral field spectrographs. We describe a method of reconstructing full spectral cubes, with two spatial and one spectral dimension, from rastered spectral mapping observations employing a single slit in a traditional slit spectrograph. When the background and image characteristics are stable, as is often achieved in space, the use of traditional long slits for integral field spectroscopy can substantially reduce instrument complexity over dedicated integral field designs, without loss of mapping efficiency -- particularly compelling when a long slit mode for single unresolved source followup is separately required. We detail a custom flux-conserving cube reconstruction algorithm, discuss issues of extended source flux calibration, and describe CUBISM, a tool which implements these methods for spectral maps obtained with ther Spitzer Space Telescope's Infrared Spectrograph.

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.

We construct a large data set of global structural parameters for 1300 field and cluster spiral galaxies and explore the joint distribution of luminosity L, optical rotation velocity V, and disk size R at I- and 2MASS K-bands. The I- and K-band velocity-luminosity (VL) relations have log-slopes of 0.29 and 0.27, respectively with sigma_ln(VL)~0.13, and show a small dependence on color and morphological type in the sense that redder, early-type disk galaxies rotate faster than bluer, later-type disk galaxies for most luminosities. The VL relation at I- and K-bands is independent of surface brightness, size and light concentration. The log-slope of the I- and K-band RL relations is a strong function of morphology and varies from 0.25 to 0.5. The average dispersion sigma_ln(RL) decreases from 0.33 at I-band to 0.29 at K, likely due to the 2MASS selection bias against lower surface brightness galaxies. Measurement uncertainties are sigma_ln(V)~0.09, sigma_ln(L)~0.14 and somewhat larger and harder to estimate for ln(R). The color dependence of the VL relation is consistent with expectations from stellar population synthesis models. The VL and RL residuals are largely uncorrelated with each other; the RV-RL residuals show only a weak positive correlation. These correlations suggest that scatter in luminosity is not a significant source of the scatter in the VL and RL relations. The observed scaling relations can be understood in the context of a model of disk galaxies embedded in dark matter halos that invokes low mean spin parameters and dark halo expansion, as we describe in our companion paper (Dutton et al. 2007). We discuss in two appendices various pitfalls of standard analytical derivations of galaxy scaling relations, including the Tully-Fisher relation with different slopes. (Abridged).

Dust in galaxies can be mapped by either the FIR/sub-mm emission, the optical or infrared reddening of starlight, or the extinction of a known background source. We compare two dust extinction measurements for a set of fifteen sections in thirteen nearby galaxies, to determine the scale of the dusty ISM responsible for disk opacity: one using stellar reddening and the other a known background source. In our earlier papers, we presented extinction measurements of 29 galaxies, based on calibrated counts of distant background objects identified though foreground disks in HST/WFPC2 images. For the 13 galaxies that overlap with the Spitzer Infrared Nearby Galaxies Survey (SINGS), we now compare these results with those obtained from an I-L color map. Our goal is to determine whether or not a detected distant galaxy indicates a gap in the dusty ISM, and hence to better understand the nature and geometry of the disk extinction. We find that distant galaxies are predominantly in low-extinction sections marked by the color maps, indicating that their number depends both on the cloud cover of {\it Spitzer}-resolved dust structures --mostly the spiral arms--and a diffuse, unresolved underlying disk. We note that our infrared color map (E[I-L]) underestimates the overall dust presence in these disks severely, because it implicitly assumes the presence of a dust screen in front of the stellar distribution.

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.

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.

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

Using the new capabilities of the Spitzer Space Telescope and extensive multiwavelength data from the Spitzer Infrared Nearby Galaxies Survey (SINGS), it is now possible to study the infrared properties of star formation in nearby galaxies down to scales equivalent to large HII regions. We are therefore able to determine what fraction of large, infrared-selected star-forming regions in normal galaxies are highly obscured and address how much of the star formation we miss by relying solely on the optical portion of the spectrum. Employing a new empirical method for deriving attenuations of infrared-selected star-forming regions we investigate the statistics of obscured star formation on 500pc scales in a sample of 38 nearby galaxies. We find that the median attenuation is 1.4 magnitudes in H-alpha and that there is no evidence for a substantial sub-population of uniformly highly-obscured star-forming regions. The regions in the highly-obscured tail of the attenuation distribution (A_H-alpha > 3) make up only ~4% of the sample of nearly 1800 regions, though very embedded infrared sources on the much smaller scales and lower luminosities of compact and ultracompact HII regions are almost certainly present in greater numbers. The highly-obscured cases in our sample are generally the bright, central regions of galaxies with high overall attenuation but are not otherwise remarkable. We also find that a majority of the galaxies show decreasing radial trends in H-alpha attenuation. The small fraction of highly-obscured regions seen in this sample of normal, star-forming galaxies suggests that on 500pc scales the timescale for significant dispersal or break up of nearby, optically-thick dust clouds is short relative to the lifetime of a typical star-forming region.

Physical dust models are presented for 65 galaxies in the SINGS survey that are strongly detected in the four IRAC bands and three MIPS bands. For each galaxy we estimate (1) the total dust mass, (2) the fraction of the dust mass contributed by PAHs, and (3) the intensity of the starlight heating the dust grains. We find that spiral galaxies have dust properties resembling the dust in the local region of the Milky Way, with similar dust-to-gas ratio, and similar PAH abundance. The observed SEDs, including galaxies with SCUBA photometry, can be reproduced by dust models that do not require "cold" (T<10K) dust. The dust-to-gas ratio is observed to be dependent on metallicity. In the interstellar media of galaxies with A_O=12+log(O/H)>8.1, grains contain a substantial fraction of interstellar Mg, Si and Fe. Galaxies with A_O<8.1 and extended HI envelopes in some cases appear to have global dust-to-gas ratios that are low for their measured oxygen abundance, but in the regions where infrared emission is detected, the dust-to-gas ratio generally appears to be consistent with a substantial fraction of interstellar Mg, Si, and Fe being contained in dust. The PAH index q_PAH -- the fraction of the dust mass in the form of PAHs -- correlates with metallicity. The nine galaxies in our sample with A_O<8.1 have a median q_PAH=1.0%, whereas galaxies with A_O>8.1 have a median q_PAH=3.55%. The derived dust masses favor a value X_CO approx 4e20 cm^{-2}(K kms)^{-1} for the CO to H_2 conversion factor. Except for some starbursting systems (Mrk33, Tolo89, NGC3049), dust in the diffuse ISM dominates the IR power.

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

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.

The ultraviolet-to-radio continuum spectral energy distributions are presented for all 75 galaxies in the Spitzer Infrared Nearby Galaxies Survey (SINGS). A principal component analysis of the sample shows that most of the sample's spectral variations stem from two underlying components, one representative of a galaxy with a low infrared-to-ultraviolet ratio and one representative of a galaxy with a high infrared-to-ultraviolet ratio. The influence of several parameters on the infrared-to-ultraviolet ratio is studied (e.g., optical morphology, disk inclination, far-infrared color, ultraviolet spectral slope, and star formation history). Consistent with our understanding of normal star-forming galaxies, the SINGS sample of galaxies in comparison to more actively star-forming galaxies exhibits a larger dispersion in the infrared-to-ultraviolet versus ultraviolet spectral slope correlation. Early type galaxies, exhibiting low star formation rates and high optical surface brightnesses, have the most discrepant infrared-to-ultraviolet correlation. These results suggest that the star formation history may be the dominant regulator of the broadband spectral variations between galaxies. Finally, a new discovery shows that the 24 micron morphology can be a useful tool for parametrizing the global dust temperature and ultraviolet extinction in nearby galaxies. The dust emission in dwarf/irregular galaxies is clumpy and warm accompanied by low ultraviolet extinction, while in spiral galaxies there is typically a much larger diffuse component of cooler dust and average ultraviolet extinction. For galaxies with nuclear 24 micron emission, the dust temperature and ultraviolet extinction are relatively high compared to disk 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.

We present Spitzer imaging of the metal-deficient (Z ~30% Z_sun) Local Group dwarf galaxy NGC 6822. On spatial scales of ~130 pc, we study the nature of IR, H alpha, HI, and radio continuum emission. Nebular emission strength correlates with IR surface brightness; however, roughly half of the IR emission is associated with diffuse regions not luminous at H alpha (as found in previous studies). The global ratio of dust to HI gas in the ISM, while uncertain at the factor of ~2 level, is ~25 times lower than the global values derived for spiral galaxies using similar modeling techniques; localized ratios of dust to HI gas are about a factor of five higher than the global value in NGC 6822. There are strong variations (factors of ~10) in the relative ratios of H alpha and IR flux throughout the central disk; the low dust content of NGC 6822 is likely responsible for the different H alpha/IR ratios compared to those found in more metal-rich environments. The H alpha and IR emission is associated with high-column density (> ~1E21 cm^-2) neutral gas. Increases in IR surface brightness appear to be affected by both increased radiation field strength and increased local gas density. Individual regions and the galaxy as a whole fall within the observed scatter of recent high-resolution studies of the radio-far IR correlation in nearby spiral galaxies; this is likely the result of depleted radio and far-IR emission strengths in the ISM of this dwarf galaxy.

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.

The recent star formation (SF) in the early-type spiral galaxy M81 is characterized using imaging observations from the far-ultraviolet (UV) to the far-infrared (IR). We compare these data with models of the stellar, gas, and dust emission for sub-galactic regions. Our results suggest the existence of a diffuse dust emission not directly linked to the recent SF. We find a radial decrease of the dust temperature and dust mass density, and in the attenuation of the stellar light. The IR emission in M81 can be modeled with three components: 1) cold dust with a temperature <T_c>=18+-2 K, concentrated near the HII regions but also presenting a diffuse distribution; 2) warm dust with T_w=53+-7 K, directly linked with the HII regions; and 3) aromatic molecules, with diffuse morphology peaking around the HII regions. We derive several relationships to obtain total IR luminosities from IR monochromatic fluxes, and we compare five different star formation rate (SFR) estimators for HII regions in M81 and M51: the UV, Halpha, and three estimators based on Spitzer data. We find that the Halpha luminosity absorbed by dust correlates tightly with the 24 microns emission. The correlation with the total IR luminosity is not as good. Important variations from galaxy to galaxy are found when estimating the total SFR with the 24 microns or the total IR emission alone. The most reliable estimations of the total SFRs are obtained by combining the Halpha emission (or the UV) and an IR luminosity (especially the 24 microns emission), which probe the unobscured and obscured SF, respectively. For the entire M81 galaxy, about 50% of the total SF is obscured by dust. The percentage of obscured SF ranges from 60% in the inner regions of the galaxy to 30% in the outer zones.

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)

We analyze extensive data on NGC1377 from the Spitzer Infrared Nearby Galaxies Survey (SINGS). Within the category of nascent starbursts, that we previously selected by their infrared to radio continuum ratios in large excess of the average and their hot dust, NGC1377 has the largest infrared excess yet measured. Optical imaging reveals a morphological distortion suggestive of a recent accretion event. Infrared spectroscopy reveals a compact and opaque source dominated by a hot, self-absorbed continuum (tau ~ 20 in the 10 micron silicate band). We provide physical evidence against non-stellar activity being the heating source. HII regions are detected through the single [NeII] line, probing <1% of the ionizing radiation. Not only is the optical depth very high, but >85% of ionizing photons are suppressed by dust. The only other detected emission features are molecular hydrogen lines, arguably excited mainly by shocks, besides photodissociation regions, and weak aromatic bands. The new observations support our interpretation in terms of an extremely young starburst (<1 Myr). More generally, galaxies deficient in radio synchrotron are likely observed within a few Myr of the onset of a starburst and after a long quiescence, prior to the replenishment of the ISM with cosmic rays. The similar infrared-radio properties of NGC1377 and some infrared-luminous galaxies suggest that NGC1377 constitutes an archetype to better understand starburst evolution. Although rare locally because observed in a brief evolutionary stage, nascent starbursts may represent a non-negligible fraction of merger-induced starbursts dominating deep infrared counts. Since they differ dramatically from usual starburst templates, they bear important consequences for the interpretation of deep surveys.

The Spitzer Space Telescope is pushing into new frontiers in high redshift astronomy. Closer to home, Spitzer is making an equally large impact on our understanding of galaxy formation and evolution. In this contribution we present mid-infrared diagnostics based largely on data from the Spitzer Infrared Nearby Galaxies Survey (SINGS). Our main result is that these mid-infrared diagnostics effectively constrain a target's dominant power source. The combination of a high ionization line index and PAH strength serves as an efficient discriminant between AGN and star-forming nuclei, confirming progress made with ISO spectroscopy on starbursting and ultraluminous infrared galaxies. The sensitivity of Spitzer allows us to probe fainter nuclei and star-forming regions within galaxy disks. We find that both star-forming nuclei and extranuclear regions stand apart from nuclei that are powered by Seyfert or LINER activity. In fact, we identify areas within three diagnostic diagrams containing >90% Seyfert/LINER nuclei or >90% HII regions/HII nuclei. We also find that, compared to starbursting nuclei, extranuclear regions typically separate even further from AGN, especially for low-metallicity extranuclear environments. In addition, instead of the traditional mid-infrared approach to differentiating between AGN and star-forming sources that utilizes relatively weak high-ionization lines, we show that strong low-ionization cooling lines of X-ray dominated regions like [SiII] 34.82 micron can alternatively be used as excellent discrimants.

Mid-infrared diagnostics are presented for a large portion of the Spitzer Infrared Nearby Galaxies Survey (SINGS) sample plus archival data from the Infrared Space Observatory and the Spitzer Space Telescope. Our main result is that these mid-infrared diagnostics effectively constrain a target's dominant power source. The combination of a high ionization line index and PAH strength serves as an efficient discriminant between AGN and star-forming nuclei, confirming progress made with ISO spectroscopy on starbursting and ultraluminous infrared galaxies. The sensitivity of Spitzer allows us to probe fainter nuclei and star-forming regions within galaxy disks. We find that both star-forming nuclei and extranuclear regions stand apart from nuclei that are powered by Seyfert or LINER activity. In fact, we identify areas within four diagnostic diagrams containing >90% Seyfert/LINER nuclei or >90% HII regions/HII nuclei. We also find that, compared to starbursting nuclei, extranuclear regions typically separate even further from AGN, especially for low-metallicity extranuclear environments. In addition, instead of the traditional mid-infrared approach to differentiating between AGN and star-forming sources that utilizes relatively weak high-ionization lines, we show that strong low-ionization cooling lines of X-ray dominated regions like [SiII] 34.82 micron can alternatively be used as excellent discrimants. Finally, the typical target in this sample shows relatively modest interstellar electron density (~400 cm^{-3}) and obscuration (A_V ~ 1.0 mag for a foreground screen), consistent with a lack of dense clumps of highly obscured gas and dust residing in the emitting regions. [Abridged]

We present new images (groundbased optical and mid-infrared (MIR) from the Spitzer Space Telescope) and spectra (from Spitzer) of the archetypal starburst galaxy M 82. The Spitzer data show that the MIR emission extends at least 6 kpc along the minor axis of the galaxy. We use the optical and infrared data to demonstrate that the extended emission is dominated by emission from dust. The colors of the MIR emission and the spectra indicate that there is a strong component of aromatic feature emission (the MIR features commonly attributed to polycyclic aromatic hydrocarbons). The dust continuum and aromatic feature emission are both strong in the well-known superwind region of this galaxy; clearly the carrier of the aromatic features can survive in close proximity to the wind, far from the plane of the galaxy. We also see significant emission by dust well outside the superwind region, providing the clearest picture to date of the dust distribution in the halo of this galaxy.

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.

We look for evidence of tidal stripping in elliptical galaxies through the analysis of homogeneous CCD data corresponding to a sample of 228 elliptical galaxies belonging to 24 clusters of galaxies at $0.015<z<0.080$. We investigate departures from the standard magnitude-isophotal size relation, as a function of environmental (cluster-centric distance, local galaxy density) and structural (cluster velocity dispersion, Bautz-Morgan type) properties. We find that, for any particular galaxy luminosity, the ellipticals in the inner and denser regions of the clusters are about 5% smaller than those in the outer regions, which is in good agreement with the finding of Strom & Strom (1978) based on photographic photometry. The null hypothesis (ie., galaxy sizes are independent of the cluster-centric distance or density) is rejected at a significance level of better than 99.7%. Numericals models of Aguilar & White (1986) predict that tidal stripping can lead to changes in the whole structure of ellipticals producing shrinkage and brightening of the galaxy, qualitatively consistent with our measurements and also with the findings of Trujillo et al. (2002), that more centrally concentrated ellipticals populate denser regions. Our observational results can be interpreted as evidence for stripping of stars from ellipticals in the central/denser regions of clusters, contributing to the intra-cluster light observed in these structures.

(Abridged) We present an initial look at the far infrared-radio correlation within the star-forming disks of four nearby, nearly face-on galaxies (NGC 2403, NGC 3031, NGC 5194, and NGC 6946). Using Spitzer MIPS imaging and WSRT radio continuum data, observed as part of the Spitzer Infrared Nearby Galaxies Survey (SINGS), we are able to probe variations in the logarithmic 24mu/22cm (q_24) and 70mu/22cm (q_70) surface brightness ratios across each disk at sub-kpc scales. We find general trends of decreasing q_24 and q_70 with declining surface brightness and with increasing radius. The residual dispersion around the trend of q_24 and q_70 versus surface brightness is smaller than the residual dispersion around the trend of q_24 and q_70 versus radius, on average by ~0.1 dex, indicating that the distribution of star formation sites is more important in determining the infrared/radio disk appearance than the exponential profiles of disks. We have also performed preliminary phenomenological modeling of cosmic ray electron (CRe^-) diffusion using an image-smearing technique, and find that smoothing the infrared maps improves their correlation with the radio maps. Exponential kernels tend to work better than Gaussian kernels which suggests that additional processes besides simple random-walk diffusion in three dimensions must affect the evolution of CRe^-s. The best fit smoothing kernels for the two less active star-forming galaxies (NGC 2403 and NGC 3031) have much larger scale-lengths than those of the more active star-forming galaxies (NGC 5194 and NGC 6946). This difference may be due to the relative deficit of recent CRe^- injection into the interstellar medium (ISM) for the galaxies having largely quiescent disks.

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.

The Spitzer Infrared Nearby Galaxies Survey (SINGS) is carrying out a comprehensive multi-wavelength survey on a sample of 75 nearby galaxies. The 1-850um spectral energy distributions are presented using broadband imaging data from Spitzer, 2MASS, ISO, IRAS, and SCUBA. The infrared colors derived from the globally-integrated Spitzer data are generally consistent with the previous generation of models that were developed based on global data for normal star-forming galaxies, though significant deviations are observed. Spitzer's excellent sensitivity and resolution also allow a detailed investigation of the infrared spectral energy distributions for various locations within the three large, nearby galaxies NGC3031 (M81), NGC5194 (M51), and NGC7331. Strong correlations exist between the local star formation rate and the infrared colors f_nu(70um)/f_nu(160um) and f_nu(24um)/f_nu(160um), suggesting that the 24 and 70um emission are useful tracers of the local star formation activity level. Preliminary evidence indicates that variations in the 24um emission, and not variations in the emission from polycyclic aromatic hydrocarbons at 8um, drive the variations in the f_nu(8.0um)/f_nu(24um) colors within NGC3031, NGC5194, and NGC7331. If the galaxy-to-galaxy variations in spectral energy distributions seen in our sample are representative of the range present at high redshift then extrapolations of total infrared luminosities and star formation rates from the observed 24um flux will be uncertain at the factor-of-five level (total range). The corresponding uncertainties using the redshifted 8.0um flux (e.g. observed 24um flux for a z=2 source) are factors of 10-20. Considerable caution should be used when interpreting such extrapolated infrared luminosities.

(Abridged) Far ultraviolet to far infrared images of the nearby galaxy NGC5194, from Spitzer, GALEX, Hubble Space Telescope and ground--based data, are used to investigate local and global star formation, and the impact of dust extinction in HII-emitting knots. In the IR/UV-UV color plane, the NGC5194 HII knots show the same trend observed for normal star-forming galaxies, having a much larger dispersion than starburst galaxies. We identify the dispersion as due to the UV emission predominantly tracing the evolved, non-ionizing stellar population, up to ages 50-100 Myr. While in starbursts the UV light traces the current SFR, in NGC5194 it traces a combination of current and recent-past SFR. Unlike the UV emission, the monochromatic 24 micron luminosity is an accurate local SFR tracer for the HII knots in NGC5194; this suggests that the 24 micron emission carriers are mainly heated by the young, ionizing stars. However, preliminary results show that the ratio of the 24 micron emission to the SFR varies by a factor of a few from galaxy to galaxy. While also correlated with star formation, the 8 micron emission is not directly proportional to the number of ionizing photons. This confirms earlier suggestions that the carriers of the 8 micron emission are heated by more than one mechanism.

We examine colors from 3.6 micron to 24 micron as a function of metallicity (O/H) for a sample of 34 galaxies. The galaxies range over 2 orders of magnitude in metallicity. They display an abrupt shift in the 8 micron to 24 micron color between metallicities 1/3 to 1/5 of the solar value. The mean 8 micron to 24 micron flux density ratio below and above 12 + log (O/H) = 8.2 is 0.08 +/- 0.04 and 0.70 +/- 0.53, respectively. We use mid-infrared colors and spectroscopy to demonstrate that the shift is primarily due to a decrease in the 8 micron flux density as opposed to an increase in the 24 micron flux density. This result is most simply interpreted as due to a weakening at low metallicity of the mid-infrared emission bands usually attributed to PAHs (polycyclic aromatic hydrocarbons) relative to the small-grain dust emission. However, existing empirical spectral energy distribution models cannot account for the observed short-wavelength (i.e., below 8 micron) colors of the low-metallicity galaxies merely by reducing the strength of the PAH features; some other emission source (e.g., hot dust) is required.

New and archival interferometric 12CO(1->0) datasets from six nearby galaxies are combined with H_2 2.122um and H-alpha maps to explore in detail the interstellar medium in different star-forming galaxies. We investigate the relation between warm (H_2 at T~2000 K) and cold (CO at T~50 K) molecular gas from 100 pc to 2 kpc scales. On these scales, the ratio of warm-to-cold molecular hydrogen correlates with the fnu(60um)/fnu(100um) ratio, a ratio that tracks the star formation activity level. This result also holds for the global properties of galaxies from a much larger sample drawn from the literature. The trend persists for over three orders of magnitude in the mass ratio, regardless of source nuclear activity.

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.

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.

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.

We present the most recent results of an optical survey of the position error contours ("error boxes") of unidentified high energy gamma-ray sources.

The Infrared Space Observatory was used to search for a tracer of the warm and dense neutral interstellar medium, the [O I] 63.18 micron line, in four ultraluminous IRAS sources lying at redshifts between 0.6 and 1.4. While these sources are quasars, their infrared continuum emission suggests a substantial interstellar medium. No [O I] flux was securely detected after probing down to a 3 sigma sensitivity level sufficient for detecting line emission in starbursts with similar continuum emission. However, if the detection threshold is slightly relaxed, one target is detected with 2.7 sigma significance. For this radio-quiet quasar there is likely a substantial dense and warm interstellar medium; the upper limits for the three radio-loud sources do not preclude the same conclusion. Using a new, uniformly-processed database of the ISO extragalactic far-infrared spectroscopy observations, it is shown that nearby Seyfert galaxies typically have higher [O I]-to-far-infrared ratios than do normal star-forming galaxies, so the lack of strong [O I] 63 micron emission from these high-redshift ultraluminous sources cannot be attributed to their active cores.

(Abridged) This paper presents an exploration of two fundamental scaling relations of spiral galaxies, the luminosity-rotation speed (or Tully-Fisher; TF) relation, and the size-luminosity (SL) relation, and the dependences of their scatter, at red and IR bands. We verify that the observed virial relations of disk galaxies are given by Vobs \prop L^0.31 and Rd \prop L^0.33 using distance-redshift surveys of high surface brightness (HSB) and low surface brightness (LSB) non-interacting galaxies. Thus the galaxy surface brightness (SB) is Sigma \prop Rd \prop Vobs \prop Mvir^0.33. We derive the first extensive J-band TF and SL relations using 2MASS data products. The TF residuals are fully independent of SB and other tested observables. TF relations (TFRs) for various morphological types have different zero points, but a common slope, such that early-type (redder) disks rotate faster than later-types (bluer) of the same luminosity. This dependence of the TFR is a direct consequence of the more fundamental dependence on color, which itself is related to the galaxy star formation history. The SL relation scatter is mostly dominated by SB, but color also plays a small role. The TF and SL residuals for HSB and LSB galaxies are weakly correlated, in agreement with the earlier claim by Courteau & Rix (1999). This suggests that spiral disks of all color, brightness, and barredness are, on average, dominated by dark matter even in their inner parts. The observed scaling relations of disk galaxies are likely the result of the simplest scenario of galaxy formation, but their color dependence remains a challenge to hierarchical structure formation models. The scaling relations and their derivatives are contrasted against galaxy formation models in our companion paper.

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.

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.

We have obtained new optical imaging and spectroscopic observations of 78 galaxies in the fields of the rich clusters Abell 1413 (z = 0.14), Abell 2218 (z = 0.18) and Abell 2670 (z = 0.08). We have detected line emission from 25 cluster galaxies plus an additional six galaxies in the foreground and background, a much lower success rate than what was found (65%) for a sample of 52 lower-richness Abell clusters in the range 0.02 < z < 0.08. We have combined these data with our previous observations of Abell 2029 and Abell 2295 (both at z = 0.08), which yields a sample of 156 galaxies. We evaluate several parameters as a function of cluster environment: Tully-Fisher residuals, H-alpha equivalent width, and rotation curve asymmetry, shape and extent. Although H-alpha is more easily detectable in galaxies that are located further from the cluster cores, we fail to detect a correlation between H-alpha extent and galaxy location in those where it is detected, again in contrast with what is found in the clusters of lesser richness. We fail to detect any statistically significant trends for the other parameters in this study. The zero-point in the z~0.1 Tully-Fisher relation is marginally fainter (by 1.5 sigma) than that found in nearby clusters, but the scatter is essentially unchanged.

We present a characterization of the local luminosity-color, bi-variate distribution of IRAS galaxies from the 1.2Jy sample, selected at 60 micron. The S60/S100 infrared color is used as the best single parameter description of the IR spectral energy distribution of galaxies. We derive an analytical form of the distribution and use it to constrain the effect of the IR color distribution on evolution models for high redshift, farIR luminous galaxies. Our adopted evolution retains the locally observed correlation between luminosity and color, such that the larger characteristic luminosities at higher redshift have a warmer characteristic color. The width of the color distribution at a given luminosity remains constant for all redshifts. We demonstrate that there is the potential for both hotter and colder sources to be missed in cosmological surveys. An evolving bi-variate luminosity function coupled with the cold source bias of sub-mm selected surveys suggests the existence of a large population of cold sources appearing in such surveys. Likewise, a hot source bias for most SIRTF wavebands together with a bi-variate model suggests an excess of hot sources being selected. We test the evolutionary form against available data for higher redshift, farIR galaxies. The data do not reveal evidence for any strong evolution in the characteristic luminosity-color distribution as a function of redshift over 0<z<1. However, there is marginal evidence for a broadening of the color distribution at higher redshifts, consistent with our locally characterized trend of a broadening in the IR color distribution at the highest luminosities.

We present ISO-PHOT spectra of the regions 2.5-4.9um and 5.8-11.6um for a sample of 45 disk galaxies from the U.S. ISO Key Project on Normal Galaxies. The spectra can be decomposed into three spectral components: (1) continuum emission from stellar photospheres, which dominates the near-infrared (2.5- 4.9um; NIR) spectral region; (2) a weak NIR excess continuum, which has a color temperature of ~ 1000K, carries a luminosity of a few percent of the total far-infrared luminosity L(FIR), and most likely arises from the ISM; and (3) the well-known broad emission features at 6.2, 7.7, 8.6 and 11.3 um, which are generally attributed to aromatic carbon particles. These aromatic features in emission (AFEs) dominate the mid-infrared (5.8-11.6 um; MIR) part of the spectrum, and resemble the so-called Type-A spectra observed in many non-stellar sources and the diffuse ISM in our own Galaxy. The relative strengths of the AFEs vary by 15-25% among the galaxies. However, little correlation is seen between these variations and either IRAS 60um-to-100um flux density ratio R(60/100) or the FIR-to-blue luminosity ratio L(FIR)/L(B), suggesting that the observed variations are not a direct consequence of the radiation field differences among the galaxies. We demonstrate that the NIR excess continuum and AFE emission are correlated, suggesting that they are produced by similar mechanisms and similar (or the same) material. On the other hand, as the current star-formation activity increases, the overall strengths of the AFEs and the NIR excess continuum drop significantly with respect to that of the far-infrared emission from large dust grains. This is likely a consequence of the preferential destruction in intense radiation fields of the small carriers responsible for the NIR/AFE emission.

New far-infrared and submillimeter data are used to solidify and to extend to long wavelengths the empirical calibration of the infrared spectral energy distribution (SED) of normal star-forming galaxies. As was found by Dale et al. (2001), a single parameter family, characterized by f_nu(60 microns)/f_nu(100 microns), is adequate to describe the range of normal galaxy spectral energy distributions observed by IRAS and ISO from 3 to 100 microns. However, predictions based on the first generation models at longer wavelengths (122 to 850 microns) are increasingly overluminous compared to the data for smaller f_nu(60 microns)/f_nu(100 microns), or alternatively, for weaker global interstellar radiation fields. After slightly modifying the far-infrared/submillimeter dust emissivity in those models as a function of the radiation field intensity to better match the long wavelength data, a suite of SEDs from 3 microns to 20 cm in wavelength is presented. Results from relevant applications are also discussed, including submillimeter-based photometric redshift indicators, the infrared energy budget and simple formulae for recovering the bolometric infrared luminosity, and dust mass estimates in galaxies. Regarding the latter, since galaxy infrared SEDs are not well-described by single blackbody curves, the usual methods of estimating dust masses can be grossly inadequate. The improved model presented herein is used to provide a more accurate relation between infrared luminosity and dust mass.

(Abridged) We present the analysis of the main FIR fine structure lines emission in NGC1313 and NGC6946. We calculate that a component probably associated with the diffuse disks contributes <~40% in N6946 and ~30 % in N1313 to the total [CII] emission. The main PDR physical parameters responsible for the neutral atomic gas emission in N1313 and N6946 do not significantly differ from what Malhotra etal (2001) found by modelling the integrated emission of a sample of 60 normal galaxies,although there are evidences for a beam averaged contribution of a less active component inside NGC6946 higher than its contribution in the integrated emission of normal galaxies. CO and [CII] in N6946 are well correlated with a mean [CII]/CO ratio similar to that of the normal galaxies sample. In N1313 the [CII]/CO seems to systematically increase from the North to the south, along the S-shaped spiral arm, indicating much more inhomogeneous conditions than in N6946. HI and [CII] in N6946 are completely de-correlated, probably because they arise from different gas components. In N1313 we successfully detect two distinct gas components: a cirrus-like component where HI and [CII] are weakly correlated as observed in our Galaxy, and a component associated with dense PDRs completely de-correlated from HI as observed in N6946.Finally, we find that the HI residing in dense PDRs and presumably recently photo-dissociated, constitutes a few % of the total HI. In turn, this dense gas component produces most of the [CII] emission emitted by the atomic neutral medium.

We present Keck spectroscopy and UKIRT near-IR imaging observations of two 170micron-selected sources from the ISO-FIRBACK survey which have faint counterparts in the optical, and r-K~5. Both sources were expected to lie at z>1 based on their far-infrared, submillimeter and radio fluxes, assuming a similar spectral energy distribution to the local ultra-luminous infrared galaxy (ULIRG) Arp220. However, our spectroscopy indicates that the redshifts of these galaxies are z<1: z=0.91 for FN1-64 and z=0.45 for FN1-40. While the bolometric luminosities of both galaxies are similar to Arp220, it appears that the dust emission in these systems has a characteristic temperature of 30K much cooler than the ~50K seen in Arp220. Neither optical spectrum shows evidence of AGN activity. If these galaxies are characteristic of the optically faint FIRBACK population, then evolutionary models of the far-infrared background must include a substantial population of cold, luminous galaxies. These galaxies provide an important intermediate comparison between the local luminous IR galaxies, and the high redshift submillimeter-selected galaxies, for which there is very little information available.

New infrared and millimeter observations from Keck, Palomar, ISO, and OVRO and archival data from the NRAO VLA and IRAS are presented for the low metallicity blue compact dwarf galaxy SBS 0335-052. Mid-infrared imaging shows this young star-forming system is compact (0.31"; 80 pc) at 12.5 microns. The large Br-gamma equivalent width (235 Angstroms) measured from integral field spectroscopy is indicative of a ~5 Myr starburst. The central source appears to be optically thin in emission, containing both a warm (~80 K) and a hot (~210 K) dust component, and the overall interstellar radiation field is quite intense, about 10,000 times the intensity in the solar neighborhood. CO emission is not detected, though the galaxy shows an extremely high global H I gas-to-dust mass ratio, high even for blue compact dwarfs. Finally, the galaxy's mid-infrared-to-optical and mid-to-near-infrared luminosity ratios are quite high, whereas its far-infrared-to-radio and far-infrared-to-optical flux ratios are surprisingly similar to what is seen in normal star-forming galaxies. The relatively high bolometric infrared-to-radio ratio is more easily understood in the context of such a young system with negligible nonthermal radio continuum emission. These new lines of evidence may outline features common to primordial galaxies found at high redshift.

The most important cooling lines of the neutral interstellar medium (ISM) lie in the far-infrared (FIR). We present measurements by the Infrared Space Observatory Long Wavelength Spectrometer of seven lines from neutral and ionized ISM of 60 normal, star-forming galaxies. The galaxy sample spans a range in properties such as morphology, FIR colors (indicating dust temperature), and FIR/Blue ratios (indicating star-formation activity and optical depth). In two-thirds of the galaxies in this sample, the [CII] line is proportional to FIR dust continuum. The other one-third show a smooth decline in [CII]/FIR with increasing F60/F100 and FIR/B, spanning a range of a factor of more than 50. Two galaxies, at the warm and active extreme of the range have [CII]/FIR < 2 \times 10^{-4} (3-sigma upper limit). This is due to increased positive grain charge in the warmer and more active galaxies, which leads to less efficient heating by photoelectrons from dust grains. The ratio of the two principal photodissociation region (PDR) cooling lines [CII]/[OI] shows a tight correlation with F60/F100, indicating that both gas and dust temperatures increase together. We derive a theoretical scaling between [NII] and [CII] from ionized gas and use it to separate [CII] emission from neutral PDRs and ionized gas. Comparison of PDR models of Kaufman et al. (1999) with observed ratios of (a) [OI]/[CII] and ([CII]+[OI])/FIR and (b) [OI]/FIR and F60/F100 yields far-UV flux G0 and gas density n. The derived G0 scales as n to the power 1.4. We interpret this correlation as arising from Stromgren sphere scalings if much of the line and continuum luminosity arises near star-forming regions. The differences in G0 and n may be due to differences in the physical properties of the star-forming clouds.(Short abstract)

The environmental dependencies of the characteristics of spiral galaxy rotation curves are studied in this work. We use our large, homogeneously collected sample of 510 cluster spiral galaxy rotation curves to test the claim that the shape of a galaxy's rotation curve strongly depends on its location within the cluster, and thus presumably on the strength of the local intracluster medium and on the frequency and strength of tidal interactions with the cluster and cluster galaxies. Our data do not corroborate such a scenario, consistent with the fact that Tully-Fisher residuals are independent of galaxy location within the cluster; while the average late-type spiral galaxy shows more rise in the outer parts of its rotation curve than does the typical early-type spiral galaxy, there is no apparent trend for either subset with cluster environment. We also investigate as a function of cluster environment rotation curve asymmetry and the radial distribution of H II region tracers within galactic disks. Mild trends with projected cluster-centric distance are observed: (i) the (normalized) radial extent of optical line emission averaged over all spiral galaxy types shows a 4%+/-2% increase per Mpc of galaxy-cluster core separation, and (ii) rotation curve asymmetry falls by a factor of two between the inner and outer cluster for early-type spirals (a negligible decrease is found for late-type spirals). Such trends are consistent with spiral disk perturbations or even the stripping of the diffuse, outermost gaseous regions within the disks as galaxies pass through the dense cluster cores.

We present mid-infrared imaging and far-infrared (FIR) spectroscopy of 5 IBm galaxies observed by ISO as part of our larger study of the inter- stellar medium of galaxies. Most of the irregulars in our sample are very actively forming stars. The mid-infrared imaging was in a band centered at 6.75 microns that is dominated by polycyclic aromatic hydrocarbons (PAHs) and in a band centered at 15 microns that is dominated by small dust grains. The spectroscopy of 3 of the galaxies includes [CII]158 microns and [OI]63 microns, important coolants of photodissociation regions (PDRs), and [OIII]88 microns and [NII]122 microns, which come from ionized gas. [OI]145 microns and [OIII]52 microns were measured in one galaxy as well. These data are combined with PDR and HII region models to deduce properties of the interstellar medium of these galaxies.

We report a strong correlation between the [CII] 158 micron cooling line and the mid-infrared flux in the 5-10 micron range in a wide variety of star-forming galaxies. The mid-infrared flux is dominated by Aromatic Feature Emission (AFE), which is thought to arise from large polycyclic aromatic hydrocarbon molecules or `PAHs' and generally associated with the smallest interstellar grains. The [CII] line is the dominant gas coolant in most regions of atomic interstellar gas, and therefore reflects the heating input to the gas. The ratio of these two quantities, [CII]/AFE, remains nearly constant with the ratio of the IRAS 60 micron band flux to the 100 micron band flux, R(60/100). This is in contrast to the drop in the [CII]/FIR ratio with increasing R(60/100), which signal higher dust temperatures and more intense radiation fields. We interpret the stable [CII]/AFE ratio as evidence that gas heating is dominated by the PAHs or small grains which are also AFE carriers over a wide range of conditions. The trend of decreasing [CII]/FIR with increasing radiation field suggests a decrease in the importance of PAHs or small grains relative to large grains both in gas heating and in dust cooling. We summarize the observed trends and suggest two plausible scenarios.

We present a new phenomenological model for the spectral energy distribution of normal star-forming galaxies between 3 and 1100 microns. A sequence of realistic galaxy spectra are constructed from a family of dust emission curves assuming a power law distribution of dust mass over a wide range of interstellar radiation fields. For each interstellar radiation field heating intensity we combine emission curves for large and very small grains and aromatic feature carriers. The model is constrained by IRAS and ISOCAM broadband photometric and ISOPHOT spectrophotometric observations for our sample of 69 normal galaxies; the model reproduces well the empirical spectra and infrared color trends. These model spectra allow us to determine the infrared energy budget for normal galaxies, and in particular to translate far-infrared fluxes into total (bolometric) infrared fluxes. The 20 to 42 micron range appears to show the most significant growth in relative terms as the activity level increases, suggesting that the 20-42 micron continuum may be the best dust emission tracer of current star formation in galaxies. The redshift dependence of infrared color-color diagrams and the far-infrared to radio correlation for galaxies are also explored.

Four IRAS-detected early type galaxies were observed with ISO. With the exception of the 15 micron image of NGC1052, the mid-IR emission from NGC1052, NGC1155, NGC5866 and NGC6958 at 4.5, 7 and 15 microns show extended emission. Mid-IR emission from NGC1052, NGC1155, and NGC6958 follows a de Vaucouleurs profile. The ratio of 15/7 micron flux decreases with radius in these galaxies, approaching the values empirically observed for purely stellar systems. In NGC5866, the 7 and 15 micron emission is concentrated in the edge-on dust lane. All the galaxies are detected in the [CII] line, and the S0s NGC1155 and NGC5866 are detected in the [OI] line as well. The ISO-LWS observations of the [CII] line are more sensitive measures of cool, neutral ISM than HI and CO by about a factor of 10-100. Three of four early type galaxies, namely NGC1052, NGC6958 and NGC5866, have low ratio FIR/Blue and show a lower [CII]/FIR, which is due to a softer radiation field from old stellar populations. The low [CII]/CO ratio in NGC5866 ([CII]/CO(1-0) < 570) confirms this scenario. We estimate the UV radiation expected from the old stellar populations in these galaxies and compare it to that needed to heat the gas to account for the cooling observed [CII] and [OI] lines. In three out of four galaxies, NGC1052, NGC5866 and NGC6958, the predicted UV radiation falls short by a factor of 2-3. In view of the observed intrinsic scatter in the "UV-upturn" in elliptical galaxies and its great sensitivity to age and metallicity effects, this is not significant. However, the much larger difference (about a factor of 20) between the UV radiation from old stars and that needed to produce the FIR lines for NGC 1155 is strong evidence for the presence of young stars, in NGC1155.

We present an estimation of the expected number of arcs and arclets in a sample of nearby (z<0.1) clusters of galaxies, that takes into account the magnitude limit of the objects as well as seeing effects. We show that strong lensing effects are not common, but also they are not as rare as usually stated. Indeed, for a given cluster, they present a strong dependence with the magnitude limit adopted in the analysis and the seeing of the observations. We also describe the procedures and results of a search for lensing effects in a sample of 33 clusters spanning the redshift range of 0.014 to 0.076, representative of the local cluster distribution. This search produced two arc candidates. The first one is in A3408 (z=0.042), the same arc previously discovered by Campusano & Hardy (1996), with z=0.073 and associated to the brightest cluster galaxy. The second candidate is in the cluster A3266 (z=0.059) and is near a bright elliptical outside the cluster center, requiring the presence of a very massive sub-structure around this galaxy to be produced by gravitational lensing.

We present mid-infrared maps and preliminary analysis for 61 galaxies observed with the ISOCAM instrument aboard the Infrared Space Observatory. Many of the general features of galaxies observed at optical wavelengths---spiral arms, disks, rings, and bright knots of emission---are also seen in the mid-infrared, except the prominent optical bulges are absent at 6.75 and 15 microns. In addition, the maps are quite similar at 6.75 and 15 microns, except for a few cases where a central starburst leads to lower 6.75/15 ratios in the inner region. We also present infrared flux densities and mid-infrared sizes for these galaxies. The mid-infrared color 6.75/15 shows a distinct trend with the far-infrared color 60/100. The quiescent galaxies in our sample (60/100 < 0.6) show 6.75/15 near unity, whereas this ratio drops significantly for galaxies with higher global heating intensity levels. Azimuthally-averaged surface brightness profiles indicate the extent to which the mid-infrared flux is centrally concentrated, and provide information on the radial dependence of mid-infrared colors. The galaxies are mostly well resolved in these maps: almost half of them have < 10% of their flux in the central resolution element. A comparison of optical and mid-infrared isophotal profiles indicates that the optical flux at 4400 Angstroms near the optical outskirts of the galaxies is approximately eight (seven) times that at 6.75 microns (15 microns), comparable with observations of the diffuse quiescent regions of the Milky Way.

We investigate the rich cluster Abell 2029 (z~0.08) using optical imaging and long-slit spectral observations of 52 disk galaxies distributed throughout the cluster field. No strong emission-line galaxies are present within ~400 kpc of the cluster center, a region largely dominated by the similarly-shaped X-ray and low surface brightness optical envelopes centered on the giant cD galaxy. However, two-thirds of the galaxies observed outside the cluster core exhibit line emission. H-alpha rotation curves of 14 cluster members are used in conjunction with a deep I band image to study the environmental dependence of the Tully-Fisher relation. The Tully-Fisher zero-point of Abell 2029 matches that of clusters at lower redshifts, although we do observe a relatively larger scatter about the Tully-Fisher relation. We do not observe any systematic variation in the data with projected distance to the cluster center: we see no environmental dependence of Tully-Fisher residuals, R-I color, H-alpha equivalent width, and the shape and extent of the rotation curves.

We present Keck high-resolution near-IR (2.2 microns; FWHM~0.15") and mid-IR (12.5 microns; FWHM~0.4") images of APM08279+5255, a z=3.91 IR-luminous BALQSO with a prodigious apparent bolometric luminosity of 5x10^{15} Lsun, the largest known in the universe. The K-band image shows that this system consists of three components, all of which are likely to be the gravitationally lensed images of the same background object, and the 12.5 micron image shows a morphology consistent with such an image configuration. Our lens model suggests that the magnification factor is ~100 from the restframe UV to mid-IR, where most of the luminosity is released. The intrinsic bolometric luminosity and IR luminosity of APM08279+5255 are estimated to be 5x10^{13} Lsun and 1x10^{13} Lsun, respectively. This indicates that APM 08279+5255 is intriniscally luminous, but it is not the most luminous object known. As for its dust contents, little can be determined with the currently available data due to the uncertainties associated with the dust emissivity and the possible effects of differential magnification. We also suggest that the lensing galaxy is likely to be a massive galaxy at z~3.

We have obtained Tully-Fisher (TF) measurements for some 3000 late-type galaxies in the field and in 76 clusters distributed throughout the sky between 10 and 200\h Mpc. The cluster data are applied to the construction of an I band TF template, resulting in a relation with a scatter of 0.35 magnitudes and a zero-point accurate to 0.02 magnitudes. Peculiar motions are computed by referral to the template relation, and the distribution of line-of-sight cluster peculiar motions is presented. The dipole of the reflex motion of the Local Group of galaxies with respect to galaxies with measured peculiar velocity converges to the CMB dipole within less than 6000 km/s. The progression of this convergence is well illustrated when the reflex motion is referred to a well-distributed sample of field galaxies, and it is maintained when the reflex motion is referred to the reference frame constituted by the distant clusters in our sample. The field and cluster samples exhibit bulk motion amplitudes of order 200 km/s or smaller. Finally, we apply our sample of cluster peculiar velocities to a test of the putative Hubble bubble recently claimed by Zehavi and coworkers. In contrast to their findings, our data support a relatively quiescent Hubble flow beyond 35\h Mpc.

We report a mid-infrared color and surface brightness analysis of IC 10, NGC 1313, and NGC 6946, three of the nearby galaxies studied under the Infrared Space Observatory Key Project on Normal Galaxies. Images with < 9 arcsecond (170 pc) resolution of these nearly face-on, late-type galaxies were obtained using the LW2 (6.75 mu) and LW3 (15 mu) ISOCAM filters. Though their global I_nu(6.75 mu)/I_nu(15 mu) flux ratios are similar and typical of normal galaxies, they show distinct trends of this color ratio with mid-infrared surface brightness. We find that I_nu(6.75 mu)/I_nu(15 mu) ~< 1 only occurs for regions of intense heating activity where the continuum rises at 15 micron and where PAH destruction can play an important role. The shape of the color-surface brightness trend also appears to depend, to the second-order, on the hardness of the ionizing radiation. We discuss these findings in the context of a two-component model for the phases of the interstellar medium and suggest that star formation intensity is largely responsible for the mid-infrared surface brightness and colors within normal galaxies, whereas differences in dust column density are the primary drivers of variations in the mid-infrared surface brightness between the disks of normal galaxies.

We have obtained I band Tully-Fisher (TF) measurements for 522 late-type galaxies in the fields of 52 rich Abell clusters distributed throughout the sky between 50 and 200\h Mpc. Here we estimate corrections to the data for various forms of observational bias, most notably Malmquist and cluster population incompleteness bias. The bias-corrected data are applied to the construction of an I band TF template, resulting in a relation with a dispersion of 0.38 magnitudes and a kinematical zero-point accurate to 0.02 magnitudes. This represents the most accurate TF template relation currently available. Individual cluster TF relations are referred to the average template relation to compute cluster peculiar motions. The line-of-sight dispersion in the peculiar motions is 341+/-93 km/s, in general agreement with that found for the cluster sample of Giovanelli and coworkers.

We present Tully-Fisher observations for 35 rich Abell clusters of galaxies. Results from I band photometry and optical rotation curve work comprise the bulk of this paper. This is the third such data installment of an all-sky survey of 52 clusters in the distance range 50 to 200\h Mpc. The complete data set provides the basis for determining an accurate Tully-Fisher template relation and for estimating the amplitude and direction of the local bulk flow on a 100\h Mpc scale.

Zehavi et al. (1998) have suggested that the Hubble flow within 70/h Mpc may be accelerated by the existence of a void centered on the Local Group. Its underdensity would be ~20 %, which would result in a local Hubble distortion of about 6.5 %. We have combined the peculiar velocity data of two samples of clusters of galaxies, SCI and SCII, to investigate the amplitude of Hubble distortions to 200/h Mpc. Our results are not supportive of that conclusion. The amplitude of a possible distortion in the Hubble flow within 70/h Mpc in the SCI+SCII merged data is 0.010\pm0.022. The largest, and still quite marginal, geocentric deviation from smooth Hubble flow consistent with that data set is a shell with (Delta H)/H =0.027\pm0.023, centered at hd = 101 Mpc and extending over some 30/h Mpc. Our results are thus consistent with a Hubble flow that, on distances in excess of about 50/h Mpc, is remarkably smooth.

We have obtained new Tully-Fisher (TF) peculiar velocity measurements for 52 Abell galaxy clusters distributed throughout the sky between ~ 50 and 200 Mpc/h.The measurements are based on I band photometry and optical rotation curves for a sample of 522 spiral galaxies, from which an accurate TF template relation has been constructed. Individual cluster TF relations are referred to the template to compute cluster peculiar motions. The reflex motion of the Local Group of galaxies is measured with respect to the reference frame defined by our cluster sample and the distant portion of the Giovanelli et al. (1998) cluster set. We find the Local Group motion in this frame to be 565+/-113 km/s in the direction (l,b)=(267,26)+/-10 when peculiar velocities are weighted according to their errors. After optimizing the dipole calculation to sample equal volumes equally, the vector is 509+/-195 km/s towards (255,33)+/-22. Both solutions agree, to within 1-sigma or better, with the Local Group motion as inferred from the cosmic microwave background (CMB) dipole. Thus, the cluster sample as a whole moves slowly in the CMB reference frame, its bulk flow being at most 200 km/s.

We present Tully-Fisher (TF) observations for nine rich Abell clusters of galaxies. This is the second such data installment of an all-sky survey of 50 clusters in the redshift range 0.02 < z < 0.06. The data extends the TF study of nearby clusters of Giovanelli et al.; they will be used jointly to determine an accurate I band TF template and to establish a cluster inertial reference frame to z = 0.06.

We present first results of an all-sky observing program designed to improve the quality of the I band Tully-Fisher (TF) template and to obtain the reflex motion of the Local Group with respect to clusters to z = 0.06. We are obtaining between 5 and 15 TF measurements per cluster on a sample of 50 clusters at intermediate redshifts (0.02 < z < 0.06). Presentation of the data for seven Abell clusters of galaxies is given here. This data incorporates methods for estimating the true inclination of a spiral disk, an observational parameter undervalued for small angular-sized galaxies or for galaxies observed in poor seeing conditions.