Source author record

Michael A. Dopita

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

ResearcherUnclaimed source record

Catalog footprint

What is connected

56works
6topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

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

Building this map preview

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

Published work

56 published item(s)

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.

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.

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.

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.

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.

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.

preprint2016arXiv

Forbidden Iron Lines and Dust Destruction in Supernova Remnant Shocks: The Case of N49 in the Large Magellanic Cloud

We present results of a complete integral field survey of the bright SNR N49 in the LMC, obtained with the WiFeS instrument mounted on the ANU 2.3m telescope at Siding Spring Observatory. From theoretical shock modelling with the new MAPPINGS 5.1 code we have, for the first time, subjected the optical Fe emission line spectrum of a supernova remnant to a detailed abundance and dynamical analysis covering 8 separate stages of ionisation. This allows us to derive the dust depletion factors as a function of ionisation stage. We have shown that there is substantial (30% - 90%) destruction of Fe-bearing dust grains in these fast shocks ($v_s \sim 250$ km/s), and we have confirmed that the dominant dust destruction is through the non-thermal sputtering and grain-grain collision mechanisms developed in a number of theoretical works.

preprint2016arXiv

Problems for the WELS classification of planetary nebulae central stars: Self-consistent nebular modelling of four candidates

We present integral field unit (IFU) spectroscopy and self-consistent photoionisation modelling for a sample of four southern Galactic planetary nebulae (PNe) with supposed weak emission-line (WEL) central stars. The Wide Field Spectrograph (WiFeS) on the ANU 2.3 m telescope has been used to provide IFU spectroscopy for NGC 3211, NGC 5979, My 60, and M 4-2 covering the spectral range of 3400-7000Å. All objects are high excitation non-Type I PNe, with strong He II emission, strong [Ne V] emission, and weak low-excitation lines. They all appear to be predominantly optically-thin nebulae excited by central stars with $T_{\rm eff} > 10^5$K. Three PNe of the sample have central stars which have been previously classified as weak emission-line stars (WELS), and the fourth also shows the characteristic recombination lines of a WELS. However, the spatially-resolved spectroscopy shows that rather than arising in the central star, the C IV and N III recombination line emission is distributed in the nebula, and in some cases concentrated in discrete nebular knots. This may suggest that the WELS classification is spurious, and that, rather, these lines arise from (possibly chemically enriched) pockets of nebular gas. Indeed, from careful background subtraction we were able to identify three of the sample as being hydrogen rich O(H)-Type. We have constructed fully self-consistent photoionization models for each object. This allows us to independently determine the chemical abundances in the nebulae, to provide new model-dependent distance estimates, and to place the central stars on the H-R diagram. All four PNe have similar initial mass ($1.5 < M/M_{\odot} <2.0$) and are at a similar evolutionary stage.

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

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.

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.

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.

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.

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.

preprint2015arXiv

Probing the Physics of Narrow Line Regions in Active Galaxies II: The Siding Spring Southern Seyfert Spectroscopic Snapshot Survey (S7)

Here we describe the \emph{Siding Spring Southern Seyfert Spectroscopic Snapshot Survey} (S7) and present results on 64 galaxies drawn from the first data release. The S7 uses the Wide Field Spectrograph (WiFeS) mounted on the ANU 2.3m telescope located at the Siding Spring Observatory to deliver an integral field of $38\times25$~ arcsec at a spectral resolution of $R=7000$ in the red ($530-710$nm), and $R=3000$ in the blue ($340-560$nm). {From these data cubes we have extracted the Narrow Line Region (NLR) spectra from a 4 arc sec aperture centred on the nucleus. We also determine the H$β$ and [OIII]~$λ$5007 fluxes in the narrow lines, the nuclear reddening, the reddening-corrected relative intensities of the observed emission lines, and the H$β$ and \lOIII\ luminosities {determined from spectra for which the stellar continuum has been removed.} We present a set of images of the galaxies in [OIII]~$λ$5007, [NII]~$λ$6584 and H$α$ which serve to delineate the spatial extent of the extended narrow line region (ENLR) and {\bf also to} reveal the structure and morphology of the surrounding \HII\ regions. Finally, we provide a preliminary discussion of those Seyfert~1 and Seyfert~2 galaxies which display coronal emission lines in order to explore the origin of these lines.

preprint2015arXiv

Probing the Physics of Narrow Line Regions in Active Galaxies III: Accretion and Cocoon Shocks in the LINER NGC1052

We present Wide Field Spectrograph (WiFeS) integral field spectroscopy and HST FOS spectroscopy for the LINER galaxy NGC 1052. We infer the presence of a turbulent accretion flow forming a small-scale accretion disk. We find a large-scale outflow and ionisation cone along the minor axis of the galaxy. Part of this outflow region is photoionised by the AGN, and shares properties with the ENLR of Seyfert galaxies, but the inner ($R \lesssim 1.0$~arcsec) accretion disk and the region around the radio jet appear shock excited. The emission line properties can be modelled by a "double shock" model in which the accretion flow first passes through an accretion shock in the presence of a hard X-ray radiation, and the accretion disk is then processed through a cocoon shock driven by the overpressure of the radio jets. This model explains the observation of two distinct densities ($\sim10^4$ and $\sim10^6$ cm$^{-3}$), and provides a good fit to the observed emission line spectrum. We derive estimates for the velocities of the two shock components and their mixing fractions, the black hole mass, the accretion rate needed to sustain the LINER emission and derive an estimate for the jet power. Our emission line model is remarkably robust against variation of input parameters, and so offers a generic explanation for the excitation of LINER galaxies, including those of spiral type such as NGC 3031 (M81).

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.

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.

preprint2014arXiv

A Newly Recognized Very Young Supernova Remnant in M83

As part of a spectroscopic survey of supernova remnant candidates in M83 using the Gemini-South telescope and GMOS, we have discovered one object whose spectrum shows very broad lines at H$α$, [O~I] 6300,6363, and [O~III] 4959,5007, similar to those from other objects classified as `late time supernovae.' Although six historical supernovae have been observed in M83 since 1923, none were seen at the location of this object. Hubble Space Telescope Wide Field Camera 3 images show a nearly unresolved emission source, while Chandra and ATCA data reveal a bright X-ray source and nonthermal radio source at the position. Objects in other galaxies showing similar spectra are only decades post-supernova, which raises the possibility that the supernova that created this object occurred during the last century but was missed. Using photometry of nearby stars from the HST data, we suggest the precursor was at least 17 $\rm M_{sun}$, and the presence of broad H$α$ in the spectrum makes a type II supernova likely. The supernova must predate the 1983 VLA radio detection of the object. We suggest examination of archival images of M83 to search for evidence of the supernova event that gave rise to this object, and thus provide a precise age.

preprint2014arXiv

An Expanded HST/WFC3 Survey of M83: Project Overview and Targeted Supernova Remnant Search

We present an optical/NIR imaging survey of the face-on spiral galaxy M83, using data from the Hubble Space Telescope Wide Field Camera 3 (WFC3). Seven fields are used to cover a large fraction of the inner disk, with observations in nine broadband and narrowband filters. In conjunction with a deep Chandra survey and other new radio and optical ground-based work, these data enable a broad range of science projects to be pursued. We provide an overview of the WFC3 data and processing and then delve into one topic, the population of young supernova remnants. We used a search method targeted toward soft X-ray sources to identify 26 new supernova remnants. Many compact emission nebulae detected in [Fe II] 1.644 micron align with known remnants and this diagnostic has also been used to identify many new remnants, some of which are hard to find with optical images. We include 37 previously identified supernova remnants that the data reveal to be <0.5'' in angular size and thus are difficult to characterize from ground-based data. The emission line ratios seen in most of these objects are consistent with shocks in dense interstellar material rather than showing evidence of ejecta. We suggest that the overall high elemental abundances in combination with high interstellar medium pressures in M83 are responsible for this result. Future papers will expand on different aspects of the these data including a more comprehensive analysis of the overall supernova remnant population.

preprint2014arXiv

Choirs HI Galaxy Groups: The metallicity of dwarf galaxies

We present a recalibration of the luminosity-metallicity relation for gas-rich, star-forming dwarfs to magnitudes as faint as M$_R$ ~ -13. We use the Dopita et al. (2013) metallicity calibrations to calibrate the relation for all of the data in this analysis. In metallicity-luminosity space we find two sub-populations within a sample of high-confidence SDSS DR8 star-forming galaxies; 52% are metal-rich giants and 48% are metal-medium galaxies. Metal-rich dwarfs classified as tidal dwarf galaxy (TDG) candidates in the literature are typically of metallicity 12 + log(O/H) = 8.70 $\pm$ 0.05, while SDSS dwarfs fainter than M$_R$ = -16 have a mean metallicity of 12 + log(O/H) = 8.28 $\pm$ 0.10, regardless of their luminosity, indicating that there is an approximate floor to the metallicity of low luminosity galaxies. Our hydrodynamical simulations predict that TDGs should have metallicities elevated above the normal luminosity-metallicity relation. Metallicity can therefore be a useful diagnostic for identifying TDG candidate populations in the absence of tidal tails. At magnitudes brighter than M$_R$ ~ -16 our sample of 53 star-forming galaxies in 9 HI gas-rich groups is consistent with the normal relation defined by the SDSS sample. At fainter magnitudes there is an increase in dispersion in metallicity of our sample, suggestive of a wide range of HI content and environment. In our sample we identify three (16% of dwarfs) strong TDG candidates (12 + log(O/H) > 8.6), and four (21%) very metal poor dwarfs (12 + log(O/H) < 8.0), which are likely gas-rich dwarfs with recently ignited star formation.

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.

preprint2014arXiv

IZI: Inferring the Gas Phase Metallicity (Z) and Ionization Parameter (q) of Ionized Nebulae using Bayesian Statistics

We present a new method for inferring the metallicity (Z) and ionization parameter (q) of HII regions and star-forming galaxies using strong nebular emission lines (SEL). We use Bayesian inference to derive the joint and marginalized posterior probability density functions for Z and q given a set of observed line fluxes and an input photo-ionization model. Our approach allows the use of arbitrary sets of SELs and the inclusion of flux upper limits. The method provides a self-consistent way of determining the physical conditions of ionized nebulae that is not tied to the arbitrary choice of a particular SEL diagnostic and uses all the available information. Unlike theoretically calibrated SEL diagnostics the method is flexible and not tied to a particular photo-ionization model. We describe our algorithm, validate it against other methods, and present a tool that implements it called IZI. Using a sample of nearby extra-galactic HII regions we assess the performance of commonly used SEL abundance diagnostics. We also use a sample of 22 local HII regions having both direct and recombination line (RL) oxygen abundance measurements in the literature to study discrepancies in the abundance scale between different methods. We find that oxygen abundances derived through Bayesian inference using currently available photo-ionization models in the literature can be in good (~30%) agreement with RL abundances, although some models perform significantly better than others. We also confirm that abundances measured using the direct method are typically 0.2 dex lower than both RL and photo-ionization model based abundances.

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.

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.

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.

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.

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.

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.

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.

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.

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.

preprint2013arXiv

Nebular metallicities in two isolated Local Void dwarf galaxies

Isolated dwarf galaxies, especially those situated in voids, may provide insight into primordial conditions in the Universe and the physical processes that govern star formation in undisturbed stellar systems. The metallicity of H ii regions in such galaxies is key to investigating this possibility. From the SIGRID sample of isolated dwarf galaxies, we have identified two exceptionally isolated objects, the Local Void galaxy [KK98]246 (ESO 461-G036) and another somewhat larger dwarf irregular on the edge of the Local Void, MCG-01-41-006 (HIPASS J1609- 04). We report our measurements of the nebular metallicities in these objects. The first object has a single low luminosity Hii region, while the second is in a more vigorous star forming phase with several bright H ii regions. We find that the metallicities in both galaxies are typical for galaxies of this size, and do not indicate the presence of any primordial gas, despite (for [KK98]246) the known surrounding large reservoir of neutral Hydrogen.

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.

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.

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.

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.

preprint2012arXiv

Feedback in the cores of clusters Abell 3581, 2A 0335+096, and Sersic 159-03

The cores of massive galaxy clusters, where hot gas is cooling rapidly, appear to undergo cycles of self-regulating energy feedback, in which AGN outbursts in the central galaxies episodically provide sufficient heating to offset much of the gas cooling. We use deep integral-field spectroscopy to study the optical line emission from the extended nebulae of three nearby brightest cluster galaxies and investigate how they are related to the processes of heating and cooling in the cluster cores. Two of these systems, Abell 3581 and Sersic 159-03, appear to be experiencing phases of feedback that are dominated by the activity and output of a central AGN. Abell 3581, shows evidence for significant interaction between the radio outflows and the optical nebula, in addition to accretion flows into the nucleus of the galaxy. X-ray and radio data show that Sersic 159-03 is dominated by the feedback of energy from the central AGN, but the kinematics of the optical nebula are consistent with infall or outflow of material along its bright filaments. The third system, 2A 0335+096, is dominated by mass accretion and cooling, and so we suggest that it is in an accumulation phase of the feedback cycle. The outer nebula forms a disk-like structure, ~14 kpc in radius, that rotates about the central galaxy with a velocity amplitude of ~200 km/s. Overall, our data are consistent with ongoing AGN-driven feedback cycles occurring in these systems.

preprint2012arXiv

Hot Gas, Cold Gas and Sub-Halos in a Lyman-alpha Blob at Redshift 2.38

We present integral field spectroscopy of a Lyman-alpha blob at redshift 2.38, with a spectral resolution three times better than previous published work. As with previous observations, the blob has a chaotic velocity structure, much of which breaks up into multiple components. Our spectroscopy shows, however, that some of these multiple components are extremely narrow: they have velocity widths of less than 100 km/s. Combining these new data with previous observations, we argue that this Lyman-alpha blob resides in a dark-matter halo of around 10^13 solar masses. At the centre of this halo are two compact red massive galaxies. They are surrounded by hot gas, probably a super-wind from merger-induced nuclear starbursts. This hot gas has shut down star formation in the non-nuclear region of these galaxies, leading to their red-and-dead colours. A filament or lump of infalling cold gas is colliding with the hot gas phase and being shocked to high temperatures, while still around 30kpc from the red galaxies. The shock region is self-absorbed in Lyman-alpha but produces C IV emission. Further out still, the cold gas in a number of sub-halos is being lit up, most likely by a combination of tidally triggered star formation, bow-shocks as they plough through the hot halo medium, resonant scattering of Lyman-alpha from the filament collision, and tidal stripping of gas which enhances the Lyman-alpha escape fraction. The observed Lyman-alpha emission from the Blob is dominated by the sum of the emission from these sub-halos. On statistical grounds, we argue that Lyman-alpha blobs are not greatly elongated in shape, and that most are not powered by ionisation or scattering from a central active galactic nucleus or starburst.

preprint2012arXiv

Resolving the electron temperature discrepancies in HII Regions and Planetary Nebulae: kappa-distributed electrons

The measurement of electron temperatures and metallicities in H ii regions and Planetary Nebulae (PNe) has-for several decades-presented a problem: results obtained using different techniques disagree. What it worse, they disagree consistently. There have been numerous attempts to explain these discrepancies, but none has provided a satisfactory solution to the problem. In this paper, we explore the possibility that electrons in H ii regions and PNe depart from a Maxwell-Boltzmann equilibrium energy distribution. We adopt a "kappa-distribution" for the electron energies. Such distributions are widely found in Solar System plasmas, where they can be directly measured. This simple assumption is able to explain the temperature and metallicity discrepancies in H ii regions and PNe arising from the different measurement techniques. We find that the energy distribution does not need to depart dramatically from an equilibrium distribution. From an examination of data from Hii regions and PNe it appears that kappa ~ 10 is sufficient to encompass nearly all objects. We argue that the kappa-distribution offers an important new insight into the physics of gaseous nebulae, both in the Milky Way and elsewhere, and one that promises significantly more accurate estimates of temperature and metallicity in these regions.

preprint2012arXiv

The Resolved Stellar Population in 50 Regions of M83 from HST/WFC3 Early Release Science Observations

We present a multi-wavelength photometric study of ~15,000 resolved stars in the nearby spiral galaxy M83 (NGC5236, D=4.61Mpc) based on Hubble Space Telescope Wide Field Camera 3 observations using four filters: F336W, F438W, F555W, and F814W. We select 50 regions (an average size of 260 pc by 280 pc) in the spiral arm and inter-arm areas of M83, and determine the age distribution of the luminous stellar populations in each region. This is accomplished by correcting for extinction towards each individual star by comparing its colors with predictions from stellar isochrones. We compare the resulting luminosity weighted mean ages of the luminous stars in the 50 regions with those determined from several independent methods, including the number ratio of red-to-blue supergiants, morphological appearance of the regions, surface brightness fluctuations, and the ages of clusters in the regions. We find reasonably good agreement between these methods. We also find that young stars are much more likely to be found in concentrated aggregates along spiral arms, while older stars are more dispersed. These results are consistent with the scenario that star formation is associated with the spiral arms, and stars form primarily in star clusters and then disperse on short timescales to form the field population. The locations of Wolf-Rayet stars are found to correlate with the positions of many of the youngest regions, providing additional support for our ability to accurately estimate ages. We address the effects of spatial resolution on the measured colors, magnitudes, and age estimates. While individual stars can occasionally show measurable differences in the colors and magnitudes, the age estimates for entire regions are only slightly affected.

preprint2011arXiv

Detection of brown dwarf-like objects in the core of NGC3603

We use near-infrared data obtained with the Wide Field Camera 3 (WFC3) on the Hubble Space Telescope to identify objects having the colors of brown dwarfs (BDs) in the field of the massive galactic cluster NGC 3603. These are identified through use of a combination of narrow and medium band filters spanning the J and H bands, and which are particularly sensitive to the presence of the 1.3-1.5μm H2O molecular band - unique to BDs. We provide a calibration of the relationship between effective temperature and color for both field stars and for BDs. This photometric method provides effective temperatures for BDs to an accuracy of {\pm}350K relative to spectroscopic techniques. This accuracy is shown to be not significantly affected by either stellar surface gravity or uncertainties in the interstellar extinction. We identify nine objects having effective temperature between 1700 and 2200 K, typical of BDs, observed J-band magnitudes in the range 19.5-21.5, and that are strongly clustered towards the luminous core of NGC 3603. However, if these are located at the distance of the cluster, they are far too luminous to be normal BDs. We argue that it is unlikely that these objects are either artifacts of our dataset, normal field BDs/M-type giants or extra-galactic contaminants and, therefore, might represent a new class of stars having the effective temperatures of BDs but with luminosities of more massive stars. We explore the interesting scenario in which these objects would be normal stars that have recently tidally ingested a Hot Jupiter, the remnants of which are providing a short-lived extended photosphere to the central star. In this case, we would expect them to show the signature of fast rotation.

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.

preprint2011arXiv

Large-scale shock-ionized and photo-ionized gas in M83: the impact of star formation

We investigate the ionization structure of the nebular gas in M83 using the line diagnostic diagram, [O III](5007 \degA)/Hβ vs. [S II](6716 °A+6731 °A)/Hα with the newly available narrowband images from the Wide Field Camera 3 (WFC3) of the Hubble Space Telescope (HST). We produce the diagnostic diagram on a pixel-by-pixel (0.2" x 0.2") basis and compare it with several photo- and shock-ionization models. For the photo-ionized gas, we observe a gradual increase of the log([O III]/Hβ) ratios from the center to the spiral arm, consistent with the metallicity gradient, as the H II regions go from super solar abundance to roughly solar abundance from the center out. Using the diagnostic diagram, we separate the photo-ionized from the shock-ionized component of the gas. We find that the shock-ionized Hα emission ranges from ~2% to about 15-33% of the total, depending on the separation criteria used. An interesting feature in the diagnostic diagram is an horizontal distribution around log([O III]/Hβ) ~ 0. This feature is well fit by a shock-ionization model with 2.0 Z\odot metallicity and shock velocities in the range of 250 km/s to 350 km/s. A low velocity shock component, < 200 km/s, is also detected, and is spatially located at the boundary between the outer ring and the spiral arm. The low velocity shock component can be due to : 1) supernova remnants located nearby, 2) dynamical interaction between the outer ring and the spiral arm, 3) abnormal line ratios from extreme local dust extinction. The current data do not enable us to distinguish among those three possible interpretations. Our main conclusion is that, even at the HST resolution, the shocked gas represents a small fraction of the total ionized gas emission at less than 33% of the total. However, it accounts for virtually all of the mechanical energy produced by the central starburst in M83.

preprint2011arXiv

Re-ionizing the Universe without Stars

Recent observations show that the measured rates of star formation in the early universe are insufficient to produce re-ionization, and therefore, another source of ionizing photons is required. In this \emph{Letter}, we examine the possibility that these can be supplied by the fast accretion shocks formed around the cores of the most massive haloes ($10.5< \log M/M_{\odot} < 12$) on spatial scales of order 1 kpc. We model the detailed physics of these fast accretion shocks, and apply these to a simple 1-D spherical hydrodynamic accretion model for baryonic infall in dark matter halos with an Einasto density distribution. The escape of UV photons from these halos is delayed by the time taken to reach the critical accretion shock velocity for escape of UV photons; 220 km s$^{-1}$, and by the time it takes for these photons to ionize the surrounding baryonic matter in the accretion flow. Assuming that in the universe at large the baryonic matter tracks the dark matter, we can estimate the epoch of re-ionization in the case that accretion shocks act alone as the source of UV photons. We find that 50% of the volume (and 5-8% of the mass) of the universe can be ionized by $z \sim 7-8$. The UV production rate has an uncertainty of a factor of about 5 due to uncertainties in the cosmological parameters controlling the development of large scale structure. Because our mechanism is a steeply rising function of decreasing redshift, this uncertainty translates to a re-ionization redshift uncertainty of less than $\pm0.5$. We also find that, even without including the UV photon production of stars, re-ionization is essentially complete by $z \sim 5.8$. Thus, fast accretion shocks can provide an important additional source of ionizing photons in the early universe.

preprint2011arXiv

Star formation in 30 Doradus

Using observations obtained with the Wide Field Camera 3 (WFC3) on board the Hubble Space Telescope (HST), we have studied the properties of the stellar populations in the central regions of 30 Dor, in the Large Magellanic Cloud. The observations clearly reveal the presence of considerable differential extinction across the field. We characterise and quantify this effect using young massive main sequence stars to derive a statistical reddening correction for most objects in the field. We then search for pre-main sequence (PMS) stars by looking for objects with a strong (> 4 sigma) Halpha excess emission and find about 1150 of them over the entire field. Comparison of their location in the Hertzsprung-Russell diagram with theoretical PMS evolutionary tracks for the appropriate metallicity reveals that about one third of these objects are younger than ~4Myr, compatible with the age of the massive stars in the central ionising cluster R136, whereas the rest have ages up to ~30Myr, with a median age of ~12Myr. This indicates that star formation has proceeded over an extended period of time, although we cannot discriminate between an extended episode and a series of short and frequent bursts that are not resolved in time. While the younger PMS population preferentially occupies the central regions of the cluster, older PMS objects are more uniformly distributed across the field and are remarkably few at the very centre of the cluster. We attribute this latter effect to photoevaporation of the older circumstellar discs caused by the massive ionising members of R136.

preprint2011arXiv

The Hubble Space Telescope Wide Field Camera 3 Early Release Science data: Panchromatic Faint Object Counts for 0.2-2 microns wavelength

We describe the Hubble Space Telescope (HST) Wide Field Camera 3 (WFC3) Early Release Science (ERS) observations in the Great Observatories Origins Deep Survey (GOODS) South field. The new WFC3 ERS data provide calibrated, drizzled mosaics in the UV filters F225W, F275W, and F336W, as well as in the near-IR filters F098M (Ys), F125W (J), and F160W (H) with 1-2 HST orbits per filter. Together with the existing HST Advanced Camera for Surveys (ACS) GOODS-South mosaics in the BViz filters, these panchromatic 10-band ERS data cover 40-50 square arcmin at 0.2-1.7 μm in wavelength at 0.07-0.15" FWHM resolution and 0.090" Multidrizzled pixels to depths of AB\simeq 26.0-27.0 mag (5-σ) for point sources, and AB\simeq 25.5-26.5 mag for compact galaxies. In this paper, we describe: a) the scientific rationale, and the data taking plus reduction procedures of the panchromatic 10-band ERS mosaics; b) the procedure of generating object catalogs across the 10 different ERS filters, and the specific star-galaxy separation techniques used; and c) the reliability and completeness of the object catalogs from the WFC3 ERS mosaics. The excellent 0.07-0.15" FWHM resolution of HST/WFC3 and ACS makes star- galaxy separation straightforward over a factor of 10 in wavelength to AB\simeq 25-26 mag from the UV to the near-IR, respectively.

preprint2011arXiv

Using H-alpha Morphology and Surface Brightness Fluctuations to Age-Date Star Clusters in M83

We use new WFC3 observations of the nearby grand design spiral galaxy M83 to develop two independent methods for estimating the ages of young star clusters. The first method uses the physical extent and morphology of Halpha emission to estimate the ages of clusters younger than tau ~10 Myr. It is based on the simple premise that the gas in very young (tau < few Myr) clusters is largely coincident with the cluster stars, is in a small, ring-like structure surrounding the stars in slightly older clusters (e.g., tau ~5 Myr), and is in a larger ring-like bubble for still older clusters (i.e., ~5-10 Myr). The second method is based on an observed relation between pixel-to-pixel flux variations within clusters and their ages. This method relies on the fact that the brightest individual stars in a cluster are most prominent at ages around 10 Myr, and fall below the detection limit (i.e., M_V < -3.5) for ages older than about 100 Myr. These two methods are the basis for a new morphological classification system which can be used to estimate the ages of star clusters based on their appearance. We compare previous age estimates of clusters in M83 determined from fitting UBVI Halpha measurements using predictions from stellar evolutionary models with our new morphological categories and find good agreement at the ~95% level. The scatter within categories is ~0.1 dex in log tau for young clusters (<10 Myr) and ~0.5 dex for older (>10 Myr) clusters. A by-product of this study is the identification of 22 "single-star" HII regions in M83, with central stars having ages ~4 Myr.

preprint2010arXiv

Anatomy of a post-starburst minor merger: a multi-wavelength WFC3 study of NGC 4150

(Abridged) We present a spatially-resolved near-UV/optical study of NGC 4150, using the Wide Field Camera 3 (WFC3) on board the Hubble Space Telescope. Previous studies of this early-type galaxy (ETG) indicate that it has a large reservoir of molecular gas, exhibits a kinematically decoupled core (likely indication of recent merging) and strong, central H_B absorption (indicative of young stars). The core of NGC 4150 shows ubiquitous near-UV emission and remarkable dusty substructure. Our analysis shows this galaxy to lie in the near-UV green valley, and its pixel-by-pixel photometry exhibits a narrow range of near-UV/optical colours that are similar to those of nearby E+A (post-starburst) galaxies. We parametrise the properties of the recent star formation (age, mass fraction, metallicity and internal dust content) in the NGC 4150 pixels by comparing the observed near-UV/optical photometry to stellar models. The typical age of the recent star formation (RSF) is around 0.9 Gyrs, consistent with the similarity of the near-UV colours to post-starburst systems, while the morphological structure of the young component supports the proposed merger scenario. The RSF metallicity, representative of the metallicity of the gas fuelling star formation, is around 0.3 - 0.5 Zsun. Assuming that this galaxy is a merger and that the gas is sourced mainly from the infalling companion, these metallicities plausibly indicate the gas-phase metallicity (GPM) of the accreted satellite. Comparison to the local mass-GPM relation suggests (crudely) that the mass of the accreted system is around 3x10^8 Msun, making NGC 4150 a 1:20 minor merger. A summation of the pixel RSF mass fractions indicates that the RSF contributes about 2-3 percent of the stellar mass. This work reaffirms our hypothesis that minor mergers play a significant role in the evolution of ETGs at late epochs.

preprint2010arXiv

Hubble Space Telescope Imaging of Lyman Alpha Emission at z=4.4

We present the highest redshift detections of resolved Lyman alpha emission, using Hubble Space Telescope/ACS F658N narrowband-imaging data taken in parallel with the Wide Field Camera 3 Early Release Science program in the GOODS CDF-S. We detect Lyman alpha emission from three spectroscopically confirmed z = 4.4 Lyman alpha emitting galaxies (LAEs), more than doubling the sample of LAEs with resolved Lyman alpha emission. Comparing the light distribution between the rest-frame ultraviolet continuum and narrowband images, we investigate the escape of Lyman alpha photons at high redshift. While our data do not support a positional offset between the Lyman alpha and rest-frame ultraviolet (UV) continuum emission, the half-light radii in two out of the three galaxies are significantly larger in Lyman alpha than in the rest-frame UV continuum. This result is confirmed when comparing object sizes in a stack of all objects in both bands. Additionally, the narrowband flux detected with HST is significantly less than observed in similar filters from the ground. These results together imply that the Lyman alpha emission is not strictly confined to its indigenous star-forming regions. Rather, the Lyman alpha emission is more extended, with the missing HST flux likely existing in a diffuse outer halo. This suggests that the radiative transfer of Lyman alpha photons in high-redshift LAEs is complicated, with the interstellar-medium geometry and/or outflows playing a significant role in galaxies at these redshifts.

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.

preprint2010arXiv

Modelling the spectral energy distribution of galaxies. V. The dust and PAH emission SEDs of disk galaxies

We present a self-consistent model of the spectral energy distributions (SEDs) of spiral galaxies from the ultraviolet (UV) to the mid-infrared (MIR)/far-infrared (FIR)/submillimeter (submm) based on a full radiative transfer calculation of the propagation of starlight in galaxy disks. This model predicts not only the total integrated energy absorbed in the UV/optical and re-emitted in the infrared/submm, but also the colours of the dust emission based on an explicit calculation of the strength and colour of the UV/optical radiation fields heating the dust, and incorporating a full calculation of the stochastic heating of small dust grains and PAH molecules. The geometry of the translucent components of the model is empirically constrained using the results from the radiation transfer analysis of Xilouris et al. on spirals in the middle range of the Hubble sequence, while the geometry of the optically thick components is constrained from physical considerations with a posteriori checks of the model predictions with observational data. These geometrical constraints enable the dust emission to be predicted in terms of a minimum set of free parameters: the central face-on dust opacity in the B-band tau^f_B, a clumpiness factor F for the star-forming regions, the star-formation rate SFR, the normalised luminosity of the old stellar population old and the bulge-to-disk ratio B/D. We show that these parameters are almost orthogonal in their predicted effect on the colours of the dust/PAH emission. The results of the calculations are made available in the form of a large library of simulated dust emission SEDs spanning the whole parameter space of our model, together with the corresponding library of dust attenuation calculated using the same model. (see full abstract in the paper)

preprint2010arXiv

Supernova Remnants and the Interstellar Medium of M83: Imaging & Photometry with WFC3 on HST

We present Wide Field Camera 3 images taken with the Hubble Space Telescope within a single field in the southern grand design star-forming galaxy M83. Based on their size, morphology and photometry in continuum-subtracted H$α$, [\SII], H$β$, [\OIII] and [\OII] filters, we have identified 60 supernova remnant candidates, as well as a handful of young ejecta-dominated candidates. A catalog of these remnants, their sizes and, where possible their H$α$ fluxes are given. Radiative ages and pre-shock densities are derived from those SNR which have good photometry. The ages lie in the range $2.62 < log(τ_{\rm rad}/{\rm yr}) < 5.0$, and the pre-shock densities at the blast wave range over $0.56 < n_0/{\rm cm^{-3}} < 1680$. Two populations of SNR have been discovered. These divide into a nuclear and spiral arm group and an inter-arm population. We infer an arm to inter-arm density contrast of 4. The surface flux in diffuse X-rays is correlated with the inferred pre-shock density, indicating that the warm interstellar medium is pressurised by the hot X-ray plasma. We also find that the interstellar medium in the nuclear region of M83 is characterized by a very high porosity and pressure and infer a SNR rate of one per 70-150 yr for the nuclear ($R<300 $pc) region. On the basis of the number of SNR detected and their radiative ages, we infer that the lower mass of Type II SNe in M83 is $M_{\rm min} = 16^{+7}_ {-5}$ M$_{\odot}$. Finally we give evidence for the likely detection of the remnant of the historical supernova, SN1968L.

preprint2010arXiv

The Luminosity, Mass, and Age Distributions of Compact Star Clusters in M83 Based on HST/WFC3 Observations

The newly installed Wide Field Camera 3 (WFC3) on the Hubble Space Telescope has been used to obtain multi-band images of the nearby spiral galaxy M83. These new observations are the deepest and highest resolution images ever taken of a grand-design spiral, particularly in the near ultraviolet, and allow us to better differentiate compact star clusters from individual stars and to measure the luminosities of even faint clusters in the U band. We find that the luminosity function for clusters outside of the very crowded starburst nucleus can be approximated by a power law, dN/dL \propto L^{alpha}, with alpha = -2.04 +/- 0.08, down to M_V ~ -5.5. We test the sensitivity of the luminosity function to different selection techniques, filters, binning, and aperture correction determinations, and find that none of these contribute significantly to uncertainties in alpha. We estimate ages and masses for the clusters by comparing their measured UBVI,Halpha colors with predictions from single stellar population models. The age distribution of the clusters can be approximated by a power-law, dN/dt propto t^{gamma}, with gamma=-0.9 +/- 0.2, for M > few x 10^3 Msun and t < 4x10^8 yr. This indicates that clusters are disrupted quickly, with ~80-90% disrupted each decade in age over this time. The mass function of clusters over the same M-t range is a power law, dN/dM propto M^{beta}, with beta=-1.94 +/- 0.16, and does not have bends or show curvature at either high or low masses. Therefore, we do not find evidence for a physical upper mass limit, M_C, or for the earlier disruption of lower mass clusters when compared with higher mass clusters, i.e. mass-dependent disruption. We briefly discuss these implications for the formation and disruption of the clusters.

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.

preprint1997arXiv

M87: A Misaligned BL LAC?

The nuclear region of M87 was observed with the Faint Object Spectrograph (FOS) on the Hubble Space Telescope (HST) at 6 epochs, spanning 18 months, after the HST image quality was improved with the deployment of the corrective optics (COSTAR) in December 1993. From the FOS target acquisition data, we have established that the flux from the optical nucleus of M87 varies by a factor ~2 on time scales of ~2.5 months and by as much as 25% over 3 weeks, and remains unchanged (<= 2.5%) on time scales of ~1 day. The changes occur in an unresolved central region <= 5 pc in diameter, with the physical size of the emitting region limited by the observed time scales to a few hundred gravitational radii. The featureless continuum spectrum becomes bluer as it brightens while emission lines remain unchanged. This variability combined with the observations of the continuum spectral shape, strong relativistic boosting and the detection of significant superluminal motions in the jet, strongly suggest that M87 belongs to the class of BL Lac objects but is viewed at an angle too large to reveal the classical BL Lac properties.