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

Lisa J. Kewley

Lisa J. Kewley 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.HEastro-phastro-ph.SR

Catalog footprint

What is connected

69works
5topics
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

69 published item(s)

preprint2022arXiv

CLASSY III: The Properties of Starburst-Driven Warm Ionized Outflows

We report the results of analyses of galactic outflows in a sample of 45 low-redshift starburst galaxies in the COS Legacy Archive Spectroscopic SurveY (CLASSY), augmented by five additional similar starbursts with COS data. The outflows are traced by blueshifted absorption-lines of metals spanning a wide range of ionization potential. The high quality and broad spectral coverage of CLASSY data enable us to disentangle the absorption due to the static ISM from that due to outflows. We further use different line multiplets and doublets to determine the covering fraction, column density, and ionization state as a function of velocity for each outflow. We measure the outflow's mean velocity and velocity width, and find that both correlate in a highly significant way with the star-formation rate, galaxy mass, and circular velocity over ranges of four orders-of-magnitude for the first two properties. We also estimate outflow rates of metals, mass, momentum, and kinetic energy. We find that, at most, only about 20% of silicon created and ejected by supernovae in the starburst is carried in the warm phase we observe. The outflows' mass-loading factor increases steeply and inversely with both circular and outflow velocity (log-log slope $\sim$ -1.6), and reaches $\sim 10$ for dwarf galaxies. We find that the outflows typically carry about 10 to 100% of the momentum injected by massive stars and about 1 to 20% of the kinetic energy. We show that these results place interesting constraints on, and new insights into, models and simulations of galactic winds.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Gas-phase metallicity break radii of star-forming galaxies in IllustrisTNG

We present radial gas-phase metallicity profiles, gradients, and break radii at redshift $z = 0 - 3$ from the TNG50-1 star-forming galaxy population. These metallicity profiles are characterized by an emphasis on identifying the steep inner gradient and flat outer gradient. From this, the break radius, $r_{\rm Break}$, is defined as the region where the transition occurs. We observe the break radius having a positive trend with mass that weakens with redshift. When normalized by the stellar half-mass radius, the break radius has a weaker relation with both mass and redshift. To test if our results are dependent on the resolution or adopted physics of TNG50-1, the same analysis is performed in TNG50-2 and Illustris-1. We find general agreement between each of the simulations in their qualitative trends; however, the adopted physics between TNG and Illustris differ and therefore the breaks, normalized by galaxy size, deviate by a factor of $\sim$2. In order to understand where the break comes from, we define two relevant time-scales: an enrichment time-scale and a radial gas mixing time-scale. We find that $r_{\rm Break}$ occurs where the gas mixing time-scale is $\sim$10 times as long as the enrichment time-scale in all three simulation runs, with some weak mass and redshift dependence. This implies that galactic disks can be thought of in two-parts: a star-forming inner disk with a steep gradient and a mixing-dominated outer disk with a flat gradient, with the break radius marking the region of transition between them.

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 COS Legacy Archive Spectroscopy SurveY (CLASSY) Treasury Atlas

Far-ultraviolet (FUV; ~1200-2000 angstroms) spectra are fundamental to our understanding of star-forming galaxies, providing a unique window on massive stellar populations, chemical evolution, feedback processes, and reionization. The launch of JWST will soon usher in a new era, pushing the UV spectroscopic frontier to higher redshifts than ever before, however, its success hinges on a comprehensive understanding of the massive star populations and gas conditions that power the observed UV spectral features. This requires a level of detail that is only possible with a combination of ample wavelength coverage, signal-to-noise, spectral-resolution, and sample diversity that has not yet been achieved by any FUV spectral database. We present the COS Legacy Spectroscopic SurveY (CLASSY) treasury and its first high level science product, the CLASSY atlas. CLASSY builds on the HST archive to construct the first high-quality (S/N_1500 >~ 5/resel), high-resolution (R~15,000) FUV spectral database of 45 nearby (0.002 < z < 0.182) star-forming galaxies. The CLASSY atlas, available to the public via the CLASSY website, is the result of optimally extracting and coadding 170 archival+new spectra from 312 orbits of HST observations. The CLASSY sample covers a broad range of properties including stellar mass (6.2 < logM_star(M_sol) < 10.1), star formation rate (-2.0 < log SFR (M_sol/yr) < +1.6), direct gas-phase metallicity (7.0 < 12+log(O/H) < 8.8), ionization (0.5 < O_32 < 38.0), reddening (0.02 < E(B-V < 0.67), and nebular density (10 < n_e (cm^-3) < 1120). CLASSY is biased to UV-bright star-forming galaxies, resulting in a sample that is consistent with z~0 mass-metallicity relationship, but is offset to higher SFRs by roughly 2 dex, similar to z >~2 galaxies. This unique set of properties makes the CLASSY atlas the benchmark training set for star-forming galaxies across cosmic time.

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

Theoretically Modelling Photoionized Regions with Fractal Geometry in Three Dimensions

We create a photoionization model embedded in the turbulent ISM by using the state-of-the-art Messenger Monte-Carlo MAPPINGS~V code (M$^3$) in conjunction with the CMFGEN stellar atmosphere model. We show that the turbulent ISM causes the inhomogeneity of electron temperature and density within the nebula. The fluctuation in the turbulent ISM creates complex ionization structures seen in nearby nebulae. The inhomogeneous density distribution within the nebula creates a significant scatter on the spatially-resolved standard optical diagnostic diagrams, which cannot be represented by the spherical constant density photoionization model. We analyze the dependence of different optical emission lines on the complexity of nebular geometry, finding that the emission-lines residing on the nebular boundary are highly sensitive to the complexity of nebular geometry, while the emission-lines produced throughout the nebula are sensitive to the density distribution of the ISM within the nebula. Our fractal photoionization model demonstrates that a complex nebular geometry is required for accurate modeling of HII regions and emission-line galaxies, especially for the high-redshift galaxies, where the ISM is highly turbulent based on the increasing observational evidence.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

On the origin of nitrogen at low metallicity

Understanding the evolution of the N/O ratio in the interstellar medium (ISM) of galaxies is essential if we are to complete our picture of the chemical evolution of galaxies at high redshift, since most observational calibrations of O/H implicitly depend upon the intrinsic N/O ratio. The observed N/O ratio, however, shows large scatter at low O/H, and is strongly dependent on galactic environment. We show that several heretofore unexplained features of the N/O distribution at low O/H can be explained by the N seen in metal-poor galaxies being mostly primary nitrogen that is returned to the ISM via pre-supernova winds from rapidly rotating massive stars ($M \gtrsim 10$ M$_\odot$, $v/v_{\rm crit} \gtrsim 0.4$). This mechanism naturally produces the observed N/O plateau at low O/H. We show that the large scatter in N/O at low O/H also arises naturally from variations in star-formation efficiency. By contrast, models in which the N and O come primarily from supernovae provide a very poor fit to the observed abundance distribution. We propose that the peculiar abundance patterns we observe at low O/H are a signature that dwarf galaxies retain little of their SN ejecta, leaving them with abundance patterns typical of winds.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Tracing the Ionization Structure of the Shocked Filaments of NGC 6240

We study the ionization and excitation structure of the interstellar medium in the late-stage gas-rich galaxy merger NGC 6240 using a suite of emission line maps at $\sim$25 pc resolution from the Hubble Space Telescope, Keck NIRC2 with Adaptive Optics, and ALMA. NGC 6240 hosts a superwind driven by intense star formation and/or one or both of two active nuclei; the outflows produce bubbles and filaments seen in shock tracers from warm molecular gas (H$_2$ 2.12$μ$m) to optical ionized gas ([O III], [N II], [S II], [O I]) and hot plasma (Fe XXV). In the most distinct bubble, we see a clear shock front traced by high [O III]/H$β$ and [O III]/[O I]. Cool molecular gas (CO(2-1)) is only present near the base of the bubble, towards the nuclei launching the outflow. We interpret the lack of molecular gas outside the bubble to mean that the shock front is not responsible for dissociating molecular gas, and conclude that the molecular clouds are partly shielded and either entrained briefly in the outflow, or left undisturbed while the hot wind flows around them. Elsewhere in the galaxy, shock-excited H$_2$ extends at least $\sim$4 kpc from the nuclei, tracing molecular gas even warmer than that between the nuclei, where the two galaxies' interstellar media are colliding. A ridgeline of high [O III]/H$β$ emission along the eastern arm aligns with the south nucleus' stellar disk minor axis; optical integral field spectroscopy from WiFeS suggests this highly ionized gas is centered at systemic velocity and likely photoionized by direct line-of-sight to the south AGN.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

MOSEL Survey: Tracking the Growth of Massive Galaxies at 2<z<4 using Kinematics and the IllustrisTNG Simulation

We use K-band spectroscopic data from the Multi-Object Spectroscopic Emission Line (MOSEL) survey to analyze the kinematic properties of galaxies at z>3. Our sample consists of 34 galaxies at 3.0<zspec<3.8 between 9.0<log(M_star)<11.0. We find that galaxies with log(M_star) > 10.2 at z > 3 have 56 +/- 21 km/s lower integrated velocity dispersion compared to galaxies at z ~ 2 of similar stellar mass. Massive galaxies at z > 3 have either a flat or declining star formation history (SFH), whereas similar stellar mass galaxies at z~2.0 exhibit a slight peak in the past 500 Myrs. Comparing with the IllustrisTNG cosmological simulation, we find that (i) the dynamical mass of massive galaxies in simulations (log(M_star) > 10.0) increases by ~0.1 dex at a fixed stellar mass between z=2.0-3.0, and (ii) dynamical mass growth is coupled with a rapid rise in the ex-situ stellar mass fraction (stars accreted from other galaxies) for massive galaxies at z < 3.5. We speculate that the rising contribution of ex-situ stellar mass to the total stellar mass growth of massive galaxies is driving the higher integrated velocity dispersion and rising SFHs of massive galaxies at z~2.0 compared to galaxies of similar stellar masses at z > 3.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Oxygen Loss from Simulated Galaxies and the Metal Flow Main Sequence: Predicting the Dependence on Mass and Environment

We predict the mass fraction of oxygen lost from galaxies in a cosmological simulation as a function of stellar mass and environment at the present day. The distribution with stellar mass is bimodal, separating star-forming and quenched galaxies. The metallicity of gas and stars is self-consistently calculated using a chemical evolution model that includes supernovae type II and Ia, hypernovae, and asymptotic giant branch stars. The mass of oxygen lost from each galaxy is calculated by comparing the existing oxygen in gas and stars in the galaxy to the oxygen that should have been produced by the present-day population of stars. More massive galaxies are able to retain a greater fraction of their metals ($\sim 100$ per cent) than low-mass galaxies ($\sim 40 - 70$ per cent). As in the star formation main sequence, star-forming galaxies follow a tight relationship also in terms of oxygen mass lost -- a metal flow main sequence, ZFMS -- whereas massive quenched galaxies tend to have lost a greater fraction of oxygen (up to 20 per cent), due to AGN-driven winds. The amount of oxygen lost by satellite galaxies depends on the details of their interaction history, and those in richer groups tend to have lost a greater fraction of their oxygen. Observational estimates of metal retention in galaxies will provide a strong constraint on models of galaxy evolution.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Quantifying the effects of spatial resolution and noise on galaxy metallicity gradients

Metallicity gradients are important diagnostics of galaxy evolution, because they record the history of events such as mergers, gas inflow and star-formation. However, the accuracy with which gradients can be measured is limited by spatial resolution and noise, and hence measurements need to be corrected for such effects. We use high resolution (~20 pc) simulation of a face-on Milky Way mass galaxy, coupled with photoionisation models, to produce a suite of synthetic high resolution integral field spectroscopy (IFS) datacubes. We then degrade the datacubes, with a range of realistic models for spatial resolution (2 to 16 beams per galaxy scale length) and noise, to investigate and quantify how well the input metallicity gradient can be recovered as a function of resolution and signal-to-noise ratio (SNR) with the intention to compare with modern IFS surveys like MaNGA and SAMI. Given appropriate propagation of uncertainties and pruning of low SNR pixels, we show that a resolution of 3-4 telescope beams per galaxy scale length is sufficient to recover the gradient to ~10-20% uncertainty. The uncertainty escalates to ~60% for lower resolution. Inclusion of the low SNR pixels causes the uncertainty in the inferred gradient to deteriorate. Our results can potentially inform future IFS surveys regarding the resolution and SNR required to achieve a desired accuracy in metallicity gradient measurements.

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

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
preprint2020arXiv

ZFIRE: Measuring Electron Density with [OII] as a function of environment at z = 1.62

The global star formation rates (SFR) of galaxies at fixed stellar masses increase with redshift and are known to vary with environment unto z~2. We explore here whether the changes in the star formation rates can also apply to the electron densities of the inter-stellar medium (ISM) by measuring the [OII] (3727A/3729A) ratio for cluster and field galaxies at z~2. We measure a median electron density of ne = 366+/-84 cm-3 for six galaxies (with 1-sigma scatter = 163 cm-3) in the UDS proton-cluster at z=1.62. We find that the median electron density of galaxies in the UDS photo-cluster environment is three times higher compared to the median electron density of field galaxies (ne = 113+/- 63 cm-3 and 1-sigma scatter = 79 cm-3) at comparable redshifts, stellar mass and SFR. However, we note that a sample of six photo-cluster galaxies is insufficient to reliably measure the electron density in the average porto-cluster environment at z~2. We conclude that the electron density increases with redshift in both cluster and field environments up to z~2 (ne = 30 +/- 1 cm-3 for z ~ 0 to ne =254+/- 76 cm-3 for z~1.5). We find tentative evidence (~2.6 sigma ) for a possible dependence of electron density on environment, but the results require confirmation with larger sample sizes.

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
preprint2016arXiv

A physically-based model of the ionizing radiation from active galaxies for photoionization modeling

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

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Chemical abundances in high-redshift galaxies: A powerful new emission line diagnostic

This Letter presents a new, remarkably simple diagnostic specifically designed to derive chemical abundances for high redshift galaxies. It uses only the H α, [N II] and [S II] emission lines, which can usually be observed in a single gating stetting, and is almost linear up to an abundance of 12+log (O/H) = 9.05. It can be used over the full abundance range encountered in high redshift galaxies. By its use of emission lines located close together in wavelength, it is also independent of reddening. Our diagnostic depends critically on the calibration of the N/O ratio. However, by using realistic stellar atmospheres combined with the N/O vs. O/H abundance calibration derived locally from stars and H II regions, and allowing for the fact that high-redshift H II regions have both high ionisation parameters \emph{and} high gas pressures, we find that the observations of high-redshift galaxies can be simply explained by the models without having to invoke arbitrary changes in N/O ratio, or the presence of unusual quantities of Wolf-Rayet stars in these galaxies.

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

Quenching star formation: insights from the local main sequence

The so-called star-forming main sequence of galaxies is the apparent tight relationship between the star formation rate and stellar mass of a galaxy. Previous studies exclude galaxies which are not strictly 'star forming' from the main sequence, because they do not lie on the same tight relation. Using local galaxies in the Sloan Digital Sky Survey we have classified galaxies according to their emission line ratios, and studied their location on the star formation rate - stellar mass plane. We find that galaxies form a sequence from the `blue cloud' galaxies which are actively forming stars, through a combination of composite, Seyfert, and LINER (Low-ionization nuclear emission-line region) galaxies, ending as 'red-and-dead' galaxies. The sequence supports an evolutionary pathway for galaxies in which star formation quenching by active galactic nuclei (AGN) plays a key role.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Radial Distribution Of ISM Gas-phase Metallicity In CLASH Clusters at z~0.35: A New Outlook On Environmental Impact On Galaxy Evolution

We present the first observation of cluster-scale radial metallicity gradients from star-forming galaxies. We use the DEIMOS spectrograph on the Keck II telescope to observe two CLASH clusters at z~0.35: MACS1115+0129 and RXJ1532+3021. Based on our measured interstellar medium (ISM) properties of star-forming galaxies out to a radius of 2.5 Mpc from the cluster centre, we find that the galaxy metallicity decreases as a function of projected cluster-centric distance (-0.15+/-0.08 dex/Mpc) in MACS1115+01. On the mass-metallicity relation (MZR), star-forming galaxies in MACS1115+01 are offset to higher metallicity (~0.2 dex) than the local SDSS galaxies at a fixed mass range. In contrast, the MZR of RXJ1532+30 is consistent with the local comparison sample. RXJ1532+30 exhibits a bimodal radial metallicity distribution, with one branch showing a similar negative gradient as MACS1115+01 (-0.14+/-0.05 dex/Mpc) and the other branch showing a positive radial gradient. The positive gradient branch in RXJ1532+30 is likely caused by either interloper galaxies or an in-plane merger, indicating that cluster-scale abundance gradients probe cluster substructures and thus the dynamical state of a cluster. Most strikingly, we discover that neither the radial metallicity gradient nor the offset from the MZR is driven by the stellar mass. We compare our observations with Rhapsody-G cosmological hydrodynamical zoom-in simulations of relaxed galaxy clusters and find that the simulated galaxy cluster also exhibits a negative abundance gradient, albeit with a shallower slope (-0.04+/-0.03 dex/Mpc). Our observations suggest that the negative radial gradient originates from ram-pressure stripping and/or strangulation processes in the cluster environments.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Shocked POststarburst Galaxy Survey II: The Molecular Gas Content and Properties of a Subset of SPOGs

We present CO(1-0) observations of objects within the Shocked POststarburst Galaxy Survey taken with the Institut de Radioastronomie Millimetrique (IRAM) 30m single dish and the Combined Array for Research for Millimeter Astronomy (CARMA) interferometer. Shocked Poststarburst Galaxies (SPOGs) represent a transitioning population of galaxies, with deep Balmer absorption (Hdelta>5A), consistent with an intermediate-age (A-star) stellar population, and ionized gas line ratios inconsistent with pure star formation. The CO(1-0) subsample was selected from SPOGs detected by the Wide-field Infrared Survey Explorer with 22um flux detected at a signal-to-noise (S/N)>3. Of the 52 objects observed in CO(1-0), 47 are detected with S/N>3. A large fraction (37-46%) of our CO-SPOG sample were visually classified as morphologically disrupted. The H2 masses detected were between 10^(8.7-10.8) Msuns, consistent with the gas masses found in normal galaxies, though approximately an order of magnitude larger than the range seen in poststarburst galaxies. When comparing the 22um and CO(1-0) fluxes, SPOGs diverge from the normal star-forming relation, having 22um fluxes in excess by a factor of <EMIR,SPOG>=4.91+0.42-0.39. The Na I D characteristics of CO-SPOGs show that it is likely that many of these objects host interstellar winds. Objects with the large Na I D enhancements also tend to emit in the radio, suggesting possible AGN-driving of neutral winds.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Shocked POststarbust Galaxy Survey I: Candidate Poststarbust Galaxies with Emission Line Ratios Consistent with Shocks

[Abridged] The Shocked POststarburst Galaxy Survey (SPOGS) aims to identify transforming galaxies, in which the nebular lines are excited via shocks instead of through star formation processes. Utilizing the OSSY measurements on the Sloan Digital Sky Survey Data Release 7 catalog, we applied Balmer absorption and shock boundary criteria to identify 1,067 SPOG candidates (SPOGs*) within z=0.2. SPOGs* represent 0.2% of the OSSY sample galaxies that exceed the continuum signal-to-noise cut (and 0.7% of the emission line galaxy sample). SPOGs* colors suggest that they are in an earlier phase of transition than OSSY galaxies that meet an E+A selection. SPOGs* have a 13% 1.4GHz detection rate from the Faint Images of the Radio Sky at Twenty centimeters survey, higher than most other subsamples, and comparable only to low-ionization nuclear emission line region hosts, suggestive of the presence of active galactic nuclei. SPOGs* also have stronger NaD absorption than predicted from the stellar population, suggestive of cool gas being driven out in galactic winds. It appears that SPOGs* represent an earlier phase in galaxy transformation than traditionally selected poststarburst galaxies, and that a large proportion of SPOGs* also have properties consistent with disruption of their interstellar media, a key component to galaxy transformation. It is likely that many of the known pathways to transformation undergo a SPOG phase. Studying this sample of SPOGs* further, including their morphologies, active galactic nuclei properties, and environments, has the potential for us to build a more complete picture of the initial conditions that can lead to a galaxy evolving.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

The Discovery of Raman Scattering in HII Regions

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

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

ZFIRE: The Kinematics of Star-Forming Galaxies as a Function of Environment at z~2

We perform a kinematic analysis of galaxies at $z\sim2$ in the COSMOS legacy field using near-infrared (NIR) spectroscopy from Keck/MOSFIRE as part of the ZFIRE survey. Our sample consists of 75 Ks-band selected star-forming galaxies from the ZFOURGE survey with stellar masses ranging from log(M$_{\star}$/M$_{\odot}$)$=9.0-11.0$, 28 of which are members of a known overdensity at $z=2.095$. We measure H$α$ emission-line integrated velocity dispersions ($σ_{\rm int}$) from 50$-$230 km s$^{-1}$, consistent with other emission-line studies of $z\sim2$ field galaxies. From these data we estimate virial, stellar, and gas masses and derive correlations between these properties for cluster and field galaxies at $z\sim2$. We find evidence that baryons dominate within the central effective radius. However, we find no statistically significant differences between the cluster and the field, and conclude that the kinematics of star-forming galaxies at $z\sim2$ are not significantly different between the cluster and field environments.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

A rise in the ionizing photons in star-forming galaxies over the past 5 billion years

We investigate the change in ionizing photons in galaxies between 0.2<z<0.6 using the F2 field of the SHELS complete galaxy redshift survey. We show, for the first time, that while the [OIII]/Hb and [OIII]/[OII] ratios rise, the [NII]/H-alpha and [SII]/H-alpha ratios fall significantly over the 0.2<z<0.35 redshift range for stellar masses between 9.2<log(M/Msun)<10.6. The [OIII]/H-beta and [OIII]/[OII] ratios continue to rise across the full 0.2<z<0.6 redshift range for stellar masses between 9.8<log(M/Msun)<10.0. We conclusively rule out AGN contamination, a changing ISM pressure, and a change in the hardness of the EUV radiation field as the cause of the change in the line ratios between 0.2<z<0.35. We find that the ionization parameter rises significantly with redshift (by 0.1 to 0.25 dex depending on the stellar mass of the sample). We show that the ionization parameter is strongly correlated with the fraction of young-to-old stars, as traced by the H-beta equivalent width. We discuss the implications of this result on higher redshift studies, and we consider the implications on the use of standard optical metallicity diagnostics at high redshift.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

A Turnover in the Galaxy Main Sequence of Star Formation at $M_{*} \sim 10^{10} M_{\odot}$ for Redshifts $z < 1.3$

The relationship between galaxy star formation rates (SFR) and stellar masses ($M_\ast$) is re-examined using a mass-selected sample of $\sim$62,000 star-forming galaxies at $z \le 1.3$ in the COSMOS 2-deg$^2$ field. Using new far-infrared photometry from $Herschel$-PACS and SPIRE and $Spitzer$-MIPS 24 $μ$m, along with derived infrared luminosities from the NRK method based on galaxies' locations in the restframe color-color diagram $(NUV - r)$ vs. $(r - K)$, we are able to more accurately determine total SFRs for our complete sample. At all redshifts, the relationship between median $SFR$ and $M_\ast$ follows a power-law at low stellar masses, and flattens to nearly constant SFR at high stellar masses. We describe a new parameterization that provides the best fit to the main sequence and characterizes the low mass power-law slope, turnover mass, and overall scaling. The turnover in the main sequence occurs at a characteristic mass of about $M_{0} \sim 10^{10} M_{\odot}$ at all redshifts. The low mass power-law slope ranges from 0.9-1.3 and the overall scaling rises in SFR as a function of $(1+z)^{4.12 \pm 0.10}$. A broken power-law fit below and above the turnover mass gives relationships of $SFR \propto M_{*}^{0.88 \pm 0.06}$ below the turnover mass and $SFR \propto M_{*}^{0.27 \pm 0.04}$ above the turnover mass. Galaxies more massive than $M_\ast \gtrsim 10^{10}\ M_{\rm \odot}$ have on average, a much lower specific star formation rate (sSFR) than would be expected by simply extrapolating the traditional linear fit to the main sequence found for less massive galaxies.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

A universal, turbulence-regulated star formation law: from Milky Way clouds to high-redshift disk and starburst galaxies

Whilst the star formation rate (SFR) of molecular clouds and galaxies is key in understanding galaxy evolution, the physical processes which determine the SFR remain unclear. This uncertainty about the underlying physics has resulted in various different star formation laws, all having substantial intrinsic scatter. Extending upon previous works that define the column density of star formation (Sigma_SFR) by the gas column density (Sigma_gas), we develop a new universal star formation (SF) law based on the multi-freefall prescription of gas. This new SF law relies predominantly on the probability density function (PDF) and on the sonic Mach number of the turbulence in the star-forming clouds. By doing so we derive a relation where the star formation rate (SFR) correlates with the molecular gas mass per multi-freefall time, whereas previous models had used the average, single-freefall time. We define a new quantity called maximum (multi-freefall) gas consumption rate (MGCR) and show that the actual SFR is only about 0.4% of this maximum possible SFR, confirming the observed low efficiency of star formation. We show that placing observations in this new framework (Sigma_SFR vs. MGCR) yields a significantly improved correlation with 3-4 times reduced scatter compared to previous SF laws and a goodness-of-fit parameter R^2=0.97. By inverting our new relationship, we provide sonic Mach number predictions for kpc-scale observations of Local Group galaxies as well as unresolved observations of local and high-redshift disk and starburst galaxies that do not have independent, reliable estimates for the turbulent cloud Mach number.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

HST/WFC3 Observations of an Off-Nuclear Superbubble in Arp 220

We present a high spatial resolution optical and infrared study of the circumnuclear region in Arp 220, a late-stage galaxy merger. Narrowband imaging using HST/WFC3 has resolved the previously observed peak in H$α$+[NII] emission into a bubble-shaped feature. This feature measures 1.6" in diameter, or 600 pc, and is only 1" northwest of the western nucleus. The bubble is aligned with the western nucleus and the large-scale outflow axis seen in X-rays. We explore several possibilities for the bubble origin, including a jet or outflow from a hidden active galactic nucleus (AGN), outflows from high levels of star formation within the few hundred pc nuclear gas disk, or an ultraluminous X-ray source. An obscured AGN or high levels of star formation within the inner $\sim$100 pc of the nuclei are favored based on the alignment of the bubble and energetics arguments.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Kinematic classifications of local interacting galaxies: implications for the merger/disk classifications at high-z

The classification of galaxy mergers and isolated disks is key for understanding the relative importance of galaxy interactions and secular evolution during the assembly of galaxies. The kinematic properties of galaxies as traced by emission lines have been used to suggest the existence of a significant population of high-z star-forming galaxies consistent with isolated rotating disks. However, recent studies have cautioned that post-coalescence mergers may also display disk-like kinematics. To further investigate the robustness of merger/disk classifications based on kinematic properties, we carry out a systematic classification of 24 local (U)LIRGs spanning a range of galaxy morphologies: from isolated spiral galaxies, ongoing interacting systems, to fully merged remnants. We artificially redshift the WiFeS observations of these local (U)LIRGs to z=1.5 to make a realistic comparison with observations at high-z, and also to ensure that all galaxies have the same spatial sampling of ~900 pc. Using both kinemetry-based and visual classifications, we find that the reliability of kinematic classification shows a strong trend with the interaction stage of galaxies. Mergers with two nuclei and tidal tails have the most distinct kinematic properties compared to isolated disks, whereas a significant population of the interacting disks and merger remnants are indistinguishable from isolated disks. The high fraction of late-stage mergers showing disk-like kinematics reflects the complexity of the dynamics during galaxy interactions. However, the exact fractions of misidentified disks and mergers depend on the definition of kinematic asymmetries and the classification threshold when using kinemetry-based classifications. Our results suggest that additional indicators such as morphologies traced by stars or molecular gas are required to further constrain the merger/disk classifications at high-z.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Metallicity gradients in local field star-forming galaxies: Insights on inflows, outflows, and the coevolution of gas, stars and metals

We present metallicity gradients in 49 local field star-forming galaxies. We derive gas-phase oxygen abundances using two widely adopted metallicity calibrations based on the [OIII]/Hbeta, [NII]/Halpha and [NII]/[OII] line ratios. The two derived metallicity gradients are usually in good agreement within +/-0.14 dex/R25 (R25 is the B-band iso-photoal radius), but the metallicity gradients can differ significantly when the ionisation parameters change systematically with radius. We investigate the metallicity gradients as a function of stellar mass (8<log(M*/Msun)<11) and absolute B-band luminosity (-16 > MB > -22). When the metallicity gradients are expressed in dex/kpc, we show that galaxies with lower mass and luminosity, on average, have steeper metallicity gradients. When the metallicity gradients are expressed in dex/R25, we find no correlation between the metallicity gradients, and stellar mass and luminosity. We provide a local benchmark metallicity gradient of field star-forming galaxies useful for comparison with studies at high redshifts. We investigate the origin of the local benchmark gradient using simple chemical evolution models and observed gas and stellar surface density profiles in nearby field spiral galaxies. Our models suggest that the local benchmark gradient is a direct result of the coevolution of gas and stellar disk under virtually closed-box chemical evolution when the stellar-to-gas mass ratio becomes high (>>0.3). These models imply low current mass accretion rates (<0.3xSFR), and low mass outflow rates (<3xSFR) in local field star-forming galaxies.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Shocked Gas in IRAS F17207-0014: ISM Collisions and Outflows

We combine optical and near-infrared AO-assisted integral field observations of the merging ULIRG IRAS F17207-0014 from the Wide-Field Spectrograph (WiFeS) and Keck/OSIRIS. The optical emission line ratios [N II]/H$α$, [S II]/H$α$, and [O I]/H$α$ reveal a mixing sequence of shocks present throughout the galaxy, with the strongest contributions coming from large radii (up to 100% at $\sim$5 kpc in some directions), suggesting galactic-scale winds. The near-infrared observations, which have approximately 30 times higher spatial resolution, show that two sorts of shocks are present in the vicinity of the merging nuclei: low-level shocks distributed throughout our field-of-view evidenced by an H$_{2}$/Br$γ$ line ratio of $\sim$0.6-4, and strong collimated shocks with a high H$_{2}$/Br$γ$ line ratio of $\sim$4-8, extending south from the two nuclear disks approximately 400 pc ($\sim$0.5 arcsec). Our data suggest that the diffuse shocks are caused by the collision of the interstellar media associated with the two progenitor galaxies and the strong shocks trace the base of a collimated outflow coming from the nucleus of one of the two disks.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

The Absence of an Environmental Dependence in the Mass-Metallicity Relation at z=2

We investigate the environmental dependence of the mass-metallicity relation at z=2 with MOSFIRE/Keck as part of the ZFIRE survey. Here, we present the chemical abundance of a Virgo-like progenitor at z=2.095 that has an established red sequence. We identified 43 cluster ($<z>=2.095\pm0.004$) and 74 field galaxies ($<z>=2.195\pm0.083$) for which we can measure metallicities. For the first time, we show that there is no discernible difference between the mass-metallicity relation of field and cluster galaxies to within 0.02dex. Both our field and cluster galaxy mass-metallicity relations are consistent with recent field galaxy studies at z~2. We present hydrodynamical simulations for which we derive mass-metallicity relations for field and cluster galaxies. We find at most a 0.1dex offset towards more metal-rich simulated cluster galaxies. Our results from both simulations and observations are suggestive that environmental effects, if present, are small and are secondary to the ongoing inflow and outflow processes that are governed by galaxy halo mass.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

The SAMI Galaxy Survey: Shocks and Outflows in a normal star-forming galaxy

We demonstrate the feasibility and potential of using large integral field spectroscopic surveys to investigate the prevalence of galactic-scale outflows in the local Universe. Using integral field data from SAMI and the Wide Field Spectrograph, we study the nature of an isolated disk galaxy, SDSS J090005.05+000446.7 (z = 0.05386). In the integral field datasets, the galaxy presents skewed line profiles changing with position in the galaxy. The skewed line profiles are caused by different kinematic components overlapping in the line-of-sight direction. We perform spectral decomposition to separate the line profiles in each spatial pixel as combinations of (1) a narrow kinematic component consistent with HII regions, (2) a broad kinematic component consistent with shock excitation, and (3) an intermediate component consistent with shock excitation and photoionisation mixing. The three kinematic components have distinctly different velocity fields, velocity dispersions, line ratios, and electron densities. We model the line ratios, velocity dispersions, and electron densities with our MAPPINGS IV shock and photoionisation models, and we reach remarkable agreement between the data and the models. The models demonstrate that the different emission line properties are caused by major galactic outflows that introduce shock excitation in addition to photoionisation by star-forming activities. Interstellar shocks embedded in the outflows shock-excite and compress the gas, causing the elevated line ratios, velocity dispersions, and electron densities observed in the broad kinematic component. We argue from energy considerations that, with the lack of a powerful active galactic nucleus, the outflows are likely to be driven by starburst activities. Our results set a benchmark of the type of analysis that can be achieved by the SAMI Galaxy Survey on large numbers of galaxies.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Z-FIRE: ISM properties of the z = 2.095 COSMOS Cluster

We investigate the ISM properties of 13 star-forming galaxies within the z~2 COSMOS cluster. We show that the cluster members have [NII]/Ha and [OIII]/Hb emission-line ratios similar to z~2 field galaxies, yet systematically different emission-line ratios (by ~0.17 dex) from the majority of local star-forming galaxies. We find no statistically significant difference in the [NII]/Ha and [OIII]/Hb line ratios or ISM pressures among the z~2 cluster galaxies and field galaxies at the same redshift. We show that our cluster galaxies have significantly larger ionization parameters (by up to an order of magnitude) than local star-forming galaxies. We hypothesize that these high ionization parameters may be associated with large specific star formation rates (i.e. a large star formation rate per unit stellar mass). If this hypothesis is correct, then this relationship would have important implications for the geometry and/or the mass of stars contained within individual star clusters as a function of redshift.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

ZFIRE: Galaxy Cluster Kinematics, H$α$ Star Formation Rates, and Gas-Phase Metallicities of XMM-LSS J02182-05102 at z=1.6233

We spectroscopically survey the galaxy cluster XMM-LSS J02182-05102 (hereafter IRC 0218) using LRIS (optical) and MOSFIRE (near-infrared) on Keck I as part of the ZFIRE survey. IRC 0218 has a narrow redshift range of $1.612<z_{\rm spec}<1.635$ defined by 33 members of which 20 are at R$_{\rm proj}<1$ Mpc. The cluster redshift and velocity dispersion are $z_{\rm cl}=1.6233\pm0.0003$ and $σ_{\rm cl}=254\pm50$ km s$^{-1}$. We reach NIR line sensitivities of $\sim0.3\times10^{-17}$ erg s$^{-1}$ cm$^{-2}$ that, combined with multi-wavelength photometry, provide extinction-corrected H$α$ star formation rates (SFR), gas phase metallicities from [NII]/H$α$, and stellar masses. We measure an integrated H$α$ SFR of $\sim325{\rm M}_{\odot}$ yr$^{-1}$ (26 members; R$_{\rm proj}<2$ Mpc) and show that the elevated star formation in the cluster core (R$_{\rm proj}<0.25$ Mpc) is driven by the concentration of star-forming members, but the average SFR per H$α$-detected galaxy is half that of members at R$_{\rm proj}\sim1$ Mpc. However, we do not detect any environmental imprint when comparing attenuation and gas phase metallicities: the cluster galaxies show similar trends with M$_{\star}$ as to the field, e.g. more massive galaxies have larger stellar attenuation. IRC 0218's gas phase metallicity-M$_{\star}$ relation (MZR) is offset to lower metallicities relative to $z\sim0$ and has a slope of $0.13\pm0.10$. Comparing the MZR in IRC 0218 to the COSMOS cluster at $z=2.1$ shows no evolution ($Δt\sim1$ Gyr): the MZR for both galaxy clusters are remarkably consistent with each other and virtually identical to several field surveys at $z\sim2$.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Catching Quenching Galaxies: The Nature of the WISE Infrared Transition Zone

We present the discovery of a prominent bifurcation between early-type galaxies and late-type galaxies, in [4.6]-[12] micron colors from the Wide Field Infrared Survey Explorer (WISE). We then use an emission-line diagnostic comparison sample to explore the nature of objects found both within, and near the edges of, this WISE infrared transition zone (IRTZ). We hypothesize that this birfurcation might be due to the presence of hot dust and PAH emission features in late-type galaxies. Using a sample of galaxies selected through the Shocked Poststarburst Galaxy Survey (SPOGS), we are able to identify galaxies with strong Balmer absorption (EW(Hdelta)>5 Angstroms) as well as emission lines inconsistent with star formation (deemed SPOG candidates, or SPOGs*) that lie within the optical green valley. Seyferts and low ionization nuclear emission line regions, whose u-r colors tend to be red, are strongly represented within the IRTZ, whereas SPOGs* tend to sit near the star-forming edge. Although AGN are well-represented in the IRTZ, we argue that the dominant IRTZ population are galaxies that are in late stages of transitioning across the optical green valley, shedding the last of their remnant interstellar media.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Galaxy emission line classification using 3D line ratio diagrams

Two-dimensional (2D) line ratio diagnostic diagrams have become a key tool in understanding the excitation mechanisms of galaxies. The curves used to separate the different regions - HII-like or else excited by an active galactic nucleus (AGN) - have been refined over time but the core technique has not evolved significantly. However, the classification of galaxies based on their emission line ratios really is a multi-dimensional problem. Here we exploit recent software developments to explore the potential of three-dimensional (3D) line ratio diagnostic diagrams. We introduce a specific set of 3D diagrams, the ZQE diagrams, which separate the oxygen abundance and the ionisation parameter of HII region-like spectra, and which also enable us to probe the excitation mechanism of the gas. By examining these new 3D spaces interactively, we define a new set of 2D diagnostics, the ZE diagnostics, which can provide the metallicity of objects excited by hot young stars, and which cleanly separate HII region-like objects from the different classes of AGNs. We show that these ZE diagnostics are consistent with the key log[NII]/H$α$ vs. log[OIII]/H$β$ diagnostic currently used by the community. They also have the advantage of attaching a probability that a given object belongs to one class or to the other. Finally, we discuss briefly why ZQE diagrams can provide a new way to differentiate and study the different classes of AGNs in anticipation of a dedicated follow-up study.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

KECK/MOSFIRE spectroscopic confirmation of a Virgo-like cluster ancestor at z=2.095

We present the spectroscopic confirmation of a galaxy cluster at $z=2.095$ in the COSMOS field. This galaxy cluster was first reported in the ZFOURGE survey as harboring evolved massive galaxies using photometric redshifts derived with deep near-infrared (NIR) medium-band filters. We obtain medium resolution ($R \sim$ 3600) NIR spectroscopy with MOSFIRE on the Keck 1 telescope and secure 180 redshifts in a $12'\times12'$ region. We find a prominent spike of 57 galaxies at $z=2.095$ corresponding to the galaxy cluster. The cluster velocity dispersion is measured to be $σ_{\rm v1D}$ = 552 $\pm$ 52 km/s. This is the first study of a galaxy cluster in this redshift range ($z \gt 2.0$) with the combination of spectral resolution ($\sim$26 km/s) and the number of confirmed members (${>}50$) needed to impose a meaningful constraint on the cluster velocity dispersion and map its members over a large field of view. Our $Λ$CDM cosmological simulation suggests that this cluster will most likely evolve into a Virgo-like cluster with ${\rm M_{vir}}{=}10^{14.4\pm0.3} {\rm M_\odot}$ ($68\%$ confidence) at $z\sim$ 0. The theoretical expectation of finding such a cluster is $\sim$ $4\%$. Our results demonstrate the feasibility of studying galaxy clusters at $z > 2$ in the same detailed manner using multi-object NIR spectrographs as has been done in the optical in lower redshift clusters.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Metal-poor dwarf galaxies in the SIGRID galaxy sample. I. HII region observations and chemical abundances

In this paper we present the results of observations of seventeen HII regions in thirteen galaxies from the SIGRID sample of isolated gas rich irregular dwarf galaxies. The spectra of all but one of the galaxies exhibit the auroral [OIII] 4363A line, from which we calculate the electron temperature, Te, and gas-phase oxygen abundance. Five of the objects are blue compact dwarf (BCD) galaxies, of which four have not previously been analysed spectroscopically. We include one unusual galaxy which exhibits no evidence of the [NII] λλ 6548,6584A lines, suggesting a particularly low metallicity (< Zsolar/30). We compare the electron temperature based abundances with those derived using eight of the new strong line diagnostics presented by Dopita et al. (2013). Using a method derived from first principles for calculating total oxygen abundance, we show that the discrepancy between the Te-based and strong line gas-phase abundances have now been reduced to within ~0.07 dex. The chemical abundances are consistent with what is expected from the luminosity-metallicity relation. We derive estimates of the electron densities and find them to be between ~5 and ~100 cm-3. We find no evidence for a nitrogen plateau for objects in this sample with metallicities 0.5 > Zsolar > 0.15.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Metal-poor dwarf galaxies in the SIGRID galaxy sample. II. The electron temperature-abundance calibration and the parameters that affect it

In this paper, we use the Mappings photoionization code to explore the physical parameters that impact on the measurement of electron temperature and abundance in HII regions. In the previous paper we presented observations and measurements of physical properties from the spectra of seventeen HII regions in fourteen isolated dwarf irregular galaxies from the SIGRID sample. Here, we analyze these observations further, together with three additional published data sets. We explore the effects of optical thickness, electron density, ionization parameter, ionization source, and non-equilibrium effects on the relation between electron temperature and metallicity. We present a standard model that fits the observed data remarkably well at metallicities between 1/10 and 1 solar. We investigate the effects of optically thin HII regions, and show that they can have a considerable effect on the measured electron temperature, and that there is evidence that some of the observed objects are optically thin. We look at the role of the ionization parameter and find that lower ionization parameter values give better fits at higher oxygen abundance. We show that higher pressures combined with low optical depth, and also kappa electron energy distributions at low kappa values, can generate the apparent high electron temperatures in low metallicity HII regions, and that the former provides the better fit to observations. We examine the effects of these parameters on the strong line diagnostic methods. We extend this to three-dimensional diagnostic grids to confirm how well the observations are described by the grids.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Probing the physics of narrow-line regions of Seyfert galaxies I: The case of NGC 5427

We have used the Wide Field Spectrograph (WiFeS) on the ANU 2.3m telescope at Siding Spring to observe the nearby, nearly face-on, Seyfert 2 galaxy, NGC 5427. We have obtained integral field spectroscopy of both the nuclear regions and the HII regions in the spiral arms. We have constrained the chemical abundance in the interstellar medium of the extended narrow line region (ENLR) by measuring the abundance gradient in the circum-nuclear \ion{H}{ii} regions to determine the nuclear chemical abundances, and to use these to in turn determine the EUV spectral energy distribution for comparison with theoretical models. We find a very high nuclear abundance, $\sim 3.0$ times solar, with clear evidence of a nuclear enhancement of N and He, possibly caused by massive star formation in the extended ($\sim 100$pc) central disk structure. The circum-nuclear narrow-line region spectrum is fit by a radiation pressure dominated photoionisation model model with an input EUV spectrum from a Black Hole with mass $5\times10^7 M_{\odot}$ radiating at $\sim 0.1$ of its Eddington luminosity. The bolometric luminosity is closely constrained to be $\log L_{\mathrm bol.} = 44.3\pm 0.1$ erg s$^{-1}$. The EUV spectrum characterised by a soft accretion disk and a harder component extending to above 15keV. The ENLR region is extended in the NW-SE direction. The line ratio variation in circum-nuclear spaxels can be understood as the result of mixing \ion{H}{ii} regions with an ENLR having a radius-invariant spectrum.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

S7 : Probing the physics of Seyfert Galaxies through their ENLR & HII Regions

Here we present the first results from the Siding Spring Southern Seyfert Spectroscopic Snapshot Survey (S7) which aims to investigate the physics of ~140 radio-detected southern active Galaxies with z<0.02 through Integral Field Spectroscopy using the Wide Field Spectrograph (WiFeS). This instrument provides data cubes of the central 38 x 25 arc sec. of the target galaxies in the waveband 340-710nm with the unusually high resolution of R=7000 in the red (530-710nm), and R=3000 in the blue (340-560nm). These data provide the morphology, kinematics and the excitation structure of the extended narrow-line region, probe relationships with the black hole characteristics and the host galaxy, measures host galaxy abundance gradients and the determination of nuclear abundances from the HII regions. From photoionisation modelling, we may determine the shape of the ionising spectrum of the AGN, discover whether AGN metallicities differ from nuclear abundances determined from HII regions, and probe grain destruction in the vicinity of the AGN. Here we present some preliminary results and modelling of both Seyfert galaxies observed as part of the survey.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Starburst-AGN mixing: I. NGC 7130

We present an integral field spectroscopic study of the Luminous Infrared Galaxy (LIRG) NGC 7130, a known starburst-AGN composite galaxy. We employ standard emission line ratio diagnostics and maps of velocity dispersion and velocity field to investigate how the dominant ionising sources change as a function of radius. From the signatures of both star formation and AGN activity we show that NGC 7130 is a remarkably clean case of starburst-AGN mixing. We find a smooth transition from AGN dominated emission in the centre to pure star forming activity further out, from which we can estimate the radius of the extended narrow line region to be 1.8 +/- 0.8 kpc. We calculate that the fraction of [OIII] luminosity due to star formation and AGN activity is 30 +/- 2% and 70 +/- 3% respectively, and that the fraction of H-alpha luminosity due to star formation and AGN activity is 65 +/- 3% and 35 +/- 2% respectively. We conclude with a discussion of the importance and potential of starburst-AGN mixing for future studies of the starburst-AGN connection.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Starburst-AGN mixing: II. Optically-selected active galaxies

We use 4 galaxies from the Calar Alto Legacy Integral Field Area (CALIFA) survey with clear signs of accretion onto supermassive black holes to investigate the relative contribution of star-formation and active galactic nucleus (AGN) activity to the line-emission of each galaxy as a function of radius. The combination of star-formation and AGN activity produces curved "mixing sequences" on standard optical diagnostic diagrams, and the fraction of emission due to AGN activity decreases smoothly with distance from the centre of the galaxy. We use the AGN activity profiles to calculate the size of the AGN narrow line regions, which have radii of ~ 6.3 kpc. We calculate the fractional contribution of the star-formation and the AGN activity to the global Halpha, [O II] $λλ$ 3727,3729 and [O III] $λ$ 5007 luminosities of each galaxy, and show that both ionization sources contribute significantly to the emission in all three lines. We use weighted combinations of stellar and AGN photoionization models to produce mixing models, which are consistent with observations for 85 percent of spaxels across the galaxies in our sample. We also present a new diagnostic for starburst-AGN mixing which requires only the [O II], [O III] and Hbeta emission lines, and can therefore be used to calculate AGN fractions up to z ~ 0.85 in the optical and z ~ 3.5 in the near-infrared. We anticipate that this diagnostic will facilitate studies of the properties of AGN ionizing radiation fields and the relative optical contribution of star-formation and AGN activity over cosmic time.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Composite Spectra in Merging U/LIRGs Caused by Shocks

We present a key result from our optical integral field spectroscopic survey of 27 nearby ultra luminous and luminous infrared galaxies (U/LIRGs) from the Great Observatory All-sky LIRG Survey. Using spatially resolved multi-component emission line fitting to trace the emission line ratios and velocity dispersion of the ionized gas, we quantify for the first time the widespread shock ionization in gas-rich merging U/LIRGs. Our results show a fractional contribution to the total observed \Ha flux from radiative shocks increasing from a few percent during early merger stages to upwards of 60% of the observed optical emission line flux in late stage mergers. We compare our resolved spectroscopy to nuclear spectra and find that 3/4 of the galaxies in our sample that would be classified as "composite" based on optical spectroscopy are primarily characterized by a combination of star formation and merger-driven shocks. Our results have important implications for the interpretation of "composite" rest-frame optical spectra of U/LIRGs as starburst+AGN, as the shock emission combined with star formation can mimic "composite" optical spectra in the absence of any contribution from an AGN.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Galaxy Interactions in Compact Groups I : The Galactic Winds of HCG16

Using the WiFeS integral field spectrograph, we have undertaken a series of observations of star-forming galaxies in Compact Groups. In this first paper dedicated to the project, we present the analysis of the spiral galaxy NGC838, a member of the Hickson Compact Group 16, and of its galactic wind. Our observations reveal that the wind forms an asymmetric, bipolar, rotating structure, powered by a nuclear starburst. Emission line ratio diagnostics indicate that photoionization is the dominant excitation mechanism at the base of the wind. Mixing from slow shocks (up to 20%) increases further out along the outflow axis. The asymmetry of the wind is most likely caused by one of the two lobes of the wind bubble bursting out of its HI envelope, as indicated by line ratios and radial velocity maps. The characteristics of this galactic wind suggest that it is caught early (a few Myr) in the wind evolution sequence. The wind is also quite different to the galactic wind in the partner galaxy NGC839 which contains a symmetric, shock-excited wind. Assuming that both galaxies have similar interaction histories, the two different winds must be a consequence of the intrinsic properties of NGC838 and NGC839 and their starbursts.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Measuring nebular temperatures: the effect of new collision strengths with equilibrium and kappa-distributed electron energies

In this paper we develop tools for observers to use when analysing nebular spectra for temperatures and metallicities, with two goals: to present a new, simple method to calculate equilibrium electron temperatures for collisionally excited line flux ratios, using the latest atomic data; and to adapt current methods to include the effects of possible non-equilibrium 'κ' electron energy distributions. Adopting recent collision strength data for [O iii], [S iii], [O ii], [S ii], and [N ii], we find that existing methods based on older atomic data seriously overestimate the electron temperatures, even when considering purely Maxwellian statistics. If κ distributions exist in H ii regions and planetary nebulae as they do in solar system plasmas, it is important to investigate the observational consequences. This paper continues our previous work on the κ distribution (Nicholls et al. 2012). We present simple formulaic methods that allow observers to (a) measure equilibrium electron temperatures and atomic abundances using the latest atomic data, and (b) to apply simple corrections to existing equilibrium analysis techniques to allow for possible non-equilibrium effects. These tools should lead to better consistency in temperature and abundance measurements, and a clearer understanding of the physics of H ii regions and planetary nebulae.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

New Strong Line Abundance Diagnostics for \HII Regions: Effects of $κ$-Distributed Electron Energies and New Atomic Data

Recently, \citet{Nicholls12}, inspired by \emph{in situ} observations of solar system astrophysical plasmas, suggested that the electrons in \HII regions are characterised by a $κ$-distribution of electron energies rather than by a simple Maxwell-Boltzmann distribution. Here we have collected together the new atomic data within a modified photoionisation code to explore the effects of both the new atomic data and the $κ$-distribution on the strong-line techniques used to determine chemical abundances in \HII regions. By comparing the recombination temperatures ($T_{\rm rec}$) with the forbidden line temperatures ($T_{\rm FL}$) we conclude that $ κ\sim 20$. While representing only a mild deviation from equilibrium, this is sufficient to strongly influence abundances determined using methods which depend on measurements of the electron temperature from forbidden lines. We present a number of new emission line ratio diagnostics which cleanly separate the two parameters determining the optical spectrum of \HII regions - the ionisation parameter $q$ or $\cal{U}$ and the chemical abundance; 12+log(O/H). An automated code to extract these parameters is presented. Using the homogeneous dataset from \citet{vanZee98}, we find self-consistent results between all these different diagnostics. The systematic errors between different line ratio diagnostics are much smaller than was found in the earlier strong line work. Overall the effect of the $κ$-distribution on the strong line abundances derived solely on the basis of theoretical models is rather small.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Spaxel Analysis: Probing the Physics of Star Formation in Ultraluminous Infrared Galaxies

This paper presents a detailed spectral pixel (spaxel) analysis of the ten Luminous Infrared Galaxies (LIRGs) previously observed with the Wide Field Spectrograph (WiFeS), an integral field spectrograph mounted on the ANU 2.3m telescope, and for which an abundance gradient analysis has already been presented by Rich et al. (2012). Here we use the strong emission line analysis techniques developed by Dopita et al. (2013) to measure the ionisation parameter and the oxygen abundance in each spaxel. In addition, we use the observed H$α$ flux to determine the surface rate of star formation ($M_{\odot}$ yr$^{-1}$ kpc$^{-2}$) and use the [\ion{S}{2}] $λ\lambda6717/6731$ ratio to estimate the local pressure in the ionised plasma. We discuss the correlations discovered between these physical quantities, and use them to infer aspects of the physics of star formation in these extreme star forming environments. In particular, we find a correlation between the star formation rate and the inferred ionisation parameter. We examine the possible reasons for this correlation, and determine that the most likely explanation is that the more active star forming regions have a different distribution of molecular gas which favour higher ionisation parameters in the ionised plasma.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

The Cosmic BPT Diagram: Confronting Theory with Observations

We compare a large sample of galaxies between 0.5<z<2.6 with theoretical predictions for how the optical diagnostic line ratios in galaxy ensembles change as a function of cosmic time. We show that star forming galaxies at high redshift (z>1.5) are consistent with a model in which the ISM conditions are more extreme at high redshift than seen in the global spectra of local galaxies. We speculate that global spectra of our high redshift galaxies may be dominated by HII regions similar to the extreme clumpy, dense star-forming complexes in the Antennae and M82. The transition to local-type conditions occurs between 0.8<z<1.5. We conclude that classification schemes developed for local samples should not be applied at high redshift (z > 1.5). We use our theoretical models to derive a new redshift-dependent classification line that utilizes the standard optical diagnostic line ratios [OIII]/H-beta and [NII]/H-alpha. Our new line can be used to separate star-forming galaxies from AGN between z=0 to z~3.5. We anticipate that our redshift-dependent optical classification line will be useful for future large surveys with near-infrared multi-object spectrographs. We apply our classification line to a sample of gravitationally lensed galaxies at z~2.5. Although limited by small numbers, we show that our classification line is consistent with the position of AGN that have been independently confirmed via other methods.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

The Observed Relation between Stellar Mass, Dust Extinction and Star Formation Rate in Local Galaxies

In this study we investigate the relation between stellar mass, dust extinction and star formation rate (SFR) using ~150,000 star-forming galaxies from the SDSS DR7. We show that the relation between dust extinction and SFR changes with stellar mass. For galaxies at the same stellar mass dust extinction is anti-correlated with the SFR at stellar masses <10^10 M_solar. There is a sharp transition in the relation at a stellar mass of 10^10 M_solar. At larger stellar masses dust extinction is positively correlated with the SFR for galaxies at the same stellar mass. The observed relation between stellar mass, dust extinction and SFR presented in this study helps to confirm similar trends observed in the relation between stellar mass, metallicity and SFR. The relation reported in this study provides important new constraints on the physical processes governing the chemical evolution of galaxies. The correlation between SFR and dust extinction for galaxies with stellar masses >10^10 M_solar is shown to extend to the population of quiescent galaxies suggesting that the physical processes responsible for the observed relation between stellar mass, dust extinction and SFR may be related to the processes leading to the shut down of star formation in galaxies.

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

A Census of Oxygen in Star-Forming Galaxies: An Empirical Model Linking Metallicities, Star Formation Rates and Outflows

In this contribution we present the first census of oxygen in star-forming galaxies in the local universe. We examine three samples of galaxies with metallicities and star formation rates at z = 0.07, 0.8 and 2.26, including the SDSS and DEEP2 surveys. We infer the total mass of oxygen produced and mass of oxygen found in the gas-phase from our local SDSS sample. The star formation history is determined by requiring that galaxies evolve along the relation between stellar mass and star formation rate observed in our three samples. We show that the observed relation between stellar mass and star formation rate for our three samples is consistent with other samples in the literature. The mass-metallicity relation is well established for our three samples and from this we empirically determine the chemical evolution of star-forming galaxies. Thus, we are able to simultaneously constrain the star formation rates and metallicities of galaxies over cosmic time allowing us to estimate the mass of oxygen locked up in stars. Combining this work with independent measurements reported in the literature we conclude that the loss of oxygen from the interstellar medium of local star-forming galaxies is likely to be a ubiquitous process with the oxygen mass loss scaling (almost) linearly with stellar mass. We estimate the total baryonic mass loss and argue that only a small fraction of the baryons inferred from cosmological observations accrete onto galaxies.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Galaxy-Wide Shocks in Late-Merger Stage Luminous Infrared Galaxies

We present an integral field spectroscopic study of two nearby Luminous Infrared Galaxies (LIRGs) that exhibit evidence of widespread shock excitation induced by ongoing merger activity, IC 1623 and NGC 3256. We show the importance of carefully separating excitation due to shocks vs. excitation by HII regions and the usefulness of IFU data in interpreting the complex processes in LIRGs. Our analysis focuses primarily on the emission line gas which is extensive in both systems and is a result of the abundant ongoing star formation as well as widespread LINER-like excitation from shocks. We use emission-line ratio maps, line kinematics, line-ratio diagnostics and new models as methods for distinguishing and analyzing shocked gas in these systems. We discuss how our results inform the merger sequence associated with local U/LIRGs and the impact that widespread shock excitation has on the interpretation of emission-line spectra and derived quantities of both local and high-redshift galaxies.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

A High-Metallicity Host Environment for the Long-Duration GRB 020819

We present spectroscopic observations of the host galaxy and explosion site of the long-duration gamma-ray burst GRB 020819. We determine a metallicity for this host environment of log(O/H) + 12 = 9.0 +/- 0.1, by far the highest metallicity determined for a long-duration GRB to date. We compare the metallicity and other properties of the GRB 020819 host environment to existing observations of long-duration GRB host galaxies, and consider the implications that this metallicity has for our understanding of long-duration GRB progenitor scenarios. We also consider how this unusually high metallicity may relate to the status of GRB 020819 as a "dark" burst, with no detected optical afterglow.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Galaxy Mergers and the Mass-Metallicity Relation: Evidence for Nuclear Metal Dilution and Flattened Gradients from Numerical Simulations

Recent results comparing interacting galaxies to the mass-metallicity relation show that their nuclear oxygen abundances are unexpectedly low. We present analysis of N-body/SPH numerical simulations of equal-mass mergers that confirm the hypothesis that these underabundances are accounted for by radial inflow of low-metallicity gas from the outskirts of the two merging galaxies. The underabundances arise between first and second pericenter, and the simulated abundance dilution is in good agreement with observations. The simulations further predict that the radial metallicity gradients of the disk galaxies flatten shortly after first passage, due to radial mixing of gas. These predictions will be tested by future observations of the radial metallicity distributions in interacting galaxies.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Gas-Phase Oxygen Gradients in Strongly Interacting Galaxies: I. Early-Stage Interactions

A consensus is emerging that interacting galaxies show depressed nuclear gas metallicities compared to isolated star-forming galaxies. Simulations suggest that this nuclear underabundance is caused by interaction-induced inflow of metal-poor gas, and that this inflow concurrently flattens the radial metallicity gradients in strongly interacting galaxies. We present metallicities of over 300 HII regions in a sample of 16 spirals that are members of strongly interacting galaxy pairs with mass ratio near unity. The deprojected radial gradients in these galaxies are about half of those in a control sample of isolated, late-type spirals. Detailed comparison of the gradients with simulations show remarkable agreement in gradient distributions, the relationship between gradients and nuclear underabundances, and the shape of profile deviations from a straight line. Taken together, this evidence conclusively demonstrates that strongly interacting galaxies at the present day undergo nuclear metal dilution due to gas inflow, as well as significant flattening of their gas-phase metallicity gradients, and that current simulations can robustly reproduce this behavior at a statistical level.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Metallicity Gradients and Gas Flows in Galaxy Pairs

We present the first systematic investigation into the metallicity gradients in galaxy close pairs. We determine the metallicity gradients for 8 galaxies in close pairs using HII region metallicities obtained with high signal-to-noise multi-slit observations with the Keck LRIS Spectrograph. We show that the metallicity gradients in close pairs are significantly shallower than gradients in isolated spiral galaxies such as the Milky Way, M83, and M101. These observations provide the first solid evidence that metallicity gradients in interacting galaxies are systematically different from metallicity gradients in isolated spiral galaxies. Our results suggest that there is a strong relationship between metallicity gradients and the gas dynamics in galaxy interactions and mergers.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Modeling IR Spectral Energy Distributions: A Pilot Study of Starburst Parameters and Silicate Absorption Curves for Some GOALS Galaxies

This paper describes a pilot study into the spectral energy distribution (SED) fitting and the derivation of physical parameters for 19 galaxies observed as part of the Great Observatories All-sky LIRG Survey (GOALS) survey as observed with the \emph{Spitzer Space Telescope}. For this we have used the pan-spectral fitting tools developed in a series of papers by Dopita and his co-workers. We show that the standard Lee and Draine `astronomical silicate' model cannot provide a good fit to the silicate absorption features as observed in the heavily dust-extinguished ($A_{\rm V} \sim 50$mag.) starbursts. We have derived an empirical fit to the `starburst silicate' absorption in these objects. This absorption curve is consistent with the silicate grains being systematically larger in starburst environments than in the local Galactic interstellar medium. We demonstrate the sensitivity of the SED fitting to each of the fitted parameters, and derive these parameters for those galaxies which do not have an embedded AGN. This technique is simple and provides reasonably robust and uniform parameters for the starburst, especially as far as the star formation rate, population of old stars, compactness of the starburst region and total foreground extinction are concerned. However, the chemical abundances and the optical extinction cannot be reliably determined by this analysis, and optical SEDs will also be required to provide a complete characterization of the starburst region and of the surrounding galaxy.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

No Correlation Between Host Galaxy Metallicity and Gamma-Ray Energy Release for Long-Duration Gamma-Ray Bursts

We compare the redshifts, host galaxy metallicities, and isotropic (E_gamma,iso) and beaming-corrected (E_gamma) gamma-ray energy release of 16 long-duration gamma-ray bursts (LGRBs) at z < 1. From this comparison, we find no statistically significant correlation between host metallicity and redshift, E_gamma,iso, or E_gamma. These results are at odds with previous theoretical and observational predictions of an inverse correlation between gamma-ray energy release and host metallicity, as well as the standard predictions of metallicity-driven wind effects in stellar evolutionary models. We consider the implications that these results have for LGRB progenitor scenarios, and discuss our current understanding of the role that metallicity plays in the production of LGRBs.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

The Host Galaxies of Gamma-Ray Bursts II: A Mass-Metallicity Relation for Long-Duration Gamma-Ray Burst Host Galaxies

We present a statistically robust mass-metallicity (M-Z) relation for long-duration gamma-ray burst (LGRB) host galaxies at z < 1. By comparing the LGRB host M-Z relation to samples representative of the general star-forming galaxy population, we conclude that LGRBs occur in host galaxies with lower metallicities than the general population, and that this trend extends to z ~ 1, with an average offset of -0.42 +/- 0.18 from the M-Z relation for star-forming galaxies. Our sample in this work includes new spectroscopic data for 6 LGRB host galaxies obtained at the Keck and Magellan telescopes, as well as 2 new host galaxies from the literature. Combined with data from our previous work, this yields a total sample of 6 LGRB host galaxies at z < 0.3 and 10 host galaxies at 0.3 < z < 1. We have determined a number of interstellar medium properties for our host galaxies using optical emission-line diagnostics, including metallicity, ionization parameter, young stellar population age, and star formation rate. Across our full sample of 16 LGRB hosts we find an average metallicity of log(O/H) + 12 = 8.4 +/- 0.3. Notably, we also measure a comparatively high metallicity of log(O/H) + 12 = 8.83 +/- 0.1 for the z = 0.296 host galaxy of GRB 050826. We also determine stellar masses (M*) for our LGRB host galaxy sample, finding a mean stellar mass of log(M*/Msun) = 9.25 (+0.19,-0.23).

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

GRB 090426: The Environment of a Rest-Frame 0.35-second Gamma-Ray Burst at Redshift z=2.609

We present the discovery of an absorption-line redshift of z = 2.609 for GRB 090426, establishing the first firm lower limit to a redshift for a gamma-ray burst with an observed duration of <2 s. With a rest-frame burst duration of T_90z = 0.35 s and a detailed examination of the peak energy of the event, we suggest that this is likely (at >90% confidence) a member of the short/hard phenomenological class of GRBs. From analysis of the optical-afterglow spectrum we find that the burst originated along a very low HI column density sightline, with N_HI < 3.2 x 10^19 cm^-2. Our GRB 090426 afterglow spectrum also appears to have weaker low-ionisation absorption (Si II, C II) than ~95% of previous afterglow spectra. Finally, we also report the discovery of a blue, very luminous, star-forming putative host galaxy (~2 L*) at a small angular offset from the location of the optical afterglow. We consider the implications of this unique GRB in the context of burst duration classification and our understanding of GRB progenitor scenarios.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

The Host Galaxies of Gamma-Ray Bursts I: ISM Properties of Ten Nearby Long-Duration GRB Hosts

We present the first observations from a large-scale survey of nearby (z < 1) long-duration gamma-ray burst (LGRB) host galaxies, which consist of eight rest-frame optical spectra obtained at Keck and Magellan. Along with two host galaxy observations from the literature, we use optical emission line diagnostics to determine metallicities, ionization parameters, young stellar population ages, and star formation rates. We compare the LGRB host environments to a variety of local and intermediate-redshift galaxy populations, as well as the newest grid of stellar population synthesis and photoionization models generated with the Starburst99/Mappings codes. With these comparisons we investigate whether the GRB host galaxies are consistent with the properties of the general galaxy population, and therefore whether they may be used as reliable tracers of star formation. We find that LGRB host galaxies generally have low-metallicity ISM environments out to z ~ 1. The ISM properties of our GRB hosts, including metallicity, ionization parameter, and young stellar population age, are significantly different from the general galaxy population, host galaxies of nearby broad-lined Type Ic supernovae, and nearby metal-poor galaxies.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Theoretical Modeling of Star-Forming Galaxies I. Emission Line Diagnostic Grids for Local and Low-Metallicity Galaxies

We use the newest generation of the Starburst99/Mappings code to generate an extensive suite of models to facilitate detailed studies of star-forming galaxies and their ISM properties, particularly at low metallicities. The new models used include a rigorous treatment of metal opacities in the population synthesis modeling and more detailed dust physics in the photoionization code. These models span a wide range of physical parameters including metallicity, ionization parameter, and the adoption of both a instantaneous burst and continuous star formation history. We examine the agreement between our models and local (z < 0.1) star-forming galaxy populations from several large datasets, including the Sloan Digital Sky Survey, the Nearby Field Galaxy Survey, and samples of blue compact galaxies and metal-poor galaxies. We find that models adopting a continuous star formation history reproduce the metallicity-sensitive line ratios observed in the local population of star-forming galaxies, including the low-metallicity sample. However, we find that the current codes generate an insufficiently hard ionizing radiation field, leading to deficiencies in the [SII] fluxes produced by the models. We consider the advantages and short-comings of this suite of models, and discuss future work and improvements that can be applied to the modeling of star-forming galaxies.

0 citations0 reviewsNo code linkedOpen access
preprint2008arXiv

Metallicity Calibrations and the Mass-Metallicity Relation for Star-Forming Galaxies

(Abridged) We investigate the effect of metallicity calibrations, AGN classification, and aperture covering fraction on the local mass-metallicity (MZ) relation using 27,730 star-forming galaxies from the Sloan Digital Sky Survey (SDSS) Data Release 4. We analyse the SDSS MZ relation with 10 metallicity calibrations, including theoretical and empirical methods. We show that the choice of metallicity calibration has a significant effect on the shape and y-intercept(12+log(O/H)) of the MZ relation. The absolute metallicity scale (y-int) varies up to 0.7 dex, depending on the calibration used, and the change in shape is substantial. These results indicate that it is critical to use the same metallicity calibration when comparing different luminosity-metallicity or mass-metallicity relations. We present new metallicity conversions that allow metallicities that have been derived using different strong-line calibrations to be converted to the same base calibration. These conversions facilitate comparisons between different samples, particularly comparisons between galaxies at different redshifts for which different suites of emission-lines are available.Our new conversions successfully remove the large 0.7 dex discrepancies between the metallicity calibrations, and we reach agreement in the MZ relation to within 0.03 dex on average. We find that different AGN classification methods have negligible effect on the SDSS MZ relation. We compare the SDSS MZ relation with nuclear and global relations from the Nearby Field Galaxy Survey (NFGS). The turn over of the MZ relation depends on aperture covering fraction. We find that a lower redshift limit of z<0.04 is insufficient for avoiding aperture effects in fiber spectra of the highest stellar mass galaxies.

0 citations0 reviewsNo code linkedOpen access
preprint2006arXiv

Modelling the Pan-Spectral Energy Distribution of Starburst Galaxies: III. Emission Line Diagnostics of Ensembles of Evolving HII Regions

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

0 citations0 reviewsNo code linkedOpen access