sort results by

Use logical operators AND, OR, NOT and round brackets to construct complex queries. Whitespace-separated words are treated as ANDed.

Show articles per page in mode

Meidt, Sharon E.

Normalized to: Meidt, S.

86 article(s) in total. 391 co-authors, from 1 to 61 common article(s). Median position in authors list is 6,0.

[1]  oai:arXiv.org:2003.01131  [pdf] - 2131876
When Gas Dynamics Decouples from Galactic Rotation: Characterizing ISM Circulation in Disk Galaxies
Comments: The Astrophysical Journal, 20 pages, 16 figures
Submitted: 2020-03-02, last modified: 2020-07-10
In galactic disks, galactic rotation sets the bulk motion of gas, and its energy and momentum can be transferred toward small scales. Additionally, in the interstellar medium, random and noncircular motions arise from stellar feedback, cloud-cloud interactions, and instabilities, among other processes. Our aim is to comprehend to which extent small-scale gas dynamics is decoupled from galactic rotation. We study the relative contributions of galactic rotation and local noncircular motions to the circulation of gas, $\Gamma$, a macroscopic measure of local rotation, defined as the line integral of the velocity field around a closed path. We measure the circulation distribution as a function of spatial scale in a set of simulated disk galaxies and we model the velocity field as the sum of galactic rotation and a Gaussian random field. The random field is parameterized by a broken power law in Fourier space, with a break at the scale $\lambda_c$. We define the spatial scale $\lambda_{\rm eq}$ at which galactic rotation and non-circular motions contribute equally to $\Gamma$. For our simulated galaxies, the gas dynamics at the scale of molecular clouds is usually dominated by noncircular motions, but in the center of galactic disks galactic rotation is still relevant. Our model shows that the transfer of rotation from large scales breaks at the scale $\lambda_c$ and this transition is necessary to reproduce the circulation distribution. We find that $\lambda_{\rm eq}$, and therefore the structure of the gas velocity field, is set by the local conditions of gravitational stability and stellar feedback.
[2]  oai:arXiv.org:2007.01877  [pdf] - 2128465
Ubiquitous velocity fluctuations throughout the molecular interstellar medium
Comments: Published in Nature Astronomy on July 6th 2020. This is the authors' version before final edits. Includes methods and supplementary information. Link to the NA publication: https://www.nature.com/articles/s41550-020-1126-z
Submitted: 2020-07-03
The density structure of the interstellar medium (ISM) determines where stars form and release energy, momentum, and heavy elements, driving galaxy evolution. Density variations are seeded and amplified by gas motion, but the exact nature of this motion is unknown across spatial scale and galactic environment. Although dense star-forming gas likely emerges from a combination of instabilities, convergent flows, and turbulence, establishing the precise origin is challenging because it requires quantifying gas motion over many orders of magnitude in spatial scale. Here we measure the motion of molecular gas in the Milky Way and in nearby galaxy NGC 4321, assembling observations that span an unprecedented spatial dynamic range ($10^{-1}{-}10^3$ pc). We detect ubiquitous velocity fluctuations across all spatial scales and galactic environments. Statistical analysis of these fluctuations indicates how star-forming gas is assembled. We discover oscillatory gas flows with wavelengths ranging from $0.3{-}400$ pc. These flows are coupled to regularly-spaced density enhancements that likely form via gravitational instabilities. We also identify stochastic and scale-free velocity and density fluctuations, consistent with the structure generated in turbulent flows. Our results demonstrate that ISM structure cannot be considered in isolation. Instead, its formation and evolution is controlled by nested, interdependent flows of matter covering many orders of magnitude in spatial scale.
[3]  oai:arXiv.org:2003.05384  [pdf] - 2119909
Dense Gas in a Giant Molecular Filament
Comments: 17 pages, 13 figures, accepted by A&A
Submitted: 2020-03-11, last modified: 2020-06-22
Recent surveys of the Galactic plane in the dust continuum and CO emission lines reveal that large ($\gtrsim 50$~pc) and massive ($\gtrsim 10^5$~$M_\odot$) filaments, know as giant molecular filaments (GMFs), may be linked to galactic dynamics and trace the mid-plane of the gravitational potential in the Milky Way. We have imaged one entire GMF located at $l\sim$52--54$^\circ$ longitude, GMF54 ($\sim$68~pc long), in the empirical dense gas tracers using the HCN(1--0), HNC(1--0), HCO$^+$(1--0) lines, and their $^{13}$C isotopologue transitions, as well as the N$_2$H$^+$(1--0) line. We study the dense gas distribution, the column density probability density functions (N-PDFs) and the line ratios within the GMF. The dense gas molecular transitions follow the extended structure of the filament with area filling factors between 0.06 and 0.28 with respect to $^{13}$CO(1--0). We constructed the N-PDFs of H$_2$ for each of the dense gas tracers based on their column densities and assumed uniform abundance. The N-PDFs of the dense gas tracers appear curved in log-log representation, and the HCO$^+$ N-PDF has the largest log-normal width and flattest power-law slope index. Studying the N-PDFs for sub-regions of GMF54, we found an evolutionary trend in the N-PDFs that high-mass star forming and Photon-Dominate Regions (PDRs) have flatter power-law indices. The integrated intensity ratios of the molecular lines in GMF54 are comparable to those in nearby galaxies. In particular, the N$_2$H$^+$/$^{13}$CO ratio, which traces the dense gas fraction, has similar values in GMF54 and all nearby galaxies except ULIRGs.
[4]  oai:arXiv.org:2005.11709  [pdf] - 2100641
PHANGS CO kinematics: disk orientations and rotation curves at 150 pc resolution
Comments: 19 figures, 36 pages, accepted for publication in ApJ. Table of PHANGS rotation curves available from http://phangs.org/data
Submitted: 2020-05-24
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.
[5]  oai:arXiv.org:2003.12576  [pdf] - 2077130
Reproducing the Universe: a comparison between the EAGLE simulations and the nearby DustPedia galaxy sample
Comments: 17 pages, 14 figures, accepted for publication in MNRAS
Submitted: 2020-03-27
We compare the spectral energy distributions (SEDs) and inferred physical properties for simulated and observed galaxies at low redshift. We exploit UV-submillimetre mock fluxes of ~7000 z=0 galaxies from the EAGLE suite of cosmological simulations, derived using the radiative transfer code SKIRT. We compare these to ~800 observed galaxies in the UV-submillimetre range, from the DustPedia sample of nearby galaxies. To derive global properties, we apply the SED fitting code CIGALE consistently to both data sets, using the same set of ~80 million models. The results of this comparison reveal overall agreement between the simulations and observations, both in the SEDs and in the derived physical properties, with a number of discrepancies. The optical and far-infrared regimes, and the scaling relations based upon the global emission, diffuse dust and stellar mass, show high levels of agreement. However, the mid-infrared fluxes of the EAGLE galaxies are overestimated while the far-UV domain is not attenuated enough, compared to the observations. We attribute these discrepancies to a combination of galaxy population differences between the samples, and limitations in the subgrid treatment of star-forming regions in the EAGLE-SKIRT post-processing recipe. Our findings show the importance of detailed radiative transfer calculations and consistent comparison, and provide suggestions for improved numerical models.
[6]  oai:arXiv.org:2002.08964  [pdf] - 2131871
Dynamical Equilibrium in the Molecular ISM in 28 Nearby Star-Forming Galaxies
Comments: 28 pages + 3 appendices, ApJ in press. See https://www.youtube.com/watch?v=qxkd-RXB0Ek for a short video describing the main results. Data tables available at https://www.canfar.net/storage/list/phangs/RELEASES/Sun_etal_2020 prior to publication
Submitted: 2020-02-20
We compare the observed turbulent pressure in molecular gas, $P_\mathrm{turb}$, to the required pressure for the interstellar gas to stay in equilibrium in the gravitational potential of a galaxy, $P_\mathrm{DE}$. To do this, we combine arcsecond resolution CO data from PHANGS-ALMA with multi-wavelength data that traces the atomic gas, stellar structure, and star formation rate (SFR) for 28 nearby star-forming galaxies. We find that $P_\mathrm{turb}$ correlates with, but almost always exceeds the estimated $P_\mathrm{DE}$ on kiloparsec scales. This indicates that the molecular gas is over-pressurized relative to the large-scale environment. We show that this over-pressurization can be explained by the clumpy nature of molecular gas; a revised estimate of $P_\mathrm{DE}$ on cloud scales, which accounts for molecular gas self-gravity, external gravity, and ambient pressure, agrees well with the observed $P_\mathrm{turb}$ in galaxy disks. We also find that molecular gas with cloud-scale ${P_\mathrm{turb}}\approx{P_\mathrm{DE}}\gtrsim{10^5\,k_\mathrm{B}\,\mathrm{K\,cm^{-3}}}$ in our sample is more likely to be self-gravitating, whereas gas at lower pressure appears more influenced by ambient pressure and/or external gravity. Furthermore, we show that the ratio between $P_\mathrm{turb}$ and the observed SFR surface density, $\Sigma_\mathrm{SFR}$, is compatible with stellar feedback-driven momentum injection in most cases, while a subset of the regions may show evidence of turbulence driven by additional sources. The correlation between $\Sigma_\mathrm{SFR}$ and kpc-scale $P_\mathrm{DE}$ in galaxy disks is consistent with the expectation from self-regulated star formation models. Finally, we confirm the empirical correlation between molecular-to-atomic gas ratio and kpc-scale $P_\mathrm{DE}$ reported in previous works.
[7]  oai:arXiv.org:2001.07459  [pdf] - 2076845
A model for the onset of self-gravitation and star formation in molecular gas governed by galactic forces: II. the bottleneck to collapse set by cloud-environment decoupling
Comments: Accepted for publication in ApJ, 32 pages, 11 figures
Submitted: 2020-01-21
In Meidt et al. (2018), we showed that gas kinematics on the scale of individual molecular clouds are not dominated by self-gravity but also track a component that originates with orbital motion in the potential of the host galaxy. This agrees with observed cloud line widths, which show systematic variations from virial motions with environment, pointing at the influence of the galaxy potential. In this paper, we hypothesize that these motions act to slow down the collapse of gas and so help regulate star formation. Extending the results of Meidt et al. (2018), we derive a dynamical collapse timescale that approaches the free-fall time only once the gas has fully decoupled from the galactic potential. Using this timescale we make predictions for how the fraction of free-falling, strongly self-gravitating gas varies throughout the disks of star-forming galaxies. We also use this collapse timescale to predict variations in the molecular gas star formation efficiency, which is lowered from a maximum, feedback-regulated level in the presence of strong coupling to the galactic potential. Our model implies that gas can only decouple from the galaxy to collapse and efficiently form stars deep within clouds. We show that this naturally explains the observed drop in star formation rate per unit gas mass in the Milky Way's CMZ and other galaxy centers. The model for a galactic bottleneck to star formation also agrees well with resolved observations of dense gas and star formation in galaxy disks and the properties of local clouds.
[8]  oai:arXiv.org:1911.08977  [pdf] - 2020603
Rotation of molecular clouds in M~51
Comments: Accepted by Astronomy and Astrophysics
Submitted: 2019-11-20
The grand-design spiral galaxy M~51 was observed at 40pc resolution in CO(1--0) by the PAWS project. A large number of molecular clouds were identified and we search for velocity gradients in two high signal-to-noise subsamples, containing 682 and 376 clouds. The velocity gradients are found to be systematically prograde oriented, as was previously found for the rather flocculent spiral M~33. This strongly supports the idea that the velocity gradients reflect cloud rotation, rather than more random dynamical forces, such as turbulence. Not only are the gradients prograde, but their $\frac{\partial v}{\partial x}$ and $\frac{\partial v}{\partial y}$ coefficients follow galactic shear in sign, although with a lower amplitude. No link is found between the orientation of the gradient and the orientation of the cloud. The values of the cloud angular momenta appear to be an extension of the values noted for galactic clouds despite the orders of magnitude difference in cloud mass. Roughly 30\% of the clouds show retrograde velocity gradients. For a strictly rising rotation curve, as in M~51, gravitational contraction would be expected to yield strictly prograde rotators within an axisymmetric potential. In M~51, the fraction of retrograde rotators is found to be higher in the spiral arms than in the disk as a whole. Along the leading edge of the spiral arms, a majority of the clouds are retrograde rotators. While this work should be continued on other nearby galaxies, the M~33 and M~51 studies have shown that clouds rotate and that they rotate mostly prograde, although the amplitudes are not such that rotational energy is a significant support mechanism against gravitation. In this work, we show that retrograde rotation is linked to the presence of a spiral gravitational potential.
[9]  oai:arXiv.org:1911.03479  [pdf] - 2026285
The lifecycle of molecular clouds in nearby star-forming disc galaxies
Comments: 39 pages, 14 figures; resubmitted to MNRAS after a favourable referee report (November 6, 2019)
Submitted: 2019-11-08
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.
[10]  oai:arXiv.org:1910.14311  [pdf] - 2054060
The headlight cloud in NGC 628: An extreme giant molecular cloud in a typical galaxy disk
Comments: A&A accepted
Submitted: 2019-10-31, last modified: 2019-11-04
Cloud-scale surveys of molecular gas reveal the link between molecular clouds properties and star formation (SF) across a range of galactic environments. Cloud populations in galaxy disks are considered to be representative of the `normal' SF. At high resolution, however, clouds with exceptional gas properties and SF activity may also be observed in normal disk environments. In this paper, we study the brightest cloud traced in CO emission in the disk of NGC628. The cloud is spatially coincident with an extremely bright HII region. We characterize its molecular gas properties and investigate how feedback and large-scale processes influence the properties of the molecular gas. High resolution CO ALMA observations are used to characterize its mass and dynamical state, which are compared to other clouds in NGC628. A LVG analysis is used to constrain the beam-diluted density and temperature of the molecular gas. We analyze the MUSE spectrum using Starburst99 to characterize the young stellar population associated with the HII region. The cloud is massive ($1-2\times10^7$M$_{\odot}$), with a beam-diluted density of $n_{\rm H_2}=5\times10^4$ cm$^{-3}$. It has a low virial parameter, suggesting that its CO emission may be overluminous due to heating by the HII region. A young ($2-4$ Myr), massive $3\times10^{5}$ M$_{\odot}$ stellar population is associated. We argue that the cloud is currently being destroyed by feedback from young massive stars. Due to the cloud's large mass, this phase of the cloud's evolution is long enough for the impact of feedback on the excitation of the gas to be observed. Its high mass may be related to its location at a spiral co-rotation radius, where gas experiences reduced galactic shear compared to other regions of the disk, and receives a sustained inflow of gas that can promote the cloud's mass growth.
[11]  oai:arXiv.org:1910.10520  [pdf] - 2026110
The Gas-Star Formation Cycle in Nearby Star-Forming Galaxies I. Assessment of Multi-scale Variations
Comments: Accepted to ApJ, 29 pages with 9 figures and 7 tables
Submitted: 2019-10-23
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.
[12]  oai:arXiv.org:1910.09815  [pdf] - 2026087
Mapping electron temperature variations across a spiral arm in NGC 1672
Comments: 9 pages, 5 figure, accepted to ApJ Letter
Submitted: 2019-10-22
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.
[13]  oai:arXiv.org:1910.07190  [pdf] - 2026025
Mapping metallicity variations across nearby galaxy disks
Comments: 14 pages, 9 figures, 4 Appendices including image atlases, ApJ accepted
Submitted: 2019-10-16
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.
[14]  oai:arXiv.org:1902.10437  [pdf] - 1875340
Dense gas is not enough: environmental variations in the star formation efficiency of dense molecular gas at 100pc scales in M51
Comments: 23 pages, 9 figures, accepted for publication in A&A
Submitted: 2019-02-27, last modified: 2019-04-09
It remains unclear what sets the efficiency with which molecular gas transforms into stars. Here we present a new VLA map of the spiral galaxy M51 in 33GHz radio continuum, an extinction-free tracer of star formation, at 3" scales (~100pc). We combined this map with interferometric PdBI/NOEMA observations of CO(1-0) and HCN(1-0) at matched resolution for three regions in M51 (central molecular ring, northern and southern spiral arm segments). While our measurements roughly fall on the well-known correlation between total infrared and HCN luminosity, bridging the gap between Galactic and extragalactic observations, we find systematic offsets from that relation for different dynamical environments probed in M51, e.g. the southern arm segment is more quiescent due to low star formation efficiency (SFE) of the dense gas, despite having a high dense gas fraction. Combining our results with measurements from the literature at 100pc scales, we find that the SFE of the dense gas and the dense gas fraction anti-correlate and correlate, respectively, with the local stellar mass surface density. This is consistent with previous kpc-scale studies. In addition, we find a significant anti-correlation between the SFE and velocity dispersion of the dense gas. Finally, we confirm that a correlation also holds between star formation rate surface density and the dense gas fraction, but it is not stronger than the correlation with dense gas surface density. Our results are hard to reconcile with models relying on a universal gas density threshold for star formation and suggest that turbulence and galactic dynamics play a major role in setting how efficiently dense gas converts into stars.
[15]  oai:arXiv.org:1812.06109  [pdf] - 1806286
The Spatial Relation between Young Star Clusters and Molecular Clouds in M 51 with LEGUS
Comments: 17 pages, 12 figures. Accepted for publication in MNRAS
Submitted: 2018-12-14
We present a study correlating the spatial locations of young star clusters with those of molecular clouds in NGC~5194, in order to investigate the timescale over which clusters separate from their birth clouds. The star cluster catalogues are from the Legacy ExtraGalactic UV Survey (LEGUS) and the molecular clouds from the Plateau de Bure Interefrometer Arcsecond Whirpool Survey (PAWS). We find that younger star clusters are spatially closer to molecular clouds than older star clusters. The median ages for clusters associated with clouds is 4~Myr whereas it is 50~Myr for clusters that are sufficiently separated from a molecular cloud to be considered unassociated. After $\sim$6~Myr, the majority of the star clusters lose association with their molecular gas. Younger star clusters are also preferentially located in stellar spiral arms where they are hierarchically distributed in kpc-size regions for 50-100~Myr before dispersing. The youngest star clusters are more strongly clustered, yielding a two-point correlation function with $\alpha=-0.28\pm0.04$, than the GMCs ($\alpha=-0.09\pm0.03$) within the same PAWS field. However, the clustering strength of the most massive GMCs, supposedly the progenitors of the young clusters for a star formation efficiency of a few percent, is comparable ($\alpha=-0.35\pm0.05$) to that of the clusters. We find a galactocentric-dependence for the coherence of star formation, in which clusters located in the inner region of the galaxy reside in smaller star-forming complexes and display more homogeneous distributions than clusters further from the centre. This result suggests a correlation between the survival of a cluster complex and its environment.
[16]  oai:arXiv.org:1811.07940  [pdf] - 1811219
Do Spectroscopic Dense Gas Fractions Track Molecular Cloud Surface Densities?
Comments: 11 pages, 4 figures, accepted for publication in The Astrophysical Journal Letters
Submitted: 2018-11-19
We use ALMA and IRAM 30-m telescope data to investigate the relationship between the spectroscopically-traced dense gas fraction and the cloud-scale (120 pc) molecular gas surface density in five nearby, star-forming galaxies. We estimate the dense gas mass fraction at 650 pc and 2800 pc scales using the ratio of HCN (1-0) to CO (1-0) emission. We then use high resolution (120 pc) CO (2-1) maps to calculate the mass-weighted average molecular gas surface density within 650 pc or 2770 pc beam where the dense gas fraction is estimated. On average, the dense gas fraction correlates with the mass-weighted average molecular gas surface density. Thus, parts of a galaxy with higher mean cloud-scale gas surface density also appear to have a larger fraction of dense gas. The normalization and slope of the correlation do vary from galaxy to galaxy and with the size of the regions studied. This correlation is consistent with a scenario where the large-scale environment sets the gas volume density distribution, and this distribution manifests in both the cloud-scale surface density and the dense gas mass fraction.
[17]  oai:arXiv.org:1807.02043  [pdf] - 1767504
Azimuthal variations of gas-phase oxygen abundance in NGC 2997
Comments: 13 pages, 17 figures, accepted to A&A
Submitted: 2018-07-05
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.
[18]  oai:arXiv.org:1805.00937  [pdf] - 2131835
Cloud-Scale Molecular Gas Properties in 15 Nearby Galaxies
Comments: Accepted for publication in ApJ. 45 pages, 11 figures, 8 tables, 4 Appendices; key results summarized in Figure 10. Machine-readable table can be downloaded at http://www.astronomy.ohio-state.edu/~sun.1608/datafile3.txt prior to publication. For a brief video describing the main results of this paper, please see https://www.youtube.com/watch?v=-_eL7t1PVq8&t
Submitted: 2018-05-02
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.
[19]  oai:arXiv.org:1803.10785  [pdf] - 1686808
Dense Gas, Dynamical Equilibrium Pressure, and Star Formation in Nearby Star-Forming Galaxies
Comments: 31 pages, 13 figures, accepted for publication in The Astrophysical Journal, email for access to data table before publication
Submitted: 2018-03-28
We use new ALMA observations to investigate the connection between dense gas fraction, star formation rate, and local environment across the inner region of four local galaxies showing a wide range of molecular gas depletion times. We map HCN (1-0), HCO$^+$ (1-0), CS (2-1), $^{13}$CO (1-0), and C$^{18}$O (1-0) across the inner few kpc of each target. We combine these data with short spacing information from the IRAM large program EMPIRE, archival CO maps, tracers of stellar structure and recent star formation, and recent HCN surveys by Bigiel et al. and Usero et al. We test the degree to which changes in the dense gas fraction drive changes in the SFR. $I_{HCN}/I_{CO}$ (tracing the dense gas fraction) correlates strongly with $I_{CO}$ (tracing molecular gas surface density), stellar surface density, and dynamical equilibrium pressure, $P_{DE}$. Therefore, $I_{HCN}/I_{CO}$ becomes very low and HCN becomes very faint at large galactocentric radii, where ratios as low as $I_{HCN}/I_{CO} \sim 0.01$ become common. The apparent ability of dense gas to form stars, $\Sigma_{SFR}/\Sigma_{dense}$ (where $\Sigma_{dense}$ is traced by the HCN intensity and the star formation rate is traced by a combination of H$\alpha$ and 24$\mu$m emission), also depends on environment. $\Sigma_{SFR}/\Sigma_{dense}$ decreases in regions of high gas surface density, high stellar surface density, and high $P_{DE}$. Statistically, these correlations between environment and both $\Sigma_{SFR}/\Sigma_{dense}$ and $I_{HCN}/I_{CO}$ are stronger than that between apparent dense gas fraction ($I_{HCN}/I_{CO}$) and the apparent molecular gas star formation efficiency $\Sigma_{SFR}/\Sigma_{mol}$. We show that these results are not specific to HCN.
[20]  oai:arXiv.org:1712.06364  [pdf] - 1641377
A Model for the Onset of Self-gravitation and Star Formation in Molecular Gas Governed by Galactic Forces: I. Cloud-scale Gas Motions
Comments: Accepted for publication in ApJ, 26 pages, 11 figures
Submitted: 2017-12-18
Modern extragalactic molecular gas surveys now reach the scales of star-forming giant molecular clouds (GMCs, 20-50 pc). Systematic variations in GMC properties with galaxy environment imply that clouds are not universally self-gravitating objects, decoupled from their surroundings. Here we reexamine the coupling of clouds to their environment and develop a model for 3D gas motions generated by forces arising with the galaxy gravitational potential defined by the background disk of stars and dark matter. We show that these motions can resemble or even exceed the motions needed to support gas against its own self-gravity throughout typical galaxy disks. The importance of the galactic potential in spiral arms and galaxy centers suggests that the response to self-gravity does not always dominate the motions of gas at GMC scales, with implications for observed gas kinematics, virial equilibrium and cloud morphology. We describe how a uniform treatment of gas motions in the plane and in the vertical direction synthesizes the two main mechanisms proposed to regulate star formation: vertical pressure equilibrium and shear/Coriolis forces as parameterized by Toomre Q~1. As the modeled motions are coherent and continually driven by the external potential, they represent support for the gas that is distinct from that conventionally attributed to turbulence, which decays rapidly and requires thus maintenance, e.g. via feedback from star formation. Thus our model suggests the galaxy itself can impose an important limit to star formation, as we explore in a second paper in this series.
[21]  oai:arXiv.org:1709.06649  [pdf] - 1593628
Orbital decomposition of CALIFA spiral galaxies
Comments: accepted by MNRAS
Submitted: 2017-09-19
Schwarzschild orbit-based dynamical models are widely used to uncover the internal dynamics of early-type galaxies and globular clusters. Here we present for the first time the Schwarzschild models of late-type galaxies: an SBb galaxy NGC 4210 and an S0 galaxy NGC 6278 from the CALIFA survey. The mass profiles within $2\,R_e$ are constrained well with $1\sigma$ statistical error of $\sim 10\%$. The luminous and dark mass can be disentangled with uncertainties of $\sim 20\%$ and $\sim 50\%$ respectively. From $R_e$ to $2\,R_e$, the dark matter fraction increases from $14\pm10\%$ to $18\pm10\%$ for NGC 4210 and from $15\pm10\%$ to $30\pm20\%$ for NGC 6278. The velocity anisotropy profiles of both $\sigma_r/\sigma_t$ and $\sigma_z/\sigma_R$ are well constrained. The inferred internal orbital distributions reveal clear substructures. The orbits are naturally separated into three components: a cold component with near circular orbits, a hot component with near radial orbits, and a warm component in between. The photometrically-identified exponential disks are predominantly made up of cold orbits only beyond $\sim 1\,R_e$, while they are constructed mainly with the warm orbits inside. Our dynamical hot components are concentrated in the inner regions, similar to the photometrically-identified bulges. The reliability of the results, especially the orbit distribution, are verified by applying the model to mock data.
[22]  oai:arXiv.org:1708.04966  [pdf] - 1587172
The Chemical Evolution Carousel of Spiral Galaxies: Azimuthal Variations of Oxygen Abundance in NGC1365
Comments: 19 pages, 13 figures. Accepted to ApJ
Submitted: 2017-08-16
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.
[23]  oai:arXiv.org:1706.08540  [pdf] - 1585161
Cloud Scale ISM Structure and Star Formation in M51
Comments: Accepted for publication in the Astrophysical Journal, 22 pages, 11 figures, email for access to data table before publication, for high resolution version before publication: http://www.astronomy.ohio-state.edu/~leroy.42/m51_cloudsf.pdf
Submitted: 2017-06-26
We compare the structure of molecular gas at $40$ pc resolution to the ability of gas to form stars across the disk of the spiral galaxy M51. We break the PAWS survey into $370$ pc and $1.1$ kpc resolution elements, and within each we estimate the molecular gas depletion time ($\tau_{\rm Dep}^{\rm mol}$), the star formation efficiency per free fall time ($\epsilon_{\rm ff}$), and the mass-weighted cloud-scale (40 pc) properties of the molecular gas: surface density, $\Sigma$, line width, $\sigma$, and $b\equiv\Sigma/\sigma^2\propto\alpha_{\rm vir}^{-1}$, a parameter that traces the boundedness of the gas. We show that the cloud-scale surface density appears to be a reasonable proxy for mean volume density. Applying this, we find a typical star formation efficiency per free-fall time, $\epsilon_{ff} \left( \left< \Sigma_{40pc} \right> \right) \sim 0.3{-}0.36\%$, lower than adopted in many models and found for local clouds. More, the efficiency per free fall time anti-correlates with both $\Sigma$ and $\sigma$, in some tension with turbulent star formation models. The best predictor of the rate of star formation per unit gas mass in our analysis is $b \equiv \Sigma / \sigma^2$, tracing the strength of self gravity, with $\tau_{\rm Dep}^{\rm mol} \propto b^{-0.9}$. The sense of the correlation is that gas with stronger self-gravity (higher $b$) forms stars at a higher rate (low $\tau_{\rm Dep}^{\rm mol}$). The different regions of the galaxy mostly overlap in $\tau_{\rm Dep}^{\rm mol}$ as a function of $b$, so that low $b$ explains the surprisingly high $\tau_{\rm Dep}^{\rm mol}$ found towards the inner spiral arms found by by Meidt et al. (2013).
[24]  oai:arXiv.org:1701.02184  [pdf] - 1534297
The PdBI Arcsecond Whirlpool Survey (PAWS). The Role of Spiral Arms in Cloud and Star Formation
Comments: 52 pages, 9 figures, 8 tables, accepted to ApJ
Submitted: 2017-01-09
The process that leads to the formation of the bright star forming sites observed along prominent spiral arms remains elusive. We present results of a multi-wavelength study of a spiral arm segment in the nearby grand-design spiral galaxy M51 that belongs to a spiral density wave and exhibits nine gas spurs. The combined observations of the(ionized, atomic, molecular, dusty) interstellar medium (ISM) with star formation tracers (HII regions, young <10Myr stellar clusters) suggest (1) no variation in giant molecular cloud (GMC) properties between arm and gas spurs, (2) gas spurs and extinction feathers arising from the same structure with a close spatial relation between gas spurs and ongoing/recent star formation (despite higher gas surface densities in the spiral arm), (3) no trend in star formation age either along the arm or along a spur, (4) evidence for strong star formation feedback in gas spurs: (5) tentative evidence for star formation triggered by stellar feedback for one spur, and (6) GMC associations (GMAs) being no special entities but the result of blending of gas arm/spur cross-sections in lower resolution observations. We conclude that there is no evidence for a coherent star formation onset mechanism that can be solely associated to the presence of the spiral density wave. This suggests that other (more localized) mechanisms are important to delay star formation such that it occurs in spurs. The evidence of star formation proceeding over several million years within individual spurs implies that the mechanism that leads to star formation acts or is sustained over a longer time-scale.
[25]  oai:arXiv.org:1611.01705  [pdf] - 1542784
The Radio Spectral Energy Distribution and Star Formation Rate Calibration in Galaxies
Comments: 29 pages, 15 figures, accepted for publication in Astrophysical Journal
Submitted: 2016-11-05, last modified: 2016-12-28
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.
[26]  oai:arXiv.org:1609.06448  [pdf] - 1531679
Galactic Bar/Spiral Arm Interactions in NGC3627
Comments: 14 pages, 8 figures, accepted for Astronomy & Astrophysics, a high-resolution version can be found at http://www.mpia.de/homes/beuther/papers.html
Submitted: 2016-09-21
Aims: To gain insight into the expected gas dynamics at the interface of the Galactic bar and spiral arms in our own Milky Way galaxy, we examine as an extragalactic counterpart the evidence for multiple distinct velocity components in the cold, dense molecular gas populating a comparable region at the end of the bar in the nearby galaxy NGC3627. Methods: We assemble a high resolution view of molecular gas kinematics traced by CO(2-1) emission and extract line-of-sight velocity profiles from regions of high and low gas velocity dispersion. Results: The high velocity dispersions arise with often double-peaked or multiple line-profiles. We compare the centroids of the different velocity components to expectations based on orbital dynamics in the presence of bar and spiral potential perturbations. A model of the region as the interface of two gas-populated orbits families supporting the bar and the independently rotating spiral arms provides an overall good match to the data. An extent of the bar to the corotation radius of the galaxy is favored. Conclusions: Using NGC3627 as an extragalactic example, we expect situations like this to favor strong star formation events such as observed in our own Milky Way since gas can pile up at the crossings between the orbit families. The relative motions of the material following these orbits is likely even more important for the build up of high density in the region. The surface densities in NGC3627 are also so high that shear at the bar end is unlikely to significantly weaken the star formation activity. We speculate that scenarios in which the bar and spiral rotate at two different pattern speeds may be the most favorable for intense star formation at such interfaces.
[27]  oai:arXiv.org:1608.07584  [pdf] - 1531368
Being WISE II: Reducing the Influence of Star Formation History on the Mass-to-Light Ratio of Quiescent Galaxies
Comments: 15 pages, 11 figures. Accepted for publication in ApJ
Submitted: 2016-08-26
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$\%$.
[28]  oai:arXiv.org:1607.00010  [pdf] - 1490763
AGN feedback in the nucleus of M51
Comments: 21 pages, 17 figures, accepted for publication in A&A
Submitted: 2016-06-30
AGN feedback is invoked as one of the most relevant mechanisms that shape the evolution of galaxies. Our goal is to understand the interplay between AGN feedback and the interstellar medium in M51, a nearby spiral galaxy with a modest AGN and a kpc-scale radio jet expanding through the disc of the galaxy. For that purpose, we combine molecular gas observations in the CO(1-0) and HCN(1-0) lines from the Plateau de Bure interferometer with archival radio, X-ray, and optical data. We show that there is a significant scarcity of CO emission in the ionisation cone, while molecular gas emission tends to accumulate towards the edges of the cone. The distribution and kinematics of CO and HCN line emission reveal AGN feedback effects out to r~500pc, covering the whole extent of the radio jet, with complex kinematics in the molecular gas which displays strong local variations. We propose that this is the result of the almost coplanar jet pushing on molecular gas in different directions as it expands; the effects are more pronounced in HCN than in CO emission, probably as the result of radiative shocks. Following previous interpretation of the redshifted molecular line in the central 5" as caused by a molecular outflow, we estimate the outflow rates to be Mdot_H2~0.9Msun/yr and Mdot_dense~0.6Msun/yr, which are comparable to the molecular inflow rates (~1Msun/yr); gas inflow and AGN feedback could be mutually regulated processes. The agreement with findings in other nearby radio galaxies suggests that this is not an isolated case, and probably the paradigm of AGN feedback through radio jets, at least for galaxies hosting low-luminosity active nuclei.
[29]  oai:arXiv.org:1606.07077  [pdf] - 1510269
A Portrait of Cold Gas in Galaxies at 60pc Resolution and a Simple Method to Test Hypotheses That Link Small-Scale ISM Structure to Galaxy-Scale Processes
Comments: Accepted for publication in the Astrophysical Journal, 36 pages (24+appendix), 21 figures (12+appendix), until publication high resolution version at http://www.astronomy.ohio-state.edu/~leroy.42/cloudscale.pdf
Submitted: 2016-06-22
The cloud-scale density, velocity dispersion, and gravitational boundedness of the interstellar medium (ISM) vary within and among galaxies. In turbulent models, these properties play key roles in the ability of gas to form stars. New high fidelity, high resolution surveys offer the prospect to measure these quantities across galaxies. We present a simple approach to make such measurements and to test hypotheses that link small-scale gas structure to star formation and galactic environment. Our calculations capture the key physics of the Larson scaling relations, and we show good correspondence between our approach and a traditional "cloud properties" treatment. However, we argue that our method is preferable in many cases because of its simple, reproducible characterization of all emission. Using, low-J 12CO data from recent surveys, we characterize the molecular ISM at 60pc resolution in the Antennae, the Large Magellanic Cloud, M31, M33, M51, and M74. We report the distributions of surface density, velocity dispersion, and gravitational boundedness at 60pc scales and show galaxy-to-galaxy and intra-galaxy variations in each. The distribution of flux as a function of surface density appears roughly lognormal with a 1sigma width of ~0.3 dex, though the center of this distribution varies from galaxy to galaxy. The 60pc resolution line width and molecular gas surface density correlate well, which is a fundamental behavior expected for virialized or free-falling gas. Varying the measurement scale for the LMC and M31, we show that the molecular ISM has higher surface densities, lower line widths, and more self-gravity at smaller scales.
[30]  oai:arXiv.org:1603.08009  [pdf] - 1457209
Characterizing spiral arm and interarm star formation
Comments: 10 pages, 4 figures, 1 table, accepted for publication in ApJ
Submitted: 2016-03-25, last modified: 2016-06-14
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.
[31]  oai:arXiv.org:1605.03991  [pdf] - 1457295
First survey of Wolf-Rayet star populations over the full extension of nearby galaxies observed with CALIFA
Comments: 33 pages, accepted for publication in A&A
Submitted: 2016-05-12
The search of extragalactic regions with conspicuous presence of Wolf-Rayet (WR) stars outside the Local Group is challenging task due to the difficulties in detecting their faint spectral features. In this exploratory work, we develop a methodology to perform an automated search of WR signatures through a pixel-by-pixel analysis of integral field spectroscopy (IFS) data belonging to the Calar Alto Legacy Integral Field Area survey, CALIFA. This technique allowed us to build the first catalogue of Wolf-Rayet rich regions with spatially-resolved information, allowing to study the properties of these complexes in a 2D context. The detection technique is based on the identification of the blue WR bump (around He II 4686 {\AA}, mainly associated to nitrogen-rich WR stars, WN) and the red WR bump (around C IV 5808 {\AA} and associated to carbon-rich WR stars, WC) using a pixel-by-pixel analysis. We identified 44 WR-rich regions with blue bumps distributed in 25 galaxies of a total of 558. The red WR bump was identified only in 5 of those regions. We found that the majority of the galaxies hosting WR populations in our sample are involved in some kind of interaction process. Half of the host galaxies share some properties with gamma-ray burst (GRB) hosts where WR stars, as potential candidates to being the progenitors of GRBs, are found. We also compared the WR properties derived from the CALIFA data with stellar population synthesis models, and confirm that simple star models are generally not able to reproduce the observations. We conclude that other effects, such as the binary star channel (which could extend the WR phase up to 10 Myr), fast rotation or other physical processes that causes the loss of observed Lyman continuum photons, are very likely affecting the derived WR properties, and hence should be considered when modelling the evolution of massive stars.
[32]  oai:arXiv.org:1604.08777  [pdf] - 1457258
Comparing [CII], HI, and CO dynamics of nearby galaxies
Comments: Accepted for publication in the Astronomical Journal
Submitted: 2016-04-29
The HI and CO components of the interstellar medium (ISM) are usually used to derive the dynamical mass M_dyn of nearby galaxies. Both components become too faint to be used as a tracer in observations of high-redshift galaxies. In those cases, the 158 $\mu$m line of atomic carbon [CII] may be the only way to derive M_dyn. As the distribution and kinematics of the ISM tracer affects the determination of M_dyn, it is important to quantify the relative distributions of HI, CO and [CII]. HI and CO are well-characterised observationally, however, for [CII] only very few measurements exist. Here we compare observations of CO, HI, and [CII] emission of a sample of nearby galaxies, drawn from the HERACLES, THINGS and KINGFISH surveys. We find that within R_25, the average [CII] exponential radial profile is slightly shallower than that of the CO, but much steeper than the HI distribution. This is also reflected in the integrated spectrum ("global profile"), where the [CII] spectrum looks more like that of the CO than that of the HI. For one galaxy, a spectrally resolved comparison of integrated spectra was possible; other comparisons were limited by the intrinsic line-widths of the galaxies and the coarse velocity resolution of the [CII] data. Using high-spectral-resolution SOFIA [CII] data of a number of star forming regions in two nearby galaxies, we find that their [CII] linewidths agree better with those of the CO than the HI. As the radial extent of a given ISM tracer is a key input in deriving M_dyn from spatially unresolved data, we conclude that the relevant length-scale to use in determining M_dyn based on [CII] data, is that of the well-characterised CO distribution. This length scale is similar to that of the optical disk.
[33]  oai:arXiv.org:1603.07808  [pdf] - 1636578
Nearby supernova host galaxies from the CALIFA Survey: II. SN environmental metallicity
Comments: 24 pages, 16 Figures, 13 Tables, Accepted in A&A
Submitted: 2016-03-24, last modified: 2016-04-15
The metallicity of a supernova (SN) progenitor, together with its mass, is one of the main parameters that rules their outcome. We present a metallicity study of 115 nearby SN host galaxies (0.005<z<0.03) which hosted 142 SNe using Integral Field Spectroscopy (IFS) from the CALIFA survey. Using O3N2 we found no statistically significant differences between the gas-phase metallicities at the locations of the three main SN types (Ia, Ib/c and II) all having ~8.50$\pm$0.02 dex. The total galaxy metallicities are also very similar and we argue that this is because our sample consists only of SNe discovered in massive galaxies (log(M/Msun)>10 dex) by targeted searches. We also found no evidence that the metallicity at the SN location differs from the average metallicity at the GCD of the SNe. By extending our SN sample with published metallicities at the SN location, we studied the metallicity distributions for all SN subtypes split into SN discovered in targeted and untargeted searches. We confirm a bias toward higher host masses and metallicities in the targeted searches. Combining data from targeted and untargeted searches we found a sequence from higher to lower local metallicity: SN Ia, Ic, and II show the highest metallicity, which is significantly higher than SN Ib, IIb, and Ic-BL. Our results support the picture of SN Ib resulting from binary progenitors and, at least part of, SN Ic being the result of single massive stars stripped of their outer layers by metallicity driven winds. We studied several proxies of the local metallicity frequently used in the literature and found that the total host metallicity allows for the estimation of the metallicity at the SN location with an accuracy better than 0.08 dex and very small bias. In addition, weak AGNs not seen in total spectra may only weakly bias (by 0.04 dex) the metallicity estimate from integrated spectra. (abridged)
[34]  oai:arXiv.org:1510.03440  [pdf] - 1374125
Gravitational torques imply molecular gas inflow towards the nucleus of M51
Comments: 19 pages, 12 figures. Accepted by A&A
Submitted: 2015-10-12, last modified: 2016-02-01
The transport of gas towards the centre of galaxies is critical for black hole feeding and, indirectly, it can control active galactic nucleus (AGN) feedback. We have quantified the molecular gas inflow in the central R<1kpc of M51 to be 1 Msun/yr, using a new gravitational torque map and the molecular gas traced by the PdBI Arcsecond Whirlpool Survey (PAWS). The nuclear stellar bar is responsible for this gas inflow. We have also used torque profiles to estimate the location of dynamical resonances, suggesting a corotation for the bar at R~20", and a corotation for the spiral at R~100". We demonstrate how important it is to correct 3.6um images for dust emission in order to compute gravitational torques, and we carefully examine further sources of uncertainty. Our observational measurement of gas inflow can be compared with nuclear molecular outflow rates and provide useful constraints for numerical simulations.
[35]  oai:arXiv.org:1512.03260  [pdf] - 1359050
How galactic environment regulates star formation
Comments: 18 pages, 7 figures, accepted for publication in ApJ
Submitted: 2015-12-10
In a new simple model I reconcile two contradictory views on the factors that determine the rate at which molecular clouds form stars -- internal structure vs. external, environmental influences -- providing a unified picture for the regulation of star formation in galaxies. In the presence of external pressure, the pressure gradient set up within a self-gravitating isothermal cloud leads to a non-uniform density distribution. Thus the local environment of a cloud influences its internal structure. In the simple equilibrium model, the fraction of gas at high density in the cloud interior is determined simply by the cloud surface density, which is itself inherited from the pressure in the immediate surroundings. This idea is tested using measurements of the properties of local clouds, which are found to show remarkable agreement with the simple equilibrium model. The model also naturally predicts the star formation relation observed on cloud scales and, at the same time, provides a mapping between this relation and the closer-to-linear molecular star formation relation measured on larger scales in galaxies. The key is that pressure regulates not only the molecular content of the ISM but also the cloud surface density. I provide a straightforward prescription for the pressure regulation of star formation that can be directly implemented in numerical models. Predictions for the dense gas fraction and star formation efficiency measured on large-scales within galaxies are also presented, establishing the basis for a new picture of star formation regulated by galactic environment.
[36]  oai:arXiv.org:1410.0009  [pdf] - 1260746
The Spitzer Survey of Stellar Structure in Galaxies (S$^4$G): Precise Stellar Mass Distributions from Automated Dust Correction at 3.6 microns
Comments: Published by the Astrophysical Journal (Supplements). Data available as a tar-ball from IRSA (http://irsa.ipac.caltech.edu/data/SPITZER/S4G/S4Gp5.tar.gz); searchable option as part of S4G will become available in November 2014. Readme file and table of parameters inside tar-ball (P5_README.txt, P5_table.txt). 2015ApJS..219....5Q
Submitted: 2014-09-30, last modified: 2015-08-11
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.
[37]  oai:arXiv.org:1506.05153  [pdf] - 1259028
The Odd Offset between the Galactic Disk and Its Bar in NGC 3906
Comments: 18 pages, 8 figures, ApJ in press
Submitted: 2015-06-16
We use mid-infrared 3.6 and 4.5microns imaging of NGC 3906 from the Spitzer Survey of Stellar Structure in Galaxies (S4G) to understand the nature of an unusual offset between its stellar bar and the photometric center of an otherwise regular, circular outer stellar disk. We measure an offset of ~720 pc between the center of the stellar bar and photometric center of the stellar disk; the bar center coincides with the kinematic center of the disk determined from previous HI observations. Although the undisturbed shape of the disk suggests that NGC 3906 has not undergone a significant merger event in its recent history, the most plausible explanation for the observed offset is an interaction. Given the relatively isolated nature of NGC 3906 this interaction could be with dark matter sub structure in the galaxy's halo or from a recent interaction with a fast moving neighbor which remains to be identified. Simulations aimed at reproducing the observed offset between the stellar bar / kinematic center of the system and the photometric center of the disk are necessary to confirm this hypothesis and constrain the interaction history of the galaxy.
[38]  oai:arXiv.org:1505.03534  [pdf] - 1021972
The Spitzer Survey of Stellar Structure in Galaxies (S4G): Stellar Masses, Sizes and Radial Profiles for 2352 Nearby Galaxies
Comments: Accepted for publication in ApJS
Submitted: 2015-05-13
The Spitzer Survey of Stellar Structure in Galaxies (S4G) is a volume, magnitude, and size-limited survey of 2352 nearby galaxies with deep imaging at 3.6 and 4.5um. In this paper we describe our surface photometry pipeline and showcase the associated data products that we have released to the community. We also identify the physical mechanisms leading to different levels of central stellar mass concentration for galaxies with the same total stellar mass. Finally, we derive the local stellar mass-size relation at 3.6um for galaxies of different morphologies. Our radial profiles reach stellar mass surface densities below 1 Msun pc-2. Given the negligible impact of dust and the almost constant mass-to-light ratio at these wavelengths, these profiles constitute an accurate inventory of the radial distribution of stellar mass in nearby galaxies. From these profiles we have also derived global properties such as asymptotic magnitudes (and the corresponding stellar masses), isophotal sizes and shapes, and concentration indices. These and other data products from our various pipelines (science-ready mosaics, object masks, 2D image decompositions, and stellar mass maps), can be publicly accessed at IRSA (http://irsa.ipac.caltech.edu/data/SPITZER/S4G/).
[39]  oai:arXiv.org:1504.06282  [pdf] - 1258838
H-alpha kinematics of S4G spiral galaxies-II. Data description and non-circular motions
Comments: 23 pages, 14 figures, without appendices, accepted to be published in MNRAS
Submitted: 2015-04-23, last modified: 2015-04-24
We present a kinematical study of 29 spiral galaxies included in the Spitzer Survey of Stellar Structure in Galaxies, using Halpha Fabry-Perot data obtained with the Galaxy Halpha Fabry-Perot System instrument at the William Herschel Telescope in La Palma, complemented with images in the R-band and in Halpha. The primary goal is to study the evolution and properties of the main structural components of galaxies through the kinematical analysis of the FP data, complemented with studies of morphology, star formation and mass distribution. In this paper we describe how the FP data have been obtained, processed and analysed. We present the resulting moment maps, rotation curves, velocity model maps and residual maps. Images are available in FITS format through the NASA/IPAC Extragalactic Database and the Centre de Donn\'ees Stellaires. With these data products we study the non-circular motions, in particular those found along the bars and spiral arms. The data indicate that the amplitude of the non-circular motions created by the bar does not correlate with the bar strength indicators. The amplitude of those non-circular motions in the spiral arms does not correlate with either arm class or star formation rate along the spiral arms. This implies that the presence and the magnitude of the streaming motions in the arms is a local phenomenon.
[40]  oai:arXiv.org:1504.04528  [pdf] - 1182962
Short GMC lifetimes: an observational estimate with the PdBI Arcsecond Whirlpool Survey (PAWS)
Comments: 16 pages, 8 figures, accepted for publication in ApJ
Submitted: 2015-04-17
We describe and execute a novel approach to observationally estimate the lifetimes of giant molecular clouds (GMCs). We focus on the cloud population between the two main spiral arms in M51 (the inter-arm region) where cloud destruction via shear and star formation feedback dominates over formation processes. By monitoring the change in GMC number densities and properties from one side of the inter-arm to the other, we estimate the lifetime as a fraction of the inter-arm travel time. We find that GMC lifetimes in M51's inter-arm are finite and short, 20 to 30 Myr. Such short lifetimes suggest that cloud evolution is influenced by environment, in which processes can disrupt GMCs after a few free-fall times. Over most of the region under investigation shear appears to regulate the lifetime. As the shear timescale increases with galactocentric radius, we expect cloud destruction to switch primarily to star formation feedback at larger radii. We identify a transition from shear- to feedback-dominated disruption through a change in the behavior of the GMC number density. The signature suggests that shear is more efficient at completely dispersing clouds, whereas feedback transforms the population, e.g. by fragmenting high mass clouds into lower mass pieces. Compared to the characteristic timescale for molecular hydrogen in M51, our short lifetimes suggest that gas can remain molecular while clouds disperse and reassemble. We propose that galaxy dynamics regulates the cycling of molecular material from diffuse to bound (and ultimately star-forming) objects, contributing to long observed molecular depletion times in normal disk galaxies. We also speculate that, in more extreme environments such as elliptical galaxies and concentrated galaxy centers, star formation can be suppressed when the shear timescale becomes so short that some clouds can not survive to collapse and form stars.
[41]  oai:arXiv.org:1503.06550  [pdf] - 969451
Spitzer Survey of Stellar Structure in Galaxies (S$^4$G). The Pipeline 4: Multi-component decomposition strategies and data release
Comments: Accepted to ApJS
Submitted: 2015-03-23
The Spitzer Survey of Stellar Structure in Galaxies (S$^4$G, Sheth et. al. 2010) is a deep 3.6 and 4.5 $\mu$m imaging survey of 2352 nearby ($< 40$ Mpc) galaxies. We describe the S$^4$G data analysis pipeline 4, which is dedicated to 2-dimensional structural surface brightness decompositions of 3.6 $\mu$m images, using GALFIT3.0 \citep{peng2010}. Besides automatic 1-component S\'ersic fits, and 2-component S\'ersic bulge + exponential disk fits, we present human supervised multi-component decompositions, which include, when judged appropriate, a central point source, bulge, disk, and bar components. Comparison of the fitted parameters indicates that multi-component models are needed to obtain reliable estimates for the bulge S\'ersic index and bulge-to-total light ratio ($B/T$), confirming earlier results \citep{laurikainen2007, gadotti2008, weinzirl2009}. In this first paper, we describe the preparations of input data done for decompositions, give examples of our decomposition strategy, and describe the data products released via IRSA and via our web page ({\bf \tt www.oulu.fi/astronomy/S4G\_PIPELINE4/MAIN}). These products include all the input data and decomposition files in electronic form, making it easy to extend the decompositions to suit specific science purposes. We also provide our IDL-based visualization tools (GALFIDL) developed for displaying/running GALFIT-decompositions, as well as our mask editing procedure (MASK\_EDIT) used in data preparation. In the second paper we will present a detailed analysis of the bulge, disk, and bar parameter derived from multi-component decompositions.
[42]  oai:arXiv.org:1409.8302  [pdf] - 982073
CALIFA, the Calar Alto Legacy Integral Field Area survey. III. Second public data release
Comments: 33 pages, 20 figures, 7 tables, submitted to A&A. The CALIFA DR2 homepage http://califa.caha.es/DR2 is open for the public
Submitted: 2014-09-29
This paper describes the Second Public Data Release (DR2) of the Calar Alto Legacy Integral Field Area (CALIFA) survey. The data for 200 objects are made public, including the 100 galaxies of the First Public Data Release (DR1). Data were obtained with the integral-field spectrograph PMAS/PPak mounted on the 3.5 m telescope at the Calar Alto observatory. Two different spectral setups are available for each galaxy, (i) a low-resolution V500 setup covering the wavelength range 3745-7500 \AA\ with a spectral resolution of 6.0 \AA\ (FWHM), and (ii) a medium-resolution V1200 setup covering the wavelength range 3650-4840 \AA\ with a spectral resolution of 2.3 \AA\ (FWHM). The sample covers a redshift range between 0.005 and 0.03, with a wide range of properties in the Color-Magnitude diagram, stellar mass, ionization conditions, and morphological types. All released cubes were reduced with the latest pipeline, including improved spectrophotometric calibration, spatial registration and spatial resolution. The spectrophotometric calibration is better than 6% and the median spatial resolution is 2.5". Altogether the second data release contains over 1.5 million spectra. It is available at http://califa.caha.es/DR2.
[43]  oai:arXiv.org:1409.7123  [pdf] - 1222236
[CII] 158 $\mu$m Emission as a Star Formation Tracer
Comments: 25 pages, 14 figures, accepted for publication in ApJ
Submitted: 2014-09-24
The [CII] 157.74 $\mu$m transition is the dominant coolant of the neutral interstellar gas, and has great potential as a star formation rate (SFR) tracer. Using the Herschel KINGFISH sample of 46 nearby galaxies, we investigate the relation of [CII] surface brightness and luminosity with SFR. We conclude that [CII] can be used for measurements of SFR on both global and kiloparsec scales in normal star-forming galaxies in the absence of strong active galactic nuclei (AGN). The uncertainty of the $\Sigma_{\rm [CII]}-\Sigma_{\rm SFR}$ calibration is $\pm$0.21 dex. The main source of scatter in the correlation is associated with regions that exhibit warm IR colors, and we provide an adjustment based on IR color that reduces the scatter. We show that the color-adjusted $\Sigma_{\rm[CII]}-\Sigma_{\rm SFR}$ correlation is valid over almost 5 orders of magnitude in $\Sigma_{\rm SFR}$, holding for both normal star-forming galaxies and non-AGN luminous infrared galaxies. Using [CII] luminosity instead of surface brightness to estimate SFR suffers from worse systematics, frequently underpredicting SFR in luminous infrared galaxies even after IR color adjustment (although this depends on the SFR measure employed). We suspect that surface brightness relations are better behaved than the luminosity relations because the former are more closely related to the local far-UV field strength, most likely the main parameter controlling the efficiency of the conversion of far-UV radiation into gas heating. A simple model based on Starburst99 population-synthesis code to connect SFR to [CII] finds that heating efficiencies are $1\%-3\%$ in normal galaxies.
[44]  oai:arXiv.org:1409.1623  [pdf] - 1636570
Nearby supernova host galaxies from the CALIFA Survey: I. Sample, data analysis, and correlation to star-forming regions
Comments: (Accepted in A&A, 24 pages, 19 Figures, 9 Tables)
Submitted: 2014-09-04, last modified: 2014-09-17
[Abridged] We use optical IFS of nearby SN host galaxies provided by the CALIFA Survey with the goal of finding correlations in the environmental parameters at the location of different SN types. We recover the sequence in association of different SN types to the star-forming regions by using several indicators of the ongoing and recent SF related to both the ionized gas and the stellar populations. While the total ongoing SF is on average the same for the three SN types, SNe Ibc/IIb tend to happen closer to star-forming regions and occur in higher SF density locations compared to SNe II and SNe~Ia, the latter showing the weakest correlation. SNe~Ia host galaxies have on average masses that are $\sim$0.3-0.8~dex higher than CC SNe hosts due to a larger fraction of old stellar populations in the SNe~Ia hosts. Using the recent SN~Ia delay-time distribution and the SFHs of the galaxies, we show that the SN~Ia hosts in our sample should presently produce a factor 2 more SNe~Ia than the CC~SN hosts. Since both types are in hosts with similar SF rate and hence similar CC~SN rate, this can explain the mass difference between the SN~Ia and CC~SN hosts, and reinforce the finding that at least part of SNe~Ia should originate from very old progenitors. Comparing the mean SFH of the eight least massive galaxies to that of the massive SF SN~Ia hosts we find that the low-mass galaxies formed their stars over more extended time than the massive SN~Ia hosts. We estimate that the low-mass galaxies should produce by a factor of 10 less SNe~Ia, and a factor of 3 less CC~SNe than the high-mass group. Therefore the ratio between the number of CC~SNe and SNe~Ia is expected to increase with decreasing the galaxy mass. CC~SNe tend to explode at positions with younger stellar populations than the galaxy average, but the galaxy properties at SNe~Ia locations are one average the same as the global ones.
[45]  oai:arXiv.org:1407.6005  [pdf] - 1215845
Being WISE I: Validating Stellar Population Models and M/L ratios at 3.4 and 4.6 microns
Comments: 11 pages, 6 figures, submitted to ApJ
Submitted: 2014-07-22
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.
[46]  oai:arXiv.org:1406.3351  [pdf] - 1210129
The Nature of Extreme Emission Line Galaxies at z=1-2: Kinematics and Metallicities from Near-Infrared Spectroscopy
Comments: Resubmitted to ApJ after first referee report, 18 pages, 14 figures, 4 tables
Submitted: 2014-06-12
We present near-infrared spectroscopy of a sample of 22 Extreme Emission Line Galaxies at redshifts 1.3 < z < 2.3, confirming that these are low-mass (M* = 10^8 - 10^9 M_sun) galaxies undergoing intense starburst episodes (M*/SFR ~ 10-100 Myr). The sample is selected by [O III] or H{\alpha} emission line flux and equivalent width using near-infrared grism spectroscopy from the 3D-HST survey. High-resolution NIR spectroscopy is obtained with LBT/LUCI and VLT/X-SHOOTER. The [O III]/H{\beta} line ratio is high (> 5) and [N II]/H{\alpha} is always significantly below unity, which suggests a low gas-phase metallicity. We are able to determine gas-phase metallicities for 7 of our objects using various strong-line methods, with values in the range 0.05-0.30 Z_sun and with a median of 0.15 Z_sun; for 3 of these objects we detect [O III]{\lambda}4363 which allows for a direct constraint on the metallicity. The velocity dispersion, as measured from the nebular emission lines, is typically ~50 km/s. Combined with the observed star-forming activity, the Jeans and Toomre stability criteria imply that the gas fraction must be large (> 2/3), consistent with the difference between our dynamical and stellar mass estimates. The implied gas depletion time scale (several hundred Myr) is substantially longer than the inferred mass-weighted ages (~50 Myr), which further supports the emerging picture that most stars in low-mass galaxies form in short, intense bursts of star formation.
[47]  oai:arXiv.org:1403.2381  [pdf] - 1208288
A far-IR view of the starburst driven superwind in NGC 2146
Comments: 16 pages, 13 figures, accepted for publication in ApJ
Submitted: 2014-03-10, last modified: 2014-05-05
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.
[48]  oai:arXiv.org:1402.5210  [pdf] - 1203485
Reconstructing the stellar mass distributions of galaxies using S4G IRAC 3.6 and 4.5 micron images: II. The conversion from light to mass
Comments: 14 pages, 6 figures, accepted for publication in ApJ
Submitted: 2014-02-21
We present a new approach for estimating the 3.6 micron stellar mass-to-light ratio in terms of the [3.6]-[4.5] colors of old stellar populations. Our approach avoids several of the largest sources of uncertainty in existing techniques. By focusing on mid-IR wavelengths, we gain a virtually dust extinction-free tracer of the old stars, avoiding the need to adopt a dust model to correctly interpret optical or optical/NIR colors normally leveraged to assign M/L. By calibrating a new relation between NIR and mid-IR colors of GLIMPSE giant stars we also avoid discrepancies in model predictions for the [3.6]-[4.5] colors of old stellar populations due to uncertainties in molecular line opacities. We find that the [3.6]-[4.5] color, which is driven primarily by metallicity, provides a tight constraint on M/L_3.6, which varies intrinsically less than at optical wavelengths. The uncertainty on M/L_3.6 of ~0.07 dex due to unconstrained age variations marks a significant improvement on existing techniques for estimating the stellar M/L with shorter wavelength data. A single M/L_3.6=0.6 (assuming a Chabrier IMF), independent of [3.6]-[4.5] color, is also feasible as it can be applied simultaneously to old, metal-rich and young, metal-poor populations, and still with comparable (or better) accuracy (~0.1 dex) as alternatives. We expect our M/L_3.6 to be optimal for mapping the stellar mass distributions in S4G galaxies, for which we have developed an Independent Component Analysis technique to first isolate the old stellar light at 3.6 micron from non-stellar emission (e.g. hot dust and the 3.3 PAH feature). Our estimate can also be used to determine the fractional contribution of non-stellar emission to global (rest-frame) 3.6 micron fluxes, e.g. in WISE imaging, and establishes a reliable basis for exploring variations in the stellar IMF.
[49]  oai:arXiv.org:1401.3759  [pdf] - 791894
The PdBI Arcsecond Whirlpool Survey (PAWS): Multi-phase cold gas kinematic of M51
Comments: 42 pages, 12 figures, accepted for publication in ApJ
Submitted: 2014-01-15
The kinematic complexity and the favorable position of M51 on the sky make this galaxy an ideal target to test different theories of spiral arm dynamics. Taking advantage of the new high resolution PdBI Arcsecond Whirlpool Survey (PAWS) data, we undertake a detailed kinematic study of M51 to characterize and quantify the origin and nature of the non-circular motions. Using a tilted-ring analysis supported by several other archival datasets we update the estimation of M51's position angle (PA=(173 +/- 3) deg) and inclination (i=(22 +/- 5) deg). Harmonic decomposition of the high resolution (40 pc) CO velocity field shows the first kinematic evidence of an m=3 wave in the inner disk of M51 with a corotation at R(CR,m=3)=1.1 +/- 0.1 kpc and a pattern speed of Omega_p(m=3) = 140 km/(s kpc). This mode seems to be excited by the nuclear bar, while the beat frequencies generated by the coupling between the m=3 mode and the main spiral structure confirm its density-wave nature. We observe also a signature of an m=1 mode that is likely responsible for the lopsidedness of M51 at small and large radii. We provide a simple method to estimate the radial variation of the amplitude of the spiral perturbation (Vsp) attributed to the different modes. The main spiral arm structure has <Vsp>=50-70 km/s, while the streaming velocity associated with the m=1 and m=3 modes is, in general, 2 times lower. Our joint analysis of HI and CO velocity fields at low and high spatial resolution reveals that the atomic and molecular gas phases respond differently to the spiral perturbation due to their different vertical distribution and emission morphology.
[50]  oai:arXiv.org:1401.1505  [pdf] - 791862
The PdBI Arcsecond Whirlpool Survey (PAWS): Environmental Dependence of Giant Molecular Cloud Properties in M51
Comments: 73 pages, 18 figures, 14 tables, accepted for publication in ApJ
Submitted: 2014-01-07
Using data from the PdBI Arcsecond Whirlpool Survey (PAWS), we have generated the largest extragalactic Giant Molecular Cloud (GMC) catalog to date, containing 1,507 individual objects. GMCs in the inner M51 disk account for only 54% of the total 12CO(1-0) luminosity of the survey, but on average they exhibit physical properties similar to Galactic GMCs. We do not find a strong correlation between the GMC size and velocity dispersion, and a simple virial analysis suggests that 30% of GMCs in M51 are unbound. We have analyzed the GMC properties within seven dynamically-motivated galactic environments, finding that GMCs in the spiral arms and in the central region are brighter and have higher velocity dispersions than inter-arm clouds. Globally, the GMC mass distribution does not follow a simple power law shape. Instead, we find that the shape of the mass distribution varies with galactic environment: the distribution is steeper in inter-arm region than in the spiral arms, and exhibits a sharp truncation at high masses for the nuclear bar region. We propose that the observed environmental variations in the GMC properties and mass distributions are a consequence of the combined action of large-scale dynamical processes and feedback from high mass star formation. We describe some challenges of using existing GMC identification techniques for decomposing the 12CO(1-0) emission in molecule-rich environments, such as M51's inner disk.
[51]  oai:arXiv.org:1309.3453  [pdf] - 1179206
A Comparative Study of Giant Molecular Clouds in M51, M33 and the Large Magellanic Cloud
Comments: 25 pages, 11 figures, accepted for publication in ApJ. Uses emulateapj LaTeX macros. For more information on PAWS, further papers and data, see http://www.mpia.de/PAWS/
Submitted: 2013-09-13
We compare the properties of giant molecular clouds (GMCs) in M51 identified by the Plateau de Bure Interferometer Whirlpool Arcsecond Survey (PAWS) with GMCs identified in wide-field, high resolution surveys of CO emission in M33 and the Large Magellanic Cloud (LMC). We find that GMCs in M51 are larger, brighter and have higher velocity dispersions relative to their size than equivalent structures in M33 and the LMC. These differences imply that there are genuine variations in the average mass surface density of the different GMC populations. To explain this, we propose that the pressure in the interstellar medium surrounding the GMCs plays a role in regulating their density and velocity dispersion. We find no evidence for a correlation between size and linewidth in any of M51, M33 or the LMC when the CO emission is decomposed into GMCs, although moderately robust correlations are apparent when regions of contiguous CO emission (with no size limitation) are used. Our work demonstrates that observational bias remains an important obstacle to the identification and study of extragalactic GMC populations using CO emission, especially in molecule-rich galactic environments.
[52]  oai:arXiv.org:1309.1444  [pdf] - 1179038
Morphological Parameters of Spitzer Survey of Stellar Structure in Galaxies (S4G)
Comments: 22 pages, 19 figures, 7 tables. Full catalogs available on request. Accepted by ApJ
Submitted: 2013-09-05
The morphology of galaxies can be quantified to some degree using a set of scale-invariant parameters. Concentration (C), Asymmetry (A), Smoothness (S), the Gini index (G), relative contribution of the brightest pixels to the second order moment of the flux (M20), ellipticity (E), and the Gini index of the second order moment (GM) have all been applied to morphologically classify galaxies at various wavelengths. Here we present a catalog of these parameters for the Spitzer Survey of Stellar Structure in Galaxies (S4G), a volume-limited near-infrared imaging survey of nearby galaxies using the 3.6 and 4.5 micron channels of the IRAC camera. Our goal is to provide a reference catalog of near-infrared quantified morphology for high-redshift studies and galaxy evolution models with enough detail to resolve stellar mass morphology. We explore where normal, non-interacting galaxies -those typically found on the Hubble tuning fork- lie in this parameter space and show that there is a tight relation between Concentration and M20 for normal galaxies. M20 can be used to classify galaxies into earlier and later types (e.g., to separate spirals from irregulars). Several criteria using these parameters exist to select systems with a disturbed morphology, i.e., those that appear to be undergoing a tidal interaction. We examine the applicability of these criteria to Spitzer near-infrared imaging. We find that four relations, based on the parameters A&S, G&M20, GM, and C&M20, respectively, select outliers in the morphological parameter space, but each selects different subsets of galaxies. Two criteria (GM > 0.6, G > -0.115 x M20 + 0.384) seem most appropriate to identify possible mergers and the merger fraction in near-infrared surveys. We find no strong relation between lopsidedness and most of these morphological parameters, except for a weak dependence of lopsidedness on Concentration and M20.
[53]  oai:arXiv.org:1212.1208  [pdf] - 1158309
The CO-to-H2 Conversion Factor and Dust-to-Gas Ratio on Kiloparsec Scales in Nearby Galaxies
Comments: accepted to ApJ July 9, 2013, 63 pages, 2-column ApJ format, p38-63 are an image atlas
Submitted: 2012-12-05, last modified: 2013-08-18
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.
[54]  oai:arXiv.org:1307.7709  [pdf] - 1173042
X-ray nuclear activity in S4G barred galaxies: No link between bar strength and co-occurrent supermassive black hole fueling
Comments: 16 pages, 6 figures, accepted for publication in ApJ
Submitted: 2013-07-29
Stellar bars can lead to gas inflow toward the center of a galaxy and stimulate nuclear star formation. However, there is no compelling evidence on whether they also feed a central supermassive black hole: by measuring the fractions of barred active and inactive galaxies, previous studies have yielded conflicting results. In this paper, we aim to understand the lack of observational evidence for bar-driven active galactic nucleus (AGN) activity by studying a sample of 41 nearby (d < 35 Mpc) barred galaxies from the Spitzer Survey for Stellar Structure in Galaxies. We use Chandra observations to measure nuclear 2--10 keV X-ray luminosities and estimate Eddington ratios, together with Spitzer 3.6um imaging to quantify the strength of the stellar bar in two independent ways: (1) from its structure, as traced by its ellipticity and boxiness, and (2) from its gravitational torque Q_b, taken as the maximum ratio of the tangential force to the mean background radial force. In this way, rather than discretizing the presence of both stellar bars and nuclear activity, we are able to account for the continuum of bar strengths and degrees of AGN activity. We find nuclear X-ray sources in 31 out of 41 galaxies with median X-ray luminosity and Eddington ratio of L_X=4.3x10^{38} erg/s and L_bol/L_Edd=6.9x10^{-6} respectively, consistent with low-luminosity AGN activity. Including upper limits for those galaxies without nuclear detections, we find no significant correlation between any of the bar strength indicators and the degree of nuclear activity, irrespective of galaxy luminosity, stellar mass, Hubble type, or bulge size. Strong bars do not favor brighter or more efficient nuclear activity, implying that at least for the low-luminosity regime, supermassive black hole fueling is not closely connected to large scale features.
[55]  oai:arXiv.org:1304.1801  [pdf] - 1165778
The PdBI Arcsecond Whirlpool Survey (PAWS). I. A Cloud-Scale/Multi-Wavelength View of the Interstellar Medium in a Grand-Design Spiral Galaxy
Comments: 52 pages, 14 figures, 2 tables, (several minor typos corrected) accepted by ApJ, high resolution version available, see http://www.mpia.de/PAWS/pub/paws_schinnerer.pdf ; for more information on PAWS, further papers and the data, see http://www.mpia.de/PAWS
Submitted: 2013-04-05, last modified: 2013-07-22
The PdBI (Plateau de Bure Interferometer) Arcsecond Whirlpool Survey (PAWS) has mapped the molecular gas in the central ~9kpc of M51 in its 12CO(1-0) line emission at cloud-scale resolution of ~40pc using both IRAM telescopes. We utilize this dataset to quantitatively characterize the relation of molecular gas (or CO emission) to other tracers of the interstellar medium (ISM), star formation and stellar populations of varying ages. Using 2-dimensional maps, a polar cross-correlation technique and pixel-by-pixel diagrams, we find: (a) that (as expected) the distribution of the molecular gas can be linked to different components of the gravitational potential, (b) evidence for a physical link between CO line emission and radio continuum that seems not to be caused by massive stars, but rather depend on the gas density, (c) a close spatial relation between the PAH and molecular gas emission, but no predictive power of PAH emission for the molecular gas mass,(d) that the I-H color map is an excellent predictor of the distribution (and to a lesser degree the brightness) of CO emission, and (e) that the impact of massive (UV-intense) young star-forming regions on the bulk of the molecular gas in central ~9kpc can not be significant due to a complex spatial relation between molecular gas and star-forming regions that ranges from co-spatial to spatially offset to absent. The last point, in particular, highlights the importance of galactic environment -- and thus the underlying gravitational potential -- for the distribution of molecular gas and star formation.
[56]  oai:arXiv.org:1306.5758  [pdf] - 1172264
Explaining two circumnuclear star forming rings in NGC5248
Comments: Accepted for publication in A&A, 14pages + long table
Submitted: 2013-06-24
The distribution of gas in the central kiloparsec of a galaxy has a dynamically rapid evolution. Nonaxisymmetries in the gravitational potential of the galactic disk, such as a large scale stellar bar or spiral, can lead to significant radial motion of gaseous material from larger radii to the central region. The large influx of gas and the subsequent star formation keep the central region constantly changing. However, the ability of gas to reach the nucleus proper to fuel an AGN phase is not guaranteed. Gas inflow can be halted at a circumnuclear star forming ring several hundred parsec away. The nearby galaxy NGC5248 is especially interesting in this sense since it is said to host 2 circumnuclear star forming rings at 100pc and 370pc from its quiescent nucleus. Here we present new subarcsecond PdBI+30m CO(2-1) emission line observations of the central region. For the first time the molecular gas distribution at the smallest stellar ring is resolved into a gas ring, consistent with the presence of a quiescent nucleus. However, the molecular gas shows no ring structure at the larger ring. We combine analyses of the gaseous and stellar content in the central kiloparsec of this galaxy to understand the gas distribution and dynamics of this star forming central region. We discuss the probability of two scenarios leading to the current observations, given our full understanding of this system, and discuss whether there are really two circumnuclear star forming rings in this galaxy.
[57]  oai:arXiv.org:1305.2940  [pdf] - 1166600
On the Origin of Lopsidedness in Galaxies as Determined from the Spitzer Survey of Stellar Structure in Galaxies (S4G)
Comments: Accepted for Publication in ApJ (revised version only includes corrections to affiliations)
Submitted: 2013-05-13, last modified: 2013-05-28
We study the m=1 distortions (lopsidedness) in the stellar components of 167 nearby galaxies that span a wide range of morphologies and luminosities. We confirm the previous findings of 1) a high incidence of lopsidedness in the stellar distributions, 2) increasing lopsidedness as a function of radius out to at least 3.5 exponential scale lengths, and 3) greater lopsidedness, over these radii, for galaxies of later type and lower surface brightness. Additionally, the magnitude of the lopsidedness 1) correlates with the character of the spiral arms (stronger arm patterns occur in galaxies with less lopsidedness), 2) is not correlated with the presence or absence of a bar, or the strength of the bar when one is present, 3) is inversely correlated to the stellar mass fraction, f_*, within one radial scale length, and 4) correlates directly with f_* measured within the radial range over which we measure lopsidedness. We interpret these findings to mean that lopsidedness is a generic feature of galaxies and does not, generally, depend on a rare event, such as a direct accretion of a satellite galaxy onto the disk of the parent galaxy. While lopsidedness may be caused by several phenomena, moderate lopsidedness (<A_1>_i + <A_1>_o)/2 < 0.3) is likely to reflect halo asymmetries to which the disk responds or a gravitationally self-generated mode . We hypothesize that the magnitude of the stellar response depends both on how centrally concentrated the stars are with respect to the dark matter and whether there are enough stars in the region of the lopsidedness that self-gravity is dynamically important.
[58]  oai:arXiv.org:1305.2923  [pdf] - 1166594
Mapping dust through emission and absorption in nearby galaxies
Comments: 17 pages, 10 figures, accepted for publication in ApJ
Submitted: 2013-05-13
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.
[59]  oai:arXiv.org:1304.7910  [pdf] - 1166314
Gas Kinematics on GMC scales in M51 with PAWS: cloud stabilization through dynamical pressure
Comments: 28 pages, 14 figures, accepted for publication in ApJ
Submitted: 2013-04-30
We use the high spatial and spectral resolution of the PAWS CO(1-0) survey of the inner 9 kpc of the iconic spiral galaxy M51 to examine the effect of gas streaming motions on the star-forming properties of individual GMCs. We compare our view of gas flows in M51 -- which arise due to departures from axi-symmetry in the gravitational potential (i.e. the nuclear bar and spiral arms) -- with the global pattern of star formation as traced by Halpha and 24\mu m emission. We find that the dynamical environment of GMCs strongly affects their ability to form stars, in the sense that GMCs situated in regions with large streaming motions can be stabilized, while similarly massive GMCs in regions without streaming go on to efficiently form stars. We argue that this is the result of reduced surface pressure felt by clouds embedded in an ambient medium undergoing large streaming motions, which prevents collapse. Indeed, the variation in gas depletion time expected based on the observed streaming motions throughout the disk of M51 quantitatively agrees with the variation in observed gas depletion time scale. The example of M51 shows that streaming motions, triggered by gravitational instabilities in the form of bars and spiral arms, can alter the star formation law; this can explain the variation in gas depletion time among galaxies with different masses and morphologies. In particular, we can explain the long gas depletion times in spiral galaxies compared to dwarf galaxies and starbursts. We suggest that adding a dynamical pressure term to the canonical free-fall time produces a single star formation law that can be applied to all star-forming regions and galaxies, across cosmic time.
[60]  oai:arXiv.org:1304.6083  [pdf] - 1166150
The impact of bars on disk breaks as probed by S4G imaging
Comments: Accepted for publication in ApJ. 32 pages, 14 figures. The whole figureset of Fig. 14 can be downloaded from http://www.cv.nrao.edu/~jmunoz/f14_large.pdf
Submitted: 2013-04-22
We have analyzed the radial distribution of old stars in a sample of 218 nearby face-on disks, using deep 3.6um images from the Spitzer Survey of Stellar Structure in Galaxies (S4G). In particular, we have studied the structural properties of those disks with a broken or down-bending profile. We find that, on average, disks with a genuine single exponential profile have a scale-length and a central surface brightness which are intermediate to those of the inner and outer components of a down-bending disk with the same total stellar mass. In the case of barred galaxies, the ratio between the break and the bar radii (Rbr/Rbar) depends strongly on the total stellar mass of the galaxy. For galaxies more massive than 10^10 Msun, the distribution is bimodal, peaking at Rbr/Rbar~2 and ~3.5. The first peak, which is the most populated one, is linked to the Outer Lindblad Resonance of the bar, whereas the second one is consistent with a dynamical coupling between the bar and the spiral pattern. For galaxies below 10^10 Msun, breaks are found up to ~10 Rbar, but we show that they could still be caused by resonances given the rising nature of rotation curves in these low-mass disks. While not ruling out star formation thresholds, our results imply that radial stellar migration induced by non-axysymmetric features can be responsible not only for those breaks at 2 Rbar, but also for many of those found at larger radii.
[61]  oai:arXiv.org:1304.1396  [pdf] - 1165741
The Plateau de Bure + 30m Arcsecond Whirlpool Survey reveals a thick disk of diffuse molecular gas in the M51 galaxy
Comments: 42 pages, 35 PostScript figures. Accepted for publication in ApJ. Uses emulateapj LaTeX macros
Submitted: 2013-04-04
We present the data of the Plateau de Bure Arcsecond Whirlpool Survey (PAWS), a high spatial and spectral resolution 12CO(1-0) line survey of the inner 10 x 6 kpc of the M51 system, and the first wide-field imaging of molecular gas in a star-forming spiral galaxy with resolution matched to the typical size of Giant Molecular Clouds (40 pc). We describe the observation, reduction, and combination of the Plateau de Bure Interferometer (PdBI) and IRAM-30m "short spacing" data. The final data cube attains 1.1"-resolution over the 270" x 170" field of view, with sensitivity to all spatial scales from the combination of PdBI and IRAM-30m data, and brightness sensitivity of 0.4 K (1 sigma) in each 5 km/s-wide channel map. We find a CO-luminosity of 9.10^8 K km/s pc^2, corresponding to a molecular gas mass of 4.10^9 Msol for a standard CO-to-H2 conversion factor. Unexpectedly, we find that a large fraction, (50 +/- 10)%, of this emission arises mostly from spatial scales larger than 36" or about 1.3 kpc. Through a series of tests, we demonstrate that this extended emission does not result from a processing artifact. We discuss its origin in light of the stellar component, the 12CO/13CO ratio, and the difference between the kinematics and structure of the PdBI-only and hybrid synthesis (PdBI + IRAM-30m) images. The extended emission is consistent with a thick, diffuse disk of molecular gas with a typical scale height of about 200 pc, substructured in unresolved filaments which fills about 0.1% of the volume.
[62]  oai:arXiv.org:1304.1219  [pdf] - 1165723
Probability Distribution Functions OF 12CO(J = 1-0) Brightness and Integrated Intensity in M51: The PAWS View
Comments: 30 pages, 18 figures, accepted to ApJ PAWS special issue
Submitted: 2013-04-03
We analyse the distribution of CO brightness temperature and integrated intensity in M51 at ~40 pc resolution using new CO data from the Plateau de Bure Arcsecond Whirlpool Survey (PAWS). We present probability distribution functions (PDFs) of the CO emission within the PAWS field, which covers the inner 11 x 7 kpc of M51. We find variations in the shape of CO PDFs within different M51 environments, and between M51 and M33 and the Large Magellanic Cloud (LMC). Globally, the PDFs for the inner disk of M51 can be represented by narrow lognormal functions that cover 1 to 2 orders of magnitude in CO brightness and integrated intensity. The PDFs for M33 and the LMC are narrower and peak at lower CO intensities. However, the CO PDFs for different dynamical environments within the PAWS field depart from the shape of the global distribution. The PDFs for the interarm region are approximately lognormal, but in the spiral arms and central region of M51, they exhibit diverse shapes with a significant excess of bright CO emission. The observed environmental dependence of the shape of the CO PDFs is qualitatively consistent with changes that would be expected if molecular gas in the spiral arms has a larger range of average densities, gas temperatures and velocity fluctuations, though further work is required to disentangle the importance of large-scale dynamical effects versus star formation feedback in regulating these properties. We show that the shape of the CO PDFs for different M51 environments is only weakly related to global properties of the CO emission, but is strongly correlated with some properties of the local giant molecular cloud (GMC) and young stellar cluster populations. For galaxies with strong spiral structure such as M51, our results indicate that galactic-scale dynamical processes play a significant role in the formation and evolution of GMCs and stellar clusters.(abridged)
[63]  oai:arXiv.org:1301.6884  [pdf] - 1159343
A Detailed Study of the Radio--FIR Correlation in NGC6946 with Herschel-PACS/SPIRE from KINGFISH
Comments: 23 pages, 13 figures, accepted for publication in Astronomy & Astrophysics Journal
Submitted: 2013-01-29
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.
[64]  oai:arXiv.org:1210.8150  [pdf] - 619031
CALIFA, the Calar Alto Legacy Integral Field Area survey. II. First public data release
Comments: 27 pages, 25 figures, 6 tables, final version including corrections after proof reading, the CALIFA DR1 homepage at http://califa.caha.es/DR1 is open for the public
Submitted: 2012-10-30, last modified: 2013-01-25
We present the first public data release of the CALIFA survey. It consists of science-grade optical datacubes for the first 100 of eventually 600 nearby (0.005<z<0.03) galaxies, obtained with the integral-field spectrograph PMAS/PPak mounted on the 3.5m telescope at the Calar Alto observatory. The galaxies in DR1 already cover a wide range of properties in color-magnitude space, morphological type, stellar mass, and gas ionization conditions. This offers the potential to tackle a variety of open questions in galaxy evolution using spatially resolved spectroscopy. Two different spectral setups are available for each galaxy, (i) a low-resolution V500 setup covering the nominal wavelength range 3745-7500A with a spectral resolution of 6.0A (FWHM), and (ii) a medium-resolution V1200 setup covering the nominal wavelength range 3650-4840A with a spectral resolution of 2.3A (FWHM). We present the characteristics and data structure of the CALIFA datasets that should be taken into account for scientific exploitation of the data, in particular the effects of vignetting, bad pixels and spatially correlated noise. The data quality test for all 100 galaxies showed that we reach a median limiting continuum sensitivity of 1.0x10^-18erg/s/cm^2/A/arcsec^2 at 5635A and 2.2x10^-18erg/s/cm^2/A/arcsec^2 at 4500A for the V500 and V1200 setup respectively, which corresponds to limiting r and g band surface brightnesses of 23.6mag/arcsec^2 and 23.4mag/arcsec^2, or an unresolved emission-line flux detection limit of roughly 1x10^-17erg/s/cm^2/arcsec^2 and 0.6x10^-17erg/s/cm^2/arcsec^2, respectively. The median spatial resolution is 3.7", and the absolute spectrophotometric calibration is better than 15% (1sigma). We also describe the available interfaces and tools that allow easy access to this first public CALIFA data.
[65]  oai:arXiv.org:1301.2328  [pdf] - 1158952
Molecular Gas and Star Formation in Nearby Disk Galaxies
Comments: 39 pages (incl. 10 page appendix), 19 figures, figures degraded for astro-ph. Until publication a higher resolution version available at http://www.cv.nrao.edu/~aleroy/heracles_sfco.pdf - we suggest that one
Submitted: 2013-01-10
We compare molecular gas traced by 12CO(2-1) maps from the HERACLES survey, with tracers of the recent star formation rate (SFR) across 30 nearby disk galaxies. We demonstrate a first-order linear correspondence between Sig_mol and Sig_SFR but also find important second-order systematic variations in the apparent molecular gas depletion time, t_dep^mol = Sig_mol / Sig_SFR. At our 1 kpc common resolution, CO correlates closely with many tracers of the recent SFR. Weighting each line of sight equally and using a fixed, Milky Way alpha_CO, our data yield a molecular gas depletion time, t_dep^mol=Sig_mol/Sig_SFR ~ 2.2 Gyr with 0.3 dex scatter, in good agreement with literature data. We apply a forward-modeling approach to constrain the power-law index, N, that relates the SFR surface density and the molecular gas surface density and find N=1+/-0.15 for our full data set with some variation from galaxy to galaxy. However, we caution that a power law treatment oversimplifies the topic given that we observe correlations between t_dep^mol and other local and global quantities. The strongest of these are a decreased t_dep^mol in low-mass, low-metallicity galaxies and a correlation of the kpc-scale t_dep^mol with dust-to-gas ratio, D/G. These correlations can be explained by a CO-to-H2 conversion factor that depends on D/G in the theoretically expected way. This is not a unique interpretation, but external evidence of conversion factor variations makes it a conservative one. After applying a D/G-dependent alpha_CO, some weak correlations between t_dep^mol and local conditions persist. In particular, we observe lower t_dep^mol and enhanced CO excitation associated with some nuclear gas concentrations. These appear to reflect real enhancements in the SFR/H2 and t_dep appears multivalued at fixed Sig_mol, supporting the the idea of "disk" and "starburst" modes driven by environmental factors.
[66]  oai:arXiv.org:1207.2777  [pdf] - 1124768
Cool Dust in the Outer Ring of NGC 1291
Comments: Accepted to ApJ, 9 pages, 6 figures
Submitted: 2012-07-11
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.
[67]  oai:arXiv.org:1204.0554  [pdf] - 1117726
A Study of Heating and Cooling of the ISM in NGC 1097 with Herschel-PACS and Spitzer-IRS
Comments: 17 pages, 14 figures, 5 tables; accepted for publication in ApJ
Submitted: 2012-04-02
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.
[68]  oai:arXiv.org:1204.0524  [pdf] - 1117721
Converting From 3.6 and 4.5 Micron Fluxes to Stellar Mass
Comments: 7 pages, 6 figures, to be published in the Astronomical Journal
Submitted: 2012-04-02
We use high spatial resolution maps of stellar mass and infrared flux of the Large Magellanic Cloud (LMC) to calibrate a conversion between 3.6 and 4.5 micron fluxes and stellar mass, M_* = 10^{5.65} * F_{3.6}^{2.85} * F_{4.5}^{-1.85} * (D/0.05)^2 M_solar, where fluxes are in Jy and D is the luminosity distance to the source in Mpc, and to provide an approximate empirical estimate of the fractional internal uncertainty in M_* of 0.3*sqrt{N/10^6}, where N is the number of stars in the region. We find evidence that young stars and hot dust contaminate the measurements, but attempts to remove this contamination using data that is far superior than what is generally available for unresolved galaxies resulted in marginal gains in accuracy. The scatter among mass estimates for regions in the LMC is comparable to that found by previous investigators when modeling composite populations, and so we conclude that our simple conversion is as precise as possible for the data and models currently available. Our results allow for a reasonably bottom-heavy initial mass function, such as Salpeter or heavier, and moderately disfavor lighter versions such as a diet-Salpeter or Chabrier initial mass function.
[69]  oai:arXiv.org:1203.0467  [pdf] - 1117034
The S4G perspective on circumstellar dust extinction of AGB stars in M100
Comments: 6 pages, 4 figures, 1 table, accepted for publication in ApJL
Submitted: 2012-03-02
We examine the effect of circumstellar dust extinction on the near-IR contribution of asymptotic giant branch (AGB) stars in intermediate-age clusters throughout the disk of M100. For our sample of 17 AGB-dominated clusters we extract optical-to-mid-IR SEDs and find that NIR brightness is coupled to the mid-IR dust emission in such a way that a significant reduction of AGB light, of up to 1 mag in K-band, follows from extinction by the dust shell formed during this stage. Since the dust optical depth varies with AGB chemistry (C-rich or O-rich), our results suggest that the contribution of AGB stars to the flux from their host clusters will be closely linked to the metallicity and the progenitor mass of the AGB star, to which dust chemistry and mass-loss rate are sensitive. Our sample of clusters--each the analogue of a ~1 Gyr old post-starburst galaxy--has implications within the context of mass and age estimation via SED modelling at high z: we find that the average ~0.5 mag extinction estimated here may be sufficient to reduce the AGB contribution in (rest-frame) K-band from ~70%, as predicted in the latest generation of synthesis models, to ~35%. Our technique for selecting AGB-dominated clusters in nearby galaxies promises to be effective for discriminating the uncertainties associated with AGB stars in intermediate-age populations that plague age and mass estimation in high-z galaxies.
[70]  oai:arXiv.org:1201.1016  [pdf] - 1092738
Resolving the Far-IR Line Deficit: Photoelectric Heating and Far-IR Line Cooling in NGC 1097 and NGC 4559
Comments: 19 pages, 14 pages, accepted for publication in ApJ
Submitted: 2012-01-04
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.
[71]  oai:arXiv.org:1112.1093  [pdf] - 1092174
Herschel Far-Infrared and Sub-millimeter Photometry for the KINGFISH Sample of Nearby Galaxies
Comments: Accepted for publication in ApJ
Submitted: 2011-12-05
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.
[72]  oai:arXiv.org:1111.4438  [pdf] - 1091782
KINGFISH -- Key Insights on Nearby Galaxies: A Far-Infrared Survey with Herschel: Survey Description and Image Atlas
Comments: 52 pages, 10 figures, 2 tables. Accepted for publication in PASP. A pdf copy of the image atlas can be downloaded separately at the KINGFISH public website, http://www.ast.cam.ac.uk/sites/default/files/Kingfish_reduced.pdf
Submitted: 2011-11-18
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.
[73]  oai:arXiv.org:1111.0962  [pdf] - 1453893
CALIFA, the Calar Alto Legacy Integral Field Area survey: I. Survey presentation
Comments: 32 pages, 29 figures, Accepted for publishing in Astronomy and Astrophysics
Submitted: 2011-11-03, last modified: 2011-11-04
We present here the Calar Alto Legacy Integral Field Area (CALIFA) survey, which has been designed to provide a first step in this direction.We summarize the survey goals and design, including sample selection and observational strategy.We also showcase the data taken during the first observing runs (June/July 2010) and outline the reduction pipeline, quality control schemes and general characteristics of the reduced data. This survey is obtaining spatially resolved spectroscopic information of a diameter selected sample of $\sim600$ galaxies in the Local Universe (0.005< z <0.03). CALIFA has been designed to allow the building of two-dimensional maps of the following quantities: (a) stellar populations: ages and metallicities; (b) ionized gas: distribution, excitation mechanism and chemical abundances; and (c) kinematic properties: both from stellar and ionized gas components. CALIFA uses the PPAK Integral Field Unit (IFU), with a hexagonal field-of-view of $\sim1.3\sq\arcmin'$, with a 100% covering factor by adopting a three-pointing dithering scheme. The optical wavelength range is covered from 3700 to 7000 {\AA}, using two overlapping setups (V500 and V1200), with different resolutions: R\sim850 and R\sim1650, respectively. CALIFA is a legacy survey, intended for the community. The reduced data will be released, once the quality has been guaranteed. The analyzed data fulfill the expectations of the original observing proposal, on the basis of a set of quality checks and exploratory analysis. We conclude from this first look at the data that CALIFA will be an important resource for archaeological studies of galaxies in the Local Universe.
[74]  oai:arXiv.org:1110.2683  [pdf] - 1084820
Reconstructing the stellar mass distributions of galaxies using S4G IRAC 3.6 and 4.5 micron images: I. Correcting for contamination by PAH, hot dust, and intermediate age stars
Comments: 19 pages, 5 figures, accepted for publication in ApJ
Submitted: 2011-10-12
With the aim of constructing accurate 2D maps of the stellar mass distribution in nearby galaxies from S4G 3.6 and 4.5 micron images, we report on the separation of the light from old stars from the emission contributed by contaminants (e.g. hot dust and the 3.3 micron PAH feature). Results for a small sample of six disk galaxies (NGC 1566, NGC 2976, NGC 3031, NGC 3184, NGC 4321, and NGC 5194) with a range of morphological properties, dust contents and star formation histories are presented to demonstrate our approach. We use an Independent Component Analysis (ICA) technique designed to separate statistically independent source distributions, maximizing the distinction in the [3.6]-[4.5] colors of the sources. The technique also removes emission from intermediate-age evolved red objects with a low mass-to-light ratio, such as asymptotic giant branch (AGB) and red supergiant (RSG) stars, revealing maps of the underlying old distribution of light with [3.6]-[4.5] colors consistent with the colors of K and M giants. Contaminants are identified via comparison to the non-stellar emission imaged at 8 microns, which is dominated by the broad PAH feature. Using the measured 3.6/8 micron ratio to select the individual contaminants, we find that hot dust and PAH together contribute between ~5-15% to the integrated light at 3.6 microns, while light from regions dominated by intermediate-age stars accounts for only 1-5%. Locally, however, the contribution from either contaminant can reach much higher levels; dust contributes on average 22% to the emission in star-forming regions throughout the sample, while intermediate age-stars contribute upwards of 50% in localized knots. The removal of these contaminants with ICA leaves maps of the old stellar disk that retain a high degree of structural information and are ideally suited for tracing the stellar mass, as will be the focus in a companion paper.
[75]  oai:arXiv.org:1106.4022  [pdf] - 1077423
The Emission by Dust and Stars of Nearby Galaxies in the Herschel KINGFISH Survey
Comments: 19 pages, 12 figures; published in ApJ
Submitted: 2011-06-20, last modified: 2011-08-16
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.
[76]  oai:arXiv.org:1101.4023  [pdf] - 1051578
The Displaced Dusty ISM of NGC 3077: Tidal Stripping in the M 81 Triplet
Comments: ApJ Letters, in print
Submitted: 2011-01-20
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.
[77]  oai:arXiv.org:1101.1491  [pdf] - 1042956
Infrared Spectroscopy of Halos of Edge-on Galaxies
Comments: Accepted for publication in The Astrophysical Journal
Submitted: 2011-01-07
We present a study of ionized gas, PAHs, and molecular hydrogen emission in the halos of three edge-on galaxies, NGC 891, NGC 5775 and NGC 3044, based on 10-20 micron Spitzer Space Telescope spectra. The [Ne III]/[Ne II] ratio, an excellent measure of radiation hardness, rises with z in the halo of NGC 891. It is also higher in the halo of NGC 5775 than in the disk. NGC 3044 presents a more confusing situation. To explain the [Ne III]/[Ne II] as well as optical line ratio behavior in NGC 891, we carry out a simple exploration of parameter space with CLOUDY, which indicates a large increase in radiation temperature with height. Illustrative examples of physical models using a Monte Carlo radiative transfer code show that the rising neon ratio may be explained by adding a vertically extended, hot stellar source to a thin disk of massive stars. However, several other sources of hard spectra may be relevant. PAH features have scale heights of 430--530 pc in NGC 891 and 720--1080 pc in NGC 5775, suggesting they can be transported by disk-halo flows. Within NGC 891 and NGC 5775, scale heights are similar for all PAHs. For NGC 891, the scale heights exceed that of 8 micron emission, indicating a transition from more ionized to more neutral PAHs with height. Most PAH equivalent widths are higher in the halos. Molecular hydrogen 17.03 micron emission with scale heights of 550-580 pc in NGC 891 and 850 pc in NGC 5775 suggests a molecular component in a surprisingly thick layer.
[78]  oai:arXiv.org:1010.1592  [pdf] - 1041155
The Spitzer Survey of Stellar Structure in Galaxies (S^4G)
Comments: Accepted for Publication in PASP, 14 pages, 13 figures
Submitted: 2010-10-08, last modified: 2010-10-12
The Spitzer Survey of Stellar Structure in Galaxies S^4G is an Exploration Science Legacy Program approved for the Spitzer post-cryogenic mission. It is a volume-, magnitude-, and size-limited (d < 40 Mpc, |b| > 30 degrees, m_(Bcorr) < 15.5, D25>1') survey of 2,331 galaxies using IRAC at 3.6 and 4.5 microns. Each galaxy is observed for 240 s and mapped to > 1.5 x D25. The final mosaicked images have a typical 1 sigma rms noise level of 0.0072 and 0.0093 MJy / sr at 3.6 and 4.5 microns, respectively. Our azimuthally-averaged surface brightness profile typically traces isophotes at mu_3.6 (AB) (1 sigma) ~ 27 mag arcsec^-2, equivalent to a stellar mass surface density of ~ 1 Msun pc^-2. S^4G thus provides an unprecedented data set for the study of the distribution of mass and stellar structures in the local Universe. This paper introduces the survey, the data analysis pipeline and measurements for a first set of galaxies, observed in both the cryogenic and warm mission phase of Spitzer. For every galaxy we tabulate the galaxy diameter, position angle, axial ratio, inclination at mu_3.6 (AB) = 25.5 and 26.5 mag arcsec^-2 (equivalent to ~ mu_B (AB) =27.2 and 28.2 mag arcsec^-2, respectively). These measurements will form the initial S^4G catalog of galaxy properties. We also measure the total magnitude and the azimuthally-averaged radial profiles of ellipticity, position angle, surface brightness and color. Finally, we deconstruct each galaxy using GALFIT into its main constituent stellar components: the bulge/spheroid, disk, bar, and nuclear point source, where necessary. Together these data products will provide a comprehensive and definitive catalog of stellar structures, mass and properties of galaxies in the nearby Universe.
[79]  oai:arXiv.org:1005.1959  [pdf] - 1026817
Mapping far-IR emission from the central kiloparsec of NGC 1097
Comments: Accepted for publication in the A&A Herschel Special Edition
Submitted: 2010-05-11
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.
[80]  oai:arXiv.org:1005.1887  [pdf] - 1026799
Enhanced dust heating in the bulges of early-type spiral galaxies
Comments: Accepted for publication on the A&A Herschel Special Issue
Submitted: 2010-05-11
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.
[81]  oai:arXiv.org:1005.1611  [pdf] - 1026755
Far-Infrared Line Imaging of the Starburst Ring in NGC 1097 with the Herschel/PACS Spectrometer
Comments: Accepted for publication on the A&A Herschel Special Issue
Submitted: 2010-05-10
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}$.
[82]  oai:arXiv.org:0907.3443  [pdf] - 1002985
Uncovering the Origins of Spiral Structure by Measuring Radial Variation in Pattern Speeds
Comments: 15 pages in emulateapj format, 12 figures, accepted for publication in ApJ
Submitted: 2009-07-20
Current theories of spiral and bar structure predict a variety of pattern speed behaviors, calling for detailed, direct measurement of the radial variation of pattern speeds. Our recently developed Radial Tremaine-Weinberg (TWR) method allows this goal to be achieved for the first time. Here we present TWR spiral pattern speed estimates for M101, IC 342, NGC 3938 and NGC 3344 in order to investigate whether spiral structure is steady or winding, whether spirals are described by multiple pattern speeds, and the relation between bar and spiral speeds. Where possible, we interpret our pattern speeds estimates according to the resonance radii associated with each (established with the disk angular rotation), and compare these to previous determinations. By analyzing the high-quality HI and CO data cubes available for these galaxies, we show that it is possible to determine directly multiple pattern speeds within these systems, and hence identify the characteristic signatures of the processes that drive the spiral structure. Even this small sample of galaxies reveals a surprisingly complex taxonomy, with the first direct evidence for the presence of resonant coupling of multiple patterns found in some systems, and the measurement of a simple single pattern speed in others. Overall, this study demonstrates that we are now in a position to uncover more of the apparently complex physics that lies behind spiral structure.
[83]  oai:arXiv.org:0807.1902  [pdf] - 14440
Radial Dependence of the Pattern Speed of M51
Comments: 14 pages in emulateapj format, 12 figures, accepted for publication in ApJ
Submitted: 2008-07-11
The grand-design spiral galaxy M51 has long been a crucial target for theories of spiral structure. Studies of this iconic spiral can address the question of whether strong spiral structure is transient (e.g. interaction-driven) or long-lasting. As a clue to the origin of the structure in M51, we investigate evidence for radial variation in the spiral pattern speed using the radial Tremaine-Weinberg (TWR) method. We implement the method on CO observations tracing the ISM-dominant molecular component. Results from the method's numerical implementation--combined with regularization, which smooths intrinsically noisy solutions--indicate two distinct patterns speeds inside 4 kpc at our derived major axis PA=170 deg., both ending at corotation and both significantly higher than the conventionally adopted global value. Inspection of the rotation curve suggests that the pattern speed interior to 2 kpc lacks an ILR, consistent with the leading structure seen in HST near-IR observations. We also find tentative evidence for a lower pattern speed between 4 and 5.3 kpc measured by extending the regularized zone. As with the original TW method, uncertainty in major axis position angle (PA) is the largest source of error in the calculation; in this study, where \delta PA=+/-5 deg. a ~20% error is introduced to the parameters of the speeds at PA=170 deg. Accessory to this standard uncertainty, solutions with PA=175 deg. (also admitted by the data) exhibit only one pattern speed inside 4 kpc, and we consider this circumstance under the semblance of a radially varying PA.
[84]  oai:arXiv.org:0711.4104  [pdf] - 7394
Tests of the Radial Tremaine-Weinberg Method
Comments: 22 pages in emulateapj format, 20 figures, accepted for publication in ApJ
Submitted: 2007-11-26
At the intersection of galactic dynamics, evolution and global structure, issues such as the relation between bars and spirals and the persistence of spiral patterns can be addressed through the characterization of the angular speeds of the patterns and their possible radial variation. The Radial Tremaine-Weinberg (TWR) Method, a generalized version of the Tremaine-Weinberg method for observationally determining a single, constant pattern speed, allows the pattern speed to vary arbitrarily with radius. Here, we perform tests of the TWR method with regularization on several simulated galaxy data sets. The regularization is employed as a means of smoothing intrinsically noisy solutions, as well as for testing model solutions of different radial dependence (e.g. constant, linear or quadratic). We test these facilities in studies of individual simulations, and demonstrate successful measurement of both bar and spiral pattern speeds in a single disk, secondary bar pattern speeds, and spiral winding (in the first application of a TW calculation to a spiral simulation). We also explore the major sources of error in the calculation and find uncertainty in the major axis position angle most dominant. In all cases, the method is able to extract pattern speed solutions where discernible patterns exist to within 20% of the known values, suggesting that the TWR method should be a valuable tool in the area of galactic dynamics. For utility, we also discuss the caveats in, and compile a prescription for, applications to real galaxies.
[85]  oai:arXiv.org:astro-ph/0511069  [pdf] - 77427
The Lifetime of Grand Design
Comments: 6 pages, 3 figures, accepted for publication in MNRAS
Submitted: 2005-11-02
The lifetime of the structure in grand design spiral galaxies is observationally ill-determined, but is essentially set by how accurately the pattern's rotation can be characterized by a single angular pattern speed. This paper derives a generalized version of the Tremaine-Weinberg method for observationally determining pattern speeds, in which the pattern speed is allowed to vary arbitrarily with radius. The departures of the derived pattern speed from a constant then provides a simple metric of the lifetime of the spiral structure. Application of this method to CO observations of NGC 1068 reveal that the pattern speed of the spiral structure in this galaxy varies rapidly with radius, and that the lifetime of the spiral structure is correspondingly very short. If this result turns out to be common in grand-design spiral galaxies, then these features will have to be viewed as highly transient phenomena.
[86]  oai:arXiv.org:astro-ph/0509763  [pdf] - 1468891
Pattern Speeds of Spiral Galaxies Using the Tremaine-Weinberg Method
Comments: 4 pages, 3 figures. To appear in the proceedings "Island Universes: Structure and Evolution of Disk Galaxies", ed. R. de Jong (Springer: Dordrecht)
Submitted: 2005-09-26
Spiral galaxies displaying molecule-dominated ISMs are potential candidates for pattern speed determination by the Tremaine-Weinberg (TW) method (1984); the molecular gas of such galaxies, as traced by CO emission, is argued to approximately obey the continuity equation on orbital timescales--a key requirement of the method. We continue our work on this application of the Tremaine-Weinberg method with pattern speed determinations of several galaxies using CO emission data cubes from the BIMA SONG survey, wherein the molecular dominance of each of the three galaxies in our sample (NGC 5457, 5033, and 5055) is established using radial profiles of Wong and Blitz (2002).