BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

Brent Groves

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

ResearcherUnclaimed source record
astro-ph.GAastro-ph.COastro-ph.IMastro-ph.EPastro-ph.HEastro-ph.SR

Catalog footprint

What is connected

59works
6topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

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

59 published item(s)

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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

H$α$ Morphologies of Star Clusters in 16 LEGUS Galaxies: Constraints on HII region evolution timescales

The analysis of star cluster ages in tandem with the morphology of their HII regions can provide insight into the processes that clear a cluster's natal gas, as well as the accuracy of cluster ages and dust reddening derived from Spectral Energy Distribution (SED) fitting. We classify 3757 star clusters in 16 nearby galaxies according to their H$α$ morphology (concentrated, partially exposed, no emission), using Hubble Space Telescope (HST) imaging from the Legacy ExtraGalactic Ultraviolet Survey (LEGUS). We find: 1) The mean SED ages of clusters with concentrated (1-2 Myr) and partially exposed HII region morphologies (2-3 Myr) indicate a relatively early onset of gas clearing and a short (1-2 Myr) clearing timescale. 2) The reddening of clusters can be overestimated due to the presence of red supergiants, which is a result of stochastic sampling of the IMF in low mass clusters. 3) The age-reddening degeneracy impacts the results of the SED fitting - out of 1408 clusters with M$_*$ $\geq$ 5000 M$_{\odot}$, we find that at least 46 (3%) have SED ages which appear significantly underestimated or overestimated based on H$α$ and their environment, while the total percentage of poor age estimates is expected to be several times larger. 4) Lastly, we examine the dependence of the morphological classifications on spatial resolution. At HST resolution, our conclusions are robust to the distance range spanned by the sample (3-10 Mpc). However, analysis of ground-based H$α$ images shows that compact and partially exposed morphologies frequently cannot be distinguished from each other.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

The Black Hole-Galaxy Connection: Interplay between Feedback, Obscuration, and Host Galaxy Substructure

There is growing evidence for physical influence between supermassive black holes and their host galaxies. We present a case study of nearby galaxy NGC 7582, for which we find evidence that galactic substructure plays an important role in affecting the collimation of ionized outflows as well as contributing to the heavy active galactic nucleus (AGN) obscuration. This result contrasts with a simple, small-scale AGN torus model, according to which AGN wind collimation may take place inside the torus itself, at subparsec scale. Using 3D spectroscopy with the MUSE instrument, we probe the kinematics of the stellar and ionized gas components as well as the ionization state of the gas from a combination of emission line ratios. We report for the first time a kinematically distinct core (KDC) in NGC 7582, on a scale of ~600pc. This KDC coincides spatially with dust lanes and starbursting complexes previously observed. We interpret it as a circumnuclear ring of stars and dusty, gas-rich material. We obtain a clear view of the outflowing cones over kpc scales, and demonstrate that they are predominantly photoionized by the central engine. We detect the back cone (behind the galaxy), and confirm previous results of a large nuclear obscuration of both the stellar continuum and HII regions. While we tentatively associate the presence of the KDC to a large-scale bar and/or a minor galaxy merger, we stress the importance of gaining a better understanding of the role of galaxy substructure in controlling the fueling, feedback and obscuration of AGN.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

The SAMI Galaxy Survey: the difference between ionised gas and stellar velocity dispersions

We investigate the mean locally-measured velocity dispersions of ionised gas ($σ_{\rm gas}$) and stars ($σ_*$) for 1090 galaxies with stellar masses $\log\,(M_*/M_{\odot}) \geq 9.5$ from the SAMI Galaxy Survey. For star-forming galaxies, $σ_*$ tends to be larger than $σ_{\rm gas}$, suggesting that stars are in general dynamically hotter than the ionised gas (asymmetric drift). The difference between $σ_{\rm gas}$ and $σ_*$ ($Δσ$) correlates with various galaxy properties. We establish that the strongest correlation of $Δσ$ is with beam smearing, which inflates $σ_{\rm gas}$ more than $σ_*$, introducing a dependence of $Δσ$ on both the effective radius relative to the point spread function and velocity gradients. The second-strongest correlation is with the contribution of active galactic nuclei (AGN) (or evolved stars) to the ionised gas emission, implying the gas velocity dispersion is strongly affected by the power source. In contrast, using the velocity dispersion measured from integrated spectra ($σ_{\rm aper}$) results in less correlation between the aperture-based $Δσ$ ($Δσ_{\rm aper}$) and the power source. This suggests that the AGN (or old stars) dynamically heat the gas without causing significant deviations from dynamical equilibrium. Although the variation of $Δσ_{\rm aper}$ is much smaller than that of $Δσ$, a correlation between $Δσ_{\rm aper}$ and gas velocity gradient is still detected, implying there is a small bias in dynamical masses derived from stellar and ionised gas velocity dispersions.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

The SAMI Galaxy Survey: The Link Between [$α$/Fe] and Kinematic Morphology

We explore a sample of 1492 galaxies with measurements of the mean stellar population properties and the spin parameter proxy, $λ_{R_{\rm{e}}}$, drawn from the SAMI Galaxy Survey. We fit a global $\left[α/\rm{Fe}\right]$-$σ$ relation, finding that $\left[α/\rm{Fe}\right]=(0.395\pm0.010)\rm{log}_{10}\left(σ\right)-(0.627\pm0.002)$. We observe an anti-correlation between the residuals $Δ\left[α/\rm{Fe}\right]$ and the inclination-corrected $λ_{\,R_{\rm{e}}}^{\rm{\,eo}}$, which can be expressed as $Δ\left[α/\rm{Fe}\right]=(-0.057\pm0.008)λ_{\,R_{\rm{e}}}^{\rm{\,eo}}+(0.020\pm0.003)$. The anti-correlation appears to be driven by star-forming galaxies, with a gradient of $Δ\left[α/\rm{Fe}\right]\sim(-0.121\pm0.015)λ_{\,R_{\rm{e}}}^{\rm{\,eo}}$, although a weak relationship persists for the subsample of galaxies for which star formation has been quenched. We take this to be confirmation that disk-dominated galaxies have an extended duration of star formation. At a reference velocity dispersion of 200 km s$^{-1}$, we estimate an increase in half-mass formation time from $\sim$0.5 Gyr to $\sim$1.2 Gyr from low- to high-$λ_{\,R_{\rm{e}}}^{\rm{\,eo}}$ galaxies. Slow rotators do not appear to fit these trends. Their residual $α$-enhancement is indistinguishable from other galaxies with $λ_{\,R_{\rm{e}}}^{\rm{\,eo}}\lessapprox0.4$, despite being both larger and more massive. This result shows that galaxies with $λ_{\,R_{\rm{e}}}^{\rm{\,eo}}\lessapprox0.4$ experience a similar range of star formation histories, despite their different physical structure and angular momentum.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

The SAMI Galaxy Survey: The relationship between galaxy rotation and the motion of neighbours

Using data from the SAMI Galaxy Survey, we investigate the correlation between the projected stellar kinematic spin vector of 1397 SAMI galaxies and the line-of-sight motion of their neighbouring galaxies. We calculate the luminosity-weighted mean velocity difference between SAMI galaxies and their neighbours in the direction perpendicular to the SAMI galaxies angular momentum axes. The luminosity-weighted mean velocity offsets between SAMI and neighbours, which indicates the signal of coherence between the rotation of the SAMI galaxies and the motion of neighbours, is 9.0 $\pm$ 5.4 km s$^{-1}$ (1.7 $σ$) for neighbours within 1 Mpc. In a large-scale analysis, we find that the average velocity offsets increase for neighbours out to 2 Mpc. However, the velocities are consistent with zero or negative for neighbours outside 3 Mpc. The negative signals for neighbours at distance around 10 Mpc are also significant at $\sim 2$ $σ$ level, which indicate that the positive signals within 2 Mpc might come from the variance of large-scale structure. We also calculate average velocities of different subsamples, including galaxies in different regions of the sky, galaxies with different stellar masses, galaxy type, $λ_{Re}$ and inclination. Although low-mass, high-mass, early-type and low-spin galaxies subsamples show 2 - 3 $σ$ signal of coherence for the neighbours within 2 Mpc, the results for different inclination subsamples and large-scale results suggest that the $\sim 2 σ$ signals might result from coincidental scatter or variance of large-scale structure. Overall, the modest evidence of coherence signals for neighbouring galaxies within 2 Mpc needs to be confirmed by larger samples of observations and simulation studies.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

The SAMI Galaxy Survey: Using concentrated star-formation and stellar population ages to understand environmental quenching

We study environmental quenching using the spatial distribution of current star-formation and stellar population ages with the full SAMI Galaxy Survey. By using a star-formation concentration index [C-index, defined as log10(r_{50,Halpha}/r_{50,cont})], we separate our sample into regular galaxies (C-index>-0.2) and galaxies with centrally concentrated star-formation (SF-concentrated; C-index<-0.2). Concentrated star-formation is a potential indicator of galaxies currently undergoing `outside-in' quenching. Our environments cover ungrouped galaxies, low-mass groups (M_200<10^12.5 M_sun), high-mass groups (M_200 in the range 10^{12.5-14} M_sun) and clusters (M_200>10^14 M_sun). We find the fraction of SF-concentrated galaxies increases as halo mass increases with 9\pm2 per cent, 8\pm3 per cent, 19\pm4 per cent and 29\pm4 per cent for ungrouped galaxies, low-mass groups, high-mass groups and clusters, respectively. We interpret these results as evidence for `outside-in' quenching in groups and clusters. To investigate the quenching time-scale in SF-concentrated galaxies, we calculate light-weighted age (Age_L) and mass-weighted age (Age_M) using full spectral fitting, as well as the Dn4000 and Hdelta_A indices. We assume that the average galaxy age radial profile before entering a group or cluster is similar to ungrouped regular galaxies. At large radius (1-2 R_e), SF-concentrated galaxies in high-mass groups have older ages than ungrouped regular galaxies with an age difference of 1.83\pm0.38 Gyr for Age_L and 1.34\pm0.56 Gyr for Age_M. This suggests that while `outside-in' quenching can be effective in groups, the process will not quickly quench the entire galaxy. In contrast, the ages at 1-2 R_e of cluster SF-concentrated galaxies and ungrouped regular galaxies are consistent (0.19\pm0.21 Gyr for Age_L, 0.40\pm0.61 Gyr for Age_M), suggesting the quenching process must be rapid.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

The SAMI Galaxy Survey: the third and final data release

We have entered a new era where integral-field spectroscopic surveys of galaxies are sufficiently large to adequately sample large-scale structure over a cosmologically significant volume. This was the primary design goal of the SAMI Galaxy Survey. Here, in Data Release 3 (DR3), we release data for the full sample of 3068 unique galaxies observed. This includes the SAMI cluster sample of 888 unique galaxies for the first time. For each galaxy, there are two primary spectral cubes covering the blue (370-570nm) and red (630-740nm) optical wavelength ranges at spectral resolving power of R=1808 and 4304 respectively. For each primary cube, we also provide three spatially binned spectral cubes and a set of standardized aperture spectra. For each galaxy, we include complete 2D maps from parameterized fitting to the emission-line and absorption-line spectral data. These maps provide information on the gas ionization and kinematics, stellar kinematics and populations, and more. All data are available online through Australian Astronomical Optics (AAO) Data Central.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

SH$α$DE: Survey description and mass-kinematics scaling relations for dwarf galaxies

The Study of H$α$ from Dwarf Emissions (SH$α$DE) is a high spectral resolution (R=13500) H$α$ integral field survey of 69 dwarf galaxies with stellar masses $10^6<M_\star<10^9 \,\rm{M_\odot}$. The survey used FLAMES on the ESO Very Large Telescope. SH$α$DE is designed to study the kinematics and stellar populations of dwarf galaxies using consistent methods applied to massive galaxies and at matching level of detail, connecting these mass ranges in an unbiased way. In this paper we set out the science goals of SH$α$DE, describe the sample properties, outline the data reduction and analysis processes. We investigate the $\log{M_{\star}}-\log{S_{0.5}}$ mass-kinematics scaling relation, which have previously shown potential for combining galaxies of all morphologies in a single scaling relation. We extend the scaling relation from massive galaxies to dwarf galaxies, demonstrating this relation is linear down to a stellar mass of $M_{\star}\sim10^{8.6}\,\rm{M_\odot}$. Below this limit, the kinematics of galaxies inside one effective radius appear to be dominated by the internal velocity dispersion limit of the H$α$-emitting gas, giving a bend in the $\log{M_{\star}}-\log{S_{0.5}}$ relation. Replacing stellar mass with total baryonic mass using gas mass estimate reduces the severity but does not remove the linearity limit of the scaling relation. An extrapolation to estimate the galaxies' dark matter halo masses, yields a $\log{M_{h}}-\log{S_{0.5}}$ scaling relation that is free of any bend, has reduced curvature over the whole mass range, and brings galaxies of all masses and morphologies onto the virial relation.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Spinning Bar and a Star-formation Inefficient Repertoire: Turbulence in Hickson Compact Group NGC7674

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

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

The SAMI Galaxy Survey: Decomposed Stellar Kinematics of Galaxy Bulges and Disks

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

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

The SAMI Galaxy Survey: Gas velocity dispersions in low-$z$ star-forming galaxies and the drivers of turbulence

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

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

What Drives the Redshift Evolution of Strong Emission Line Ratios?

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

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

The large-scale ionization cones in the Galaxy

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

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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$\%$.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Dissecting Galaxies: Spatial and Spectral Separation of Emission Excited by Star Formation and AGN Activity

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

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

LZIFU: an emission-line fitting toolkit for integral field spectroscopy data

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

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

The COSMOS [OII] Survey: Evolution of Electron Density with Star Formation Rate

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

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

The Role of Radiation Pressure in the Narrow Line Regions of Seyfert Host Galaxies

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

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

The SAMI Galaxy Survey: extraplanar gas, galactic winds, and their association with star formation history

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

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

The Spatially Resolved [CII] Cooling Line Deficit in Galaxies

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

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

The nebular emission of star-forming galaxies in a hierarchical universe

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

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Variations of the ISM Compactness Across the Main Sequence of Star-Forming Galaxies: Observations and Simulations

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

0 citations0 reviewsNo code linkedOpen access
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

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Untangling the Nature of Spatial Variations of Cold Dust Properties in Star Forming Galaxies

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

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

On the effect of the Cosmic Microwave Background in high-redshift (sub-)millimeter observations

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

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Theoretical Evolution of Optical Strong Lines across Cosmic Time

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

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Ongoing Massive Star Formation in NGC 604

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

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

The dust energy balance in the edge-on spiral galaxy NGC 4565

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

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

The dust SED of dwarf galaxies I. The case of NGC 4214

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

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Resolved optical-infrared SEDs of galaxies: universal relations and their break-down on local scales

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

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

A complete census of AGN and their hosts from optical surveys?

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

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

A Physical Model for z~2 Dust Obscured Galaxies

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

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

ITERA: IDL Tool for Emission-line Ratio Analysis

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

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Optical vs. infrared studies of dusty galaxies and AGN: (I) Nebular emission lines

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

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

High-Resolution Panchromatic Spectral Models of Galaxies including Photoionisation and Dust

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

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

The Formation of High Redshift Submillimeter Galaxies

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

0 citations0 reviewsNo code linkedOpen access