We provide a coherent, uniform measurement of the evolution of the logarithmic star formation rate (SFR) - stellar mass ($M_*$) relation, called the main sequence of star-forming galaxies (MS), for galaxies out to $z\sim5$. We measure the MS using mean stacks of 3 GHz radio continuum images to derive average SFRs for $\sim$200,000 mass-selected galaxies at $z>0.3$ in the COSMOS field. We describe the MS relation adopting a new model that incorporates a linear relation at low stellar mass (log($M_*$/M$_\odot$)$<$10) and a flattening at high stellar mass that becomes more prominent at low redshift ($z<1.5$). We find that the SFR density peaks at $1.5<z<2$ and at each epoch there is a characteristic stellar mass ($M_* = 1 - 4 \times 10^{10}\mathrm{M}_\odot$) that contributes the most to the overall SFR density. This characteristic mass increases with redshift, at least to $z\sim2.5$. We find no significant evidence for variations in the MS relation for galaxies in different environments traced by the galaxy number density at $0.3<z<3$, nor for galaxies in X-ray groups at $z\sim0.75$. We confirm that massive bulge-dominated galaxies have lower SFRs than disk-dominated galaxies at a fixed stellar mass at $z<1.2$. As a consequence, the increase in bulge-dominated galaxies in the local star-forming population leads to a flattening of the MS at high stellar masses. This indicates that "mass-quenching" is linked with changes in the morphological composition of galaxies at a fixed stellar mass.

In the local (redshift z~0) Universe, collisional ring galaxies make up only ~0.01% of galaxies and are formed by head-on galactic collisions that trigger radially propagating density waves. These striking systems provide key snapshots for dissecting galactic disks and are studied extensively in the local Universe. However, not much is known about distant (z>0.1) collisional rings. Here we present a detailed study of a ring galaxy at a look-back time of 10.8 Gyr (z=2.19). Compared with our Milky Way, this galaxy has a similar stellar mass, but has a stellar half-light radius that is 1.5-2.2 times larger and is forming stars 50 times faster. The large, diffuse stellar light outside the star-forming ring, combined with a radial velocity on the ring and an intruder galaxy nearby, provides evidence for this galaxy hosting a collisional ring. If the ring is secularly evolved, the implied large bar in a giant disk would be inconsistent with the current understanding of the earliest formation of barred spirals. Contrary to previous predictions, this work suggests that massive collisional rings were as rare 11 Gyr ago as they are today. Our discovery offers a unique pathway for studying density waves in young galaxies, as well as constraining the cosmic evolution of spiral disks and galaxy groups.

We present kinematic orientations and high resolution (150 pc) rotation curves for 67 main sequence star-forming galaxies surveyed in CO (2-1) emission by PHANGS-ALMA. Our measurements are based on the application of a new fitting method tailored to CO velocity fields. Our approach identifies an optimal global orientation as a way to reduce the impact of non-axisymmetric (bar and spiral) features and the uneven spatial sampling characteristic of CO emission in the inner regions of nearby galaxies. The method performs especially well when applied to the large number of independent lines-of-sight contained in the PHANGS CO velocity fields mapped at 1'' resolution. The high resolution rotation curves fitted to these data are sensitive probes of mass distribution in the inner regions of these galaxies. We use the inner slope as well as the amplitude of our fitted rotation curves to demonstrate that CO is a reliable global dynamical mass tracer. From the consistency between photometric orientations from the literature and kinematic orientations determined with our method, we infer that the shapes of stellar disks in the mass range of log($\rm M_{\star}(M_{\odot})$)=9.0-10.9 probed by our sample are very close to circular and have uniform thickness.

We infer the intrinsic ionised gas kinematics for 383 star-forming galaxies across a range of integrated star-formation rates (SFR $\in [10^{-3}, 10^2]$ M$_\odot$ yr$^{-1}$) at $z \lesssim 0.1$ using a consistent 3D forward-modelling technique. The total sample is a combination of galaxies from the SAMI Galaxy Survey and DYNAMO survey. For typical low-$z$ galaxies taken from the SAMI Galaxy Survey, we find the vertical velocity dispersion ($\sigma_{v, z}$) to be positively correlated with measures of star-formation rate, stellar mass, HI gas mass, and rotational velocity. The greatest correlation is with star-formation rate surface density ($\Sigma_\text{SFR}$). Using the total sample, we find $\sigma_{v, z}$ increases slowly as a function of integrated star-formation rate in the range SFR $\in$ [$10^{-3}$, 1] M$_\odot$ yr$^{-1}$ from $17\pm3$ km s$^{-1}$ to $24\pm5$ km s$^{-1}$ followed by a steeper increase up to $\sigma_{v, z}$ $\sim 80$ km s$^{-1}$ for SFR $\gtrsim 1$ M$_\odot$ yr$^{-1}$. This is consistent with recent theoretical models that suggest a $\sigma_{v, z}$ floor driven by star-formation feedback processes with an upturn in $\sigma_{v, z}$ at higher SFR driven by gravitational transport of gas through the disc.

We investigate the stellar kinematics of the bulge and disk components in 826 galaxies with a wide range of morphology from the Sydney-AAO Multi-object Integral-field spectroscopy (SAMI) Galaxy Survey. The spatially-resolved rotation velocity (V) and velocity dispersion ($\sigma$) of bulge and disk components have been simultaneously estimated using the penalized pixel fitting (pPXF) method with photometrically defined weights for the two components. We introduce a new subroutine of pPXF for dealing with degeneracy in the solutions. We show that the V and $\sigma$ distributions in each galaxy can be reconstructed using the kinematics and weights of the bulge and disk components. The combination of two distinct components provides a consistent description of the major kinematic features of galaxies over a wide range of morphological types. We present Tully-Fisher and Faber-Jackson relations showing that the galaxy stellar mass scales with both V and $\sigma$ for both components of all galaxy types. We find a tight Faber-Jackson relation even for the disk component. We show that the bulge and disk components are kinematically distinct: (1) the two components show scaling relations with similar slopes, but different intercepts; (2) the spin parameter $\lambda_R$ indicates bulges are pressure-dominated systems and disks are supported by rotation; (3) the bulge and disk components have, respectively, low and high values in intrinsic ellipticity. Our findings suggest that the relative contributions of the two components explain, at least to first order, the complex kinematic behaviour of galaxies.

The physics regulating star formation (SF) in Hickson Compact Groups (HCG) has thus far been difficult to describe, due to their unique kinematic properties. In this study we expand upon previous works to devise a more physically meaningful SF relation able to better encompass the physics of these unique systems. We combine CO(1--0) data from the Combined Array from Research in Millimeter Astronomy (CARMA) to trace the column density of molecular gas $\Sigma_\mathrm{gas}$ and deep H$\alpha$ imaging taken on the Southern Astrophysical Research (SOAR) Telescope tracing $\Sigma_\mathrm{SFR}$ to investigate star formation efficiency across face-on HCG, NGC7674. We find a lack of universality in star formation, with two distinct sequences present in the $\Sigma_\mathrm{gas}-\Sigma_\mathrm{SFR}$ plane; one for inside and one for outside the nucleus. We devise a SF relation based on the multi-freefall nature of gas and the critical density, which itself is dependent on the virial parameter $\alpha_{\mathrm{vir}}$, the ratio of turbulent to gravitational energy. We find that our modified SF relation fits the data and describes the physics of this system well with the introduction of a virial parameter of about 5--10 across the galaxy. This $\alpha_{\mathrm{vir}}$ leads to an order-of-magnitude reduction in SFR compared to $\alpha_{\mathrm{vir}}\approx 1$ systems.

We present the detection of CO(5-4) with S/N> 7 - 13 and a lower CO transition with S/N > 3 (CO(4-3) for 4 galaxies, and CO(3-2) for one) with ALMA in band 3 and 4 in five main sequence star-forming galaxies with stellar masses 3-6x10^10 M/M_sun at 3 < z < 3.5. We find a good correlation between the total far-infrared luminosity LFIR and the luminosity of the CO(5-4) transition L'CO(5-4), where L'CO(5-4) increases with SFR, indicating that CO(5-4) is a good tracer of the obscured SFR in these galaxies. The two galaxies that lie closer to the star-forming main sequence have CO SLED slopes that are comparable to other star-forming populations, such as local SMGs and BzK star-forming galaxies; the three objects with higher specific star formation rates (sSFR) have far steeper CO SLEDs, which possibly indicates a more concentrated episode of star formation. By exploiting the CO SLED slopes to extrapolate the luminosity of the CO(1-0) transition, and using a classical conversion factor for main sequence galaxies of alpha_CO = 3.8 M_sun(K km s^-1 pc^-2)^-1, we find that these galaxies are very gas rich, with molecular gas fractions between 60 and 80%, and quite long depletion times, between 0.2 and 1 Gyr. Finally, we obtain dynamical masses that are comparable with the sum of stellar and gas mass (at least for four out of five galaxies), allowing us to put a first constraint on the alpha_CO parameter for main sequence galaxies at an unprecedented redshift.

We study the physical mechanisms that cause the offset between low-redshift and high-redshift galaxies on the [OIII]/H$\beta$ versus [NII]/H$\alpha$ ``Baldwin, Phillips & Terlevich'' (BPT) diagram using a sample of local analogues of high-redshift galaxies. These high-redshift analogue galaxies are selected from the Sloan Digital Sky Survey. Located in the same region on the BPT diagram as the ultra-violet selected galaxies at $z\sim2$, these high-redshift analogue galaxies provide an ideal local benchmark to study the offset between the local and high-redshift galaxies on the BPT diagram. We compare the nitrogen-to-oxygen ratio (N/O), the shape of the ionising radiation field, and ionisation parameters between the high-redshift analogues and a sample of local reference galaxies. The higher ionisation parameter in the high-redshift analogues is the dominant physical mechanism driving the BPT offset from low- to high-redshift, particularly at high {\nii/\ha}. Furthermore, the N/O ratio enhancement also plays a minor role to cause the BPT offset. However, the shape of the ionising radiation field is unlikely to cause the BPT offset because the high-redshift analogues have a similar hard ionising radiation field as local reference galaxies. This hard radiation field cannot be produced by the current standard stellar synthesis models. The stellar rotation and binarity may help solve the discrepancy.

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.

Recent integral field spectroscopy observations have found that about 11\% of galaxies show star-gas misalignment. The misalignment possibly results from external effects such as gas accretion, interaction with other objects, and other environmental effects, hence providing clues to these effects. We explore the properties of misaligned galaxies using Horizon-AGN, a large-volume cosmological simulation, and compare the result with the result of the Sydney-AAO Multi-object integral field spectrograph (SAMI) Galaxy Survey. Horizon-AGN can match the overall misalignment fraction and reproduces the distribution of misalignment angles found by observations surprisingly closely. The misalignment fraction is found to be highly correlated with galaxy morphology both in observations and in the simulation: early-type galaxies are substantially more frequently misaligned than late-type galaxies. The gas fraction is another important factor associated with misalignment in the sense that misalignment increases with decreasing gas fraction. However, there is a significant discrepancy between the SAMI and Horizon-AGN data in the misalignment fraction for the galaxies in dense (cluster) environments. We discuss possible origins of misalignment and disagreement.

It remains a major challenge to derive a theory of cloud-scale ($\lesssim100$ pc) star formation and feedback, describing how galaxies convert gas into stars as a function of the galactic environment. Progress has been hampered by a lack of robust empirical constraints on the giant molecular cloud (GMC) lifecycle. We address this problem by systematically applying a new statistical method for measuring the evolutionary timeline of the GMC lifecycle, star formation, and feedback to a sample of nine nearby disc galaxies, observed as part of the PHANGS-ALMA survey. We measure the spatially-resolved ($\sim100$ pc) CO-to-H$\alpha$ flux ratio and find a universal de-correlation between molecular gas and young stars on GMC scales, allowing us to quantify the underlying evolutionary timeline. GMC lifetimes are short, typically 10-30 Myr, and exhibit environmental variation, between and within galaxies. At kpc-scale molecular gas surface densities $\Sigma_{\rm H_2}\geqslant8$M$_{\odot}$pc$^{-2}$, the GMC lifetime correlates with time-scales for galactic dynamical processes, whereas at $\Sigma_{\rm H_2}\leqslant8$M$_{\odot}$pc$^{-2}$ GMCs decouple from galactic dynamics and live for an internal dynamical time-scale. After a long inert phase without massive star formation traced by H$\alpha$ (75-90% of the cloud lifetime), GMCs disperse within just 1-5 Myr once massive stars emerge. The dispersal is most likely due to early stellar feedback, causing GMCs to achieve integrated star formation efficiencies of 4-10% These results show that galactic star formation is governed by cloud-scale, environmentally-dependent, dynamical processes driving rapid evolutionary cycling. GMCs and HII regions are the fundamental units undergoing these lifecycles, with mean separations of 100-300 pc in star-forming discs. Future work should characterise the multi-scale physics and mass flows driving these lifecycles.

We present new measurements of the cosmic cold molecular gas evolution out to redshift 6 based on systematic mining of the ALMA public archive in the COSMOS deep field (A3COSMOS). Our A3COSMOS dataset contains ~700 galaxies (0.3 < z < 6) with high-confidence ALMA detections in the (sub-)millimeter continuum and multi-wavelength spectral energy distributions (SEDs). Multiple gas mass calibration methods are compared and biases in band conversions (from observed ALMA wavelength to rest-frame Rayleigh-Jeans(RJ)-tail continuum) have been tested. Combining our A3COSMOS sample with ~1,000 CO-observed galaxies at 0 < z < 4 (75% at z < 0.1), we parameterize galaxies' molecular gas depletion time and molecular gas to stellar mass ratio (gas fraction) each as a function of the stellar mass, offset from the star-forming main sequence (Delta MS) and cosmic age (or redshift). Our proposed functional form provides a statistically better fit to current data (than functional forms in the literature), and implies a "downsizing" effect (i.e., more-massive galaxies evolve earlier than less-massive ones) and "mass-quenching" (gas consumption slows down with cosmic time for massive galaxies but speeds up for low-mass ones). Adopting galaxy stellar mass functions and applying our function for gas mass calculation, we for the first time infer the cosmic cold molecular gas density evolution out to redshift 6 and find agreement with CO blind surveys as well as semi-analytic modeling. These together provide a coherent picture of cold molecular gas, SFR and stellar mass evolution in galaxies across cosmic time.

The rich information on (sub)millimeter dust continuum emission from distant galaxies in the public Atacama Large Millimeter/submillimeter Array (ALMA) archive is contained in thousands of inhomogeneous observations from individual PI-led programs. To increase the usability of these data for studies deepening our understanding of galaxy evolution, we have developed automated mining pipelines for the ALMA archive in the COSMOS field (A3COSMOS) that efficiently exploit the available information for large numbers of galaxies across cosmic time and keep the data products in sync with the increasing public ALMA archive: (a) a dedicated ALMA continuum imaging pipeline, (b) two complementary photometry pipelines for both blind source extraction and prior source fitting, (c) a counterpart association pipeline utilizing the multiwavelength data available (including quality assessment based on machine-learning techniques), (d) an assessment of potential (sub)millimeter line contribution to the measured ALMA continuum, and (e) extensive simulations to provide statistical corrections to biases and uncertainties in the ALMA continuum measurements. Application of these tools yields photometry catalogs with ~1000 (sub)millimeter detections (spurious fraction ~8%-12%) from over 1500 individual ALMA continuum images. Combined with ancillary photometric and redshift catalogs and the above quality assessments, we provide robust information on redshift, stellar mass, and star formation rate for ~700 galaxies at redshifts 0.5-6 in the COSMOS field (with undetermined selection function). The ALMA photometric measurements and galaxy properties are released publicly within our blind extraction, prior fitting, and galaxy property catalogs, plus the images. These products will be updated on a regular basis in the future.

We present an analysis of the dust attenuation of star forming galaxies at $z=2.5-4.0$ through the relationship between the UV spectral slope ($\beta$), stellar mass ($M_{\ast}$) and the infrared excess (IRX$=L_{\rm{IR}}/L_{\rm{UV}}$) based on far-infrared continuum observations from the Atacama Large Millimeter/sub-millimeter Array (ALMA). Our study exploits the full ALMA archive over the COSMOS field processed by the A$^3$COSMOS team, which includes an unprecedented sample of $\sim1500$ galaxies at $z\sim3$ as primary or secondary targets in ALMA band 6 or 7 observations with a median continuum sensitivity of 126 $\rm{\mu Jy/beam}$ (1$\sigma$). The detection rate is highly mass dependent, decreasing drastically below $\log (M_{\ast}/M_{\odot})=10.5$. The detected galaxies show that the IRX-$\beta$ relationship of massive ($\log M_{\ast}/M_{\odot} > 10$) main sequence galaxies at $z=2.5-4.0$ is consistent with that of local galaxies, while starbursts are generally offset by $\sim0.5\,{\rm dex}$ to larger IRX values. At the low mass end, we derive upper limits on the infrared luminosities through stacking of the ALMA data. The combined IRX-$M_{\ast}$ relation at $\rm{log\,(M_{\ast}/M_{\odot})>9}$ exhibits a significantly steeper slope than reported in previous studies at similar redshifts, implying little dust obscuration at $\log M_{\ast}/M_{\odot}<10$. However, our results are consistent with early measurements at $z\sim5.5$, indicating a potential redshift evolution between $z\sim2$ and $z\sim6$. Deeper observations targeting low mass galaxies will be required to confirm this finding.

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 distribution of metals within a galaxy traces the baryon cycle and the buildup of galactic disks, but the detailed gas phase metallicity distribution remains poorly sampled. We have determined the gas phase oxygen abundances for 7,138 HII regions across the disks of eight nearby galaxies using VLT/MUSE optical integral field spectroscopy as part of the PHANGS-MUSE survey. After removing the first order radial gradients present in each galaxy, we look at the statistics of the metallicity offset (Delta O/H) and explore azimuthal variations. Across each galaxy, we find low (sigma=0.03-0.05 dex) scatter at any given radius, indicative of efficient mixing. We compare physical parameters for those HII regions that are 1 sigma outliers towards both enhanced and reduced abundances. Regions with enhanced abundances have high ionization parameter, higher Halpha luminosity, lower Halpha velocity dispersion, younger star clusters and associated molecular gas clouds show higher molecular gas densities. This indicates recent star formation has locally enriched the material. Regions with reduced abundances show increased Halpha velocity dispersions, suggestive of mixing introducing more pristine material. We observe subtle azimuthal variations in half of the sample, but can not always cleanly associate this with the spiral pattern. Regions with enhanced and reduced abundances are found distributed throughout the disk, and in half of our galaxies we can identify subsections of spiral arms with clearly associated metallicity gradients. This suggests spiral arms play a role in organizing and mixing the ISM.

There is compelling evidence for a highly energetic Seyfert explosion (10^{56-57} erg) that occurred in the Galactic Centre a few million years ago. The clearest indications are the x-ray/gamma-ray "10 kpc bubbles" identified by the Rosat and Fermi satellites. In an earlier paper, we suggested another manifestation of this nuclear activity, i.e. elevated H-alpha emission along a section of the Magellanic Stream due to a burst (or flare) of ionizing radiation from Sgr A*. We now provide further evidence for a powerful flare event: UV absorption line ratios (in particular CIV/CII, SiIV/SiII) observed by the Hubble Space Telescope reveal that some Stream clouds towards both galactic poles are highly ionized by a source capable of producing ionization energies up to at least 50 eV. We show how these are clouds caught in a beam of bipolar, radiative "ionization cones" from a Seyfert nucleus associated with Sgr A*. In our model, the biconic axis is tilted by about 15 deg from the South Galactic Pole with an opening angle of roughly 60 deg. For the Stream at such large Galactic distances (D > 75 kpc), nuclear activity is a plausible explanation for all of the observed signatures: elevated H-alpha emission and H ionization fraction (X_e > 0.5), enhanced CIV/CII and SiIV/SiII ratios, and high CIV and SiIV column densities. Wind-driven "shock cones" are ruled out because the Fermi bubbles lose their momentum and energy to the Galactic corona long before reaching the Stream. The nuclear flare event must have had a radiative UV luminosity close to the Eddington limit (f_E ~ 0.1-1). Our time-dependent Seyfert flare models adequately explain the observations and indicate the Seyfert flare event took place T_o = 3.5 +/- 1 Myr ago. The timing estimates are consistent with the mechanical timescales needed to explain the x-ray/gamma-ray bubbles in leptonic jet/wind models (2-8 Myr).

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.

Hydrodynamical simulations predict that the jets of young radio sources can inhibit star formation in their host galaxies by injecting heat and turbulence into the interstellar medium (ISM). To investigate jet-ISM interactions in a galaxy with a young radio source, we have carried out a multi-wavelength study of the $z = 0.025$ Compact Steep Spectrum radio source hosted by the early-type galaxy UGC 05771. Using Keck/OSIRIS observations, we detected H\textsubscript{2} 1--0 S(1) and [Fe \textsc{ii}] emission at radii of 100s of pc, which traces shocked molecular and ionised gas being accelerated outwards by the jets to low velocities, creating a `stalling wind'. At kpc radii, we detected shocked ionised gas using observations from the CALIFA survey, covering an area much larger than the pc-scale radio source. We found that existing interferometric radio observations fail to recover a large fraction of the source's total flux, indicating the likely existence of jet plasma on kpc scales, which is consistent with the extent of shocked gas in the host galaxy. To investigate the star formation efficiency in UGC 05771, we obtained IRAM CO observations to analyse the molecular gas properties. We found that UGC 05771 sits below the Kennicutt-Schmidt relation, although we were unable to definitively conclude if direct interactions from the jets are inhibiting star formation. This result shows that jets may be important in regulating star formation in the host galaxies of compact radio sources.

In the optical spectra of galaxies, the separation of line emission from gas ionised by star formation and an AGN, or by star formation and shocks, are very well-understood problems. However, separating line emission between AGN and shocks has proven difficult. With the aid of a new three-dimensional diagnostic diagram, we show the simultaneous separation of line emission from star formation, shocks, and AGN in NGC 1068, and quantify the ratio of star formation, shocks, and AGN in each spaxel. The AGN, shock, and star formation luminosity distributions across the galaxy accurately align with X-ray, radio, and CO(3-2) observations, respectively. Comparisons with previous separation methods show that the shocked emission heavily mixes with the AGN emission. We also show that if the H$\alpha$ flux is to be used as a star formation rate indicator, separating line emission from as many sources as possible should be attempted to ensure accurate results.

We study and compare different examples of stellar evolutionary synthesis input parameters used to produce photoionisation model grids using the MAPPINGS V modelling code. The aim of this study is to (a) explore the systematic effects of various stellar evolutionary synthesis model parameters on the interpretation of emission lines in optical strong-line diagnostic diagrams, (b) characterise the combination of parameters able to reproduce the spread of local galaxies located in the star-forming region in the Sloan Digital Sky Survey, and (c) investigate the emission from extremely metal-poor galaxies using photoionisation models. We explore and compare the stellar input ionising spectrum (stellar population synthesis code [Starburst99, SLUG, BPASS], stellar evolutionary tracks, stellar atmospheres, star-formation history, sampling of the initial mass function) as well as parameters intrinsic to the H II region (metallicity, ionisation parameter, pressure, H II region boundedness). We also perform a comparison of the photoionisation codes MAPPINGS and CLOUDY. On the variations in the ionising spectrum model parameters, we find that the differences in strong emission-line ratios between varying models for a given input model parameter are small, on average ~0.1 dex. An average difference of ~0.1 dex in emission-line ratio is also found between models produced with MAPPINGS and CLOUDY. Large differences between the emission-line ratios are found when comparing intrinsic H II region parameters. We find that low-metallicity galaxies are better explained by a density-bounded H II region and higher pressures better encompass the spread of galaxies at high redshift.

We use data from the Sydney-AAO Multi-object Integral-field spectroscopy (SAMI) Galaxy Survey to study the dynamical scaling relation between galaxy stellar mass $M_*$ and the general kinematic parameter $S_K = \sqrt{K V_{rot}^2 + \sigma^2}$ that combines rotation velocity $V_{rot}$ and velocity dispersion $\sigma$. We show that the $\log M_* - \log S_K$ relation: (1)~is linear above limits set by properties of the samples and observations; (2)~has slightly different slope when derived from stellar or gas kinematic measurements; (3)~applies to both early-type and late-type galaxies and has smaller scatter than either the Tully-Fisher relation ($\log M_* - \log V_{rot}$) for late types or the Faber-Jackson relation ($\log M_* - \log\sigma$) for early types; and (4)~has scatter that is only weakly sensitive to the value of $K$, with minimum scatter for $K$ in the range 0.4 and 0.7. We compare $S_K$ to the aperture second moment (the `aperture velocity dispersion') measured from the integrated spectrum within a 3-arcsecond radius aperture ($\sigma_{3^{\prime\prime}}$). We find that while $S_{K}$ and $\sigma_{3^{\prime\prime}}$ are in general tightly correlated, the $\log M_* - \log S_K$ relation has less scatter than the $\log M_* - \log \sigma_{3^{\prime\prime}}$ relation.

The [CII]158$\mu$m line is one of the strongest far-infrared (FIR) lines and an important coolant in the interstellar medium of galaxies that is accessible out to high redshifts. The excitation of [CII] is complex and can best be studied in detail at low redshifts. Here we report the discovery of the highest global [CII] excess with respect to the FIR luminosity in the nearby AGN host galaxy HE1353-1917. This galaxy is exceptional among a sample of five targets because the AGN ionization cone and radio jet directly intercept the cold galactic disk. As a consequence, a massive multiphase gas outflow on kiloparsec scales is embedded in an extended narrow-line region. Because HE1353-1917 is distinguished by these special properties from our four bright AGN, we propose that a global [CII] excess in AGN host galaxies could be a direct signature of a multiphase AGN-driven outflow with a high mass-loading factor.

We present a systematic study of the diffuse ionized gas (DIG) in M83 and its effects on the measurement of metallicity gradients at varying resolution scales. Using spectrophotometric data cubes of M83 obtained at the 2.5m duPont telescope at Las Campanas Observatory as part of the TYPHOON program, we separate the HII regions from the DIG using the [SII]/H$\alpha$ ratio, HIIphot (HII finding algorithm) and the H$\alpha$ surface brightness. We find that the contribution to the overall H$\alpha$ luminosity is approximately equal for the HII and DIG regions. The data is then rebinned to simulate low-resolution observations at varying resolution scales from 41 pc up to 1005 pc. Metallicity gradients are measured using five different metallicity diagnostics at each resolution. We find that all metallicity diagnostics used are affected by the inclusion of DIG to varying degrees. We discuss the reasons of why the metallicity gradients are significantly affected by DIG using the HII dominance and emission line ratio radial profiles. We find that applying the [SII]/H$\alpha$ cut will provide a closer estimate of the true metallicity gradient up to a resolution of 1005 pc for all metallicity diagnostics used in this study.

SN1991bg-like supernovae are a distinct subclass of thermonuclear supernovae (SNe Ia). Their spectral and photometric peculiarities indicate their progenitors and explosion mechanism differ from `normal' SNe Ia. One method of determining information about supernova progenitors we cannot directly observe is to observe the stellar population adjacent to the apparent supernova explosion site to infer the distribution of stellar population ages and metallicities. We obtain integral field observations and analyse the spectra extracted from regions of projected radius $\sim\,\mathrm{kpc}$ about the apparent SN explosion site for 11 91bg-like SNe in both early- and late-type galaxies. We utilize full-spectrum spectral fitting to determine the ages and metallicities of the stellar population within the aperture. We find that the majority of the stellar populations that hosted 91bg-like supernovae have little recent star formation. The ages of the stellar populations suggest that that 91bg-like SN progenitors explode after delay times of $>6\,\mathrm{Gyr}$, much longer than the typical delay time of normal SNe Ia, which peaks at $\sim 1\,\mathrm{Gyr}$.

We present a novel Bayesian method, referred to as Blobby3D, to infer gas kinematics that mitigates the effects of beam smearing for observations using Integral Field Spectroscopy (IFS). The method is robust for regularly rotating galaxies despite substructure in the gas distribution. Modelling the gas substructure within the disk is achieved by using a hierarchical Gaussian mixture model. To account for beam smearing effects, we construct a modelled cube that is then convolved per wavelength slice by the seeing, before calculating the likelihood function. We show that our method can model complex gas substructure including clumps and spiral arms. We also show that kinematic asymmetries can be observed after beam smearing for regularly rotating galaxies with asymmetries only introduced in the spatial distribution of the gas. We present findings for our method applied to a sample of 20 star-forming galaxies from the SAMI Galaxy Survey. We estimate the global H$\alpha$ gas velocity dispersion for our sample to be in the range $\bar{\sigma}_v \sim $[7, 30] km s$^{-1}$. The relative difference between our approach and estimates using the single Gaussian component fits per spaxel is $\Delta \bar{\sigma}_v / \bar{\sigma}_v = - 0.29 \pm 0.18$ for the H$\alpha$ flux-weighted mean velocity dispersion.

In the optical spectra of galaxies, methods for the separation of line emission arising from star formation and an additional hard component, such as shocks or AGN, is well-understood and possible with current diagnostics. However, such diagnostics fail when attempting to separate and define line emission which arises from shocked gas, and that arising from AGN. We present a new three-dimensional diagnostic diagram for IFU data which can simultaneously separate the line emission amongst star formation, shocks, and AGN within a galaxy. We show that regions we define as AGN-dominated correlate well with the hard X-ray distribution in our test case NGC 1068, as well as with known regions of AGN activity in NGC 1068. Similarly, spaxels defined as shock-dominated correlate strongly with regions of high velocity dispersion within the galaxy.

We report the discovery of shocked molecular and ionized gas resulting from jet-driven feedback in the compact radio galaxy 4C 31.04 using near-IR imaging spectroscopy. 4C 31.04 is a $\sim 100$ pc double-lobed Compact Steep Spectrum source believed to be a very young AGN. It is hosted by a giant elliptical with a $\sim 10^{9}~\rm M_\odot$ multi-phase gaseous circumnuclear disc. We used high spatial resolution, adaptive optics-assisted $H$- and $K$-band integral field Gemini/NIFS observations to probe (1) the warm ($\sim 10^3~\rm K$) molecular gas phase, traced by ro-vibrational transitions of $\rm H_2$, and (2), the warm ionized medium, traced by the [Fe II]$_{1.644~\rm \mu m}$ line. The [Fe II] emission traces shocked gas ejected from the disc plane by a jet-blown bubble $300-400~\rm pc$ in diameter, whilst the $\rm H_2$ emission traces shock-excited molecular gas in the interior $\sim 1~\rm kpc$ of the circumnuclear disc. Hydrodynamical modelling shows that the apparent discrepancy between the extent of the shocked gas and the radio emission can occur when the brightest regions of the synchrotron-emitting plasma are temporarily halted by dense clumps, whilst less bright plasma can percolate through the porous ISM and form an energy-driven bubble that expands freely out of the disc plane. This bubble is filled with low surface-brightness plasma not visible in existing VLBI observations of 4C 31.04 due to insufficient sensitivity. Additional radial flows of jet plasma may percolate to $\sim \rm kpc$ radii in the circumnuclear disc, driving shocks and accelerating clouds of gas, giving rise to the $\rm H_2$ emission.

We systematically measure the gas-phase metallicities and the mass-metallicity relation of a large sample of local active galaxies for the first time. Observed emission-line fluxes from the Sloan Digital Sky Survey (SDSS) are compared to a four-dimensional grid of photoionization models using the Bayesian parameter estimation code NebulaBayes. For the first time we take into account arbitrary mixing between HII region and narrow-line region (NLR) emission, and the models are also varied with metallicity, ionization parameter in the NLR, and the gas pressure. The active galactic nucleus (AGN) oxygen abundance is found to increase by $\Delta {\rm O/H} \sim 0.1$ dex as a function of host galaxy stellar mass over the range $10.1 < \log M_* / M_\odot < 11.3$. We also measure the metallicity and ionization parameter of 231000 star-forming galaxies for comparison with the sample of 7670 Seyfert 2 galaxies. A systematic offset in oxygen abundance of 0.09 dex is observed between the mass-metallicity relations of the star-forming and active galaxies. We investigate potential causes of the offset, including sample selection and the treatment in the models of diffuse ionized gas, pressure, and ionization parameter. We cannot identify the major cause(s), but suspect contributions due to deficiencies in modeling the ionizing spectra and the treatment of dust physics. Optical diagnostic diagrams are presented with the star-forming and Seyfert data colored by the inferred oxygen abundance, ionization parameter and gas pressure, clearly illustrating the trends in these quantities.

We calibrate commonly used star formation rate (SFR) prescriptions using observations in five kpc-sized fields in the nearby galaxy Andromeda (M31) at 10\,pc spatial resolution. Our observations at different scales enable us to resolve the star-forming regions and to distinguish them from non star-forming components. We use extinction corrected H$\alpha$ from optical integral field spectroscopy as our reference tracer and have verified its reliability via tests. It is used to calibrate monochromatic and hybrid (H$\alpha$+a$\times$IR and FUV+b$\times$IR) SFR prescriptions, which use FUV (GALEX), 22\,$\mu$m (WISE) and 24\,$\mu$m (MIPS). Additionally, we evaluate other multi-wavelength infra-red tracers. Our results indicate that the SFR prescriptions do not change (in M31) with spatial scales or with subtraction of the diffuse component. For the calibration factors in the hybrid SFR prescriptions, we find a$\approx$0.2 and b$\approx$22 in M31, which are a factor of 5 higher than in the literature. As the fields in M31 exhibit high attenuation and low dust temperatures, lie at large galacto-centric distances, and suffer from high galactic inclination compared to measurements in other galaxies, we propose that the fields probe a dust layer extended along the line of sight that is not directly spatially associated with star-forming regions. This (vertically) extended dust component increases the attenuation and alters the SFR prescriptions in M31 compared to literature measurements. We recommend that SFR prescriptions should be applied with caution at large galacto-centric distances and in highly inclined galaxies, due to variations in the relative (vertical) distribution of dust and gas.

We use integral field spectroscopy from the SAMI Galaxy Survey to identify galaxies that show evidence for recent quenching of star formation. The galaxies exhibit strong Balmer absorption in the absence of ongoing star formation in more than 10% of their spectra within the SAMI field of view. These $\rm{H}{\delta}$-strong galaxies (HDSGs) are rare, making up only $\sim 2$% (25/1220) of galaxies with stellar mass ${\rm log(}M_*/M_{\odot})>10$. The HDSGs make up a significant fraction of non-passive cluster galaxies (15%; 17/115) and a smaller fraction (2.0%; 8/387) of the non-passive population in low-density environments. The majority (9/17) of cluster HDSGs show evidence for star formation at their centers, with the HDS regions found in the outer parts of the galaxy. Conversely, the $\rm{H}{\delta}$-strong signal is more evenly spread across the galaxy for the majority (6/8) of HDSGs in low-density environments, and is often associated with emission lines that are not due to star formation. We investigate the location of the HDSGs in the clusters, finding that they are exclusively within 0.6$R_{200}$ of the cluster centre, and have a significantly higher velocity dispersion relative to the cluster population. Comparing their distribution in projected-phase-space to those derived from cosmological simulations indicates that the cluster HDSGs are consistent with an infalling population that have entered the central 0.5$r_{200, 3D}$ cluster region within the last $\sim 1\,$Gyr. In the 8/9 cluster HDSGs with central star formation, the extent of star formation is consistent with that expected of outside-in quenching by ram-pressure stripping. Our results indicate that the cluster HDSGs are currently being quenched by ram-pressure stripping on their first passage through the cluster.

We explore the radial distribution of star formation in galaxies in the SAMI Galaxy Survey as a function of their local group environment. Using a sample of galaxies in groups (with halo masses less than $ \simeq 10^{14} \, \mathrm{M_{\odot}}$) from the Galaxy And Mass Assembly Survey, we find signatures of environmental quenching in high-mass groups ($M_{G} > 10^{12.5} \, \mathrm{M_{\odot}}$). The mean integrated specific star formation rate of star-forming galaxies in high-mass groups is lower than for galaxies in low-mass groups or that are ungrouped, with $\Delta \log(sSFR/\mathrm{yr^{-1}}) = 0.45 \pm 0.07$. This difference is seen at all galaxy stellar masses. In high-mass groups, star-forming galaxies more massive than $M_{*} \sim 10^{10} \, \mathrm{M_{\odot}}$ have centrally-concentrated star formation. These galaxies also lie below the star-formation main sequence, suggesting they may be undergoing outside-in quenching. Lower mass galaxies in high-mass groups do not show evidence of concentrated star formation. In groups less massive than $M_{G} = 10^{12.5} \, \mathrm{M_{\odot}}$ we do not observe these trends. In this regime we find a modest correlation between centrally-concentrated star formation and an enhancement in total star formation rate, consistent with triggered star formation in these galaxies.

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.

We evaluate dust-corrected far ultraviolet (FUV) star formation rates (SFRs) for samples of star-forming galaxies at $z\sim0$ and $z\sim0.7$ and find significant differences between values obtained through corrections based on UV colour, from a hybrid mid-infrared (MIR) plus FUV relation, and from a radiative transfer based attenuation correction method. The performances of the attenuation correction methods are assessed by their ability to remove the dependency of the corrected SFR on inclination, as well as returning, on average, the expected population mean SFR. We find that combining MIR (rest-frame $\sim$13$\mu$m) and FUV luminosities gives the most inclination independent SFRs and reduces the intrinsic SFR scatter out of the methods tested. However, applying the radiative transfer based method of Tuffs et al. gives corrections to the FUV SFR that are inclination independent and in agreement with the expected SFRs at both $z\sim0$ and $z\sim0.7$. SFR corrections based on the UV-slope perform worse than the other two methods tested. For our local sample, the UV-slope method works on average but does not remove inclination biases. At $z\sim$0.7 we find that the UV-slope correction used locally flattens the inclination dependence compared to the raw FUV measurements but was not sufficient to correct for the large attenuation observed at $z\sim$0.7.

Disk galaxies at intermediate redshift ($z\sim0.7$) have been found in previous work to display more optically thick behaviour than their local counterparts in the rest-frame B-band surface brightness, suggesting an evolution in dust properties over the past $\sim$6 Gyr. We compare the measured luminosities of face-on and edge-on star-forming galaxies at different wavelengths (Ultraviolet (UV), mid-infrared (MIR), far-infrared (FIR), and radio) for two well-matched samples of disk-dominated galaxies: a local Sloan Digital Sky Survey (SDSS)-selected sample at $z\sim0.07$ and a sample of disks at $z\sim0.7$ drawn from Cosmic Evolution Survey (COSMOS). We have derived correction factors to account for the inclination dependence of the parameters used for sample selection. We find that typical galaxies are transparent at MIR wavelengths at both redshifts and that the FIR and radio emission is also transparent as expected. However, reduced sensitivity at these wavelengths limits our analysis; we cannot rule out opacity in the FIR or radio. Ultra-violet attenuation has increased between $z\sim0$ and $z\sim0.7$, with the $z\sim0.7$ sample being a factor of $\sim$3.4 more attenuated. The larger UV attenuation at $z\sim0.7$ can be explained by more clumpy dust around nascent star-forming regions. There is good agreement between the fitted evolution of the normalisation of the SFR$_{\text{UV}}$ versus 1-cos(i) trend (interpreted as the clumpiness fraction) and the molecular gas fraction/dust fraction evolution of galaxies found out to $z<1$.

We present the second major release of data from the SAMI Galaxy Survey. Data Release Two includes data for 1559 galaxies, about 50% of the full survey. Galaxies included have a redshift range 0.004 < z < 0.113 and a large stellar mass range 7.5 < log (M_star/M_sun) < 11.6. The core data for each galaxy consist of two primary spectral cubes covering the blue and red optical wavelength ranges. For each primary cube we also provide three spatially binned spectral cubes and a set of standardised aperture spectra. For each core data product we provide a set of value-added data products. This includes all emission line value-added products from Data Release One, expanded to the larger sample. In addition we include stellar kinematic and stellar population value-added products derived from absorption line measurements. The data are provided online through Australian Astronomical Optics' Data Central. We illustrate the potential of this release by presenting the distribution of ~350,000 stellar velocity dispersion measurements from individual spaxels as a function of R/R_e, divided in four galaxy mass bins. In the highest stellar mass bin (log (M_star/M_sun)>11), the velocity dispersion strongly increases towards the centre, whereas below log (M_star/M_sun)<10 we find no evidence for a clear increase in the central velocity dispersion. This suggests a transition mass around log (M_star/M_sun) ~10 for galaxies with or without a dispersion-dominated bulge.

We demonstrate a robust method of resolving the star-formation and AGN contributions to emission lines using two very well known AGN systems: NGC 1365, and NGC 1068, using the high spatial resolution data from the TYPHOON/PrISM survey. We expand the previous method of calculating the AGN fraction by using theoretical-based model grids rather than empirical points. The high spatial resolution of the TYPHOON/PrISM observations show evidence of both star formation and AGN activity occurring in the nuclei of the two galaxies. We rebin the data to the lower resolutions, typically found in other integral field spectroscopy surveys such as SAMI, MaNGA, and CALIFA. The results show that when rebinned from the native resolution of TYPHOON (< 200 pc/pixel) to 1 kpc/pixel, the effects include a roughly 3 kpc increase in the radius of measured AGN activity, and a factor of 2 to 7 increase in the detection of low surface brightness features such as shocks. All of this information is critical, because information on certain physical processes may be lost at varying resolutions. We make recommendations for analysing data at current IFU survey resolutions.

Star formation is a multi-scale process that requires tracing cloud formation and stellar feedback within the local (<kpc) and global galaxy environment. We present first results from two large observing programs on ALMA and VLT/MUSE, mapping cloud scales (1arcsec = 47pc) in both molecular gas and star forming tracers across 90 kpc^2 of the central disk of NGC 628 to probe the physics of star formation. Systematic spatial offsets between molecular clouds and HII regions illustrate the time evolution of star-forming regions. Using uniform sampling of both maps on 50-500 pc scales, we infer molecular gas depletion times of 1-3 Gyr, but also find the increase of scatter in the star formation relation on small scales is consistent with gas and HII regions being only weakly correlated at the cloud (50 pc) scale. This implies a short overlap phase for molecular clouds and HII regions, which we test by directly matching our catalog of 1502 HII regions and 738 GMCs. We uncover only 74 objects in the overlap phase, and we find depletion times >1 Gyr, significantly longer than previously reported for individual star-forming clouds in the Milky Way. Finally, we find no clear trends that relate variations in the depletion time observed on 500 pc scales to physical drivers (metallicity, molecular and stellar mass surface density, molecular gas boundedness) on 50 pc scales.

Evolved asymptotic giant branch (AGB) stars and Type Ia supernovae (SNe) are important contributors to the elements that form dust in the interstellar medium of galaxies, in particular, carbon and iron. However, they require at least a Gyr to start producing these elements, therefore, a change in dust quantity or properties may appear at high redshifts. In this work, we use extinction of gamma-ray burst (GRB) afterglows as a tool to look for variations in dust properties at z>3. We use a spectroscopically selected sample of GRB afterglows observed with the VLT/X-shooter instrument to determine extinction curves out to high redshifts. We present ten new z>3 X-shooter GRBs of which six are dusty. Combining these with individual extinction curves of three previously known z>3 GRBs, we find an average extinction curve consistent with the SMC-Bar. A comparison with spectroscopically selected GRBs at all redshifts indicates a drop in visual extinction (A_V) at z>3.5 with no moderate or high extinction bursts. We check for observational bias using template spectra and find that GRBs up to z~8 are detectable with X-shooter up to A_V~0.3 mag. Although other biases are noted, a uniformly low dust content above z>3.5 indicates a real drop, suggesting a transition in dust properties and/or available dust building blocks. The remarkable increase in dust content at z<3.5 could arise due to carbon and possibly iron production by the first carbon-rich AGB and Type Ia SNe, respectively. Alternatively, z>3.5 dust drop could be the result of low stellar masses of GRB host galaxies.

The azimuthal variation of the HII region oxygen abundance in spiral galaxies is a key observable for understanding how quickly oxygen produced by massive stars can be dispersed within the surrounding interstellar medium. Observational constraints on the prevalence and magnitude of such azimuthal variations remain rare in the literature. Here, we report the discovery of pronounced azimuthal variations of HII region oxygen abundance in NGC 2997, a spiral galaxy at approximately 11.3 Mpc. Using 3D spectroscopic data from the TYPHOON Program, we study the HII region oxygen abundance at a physical resolution of 125 pc. Individual HII regions or complexes are identified in the 3D optical data and their strong emission line fluxes measured to constrain their oxygen abundances. We find 0.06 dex azimuthal variations in the oxygen abundance on top of a radial abundance gradient that is comparable to those seen in other star-forming disks. At a given radial distance, the oxygen abundances are highest in the spiral arms and lower in the inter-arm regions, similar to what has been reported in NGC 1365 using similar observations. We discuss whether the azimuthal variations could be recovered when the galaxy is observed at worse physical resolutions and lower signal-to-noise ratios.

We present gas-phase metallicity and ionization parameter maps of 25 star-forming face-on spiral galaxies from the SAMI Galaxy Survey Data Release 1. Self-consistent metallicity and ionization parameter maps are calculated simultaneously through an iterative process to account for the interdependence of the strong emission line diagnostics involving ([OII]+[OIII])/H$\beta$ (R23) and [OIII]/[OII] (O32). The maps are created on a spaxel-by-spaxel basis because HII regions are not resolved at the SAMI spatial resolution. We combine the SAMI data with stellar mass, star formation rate (SFR), effective radius (R$_e$), ellipticity, and position angles (PA) from the GAMA survey to analyze their relation to the metallicity and ionization parameter. We find a weak trend of steepening metallicity gradient with galaxy stellar mass, with values ranging from -0.03 to -0.20 dex/R$_e$. Only two galaxies show radial gradients in ionization parameter. We find that the ionization parameter has no significant correlation with either SFR, sSFR (specific star formation rate), or metallicity. For several individual galaxies we find structure in the ionization parameter maps suggestive of spiral arm features. We find a typical ionization parameter range of $7.0 < \log(q) < 7.8$ for our galaxy sample with no significant overall structure. An ionization parameter range of this magnitude is large enough to caution the use of metallicity diagnostics which have not considered the effects of a varying ionization parameter distribution.

We fit theoretical models to the emission-line spectra of 2766 Seyfert galaxies from the Sloan Digital Sky Survey (SDSS). The model line fluxes are derived by `mixing' photoionization model predictions of active galactic nucleus (AGN) narrow line region (NLR) emission and HII region emission. The observed line fluxes in each spectrum were directly compared to the grid of mixed models using the Bayesian parameter estimation code NebulaBayes, thereby measuring the degree of mixing in each spectrum for the first time. We find that the majority of the Balmer line emission in the majority of Seyfert-classified SDSS spectra arises from contaminating HII regions within the fixed-size aperture. Even for spectra with log [OIII]/Hb $\gtrsim 0.9$, a median of ~30% of the Balmer flux arises in HII regions. We show that the results are qualitatively insensitive to the assumed Seyfert ionizing continuum, and that ionizing spectra with a peak energy of E_peak ~ 40-50 eV produce the most plausible distributions of mixing fractions. The analysis cleanly quantifies how the starburst - AGN `mixing fraction' increases on the BPT diagram for SDSS galaxies. Apart from the mixing fraction, the models also vary with the ionization parameter in the NLR, the gas pressure, and the metallicity. Measurements for the metallicity in particular will be presented in a companion paper.

Although studying outflows in the host galaxies of AGN have become the forefront of extra-galactic astronomy in recent years, estimating the energy associated with these outflows have been a major challenge. Determination of the energy associated with an outflow often involves an assumption of uniform density in the NLR, which span a wide range in literature leading to large systematic uncertainties in energy estimation. In this paper, we present electron density maps for a sample of outflowing and non-outflowing Seyfert galaxies at z<0.02 drawn from the S7 survey and understand the origin and values of the observed density structures to reduce the systematic uncertainties in outflow energy estimation. We use the ratio of the [SII]6716,6731 emission lines to derive spatially resolved electron densities (<50-2000 cm$^{-3}$). Using optical IFU observations, we are able to measure densities across the central 2-5 kpc of the selected AGN host galaxies. We compare the density maps with the positions of the HII regions derived from the narrow H$\alpha$ component, ionization maps from [OIII], and spatially resolved BPT diagrams, to infer the origin of the observed density structures. We also use the electron density maps to construct density profiles as a function of distance from the central AGN. We find a spatial correlation between the sites of high star formation and high electron density for targets without an active ionized outflow. The non-outflowing targets also show an exponential drop in the electron density as a function of distance from the center, with a mean exponential index of ~0.15. The correlation between the star forming sites and electron density ceases for targets with an outflow. The density within the outflowing medium is not uniform and shows both low and high density sites, most likely due to the presence of shocks and highly turbulent medium.

We present a comparative study of the relation between the aperture-based gas-phase metallicity and three structural parameters of star-forming galaxies: mass ($\mathrm{M \equiv M_*}$), average potential ($\Phi \equiv \mathrm{M_*/R_e}$) and average surface mass density ($\Sigma \equiv \mathrm{M_*/R_e^2}$; where $\mathrm{R_e}$ is the effective radius). We use a volume-limited sample drawn from the publicly available SDSS DR7, and base our analysis on aperture-matched sampling by selecting sets of galaxies where the SDSS fibre probes a fixed fraction of $\mathrm{R_e}$. We find that between 0.5 and 1.5 $\mathrm{R_e}$, the gas-phase metallicity correlates more tightly with $\Phi$ than with either $\mathrm{M}$ or $\Sigma$, in that for all aperture-matched samples, the potential-metallicity relation has (i) less scatter, (ii) higher Spearman rank correlation coefficient and (iii) less residual trend with $\mathrm{R_e}$ than either the mass-metallicity relation and the average surface density-metallicity relation. Our result is broadly consistent with the current models of gas enrichment and metal loss. However, a more natural explanation for our findings is a local relation between the gas-phase metallicity and escape velocity.

We measure the velocity dispersion, $\sigma$, and surface density, $\Sigma$, of the molecular gas in nearby galaxies from CO spectral line cubes with spatial resolution $45$-$120$ pc, matched to the size of individual giant molecular clouds. Combining $11$ galaxies from the PHANGS-ALMA survey with $4$ targets from the literature, we characterize ${\sim}30,000$ independent sightlines where CO is detected at good significance. $\Sigma$ and $\sigma$ show a strong positive correlation, with the best-fit power law slope close to the expected value for resolved, self-gravitating clouds. This indicates only weak variation in the virial parameter $\alpha_\mathrm{vir}\propto\sigma^2/\Sigma$, which is ${\sim}1.5$-$3.0$ for most galaxies. We do, however, observe enormous variation in the internal turbulent pressure $P_\mathrm{turb}\propto\Sigma\,\sigma^2$, which spans ${\sim}5\rm\;dex$ across our sample. We find $\Sigma$, $\sigma$, and $P_\mathrm{turb}$ to be systematically larger in more massive galaxies. The same quantities appear enhanced in the central kpc of strongly barred galaxies relative to their disks. Based on sensitive maps of M31 and M33, the slope of the $\sigma$-$\Sigma$ relation flattens at $\Sigma\lesssim10\rm\;M_\odot\,pc^{-2}$, leading to high $\sigma$ for a given $\Sigma$ and high apparent $\alpha_\mathrm{vir}$. This echoes results found in the Milky Way, and likely originates from a combination of lower beam filling factors and a stronger influence of local environment on the dynamical state of molecular gas in the low density regime.

We investigate the evolution of the ionisation parameter of star-forming galaxies using a high-redshift ($z\sim 1.5$) sample from the FMOS-COSMOS survey and matched low-redshift samples from the Sloan Digital Sky Survey. By constructing samples of low-redshift galaxies for which the stellar mass ($\mathrm{M}_{*}$), star formation rate (SFR) and specific star formation rate (sSFR) are matched to the high-redshift sample we remove the effects of an evolution in these properties. We also account for the effect of metallicity by jointly constraining the metallicity and ionisation parameter of each sample. We find an evolution in the ionisation parameter for main-sequence, star-forming galaxies and show that this evolution is driven by the evolution of sSFR. By analysing the matched samples as well as a larger sample of $z<0.3$, star-forming galaxies we show that high ionisation parameters are directly linked to high sSFRs and are not simply the byproduct of an evolution in metallicity. Our results are physically consistent with the definition of the ionisation parameter, a measure of the hydrogen ionising photon flux relative to the number density of hydrogen atoms.

NebulaBayes is a new Bayesian code that implements a general method of comparing observed emission-line fluxes to photoionization model grids. The code enables us to extract robust, spatially resolved measurements of abundances in the extended narrow line regions (ENLRs) produced by Active Galactic Nuclei (AGN). We observe near-constant ionization parameters but steeply radially-declining pressures, which together imply that radiation pressure regulates the ENLR density structure on large scales. Our sample includes four `pure Seyfert' galaxies from the S7 survey that have extensive ENLRs. NGC2992 shows steep metallicity gradients from the nucleus into the ionization cones. An {\it inverse} metallicity gradient is observed in ESO138-G01, which we attribute to a recent gas inflow or minor merger. A uniformly high metallicity and hard ionizing continuum are inferred across the ENLR of Mrk573. Our analysis of IC5063 is likely affected by contamination from shock excitation, which appears to soften the inferred ionizing spectrum. The peak of the ionizing continuum E_peak is determined by the nuclear spectrum and the absorbing column between the nucleus and the ionized nebula. We cannot separate variation in this intrinsic E_peak from the effects of shock or HII region contamination, but E_peak measurements nevertheless give insights into ENLR excitation. We demonstrate the general applicability of NebulaBayes by analyzing a nuclear spectrum from the non-active galaxy NGC4691 using a HII region grid. The NLR and HII region model grids are provided with NebulaBayes for use by the astronomical community.

We present the ~800 star formation rate maps for the SAMI Galaxy Survey based on H{\alpha} emission maps, corrected for dust attenuation via the Balmer decrement, that are included in the SAMI Public Data Release 1. We mask out spaxels contaminated by non-stellar emission using the [O III]/H{\beta}, [N II]/H{\alpha}, [S II]/H{\alpha}, and [O I]/H{\alpha} line ratios. Using these maps, we examine the global and resolved star-forming main sequences of SAMI galaxies as a function of morphology, environmental density, and stellar mass. Galaxies further below the star-forming main sequence are more likely to have flatter star formation profiles. Early-type galaxies split into two populations with similar stellar masses and central stellar mass surface densities. The main sequence population has centrally-concentrated star formation similar to late-type galaxies, while galaxies >3{\sigma} below the main sequence show significantly reduced star formation most strikingly in the nuclear regions. The split populations support a two-step quenching mechanism, wherein halo mass first cuts off the gas supply and remaining gas continues to form stars until the local stellar mass surface density can stabilize the reduced remaining fuel against further star formation. Across all morphologies, galaxies in denser environments show a decreased specific star formation rate from the outside in, supporting an environmental cause for quenching, such as ram-pressure stripping or galaxy interactions.

We report an integral-field spectroscopic (IFS) observation of a gravitationally lensed spiral galaxy A1689B11 at redshift $z=2.54$. It is the most ancient spiral galaxy discovered to date and the second kinematically confirmed spiral at $z\gtrsim2$. Thanks to gravitational lensing, this is also by far the deepest IFS observation with the highest spatial resolution ($\sim$ 400 pc) on a spiral galaxy at a cosmic time when the Hubble sequence is about to emerge. After correcting for a lensing magnification of 7.2 $\pm$ 0.8, this primitive spiral disk has an intrinsic star formation rate of 22 $\pm$ 2 $M_{\odot}$ yr$^{-1}$, a stellar mass of 10$^{9.8 \pm 0.3}$$M_{\odot}$ and a half-light radius of $r_{1/2}=2.6 \pm 0.7$ kpc, typical of a main-sequence star-forming (SF) galaxy at $z\sim2$. However, the H\alpha\ kinematics show a surprisingly tranquil velocity field with an ordered rotation ($V_{\rm c}$ = 200 $\pm$ 12 km/s) and uniformly small velocity dispersions ($V_{\rm \sigma, mean}$ = 23 $\pm$ 4 km/s and $V_{\rm \sigma, outer-disk}$ = 15 $\pm$ 2 km/s). The low gas velocity dispersion is similar to local spiral galaxies and is consistent with the classic density wave theory where spiral arms form in dynamically cold and thin disks. We speculate that A1689B11 belongs to a population of rare spiral galaxies at $z\gtrsim2$ that mark the formation epoch of thin disks. Future observations with JWST will greatly increase the sample of these rare galaxies and unveil the earliest onset of spiral arms.

We present an analysis of the global stellar populations of galaxies in the SAMI Galaxy Survey. Our sample consists of 1319 galaxies spanning four orders of magnitude in stellar mass and includes all morphologies and environments. We derive luminosity-weighted, single stellar population equivalent stellar ages, metallicities and alpha enhancements from spectra integrated within one effective radius apertures. Variations in galaxy size explain the majority of the scatter in the age--mass and metallicity--mass relations. Stellar populations vary systematically in the plane of galaxy size and stellar mass, such that galaxies with high stellar surface mass density are older, more metal-rich and alpha-enhanced than less dense galaxies. Galaxies with high surface mass densities have a very narrow range of metallicities, however, at fixed mass, the spread in metallicity increases substantially with increasing galaxy size (decreasing density). We identify residual correlations with morphology and environment. At fixed mass and size, galaxies with late-type morphologies, small bulges and low Sersic n are younger than early-type, high n, high bulge-to-total galaxies. Age and metallicity both show small residual correlations with environment; at fixed mass and size, galaxies in denser environments or more massive halos are older and somewhat more metal rich than those in less dense environments. We connect these trends to evolutionary tracks within the size--mass plane.

The spatial distribution of oxygen in the interstellar medium of galaxies is the key to understanding how efficiently metals that are synthesized in massive stars can be redistributed across a galaxy. We present here a case study in the nearby spiral galaxy NGC1365 using 3D optical data obtained in the TYPHOON Program. We find systematic azimuthal variations of the HII region oxygen abundance imprinted on a negative radial gradient. The 0.2 dex azimuthal variations occur over a wide radial range of 0.3 to 0.7 R25 and peak at the two spiral arms in NGC1365. We show that the azimuthal variations can be explained by two physical processes: gas undergoes localized, sub-kpc scale self-enrichment when orbiting in the inter-arm region, and experiences efficient, kpc scale mixing-induced dilution when spiral density waves pass through. We construct a simple chemical evolution model to quantitatively test this picture and find that our toy model can reproduce the observations. This result suggests that the observed abundance variations in NGC1365 are a snapshot of the dynamical local enrichment of oxygen modulated by spiral-driven, periodic mixing and dilution.

We present the second and final data release of the Siding Spring Southern Seyfert Spectroscopic Snapshot Survey (S7). Data are presented for 63 new galaxies not included in the first data release, and we provide 2D emission-line fitting products for the full S7 sample of 131 galaxies. The S7 uses the WiFeS instrument on the ANU 2.3m telescope to obtain spectra with a spectral resolution of R = 7000 in the red (540 - 700 nm) and R = 3000 in the blue (350 - 570 nm), over an integral field of 25x38 arcsec^2 with 1x1 arcsec^2 spatial pixels. The S7 contains both the largest sample of active galaxies and the highest spectral resolution of any comparable integral field survey to date. The emission-line fitting products include line fluxes, velocities and velocity dispersions across the WiFeS field of view, and an artificial neural network has been used to determine the optimal number of Gaussian kinematic components for emission-lines in each spaxel. Broad Balmer lines are subtracted from the spectra of nuclear spatial pixels in Seyfert 1 galaxies before fitting the narrow lines. We bin nuclear spectra and measure reddening-corrected nuclear fluxes of strong narrow lines for each galaxy. The nuclear spectra are classified on optical diagnostic diagrams, where the strength of the coronal line [FeVII]6087 is shown to be correlated with [OIII]/Hbeta. Maps revealing gas excitation and kinematics are included for the entire sample, and we provide notes on the newly observed objects.

We present the first major release of data from the SAMI Galaxy Survey. This data release focuses on the emission-line physics of galaxies. Data Release One includes data for 772 galaxies, about 20% of the full survey. Galaxies included have the redshift range 0.004 < z < 0.092, a large mass range (7.6 < log(Mstellar/M$_\odot$) < 11.6), and star-formation rates of 10^-4 to 10^1\ M$_\odot$/yr. For each galaxy, we include two spectral cubes and a set of spatially resolved 2D maps: single- and multi-component emission-line fits (with dust extinction corrections for strong lines), local dust extinction and star-formation rate. Calibration of the fibre throughputs, fluxes and differential-atmospheric-refraction has been improved over the Early Data Release. The data have average spatial resolution of 2.16 arcsec (FWHM) over the 15~arcsec diameter field of view and spectral (kinematic) resolution R=4263 (sigma=30km/s) around Halpha. The relative flux calibration is better than 5\% and absolute flux calibration better than $\pm0.22$~mag, with the latter estimate limited by galaxy photometry. The data are presented online through the Australian Astronomical Observatory's Data Central.

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.

In this paper, we compare the radio emission at 1.4 GHz with optical outflow signatures of edge-on galaxies. We report observations of six edge-on star-forming galaxies in the Sydney-AAO Multi-object Integral-field spectrograph (SAMI) Galaxy Survey with 1.4 GHz luminosities $> 1\times10^{21}$ W Hz$^{-1}$. Extended minor axis optical emission is detected with enhanced \nii/H$\alpha$ line ratios and velocity dispersions consistent with galactic winds in three of six galaxies. These galaxies may host outflows driven by a combination of thermal and cosmic ray processes. We find that galaxies with the strongest wind signatures have extended radio morphologies. Our results form a baseline for understanding the driving mechanisms of galactic winds.

The most rapidly evolving regions of galaxies often display complex optical spectra with emission lines excited by massive stars, shocks and accretion onto supermassive black holes. Standard calibrations (such as for the star formation rate) cannot be applied to such mixed spectra. In this paper we isolate the contributions of star formation, shock excitation and active galactic nucleus (AGN) activity to the emission line luminosities of individual spatially resolved regions across the central 3 $\times$ 3 kpc$^2$ region of the active barred spiral galaxy NGC$\sim$613. The star formation rate and AGN luminosity calculated from the decomposed emission line maps are in close agreement with independent estimates from data at other wavelengths. The star formation component traces the B-band stellar continuum emission, and the AGN component forms an ionization cone which is aligned with the nuclear radio jet. The optical line emission associated with shock excitation is cospatial with strong $H_2$ and [Fe II] emission and with regions of high ionized gas velocity dispersion ($\sigma > 100$ km s$^{-1}$). The shock component also traces the outer boundary of the AGN ionization cone and may therefore be produced by outflowing material interacting with the surrounding interstellar medium. Our decomposition method makes it possible to determine the properties of star formation, shock excitation and AGN activity from optical spectra, without contamination from other ionization mechanisms.

The spatial distribution of dust in galaxies affects the global attenuation, and hence inferred properties, of galaxies. We trace the spatial distribution of dust in five fields (at 0.6-0.9 kpc scale) of M31 by comparing optical attenuation with the total dust mass distribution. We measure the attenuation from the Balmer decrement using Integral Field Spectroscopy and the dust mass from Herschel far-IR observations. Our results show that M31's dust attenuation closely follows a foreground screen model, contrary to what was previously found in other nearby galaxies. By smoothing the M31 data we find that spatial resolution is not the cause for this difference. Based on the emission line ratios and two simple models, we conclude that previous models of dust/gas geometry need to include a weakly or non-attenuated diffuse ionized gas (DIG) component. Due to the variation of dust and DIG scale heights with galactic radius, we conclude that different locations in galaxies will have different vertical distributions of gas and dust and therefore different measured attenuation. The difference between our result in M31 with that found in other nearby galaxies can be explained by our fields in M31 lying at larger galactic radii than the previous studies that focused on the centres of galaxies.

We present results on the dust attenuation of galaxies at redshift $\sim3-6$ by studying the relationship between the UV spectral slope ($\beta_{\rm UV}$) and the infrared excess (IRX; $L_{\rm IR}$/$L_{\rm UV}$) using ALMA far-infrared continuum observations. Our study is based on a sample of 67 massive, star-forming galaxies with a median mass of $M_{\ast}\sim 10^{10.7}\,M_{\rm \odot}$ spanning a redshift range $z=2.6-3.7$ (median $z=3.2$) that were observed with ALMA at $\lambda_{rest}=300\,{\rm \mu m}$. Both the individual ALMA detections (41 sources) and stacks including all galaxies show the IRX-$\beta_{\rm UV}$ relationship at $z\sim3$ is mostly consistent with that of local starburst galaxies on average. However, we find evidence for a large dispersion around the mean relationship by up to $\pm0.5$ dex. Nevertheless, the locally calibrated dust correction factors based on the IRX-$\beta_{\rm UV}$ relation are on average applicable to main-sequence $z\sim3$ galaxies. This does not appear to be the case at even higher redshifts, however. Using public ALMA observations of $z\sim4-6$ galaxies we find evidence for a significant evolution in the IRX-$\beta_{\rm UV}$ and the IRX-$M_{\ast}$ relations beyond $z\sim3$ toward lower IRX values. We discuss several caveats that could affect these results, including the assumed dust temperature. ALMA observations of larger $z>3$ galaxy samples will be required to confirm this intriguing redshift evolution.

The ratio of the [CII] 158$\,\mu$m emission line over the total infrared emission (TIR) is often used as a proxy for the photoelectric (PE) heating efficiency ($\epsilon_{\rm PE}$) of the far-ultraviolet (FUV) photons absorbed by dust in the interstellar medium. In the nearby galaxy M31, we measure a strong radial variation of [CII]/TIR that we rule out as being due to an intrinsic variation in $\epsilon_{\rm PE}$. [CII]/TIR fails as a proxy for $\epsilon_{\rm PE}$, because the TIR measures all dust heating, not just the contribution from FUV photons capable of ejecting electrons from dust grains. Using extensive multiwavelength coverage from the FUV to far-infrared (FIR), we infer the attenuated FUV emission ($\rm UV_{att}$), and the total attenuated flux ($\rm TOT_{att}$). We find [CII]/TIR to be strongly correlated with $\rm UV_{att}$/$\rm TOT_{att}$, indicating that, in M31 at least, one of the dominant drivers for [CII]/TIR variation is the relative hardness of the absorbed stellar radiation field. We define $\rm{ \epsilon_{PE}^{UV}}$, [CII]/$\rm{ UV_{att}}$ which should be more closely related to the actual PE efficiency, which we find to be essentially constant ($1.85 \pm 0.8 \%$) in all explored fields in M31. This suggests that part of the observed variation of [CII]/TIR in other galaxies is likely due to a change in the relative hardness of the absorbed stellar radiation field, caused by a combination of variations in the stellar population, dust opacity and galaxy metallicity, although PE efficiency may also vary across a wider range of environments.

Recently, large samples of visually classified early-type galaxies (ETGs) containing dust have been identified using space-based infrared observations with the Herschel Space Telescope. The presence of large quantities of dust in massive ETGs is peculiar as X-ray halos of these galaxies are expected to destroy dust in 10 Myr (or less). This has sparked a debate regarding the origin of the dust: is it internally produced by asymptotic giant branch (AGB) stars, or is it accreted externally through mergers? We examine the 2D stellar and ionised gas kinematics of dusty ETGs using IFS observations from the SAMI galaxy survey, and integrated star-formation rates, stellar masses, and dust masses from the GAMA survey. Only 8% (4/49) of visually-classified ETGs are kinematically consistent with being dispersion-supported systems. These "dispersion-dominated galaxies" exhibit discrepancies between stellar and ionised gas kinematics, either offsets in the kinematic position angle or large differences in the rotational velocity, and are outliers in star-formation rate at a fixed dust mass compared to normal star-forming galaxies. These properties are suggestive of recent merger activity. The remaining 90% of dusty ETGs have low velocity dispersions and/or large circular velocities, typical of "rotation-dominated galaxies". These results, along with the general evidence of published works on X-ray emission in ETGs, suggest that they are unlikely to host hot, X-ray gas consistent with their low stellar mass when compared to dispersion-dominated galaxies. This means dust will be long lived and thus these galaxies do not require external scenarios for the origin of their dust content.

Stars form in cold molecular clouds. However, molecular gas is difficult to observe because the most abundant molecule (H2) lacks a permanent dipole moment. Rotational transitions of CO are often used as a tracer of H2, but CO is much less abundant and the conversion from CO intensity to H2 mass is often highly uncertain. Here we present a new method for estimating the column density of cold molecular gas (Sigma_gas) using optical spectroscopy. We utilise the spatially resolved H-alpha maps of flux and velocity dispersion from the Sydney-AAO Multi-object Integral-field spectrograph (SAMI) Galaxy Survey. We derive maps of Sigma_gas by inverting the multi-freefall star formation relation, which connects the star formation rate surface density (Sigma_SFR) with Sigma_gas and the turbulent Mach number (Mach). Based on the measured range of Sigma_SFR = 0.005-1.5 M_sol/yr/kpc^2 and Mach = 18-130, we predict Sigma_gas = 7-200 M_sol/pc^2 in the star-forming regions of our sample of 260 SAMI galaxies. These values are close to previously measured Sigma_gas obtained directly with unresolved CO observations of similar galaxies at low redshift. We classify each galaxy in our sample as 'Star-forming' (219) or 'Composite/AGN/Shock' (41), and find that in Composite/AGN/Shock galaxies the average Sigma_SFR, Mach, and Sigma_gas are enhanced by factors of 2.0, 1.6, and 1.3, respectively, compared to Star-forming galaxies. We compare our predictions of Sigma_gas with those obtained by inverting the Kennicutt-Schmidt relation and find that our new method is a factor of two more accurate in predicting Sigma_gas, with an average deviation of 32% from the actual Sigma_gas.

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 a new basis for scaling abundances with total metallicity in nebular photoionisation models, based on extensive Milky Way stellar abundance data, to replace the uniform scaling normally used in the analysis of HII regions. Our goal is to provide a single scaling method and local abundance reference standard for use in nebular modelling and its key inputs, the stellar atmosphere and evolutionary track models. We introduce a parametric enrichment factor, zeta, to describe how atomic abundances scale with total abundance, and which allows for a simple conversion between scales based on different reference elements (usually oxygen or iron) . The models and parametric description provide a more physically realistic approach than simple uniform abundance scaling. With appropriate parameters, the methods described here may be applied to HII regions in the Milky Way, large and dwarf galaxies in the local universe, Active Galactic Nuclei (AGNs), and to star forming regions at high redshift.

Distance uncertainties plague our understanding of the physical scales relevant to the physics of star formation in extragalactic studies. The planetary nebulae luminosity function (PNLF) is one of very few techniques that can provide distance estimates to within ~10%, however it requires a planetary nebula (PN) sample that is uncontaminated by other ionizing sources. We employ optical IFU spectroscopy using MUSE on the VLT to measure [OIII] line fluxes for sources unresolved on 50 pc scales within the central star-forming galaxy disk of NGC 628. We use diagnostic line ratios to identify 62 PNe, 30 supernova remnants and 87 HII regions within our fields. Using the 36 brightest PNe we determine a new PNLF distance modulus of 29.91^{+0.08}_{-0.13} mag (9.59^{+0.35}_{-0.57} Mpc), in good agreement with literature values but significantly larger than the previously reported PNLF distance. We are able to explain the discrepancy and recover the previous result when we reintroduce SNR contaminants to our sample. This demonstrates the power of full spectral information over narrowband imaging in isolating PNe. Given our limited spatial coverage within the galaxy, we show that this technique can be used to refine distance estimates even when IFU observations cover only a fraction of a galaxy disk.

Recent cosmological hydrodynamical simulations suggest that integral field spectroscopy can connect the high-order stellar kinematic moments h3 (~skewness) and h4 (~kurtosis) in galaxies to their cosmological assembly history. Here, we assess these results by measuring the stellar kinematics on a sample of 315 galaxies, without a morphological selection, using 2D integral field data from the SAMI Galaxy Survey. A proxy for the spin parameter ($\lambda_{R_e}$) and ellipticity ($\epsilon_e$) are used to separate fast and slow rotators; there exists a good correspondence to regular and non-regular rotators, respectively, as also seen in earlier studies. We confirm that regular rotators show a strong h3 versus $V/\sigma$ anti-correlation, whereas quasi-regular and non-regular rotators show a more vertical relation in h3 and $V/\sigma$. Motivated by recent cosmological simulations, we develop an alternative approach to kinematically classify galaxies from their individual h3 versus $V/\sigma$ signatures. We identify five classes of high-order stellar kinematic signatures using Gaussian mixture models. Class 1 corresponds to slow rotators, whereas Classes 2-5 correspond to fast rotators. We find that galaxies with similar $\lambda_{R_e}-\epsilon_e$ values can show distinctly different h3-$V/\sigma$ signatures. Class 5 objects are previously unidentified fast rotators that show a weak h3 versus $V/\sigma$ anti-correlation. These objects are predicted to be disk-less galaxies formed by gas-poor mergers. From morphological examination, however, there is evidence for large stellar disks. Instead, Class 5 objects are more likely disturbed galaxies, have counter-rotating bulges, or bars in edge-on galaxies. Finally, we interpret the strong anti-correlation in h3 versus $V/\sigma$ as evidence for disks in most fast rotators, suggesting a dearth of gas-poor mergers among fast rotators.

We present a simplified model of Active Galactic Nucleus (AGN) continuum emission designed for photoionization modeling. The new model {\sc oxaf} reproduces the diversity of spectral shapes that arise in physically-based models. We identify and explain degeneracies in the effects of AGN parameters on model spectral shapes, with a focus on the complete degeneracy between the black hole mass and AGN luminosity. Our re-parametrized model {\sc oxaf} removes these degeneracies and accepts three parameters which directly describe the output spectral shape: the energy of the peak of the accretion disk emission $E_\mathrm{peak}$, the photon power-law index of the non-thermal emission $\Gamma$, and the proportion of the total flux which is emitted in the non-thermal component $p_\mathrm{NT}$. The parameter $E_\mathrm{peak}$ is presented as a function of the black hole mass, AGN luminosity, and `coronal radius' of the {\sc optxagnf} model upon which {\sc oxaf} is based. We show that the soft X-ray excess does not significantly affect photoionization modeling predictions of strong emission lines in Seyfert narrow-line regions. Despite its simplicity, {\sc oxaf} accounts for opacity effects where the accretion disk is ionized because it inherits the `color correction' of {\sc optxagnf}. We use a grid of {\sc mappings} photoionization models with {\sc oxaf} ionizing spectra to demonstrate how predicted emission-line ratios on standard optical diagnostic diagrams are sensitive to each of the three {\sc oxaf} parameters. The {\sc oxaf} code is publicly available in the Astrophysics Source Code Library.

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.

Star-forming galaxies at $z > 1$ exhibit significantly different properties to local galaxies of equivalent stellar mass. Not only are high-redshift star-forming galaxies characterized by higher star formation rates and gas fractions than their local counterparts, they also appear to host star-forming regions with significantly different physical conditions, including greater electron densities. To understand what physical mechanisms are responsible for the observed evolution of star-forming conditions we have assembled the largest sample of star-forming galaxies at $z\sim 1.5$ with emission-line measurements of the $\mathrm{[OII]} \lambda \lambda 3726,3729$ doublet. By comparing our $z\sim 1.5$ sample to local galaxy samples with equivalent distributions of stellar mass, star formation rate and specific star formation rate we investigate the proposed evolution in electron density and its dependence on global properties. We measure an average electron density of $114_{-27}^{+28} \, \mathrm {cm}^{-3} $ for our $z\sim 1.5$ sample, a factor of five greater than the typical electron density of local star-forming galaxies. However, we find no offset between the typical electron densities of local and high-redshift galaxies with equivalent star-formation rates. Our work indicates that the average electron density of a sample is highly sensitive to the star formation rates, implying that the previously observed evolution is mainly the result of selection effects.

The observed evolution of the gas fraction and its associated depletion time in main sequence (MS) galaxies provides insights on how star formation proceeds over cosmic time. We report ALMA detections of the rest-frame $\sim$300$\mu$m continuum observed at 240 GHz for 45 massive ($\rm \langle log(M_{\star}(M_{\odot}))\rangle=10.7$), normal star forming ($\rm \langle log(sSFR(yr^{-1}))\rangle=-8.6$), i.e. MS, galaxies at $\rm z\approx3.2$ in the COSMOS field. From an empirical calibration between cold neutral, i.e. molecular and atomic, gas mass $\rm M_{gas}$ and monochromatic (rest-frame) infrared luminosity, the gas mass for this sample is derived. Combined with stellar mass $\rm M_{\star}$ and star formation rate (SFR) estimates (from {\sc MagPhys} fits) we obtain a median gas fraction of $\rm \mu_{gas}=M_{gas}/M_{\star}=1.65_{-0.19}^{+0.18}$ and a median gas depletion time $\rm t_{depl.}(Gyr)=M_{gas}/SFR=0.68_{-0.08}^{+0.07}$; correction for the location on the MS will only slightly change the values. The reported uncertainties are the $\rm 1\sigma$ error on the median. Our results are fully consistent with the expected flattening of the redshift evolution from the 2-SFM (2 star formation mode) framework that empirically prescribes the evolution assuming a universal, log-linear relation between SFR and gas mass coupled to the redshift evolution of the specific star formation rate (sSFR) of main sequence galaxies. While $\rm t_{dep.}$ shows only a mild dependence on location within the MS, a clear trend of increasing $\rm \mu_{gas}$ across the MS is observed (as known from previous studies). Further we comment on trends within the MS and (in)consistencies with other studies.

Stellar population synthesis models can now reproduce the photometry of old stellar systems (age $>$ 2 Gyr) in the near-infrared (NIR) bands at 3.4 and 4.6$\mu$m (WISE W1 $\&$ W2 or IRAC 1 $\&$ 2). In this paper we derive stellar mass-to-light ratios for these and optical bands, and confirm that the NIR M/L shows dramatically reduced sensitivity to both age and metallicity compared to optical bands, and further, that this behavior leads to significantly more robust stellar masses for quiescent galaxies with [Fe/H] > -0.5 regardless of star formation history (SFH). We then use realistic early-type galaxy SFHs and metallicity distributions from the EAGLE simulations of galaxy formation to investigate two methods to determine the appropriate M/L for a galaxy: 1) We show that the uncertainties introduced by an unknown SFH can be largely removed using a spectroscopically inferred luminosity-weighted age and metallicity for the population to select the appropriate single stellar population (SSP) equivalent M/L. Using this method, the maximum systematic error due to SFH on the M/L of an early-type galaxy is $<$ 4$\%$ at 3.4 $\mu$m and typical uncertainties due to errors in the age and metallicity create scatter of $\lesssim$13$\%$. The equivalent values for optical bands are more than 2-3 times greater, even before considering uncertainties associated with internal dust extinction. 2) We demonstrate that if the EAGLE SFHs and metallicities accurately reproduce the true properties of early-type galaxies, the use of an iterative approach to select a mass dependent M/L can provide even more accurate stellar masses for early-type galaxies, with typical uncertainties $<$ 9$\%$.

We present LZIFU (LaZy-IFU), an IDL toolkit for fitting multiple emission lines simultaneously in integral field spectroscopy (IFS) data. LZIFU is useful for the investigation of the dynamical, physical and chemical properties of gas in galaxies. LZIFU has already been applied to many world-class IFS instruments and large IFS surveys, including the Wide Field Spectrograph, the new Multi Unit Spectroscopic Explorer (MUSE), the Calar Alto Legacy Integral Field Area (CALIFA) survey, the Sydney-Australian-astronomical-observatory Multi-object Integral-field spectrograph (SAMI) Galaxy Survey. Here we describe in detail the structure of the toolkit, and how the line fluxes and flux uncertainties are determined, including the possibility of having multiple distinct kinematic components. We quantify the performance of LZIFU, demonstrating its accuracy and robustness. We also show examples of applying LZIFU to CALIFA and SAMI data to construct emission line and kinematic maps, and investigate complex, skewed line profiles presented in IFS data. The code is made available to the astronomy community through github. LZIFU will be further developed over time to other IFS instruments, and to provide even more accurate line and uncertainty estimates.

The optical spectra of Seyfert galaxies are often dominated by emission lines excited by both star formation and AGN activity. Standard calibrations (such as for the star formation rate) are not applicable to such composite (mixed) spectra. In this paper, we describe how integral field data can be used to spectrally and spatially separate emission associated with star formation from emission associated with accretion onto an active galactic nucleus (AGN). We demonstrate our method using integral field data for two AGN host galaxies (NGC 5728 and NGC 7679) from the Siding Spring Southern Seyfert Spectroscopic Snapshot Survey (S7). The spectra of NGC 5728 and NGC 7679 form clear sequences of AGN fraction on standard emission line ratio diagnostic diagrams. We show that the emission line luminosities of the majority (> 85 per cent) of spectra along each AGN fraction sequence can be reproduced by linear superpositions of the emission line luminosities of one AGN dominated spectrum and one star formation dominated spectrum. We separate the Halpha, Hbeta, [N II]$\lambda$6583, [S II]$\lambda \lambda$6716, 6731, [O III]$\lambda$5007 and [O II]$\lambda \lambda$3726, 3729 luminosities of every spaxel into contributions from star formation and AGN activity. The decomposed emission line images are used to derive the star formation rates and AGN bolometric luminosities for NGC 5728 and NGC 7679. Our calculated values are mostly consistent with independent estimates from data at other wavelengths. The recovered star forming and AGN components also have distinct spatial distributions which trace structures seen in high resolution imaging of the galaxies, providing independent confirmation that our decomposition has been successful.

We present The Machine, an artificial neural network (ANN) capable of differentiating between the numbers of Gaussian components needed to describe the emission lines of Integral Field Spectroscopic (IFS) observations. Here we show the preliminary results of the S7 first data release (Siding Spring Southern Seyfert Spectro- scopic Snapshot Survey, Dopita et al. 2015) and SAMI Galaxy Survey (Sydney-AAO Multi-object Integral Field Unit, Croom et al. 2012) to classify whether the emission lines in each spatial pixel are composed of 1, 2, or 3 different Gaussian components. Previously this classification has been done by individual people, taking an hour per galaxy. This time investment is no longer feasible with the large spectroscopic surveys coming online.

Interarm star formation contributes significantly to a galaxy's star formation budget, and provides an opportunity to study stellar birthplaces unperturbed by spiral arm dynamics. Using optical integral field spectroscopy of the nearby galaxy NGC 628 with VLT/MUSE, we construct Halpha maps including detailed corrections for dust extinction and stellar absorption to identify 391 HII regions at 35pc resolution over 12 kpc^2. Using tracers sensitive to the underlying gravitational potential, we associate HII regions with either arm (271) or interarm (120) environments. Using our full spectral coverage of each region, we find that most HII region physical properties (luminosity, size, metallicity, ionization parameter) are independent of environment. We calculate the fraction of Halpha luminosity due to the diffuse ionized gas (DIG) background contaminating each HII region, and find the DIG surface brightness to be higher within HII regions compared to the surroundings, and slightly higher within arm HII regions. Use of the temperature sensitive [SII]/Halpha line ratio map instead of the Halpha surface brightness to identify HII region boundaries does not change this result. Using the dust attenuation as a tracer of the gas, we find depletion times consistent with previous work (2 x 10^9 yr) with no differences between the arm and interarm, however this is very sensitive to the DIG correction. Unlike molecular clouds, which can be dynamically affected by the galactic environment, we see fairly consistent HII region properties in both arm and interarm environments. This suggests either a difference in arm star formation and feedback, or a decoupling of dense star forming clumps from the more extended surrounding molecular gas.

We report here on the discovery of faint extended wings of H\alpha\ observed out to an apparent velocity of ~ 7600 km/s in the Orion Nebula (M42) and in five HII regions in the Large and the Small Magellanic Clouds. We show that, these wings are caused by Raman scattering of both the O I and Si II resonance lines and stellar continuum UV photons with H I followed by radiative decay to the H I n=2 level. The broad wings also seen in H\beta\ and in H\gamma\ result from Raman scattering of the UV continuum in the H I n=4 and n=5 levels respectively.The Raman scattering fluorescence is correlated with the intensity of the narrow permitted lines of O I and Si II. In the case of Si II, this is explained by radiative pumping of the same 1023.7\AA\ resonance line involved in the Raman scattering by the Ly\beta\ radiation field. The subsequent radiative cascade produces enhanced Si II 5978.9, 6347.1$ and 6371.4\AA\ permitted transitions. Finally we show that in O I, radiative pumping of the 1025.76\AA\ resonance line by the Lyman series radiation field is also the cause of the enhancement in the permitted lines of this species lying near H\alpha\ in wavelength, but here the process is a little more complex. We argue that all these processes are active in the zone of the HII region near the ionisation front

The [NII] 122 and 205 \mu m transitions are powerful tracers of the ionized gas in the interstellar medium. By combining data from 21 galaxies selected from the Herschel KINGFISH and Beyond the Peak surveys, we have compiled 141 spatially resolved regions with a typical size of ~1 kiloparsec, with observations of both [NII] far-infrared lines. We measure [NII] 122/205 line ratios in the ~0.6-6 range, which corresponds to electron gas densities $n_e$~1-300 cm$^{-3}$, with a median value of $n_e$=30 cm$^{-3}$. Variations in the electron density within individual galaxies can be as a high as a factor of ~50, frequently with strong radial gradients. We find that $n_e$ increases as a function of infrared color, dust-weighted mean starlight intensity, and star formation rate surface density ($\Sigma_{SFR}$). As the intensity of the [NII] transitions is related to the ionizing photon flux, we investigate their reliability as tracers of the star formation rate (SFR). We derive relations between the [NII] emission and SFR in the low-density limit and in the case of a log-normal distribution of densities. The scatter in the correlation between [NII] surface brightness and $\Sigma_{SFR}$ can be understood as a property of the $n_e$ distribution. For regions with $n_e$ close to or higher than the [NII] line critical densities, the low-density limit [NII]-based SFR calibration systematically underestimates the SFR since [NII] emission is collisionally quenched. Finally, we investigate the relation between [NII] emission, SFR, and $n_e$ by comparing our observations to predictions from the MAPPINGS-III code.

We investigate the relative significance of radiation pressure and gas pressure in the extended narrow line regions (ENLRs) of four Seyfert galaxies from the integral field Siding Spring Southern Seyfert Spectroscopic Snapshot Survey (S7). We demonstrate that there exist two distinct types of starburst-AGN mixing curves on standard emission line diagnostic diagrams which reflect the balance between gas pressure and radiation pressure in the ENLR. In two of the galaxies the ENLR is radiation pressure dominated throughout and the ionization parameter remains constant (log U ~ 0). In the other two galaxies radiation pressure is initially important, but gas pressure becomes dominant as the ionization parameter in the ENLR decreases from log U ~ 0 to -3.4 <= log U <= -3.2. Where radiation pressure is dominant, the AGN regulates the density of the interstellar medium on kpc scales and may therefore have a direct impact on star formation activity and/or the incidence of outflows in the host galaxy to scales far beyond the zone of influence of the black hole. We find that both radiation pressure dominated and gas pressure dominated ENLRs are dynamically active with evidence for outflows, indicating that radiation pressure may be an important source of AGN feedback even when it is not dominant over the entire ENLR.

We map the dust distribution in the central 180" (~680 pc) region of the M31 bulge, based on HST/WFC3 and ACS observations in ten bands from near-ultraviolet (2700 A) to near-infrared (1.5 micron). This large wavelength coverage gives us great leverage to detect not only dense dusty clumps, but also diffuse dusty molecular gas. We fit a pixel-by-pixel spectral energy distributions to construct a high-dynamic-range extinction map with unparalleled angular resolution (~0.5" , i.e., ~2 pc) and sensitivity (the extinction uncertainty, \delta A_V~0.05). In particular, the data allow to directly fit the fractions of starlight obscured by individual dusty clumps, and hence their radial distances in the bulge. Most of these clumps seem to be located in a thin plane, which is tilted with respect to the M31 disk and appears face-on. We convert the extinction map into a dust mass surface density map and compare it with that derived from the dust emission as observed by Herschel . The dust masses in these two maps are consistent with each other, except in the low-extinction regions, where the mass inferred from the extinction tends to be underestimated. Further, we use simulations to show that our method can be used to measure the masses of dusty clumps in Virgo cluster early-type galaxies to an accuracy within a factor of ~2.

We investigate a sample of 40 local, main-sequence, edge-on disc galaxies using integral field spectroscopy with the Sydney-AAO Multi-object Integral field spectrograph (SAMI) Galaxy Survey to understand the link between properties of the extraplanar gas and their host galaxies. The kinematics properties of the extraplanar gas, including velocity asymmetries and increased dispersion, are used to differentiate galaxies hosting large-scale galactic winds from those dominated by the extended diffuse ionized gas. We find rather that a spectrum of diffuse gas-dominated to wind dominated galaxies exist. The wind-dominated galaxies span a wide range of star formation rates ($-1 \lesssim \log({\rm SFR/M_{\odot} yr^{-1}}) \lesssim 0.5$) across the whole stellar mass range of the sample ($8.5 \lesssim \log({\rm M_{*}/M_{\odot}}) \lesssim 11$). The wind galaxies also span a wide range in SFR surface densities ($10^{-3} \textrm{--} 10^{-1.5}\rm~M_{\odot} ~yr^{-1}~kpc^{-2}$) that is much lower than the canonical threshold of $\rm0.1~M_{\odot} ~yr^{-1}~kpc^{-2}$. The wind galaxies on average have higher SFR surface densities and higher $\rm H\delta_A$ values than those without strong wind signatures. The enhanced $\rm H\delta_A$ indicates that bursts of star formation in the recent past are necessary for driving large-scale galactic winds. We demonstrate with Sloan Digital Sky Survey data that galaxies with high SFR surface density have experienced bursts of star formation in the recent past. Our results imply that the galactic winds revealed in our study are indeed driven by bursts of star formation, and thus probing star formation in the time domain is crucial for finding and understanding galactic winds.

(abridged) The majority of star-forming galaxies follow a simple empirical correlation in the star formation rate (SFR) versus stellar mass ($M_*$) plane, usually referred to as the star formation Main Sequence (MS). Here we combine a set of hydro-dynamical simulations of interacting galactic disks with state-of-the-art radiative transfer codes to analyze how the evolution of mergers is reflected upon the properties of the MS. We present \textsc{Chiburst}, a Markov Chain Monte Carlo (MCMC) Spectral Energy Distribution (SED) code that fits the multi-wavelength, broad-band photometry of galaxies and derives stellar masses, star formation rates, and geometrical properties of the dust distribution. We apply this tool to the SEDs of simulated mergers and compare the derived results with the reference output from the simulations. Our results indicate that changes in the SEDs of mergers as they approach coalescence and depart from the MS are related to an evolution of dust geometry in scales larger than a few hundred parsecs. This is reflected in a correlation between the specific star formation rate (sSFR), and the compactness parameter $\mathcal{C}$, that parametrizes this geometry and hence the evolution of dust temperature ($T_{\rm{dust}}$) with time. As mergers approach coalescence, they depart from the MS and increase their compactness, which implies that moderate outliers of the MS are consistent with late-type mergers. By further applying our method to real observations of Luminous Infrared Galaxies (LIRGs), we show that the merger scenario is unable to explain these extreme outliers of the MS. Only by significantly increasing the gas fraction in the simulations are we able to reproduce the SEDs of LIRGs.

The mid-infrared is an optimal window to trace stellar mass in nearby galaxies and the 3.6$\mu m$ IRAC band has been exploited to this effect, but such mass estimates can be biased by dust emission. We present our pipeline to reveal the old stellar flux at 3.6$\mu m$ and obtain stellar mass maps for more than 1600 galaxies available from the Spitzer Survey of Stellar Structure in Galaxies (S$^{4}$G). This survey consists of images in two infrared bands (3.6 and 4.5$\mu m$), and we use the Independent Component Analysis (ICA) method presented in Meidt et al. (2012) to separate the dominant light from old stars and the dust emission that can significantly contribute to the observed 3.6$\mu m$ flux. We exclude from our ICA analysis galaxies with low signal-to-noise ratio (S/N < 10) and those with original [3.6]-[4.5] colors compatible with an old stellar population, indicative of little dust emission (mostly early Hubble types, which can directly provide good mass maps). For the remaining 1251 galaxies to which ICA was successfully applied, we find that as much as 10-30% of the total light at 3.6$\mu m$ typically originates from dust, and locally it can reach even higher values. This contamination fraction shows a correlation with specific star formation rates, confirming that the dust emission that we detect is related to star formation. Additionally, we have used our large sample of mass estimates to calibrate a relationship of effective mass-to-light ratio ($M/L$) as a function of observed [3.6]-[4.5] color: $\log(M/L)=-0.339 (\pm 0.057) \times ([3.6]-[4.5]) -0.336 (\pm 0.002)$. Our final pipeline products have been made public through IRSA, providing the astronomical community with an unprecedentedly large set of stellar mass maps ready to use for scientific applications.

Galaxy surveys targeting emission lines are characterising the evolution of star-forming galaxies, but there is still little theoretical progress in modelling their physical properties. We predict nebular emission from star-forming galaxies within a cosmological galaxy formation model. Emission lines are computed by combining the semi-analytical model \sag\ with the photoionisation code \mapp. We characterise the interstellar medium (ISM) of galaxies by relating the ionisation parameter of gas in galaxies to their cold gas metallicity, obtaining a reasonable agreement with the observed \ha, \oii, \oiii\ luminosity functions, and the the BPT diagram for local star-forming galaxies. The average ionisation parameter is found to increase towards low star-formation rates and high redshifts, consistent with recent observational results. The predicted link between different emission lines and their associated star-formation rates is studied by presenting scaling relations to relate them. Our model predicts that emission line galaxies have modest clustering bias, and thus reside in dark matter haloes of masses below $M_{\rm halo} \lesssim 10^{12} {[\rm h^{-1} M_{\odot}]}$. Finally, we exploit our modelling technique to predict galaxy number counts up to $z\sim 10$ by targeting far-infrared (FIR) emission lines detectable with submillimetre facilities

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.

The current Lambda CDM cosmological model predicts that galaxy evolution proceeds more slowly in lower density environments, suggesting that voids are a prime location to search for relatively pristine galaxies that are representative of the building blocks of early massive galaxies. To test the assumption that void galaxies are more pristine, we compare the evolutionary properties of a sample of dwarf galaxies selected specifically to lie in voids with a sample of similar isolated dwarf galaxies in average density environments. We measure gas-phase oxygen abundances and gas fractions for eight dwarf galaxies (M_r > -16.2), carefully selected to reside within the lowest density environments of seven voids, and apply the same calibrations to existing samples of isolated dwarf galaxies. We find no significant difference between these void dwarf galaxies and the isolated dwarf galaxies, suggesting that dwarf galaxy chemical evolution proceeds independent of the large-scale environment. While this sample is too small to draw strong conclusions, it suggests that external gas accretion is playing a limited role in the chemical evolution of these systems, and that this evolution is instead dominated mainly by the internal secular processes that are linking the simultaneous growth and enrichment of these galaxies.

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] 158 micron line is one of the strongest emission lines observed in star-forming galaxies, and has been empirically measured to correlate with the star formation rate (SFR) globally and on ~kpc scales. However, due to the multi-phase origins of [CII], one might expect this relation to break down at small scales. We investigate the origins of [CII] emission by examining high spatial resolution observations of [CII] in M31, with the Survey of Lines in M31 (SLIM). We present five ~700x700 pc (3"x3") Fields mapping the [CII] emission, Halpha emission, combined with ancillary infrared (IR) data. We spatially separate star-forming regions from diffuse gas and dust emission on ~50 pc scales. We find that the [CII] - SFR correlation holds even at these scales, although the relation typically has a flatter slope than found at larger (~kpc) scales. While the Halpha emission in M31 is concentrated in the SFR regions, we find that a significant amount (~20-90%) of the [CII] emission comes from outside star-forming regions, and that the total IR (TIR) emission has the highest diffuse fraction of all SFR tracers. We find a weak correlation of the [CII]/TIR to dust color in each Field, and find a large scale trend of increasing [CII]/TIR with galactocentric radius. The differences in the relative diffuse fractions of [CII], Halpha and IR tracers are likely caused by a combination of energetic photon leakage from HII regions and heating by the diffuse radiation field arising from older (B-star) stellar populations. However, we find that by averaging our measurements over ~kpc scales, these effects are minimized, and the relation between [CII] and SFR found in other nearby galaxy studies is retrieved.

Using data from the WISE mission, we have measured near infra-red (NIR) photometry of a diverse sample of dust-free stellar systems (globular clusters, dwarf and giant early-type galaxies) which have metallicities that span the range -2.2 < [Fe/H] (dex) < 0.3. This dramatically increases the sample size and broadens the metallicity regime over which the 3.4 (W1) and 4.6 micron (W2) photometry of stellar populations have been examined. We find that the W1 - W2 colors of intermediate and old (> 2 Gyr) stellar populations are insensitive to the age of the stellar population, but that the W1 - W2 colors become bluer with increasing metallicity, a trend not well reproduced by most stellar population synthesis (SPS) models. In common with previous studies, we attribute this behavior to the increasing strength of the CO absorption feature located in the 4.6 micron bandpass with metallicity. Having used our sample to validate the efficacy of some of the SPS models, we use these models to derive stellar mass-to-light ratios in the W1 and W2 bands. Utilizing observational data from the SAURON and ATLAS3D surveys, we demonstrate that these bands provide extremely simple, yet robust stellar mass tracers for dust free older stellar populations that are freed from many of the uncertainties common among optical estimators.

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.

We study the properties of the interstellar medium in the interacting system BR1202-0725 at z=4.7 via its [NII] and [CII] fine-structure line emission. This system consists of a QSO, a sub-mm galaxy (SMG), and two Ly-alpha emitters (LAEs). Such a diversity in galaxy properties makes BR1202-0725 a unique laboratory of star formation and galaxy evolution at high redshift. We present ionized nitrogen ([NII] 205 micron) observations of this system, obtained with the IRAM Plateau de Bure Interferometer. We find no [NII] emission at the quasar location, but tentative [NII] line detections associated with the SMG and one of the LAEs. Together with available ionized carbon ([CII] 158 micron) ALMA observations of this system, we find the following: The [CII]/[NII] luminosity ratio is >5.5 for the QSO and the SMG, but it is as low as ~2 in the LAE, suggesting that, in this source, most of the [CII] emission is associated with the ionized medium (HII regions) rather than the neutral one (PDRs). This study demonstrates the importance of combined studies of multiple fine--structure lines in order to pin down the physical properties of the interstellar medium in distant galaxies.

Spitzer Space Telescope and Herschel Space Observatory imaging of M31 is used, with a physical dust model, to construct maps of dust surface density, dust-to-gas ratio, starlight heating intensity, and PAH abundance, out to R=25kpc. The global dust mass is M_d=5.4x10^7Msol, the global dust/H mass ratio is M_d/M_H=0.0081, and the global PAH abundance is <q_PAH>=0.039. The dust surface density has an inner ring at R=5.6kpc, a maximum at R=11.2kpc, and an outer ring at R=15.1kpc. The dust/gas ratio varies from M_d/M_H=0.026 at the center to ~0.0027 at R=25kpc. From the dust/gas ratio, we estimate the ISM metallicity to vary by a factor ~10, from Z/Zsol=3 at R=0 to ~0.3 at R=25kpc. The dust heating rate parameter <U> peaks at the center, with <U> approx 35, declining to <U> approx 0.25 at R=20kpc. Within the central kpc, the starlight heating intensity inferred from the dust modeling is close to what is estimated from the stars in the bulge. The PAH abundance reaches a peak q_PAH=0.045 at R=11.2kpc. When allowance is made for the different spectrum of the bulge stars, q_PAH for the dust in the central kpc is similar to the overall value of q_PAH in the disk. The silicate-graphite-PAH dust model used here is generally able to reproduce the observed dust spectral energy distribution across M31, but overpredicts 500um emission at R=2-6kpc, suggesting that at R=2-6kpc, the dust opacity varies more steeply with frequency (with beta approx 2.3 between 200 and 600um) than in the model

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.

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

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 use the chemical evolution predictions of cosmological hydrodynamic simulations with our latest theoretical stellar population synthesis, photoionization and shock models to predict the strong line evolution of ensembles of galaxies from z=3 to the present day. In this paper, we focus on the brightest optical emission-line ratios, [NII]/H-alpha and [OIII]/H-beta. We use the optical diagnostic Baldwin-Phillips-Terlevich (BPT) diagram as a tool for investigating the spectral properties of ensembles of active galaxies. We use four redshift windows chosen to exploit new near-infrared multi-object spectrographs. We predict how the BPT diagram will appear in these four redshift windows given different sets of assumptions. We show that the position of star-forming galaxies on the BPT diagram traces the ISM conditions and radiation field in galaxies at a given redshift. Galaxies containing AGN form a mixing sequence with purely star-forming galaxies. This mixing sequence may change dramatically with cosmic time, due to the metallicity sensitivity of the optical emission-lines. Furthermore, the position of the mixing sequence may probe metallicity gradients in galaxies as a function of redshift, depending on the size of the AGN narrow line region. We apply our latest slow shock models for gas shocked by galactic-scale winds. We show that at high redshift, galactic wind shocks are clearly separated from AGN in line ratio space. Instead, shocks from galactic winds mimic high metallicity starburst galaxies. We discuss our models in the context of future large near-infrared spectroscopic surveys.

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.

Modern (sub-)millimeter interferometers enable the measurement of the cool gas and dust emission of high-redshift galaxies (z>5). However, at these redshifts the cosmic microwave background (CMB) temperature is higher, approaching, and even exceeding, the temperature of cold dust and molecular gas observed in the local Universe. In this paper, we discuss the impact of the warmer CMB on (sub-)millimeter observations of high-redshift galaxies. The CMB affects the observed (sub-)millimeter dust continuum and the line emission (e.g. carbon monoxide, CO) in two ways: (i) it provides an additional source of (both dust and gas) heating; and (ii) it is a non-negligible background against which the line and continuum emission are measured. We show that these two competing processes affect the way we interpret the dust and gas properties of high-redshift galaxies using spectral energy distribution models. We quantify these effects and provide correction factors to compute what fraction of the intrinsic dust (and line) emission can be detected against the CMB as a function of frequency, redshift and temperature. We discuss implications on the derived properties of high-redshift galaxies from (sub-)millimeter data. Specifically, the inferred dust and molecular gas masses can be severely underestimated for cold systems if the impact of the CMB is not properly taken into account.

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.

The goal of the present study is to establish the physical origin of dust heating and emission based on radiation transfer models, which self-consistently connect the emission components from diffuse dust and the dust in massive star forming regions. NGC 4214 is a nearby dwarf galaxy with a large set of ancillary data, ranging from the ultraviolet (UV) to radio, including maps from SPITZER, HERSCHEL and detections from PLANCK. We mapped this galaxy with MAMBO at 1.2 mm at the IRAM 30 m telescope. We extract separate dust emission components for the HII regions (plus their associated PDRs on pc scales) and for the diffuse dust (on kpc scales). We analyse the full UV to FIR/submm SED of the galaxy using a radiation transfer model which self-consistently treats the dust emission from diffuse and SF complexes components, considering the illumination of diffuse dust both by the distributed stellar populations, and by escaping light from the HII regions. While maintaining consistency with the framework of this model we additionally use a model that provides a detailed description of the dust emission from the HII regions and their surrounding PDRs on pc scales. Due to the large amount of available data and previous studies for NGC 4214 very few free parameters remained in the model fitting process. We achieve a satisfactory fit for the emission from HII+PDR regions on pc scales, with the exception of the emission at 8\mi, which is underpredicted by the model. For the diffuse emission we achieve a good fit if we assume that about 30-70% of the emission escaping the HII+PDR regions is able to leave the galaxy without passing through a diffuse ISM, which is not an unlikely scenario for a dwarf galaxy which has recently undergone a nuclear starburst. We determine a dust-to-gas mass ratio of 350-390 which is close to the expected value based on the metallicity.

NGC 604 is the second most massive H II region in the Local Group, thus an important laboratory for massive star formation. Using a combination of observational and analytical tools that include Spitzer spectroscopy, Herschel photometry, Chandra imaging, and Bayesian Spectral Energy Distribution fitting, we investigate the physical conditions in NGC 604, and quantify the amount of massive star formation currently taking place. We derive an average age of 4 +/- 1 Myr and a total stellar mass of 1.6 (+1.6)(-1.0) x 10^5 M_sun for the entire region, in agreement with previous optical studies. Across the region we find an effect of the X-ray field on both the abundance of aromatic molecules and the [Si II] emission. Within NGC 604 we identify several individual bright infrared sources with diameters of about 15 pc and luminosity weighted masses between 10^3 M_sun and 10^4 M_sun. Their spectral properties indicate that some of these sources are embedded clusters in process of formation, which together account for ~8% of the total stellar mass in the NGC 604 system. The variations of the radiation field strength across NGC 604 are consistent with a sequential star formation scenario, with at least two bursts in the last few million years. Our results indicate that massive star formation in NGC 604 is still ongoing, likely triggered by the earlier bursts.

We combine new dust continuum observations of the edge-on spiral galaxy NGC 4565 in all Herschel/SPIRE (250, 350, 500 micron) wavebands, obtained as part of the Herschel Reference Survey, and a large set of ancillary data (Spitzer, SDSS, GALEX) to analyze its dust energy balance. We fit a radiative transfer model for the stars and dust to the optical maps with the fitting algorithm FitSKIRT. To account for the observed UV and mid-infrared emission, this initial model was supplemented with both obscured and unobscured star-forming regions. Even though these star-forming complexes provide an additional heating source for the dust, the far-infrared/submillimeter emission long wards of 100 micron is underestimated by a factor of 3-4. This inconsistency in the dust energy budget of NGC 4565 suggests that a sizable fraction (two-thirds) of the total dust reservoir (Mdust ~ 2.9e+8 Msun) consists of a clumpy distribution with no associated young stellar sources. The distribution of those dense dust clouds would be in such a way that they remain unresolved in current far-infrared/submillimeter observations and hardly comtribute to the attenuation at optical wavelengths. More than two-thirds of the dust heating in NGC 4565 is powered by the old stellar population, with localized embedded sources supplying the remaining dust heating in NGC 4565. The results from this detailed dust energy balance study in NGC 4565 is consistent with that of similar analyses of other edge-on spirals.

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.

We use new Herschel multi-band imaging of the Andromeda galaxy to analyze how dust heating occurs in the central regions of galaxy spheroids that are essentially devoid of young stars. We construct a dust temperature map of M31 through fitting modified blackbody SEDs to the Herschel data, and find that the temperature within 2 kpc rises strongly from the mean value in the disk of 17 pm 1K to \sim35K at the centre. UV to near-IR imaging of the central few kpc shows directly the absence of young stellar populations, delineates the radial profile of the stellar density, and demonstrates that even the near-UV dust extinction is optically thin in M31's bulge. This allows the direct calculation of the stellar radiation heating in the bulge, U\ast(r), as a function of radius. The increasing temperature profile in the centre matches that expected from the stellar heating, i.e. that the dust heating and cooling rates track each other over nearly two orders of magnitude in U\ast. The modelled dust heating is in excess of the observed dust temperatures, suggesting that it is more than sufficient to explain the observed IR emission. Together with the wavelength dependent absorption cross section of the dust, this demonstrates directly that it is the optical, not UV, radiation that sets the heating rate. This analysis shows that neither young stellar populations nor stellar near-UV radiation are necessary to heat dust to warm temperatures in galaxy spheroids. Rather, it is the high densities of Gyr-old stellar populations that provide a sufficiently strong diffuse radiation field to heat the dust. To the extent which these results pertain to the tenuous dust found in the centres of early-type galaxies remains yet to be explored.

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 investigate the correlation between the far-infrared (FIR) and radio continuum emission from NGC6946 on spatial scales between 0.9 and 17 kpc. We use the Herschel PACS (70, 100, 160$\mu$m) and SPIRE (250$\mu$m) data from the KINGFISH project. Separating the free-free and synchrotron components of the radio continuum emission, we find that FIR is better correlated with the free-free than the synchrotron emission. Compared to a similar study in M33 and M31, we find that the scale dependence of the synchrotron--FIR correlation in NGC6946 is more similar to M31 than M33. The scale dependence of the synchrotron--FIR correlation can be explained by the turbulent-to-ordered magnetic field ratio or, equivalently, the diffusion length of the cosmic ray electrons in these galaxies.

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 have carried out a detailed modeling of the dust Spectral Energy Distribution (SED) of the nearby, starbursting dwarf galaxy NGC 4214. A key point of our modeling is that we distinguish the emission from (i) HII regions and their associated photodissociation regions (PDRs) and (ii) diffuse dust. For both components we apply templates from the literature calculated with a realistic geometry and including radiation transfer. The large amount of existing data from the ultraviolet (UV) to the radio allows the direct measurement of most of the input parameters of the models. We achieve a good fit for the total dust SED of NGC 4214. In the present contribution we describe the available data, the data reduction and the determination of the model parameters, whereas a description of the general outline of our work is presented in the contribution of Lisenfeld et al. in this volume.

We have carried out a detailed modelling of the dust heating and emission in the nearby, starbursting dwarf galaxy NGC 4214. Due to its proximity and the great wealth of data from the UV to the millimeter range (from GALEX, HST, {\it Spitzer}, Herschel, Planck and IRAM) it is possible to separately model the emission from HII regions and their associated photodissociation regions (PDRs) and the emission from diffuse dust. Furthermore, most model parameters can be directly determined from the data leaving very few free parameters. We can fit both the emission from HII+PDR regions and the diffuse emission in NGC 4214 with these models with "normal" dust properties and realistic parameters.

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.

High resolution spectra are necessary to distinguish and correctly measure the Balmer emission lines due to the presence of strong metal and Balmer absorption features in the stellar continuum. This accurate measurement is necessary for use in emission line diagnostics, such as the Balmer decrement (i.e. Halpha/Hbeta), used to determine the attenuation of galaxies. Yet at high redshifts obtaining such spectra becomes costly. Balmer emission line equivalent widths are much easier to measure, requiring only low resolution spectra or even simple narrow band filters and therefore shorter observation times. However a correction for the stellar continuum is still needed for this equivalent width Balmer decrement. We present here a statistical analysis of the Sloan Digital Sky Survey Data Release 7 emission line galaxy sample, using the spectrally determined Balmer emission line fluxes and equivalent widths. Using the large numbers of galaxies available in the SDSS catalogue, we determined an equivalent width Balmer decrement including a statistically-based correction for the stellar continuum. Based on this formula, the attenuation of galaxies can now be obtained from low spectral resolution observations. In addition, this investigation also revealed an error in the Hbeta line fluxes, within the SDSS DR7 MPA/JHU catalogue, with the equivalent widths underestimated by average ~0.35A in the emission line galaxy sample. This error means that Balmer decrement determined attenuations are overestimated by a systematic 0.1 magnitudes in A_V, and future analyses of this sample need to include this correction.

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.

To understand and interpret the observed Spectral Energy Distributions (SEDs) of starbursts, theoretical or semi-empirical SED models are necessary. Yet, while they are well-founded in theory, independent verification and calibration of these models, including the exploration of possible degeneracies between their parameters, are rarely made. As a consequence, a robust fitting method that leads to unique and reproducible results has been lacking. Here we introduce a novel approach based on Bayesian analysis to fit the Spitzer-IRS spectra of starbursts using the SED models proposed by Groves et al. (2008). We demonstrate its capabilities and verify the agreement between the derived best fit parameters and actual physical conditions by modelling the nearby, well-studied, giant HII region 30 Dor in the LMC. The derived physical parameters, such as cluster mass, cluster age, ISM pressure and covering fraction of photodissociation regions, are representative of the 30 Dor region. The inclusion of the emission lines in the modelling is crucial to break degeneracies. We investigate the limitations and uncertainties by modelling sub-regions, which are dominated by single components, within 30 Dor. A remarkable result for 30 Doradus in particular is a considerable contribution to its mid-infrared spectrum from hot ({\simeq} 300K) dust. The demonstrated success of our approach will allow us to derive the physical conditions in more distant, spatially unresolved starbursts.

A large body of evidence has demonstrated that the global rest-frame optical and IR colours of galaxies correlate well with each other, as well as with other galactic properties such as surface brightness and morphology. However the processes that lead to the observed correlations are contrary; the stellar light that contributes to the optical is readily absorbed by dust which emits in the IR. Thus on small scales we expect these correlations to break down. We examine seven nearby galaxies ranging from early- to late-types, on a pixel-by-pixel basis and we demonstrate that there is disconnect between the optical and IR when normalized to the near-IR (H-band). We can decompose this disconnect into two distinct components through a Principal Component Analysis of the H-band normalized SED of the pixels: one mainly correlated with variations in the IR, the other correlated with variations in the optical. By mapping these two components it is clear they arise from distinct spatial regions. The IR dominated component is strongly associated with the specific star-formation rate, while the optical-dominated component is broadly associated with the stellar mass density. However, when the pixels of all galaxies are compared, the well known optical-IR colour correlations return, demonstrating that the variance observed within galaxies is around a mean which follows the well-known trend. We also examine the extremely strong correlations between the IRAC-NIR colours and demonstrate that they are tight enough to use a single IRAC-NIR colour (i.e. 8mum-H) to determine the fluxes in the other IRAC bands. These correlations arise from the differing contribution of stellar light and dust to the IRAC bands, enabling us to determine pure stellar colours for these bands, but still demonstrating the need for dust (or stellar) corrections in these bands when being used as stellar (dust) tracers.

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.

Optical nebular emission lines are commonly used to estimate the star formation rate of galaxies and the black hole accretion rate of their central active nucleus. The accuracy of the conversion from line strengths to physical properties depends upon the accuracy to which the lines can be corrected for dust attenuation. For studies of single galaxies with normal amounts of dust, most dust corrections result in the same derived properties within the errors. However, for statistical studies of populations of galaxies, or for studies of galaxies with higher dust contents such as might be found in some classes of "transition" galaxies, significant uncertainty arises from the dust attenuation correction. We compare the strength of the predominantly unobscured mid-IR [NeII]15.5um + [NeIII]12.8um emission lines to the optical H alpha emission lines in four samples of galaxies: (i) ordinary star forming galaxies, (ii) optically selected dusty galaxies, (iii) ULIRGs, (iv) Seyfert 2 galaxies. We show that a single dust attenuation curve applied to all samples can correct H alpha emission for dust attenuation to a factor better than 2. Similarly, we compare mid-IR [OIV] and optical [OIII] luminosities to find that [OIII] can be corrected to a factor better than 3. This shows that the total dust attenuation suffered by the AGN narrow line region is not significantly different to that suffered by the starforming HII regions in the galaxy. We provide explicit dust attenuation corrections, together with errors, for [OII], [OIII] and H alpha. The best-fit average attenuation curve is slightly greyer than the Milky-Way extinction law, indicating either that external galaxies have slightly different typical dust properties to the Milky Way, or that there is a significant contribution from scattering. Finally, we uncover an intriguing correlation between Silicate absorption and Balmer decrement.

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 present a physical model for the origin of z~2 Dust-Obscured Galaxies (DOGs), a class of high-redshift ULIRGs selected at 24 micron which are particularly optically faint (24/R>1000). By combining N-body/SPH simulations of high redshift galaxy evolution with 3D polychromatic dust radiative transfer models, we find that luminous DOGs (with F24 > 0.3 mJy at z~2 are well-modeled as extreme gas-rich mergers in massive (~5x10^12-10^13 Msun) halos, with elevated star formation rates (~500-1000 Msun/yr) and/or significant AGN growth (Mdot > 0.5 Msun/yr), whereas less luminous DOGs are more diverse in nature. At final coalescence, merger-driven DOGs transition from being starburst dominated to AGN dominated, evolving from a "bump" to a power-law shaped mid-IR (IRAC) spectral energy distribution (SED). After the DOG phase, the galaxy settles back to exhibiting a "bump" SED with bluer colors and lower star formation rates. While canonically power-law galaxies are associated with being AGN-dominated, we find that the power-law mid-IR SED can owe both to direct AGN contribution, as well as to a heavily dust obscured stellar bump at times that the galaxy is starburst dominated. Thus power-law galaxies can be either starburst or AGN dominated. Less luminous DOGs can be well-represented either by mergers, or by massive ($M_{\rm baryon} ~5x10^11 Msun) secularly evolving gas-rich disc galaxies (with SFR > 50 Msun/yr). By utilising similar models as those employed in the SMG formation study of Narayanan et al. (2010), we investigate the connection between DOGs and SMGs. We find that the most heavily star-forming merger driven DOGs can be selected as Submillimetre Galaxies (SMGs), while both merger-driven and secularly evolving DOGs typically satisfy the BzK selection criteria.

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.

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.

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

We present a new software tool to enable astronomers to easily compare observations of emission line ratios with those determined by photoionization and shock models, ITERA, the IDL Tool for Emission-line Ratio Analysis. This tool can plot ratios of emission lines predicted by models and allows for comparison of observed line ratios against grids of these models selected from model libraries associated with the tool. We provide details of the libraries of standard photoionization and shock models available with ITERA, and, in addition, present three example emission line ratio diagrams covering a range of wavelengths to demonstrate the capabilities of ITERA. ITERA, and associated libraries, is available from \url{http://www.brentgroves.net/itera.html}

Large optical surveys provide an unprecedented census of galaxies in the local Universe, forming an invaluable framework into which more detailed studies of objects can be placed. But how useful are optical surveys for understanding the co-evolution of black holes and galaxies, given their limited wavelength coverage, selection criteria, and depth? In this conference paper I present work-in-progress comparing optical and mid-IR diagnostics of three "unusual" low redshift populations (luminous Seyferts, dusty Balmer-strong AGN, ULIRGs) with a set of ordinary star-forming galaxies from the SDSS. I address the questions: How well do the mid-infrared and optical diagnostics of star formation and AGN strength agree? To what extent do optical surveys allow us to include extreme, dusty, morphologically disturbed galaxies in our "complete" census of black hole-galaxy co-evolution?

We describe a model for the formation of \zsim 2 Submillimeter Galaxies (SMGs) which simultaneously accounts for both average and bright SMGs while providing a reasonable match to their mean observed spectral energy distributions (SEDs). By coupling hydrodynamic simulations of galaxy mergers with the high resolution 3D polychromatic radiative transfer code Sunrise, we find that a mass sequence of merger models which use observational constraints as physical input naturally yield objects which exhibit black hole, bulge, and H2 gas masses similar to those observed in SMGs. The dominant drivers behind the 850 micron flux are the masses of the merging galaxies and the stellar birthcloud covering fraction. The most luminous (S850 ~ 15 mJy) sources are recovered by ~10^13 Msun 1:1 major mergers with a birthcloud covering fraction close to unity, whereas more average SMGs ~5-7 mJy) may be formed in lower mass halos ~5x10^12 Msun. These models demonstrate the need for high spatial resolution hydrodynamic and radiative transfer simulations in matching both the most luminous sources as well as the full SEDs of SMGs. While these models suggest a natural formation mechanism for SMGs, they do not attempt to match cosmological statistics of galaxy populations; future efforts along this line will help ascertain the robustness of these models.

An updated version of the dust radiation transfer code Sunrise, including models for star-forming regions and a self-consistent calculation of the spatially dependent dust and PAH emission, is presented. Given a hydrodynamic simulation of a galaxy, this model can calculate a realistic 2-dimensional ultraviolet--submillimeter spectral energy distribution of the galaxy, including emission lines from HII regions, from any viewpoint. To model the emission from star-forming regions, the MAPPINGSIII photoionization code is used. The high wavelength resolution (~ 1000 wavelengths) is made possible by the polychromatic Monte-Carlo algorithm employed by Sunrise. From the 2-D spectral energy distributions, images in any filter bands or integrated galaxy SEDs can be created. Using a suite of hydrodynamic simulations of disc galaxies, the output broad-band images and spectral energy distributions are compared with observed galaxies from the multiwavelength SINGS and SLUGS galaxy surveys. Overall, the output spectral energy distributions show a good match with observed galaxies in colours ranging from GALEX far-UV to SCUBA submillimeter wavelengths. The only possible exception is the 160 micron/850 micron colour, which the simulations underestimate by a factor "of order 5" compared to the SINGS sample. However, the simulations here agree with the SLUGS galaxies, which consistently have significantly larger amounts of cold dust than the SINGS galaxies. The Sunrise model can be used to generate simulated observations of arbitrary hydrodynamic galaxy simulations. In this way, predictions of galaxy formation theories can be directly tested against observations of galaxies.

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

We present an analysis of deep mid- to far-infrared (MFIR) Spitzer photometric observations of the southern 2Jy sample of powerful radio sources (0.05 < z < 0.7), conducting a statistical investigation of the links between radio jet, AGN, starburst activity and MFIR properties. This is part of an ongoing extensive study of powerful radio galaxies that benefits from both complete optical emission line information and a uniquely high detection rate in the far-infrared (far-IR). We find tight correlations between the MFIR and [OIII] emission luminosities, which are significantly better than those between MFIR and extended radio luminosities, or between radio and [OIII] luminosities. Since [OIII] is a known indicator of intrinsic AGN power, these correlations confirm AGN illumination of the circum-nuclear dust as the primary heating mechanism for the dust producing thermal MFIR emission at both 24 and 70 microns. We demonstrate that AGN heating is energetically feasible, and identify the narrow line region clouds as the most likely location of the cool, far-IR emitting dust. Starbursts make a major contribution to the heating of the cool dust in only 15-28% of our targets. We also investigate the orientation dependence of the continuum properties, finding that the broad- and narrow-line objects in our sample with strong emission lines have similar distributions of MFIR luminosities and colours. Therefore our results are entirely consistent with the orientation-based unified schemes for powerful radio galaxies. However, the weak line radio galaxies (WLRG) form a separate class of objects with intrinsically low luminosity AGN in which both the optical emission lines and the MFIR continuum are weak.

The mid-infrared ratio [NeIII]15.6mum/[NeII]12.8mum is a strong diagnostic of the ionization state of emission line objects, due to its use of only strong neon emission lines only weakly affected by extinction. However this ratio is not available to ground-based telescopes as only a few spectroscopic windows are available in the MIR. To deal with this problem we aimed to verify if there exists a conversion law between ground-accessible, strong MIR line ratio [SIV]/[NeII] and the diagnostic [NeIII]/[NeII] ratio that can serve as a reference for future ground-based observations. We collated the [SIV]10.5mum, [NeII]12.8mum, [NeIII]15.6\mum and [SIII]18.7mum emission line fluxes from a wide range of sources in the rich Spitzer and ISO archives, and compared the [NeIII]/[NeII], [SIV]/[SIII], and [SIV]/[NeII] ratios. We find a strong correlation between the [SIV]/[NeII] and [\neiii]/[\neii] ratio, with a linear fit of log([NeIII]/[NeII]) = 0.81log([SIV]/[NeII])+0.36, accurate to a factor of ~2 over four orders of magnitude in the line ratios. This demonstrates clearly the ability of ground-based infrared spectrographs to do ionization studies of nebulae.

We present a new library of fully-radiative shock models calculated with the MAPPINGS III shock and photoionization code. The library consists of grids of models with shock velocities in the range v=100-1000 km/s and magnetic parameters B/sqrt(n) of 10^-4 - 10 muG cm^(3/2) for five different atomic abundance sets, and for a pre-shock density of 1.0 cm^(-3). Additionally, Solar abundance model grids have been calculated for densities of 0.01, 0.1, 10, 100, and 1000 cm^(-3) with the same range in v and B/sqrt(n). Each model includes components of both the radiative shock and its photoionized precursor, ionized by the EUV and soft X-ray radiation generated in the radiative gas. We present the details of the ionization structure, the column densities, and the luminosities of the shock and its precursor. Emission line ratio predictions are separately given for the shock and its precursor as well as for the composite shock+precursor structure to facilitate comparison with observations in cases where the shock and its precursor are not resolved. Emission line ratio grids for shock and shock+precursor are presented on standard line ratio diagnostic diagrams, and we compare these grids to observations of radio galaxies and a sample of AGN and star forming galaxies from the Sloan Digital Sky Survey. This library is available online, along with a suite of tools to enable the analysis of the shocks and the easy creation of emission line ratio diagnostic diagrams. These models represent a significant increase in parameter space coverage over previously available models, and therefore provide a unique tool in the diagnosis of emission by shocks.

This paper addresses the origin of the silicate emission observed in PG QSOs, based on observations with the Spitzer Space Telescope. Scenarios based on the unified model suggest that silicate emission in AGN arises mainly from the illuminated faces of the clouds in the torus at temperatures near sublimation. However, detections of silicate emission in Type 2 QSOs, and the estimated cool dust temperatures, argue for a more extended emission region.To investigate this issue we present the mid-infrared spectra of 23 QSOs. These spectra, and especially the silicate emission features at ~10 and ~18 mu can be fitted using dusty narrow line region (NLR) models and a combination of black bodies. The bolometric luminosities of the QSOs allow us to derive the radial distances and covering factors for the silicate-emitting dust. The inferred radii are 100-200 times larger than the dust sublimation radius, much larger than the expected dimensions of the inner torus. Our QSO mid-IR spectra are consistent with the bulk of the silicate dust emission arising from the dust in the innermost parts of the NLR.

We combine the the stellar spectral synthesis code Starburst99, the nebular modelling code MAPPINGSIII, and a 1-D dynamical evolution model of HII regions around massive clusters of young stars to generate improved models of the spectral energy distribution (SED) of starburst galaxies. We introduce a compactness parameter, C, which characterizes the specific intensity of the radiation field at ionization fronts in HII regions, and which controls the shape of the far-IR dust re-emission, often referred to loosely as the dust ``temperature''. We also investigate the effect of metallicity on the overall SED and in particular, on the strength of the PAH features. We provide templates for the mean emission produced by the young compact HII regions, the older (10 - 100 Myr) stars and for the wavelength-dependent attenuation produced by a foreground screen of the dust used in our model. We demonstrate that these components may be combined to produce a excellent fit to the observed SEDs of star formation dominated galaxies which are often used as templates (Arp 220 and NGC 6240). This fit extends from the Lyman Limit to wavelengths of about one mm. The methods presented in both this paper and in the previous papers of this series allow the extraction of the physical parameters of the starburst region (star formation rates, star formation rate history, mean cluster mass, metallicity, dust attenuation and pressure) from the analysis of the pan-spectral SED.

We present high spatial resolution (~ 35 parsec) 5-38 um spectra of the central region of M82, taken with the Spitzer Infrared Spectrograph. From these spectra we determined the fluxes and equivalent widths of key diagnostic features, such as the [NeII]12.8um, [NeIII]15.5um, and H_2 S(1)17.03um lines, and the broad mid-IR polycyclic aromatic hydrocarbon (PAH) emission features in six representative regions and analysed the spatial distribution of these lines and their ratios across the central region. We find a good correlation of the dust extinction with the CO 1-0 emission. The PAH emission follows closely the ionization structure along the galactic disk. The observed variations of the diagnostic PAH ratios across M82 can be explained by extinction effects, within systematic uncertainties. The 16-18um PAH complex is very prominent, and its equivalent width is enhanced outwards from the galactic plane. We interpret this as a consequence of the variation of the UV radiation field. The EWs of the 11.3um PAH feature and the H_2 S(1) line correlate closely, and we conclude that shocks in the outflow regions have no measurable influence on the H_2 emission. The [NeIII]/[NeII] ratio is on average low at ~0.18, and shows little variations across the plane, indicating that the dominant stellar population is evolved (5 - 6 Myr) and well distributed. There is a slight increase of the ratio with distance from the galactic plane of M82 which we attribute to a decrease in gas density. Our observations indicate that the star formation rate has decreased significantly in the last 5 Myr. The quantities of dust and molecular gas in the central area of the galaxy argue against starvation and for negative feedback processes, observable through the strong extra-planar outflows.

Using the large emission line galaxy sample from the Sloan Digital Sky Survey we show that Star forming galaxies, Seyferts, and low-ionization nuclear emission-line regions (LINERs) form clearly separated branches on the standard optical diagnostic diagrams. We derive a new empirical classification scheme which cleanly separates these emission-line galaxies, using strong optical emission lines. Using this classification we identify a few distinguishing host galaxy properties of each class, which, along with the emission line analysis, suggest continuous evolution from one class to another. As a final note, we introduce models of both Starforming galaxies and AGN narrow line regions which can explain the distribution of galaxies on standard emission line ratio diagrams, and possibly suggest new diagnostics across the emission spectrum.

I present a broad overview of modelling of the Narrow Line Region (NLR) of active galaxies, and discuss some of the more recent models we currently have for the emission from the NLR. I show why the emission line ratios from the NLR are constrained to certain observed values, and describe what physical parameters we can derive from observations using emission line models. Also presented are some examples of this, looking at the metallicity and excitation mechanism of active galaxies. As a final point, the limitations of the current models are discussed, and how how the combination of modelling and theory can help us solve some of the questions that still remain within the NLR.

We present models for the mid- and far- infrared emission from the Narrow Line Region (NLR) of Active Galactic Nuclei (AGN). Using the MAPPINGS III code we explore the effect of typical NLR parameters on the spectral characteristics of the IR emission. These include useful IR emission line ratio diagnostic diagrams for the determination of these parameters, as well as Star formation--AGN mixing diagnostics. We also examine emission line to continuum correlations which would assist in separating the IR emission arising from the NLR from that coming from the inner torus. We find for AGN like NGC 1068 and NGC 4151 that the NLR only contributes ~10% to the total IRAS 25 mum flux, and that other components such as a dusty torus are necessary to explain the total AGN IR emission.

We build, as far as theory will permit, self consistent model HII regions around central clusters of aging stars. These produce strong emission line diagnostics applicable to either individual HII regions in galaxies, or to the integrated emission line spectra of disk or starburst galaxies. The models assume that the expansion and internal pressure of individual HII regions is driven by the net input of mechanical energy from the central cluster, be it through winds or supernova events. This eliminates the ionization parameter as a free variable, replacing it with a parameter which depends on the ratio of the cluster mass to the pressure in the surrounding interstellar medium. These models explain why HII regions with low abundances have high excitation, and demonstrate that at least part of the warm ionized medium is the result of overlapping faint, old, large, and low pressure HII regions. We present line ratios (at both optical and IR wavelengths) which provide reliable abundance diagnostics for both single HII regions or for integrated galaxy spectra, and we find a number that can be used to estimate the mean age of the cluster stars exciting individual HII regions.

Current emission-line based estimates of the metallicity of active galactic nuclei (AGN) at both high and low redshifts indicate that AGN have predominantly solar to supersolar metallicities. This leads to the question: do low metallicity AGN exist? In this paper we use photoionization models to examine the effects of metallicity variations on the narrow emission lines from an AGN. We explore a variety of emission-line diagnostics that are useful for identifying AGN with low metallicity gas. We find that line ratios involving [NII] are the most robust metallicity indicators in galaxies where the primary source of ionization is from the active nucleus. Ratios involving [SII] and [OI] are strongly affected by uncertainties in modelling the density structure of the narrow line clouds. To test our diagnostics, we turn to an analysis of AGN in the Sloan Digital Sky Survey (SDSS). We find a clear trend in the relative strength of [NII] with the mass of the AGN host galaxy. The metallicity of the ISM is known to be correlated with stellar mass in star-forming galaxies; our results indicate that a similar trend exists for AGN. We also find that the best-fit models for typical Seyfert narrow line regions have supersolar abundances. Although there is a mass-dependent range of a factor of 2-3 in the NLR metallicities of the AGN in our sample, AGN with sub-solar metallicities are very rare in the SDSS. Out of a sample of ~23000 Seyfert 2 galaxies we find only ~40 clear candidates for AGN with NLR abundances that are below solar.

We examine, from a theoretical viewpoint, how the physical parameters of HII regions are controlled in both normal galaxies and in starburst environments. These parameters are the HII region luminosity function, the time-dependent size, the covering fraction of molecular clouds, the pressure in the ionized gas and the ionization parameter. The factors which control them are the initial mass function of the exciting stars, the cluster mass function, the metallicity and the mean pressure in the surrounding interstellar medium. We investigate the sensitivity of the H$\alpha$ luminosity to the IMF, and find that this can translate to more than a factor two variation in derived star formation rates. The molecular cloud dissipation timescale is estimated from a case study of M17 to be $\sim1$ Myr for this object. Based upon HII luminosity function fitting for nearby galaxies, we suggest that the \HII region cluster mass function is fitted by a log-normal form peaking at $\sim 100 M_{\odot}$. The cluster mass function continues the stellar IMF to higher mass. The pressure in the HII regions is controlled by the mechanical luminosity flux from the central cluster. Since this is closely related to the ionizing photon flux, we show that the ionization parameter is not a free variable, and that the diffuse ionized medium may be composed of many large, faint and old HII regions. Finally, we derive theoretical probability distributions for the ionization parameter as a function of metallicity and compare these to those derived for SDSS galaxies.

We present an analysis of the host properties of 85224 emission-line galaxies from the Sloan Digital Sky Survey. We show that Seyferts and LINERs form clearly separated branches on the standard optical diagnostic diagrams. We derive a new empirical classification scheme which cleanly separates star-forming galaxies, composite AGN-HII galaxies, Seyferts and LINERs and we study the host galaxy properties of these different classes of objects. LINERs are older, more massive, less dusty and more concentrated, and have higher velocity dispersions and lower [OIII] luminosities than Seyfert galaxies. We consider the quantity L[OIII]/sigma^4, which is an indicator of the black hole accretion rate relative to the Eddington rate. Remarkably, we find that at fixed L[OIII]/sigma^4, all differences between Seyfert and LINER host properties disappear. LINERs and Seyferts form a continuous sequence, with LINERs dominant at low L/L(EDD) and Seyferts dominant at high L/L(EDD). These results suggest that the majority of LINERs are AGN and that the Seyfert/LINER dichotomy is analogous to the high/low-state transition for X-ray binary systems. We apply theoretical photoionization models and show that pure LINERs require a harder ionizing radiation field with lower ionization parameter than Seyfert galaxies, consistent with the low and high X-ray binary states.

In this paper, we combine the stellar spectral synthesis code STARBURST 99, the nebular modelling code MAPPINGS IIIq, a 1-D dynamical evolution model of \HII regions around massive clusters of young stars and a simplified model of synchrotron emissivity to produce purely theoretical self-consistent synthetic spectral energy distributions (SEDs) for (solar metallicity) starbursts lasting some $10^8$ years. These SEDs extend from the Lyman Limit to beyond 21 cm. We find that two ISM parameters control the form of the SED; the pressure in the diffuse phase of the ISM (or, equivalently, its density), and the molecular cloud dissipation timescale. We present detailed SED fits to Arp 220 and NGC 6240, and we give the predicted colors for starburst galaxies derived from our models for the IRAS and the Spitzer Space Observatory MIPS and IRAC instruments. Our models reproduce the spread in observed colors of starburst galaxies. Finally, we present absolute calibrations to convert observed fluxes into star formation rates in the UV (GALEX), at optical wavelengths (H$\alpha$), and in the IR (IRAS or the Spitzer Space Observatory). (Abstract Truncated)

We present ultraviolet emission-line maps of the narrow-line region (NLR) of NGC 1068. The maps span 115--318 nm, the biconical ionization cone, several posited jet/ISM interactions, and the compact knots whose optical spectra we reported previously resemble kinematically the quasar Associated Absorption Line systems. Across the NLR, we find that ultraviolet flux ratios are consistent with photoionization, not shock excitation, even for gas blueshifted abruptly to 3000 \kms relative to galaxy systemic velocity or for gas projected near the radio jet. The knots may be radiatively accelerated, photoablated fragments of molecular clouds.

Seyfert narrow line regions (NLR) have emission line ratios which are remarkably uniform, displaying only ~0.5 dex variation between different galaxies. Existing models have been unable to explain these observations without the introduction of ad hoc assumptions, geometrical restrictions or new parameters. Here we introduce a new model: dusty radiation pressure dominated photoionization, which provides a natural self-regulating characteristic leading to an invariance of the spectrum over a very wide range (>100) of ionization parameter. The dusty model is able to reproduce both the range and the absolute value of the observational line ratios not only in the standard optical diagnostic diagrams but also in UV diagnostic plots, providing an explanation to the problem in NLR observations.

Seyfert narrow line region (NLR) emission line ratios are remarkably uniform, displaying only ~0.5 dex variation between galaxies, and even less within an individual object. Previous photoionization and shock models of this region were unable to explain this observation without the introduction of arbitrary assumptions or additional parameters. Dusty, radiation pressure dominated photoionization models provide a simple physical mechanism which can reproduce this spectral uniformity between different objects. In the first paper of this series we described this model and its implementation in detail, as well as presenting grids of model emission lines and examining the model structures. Here we explore these models further, demonstrating their ability to reproduce the observed Seyfert line ratios on standard line diagnostic diagrams in both the optical and UV. We also investigate the effects that the variation of metallicity, density and ionizing spectrum have upon both the new paradigm and the standard photoionization models used hitherto. Along with the standard diagnostic diagrams we provide several new diagnostic diagrams in the UV, Optical and IR. These new diagrams can provide further tests of the dusty, radiation pressure photoionization paradigm as well as being used as diagnostics of the metallicity, density and ionizing spectrum of the emission line clouds.

We present the implementation of dusty, radiation pressure dominated photoionization models applicable to the Narrow Line Regions (NLRs) of Active Galactic Nuclei, using the \mapiii code. We give a grid of the predicted intensities of the most commonly used diagnostic spectral lines in the UV, optical and IR, covering a wide range of density, metallicity, the power-law index characterizing the photoionizing source and photoionization parameter, for use in the diagnosis of NLRs. We examine the temperature, density and ionization structure of these models, investigating the effect of variation of these parameters in order to gain a better understanding of NLR clouds themselves.

The radio - far infrared correlation is one of the tightest correlations found in astronomy. Many of the models explaining this correlation rely on the association of of global magnetic field strength with gas density. In this letter we put forward that the physical reason for this association lies within the processes of magnetohydrodynamic turbulence.

Star formation rates in galaxies are frequently estimated using the Balmer line fluxes. However, these can be systematically underestimated because dust competes for the absorption of Lyman continuum photons in the ionized gas. Here we present theoretical correction factors in a simple analytic form. T These factors scale as the product of the ionization parameter, ${\cal U}$, and the nebular O/H abundance ratio, both of which can now be derived from the observation of bright nebular line ratios. The correction factors are only somewhat dependent upon the photoelectron production by grains, but are very sensitive to the presence of complex PAH-like carbonaceous molecules in the ionized gas, providing that these can survive in such an environment.

In this paper we present a new method for obtaining the optical wavelength-dependent reddening function of planetary nebulae, using the nebular and stellar continuum. The data used was a spectrum of NGC 6302 obtained with a mean signal to noise of >10^2 A^-1 in the nebular continuum. With such a high S/N the continuum can be accurately compared with a theoretical model nebular plus stellar continuum. The nebular electron temperature and density used in the model are determined using ratios of prominent emission lines. The reddening function can then be obtained from the ratio of the theoretical and the observed continuum. We find that for NGC 6302, the visible to IR extinction law is indistinguishable from `standard' interstellar reddening, but that the UV extinction curve is much steeper than normal, suggesting that more small dust grains had been ejected into the nebula by the PN central star. Finally, using the extinction law that we have determined, we present a complete de--reddened line list of nearly 600 emission lines, and report on the detection of the He(2-10) and He(2-8) Raman Features at 4331 A and 4852 A, and the detection of Raman-Scattered OVI features at 6830 and 7087 AA.

The remarkable similarity between emission spectra of narrow line regions (NLR) in Seyfert Galaxies has long presented a mystery. In photoionization models, this similarity implies that the ionization parameter is nearly always the same, about U ~ 0.01. Here we present dusty, radiation-pressure dominated photoionization models that can provide natural physical insight into this problem. In these models, dust and the radiation pressure acting on it provide the controlling factor in moderating the density, excitation and surface brightness of photoionized NLR structures. Additionally, photoelectric heating by the dust is important in determining the temperature structure of the models. These models can also explain the coexistence of the low-, intermediate- and coronal ionization zones within a single self-consistent physical structure. The radiation pressure acting on dust may also be capable of driving the fast (~3000 km/s) outflows such as are seen in the HST observations of NGC 1068.

Using the HST STIS spectrograph we have obtained a grid of [O III] and H-beta emission-line spectra at 0"05x0"19 and 60 km/s (FWHM) resolution that covers much of the NLR of NGC 1068. We find emitting knots that have blueshifted radial velocities up to 3200 km/s relative to galaxy systemic, are 70-150 pc NE of the nucleus and up to 40 pc from the radio jet, emit several percent of the NLR line flux but no significant continuum, span velocity extents of up to 1250 km/s but a small fraction of the sky seen from the nucleus, coincide with a region of enhanced IR coronal-line emission, and have ionized masses $\sim$200 Msun/ne4 (ne4=10^4 cm^{-3}). We argue that the blueshifted knots are ablata from disintegrating molecular clouds that are being photoionized by the AGN, and are being accelerated readiatively by the AGN or mechanically by the radio jet. In their kinematic properties, the knots resemble the associated absorbers seen projected on the UV continua of some AGN. Between 2"5-4"5 from the nucleus, emission is redshifted relative to systemic, a pattern that we interpret as gas in the galaxy disk being pushed away from us by the NE radio lobe.

The ambient interstellar environment of wind- and supernova-driven superbubbles strongly affects their evolution, but its properties are rarely well-determined. We have therefore obtained HI aperture synthesis imaging of the environment around three similar, optically-selected superbubble nebulae in the Large Magellanic Cloud. The resulting HI maps show that the ambient gas distribution around these superbubbles differ to an extreme: DEM L25 shows no neutral shell component, but is nestled within an HI hole; DEM L50 shows a massive neutral shell component, but is otherwise within an HI void; and DEM L301 shows no correspondence at all between the optical nebula and HI distribution. There is also poor correspondence between the HI and optical kinematics. These results strongly caution against inferring properties of the ambient neutral environment of individual superbubbles without direct observations. Finally, all three objects show some evidence of shock activity.

Gamma-ray emission from AGN provides us with unprecedented insights into the physics of extragalactic jets. The emission from these jets fits naturally into the existing theoretical framework of relativistic jets as inferred from parsec scale and kiloparsec scale observations. Models of the $\gamma$-ray fluxes give us important knowledge of jet parameters, in particular magnetic fields and the size of the emitting region. For example, in the case of Markarian~501, the jet plasma is subequipartiion and the energy flux is particle dominated. Discussion of the physics of jets on the sub-parsec scale leads naturally to questions of jet composition. The dynamics of jets are consistent with conventional electron-proton composition if the minimum Lorentz factors, $\gamma_{\rm min} > 100$. Jets can be electron-positron if $\gamma_{\rm min} \approx 100$.