Source author record

Eva Schinnerer

Eva Schinnerer appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

117works
5topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

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

Building this map preview

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

Published work

117 published item(s)

preprint2022arXiv

A 2-3 mm high-resolution molecular line survey towards the centre of the nearby spiral galaxy NGC 6946

The complex physical, kinematic, and chemical properties of galaxy centres make them interesting environments to examine with molecular line emission. We present new $2-4$" (${\sim}75{-}150$ pc at $7.7$ Mpc) observations at 2 and 3 mm covering the central $50$" (${\sim}1.9$ kpc) of the nearby double-barred spiral galaxy NGC 6946 obtained with the IRAM Plateau de Bure Interferometer. We detect spectral lines from ten molecules: CO, HCN, HCO$^+$, HNC, CS, HC$_3$N, N$_2$H$^+$, C$_2$H, CH$_3$OH, and H$_2$CO. We complemented these with published 1mm CO observations and 33 GHz continuum observations to explore the star formation rate surface density ${Σ_{\mathrm{SFR}}}$ on 150 pc scales. In this paper, we analyse regions associated with the inner bar of NGC 6946 $-$ the nuclear region (NUC), the northern (NBE), and southern inner bar end (SBE) and we focus on short-spacing corrected bulk (CO) and dense gas tracers (HCN, HCO$^+$, and HNC). We find that HCO$^+$ correlates best with ${Σ_{\mathrm{SFR}}}$, but the dense gas fraction ($f_{\mathrm{dense}}$) and star formation efficiency of the dense gas (${\mathrm{SFE_{dense}}}$) fits show different behaviours than expected from large-scale disc observations.The SBE has a higher ${Σ_{\mathrm{SFR}}}$, $f_{\mathrm{dense}}$, and shocked gas fraction than the NBE. We examine line ratio diagnostics and find a higher CO(2-1)/CO(1-0) ratio towards NBE than for the NUC. Moreover, comparison with existing extragalactic datasets suggests that using the HCN/HNC ratio to probe kinetic temperatures is not suitable on kiloparsec and sub-kiloparsec scales in extragalactic regions. Lastly, our study shows that the HCO$^+$/HCN ratio might not be a unique indicator to diagnose AGN activity in galaxies.

preprint2022arXiv

A CO isotopologue Line Atlas within the Whirlpool galaxy Survey (CLAWS)

We present the CO isotopologue Line Atlas within the Whirpool galaxy Survey (CLAWS) based on an IRAM 30-m large programme which provides a benchmark study of numerous, faint CO isotopologues in the mm-wavelength regime across the full disc of M51 (NGC 5194). The survey's core goal is to use the low-J CO isotopologue lines to constrain CO excitation and chemistry, and therefrom the local physical conditions of the gas. In this survey paper, we describe the CLAWS observing and data reduction strategies. We map the J=1-0 and 2-1 transitions of the CO isotopologues $^{12}$CO,$^{13}$CO, C$^{18}$O and C$^{17}$O, as well as several supplementary lines within the 1 mm and 3 mm window (CN(1-0), CS(2-1), CH$_3$OH(2-1), N$_2$H$^+$(1-0), HC$_3$N(10-9)) at ~1 kpc resolution. A total observation time of 149 h offers unprecedented sensitivity. We use these data to explore several CO isotopologue line ratios in detail, study their radial (and azimuthal) trends and investigate whether changes in line ratios stem from changes in ISM properties such as gas temperatures, densities or chemical abundances. For example, we find negative radial trends for the $^{13}$CO}/$^{12}$CO, C$^{18}$O/$^{12}$CO and C$^{18}$O/$^{13}$CO line ratios in their J=1-0 transitions. We also find variations with local environment, such as higher $^{12}$CO(2-1)/(1-0) or $^{13}$CO/$^{12}$CO(1-0) line ratios in interarm regions compared to spiral arm regions. We propose that these aforementioned variations of CO line ratios are most likely due to a variation of the optical depth, while abundance variations due to selective nucleosynthesis on a galaxy-wide scale could also play a role. We also study the CO spectral line energy distribution (SLED) using archival JCMT $^{12}$CO(3-2) data and find a variation of the SLED shape with local environmental parameters further underlying changes in optical depth, gas temperatures or densities.

preprint2022arXiv

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

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

preprint2022arXiv

A$^{3}$COSMOS: A census on the molecular gas mass and extent of main-sequence galaxies across cosmic time

To constrain for the first time the mean mass and extent of the molecular gas of a mass-complete sample of $>10^{10}$M$_{\odot}$ main-sequence (MS) galaxies at $0.4<z<3.6$. We apply an innovative $uv$-based stacking analysis to a large set of archival Atacama Large Millimeter/submillimeter Array (ALMA) observations. This stacking analysis provides measurements of the mean mass and extent of the molecular gas of galaxy populations. The molecular gas mass of MS galaxies evolves with redshift and stellar mass. At all stellar masses, the molecular gas fraction decreases by a factor of 24 from $z\sim3.2$ to $z\sim0$. At a given redshift, the molecular gas fraction of MS galaxies decreases with stellar mass, at roughly the same rate as their specific star formation rate decreases. The molecular gas depletion time of MS galaxies remains roughly constant at $z>0.5$ with a value of 300--500 Myr, but increases by a factor of 3 from $z\sim0.5$ to $z\sim0$. This evolution of the molecular gas depletion time of MS galaxies can be predicted from the evolution of their molecular gas surface density and a seemingly universal MS-only $Σ_{M_{\rm mol}}-Σ_{\rm SFR}$ relation with an inferred slope of 1.13, i.e., the so-called KS relation. The far-infrared size of MS galaxies shows no significant evolution with redshift or stellar mass, with a mean circularized half-light radius of 2.2 kpc. Finally, our mean molecular gas masses are lower than previous estimates, likely caused by the fact that literature studies were biased towards individually-detected MS galaxies with massive gas reservoirs. To first order, the molecular gas content of MS galaxies regulates their star formation across cosmic time, while variation of their star formation efficiency plays a secondary role. Despite a large evolution of their gas content and SFRs, MS galaxies evolved along a seemingly universal MS-only KS relation.

preprint2022arXiv

AGN-driven Cold Gas Outflow of NGC 1068 Characterized by Dissociation-Sensitive Molecules

Recent developments in (sub-)millimeter facilities have drastically changed the amount of information obtained from extragalactic spectral scans. In this paper, we present a feature extraction technique using principal component analysis (PCA) applied to arcsecond-resolution (1.0-2.0 arcsec = 72-144 pc) spectral scan datasets for the nearby type-2 Seyfert galaxy, NGC 1068, using Band 3 of the Atacama Large Millimeter/submillimeter Array. We apply PCA to 16 well-detected molecular line intensity maps convolved to a common 150 pc resolution. In addition, we include the [SIII]/[SII] line ratio and [CI] $^3P_1$-$^3P_0$ maps in the literature, both of whose distributions show remarkable resemblance with that of a kpc-scale biconical outflow from the central AGN. We identify two prominent features: (1) central concentration at the circumnuclear disk (CND) and (2) two peaks across the center that coincide with the biconical outflow peaks. The concentrated molecular lines in the CND are mostly high-dipole molecules (e.g., H$^{13}$CN, HC$_3$N, and HCN). Line emissions from molecules known to be enhanced in irradiated interstellar medium, CN, C$_2$H, and HNC, show similar concentrations and extended components along the bicone, suggesting that molecule dissociation is a dominant chemical effect of the cold molecular outflow of this galaxy. Although further investigation should be made, this scenario is consistent with the faintness or absence of the emission lines from CO isotopologues, CH$_3$OH, and N$_2$H$^+$, in the outflow, which are easily destroyed by dissociating photons and electrons.

preprint2022arXiv

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

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

preprint2022arXiv

Evolution of thermal and nonthermal radio continuum emission on kpc scales--Predictions for SKA

Resolved maps of the thermal and nonthermal radio continuum (RC) emission of distant galaxies are a powerful tool for understanding the role of the interstellar medium (ISM) in the evolution of galaxies. We simulate the RC surface brightness of present-day star forming galaxies in the past at $0.15<z<3$ considering two cases of radio size evolution: (1)~no evolution, and (2)~same evolution as in the optical. We aim to investigate the a)~structure of the thermal and nonthermal emission on kpc scales, b)~evolution of the thermal fraction and synchrotron spectrum at mid-radio frequencies ($\simeq$1-10\,GHz), and c)~capability of the proposed SKA1-MID reference surveys in detecting the RC emitting structures. The synchrotron spectrum flattens with $z$ causing curvature in the observed mid-radio SEDs of galaxies at higher $z$. The spectral index reported in recent observational studies agrees better with the no size evolution scenario. In this case, the mean thermal fraction observed at 1.4\,GHz increases with redshift by more than 30\% from $z=0.15$ to $z=2$ because of the drop of the synchrotron emission at higher rest-frame~frequencies. More massive galaxies have lower thermal fractions and experience a faster flattening of the nonthermal spectrum. The proposed SKA1-MID band~2 reference survey, unveils the ISM in M51- and NGC6946-like galaxies (with ${\rm M_{\star}}\simeq10^{10}\,{\rm M}_{\odot}$) up to $z=3$. This survey detects lower-mass galaxies like M33 (${\rm M_{\star}}\simeq10^{9}\,{\rm M}_{\odot}$) only at low redshifts $z\lesssim 0.5$. For a proper separation of the RC emitting processes at the peak of star formation, it is vital to include band~1 into the SKA1-MID reference surveys.

preprint2022arXiv

Planetary Nebula Luminosity Function distances for 19 galaxies observed by PHANGS-MUSE

We provide new planetary nebula luminosity function (PNLF) distances to 19 nearby spiral galaxies that were observed with VLT/MUSE by the PHANGS collaboration. Emission line ratios are used to separate planetary nebulae (PNe) from other bright [OIII] emitting sources like compact supernovae remnants (SNRs) or HII regions. While many studies have used narrowband imaging for this purpose, the detailed spectral line information provided by integral field unit (IFU) spectroscopy grants a more robust way of categorising different [OIII] emitters. We investigate the effects of SNR contamination on the PNLF and find that we would fail to classify all objects correctly, when limited to the same data narrowband imaging provides. However, the few misclassified objects usually do not fall on the bright end of the luminosity function, and only in three cases does the distance change by more than $1σ$. We find generally good agreement with literature values from other methods. Using metallicity constraints that have also been derived from the same IFU data, we revisit the PNLF zero point calibration. Over a range of $8.34 < 12 + \log(\mathrm{O}/\mathrm{H}) < 8.59$, our sample is consistent with a constant zero point and yields $M^*=-4.542^{+0.103}_{-0.059}\, \mathrm{mag}$, within $1σ$ of other literature values. MUSE pushes the limits of PNLF studies and makes galaxies beyond $20\, \mathrm{Mpc}$ accessible for this kind of analysis. This approach to the PNLF shows great promise for leveraging existing archival IFU data on nearby galaxies.

preprint2022arXiv

The Column Densities of Molecular Gas across Cosmic Time: Bridging Observations and Simulations

Observations of the cosmic evolution of different gas phases across time indicate a marked increase in the molecular gas mass density towards $z\sim 2-3$. Such a transformation implies an accompanied change in the global distribution of molecular hydrogen column densities ($N_{\rm{H_2}}$). Using observations by PHANGS-ALMA/SDSS and simulations by GRIFFIN/IllustrisTNG we explore the evolution of this H$_2$ column density distribution function [$f(N_{\rm{H}_2})$]. The H$_2$ (and HI) column density maps for TNG50 and TNG100 are derived in post-processing and are made available through the IllustrisTNG online API. The shape and normalization of $f(N_{\rm{H}_2})$ of individual main-sequence star-forming galaxies are correlated with the star formation rate (SFR), stellar mass (${M_*}$), and H$_2$ mass ($M_{\rm{H}_2}$) in both observations and simulations. TNG100, combined with H$_2$ post-processing models, broadly reproduces observations, albeit with differences in slope and normalization. Also, an analytically modelled $f(N)$, based on exponential gas disks, matches well with the simulations. The GRIFFIN simulation gives first indications that the slope of $f(N_{\rm{H}_2})$ might not majorly differ when including non-equilibrium chemistry in simulations. The $f(N_{\rm{H}_2})$ by TNG100 implies that higher molecular gas column densities are reached at $z=3$ than at $z=0$. Further, denser regions contribute more to the molecular mass density at $z=3$. Finally, H$_2$ starts dominating compared to HI only at column densities above log($N_{\rm{H}_2} / \rm{cm}^{-2}) \sim 21.8-22$ at both redshifts. These results imply that neutral atomic gas is an important contributor to the overall cold gas mass found in the ISM of galaxies including at densities typical for molecular clouds at $z=0$ and $z=3$.

preprint2022arXiv

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

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

preprint2022arXiv

The Kiloparsec-scale Neutral Atomic Carbon Outflow in the Nearby Type-2 Seyfert Galaxy NGC 1068: Evidence for Negative AGN Feedback

Active galactic nucleus (AGN) feedback is postulated as a key mechanism for regulating star formation within galaxies. Studying the physical properties of the outflowing gas from AGN is thus crucial for understanding the co-evolution of galaxies and supermassive black holes. Here we report 55 pc resolution ALMA neutral atomic carbon [CI] $^3P_1\text{-}^3P_0$ observations toward the central 1 kpc of the nearby type-2 Seyfert galaxy NGC 1068, supplemented by 55 pc resolution CO($J=1\text{-}0$) observations. We find that [CI] emission within the central kpc is strongly enhanced by a factor of $>$5 compared to the typical [CI]/CO intensity ratio of $\sim$0.2 for nearby starburst galaxies (in units of brightness temperature). The most [CI]-enhanced gas (ratio $>$ 1) exhibits a kpc-scale elongated structure centered at the AGN that matches the known biconical ionized gas outflow entraining molecular gas in the disk. A truncated, decelerating bicone model explains well the kinematics of the elongated structure, indicating that the [CI] enhancement is predominantly driven by the interaction between the ISM in the disk and the highly inclined ionized gas outflow (which is likely driven by the radio jet). Our results strongly favor the "CO dissociation scenario" rather than the "in-situ C formation" one which prefers a perfect bicone geometry. We suggest that the high [CI]/CO intensity ratio gas in NGC 1068 directly traces ISM in the disk that is currently dissociated and entrained by the jet and the outflow, i.e., the "negative" effect of the AGN feedback.

preprint2022arXiv

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

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

preprint2021arXiv

CO excitation, molecular gas density and interstellar radiation field in local and high-redshift galaxies

We study the Carbon Monoxide (CO) excitation, mean molecular gas density and interstellar radiation field (ISRF) intensity in a comprehensive sample of 76 galaxies from local to high redshift (z~0-6), selected based on detections of their CO transitions J=2-1 and 5-4 and their optical/infrared/(sub-)millimeter spectral energy distributions (SEDs). We confirm the existence of a tight correlation between CO excitation as traced by the CO(5-4)/(2-1) line ratio (R52), and the mean ISRF intensity U as derived from infrared SED fitting using dust SED templates. By modeling the molecular gas density probability distribution function (PDF) in galaxies and predicting CO line ratios with large velocity gradient radiative transfer calculations, we present a framework linking global CO line ratios to the mean molecular hydrogen gas density nH2 and kinetic temperature Tkin. Mapping in this way observed R52 ratios to nH2 and Tkin probability distributions, we obtain positive U-nH2 and U-Tkin correlations, which imply a scenario in which the ISRF in galaxies is mainly regulated by Tkin and (non-linearly) by nH2. A small fraction of starburst galaxies showing enhanced nH2 could be due to merger-driven compaction. Our work demonstrates that ISRF and CO excitation are tightly coupled, and that density-PDF modeling is a promising tool for probing detailed ISM properties inside galaxies.

preprint2021arXiv

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

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

preprint2021arXiv

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

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

preprint2021arXiv

Molecular Gas Excitation of the Massive Dusty Starburst CRLE and the Main-Sequence Galaxy HZ10 at z=5.7 in the COSMOS Field

We report CO(5$\rightarrow$4) and CO(6$\rightarrow$5) line observations in the dusty starbursting galaxy CRLE ($z = 5.667$) and the main-sequence (MS) galaxy HZ10 ($z = 5.654$) with the Northern Extended Millimeter Array (NOEMA). CRLE is the most luminous $z>5$ starburst in the COSMOS field and HZ10 is the most gas-rich "normal" galaxy currently known at $z>5$. We find line luminosities for CO(5$\rightarrow$4) and CO(6$\rightarrow$5) of (4.9 $\pm$ 0.5) and (3.8 $\pm$ 0.4) $\times$ 10$^{10}$ K km s$^{-1}$ pc$^{2}$ for CRLE and upper limits of $< 0.76$ and $< 0.60$ $\times$ 10$^{10}$ K km s$^{-1}$ pc$^{2}$ for HZ10, respectively. The CO excitation of CRLE appears comparable to other $z>5$ dusty star-forming galaxies (DSFGs). For HZ10, these line luminosity limits provide the first significant constraints of this kind for a MS galaxy at $z > 5$. We find the upper limit of $L'_{5\rightarrow4}/L'_{2\rightarrow1}$ in HZ10 could be similar to the average value for MS galaxies around $z\approx 1.5$, suggesting that MS galaxies with comparable gas excitation may already have existed one billion years after the Big Bang. For CRLE we determine the most likely values for the H$_2$ density, kinetic temperature and dust temperature based on excitation modeling of the CO line ladder. We also derive a total gas mass of $(7.1 \pm 1.3) \times 10^{10} M_\odot$. Our findings provide some of the currently most detailed constraints on the gas excitation that sets the conditions for star formation in a galaxy protocluster environment at $z > 5$.

preprint2021arXiv

Molecular Gas Properties and CO-to-H2 Conversion Factors in the Central Kiloparsec of NGC 3351

The CO-to-H$_2$ conversion factor ($α_\rm{CO}$) is critical to studying molecular gas and star formation in galaxies. The value of $α_\rm{CO}$ has been found to vary within and between galaxies, but the specific environmental conditions that cause these variations are not fully understood. Previous observations on $\sim$kpc scales revealed low values of $α_\rm{CO}$ in the centers of some barred spiral galaxies, including NGC 3351. We present new ALMA Band 3, 6, and 7 observations of $^{12}$CO, $^{13}$CO, and C$^{18}$O lines on 100 pc scales in the inner $\sim$2 kpc of NGC 3351. Using multi-line radiative transfer modeling and a Bayesian likelihood analysis, we infer the H$_2$ density, kinetic temperature, CO column density per line width, and CO isotopologue abundances on a pixel-by-pixel basis. Our modeling implies the existence of a dominant gas component with a density of $2{-}3\times10^3$ $\rm{cm^{-3}}$ in the central ${\sim}$1 kpc and a high temperature of 30$-$60 K near the nucleus and near the contact points that connect to the bar-driven inflows. Assuming a CO/H$_2$ abundance of $3\times10^{-4}$, our analysis yields $α_\rm{CO}{\sim}0.5{-}2.0$ $\rm{M_\odot\,(K~km~s^{-1}~pc^2)^{-1}}$ with a decreasing trend with galactocentric radius in the central $\sim$1 kpc. The inflows show a substantially lower $α_\rm{CO} < 0.1$ $\rm{M_\odot\,(K~km~s^{-1}~pc^2)^{-1}}$, likely due to lower optical depths caused by turbulence or shear in the inflows. Over the whole region, this gives an intensity-weighted $α_\rm{CO}$ of ${\sim}1.5$ $\rm{M_\odot\,(K~km~s^{-1}~pc^2)^{-1}}$, which is similar to previous dust modeling based results at kpc scales. This suggests that low $α_\rm{CO}$ on kpc scales in the centers of some barred galaxies may be due to the contribution of low optical depth CO emission in bar-driven inflows.

preprint2021arXiv

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

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

preprint2021arXiv

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

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

preprint2021arXiv

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

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

preprint2020arXiv

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

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

preprint2020arXiv

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

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.

preprint2020arXiv

ALMA reveals the molecular gas properties of 5 star-forming galaxies across the main sequence at 3 < z < 3.5

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

preprint2020arXiv

Bringing high spatial resolution to the Far-infrared -- A giant leap for astrophysics

The far-infrared (FIR) regime is one of the few wavelength ranges where no astronomical data with sub-arcsecond spatial resolution exist. Neither of the medium-term satellite projects like SPICA, Millimetron nor O.S.T. will resolve this malady. For many research areas, however, information at high spatial and spectral resolution in the FIR, taken from atomic fine-structure lines, from highly excited carbon monoxide (CO), light hydrids, and especially from water lines would open the door for transformative science. A main theme will be to trace the role of water in proto-planetary disks, to observationally advance our understanding of the planet formation process and, intimately related to that, the pathways to habitable planets and the emergence of life. Furthermore, key observations will zoom into the physics and chemistry of the star-formation process in our own Galaxy, as well as in external galaxies. The FIR provides unique tools to investigate in particular the energetics of heating, cooling and shocks. The velocity-resolved data in these tracers will reveal the detailed dynamics engrained in these processes in a spatially resolved fashion, and will deliver the perfect synergy with ground-based molecular line data for the colder dense gas.

preprint2020arXiv

Distances to PHANGS Galaxies: New Tip of the Red Giant Branch Measurements and Adopted Distances

PHANGS-HST is an ultraviolet-optical imaging survey of 38 spiral galaxies within ~20 Mpc. Combined with the PHANGS-ALMA, PHANGS-MUSE surveys and other multiwavelength data, the dataset will provide an unprecedented look into the connections between young stars, HII regions, and cold molecular gas in these nearby star-forming galaxies. Accurate distances are needed to transform measured observables into physical parameters (e.g., brightness to luminosity, angular to physical sizes of molecular clouds, star clusters and associations). PHANGS-HST has obtained parallel ACS imaging of the galaxy halos in the F606W and F814W bands. Where possible, we use these parallel fields to derive tip of the red giant branch (TRGB) distances to these galaxies. In this paper, we present TRGB distances for 11 galaxies from ~4 to ~15 Mpc, based on the first year of PHANGS-HST observations. Five of these represent the first published TRGB distance measurements (IC 5332, NGC 2835, NGC 4298, NGC 4321, and NGC 4328), and eight of which are the best available distances to these targets. We also provide a compilation of distances for the 118 galaxies in the full PHANGS sample, which have been adopted for the first PHANGS-ALMA public data release.

preprint2020arXiv

Dynamical Equilibrium in the Molecular ISM in 28 Nearby Star-Forming Galaxies

We compare the observed turbulent pressure in molecular gas, $P_\mathrm{turb}$, to the required pressure for the interstellar gas to stay in equilibrium in the gravitational potential of a galaxy, $P_\mathrm{DE}$. To do this, we combine arcsecond resolution CO data from PHANGS-ALMA with multi-wavelength data that traces the atomic gas, stellar structure, and star formation rate (SFR) for 28 nearby star-forming galaxies. We find that $P_\mathrm{turb}$ correlates with, but almost always exceeds the estimated $P_\mathrm{DE}$ on kiloparsec scales. This indicates that the molecular gas is over-pressurized relative to the large-scale environment. We show that this over-pressurization can be explained by the clumpy nature of molecular gas; a revised estimate of $P_\mathrm{DE}$ on cloud scales, which accounts for molecular gas self-gravity, external gravity, and ambient pressure, agrees well with the observed $P_\mathrm{turb}$ in galaxy disks. We also find that molecular gas with cloud-scale ${P_\mathrm{turb}}\approx{P_\mathrm{DE}}\gtrsim{10^5\,k_\mathrm{B}\,\mathrm{K\,cm^{-3}}}$ in our sample is more likely to be self-gravitating, whereas gas at lower pressure appears more influenced by ambient pressure and/or external gravity. Furthermore, we show that the ratio between $P_\mathrm{turb}$ and the observed SFR surface density, $Σ_\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 $Σ_\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.

preprint2020arXiv

Pa$β$, H$α$, and Attenuation in NGC 5194 and NGC 6946

We combine Hubble Space Telescope (HST) Paschen $β$ (Pa$β$) imaging with ground-based, previously published H$α$ maps to estimate the attenuation affecting H$α$, A(H$α$), across the nearby, face-on galaxies NGC 5194 and NGC 6946. We estimate A(H$α$) in ~ 2,000 independent 2" ~75 pc diameter apertures in each galaxy, spanning out to a galactocentric radius of almost 10 kpc. In both galaxies, A(H$α$) drops with radius, with a bright, high attenuation inner region, though in detail the profiles differ between the two galaxies. Regions with the highest attenuation-corrected H$α$ luminosity show the highest attenuation, but the observed H$α$ luminosity of a region is not a good predictor of attenuation in our data. Consistent with much previous work, the IR-to-H$α$ color does a good job of predicting A(H$α$). We calculate the best-fit empirical coefficients for use combining H$α$ with 8, 12, 24, 70, or 100 $μ$m to correct for attenuation. These agree well with previous work but we also measure significant scatter around each of these linear relations. The local atomic plus molecular gas column density, N(H), also predicts A(H$α$) well. We show that a screen with magnitude ~ 0.2 times the expected for a Milky Way gas-to-dust value does a reasonable job of explaining A(H$α$) as a function of N(H). This could be expected if only ~ 40% of gas and dust directly overlap regions of H$α$ emission.

preprint2020arXiv

PHANGS CO kinematics: disk orientations and rotation curves at 150 pc resolution

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.

preprint2020arXiv

The 300-pc Scale ALMA View of [CI] $^3P_1$-$^3P_0$, CO $J$=1-0, and 609 $μ$m Dust Continuum in A Luminous Infrared Galaxy

We present high-quality ALMA Band 8 observations of the [CI] $^3P_1$-$^3P_0$ line and 609 $μ$m dust continuum emission toward the nearby luminous infrared galaxy (LIRG) IRAS F18293-3413, as well as matched resolution (300-pc scale) Band 3 CO $J=$1-0 data, which allow us to assess the use of the [CI] $^3P_1$-$^3P_0$ line as a total gas mass estimator. We find that the [CI] line basically traces structures detected in CO (and dust), and a mean (median) [CI]/CO luminosity ($L'_{\rm [CI]}$/$L'_{\rm CO}$) ratio of 0.17 (0.16) with a scatter of 0.04. However, a pixel-by-pixel comparison revealed that there is a radial $L'_{\rm [CI]}$/$L'_{\rm CO}$ gradient and a superlinear $L'_{\rm CO}$ vs. $L'_{\rm [CI]}$ relation (slope = 1.54 $\pm$ 0.02) at this spatial scale, which can be explained by radial excitation and/or line opacity gradients. Based on the molecular gas masses converted from the dust continuum emission, we found that the CO-to-H$_2$ and [CI]-to-H$_2$ conversion factors are relatively flat across the molecular gas disk with a median value of 3.5$^{+1.9}_{-1.3}$ and 20.7$^{+9.2}_{-4.9}$ $M_{\odot}$ (K km s$^{-1}$ pc$^2$)$^{-1}$, respectively. A non-LTE calculation yields that typical molecular gas properties seen in nearby (U)LIRGs ($n_{\rm H_2}$ = 10$^{3-4}$ cm$^{-3}$, $T_{\rm kin}$ $\sim$ 50 K, and $X_{\rm CI}$ = (0.8-2.3) $\times$ 10$^{-5}$) can naturally reproduce the derived [CI]-to-H$_2$ conversion factor. However, we caution that a careful treatment of the physical gas properties is required in order to measure H$_2$ gas mass distributions in galaxies using a single [CI] line. Otherwise, a single [CI] line is not a good molecular gas estimator in a spatially resolved manner.

preprint2020arXiv

The VLA-COSMOS 3 GHz Large Project: Evolution of specific star formation rates out to $z\sim5$

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

preprint2020arXiv

Volumetric Star Formation Prescriptions in Vertically Resolved Edge-on Galaxies

We measure the gas disc thicknesses of the edge-on galaxy NGC 4013 and the less edge-on galaxies (NGC 4157 and 5907) using CO (CARMA/OVRO) and/or HI (EVLA) observations. We also estimate the scale heights of stars and/or the star formation rate (SFR) for our sample of five galaxies using Spitzer IR data (3.6 $μ$m and 24 $μ$m). We derive the average volume densities of the gas and the SFR using the measured scale heights along with radial surface density profiles. Using the volume density that is more physically relevant to the SFR than the surface density, we investigate the existence of a volumetric star formation law (SFL), how the volumetric SFL is different from the surface-density SFL, and how the gas pressure regulates the SFR based on our galaxy sample. We find that the volumetric and surface SFLs in terms of the total gas have significantly different slopes, while the volumetric and surface SFLs in terms of the molecular gas do not show any noticeable difference. The volumetric SFL for the total gas has a flatter power-law slope of 1.26 with a smaller scatter of 0.19 dex compared to the slope (2.05) and the scatter (0.25 dex) of the surface SFL. The molecular gas SFLs have similar slopes of 0.78 (volume density) and 0.77 (surface density) with the same rms scatter. We show that the interstellar gas pressure is strongly correlated with the SFR but find no significant difference between the correlations based on the volume and surface densities.

preprint2020arXiv

When Gas Dynamics Decouples from Galactic Rotation: Characterizing ISM Circulation in Disk Galaxies

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, $Γ$, 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 $λ_c$. We define the spatial scale $λ_{\rm eq}$ at which galactic rotation and non-circular motions contribute equally to $Γ$. 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 $λ_c$ and this transition is necessary to reproduce the circulation distribution. We find that $λ_{\rm eq}$, and therefore the structure of the gas velocity field, is set by the local conditions of gravitational stability and stellar feedback.

preprint2019arXiv

Mapping electron temperature variations across a spiral arm in NGC 1672

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

preprint2019arXiv

The Gas-Star Formation Cycle in Nearby Star-Forming Galaxies I. Assessment of Multi-scale Variations

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.

preprint2019arXiv

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

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

preprint2019arXiv

The lifecycle of molecular clouds in nearby star-forming disc galaxies

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$α$ 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 $Σ_{\rm H_2}\geqslant8$M$_{\odot}$pc$^{-2}$, the GMC lifetime correlates with time-scales for galactic dynamical processes, whereas at $Σ_{\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$α$ (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.

preprint2018arXiv

Cloud-Scale Molecular Gas Properties in 15 Nearby Galaxies

We measure the velocity dispersion, $σ$, and surface density, $Σ$, 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. $Σ$ and $σ$ 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 $α_\mathrm{vir}\proptoσ^2/Σ$, 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Σ\,σ^2$, which spans ${\sim}5\rm\;dex$ across our sample. We find $Σ$, $σ$, 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 $σ$-$Σ$ relation flattens at $Σ\lesssim10\rm\;M_\odot\,pc^{-2}$, leading to high $σ$ for a given $Σ$ and high apparent $α_\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.

preprint2016arXiv

(Sub)millimetre interferometric imaging of a sample of COSMOS/AzTEC submillimetre galaxies IV. Physical properties derived from spectral energy distributions

We characterise the physical nature of a 1.1 mm-selected, flux-limited, and interferometrically followed up sample of SMGs in COSMOS. We used the MAGPHYS code to fit the multiwavelength (UV-radio) SEDs of 16 of the target SMGs. We also constructed the pure radio SEDs of our SMGs using three different radio bands (325 MHz, 1.4 GHz, and 3 GHz). Moreover, since two SMGs in our sample, AzTEC1 and AzTEC3, benefit from previous CO line observations, we studied their properties in more detail. We found that 63% of our target SMGs lie above the galaxy main-sequence by more than a factor of 3, and hence are starbursts. The 3 GHz radio sizes we have previously measured for the target SMGs were compared with the present stellar mass estimates, and we found that the z>3 SMGs are fairly consistent with the mass-size relationship of z~2 compact, quiescent galaxies (cQGs). The median IR-radio correlation parameter is found to be q=2.27, which is lower than measured locally (median q=2.64). AzTEC1 is found to have a sub-Eddington SFR surface density (by a factor of 2.6), while AzTEC3 appears to be an Eddington-limited starburster. The gas reservoir in these two high-z SMGs would be exhausted in only ~86 and 19 Myr at the current SFR, respectively. A comparison of the MAGPHYS-based properties of our SMGs with those of equally bright ALESS SMGs suggests that the two populations share fairly similar physical characteristics, including the q parameter. A hint of negative correlation is found between the 3 GHz size and the level of starburstiness, and hence cosmic-ray electrons in more compact starbursts might be more susceptible to free-free absorption. Some of the derived low and high q values (compared to the local median) could be the result of a specific merger/post-starburst phase of galaxy evolution. Overall, our results support the scenario where z>3 SMGs evolve into today's giant ellipticals.

preprint2016arXiv

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

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.

preprint2016arXiv

Being WISE II: Reducing the Influence of Star Formation History on the Mass-to-Light Ratio of Quiescent Galaxies

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$μ$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 $μ$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$\%$.

preprint2015arXiv

(Sub)millimetre interferometric imaging of a sample of COSMOS/AzTEC submillimetre galaxies - I. Multiwavelength identifications and redshift distribution

We used the PdBI to map a sample of 15 SMGs in the COSMOS field at the wavelength of 1.3 mm. The target SMGs were originally discovered in the JCMT/AzTEC 1.1 mm continuum survey at S/N=4-4.5. This paper presents, for the first time, interferometric millimetre-wavelength observations of these sources. The angular resolution of our observations, 1.8", allowed us to accurately determine the positions of the target SMGs. Using a detection threshold of S/N>4.5 regardless of multiwavelength counterpart association, and 4<S/N<=4.5 if a multiwavelength counterpart within 1.5" is also present, the total number of detections in our survey is 22. Three of our detected SMGs (AzTEC21, 27, and 28; which corresponds to 20%) are marginally resolved at our angular resolution, and these sources are found to have elongated or clumpy morphologies and/or multiple components. Using optical to NIR photometric redshifts, available spectroscopic redshifts, and redshifts estimated from the radio-to-submm spectral index we infer a median redshift of $\tilde{z}=3.20\pm0.25$ for our sample. To study the overall multiplicity and redshift distribution of flux-limited samples of SMGs we combined these sources with the 15 brightest AzTEC SMGs detected at 1.1 mm, AzTEC1-15, and studied previously. This constitutes a complete, flux- and S/N-limited 1.1-mm selected sample. We find that the median redshift for the 15 brightest AzTEC SMGs ($\tilde{z}=3.05\pm0.44$) is consistent with that for AzTEC16-30. This conforms to recent observational findings that SMGs do not exhibit any significant trend between the redshift and (sub)mm flux density. For the combined AzTEC1-30 sample we derive a median redshift of $\tilde{z}=3.17\pm0.27$, consistent with previous results based on mm-selected samples. We further infer that within the combined AzTEC1-30 sample $\sim25\pm9\%$ of sources separate into multiple components.

preprint2015arXiv

(Sub)millimetre interferometric imaging of a sample of COSMOS/AzTEC submillimetre galaxies - II. The spatial extent of the radio-emitting regions

Radio emission at cm wavelengths from highly star-forming galaxies, such as SMGs, is dominated by synchrotron radiation arising from supernova activity. Using deep, high-resolution ($1σ=2.3$ $μ$Jy beam$^{-1}$; $0.75^{"}$) cm radio-continuum observations taken by the VLA-COSMOS 3 GHz Large Project, we studied the radio-emitting sizes of a flux-limited sample of SMGs in the COSMOS field. Of the 39 SMGs studied here, 3 GHz emission was detected towards 18 of them ($\sim46\pm11\%$) with S/N ratios in the range of ${\rm S/N=4.2-37.4}$. Using 2D elliptical Gaussian fits, we derived a median deconvolved major axis FWHM size of $0.54^{"}\pm 0.11^{"}$ for our 18 SMGs detected at 3 GHz. For the 15 SMGs with known redshift we derived a median linear major axis FWHM of $4.2\pm0.9$ kpc. No clear correlation was found between the radio-emitting size and the 3 GHz or submm flux density, or the redshift of the SMG. However, there is a hint of larger radio sizes at $z\sim2.5-5$ compared to lower redshifts. The sizes we derived are consistent with previous SMG sizes measured at 1.4 GHz and in mid-$J$ CO emission, but significantly larger than those seen in the (sub)mm continuum emission. One possible scenario is that SMGs have i) an extended gas component with a low dust temperature, and which can be traced by low- to mid-$J$ CO line emission and radio continuum emission, and ii) a warmer, compact starburst region giving rise to the high-$J$ line emission of CO, which could dominate the dust continuum size measurements. Because of the rapid cooling of CR electrons in dense starburst galaxies ($\sim10^4-10^5$ yr), the more extended synchrotron radio-emitting size being a result of CR diffusion seems unlikely. Instead, if SMGs are driven by galaxy mergers the radio synchrotron emission might arise from more extended magnetised ISM around the starburst region.

preprint2015arXiv

ALMA Observations of the Submillimeter Dense Molecular Gas Tracers in the Luminous Type-1 Active Nucleus of NGC 7469

We present ALMA Cycle 1 observations of the central kpc region of the luminous type-1 Seyfert galaxy NGC 7469 with unprecedented high resolution (0.5$"$ $\times$ 0.4$"$ = 165 pc $\times$ 132 pc) at submillimeter wavelengths. Utilizing the wide-bandwidth of ALMA, we simultaneously obtained HCN(4-3), HCO$^+$(4-3), CS(7-6), and partially CO(3-2) line maps, as well as the 860 $μ$m continuum. The region consists of the central $\sim$ 1$"$ component and the surrounding starburst ring with a radius of $\sim$ 1.5$"$-2.5$"$. Several structures connect these components. Except for CO(3-2), these dense gas tracers are significantly concentrated towards the central $\sim$ 1$"$, suggesting their suitability to probe the nuclear regions of galaxies. Their spatial distribution resembles well those of centimeter and mid-infrared continuum emissions, but it is anti-correlated with the optical one, indicating the existence of dust obscured star formation. The integrated intensity ratios of HCN(4-3)/HCO$^+$(4-3) and HCN(4-3)/CS(7-6) are higher at the AGN position than at the starburst ring, which is consistent to our previous findings (submm-HCN enhancement). However, the HCN(4-3)/HCO$^+$(4-3) ratio at the AGN position of NGC 7469 (1.11$\pm$0.06) is almost half of the corresponding value of the low-luminosity type-1 Seyfert galaxy NGC 1097 (2.0$\pm$0.2), despite the more than two orders of magnitude higher X-ray luminosity of NGC 7469. But the ratio is comparable to that of the close vicinity of the AGN of NGC 1068 ($\sim$ 1.5). Based on these results, we speculate that some other heating mechanisms than X-ray (e.g., mechanical heating due to AGN jet) can contribute significantly for shaping the chemical composition in NGC 1097.

preprint2015arXiv

Discovery of Massive, Mostly Star-formation Quenched Galaxies with Extremely Large Lyman-alpha Equivalent Widths at z ~ 3

We report a discovery of 6 massive galaxies with both extremely large Lya equivalent width and evolved stellar population at z ~ 3. These MAssive Extremely STrong Lya emitting Objects (MAESTLOs) have been discovered in our large-volume systematic survey for strong Lya emitters (LAEs) with twelve optical intermediate-band data taken with Subaru/Suprime-Cam in the COSMOS field. Based on the SED fitting analysis for these LAEs, it is found that these MAESTLOs have (1) large rest-frame equivalent width of EW_0(Lya) ~ 100--300 A, (2) M_star ~ 10^10.5--10^11.1 M_sun, and (3) relatively low specific star formation rates of SFR/M_star ~ 0.03--1 Gyr^-1. Three of the 6 MAESTLOs have extended Ly$α$ emission with a radius of several kpc although they show very compact morphology in the HST/ACS images, which correspond to the rest-frame UV continuum. Since the MAESTLOs do not show any evidence for AGNs, the observed extended Lya emission is likely to be caused by star formation process including the superwind activity. We suggest that this new class of LAEs, MAESTLOs, provides a missing link from star-forming to passively evolving galaxies at the peak era of the cosmic star-formation history.

preprint2015arXiv

Microwave Continuum Emission and Dense Gas Tracers in NGC 3627: Combining Jansky VLA and ALMA Observations

We present Karl G. Jansky Very Large Array (VLA) Ka band (33 GHz) and Atacama Large Millimeter Array (ALMA) Band 3 (94.5 GHz) continuum images covering the nucleus and two extranuclear star-forming regions within the nearby galaxy NGC 3627 (M 66), observed as part of the Star Formation in Radio Survey (SFRS). Both images achieve an angular resolution of $\lesssim$2\arcsec, allowing us to map the radio spectral indices and estimate thermal radio fractions at a linear resolution of $\lesssim$90 pc at the distance of NGC 3627. The thermal fraction at 33 GHz reaches unity at and around the peaks of each HII region; we additionally observed the spectral index between 33 and 94.5 GHz to become both increasingly negative and positive away from the peaks of the HII regions, indicating an increase of non-thermal extended emission from diffusing cosmic-ray electrons and the possible presence of cold dust, respectively. While the ALMA observations were optimized for collecting continuum data, they also detected line emission from the $J=1\rightarrow0$ transitions of HCN and HCO$^{+}$. The peaks of dense molecular gas traced by these two spectral lines are spatially offset from the peaks of the 33 and 94.5 GHz continuum emission for the case of the extranuclear star-forming regions, indicating that our data reach an angular resolution at which one can spatially distinguish sites of recent star formation from the sites of future star formation. Finally, we find trends of decreasing dense gas fraction and velocity dispersion with increasing star formation efficiency among the three regions observed, indicating that the dynamical state of the dense gas, rather than its abundance, plays a more significant role in the star formation process.

preprint2015arXiv

Next Generation Very Large Array Memo No. 7 Science Working Group 2: "Galaxy Ecosystems": The Matter Cycle in and Around Galaxies

This white paper discusses how a "next-generation" Very Large Array (ngVLA) operating in the frequency range 1-116 GHz could be a groundbreaking tool to study the detailed astrophysics of the "matter cycle" in the Milky Way and other galaxies. If optimized for high brightness sensitivity, the ngVLA would bring detailed microwave spectroscopy and modeling of the full radio spectral energy distribution into regular use as survey tools at resolutions of 0.1- 1 arcseconds. This wavelength range includes powerful diagnostics of density, excitation, and chemistry in the cold ISM, as well as multiple tracers of the rate of recent star formation, the magnetic field, shocks, and properties of the ionized ISM. We highlight design considerations that would make this facility revolutionary in this area, the foremost of which is a large amount of collecting area on ~km-length baselines. We also emphasize the strong case for harnessing the large proposed collecting area of the ngVLA for very long baseline applications as part of the core design. This would allow measurements of the three dimensional space motions of galaxies to beyond the Local Group and mapping of the Milky Way out to the far side of the disk. Finally, we discuss the gains from the proposed combination of very high resolution and sensitivity to thermal emission, which include observing the feeding of black holes and resolving forming protoclusters.

preprint2015arXiv

Short GMC lifetimes: an observational estimate with the PdBI Arcsecond Whirlpool Survey (PAWS)

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.

preprint2015arXiv

Structure, dynamical impact and origin of magnetic fields in nearby galaxies in the SKA era

Magnetic fields are an important ingredient of the interstellar medium (ISM). Besides their importance for star formation, they govern the transport of cosmic rays, relevant to the launch and regulation of galactic outflows and winds, which in turn are pivotal in shaping the structure of halo magnetic fields. Mapping the small-scale structure of interstellar magnetic fields in many nearby galaxies is crucial to understand the interaction between gas and magnetic fields, in particular how gas flows are affected. Elucidation of the magnetic role in, e.g., triggering star formation, forming and stabilising spiral arms, driving outflows, gas heating by reconnection and magnetising the intergalactic medium has the potential to revolutionise our physical picture of the ISM and galaxy evolution in general. Radio polarisation observations in the very nearest galaxies at high frequencies (>= 3 GHz) and with high spatial resolution (<= 5") hold the key here. The galaxy survey with SKA1 that we propose will also be a major step to understand the galactic dynamo, which is important for models of galaxy evolution and for astrophysical magnetohydrodynamics in general. Field amplification by turbulent gas motions, which is crucial for efficient dynamo action, has been investigated so far only in simulations, while compelling evidence of turbulent fields from observations is still lacking.

preprint2015arXiv

Submillimeter-HCN Diagram for an Energy Diagnostics in the Centers of Galaxies

Compiling data from literature and the ALMA archive, we show enhanced HCN(4-3)/HCO$^+$(4-3) and/or HCN(4-3)/CS(7-6) integrated intensity ratios in circumnuclear molecular gas around active galactic nuclei (AGNs) compared to those in starburst (SB) galaxies (submillimeter HCN-enhancement). The number of sample galaxies is significantly increased from our previous work. We expect this feature could potentially be an extinction-free energy diagnostic tool of nuclear regions of galaxies. Non-LTE radiative transfer modelings of the above molecular emission lines involving both collisional and radiative excitation, as well as a photon trapping effect were conducted to investigate the cause of the high line ratios in AGNs. As a result, we found that enhanced abundance ratios of HCN-to-HCO$^+$ and HCN-to-CS in AGNs as compared to SB galaxies by a factor of a few to even $>$ 10 is a plausible explanation for the submillimeter HCN-enhancement. However, a counter argument of a systematically higher gas density in AGNs than in SB galaxies can also be a plausible scenario. Although we could not fully discriminate these two scenarios at this moment due to insufficient amount of multi-transition, multi-species data, the former scenario equivalently claims for abnormal chemical composition in AGNs. Regarding the actual mechanism to realize the composition, we suggest it is difficult with conventional gas phase X-ray dominated region (XDR) ionization models to reproduce the observed high line ratios. We might have to take into account other mechanisms such as neutral-neutral reactions that are efficiently activated at high temperature environments and/or mechanically heated regions to further understand the high line ratios in AGNs.

preprint2015arXiv

Testing the Molecular-Hydrogen Kennicutt-Schmidt Law in the Low-Density Environments of Extended Ultraviolet Disk Galaxies

Studying star formation beyond the optical radius of galaxies allows us to test empirical relations in extreme conditions with low average gas density and low molecular fraction. Previous studies discovered galaxies with extended ultraviolet (XUV) disks, which often contain star forming regions with lower Halpha-to-far-UV (FUV) flux ratios compared to inner disk star forming regions. However, most previous studies lack measurements of molecular gas, which is presumably the component of the interstellar medium out of which stars form. We analyzed published CO measurements and upper limits for fifteen star forming regions in the XUV or outer disk of three nearby spiral galaxies and a new CO upper limit from the IRAM 30 m telescope in one star forming region at r = 3.4 r_25 in the XUV disk of NGC 4625. We found that the star forming regions are in general consistent with the same molecular-hydrogen Kennicutt-Schmidt law that applies within the optical radius, independent of whether we used Halpha or FUV as the star formation rate (SFR) tracer. However, a number of the CO detections are significantly offset towards higher SFR surface density for their molecular hydrogen surface density. Deeper CO data may enable us to use the presence or absence of molecular gas as an evolutionary probe to break the degeneracy between age and stochastic sampling of the initial mass function as the explanation for the low Halpha-to-FUV flux ratios in XUV disks.

preprint2015arXiv

The Odd Offset between the Galactic Disk and Its Bar in NGC 3906

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.

preprint2015arXiv

The Spitzer Survey of Stellar Structure in Galaxies (S4G): Stellar Masses, Sizes and Radial Profiles for 2352 Nearby Galaxies

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

preprint2015arXiv

Understanding the central kinematics of globular clusters with simulated integrated-light IFU observations

The detection of intermediate mass black holes in the centres of globular clusters is highly controversial, as complementary observational methods often deliver significantly different results. In order to understand these discrepancies, we develop a procedure to simulate integral field unit (IFU) observations of globular clusters: Simulating IFU Star Cluster Observations (SISCO). The input of our software are realistic dynamical models of globular clusters that are then converted in a spectral data cube. We apply SISCO to Monte Carlo cluster simulations from Downing et al. (2010), with a realistic number of stars and concentrations. Using independent realisations of a given simulation we are able to quantify the stochasticity intrinsic to the problem of observing a partially resolved stellar population with integrated-light spectroscopy. We show that the luminosity-weighted IFU observations can be strongly biased by the presence of a few bright stars that introduce a scatter in the velocity dispersion measurements up to $\simeq$40% around the expected value, preventing any sound assessment of the central kinematic and a sensible interpretation of the presence/absence of an intermediate mass black hole. Moreover, we illustrate that, in our mock IFU observations, the average kinematic tracer has a mass of $\simeq$0.75 solar masses, only slightly lower than the mass of the typical stars examined in studies of resolved line-of-sight velocities of giant stars. Finally, in order to recover unbiased kinematic measurements we test different masking techniques that allow us to remove the spaxels dominated by bright stars, bringing the scatter down to a level of only a few percent. The application of SISCO will allow to investigate state-of-the-art simulations as realistic observations.

preprint2014arXiv

A New Method for Measuring Metallicities of Young Super Star Clusters

We demonstrate how the metallicities of young super star clusters can be measured using novel spectroscopic techniques in the J-band. The near-infrared flux of super star clusters older than ~6 Myr is dominated by tens to hundreds of red supergiant stars. Our technique is designed to harness the integrated light of that population and produces accurate metallicities for new observations in galaxies above (M83) and below (NGC 6946) solar metallicity. In M83 we find [Z]= +0.28 +/- 0.14 dex using a moderate resolution (R~3500) J-band spectrum and in NGC 6496 we report [Z]= -0.32 +/- 0.20 dex from a low resolution spectrum of R~1800. Recently commissioned low resolution multiplexed spectrographs on the VLT (KMOS) and Keck (MOSFIRE) will allow accurate measurements of super star cluster metallicities across the disks of star-forming galaxies up to distances of 70 Mpc with single night observation campaigns using the method presented in this letter.

preprint2014arXiv

ALMA Imaging of Gas and Dust in a Galaxy Protocluster at Redshift 5.3: [CII] Emission in "Typical" Galaxies and Dusty Starbursts ~1 Billion Years after the Big Bang

We report interferometric imaging of [CII] and OH emission toward the center of the galaxy protocluster associated with the z=5.3 submillimeter galaxy (SMG) AzTEC-3, using the Atacama Large (sub)Millimeter Array (ALMA). We detect strong [CII], OH, and rest-frame 157.7 um continuum emission toward the SMG. The [CII] emission is distributed over a scale of 3.9 kpc, implying a dynamical mass of 9.7 x 10^10 Msun, and a star formation rate (SFR) surface density of Sigma_SFR = 530 Msun/yr/kpc2. This suggests that AzTEC-3 forms stars at Sigma_SFR approaching the Eddington limit for radiation pressure supported disks. We find that the OH emission is slightly blueshifted relative to the [CII] line, which may indicate a molecular outflow associated with the peak phase of the starburst. We also detect and dynamically resolve [CII] emission over a scale of 7.5 kpc toward a triplet of Lyman-break galaxies with moderate UV-based SFRs in the protocluster at ~95kpc projected distance from the SMG. These galaxies are not detected in the continuum, suggesting far-infrared SFRs of <18-54 Msun/yr, consistent with a UV-based estimate of 22 Msun/yr. The spectral energy distribution of these galaxies is inconsistent with nearby spiral and starburst galaxies, but resembles those of dwarf galaxies. This is consistent with expectations for young starbursts without significant older stellar populations. This suggests that these galaxies are significantly metal-enriched, but not heavily dust-obscured, "normal" star-forming galaxies at z>5, showing that ALMA can detect the interstellar medium in "typical" galaxies in the very early universe.

preprint2014arXiv

An ALMA Survey of Submillimeter Galaxies in the Extended Chandra Deep Field South: The Redshift Distribution and Evolution of Submillimeter Galaxies

We present the first photometric redshift distribution for a large unbiased sample of 870um selected submillimeter galaxies (SMGs) with robust identifications based on observations with the Atacama Large Millimeter Array (ALMA). In our analysis we consider 96 SMGs in the Extended Chandra Deep Field South, 77 of which have 4-19 band, optical-near-infrared, photometry. We model the Spectral Energy Distributions (SEDs) for these 77 SMGs, deriving a median photometric redshift of z=2.3+/-0.1. The remaining 19 SMGs have insufficient optical or near-infrared photometry to derive photometric redshifts, but a stacking analysis of IRAC and Herschel observations confirms they are not spurious. Assuming these sources have an absolute H-band magnitude distribution comparable to that of a complete sample of z~1-2 SMGs, we demonstrate that the undetected SMGs lie at higher redshifts, raising the median redshift for SMGs to z=2.5+/-0.2. More critically we show that the proportion of galaxies undergoing an SMG phase at z>3 is 35+/-5% of the total population. We derive a median stellar mass for SMGs of Mstar=(8+/-1)x10^10Mo, but caution that there are significant systematic uncertainties in our stellar mass estimate, up to x5 for individual sources. We compare our sample of SMGs to a volume-limited, morphologically classified sample of ellipticals in the local Universe. Assuming the star formation activity in SMGs has a timescale of ~100Myr we show that their descendants at z~0 would have a space density and M_H distribution which are in good agreement with those of local ellipticals. In addition the inferred mass-weighted ages of the local ellipticals broadly agree with the look-back times of the SMG events. Taken together, these results are consistent with a simple model that identifies SMGs as events that form most of the stars seen in the majority of luminous elliptical galaxies at the present day.

preprint2014arXiv

An ALMA Survey of Submillimetre Galaxies in the Extended Chandra Deep Field South: The Far-Infrared Properties of SMGs

We exploit ALMA 870um (345GHz) observations of submillimetre sources in the Extended Chandra Deep Field South to investigate the far-infrared properties of high-redshift submillimetre galaxies (SMGs). Using the precisely located 870um ALMA positions of 99 SMGs, together with 24um and radio imaging of this field, we deblend the Herschel/SPIRE imaging of this region to extract their far-infrared fluxes and colours. The median photometric redshifts for ALMA LESS (ALESS) SMGs which are detected in at least two SPIRE bands increases with wavelength of the peak in their SEDs, with z=2.3+/-0.2, 2.5+/-0.3 and 3.5+/-0.5 for the 250, 350 and 500-um peakers respectively. We find that 34 ALESS SMGs do not have a >3-sigma counterpart at 250, 350 or 500-um. These galaxies have a median photometric redshift of z=3.3+/-0.5, which is higher than the full ALESS SMG sample; z=2.5+/-0.2. Using the photometric redshifts together with the 250-870um photometry, we estimate the far-infrared luminosities and characteristic dust temperature of each SMG. The median infrared luminosity of the S_870um>2mJy SMGs is L_IR=(3.0+/-0.3)x10^{12}Lo(SFR=300+/-30Mo/yr). At a fixed luminosity, the characteristic dust temperature of these high-redshift SMGs is 2-3K lower than comparably luminous galaxies at z=0, reflecting the more extended star formation occurring in these systems. By extrapolating the 870um number counts to S_ 870um=1mJy, we show that the contribution of S_870um>1mJy SMGs to the cosmic star formation budget is 20% of the total over the redshift range z~1-4. We derive a median dust mass for these SMGs of M_d=(3.6+/-0.3)x10^8Mo and by adopting an appropriate gas-to-dust ratio, we estimate an average molecular mass of M_H2=(4.2+/-0.4)x10^{10}Mo. Finally, we use our estimates of the H2 masses to show that SMGs with S_870um>1mJy contain ~10% of the z~2 volume-averaged H2 mass density at this epoch.

preprint2014arXiv

Andromeda's Dust

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

preprint2014arXiv

Being WISE I: Validating Stellar Population Models and M/L ratios at 3.4 and 4.6 microns

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.

preprint2014arXiv

Confirming the Quiescent Galaxy Population out to $z=3$: A Stacking Analysis of Mid-, Far-Infrared and Radio Data

We present stringent constraints on the average mid-, far-infrared and radio emissions of $\sim$14200 quiescent galaxies (QGs), identified out to $z=3$ in the COSMOS field via their rest-frame NUV$-$r and r$-$J colors, and with stellar masses $M_{\star}=10^{9.8-12.2} \,M_{\odot} $. Stacking in deep Spitzer (MIPS $24\,μ$m), Herschel (PACS and SPIRE), and VLA (1.4 GHz) maps reveals extremely low dust-obscured star formation rates for QGs (SFR $<0.1-3\,M_{\odot}$yr$^{-1}$ at $z \leqslant 2$ and $<6-18\,M_{\odot}$yr$^{-1}$ at $z > 2$), consistent with the low unobscured SFRs ($<0.01-1.2\,M_{\odot}$yr$^{-1}$) inferred from modeling their ultraviolet-to-near-infrared photometry. The average SFRs of QGs are $>10\times$ below those of star-forming galaxies (SFGs) within the $M_{\star}$- and $z$-ranges considered. The stacked 1.4 GHz signals (S/N $> 5$) are, if attributed solely to star formation, in excess of the total (obscured plus unobscured) SFR limits, suggestive of a widespread presence of low-luminosity active galactic nuclei (AGN) among QGs. Our results reaffirm the existence of a significant population QGs out to $z = 3$, thus corroborating the need for powerful quenching mechanism(s) to terminate star formation in galaxies at earlier epochs.

preprint2014arXiv

Dust continuum emission as a tracer of gas mass in galaxies

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

preprint2014arXiv

High-Resolution Radio Continuum Measurements of the Nuclear Disks of Arp 220

We present new Karl G. Jansky Very Large Array radio continuum images of the nuclei of Arp 220, the nearest ultra-luminous infrared galaxy. These images have both the angular resolution to study detailed morphologies of the two nuclei that power the system and sensitivity to a wide range of spatial scales. At 33 GHz, and with a resolution of 0".081 x 0".063 (29.9 x 23.3 pc), we resolve the emission surrounding both nuclei and conclude that is mostly synchrotron in nature. The spatial distributions of radio emission in both nuclei are well described by exponential profiles. These have deconvolved half-light radii of 51 and 35 pc for the eastern and western nuclei, and they match the number density profile of radio supernovae observed with very long baseline interferometry. This similarity might be due to the fast cooling of cosmic rays electrons caused by the presence of a strong (~ mG) magnetic field in this system. We estimate high luminosity surface densities of $\mathrm{Σ_{IR} \sim 4.2^{+1.6}_{-0.7} \times 10^{13}}$ (east) and $\mathrm{\sim 9.7^{+3.7}_{-2.4} \times 10^{13}~(west)~L_{\odot}~kpc^{-2}}$, and star formation rate surface densities of $\mathrm{Σ_{SFR} \sim 10^{3.7\pm0.1}}$ (east) and $\mathrm{\sim 10^{4.1\pm0.1}~(west)~M_{\odot}~yr^{-1}~kpc^{-2}}$. These values, especially for the western nucleus are, to our knowledge, the highest luminosity and star formation rate surface densities measured for any star-forming system. Despite these high values, the nuclei lie below the dusty Eddington limit in which radiation pressure is balanced only by self-gravity. The small measured sizes also imply that the nuclei of Arp 220 are only transparent in the frequency range ~ 5 to 350 GHz. Our results offer no clear evidence that an active galactic nucleus dominates the emission from either nucleus at 33 GHz.

preprint2014arXiv

Inter-comparison of Radio-Loudness Criteria for Type 1 AGNs in the XMM-COSMOS Survey

Limited studies have been performed on the radio-loud fraction in X-ray selected type 1 AGN samples. The consistency between various radio-loudness definitions also needs to be checked. We measure the radio-loudness of the 407 type 1 AGNs in the XMM-COSMOS quasar sample using nine criteria from the literature (six defined in the rest-frame and three defined in the observed frame): $R_L=\log(L_{5GHz}/L_B)$, $q_{24}=\log(L_{24μm}/L_{1.4GHz})$, $R_{uv}=\log(L_{5GHz}/L_{2500Å})$, $R_{i}=\log(L_{1.4GHz}/L_i)$, $R_X=\log(νL_ν(5GHz)/L_X)$, $P_{5GHz}=\log(P_{5GHz}(W/Hz/Sr))$, $R_{L,obs}=\log(f_{1.4GHz}/f_B)$ (observed frame), $R_{i,obs}=\log(f_{1.4GHz}/f_i)$ (observed frame), and $q_{24, obs}=\log(f_{24μm}/f_{1.4GHz})$ (observed frame). Using any single criterion defined in the rest-frame, we find a low radio-loud fraction of $\lesssim 5\%$ in the XMM-COSMOS type 1 AGN sample, except for $R_{uv}$. Requiring that any two criteria agree reduces the radio-loud fraction to $\lesssim 2\%$ for about 3/4 of the cases. The low radio-loud fraction cannot be simply explained by the contribution of the host galaxy luminosity and reddening. The $P_{5GHz}=\log(P_{5GHz}(W/Hz/Sr))$ gives the smallest radio-loud fraction. Two of the three radio-loud fractions from the criteria defined in the observed frame without k-correction ($R_{L,obs}$ and $R_{i,obs}$) are much larger than the radio-loud fractions from other criteria.

preprint2014arXiv

Reconstructing the stellar mass distributions of galaxies using S4G IRAC 3.6 and 4.5 micron images: II. The conversion from light to mass

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.

preprint2014arXiv

The Astrophysics of Star Formation Across Cosmic Time at $\gtrsim$10 GHz with the Square Kilometre Array

In this chapter, we highlight a number of science investigations that are enabled by the inclusion of Band~5 ($4.6-13.8$ GHz) for SKA1-MID science operations, while focusing on the astrophysics of star formation over cosmic time. For studying the detailed astrophysics of star formation at high-redshift, surveys at frequencies $\gtrsim$10 GHz have the distinct advantage over traditional $\sim$1.4 GHz surveys as they are able to yield higher angular resolution imaging while probing higher rest frame frequencies of galaxies with increasing redshift, where emission of star-forming galaxies becomes dominated by thermal (free-free) radiation. In doing so, surveys carried out at $\gtrsim$10 GHz provide a robust, dust-unbiased measurement of the massive star formation rate by being highly sensitive to the number of ionizing photons that are produced. To access this powerful star formation rate diagnostic requires that Band~5 be available for SKA1-MID. We additionally present a detailed science case for frequency coverage extending up to 30 GHz during full SKA2 operations, as this allows for highly diverse science while additionally providing contiguous frequency coverage between the SKA and ALMA, which will likely be the two most powerful interferometers for the coming decades. To enable this synergy, it is crucial that the dish design of the SKA be flexible enough to include the possibility of being fit with receivers operating up to 30 GHz.

preprint2014arXiv

The PdBI Arcsecond Whirlpool Survey (PAWS): Environmental Dependence of Giant Molecular Cloud Properties in M51

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.

preprint2014arXiv

The PdBI Arcsecond Whirlpool Survey (PAWS): Multi-phase cold gas kinematic of M51

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.

preprint2014arXiv

The star-formation history of the Universe with the SKA

Radio wavelengths offer the unique possibility of tracing the total star-formation rate in galaxies, both obscured and unobscured. As such, they may provide the most robust measurement of the star-formation history of the Universe. In this chapter we highlight the constraints that the SKA can place on the evolution of the star-formation history of the Universe, the survey area required to overcome sample variance, the spatial resolution requirements, along with the multi-wavelength ancillary data that will play a major role in maximising the scientific promise of the SKA. The required combination of depth and resolution means that a survey to trace the star formation in the Universe should be carried out with a facility that has a resolution of at least ~0.5arcsec, with high sensitivity at < 1 GHz. We also suggest a strategy that will enable new parameter space to be explored as the SKA expands over the coming decade.

preprint2014arXiv

The VLA-COSMOS Survey: V. 324 MHz continuum observations

We present 90 cm VLA imaging of the COSMOS field, comprising a circular area of 3.14 square degrees at 8.0"x6.0" angular resolution with an average rms of 0.5 mJy/beam. The extracted catalog contains 182 sources (down to 5.5sigma), 30 of which are multi-component sources. Using Monte Carlo artificial source simulations we derive the completeness of the catalog, and we show that our 90 cm source counts agree very well with those from previous studies. Using X-ray, NUV-NIR and radio COSMOS data to investigate the population mix of our 90 cm radio sample, we find that our sample is dominated by active galactic nuclei (AGN). The average 90-20 cm spectral index (S_nu~nu**alpha, where S_nu is the flux density at frequency nu, and alpha the spectral index) of our 90 cm selected sources is -0.70, with an interquartile range of -0.90 to -0.53. Only a few ultra-steep-spectrum sources are present in our sample, consistent with results in the literature for similar fields. Our data do not show clear steepening of the spectral index with redshift. Nevertheless, our sample suggests that sources with spectral indices steeper than -1 all lie at z>1, in agreement with the idea that ultra-steep-spectrum radio sources may trace intermediate-redshift galaxies (z>1).

preprint2013arXiv

A Comparative Study of Giant Molecular Clouds in M51, M33 and the Large Magellanic Cloud

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.

preprint2013arXiv

Ammonia Thermometry of Star Forming Galaxies

With a goal toward deriving the physical conditions in external galaxies, we present a study of the ammonia (NH$_3$) emission and absorption in a sample of star forming systems. Using the unique sensitivities to kinetic temperature afforded by the excitation characteristics of several inversion transitions of NH$_3$, we have continued our characterization of the dense gas in star forming galaxies by measuring the kinetic temperature in a sample of 23 galaxies and one galaxy offset position selected for their high infrared luminosity. We derive kinetic temperatures toward 13 galaxies, 9 of which possess multiple kinetic temperature and/or velocity components. Eight of these galaxies exhibit kinetic temperatures $>100$ K, which are in many cases at least a factor of two larger than kinetic temperatures derived previously. Furthermore, the derived kinetic temperatures in our galaxy sample, which are in many cases at least a factor of two larger than derived dust temperatures, point to a problem with the common assumption that dust and gas kinetic temperatures are equivalent. As previously suggested, the use of dust emission at wavelengths greater than 160 $μ$m to derive dust temperatures, or dust heating from older stellar populations, may be skewing derived dust temperatures in these galaxies to lower values. We confirm the detection of high-excitation OH $^2Π_{3/2}$ J=9/2 absorption toward Arp220 (Ott et. al. 2011). We also report the first detections of non-metastable NH$_3$ inversion transitions toward external galaxies in the (2,1) (NGC253, NGC660, IC342, and IC860), (3,1), (3,2), (4,3), (5,4) (all in NGC660) and (10,9) (Arp220) transitions.

preprint2013arXiv

An ALMA survey of submillimetre galaxies in the Extended Chandra Deep Field South: High resolution 870um source counts

We report the first counts of faint submillimetre galaxies (SMG) in the 870-um band derived from arcsecond resolution observations with the Atacama Large Millimeter Array (ALMA). We have used ALMA to map a sample of 122 870-um-selected submillimetre sources drawn from the (0.5x0.5)deg^2 LABOCA Extended Chandra Deep Field South Submillimetre Survey (LESS). These ALMA maps have an average depth of sigma(870um)~0.4mJy, some ~3x deeper than the original LABOCA survey and critically the angular resolution is more than an order of magnitude higher, FWHM of ~1.5" compared to ~19" for the LABOCA discovery map. This combination of sensitivity and resolution allows us to precisely pin-point the SMGs contributing to the submillimetre sources from the LABOCA map, free from the effects of confusion. We show that our ALMA-derived SMG counts broadly agree with the submillimetre source counts from previous, lower-resolution single-dish surveys, demonstrating that the bulk of the submillimetre sources are not caused by blending of unresolved SMGs. The difficulty which well-constrained theoretical models have in reproducing the high-surface densities of SMGs, thus remains. However, our observations do show that all of the very brightest sources in the LESS sample, S(870um)>12mJy, comprise emission from multiple, fainter SMGs, each with 870-um fluxes of <9mJy. This implies a natural limit to the star-formation rate in SMGs of <10^3 M_Sun/yr, which in turn suggests that the space densities of z>1 galaxies with gas masses in excess of ~5x10^10 M_Sun is <10^-5 Mpc^-3. We also discuss the influence of this blending on the identification and characterisation of the SMG counterparts to these bright submillimetre sources and suggest that it may be responsible for previous claims that they lie at higher redshifts than fainter SMGs.

preprint2013arXiv

Clumping and the Interpretation of kpc-Scale Maps of the Interstellar Medium: Smooth HI and Clumpy, Variable H2 Surface Density

Many recent models consider the structure of individual interstellar medium (ISM) clouds as a way to explain observations of large parts of galaxies. To compare such models to observations, one must understand how to translate between surface densities observed averaging over large (~kpc) scales and surface densities on the scale of individual clouds (~pc scale), which are treated by models. We define a "clumping factor" that captures this translation as the ratio of the mass-weighted surface density, which is often the quantity of physical interest, to the area-weighted surface density, which is observed. We use high spatial resolution (sub-kpc) maps of CO and HI emission from nearby galaxies to measure the clumping factor of both atomic and molecular gas. The molecular and atomic ISM exhibit dramatically different degrees of clumping. As a result, the ratio H2/HI measured at ~kpc resolution cannot be trivially interpreted as a cloud-scale ratio of surface densities. HI emission appears very smooth, with a clumping factor of only ~1.3. Based on the scarce and heterogeneous high resolution data available, CO emission is far more clumped with a widely variable clumping factor, median ~7 for our heterogeneous data. Our measurements do not provide evidence for a universal mass-weighted surface density of molecular gas, but also cannot conclusively rule out such a scenario. We suggest that a more sophisticated treatment of molecular ISM structure, one informed by high spatial resolution CO maps, is needed to link cloud-scale models to kpc-scale observations of galaxies.

preprint2013arXiv

Gas Kinematics on GMC scales in M51 with PAWS: cloud stabilization through dynamical pressure

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

preprint2013arXiv

Molecular Gas and Star Formation in Nearby Disk Galaxies

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.

preprint2013arXiv

Probability Distribution Functions OF 12CO(J = 1-0) Brightness and Integrated Intensity in M51: The PAWS View

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)

preprint2013arXiv

Submillimeter ALMA Observations of the Dense Gas in the Low-Luminosity Type-1 Active Nucleus of NGC 1097

We present the first 100 pc scale view of the dense molecular gas in the central ~ 1.3 kpc region of the type-1 Seyfert NGC 1097 traced by HCN (J=4-3) and HCO+ (J=4-3) lines afforded with ALMA band 7. This galaxy shows significant HCN enhancement with respect to HCO+ and CO in the low-J transitions, which seems to be a common characteristic in AGN environments. Using the ALMA data, we study the characteristics of the dense gas around this AGN and search for the mechanism of HCN enhancement. We find a high HCN (J=4-3) to HCO+ (J=4-3) line ratio in the nucleus. The upper limit of the brightness temperature ratio of HCN (v2=1^{1f}, J=4-3) to HCN (J=4-3) is 0.08, which indicates that IR pumping does not significantly affect the pure rotational population in this nucleus. We also find a higher HCN (J=4-3) to CS (J=7-6) line ratio in NGC 1097 than in starburst galaxies, which is more than 12.7 on the brightness temperature scale. Combined from similar observations from other galaxies, we tentatively suggest that this ratio appears to be higher in AGN-host galaxies than in pure starburst ones similar to the widely used HCN to HCO+ ratio. LTE and non-LTE modeling of the observed HCN and HCO+ lines using J=4-3 and 1-0 data from ALMA, and J=3-2 data from SMA, reveals a high HCN to HCO+ abundance ratio (5 < [HCN]/[HCO+] < 20: non-LTE analysis) in the nucleus, and that the high-J lines (J=4-3 and 3-2) are emitted from dense (10^{4.5} < n_H2 [/cc] < 10^6), hot (70 < Tkin [K] < 550) regions. Finally we propose that the high temperature chemistry is more plausible to explain the observed enhanced HCN emission in NGC 1097 than the pure gas phase PDR/XDR chemistry.

preprint2013arXiv

The Paα Luminosity Function of HII Regions in Nearby Galaxies from HST/NICMOS

The HII region luminosity function (LF) is an important tool for deriving the birthrates and mass distribution of OB associations, and is an excellent tracer of the newly formed massive stars and associations. To date, extensive work (predominantly in Hα) has been done from the ground, which is hindered by dust extinction and the severe blending of adjacent (spatially or in projection) HII regions. Reliably measuring the properties of HII regions requires a linear resolution <40 pc, but analyses satisfying this requirement have been done only in a handful of galaxies, so far. As the first space-based work using a galaxy sample, we have selected 12 galaxies from our HST NICMOS Paα survey and studied the luminosity function and size distribution of HII regions both in individual galaxies and cumulatively, using a virtually extinction-free tracer of the ionizing photon rate. The high angular resolution and low sensitivity to diffuse emission of NICMOS also offer an advantage over ground-based imaging by enabling a higher degree of de-blending of the HII regions. We do not confirm the broken power-law LFs found in ground-based studies. Instead, we find that the LFs, both individual and co-added, follow a single power law dN(L)/dlnL ~ L^-1, consistent with the mass function of star clusters in nearby galaxies, and in agreement with the results of the existing analyses with HST data. The individual and co-added size distribution of HII regions are both roughly consistent with dN(D)/dlnD ~ D^-3, but the power-law scaling is probably contaminated by blended regions or complexes.

preprint2013arXiv

The PdBI Arcsecond Whirlpool Survey (PAWS). I. A Cloud-Scale/Multi-Wavelength View of the Interstellar Medium in a Grand-Design Spiral Galaxy

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.

preprint2012arXiv

An ALMA survey of Sub-millimetre Galaxies in the Extended Chandra Deep Field South: Detection of [C II] at z=4.4

We present ALMA 870-um (345GHz) observations of two sub-millimetre galaxies (SMGs) drawn from an ALMA study of the 126 sub-millimeter sources from the LABOCA Extended Chandra Deep Field South Survey (LESS). The ALMA data identify the counterparts to these previously unidentified sub-millimeter sources and serendipitously detect bright emission lines in their spectra which we show are most likely to be [C II]157.74um emission yielding redshifts of z=4.42 and z=4.44. This blind detection rate within the 7.5-GHz bandpass of ALMA is consistent with the previously derived photometric redshift distribution of SMGs and suggests a modest, but not dominant (<25%), tail of 870-um selected SMGs at z>4. We find that the ratio of L_CII/L_FIR in these SMGs is much higher than seen for similarly far-infrared-luminous galaxies at z~0, which is attributed to the more extended gas reservoirs in these high-redshift ULIRGs. Indeed, in one system we show that the [C II] emission shows hints of extended emission on >3kpc scales. Finally, we use the volume probed by our ALMA survey to show that the bright end of the [C II] luminosity function evolves strongly between z=0 and z~4.4, reflecting the increased ISM cooling in galaxies as a result of their higher star-formation rates. These observations demonstrate that even with short integrations, ALMA is able to detect the dominant fine structure cooling lines from high-redshift ULIRGs, measure their energetics and trace their evolution with redshift.

preprint2012arXiv

Disclosing the Radio Loudness Distribution Dichotomy in Quasars: An Unbiased Monte Carlo Approach Applied to the SDSS-FIRST Quasar Sample

We investigate the dichotomy in the radio loudness distribution of quasars by modelling their radio emission and various selection effects using a Monte Carlo approach. The existence of two physically distinct quasar populations, the radio-loud and radio-quiet quasars, is controversial and over the last decade a bimodal distribution of radio loudness of quasars has been both affirmed and disputed. We model the quasar radio luminosity distribution with simple unimodal and bimodal distribution functions. The resulting simulated samples are compared to a fiducial sample of 8,300 quasars drawn from the SDSS DR7 Quasar Catalog and combined with radio observations from the FIRST survey. Our results indicate that the SDSS-FIRST sample is best described by a radio loudness distribution which consists of two components, with 12+/-1 % of sources in the radio-loud component. On the other hand, the evidence for a local minimum in the loudness distribution (bimodality) is not strong and we find that previous claims for its existence were probably affected by the incompleteness of the FIRST survey close to its faint limit. We also investigate the redshift and luminosity dependence of the radio loudness distribution and find tentative evidence that at high redshift radio-loud quasars were rarer, on average "louder", and exhibited a smaller range in radio loudness. In agreement with other recent work, we conclude that the SDSS-FIRST sample strongly suggests that the radio loudness distribution of quasars is not a universal function, and that more complex models than presented here are needed to fully explain available observations.

preprint2012arXiv

Early type galaxies with tidal debris and their scaling relations in the Spitzer Survey of Stellar Structure in Galaxies (S4G)

Tidal debris around galaxies can yield important clues on their evolution. We have identified tidal debris in 11 early type galaxies (T \leq 0) from a sample of 65 early types drawn from the Spitzer Survey of Stellar Structure in Galaxies (S4G). The tidal debris includes features such as shells, ripples and tidal tails. A variety of techniques, including two-dimensional decomposition of galactic structures, was used to quantify the residual tidal features. The tidal debris contributes ~3 - 10% to the total 3.6 μm luminosity of the host galaxy. Structural parameters of the galaxies were estimated using two-dimensional profile fitting. We investigate the locations of galaxies with tidal debris in the Fundamental Plane and Kormendy relation. We find that galaxies with tidal debris lie within the scatter of early type galaxies without tidal features. Assuming that the tidal debris is indicative of recent gravitational interaction or merger, this suggests that these galaxies have either undergone minor merging events so that the overall structural properties of the galaxies are not significantly altered, or they have undergone a major merging events but already have experienced sufficient relaxation and phase-mixing so that their structural properties become similar to those of the non-interacting early type galaxies.

preprint2012arXiv

Extended X-ray emission from non-thermal sources in the COSMOS field: A detailed study of a large radio galaxy at z=1.168

X-ray selected galaxy group samples are usually generated by searching for extended X- ray sources that reflect the thermal radiation of the intragroup medium. On the other hand, large radio galaxies that regularly occupy galaxy groups also emit in the X-ray window, and their contribution to X-ray selected group samples is still not well understood. In order to investigate their relative importance, we have carried out a systematic search for non-thermal extended X-ray sources in the COSMOS field. Based on the morphological coincidence of X-ray and radio extensions, out of 60 radio galaxies, and \sim 300 extended X-ray sources, we find only one candidate where the observed extended X-ray emission arises from non- thermal processes related to radio galaxies. We present a detailed analysis of this source, and its environment. Our results yield that external Inverse Compton emission of the lobes is the dominant process that generates the observed X-ray emission of our extended X-ray candidate, with a minor contribution from the gas of the galaxy group hosting the radio galaxy. Finally, we show that finding only one potential candidate in the COSMOS field (in a redshift range 0 < z < 6 and with radio luminosity between 1025 and 1030 W/Hz) is consistent with expected X-ray-counts arising from synchrotron lobes. This implies that these sources are not a prominent source of contamination in samples of X-ray selected clusters/groups, but they could potentially dominate the z > 1 cluster counts at the bright end (S_X > 7 \cdot 10^-15 erg s^-1 cm^2).

preprint2012arXiv

Gas and dust in a z=2.8 obscured quasar

We present new detections of the CO(5-4), CO(7-6), [CI](1-0) and [CI](2-1) molecular and atomic line transitions towards the unlensed, obscured quasar AMS12 (z=2.7672), observed with the IRAM PdBI. This is the first unlensed, high redshift source to have both [CI] transitions detected. Continuum measurements between 70 $μ$m and 3 mm are used to constrain the FIR SED, and we find a best fit FIR luminosity of log[Lfir/Lsol] = 13.5+/-0.1, dust temperature T_d = 88+/-8 K and emissivity index β = 0.6+/-0.1. The highly-excited molecular gas probed by CO(3-2), (5-4) and (7-6), is modelled with large velocity gradient (LVG) models. The gas kinetic temperature T_g, density n(H2), and the characteristic size r0, are determined using the dust temperature from the FIR SED as a prior for the gas temperature. The best fitting parameters are T_g = 90+/-8 K, n(H2) = 10^(3.9+/-0.1) cm^(-3) and r0 = 0.8+/-0.04 kpc. The ratio of the [CI] lines gives a [CI] excitation temperature of 43+/-10 K, indicating the [CI] and the high-excitation CO are not in thermal equilibrium. The [CI] excitation temperature is below that of T_d and T_g of the high-excitation CO, perhaps because [CI] lies at a larger radius where there may also be a large reservoir of CO at a cooler temperature, perhaps detectable through the CO(1-0). Using the [CI](1-0) line we can estimate the strength of the CO(1-0) line and hence the gas mass. This suggests that a significant fraction (~30%) of the molecular gas is missed from the high-excitation line analysis. The Eddington limited black hole mass is found from the bolometric luminosity to be Mbh >~ 1.5x10^9 Msol. Along with the stellar mass of 3x10^11 Msol, these give a black hole - bulge mass ratio of Mbh/Mbulge >~ 0.005. This is in agreement with studies on the evolution of the Mbh/Mbulge relationship at high redshifts, which find a departure from the local value ~0.002.

preprint2012arXiv

Herschel-PACS observations of [OI]63um towards submillimetre galaxies at z~1

We present Herschel-PACS spectroscopy of the [OI]63um far-infrared cooling line from a sample of six unlensed and spectroscopically-confirmed 870um-selected submillimetre (submm) galaxies (SMGs) at 1.1<z<1.6 from the LABOCA Extended Chandra Deep Field South (ECDFS) Submm Survey (LESS). This is the first survey of [OI]63um, one of the main photodissociation region (PDR) cooling lines, in SMGs. New high-resolution ALMA interferometric 870um continuum imaging confirms that these six Herschel-targeted SMG counterparts are bona fide sources of submm emission. We detect [OI]63um in two SMGs with a SNR >3, tentatively detect [OI]63um in one SMG, and constrain the line flux for the non-detections. We also exploit the combination of submm continuum photometry from 250-870um and our new PACS continuum measurements to constrain the far-infrared (FIR) luminosity, L_FIR, in these SMGs to < 30%. We find that SMGs do not show a deficit in their [OI]63um-to-far-infrared continuum luminosity ratios (with ratios ranging from ~0.5-1.5%), similar to what was seen previously for the [CII]158um-to-FIR ratios in SMGs. These observed ratios are about an order of magnitude higher than what is seen typically for local ultra luminous infrared galaxies (ULIRGs), which adds to the growing body of evidence that SMGs are not simply `scaled up' versions of local ULIRGs. Rather, the PDR line-to-L_FIR ratios suggest that the star formation modes of SMGs are likely more akin to that of local normal (lower-luminosity) star-forming galaxies, with the bulk of the star formation occurring in extended regions, galaxy-scale (~kpc) in size. These observations represent the first step towards a census of the major PDR cooling lines in typical SMGs that will be attainable with ALMA, enabling detailed modelling to probe the global properties of the star formation and the evolutionary status of SMGs.

preprint2012arXiv

Properties of Bulgeless Disk Galaxies II. Star Formation as a Function of Circular Velocity

We study the relation between the surface density of gas and star formation rate in twenty moderately-inclined, bulgeless disk galaxies (Sd-Sdm Hubble types) using CO(1-0) data from the IRAM 30m telescope, HI emission line data from the VLA/EVLA, H-alpha data from the MDM Observatory, and PAH emission data derived from Spitzer IRAC observations. We specifically investigate the efficiency of star formation as a function of circular velocity (v_circ). Previous work found that the vertical dust structure and disk stability of edge-on, bulgeless disk galaxies transition from diffuse dust lanes with large scale heights and gravitationally-stable disks at v_circ < 120 km/s (M_star <~ 10^10 M_sun) to narrow dust lanes with small scale heights and gravitationally-unstable disks at v_circ > 120 km/s. We find no transition in star formation efficiency (Sigma_SFR/Sigma_HI+H2) at v_circ = 120 km/s, or at any other circular velocity probed by our sample (v_circ = 46 - 190 km/s). Contrary to previous work, we find no transition in disk stability at any circular velocity in our sample. Assuming our sample has the same dust structure transition as the edge-on sample, our results demonstrate that scale height differences in the cold interstellar medium of bulgeless disk galaxies do not significantly affect the molecular fraction or star formation efficiency. This may indicate that star formation is primarily affected by physical processes that act on smaller scales than the dust scale height, which lends support to local star formation models.

preprint2012arXiv

The heating of dust by old stellar populations in the Bulge of M31

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

preprint2012arXiv

The Resolved Radio--FIR Correlation in Nearby Galaxies with Herschel and Spitzer

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

preprint2012arXiv

The S4G perspective on circumstellar dust extinction of AGB stars in M100

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.

preprint2011arXiv

A constant characteristic mass for star forming galaxies since z~3 revealed by radio emission in the COSMOS field

We present results of our 1.4 GHz image stacking analysis of mass-selected galaxies in the COSMOS field. From the resulting median radio continuum flux density we have determined the evolution of the average star formation rate (SFR) of galaxies as a function of stellar mass, unbiased from effects of dust but also source confusion due to the 1.5" angular resolution achieved by the VLA. We find a power-law relation between specific SFR (SSFR) and stellar mass for star forming galaxies out to z=3. While higher mass systems exhibit lower SSFRs at any epoch, no differential, more rapid evolution of high mass galaxies is evident at least out to z~1.5 where our conclusions are most robust. Utilizing measured mass functions of star forming systems, the characteristic stellar mass for galaxies contributing most to the comoving SFR density appears not to evolve. These findings hence challenge 'downsizing' scenarios in which star formation has gradually shifted towards lower mass systems with cosmic time. Our analysis constitutes the to-date best determination of the cosmic star formation history (CSFH) since z=3 and yields indirect evidence for a rapid decline of the global mass density of molecular gas with time.

preprint2011arXiv

A Deep Chandra ACIS Study of NGC 4151. I. the X-ray Morphology of the 3 kpc-diameter Circum-nuclear Region and Relation to the Cold Interstellar Medium

We report on the imaging analysis of 200 ks sub-arcsecond resolution Chandra ACIS-S observations of the nearby Seyfert 1 galaxy NGC 4151. Bright, structured soft X-ray emission is observed to extend from 30 pc to 1.3 kpc in the south-west from the nucleus, much farther than seen in earlier X-ray studies. The terminus of the north-eastern X-ray emission is spatially coincident with a CO gas lane, where the outflow likely encounters dense gas in the host galactic disk. X-ray emission is also detected outside the boundaries of the ionization cone, which indicates that the gas there is not completely shielded from the nuclear continuum, as would be the case for a molecular torus collimating the bicone. In the central r<200 pc region, the subpixel processing of the ACIS data recovers the morphological details on scales of <30~pc (<0.5") first discovered in Chandra HRC images. The X-ray emission is more absorbed towards the boundaries of the ionization cone, as well as perpendicular to the bicone along the direction of a putative torus in NGC 4151. The innermost region where X-ray emission shows the highest hardness ratio, is spatially coincident with the near-infrared resolved H_2 emission and dusty spirals we find in an HST V-H color image. The agreement between the observed H_2 line flux and the value predicted from X-ray-irradiated molecular cloud models supports photo-excitation by X-rays from the active nucleus as the origin of the H_2 line, although contribution from UV fluorescence or collisional excitation cannot be fully ruled out with current data. The discrepancy between the mass of cold molecular gas inferred from recent CO and near-infrared H_2 observations may be explained by the anomalous CO abundance in this X-ray dominated region. The total H_2 mass derived from the X-ray observation agrees with measurement in Storchi-Bergmann et al.

preprint2011arXiv

A Deep Chandra ACIS Study of NGC 4151. III. the Line Emission and Spectral Analysis of the Ionization Cone

This paper is the third in a series in which we present deep Chandra ACIS-S imaging spectroscopy of the Seyfert 1 galaxy NGC 4151, devoted to study its complex circum-nuclear X-ray emission. Emission features in the soft X-ray spectrum of the bright extended emission (L[0.3-2keV]~10^40 erg/s) at r>130 pc (2") are consistent with the brighter OVII, OVIII, and NeIX lines seen in the Chandra HETGS and XMM-Newton RGS spectra below 2 keV. We construct emission line images of these features and find good morphological correlations with the narrow line region clouds mapped in [OIII]5007A. Self-consistent photoionization models provide good descriptions of the spectra of the large scale emission, as well as resolved structures, supporting the dominant role of nuclear photoionization, although displacement of optical and X-ray features implies a more complex medium. Collisionally ionized emission is estimated to be <12% of the extended emission. Presence of both low and high ionization spectral components and extended emission in the X-ray image perpendicular to the bicone indicates leakage of nuclear ionization, likely filtered through warm absorbers, instead of being blocked by a continuous obscuring torus. The ratios of [OIII]/soft X-ray flux are approximately constant (~15) for the 1.5 kpc radius spanned by these measurements, indicating a relatively constant ionization parameter, consistent with the photoionized outflow of a wind-like density profile. Using spatially resolved features, we estimate that the mass outflow rate in NGC 4151 is ~2Msun/yr at 130 pc and the kinematic power of the ionized outflow is 1.7x10^41 erg/s, approximately 0.3% of the bolometric luminosity of NGC 4151.

preprint2011arXiv

A massive proto-cluster of galaxies at a redshift of z {\approx} 5.3

Massive clusters of galaxies have been found as early as 3.9 Billion years (z=1.62) after the Big Bang containing stars that formed at even earlier epochs. Cosmological simulations using the current cold dark matter paradigm predict these systems should descend from "proto-clusters" - early over-densities of massive galaxies that merge hierarchically to form a cluster. These proto-cluster regions themselves are built-up hierarchically and so are expected to contain extremely massive galaxies which can be observed as luminous quasars and starbursts. However, observational evidence for this scenario is sparse due to the fact that high-redshift proto-clusters are rare and difficult to observe. Here we report a proto-cluster region 1 billion years (z=5.3) after the Big Bang. This cluster of massive galaxies extends over >13 Mega-parsecs, contains a luminous quasar as well as a system rich in molecular gas. These massive galaxies place a lower limit of >4x10^11 solar masses of dark and luminous matter in this region consistent with that expected from cosmological simulations for the earliest galaxy clusters.

preprint2011arXiv

Complex Radio Spectral Energy Distributions in Luminous and Ultraluminous Infrared Galaxies

We use the Expanded Very Large Array to image radio continuum emission from local luminous and ultraluminous infrared galaxies (LIRGs and ULIRGs) in 1 GHz windows centered at 4.7, 7.2, 29, and 36 GHz. This allows us to probe the integrated radio spectral energy distribution (SED) of the most energetic galaxies in the local universe. The 4-8 GHz flux densities agree well with previous measurements. They yield spectral indices α\approx -0.67 (where F_ν\propto ν^α) with \pm 0.15 (1σ) scatter, typical of nonthermal (synchrotron) emission from star-forming galaxies. The contrast of our 4-8 GHz data with literature 1.5 and 8.4 GHz flux densities gives further evidence for curvature of the radio SED of U/LIRGs. The SED appears flatter near \sim 1 GHz than near \sim 6 GHz, suggesting significant optical depth effects at the lower frequencies. The high frequency (28-37 GHz) flux densities are low compared to extrapolations from the 4-8 GHz data. We confirm and extend to higher frequency a previously observed deficit of high frequency radio emission for luminous starburst galaxies.

preprint2011arXiv

Cyanoacetylene in IC 342: An Evolving Dense Gas Component with Starburst Age

We present the first images of the J=5-4 and J=16-15 lines of the dense gas tracer, cyanoacetylene, HC_3N, in an external galaxy. The central 200 pc of the nearby star-forming spiral galaxy, IC 342, was mapped using the VLA and the Plateau de Bure Interferometer. HC_3N(5-4) line emission is found across the nuclear mini-spiral, but is very weak towards the starburst site, the location of the strongest mid-IR and radio emission. The J=16-15 and 10-9 lines are also faint near the large HII region complex, but are brighter relative to the 5-4 line, consistent with higher excitation. The brightest HC_3N emission is located in the northern arm of the nuclear minispiral, 100 pc away from the radio/IR source to the southwest of the nucleus. This location appears less affected by ultraviolet radiation, and may represent a more embedded, earlier stage of star formation. HC_3N excitation temperatures are consistent with those determined from C^{18}O; the gas is dense, 10^{4-5}/cc, and cool, T_K ~< 40 K. So as to not violate limits on the total H_2 mass determined from C^{18}O, at least two dense components are required to model IC 342's giant molecular clouds. These observations suggest that HC_3N(5-4) is an excellent probe of the dense, quiescent gas in galaxies. The high excitation combined with faint emission towards the dense molecular gas at the starburst indicates that it currently lacks large masses of very dense gas. We propose a scenario where the starburst is being caught in the act of dispersing or destroying its dense gas in the presence of the large HII region. This explains the high star formation efficiency seen in the dense component. The little remaining dense gas appears to be in pressure equilibrium with the starburst HII region.

preprint2011arXiv

Molecular gas around low-luminosity AGN in late-type spirals

We have studied the molecular gas in the vicinity of low-luminosity active galactic nuclei (AGNs) in three bulge-less spiral galaxies: NGC 1042, NGC 4178, and NGC 4395. The (1-0) and (2-1) transitions of gaseous carbon monoxide (CO) are clearly detected within the central kpc of all three galaxies. In the case of NGC 4395, this constitutes the first reported detection of CO. In general, the CO emission is faint, as may be expected from their less-than-spectacular star formation activity. Interestingly, however, both face-on galaxies in our sample (which allow an unimpeded view of their nucleus) show an elevated intensity ratio CO(2-1)/CO(1-0) when compared to similar late-type spirals without an AGN. We discuss that this is unlikely due to a very compact CO source. Instead, we speculate that even energetically weak AGN can impact the physical state of the surrounding gas. We do not detect any tracers of dense molecular gas such as HCN or HCO+, but the sensitivity of our observations allows us to establish upper limits that lie at the low end of the range observed in more energetic AGN. The derived gas density is less than n(H2)\approx 2e03 cm^-3 which is significantly lower than in most other nearby galaxies. The scarcity of dense gas suggests that the conditions for star formation are poor in these nuclei.

preprint2011arXiv

New Constraints on Mass-Dependent Disruption of Star Clusters in M51

We use UBVI,Ha images of the Whirlpool galaxy, M51, taken with the ACS and WFPC2 cameras on the Hubble Space Telescope (HST) to select star clusters, and to estimate their masses and ages by comparing their observed colors with predictions from population synthesis models. We construct the mass function of intermediate age (1-4x10^8 yr) clusters, and find that it is well described by a power law, psi(M) propto M^beta, with beta=-2.1 +/- 0.2, for clusters more massive than approximately 6x10^3 Msun. This extends the mass function of intermediate age clusters in M51 to masses lower by nearly a factor of five over previous determinations. The mass function does not show evidence for curvature at either the high or low mass end. This shape indicates that there is no evidence for the earlier disruption of lower mass clusters compared with their higher mass counterparts (i.e., no mass-dependent disruption) over the observed range of masses and ages, or for a physical upper mass limit Mc with which clusters in M51 can form. These conclusions differ from previous suggestions based on poorer-quality HST observations. We discuss their implications for the formation and disruption of the clusters. Ages of clusters in two "feathers," stellar features extending from the outer portion of a spiral arm, show that the feather with a larger pitch angle formed earlier, and over a longer period, than the other.

preprint2011arXiv

Properties of Bulgeless Disk Galaxies. I. Atomic Gas

We study the neutral hydrogen properties of a sample of 20 bulgeless disk galaxies (Sd - Sdm Hubble types), an interesting class that can be used to constrain galaxy formation and evolution, especially the role of mergers versus internal processes. Our sample is composed of nearby (within 32 Mpc), moderately inclined galaxies that bracket the circular velocity of 120 km/s, which has been found to be associated with a transition in dust scale heights in edge-on, late-type disks. Here we present HI channel maps, line profiles, and integrated intensity maps. We also derive kinematic parameters, including the circular velocity, from rotation curve analyses and calculate the integrated HI flux and HI mass for each galaxy in the sample. Three of the 20 galaxies in our sample have kinematically distinct outer components with major axes that differ by 30 - 90 degrees from the main disk. These distinct outer components may be due to a recent interaction, which would be somewhat surprising because the disks do not contain bulges. We will use the data products and derived properties in subsequent investigations into star formation and secular evolution in bulgeless disks with circular velocities above and below 120 km/s.

preprint2011arXiv

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

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.

preprint2011arXiv

Star and dust formation activities in AzTEC-3: A starburst galaxy at z = 5.3

Analyses of high-redshift ultraluminous infrared (IR) galaxies traditionally use the observed optical to submillimeter spectral energy distribution (SED) and estimates of the dynamical mass as observational constraints to derive the star formation rate (SFR), the stellar mass, and age of these objects. An important observational constraint neglected in the analysis is the mass of dust giving rise to the IR emission. In this paper we add this constraint to the analysis of AzTEC-3. Adopting an upper limit to the mass of stars and a bolometric luminosity for this object, we construct stellar and chemical evolutionary scenarios, constrained to produce the inferred dust mass and observed luminosity before the associated stellar mass exceeds the observational limit. We find that the model with a Top Heavy IMF provided the most plausible scenario consistent with the observational constraints. In this scenario the dust formed over a period of ~200 Myr, with a SFR of ~500 Msun/yr. These values for the age and SFR in AzTEC-3 are significantly higher and lower, respectively, from those derived without the dust mass constraint. However, this scenario is not unique, and others cannot be completely ruled out because of the prevailing uncertainties in the age of the galaxy, its bolometric luminosity, and its stellar and dust masses. A robust result of our models is that all scenarios require most of the radiating dust mass to have been accreted in molecular clouds. Our new procedure highlights the importance of a multiwavelength approach, and of the use of dust evolution models in constraining the age and the star formation activity and history in galaxies.

preprint2011arXiv

The Emission by Dust and Stars of Nearby Galaxies in the Herschel KINGFISH Survey

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.

preprint2011arXiv

The star formation history of mass-selected galaxies in the COSMOS field

We explore the evolution of the specific star formation rate (SSFR) for 3.6um-selected galaxies of different M_* in the COSMOS field. The average SFR for sub-sets of these galaxies is estimated with stacked 1.4GHz radio continuum emission. We separately consider the total sample and a subset of galaxies (SF) that shows evidence for substantive recent star formation in the rest-frame optical SED. At 0.2<z<3 both populations show a strong and M_*-independent decrease in their SSFR towards z=0.2, best described by a power- law (1+z)^n, where n~4.3 for all galaxies and n~3.5 for SF sources. The decrease appears to have started at z>2, at least above 4x10^10M_Sun where our conclusions are most robust. We find a tight correlation with power-law dependence, SSFR (M_*)^beta, between SSFR and M_* at all z. It tends to flatten below ~10^10M_Sun if quiescent galaxies are included; if they are excluded a shallow index beta_SFG -0.4 fits the correlation. On average, higher M_* objects always have lower SSFRs, also among SF galaxies. At z>1.5 there is tentative evidence for an upper SSFR-limit that an average galaxy cannot exceed. It is suggested by a flattening of the SSFR-M_* relation (also for SF sources), but affects massive (>10^10M_Sun) galaxies only at the highest z. Below z=1.5 there thus is no direct evidence that galaxies of higher M_* experience a more rapid waning of their SSFR than lower M_* SF systems. In this sense, the data rule out any strong 'downsizing'. We combine our results with recent measurements of the galaxy (stellar) mass function in order to determine the characteristic mass of a SF galaxy (M_*=10^(10.6\pm0.4)M_Sun). In this sense, too, there is no 'downsizing'. Our analysis constitutes the most extensive SFR density determination with a single technique to z=3. Recent Herschel results are consistent with our results, but rely on far smaller samples.

preprint2011arXiv

The unusual vertical mass distribution of NGC 4013 seen through the Spitzer Survey of Stellar Structure in Galaxies (S4G)

NGC 4013 is a nearby Sb edge-on galaxy known for its "prodigious" HI warp and its "giant" tidal stream. Previous work on this unusual object shows that it cannot be fitted satisfactorily by a canonical thin+thick disk structure. We have produced a new decomposition of NGC 4013, considering three stellar flattened components (thin+thick disk plus an extra and more extended component) and one gaseous disk. All four components are considered to be gravitationally coupled and isothermal. To do so, we have used the 3.6micron images from the Spitzer Survey of Stellar Structure in Galaxies (S4G). We find evidence for NGC4013 indeed having a thin and a thick disk and an extra flattened component. This smooth and extended component (scaleheight z_EC kpc) could be interpreted as a thick disk or as a squashed ellipsoidal halo and contains ~20% of the total mass of all three stellar components. We argue it is unlikely to be related to the ongoing merger or due to the off-plane stars from a warp in the other two disk components. Instead, we favor a scenario in which the thick disk and the extended component were formed in a two-stage process, in which an initially thick disk has been dynamically heated by a merger soon enough in the galaxy history to have a new thick disk formed within it.

preprint2011arXiv

Thick disks of edge-on galaxies seen through the Spitzer Survey of Stellar Structure in Galaxies (S4G): Lair of missing baryons?

Most, if not all, disk galaxies have a thin (classical) disk and a thick disk. In most models thick disks are thought to be a necessary consequence of the disk formation and/or evolution of the galaxy. We present the results of a study of the thick disk properties in a sample of carefully selected edge-on galaxies with types ranging from T=3 to T=8. We fitted one-dimensional luminosity profiles with physically motivated functions - the solutions of two stellar and one gaseous isothermal coupled disks in equilibrium - which are likely to yield more accurate results than other functions used in previous studies. The images used for the fits come from the Spitzer Survey of Stellar Structure in Galaxies (S4G). We found that thick disks are on average more massive than previously reported, mostly due to the selected fitting function. Typically, the thin and the thick disk have similar masses. We also found that thick disks do not flare significantly within the observed range in galactocentric radii and that the ratio of thick to thin disk scaleheights is higher for galaxies of earlier types. Our results tend to favor an in situ origin for most of the stars in the thick disk. In addition the thick disk may contain a significant amount of stars coming from satellites accreted after the initial build-up of the galaxy and an extra fraction of stars coming from the secular heating of the thin disk by its own overdensities. Assigning thick disk light to the thin disk component may lead to an underestimate of the overall stellar mass in galaxies, because of different mass to light ratios in the two disk components. On the basis of our new results, we estimate that disk stellar masses are between 10% and 50% higher than previously thought and we suggest that thick disks are a reservoir of "local missing baryons".

preprint2010arXiv

A Massive Molecular Gas Reservoir in the z=5.3 Submillimeter Galaxy AzTEC-3

We report the detection of CO 2-1, 5-4, and 6-5 emission in the highest-redshift submillimeter galaxy (SMG) AzTEC-3 at z=5.298, using the Expanded Very Large Array and the Plateau de Bure Interferometer. These observations ultimately confirm the redshift, making AzTEC-3 the most submillimeter-luminous galaxy in a massive z=5.3 protocluster structure in the COSMOS field. The strength of the CO line emission reveals a large molecular gas reservoir with a mass of 5.3e10 (alpha_CO/0.8) Msun, which can maintain the intense 1800 Msun/yr starburst in this system for at least 30 Myr, increasing the stellar mass by up to a factor of six in the process. This gas mass is comparable to `typical' z~2 SMGs, and constitutes >~80% of the baryonic mass (gas+stars) and 30%-80% of the total (dynamical) mass in this galaxy. The molecular gas reservoir has a radius of <4 kpc and likely consists of a `diffuse', low-excitation component, containing (at least) 1/3 of the gas mass (depending on the relative conversion factor alpha_CO), and a `dense', high-excitation component, containing ~2/3 of the mass. The likely presence of a substantial diffuse component besides highly-excited gas suggests different properties between the star-forming environments in z>4 SMGs and z>4 quasar host galaxies, which perhaps trace different evolutionary stages. The discovery of a massive, metal-enriched gas reservoir in a SMG at the heart of a large z=5.3 protocluster considerably enhances our understanding of early massive galaxy formation, pushing back to a cosmic epoch where the Universe was less than 1/12 of its present age.

preprint2010arXiv

Detection of molecular gas in a distant submillimetre galaxy at z=4.76 with ATCA

We have detected the CO(2-1) transition from the submillimetre galaxy (SMG) LESSJ033229.4-275619 at z=4.755 using the new Compact Array Broadband Backend system on the Australian Telescope Compact Array. These data have identified a massive gas reservoir available for star formation for the first time in an SMG at z~5. We use the luminosity and velocity width (FWHM of 160 km/s) of the CO(2--1) line emission to constrain the gas and dynamical mass of Mgas~1.6x10^10 Msun and Mdyn(<2kpc)~5x10^10 (0.25/sin^2(i)) Msun, respectively, similar to that observed for SMGs at lower redshifts of z~2-4, although we note that our observed CO FWHM is a factor of ~3 narrower than typically seen in SMGs. Together with the stellar mass we estimate a total baryonic mass of Mbary~1x10^11 Msun, consistent with the dynamical mass for this young galaxy within the uncertainties. Dynamical and baryonic mass limits of high-redshift galaxies are useful tests of galaxy formation models: using the known z~4-5 SMGs as examples of massive baryonic systems, we find that their space density is consistent with that predicted by current galaxy formation models. In addition, these observations have helped to confirm that z~4-5 SMGs possess the baryonic masses and gas consumption timescales necessary to be the progenitors of the luminous old red galaxies seen at z~3. Our results provide a preview of the science that ALMA will enable on the formation and evolution of the earliest massive galaxies in the Universe.

preprint2010arXiv

Extended X-ray Emission in the HI Cavity of NGC 4151: Galaxy-scale AGN Feedback?

We present the Chandra discovery of soft diffuse X-ray emission in NGC 4151 (L[0.5-2keV]~10^{39} erg s$^{-1}$), extending ~2 kpc from the active nucleus and filling in the cavity of the HI material. The best fit to the X-ray spectrum requires either a kT~0.25 keV thermal plasma or a photoionized component. In the thermal scenario, hot gas heated by the nuclear outflow would be confined by the thermal pressure of the HI gas and the dynamic pressure of inflowing neutral material in the galactic disk. In the case of photoionization, the nucleus must have experienced an Eddington limit outburst. For both scenarios, the AGN-host interaction in NGC 4151 must have occured relatively recently (some 10^4 yr ago). This very short timescale to the last episode of high activity phase may imply such outbursts occupy $\gtrsim$1% of AGN lifetime.

preprint2010arXiv

The Spitzer Survey of Stellar Structure in Galaxies (S^4G)

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.

preprint2010arXiv

The thick disk in the galaxy NGC 4244 from S4G imaging

If thick disks are ubiquitous and a natural product of disk galaxy formation and/or evolution processes, all undisturbed galaxies which have evolved during a significant fraction of a Hubble time should have a thick disk. The late-type spiral galaxy NGC 4244 has been reported as the only nearby edge-on galaxy without a confirmed thick disk. Using data from the Spitzer Survey of Stellar Structure in Galaxies (S4G) we have identified signs of two disk components in this galaxy. The asymmetries between the light profiles on both sides of the mid-plane of NGC 4244 can be explained by a combination of the galaxy not being perfectly edge-on and a certain degree of opacity of the thin disk. We argue that the subtlety of the thick disk is a consequence of either a limited secular evolution in NGC 4244, a small fraction of stellar material in the fragments which built the galaxy, or a high amount of gaseous accretion after the formation of the galaxy.

preprint2010arXiv

The VLA-COSMOS Survey. IV. Deep Data and Joint Catalog

In the context of the VLA-COSMOS Deep project additional VLA A array observations at 1.4 GHz were obtained for the central degree of the COSMOS field and combined with the existing data from the VLA-COSMOS Large project. A newly constructed Deep mosaic with a resolution of 2.5" was used to search for sources down to 4 sigma with 1 sigma ~ 12 microJy/beam in the central 50'x50'. This new catalog is combined with the catalog from the Large project (obtained at 1.5"x1.4" resolution) to construct a new Joint catalog. All sources listed in the new Joint catalog have peak flux densities of >5 sigma at 1.5" and/or 2.5" resolution to account for the fact that a significant fraction of sources at these low flux levels are expected to be slighty resolved at 1.5" resolution. All properties listed in the Joint catalog such as peak flux density, integrated flux density and source size are determined in the 2.5" resolution Deep image. In addition, the Joint catalog contains 43 newly identified multi-component sources.

preprint2009arXiv

BLAST: the far-infrared/radio correlation in distant galaxies

We investigate the correlation between far-infrared (FIR) and radio luminosities in distant galaxies, a lynchpin of modern astronomy. We use data from the Balloon-borne Large Aperture Submillimetre Telescope (BLAST), Spitzer, the Large Apex BOlometer CamerA (LABOCA), the Very Large Array (VLA) and the Giant Metre-wave Radio Telescope (GMRT) in the Extended Chandra Deep Field South (ECDFS). For a catalogue of BLAST 250-micron-selected galaxies, we re-measure the 70--870-micron flux densities at the positions of their most likely 24-micron counterparts, which have a median [interquartile] redshift of 0.74 [0.25, 1.57]. From these, we determine the monochromatic flux density ratio, q_250 = log_10 (S_250micron / S_1400MHz), and the bolometric equivalent, q_IR. At z ~= 0.6, where our 250-micron filter probes rest-frame 160-micron emission, we find no evolution relative to q_160 for local galaxies. We also stack the FIR and submm images at the positions of 24-micron- and radio-selected galaxies. The difference between q_IR seen for 250-micron- and radio-selected galaxies suggests star formation provides most of the IR luminosity in ~< 100-uJy radio galaxies, but rather less for those in the mJy regime. For the 24-micron sample, the radio spectral index is constant across 0 < z < 3, but q_IR exhibits tentative evidence of a steady decline such that q_IR is proportional to (1+z)^(-0.15 +/- 0.03) - significant evolution, spanning the epoch of galaxy formation, with major implications for techniques that rely on the FIR/radio correlation. We compare with model predictions and speculate that we may be seeing the increase in radio activity that gives rise to the radio background.

preprint2009arXiv

Environment of MAMBO galaxies in the COSMOS field

Submillimeter galaxies (SMG) represent a dust-obscured high-redshift population undergoing massive star formation activity. Their properties and space density have suggested that they may evolve into spheroidal galaxies residing in galaxy clusters. In this paper, we report the discovery of compact (~10"-20") galaxy overdensities centered at the position of three SMGs detected with the Max-Planck Millimeter Bolometer camera (MAMBO) in the COSMOS field. These associations are statistically significant. The photometric redshifts of galaxies in these structures are consistent with their associated SMGs; all of them are between z=1.4-2.5, implying projected physical sizes of ~170 kpc for the overdensities. Our results suggest that about 30% of the radio-identified bright SMGs in that redshift range form in galaxy density peaks in the crucial epoch when most stars formed.

preprint2009arXiv

Massive Galaxies in COSMOS: Evolution of Black hole versus bulge mass but not versus total stellar mass over the last 9 Gyrs?

We constrain the ratio of black hole (BH) mass to total stellar mass of type-1 AGN in the COSMOS survey at 1<z<2. For 10 AGN at mean redshift z~1.4 with both HST/ACS and HST/NICMOS imaging data we are able to compute total stellar mass M_(*,total), based on restframe UV-to-optical host galaxy colors which constrain mass-to-light ratios. All objects have virial BH mass-estimates available from the COSMOS Magellan/IMACS and zCOSMOS surveys. We find zero difference between the M_BH--M_(*,total)-relation at z~1.4 and the M_BH--M_(*,bulge)-relation in the local Universe. Our interpretation is: (a) If our objects were purely bulge-dominated, the M_BH--M_(*,bulge)-relation has not evolved since z~1.4. However, (b) since we have evidence for substantial disk components, the bulges of massive galaxies (logM_(*,total)=11.1+-0.25 or logM_BH~8.3+-0.2) must have grown over the last 9 Gyrs predominantly by redistribution of disk- into bulge-mass. Since all necessary stellar mass exists in the galaxy at z=1.4, no star-formation or addition of external stellar material is required, only a redistribution e.g. induced by minor and major merging or through disk instabilities. Merging, in addition to redistributing mass in the galaxy, will add both BH and stellar/bulge mass, but does not change the overall final M_BH/M_(*,bulge) ratio. Since the overall cosmic stellar and BH mass buildup trace each other tightly over time, our scenario of bulge-formation in massive galaxies is independent of any strong BH-feedback and means that the mechanism coupling BH and bulge mass until the present is very indirect.

preprint2009arXiv

Millimetre observations of a sample of high-redshift obscured quasars

We present observations at 1.2 mm with MAMBO-II of a sample of z>~2 radio-intermediate obscured quasars, as well as CO observations of two sources with the Plateau de Bure Interferometer. Five out of 21 sources (24%) are detected at a significance of >=3sigma. Stacking all sources leads to a statistical detection of <S_1.2mm>= 0.96+-0.11 mJy and stacking only the non-detections also yields a statistical detection, with <S_1.2mm>= 0.51+-0.13 mJy. This corresponds to a typical far-infrared luminosity L_FIR~4x10^12 Lsol. If the far-infrared luminosity is powered entirely by star-formation, and not by AGN-heated dust, then the characteristic inferred star-formation rate is ~700 Msol yr-1. This far-infrared luminosity implies a dust mass of M_dust~3x10^8 Msol. We estimate that such large dust masses on kpc scales can plausibly cause the obscuration of the quasars. We present dust SEDs for our sample and derive a mean SED for our sample. This mean SED is not well fitted by clumpy torus models, unless additional extinction and far-infrared re-emission due to cool dust are included. There is a hint that the host galaxies of obscured quasars must have higher far-infrared luminosities and cool-dust masses and are therefore often found at an earlier evolutionary phase than those of unobscured quasars. For one source at z=2.767, we detect the CO(3-2) transition, with S_CO Delta nu=630+-50 mJy km s-1, corresponding to L_CO(3-2)= 3.2x10^7 Lsol, or L'_CO(3-2)=2.4x10^10 K km s-1 pc2. For another source at z=4.17, the lack of detection of the CO(4-3) line yields a limit of L'_CO(4-3)<1x10^10 K km s-1 pc2. Molecular gas masses, gas depletion timescales and gas-to-dust ratios are estimated (Abridged).

preprint2008arXiv

Dynamical Evolution of AGN Host Galaxies - Gas In/Out-Flow Rates in 7 NUGA Galaxies

To examine the role of the host galaxy structure in fueling nuclear activity, we estimated gas flow rates from several kpc down to the inner few 10 pc for seven nearby spiral galaxies, selected from the NUGA sample (NUclei of GAlaxies). We calculated gravitational torques from near-IR images and determined gas in/out-flow rates as a function of radius and location within the galactic disks, based on high angular resolution interferometric observations of molecular (CO using PdBI) and atomic (HI using the VLA) gas. The results are compared with kinematic evidence for radial gas flows and the dynamical state of the galaxies (via resonances) derived from several different methods. We show that gravitational torques are very efficient at transporting gas from the outer disk all the way into the galaxies centers at ~100 pc; previously assumed dynamical barriers to gas transport, such as the Corotation Resonance of stellar bars, seem to be overcome by gravitational torque induced gas flows from other non-axisymmmetric structures. The resulting rates of gas mass inflow range from 0.01 to 50 solar masses per year and are larger for the galaxy center than for the outer disk. Our gas flow maps show the action of nested bars within larger bars for 3 galaxies. Non-circular streaming motions found in the kinematic maps are larger in the center than in the outer disk and appear to correlate only loosely with the in/out-flow rates as a function of radius. We demonstrate that spiral gas disks are very dynamic systems that undergo strong radial evolution on timescales of a few rotation periods (e.g. 5 times 10^8 yrs at a radius of 5 kpc), due to the effectiveness of gravitational torques in redistributing the cold galactic gas.

preprint2008arXiv

The COSMOS AGN Spectroscopic Survey I: XMM Counterparts

We present optical spectroscopy for an X-ray and optical flux-limited sample of 677 XMM-Newton selected targets covering the 2 deg^2 COSMOS field, with a yield of 485 high-confidence redshifts. The majority of the spectra were obtained over three seasons (2005-2007) with the IMACS instrument on the Magellan (Baade) telescope. We also include in the sample previously published Sloan Digital Sky Survey spectra and supplemental observations with MMT/Hectospec. We detail the observations and classification analyses. The survey is 90% complete to flux limits of f_{0.5-10 keV}>8 x 10^-16 erg cm^-2 s^-1 and i_AB+<22, where over 90% of targets have high-confidence redshifts. Making simple corrections for incompleteness due to redshift and spectral type allows for a description of the complete population to $i_AB+<23. The corrected sample includes 57% broad emission line (Type 1, unobscured) AGN at 0.13<z<4.26, 25% narrow emission line (Type 2, obscured) AGN at 0.07<z<1.29, and 18% absorption line (host-dominated, obscured) AGN at 0<z<1.22 (excluding the stars that made up 4% of the X-ray targets). We show that the survey's limits in X-ray and optical flux include nearly all X-ray AGN (defined by L_{0.5-10 keV}>3 x 10^42 erg s^-1) to z<1, of both optically obscured and unobscured types. We find statistically significant evidence that the obscured to unobscured AGN ratio at z<1 increases with redshift and decreases with luminosity.

preprint2007arXiv

Evolution of the Bar Fraction in COSMOS: Quantifying the Assembly of the Hubble Sequence

We have analyzed the redshift-dependent fraction of galactic bars over 0.2<z<0.84 in 2,157 luminous face-on spiral galaxies from the COSMOS 2-square degree field. Our sample is an order of magnitude larger than that used in any previous investigation, and is based on substantially deeper imaging data than that available from earlier wide-area studies of high-redshift galaxy morphology. We find that the fraction of barred spirals declines rapidly with redshift. Whereas in the local Universe about 65% of luminous spiral galaxies contain bars (SB+SAB), at z ~0.84 this fraction drops to about 20%. Over this redshift range the fraction of strong (SB) bars drops from about 30% to under 10%. It is clear that when the Universe was half its present age, the census of galaxies on the Hubble sequence was fundamentally different from that of the present day. A major clue to understanding this phenomenon has also emerged from our analysis, which shows that the bar fraction in spiral galaxies is a strong function of stellar mass, integrated color and bulge prominence. The bar fraction in very massive, luminous spirals is about constant out to z ~ 0.84 whereas for the low mass, blue spirals it declines significantly with redshift beyond z=0.3. There is also a slight preference for bars in bulge dominated systems at high redshifts which may be an important clue towards the co-evolution of bars, bulges and black holes. Our results thus have important ramifications for the processes responsible for galactic downsizing, suggesting that massive galaxies matured early in a dynamical sense, and not just as a result of the regulation of their star formation rate.