Source author record

M. S. Oey

M. S. Oey 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

55works
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

55 published item(s)

preprint2026arXiv

Inflight performance and future improvements for the INtegral Field Ultraviolet Spectrographic Experiment, the first far ultraviolet integral field spectrograph

Integral field spectroscopy allows for spectral mapping of extended sources in a time efficient manner. An integral field unit (IFU) in the ultraviolet on Habitable Worlds Observatory (HWO) could be used to quickly map extended objects like supernova remnants or galaxies and their surroundings, but there are technical challenges to an ultraviolet IFU. The INtegral Field Ultraviolet Spectrographic Experiment (INFUSE), a sounding rocket project, is the first static configuration far ultraviolet integral field spectrograph. INFUSE features an f/16, 0.49m Cassegrain telescope and a 26-element image slicer feeding 26 replica holographic gratings, with spectra imaged by the largest cross-strip microchannel plate detector flown in space. The first launch of INFUSE occurred from White Sands Missile Range on October 29th, 2023, and demonstrated spectral multiplexing, successfully detecting ionizing gas emission in the XA region of the Cygnus Loop. INFUSE will launch again in fall 2025 to observe NGC 2366, a local analog for Green Pea type galaxies, with several enhancements including a xenon-enhanced lithium fluoride + aluminum coated grating, testing the leading flight coating for HWO for the first time. The INFUSE IFU is designed as a pathfinder for a potential IFU mode on HWO, enabling rapid 3D spectroscopy of extended sources.

preprint2022arXiv

Classical OBe Stars as Post-Supernova Runaways: Confirming Binary Origins

Massive binaries play an important role in fields ranging from gravitational wave astronomy to stellar evolution. We provide several lines of evidence that classical OBe stars in the Small Magellanic Cloud (SMC) obtain their rapid rotation from mass and angular momentum transfer in massive binaries, which predicts that the subsequent supernovae should often eject OBe stars into the field. We find that (1) OBe stars have a higher field frequency than OB stars; (2) our cumulative distribution function (CDF) of stellar distances from O stars shows that OBe stars are indeed much more isolated than ordinary OB stars of corresponding spectral types; (3) the CDFs of OBe stars approach that of high-mass X-ray binaries (HMXBs), which are confirmed post-supernova objects; and (4) Oe stars are as isolated from clusters as Be stars, implying that their final masses are relatively independent of their initial masses, consistent with major mass transfer. Lastly, we also find that the spatial distribution of supergiant OBe stars differs from that of classical OBe stars, consistent with the different mechanism responsible for their emission-line spectra.

preprint2022arXiv

Haro 11 -- Untying the knots of the nuclear starburst

Star formation is a clustered process that regulates the structure and evolution of galaxies. We investigate this process in the dwarf galaxy Haro 11, forming stars in three knots (A, B, C). The exquisite resolution of HST imaging allows us to resolve the starburst into tens of bright star clusters. We derive masses between $10^5$ and $10^7\,\rm M_{\odot}$ and ages younger than 20 Myr, using photometric modeling. We observe that the clustered star formation has propagated from knot C (the oldest) through knot A (in between) towards knot B (the youngest). We use aperture-matched ultraviolet and optical spectroscopy (HST + MUSE) to independently study the stellar populations of Haro 11 and determine the physical properties of the stellar populations and their feedback in 1 kpc diameter regions. We discuss these results in light of the properties of the ionised gas within the knots. We interpret the broad blue-shifted components of the optical emission lines as outflowing gas ($v_{max} \sim 400$ km/s). The strongest outflow is detected in knot A with a mass-rate of $\dot{M}_{out}\sim 10\,\rm M_{\odot}/yr$, ten times higher than the star-formation in the same region. Knot B hosts a young and not fully developed outflow, whereas knot C has likely been already evacuated. Because Haro 11 has properties similar to high-redshift unresolved galaxies, our work can additionally aid the understanding of star formation at high redshift, a window that will be opened by upcoming facilities.

preprint2022arXiv

No correlation of the Lyman continuum escape fraction with spectral hardness

The properties that govern the production and escape of hydrogen ionizing photons (Lyman continuum, LyC; with energies >13.6 eV) in star-forming galaxies are still poorly understood, but they are key to identifying and characterizing the sources that reionized the Universe. Here we empirically explore the relationship between the hardness of ionizing radiation and the LyC leakage in a large sample of low-$z$ star-forming galaxies from the recent Hubble Space Telescope Low-$z$ Lyman Continuum Survey. Using Sloan Digital Sky Survey stacks and deep XShooter observations, we investigate the hardness of the ionizing spectra ($Q_{\rm He^+}/Q_{\rm H}$) between 54.4 eV (He$^{+}$) and 13.6 eV (H) from the optical recombination lines HeII 4686A and H$β$ 4861A for galaxies with LyC escape fractions spanning a wide range, $f_{\rm esc} \rm (LyC) \simeq 0 - 90\%$. We find that the observed intensity of HeII/H$β$ is primarily driven by variations in the metallicity, but is not correlated with LyC leakage. Both very strong ($<f_{\rm esc} \rm (LyC)> \simeq 0.5$) and nonleakers ($ < f_{\rm esc} \rm (LyC) > \simeq 0$) present similar observed intensities of HeII/H$β$ at comparable metallicity, between $\simeq 0.01$ and $\simeq 0.02$ for $12 + \log({\rm O/H}) > 8.0$ and $<8.0$, respectively. Our results demonstrate that $Q_{\rm He^+}/Q_{\rm H}$ does not correlate with $f_{\rm esc} \rm (LyC)$, which implies that strong LyC emitters do not show harder ionizing spectra than nonleakers at similar metallicity.

preprint2022arXiv

The Low-Redshift Lyman Continuum Survey I: New, Diverse Local Lyman-Continuum Emitters

The origins of Lyman continuum (LyC) photons responsible for the reionization of the universe are as of yet unknown and highly contested. Detecting LyC photons from the epoch of reionization is not possible due to absorption by the intergalactic medium, which has prompted the development of several indirect diagnostics to infer the rate at which galaxies contribute LyC photons to reionize the universe by studying lower-redshift analogs. We present the Low-redshift Lyman Continuum Survey (LzLCS) comprising measurements made with HST/COS for a z=0.2-0.4 sample of 66 galaxies. After careful processing of the FUV spectra, we obtain a total of 35 Lyman continuum emitters (LCEs) detected with 97.725% confidence, nearly tripling the number of known local LCEs. We estimate escape fractions from the detected LyC flux and upper limits on the undetected LyC flux, finding a range of LyC escape fractions up to 50%. Of the 35 LzLCS LCEs, 12 have LyC escape fractions greater than 5%, more than doubling the number of known local LCEs with cosmologically relevant LyC escape.

preprint2022arXiv

The Low-Redshift Lyman Continuum Survey II: New Insights into LyC Diagnostics

The Lyman continuum (LyC) cannot be observed at the epoch of reionization (z {\gtrsim} 6) due to intergalactic H I absorption. To identify Lyman continuum emitters (LCEs) and infer the fraction of escaping LyC, astronomers have developed various indirect diagnostics of LyC escape. Using measurements of the LyC from the Low-redshift Lyman Continuum Survey (LzLCS), we present the first statistical test of these diagnostics. While optical depth indicators based on Lyα, such as peak velocity separation and equivalent width, perform well, we also find that other diagnostics, such as the [O III]/[O II] flux ratio and star formation rate surface density, predict whether a galaxy is a LCE. The relationship between these galaxy properties and the fraction of escaping LyC flux suggests that LyC escape depends strongly on H I column density, ionization parameter, and stellar feedback. We find LCEs occupy a range of stellar masses, metallicities, star formation histories, and ionization parameters, which may indicate episodic and/or different physical causes of LyC escape.

preprint2022arXiv

The Low-Redshift Lyman Continuum Survey. Unveiling the ISM properties of low-$z$ Lyman continuum emitters

Combining 66 ultraviolet (UV) spectra and ancillary data from the Low-Redshift Lyman Continuum Survey (LzLCS) and 23 LyC observations by earlier studies, we form a statistical sample of star-forming galaxies at $z \sim 0.3$ to study the role of the cold interstellar medium (ISM) gas in the leakage of ionizing radiation. We first constrain the massive star content (ages and metallicities) and UV attenuation, by fitting the stellar continuum with a combination of simple stellar population models. The models, together with accurate LyC flux measurements, allow to determine the absolute LyC photon escape fraction for each galaxy ($f_{\rm esc}^{\rm abs}$). We measure the equivalent widths and residual fluxes of multiple HI and low-ionization state (LIS) lines, and the geometrical covering fraction adopting the picket-fence model. The $f_{\rm esc}^{\rm abs}$ spans a wide range, with a median (0.16, 0.84 quantiles) of 0.04 (0.02, 0.20), and 50 out of the 89 galaxies detected in the LyC. The HI and LIS line equivalent widths scale with the UV luminosity and attenuation, and inversely with the residual flux of the lines. The HI and LIS residual fluxes are correlated, indicating that the neutral gas is spatially traced by the LIS transitions. We find the observed trends of the absorption lines and the UV attenuation are primarily driven by the covering fraction. The non-uniform gas coverage demonstrates that LyC photons escape through low-column density channels in the ISM. The equivalent widths and residual fluxes of the UV lines strongly correlate with $f_{\rm esc}^{\rm abs}$: strong LyC leakers show weak absorption lines, low UV attenuation, and large Ly$α$ equivalent widths. We finally show that simultaneous UV absorption line and dust attenuation measurements can predict, on average, the escape fraction of galaxies and the method can be applied to galaxies across a wide redshift range.

preprint2020arXiv

Accretion-Driven Sources in Spatially Resolved Ly$α$ Emitters

Ly$α$ emission is a standard tracer of starburst galaxies at high redshift. However, a number of local Ly$α$ emitters (LAEs) are X-ray sources, suggesting a possible origin of Ly$α$ photons other than young, hot stars, and which may be active at much later ages relative to the parent starburst. Resolved, nearby LAEs offer the opportunity to discriminate between diffuse X-ray emission arising from supernova-heated gas, high-mass X-ray binaries (HMXBs), or low-luminosity active galactic nuclei (LLAGN). We examine archival X-ray imaging from Chandra and XMM-Newton for 11 galaxies with spatially resolved Ly$α$ imaging to determine the luminosity, morphology, and spectral hardness of the X-ray sources. The data are consistent with 9 of the 12, bright Ly$α$ sources being driven by luminous, $10^{40}$ erg s$^{-1}$ X-ray sources. Half of the 8 Chandra sources are unresolved. The data suggest that nuclear activity, whether from LLAGN or nuclear starbursts, may play an important role in Ly$α$ emission. Our results also suggest a significant link between Ly$α$ emission and HMXBs, ULXs, and/or LLAGN, which would imply that Ly$α$ may be generated over timescales 1 - 2 orders of magnitude longer than produced by photoionization from OB stars. This highlights a critical need to quantify the relative contributions of different sources across cosmic time, to interpret Ly$α$ observations and the resulting properties of distant galaxies.

preprint2020arXiv

Size distribution of superbubbles

We consider the size distribution of superbubbles in a star forming galaxy. Previous studies have tried to explain the distribution by using adiabatic self-similar evolution of wind driven bubbles, assuming that bubbles stall when pressure equilibrium is reached. We show, with the help of hydrodynamical numerical simulations, that this assumption is not valid. We also include radiative cooling of shells. In order to take into account non-thermal pressure in the ambient medium, we assume an equivalent higher temperature than implied by thermal pressure alone. Assuming that bubbles stall when the outer shock speed becomes comparable to the ambient sound speed (which includes non-thermal components), we recover the size distribution with a slope of $\sim -2.7$ for typical values of ISM pressure in Milky Way, which is consistent with observations. Our simulations also allow us to follow the evolution of size distribution in the case of different values of non-thermal pressure, and we show that the size distribution steepens with lower pressure, to slopes intermediate between only-growing and only-stalled cases.

preprint2019arXiv

Catastrophic Cooling in Superwinds: Line Emission and Non-equilibrium Ionization

Outflows are a pervasive feature of mechanical feedback from super star clusters (SSC) in starburst galaxies, playing a fundamental role in galaxy evolution. Observations are now starting to confirm that outflows can undergo catastrophic cooling, suppressing adiabatic superwinds. Here we present a suite of one-dimensional, hydrodynamic simulations that study the ionization structure of these outflows and the resulting line emission generated by the cooling gas. We use the non-equilibrium atomic chemistry package within MAIHEM, our modified version of FLASH, which evolves the ionization state of the gas and computes the total cooling rate on an ion-by-ion basis. We find that catastrophically cooling models produce strong nebular line emission compared to adiabatic outflows. We also show that such models exhibit non-equilibrium conditions, thereby generating more highly ionized states than equivalent equilibrium models. When including photoionization from the parent SSC, catastrophically cooling models show strong C IV λ1549 and O VI λ1037 emission. For density bounded photoionization, He II λ1640, λ4686, C III] λ1908, Si IV λ1206, and Si III λ1400 are also strongly enhanced. These lines are seen in extreme starbursts where catastrophic cooling is likely to occur, suggesting that they may serve as diagnostics of such conditions. The higher ionization generated by these flows may help to explain line emission that cannot be attributed to SSC photoionization alone.

preprint2019arXiv

New Insights on Ly-alpha and Lyman Continuum Radiative Transfer in the Greenest Peas

As some of the only Lyman continuum (LyC) emitters at z~0, Green Pea (GP) galaxies are possible analogs of the sources that reionized the universe. We present HST COS spectra of 13 of the most highly ionized GPs, with [O III]/[O II]=6-35, and investigate correlations between Ly-alpha, galaxy properties, and low-ionization UV lines. Galaxies with high [O III]/[O II] have higher H-alpha equivalent widths (EWs), and high intrinsic Ly-alpha production may explain the prevalence of high Ly-alpha EWs among GPs. While Ly-alpha escape fraction is closely linked to low gas covering fractions, implying a clumpy gas geometry, narrow Ly-alpha velocity peak separation (delta_v,LyA) correlates with the ionization state, suggesting a density-bounded geometry. We therefore suggest that delta_v,LyA may trace the residual transparency of low-column-density pathways. Metallicity is associated with both [O III]/[O II] and delta_v,LyA. This trend may result from catastrophic cooling around low-metallicity star clusters, which generates a compact geometry of dense clouds within a low-density inter-clump medium. We find that the relative strength of low-ionization UV emission to absorption correlates with Ly-alpha emission strength and is related to Ly-alpha profile shape. However, as expected for optically thin objects, the GPs with the lowest delta_v,LyA show both weak low-ionization emission and weak absorption. The strengths of the low-ionization absorption and emission lines in a stacked spectrum do not correspond to any individual spectrum. Galaxies with high [O III]/[O II] contain a high fraction of LyC emitter candidates, but [O III]/[O II] alone is an insufficient diagnostic of LyC escape.

preprint2016arXiv

Classical Oe Stars in the Field of the Small Magellanic Cloud

We present $29\pm1$ classical Oe stars from RIOTS4, a spatially complete, spectroscopic survey of Small Magellanic Cloud (SMC) field OB stars. The two earliest are O6e stars, and four are earlier than any Milky Way (MW) Oe stars. We also find ten Ope stars, showing He~\textsc{i} infill and/or emission; five appear to be at least as hot as $\sim$O7.5e stars. The hottest, star 77616, shows He~\textsc{ii} disk emission, suggesting that even the hottest O stars can form decretion disks, and offers observational support for theoretical predictions that the hottest, fastest rotators can generate He$^+$-ionizing atmospheres. Our data also demonstrate that Ope stars correspond to Oe stars earlier than O7.5e with strong disk emission. We find that in the SMC, Oe stars extend to earlier spectral types than in the MW, and our SMC Oe/O frequency, $0.26\pm0.04$, is much greater than the MW value, $0.03\pm0.01$. These results are consistent with angular momentum transport by stronger winds suppressing decretion disk formation at higher metallicity. In addition, our SMC field Oe star frequency is indistinguishable from that for clusters, which is consistent with the similarity between rotation rates in these environments, and contrary to the pattern for MW rotation rates. Thus, our findings strongly support the viscous decretion disk model and confirm that Oe stars are the high-mass extension of the Be phenomenon. Additionally, we find that Fe~\textsc{ii} emission occurs among Oe stars later than O7.5e with massive disks, and we revise a photometric criterion for identifying Oe stars to $J-[3.6] \geq 0.1$.

preprint2016arXiv

Direct Detection of Lyman Continuum Escape from Local Starburst Galaxies with the Cosmic Origins Spectrograph

We report on the detection of Lyman continuum radiation in two nearby starburst galaxies. Tol 0440-381, Tol 1247-232 and Mrk 54 were observed with the Cosmic Origins Spectrograph onboard the Hubble Space Telescopes. The three galaxies have radial velocities of ~13,000 km/s, permitting a ~35 A window on the restframe Lyman continuum shortward of the Milky Way Lyman edge at 912 A. The chosen instrument configuration using the G140L grating covers the spectral range from 912 to 2,000 Å. We developed a dedicated background subtraction method to account for temporal and spatial background variations of the detector, which is crucial at the low flux levels around 912 A. This modified pipeline allowed us to significantly improve the statistical and systematic detector noise and will be made available to the community. We detect Lyman continuum in all three galaxies. However, we conservatively interpret the emission in Tol 0440-381 as an upper limit due to possible contamination by geocoronal Lyman series lines. We determined the current star-formation properties from the far-ultraviolet continuum and spectral lines and used synthesis models to predict the Lyman continuum radiation emitted by the current population of hot stars. We discuss the various model uncertainties such as, among others, atmospheres and evolution models. Lyman continuum escape fractions were derived from a comparison between the observed and predicted Lyman continuum fluxes. Tol 1247-232, Mrk 54 and Tol 0440-381 have absolute escape fractions of (4.5 +/- 1.2)%, (2.5 +/- 0.72)% and <(7.1 +/- 1.1)%, respectively.

preprint2015arXiv

From HI to Stars: HI Depletion in Starbursts and Star-Forming Galaxies in the ALFALFA H-alpha Survey

HI in galaxies traces the fuel for future star formation and reveals the effects of feedback on neutral gas. Using a statistically uniform, HI-selected sample of 565 galaxies from the ALFALFA H-alpha survey, we explore HI properties as a function of star formation activity. ALFALFA H-alpha provides R-band and H-alpha imaging for a volume-limited subset of the 21-cm ALFALFA survey. We identify eight starbursts based on H-alpha equivalent width and six with enhanced star formation relative to the main sequence. Both starbursts and non-starbursts have similar HI to stellar mass ratios (MHI/M*), which suggests that feedback is not depleting the starbursts' HI. Consequently, the starbursts do have shorter HI depletion times (t_dep), implying more efficient HI-to-H2 conversion. While major mergers likely drive this enhanced efficiency in some starbursts, the lowest mass starbursts may experience periodic bursts, consistent with enhanced scatter in t_dep at low M*. Two starbursts appear to be pre-coalescence mergers; their elevated MHI/M* suggest that HI-to-H2 conversion is still ongoing at this stage. By comparing with the GASS sample, we find that t_dep anti-correlates with stellar surface density for disks, while spheroids show no such trend. Among early-type galaxies, t_dep does not correlate with bulge-to-disk ratio; instead, the gas distribution may determine the star formation efficiency. Finally, the weak connection between galaxies' specific star formation rates and MHI/M* contrasts with the well-known correlation between MHI/M* and color. We show that dust extinction can explain the HI-color trend, which may arise from the relationship between M*, MHI, and metallicity.

preprint2015arXiv

The Runaways and Isolated O-Type Star Spectroscopic Survey of the SMC (RIOTS4)

We present the Runaways and Isolated O-Type Star Spectroscopic Survey of the SMC (RIOTS4), a spatially complete survey of uniformly selected field OB stars that covers the entire star-forming body of the SMC. Using the IMACS multislit spectrograph and MIKE echelle spectrograph on the Magellan telescopes, we obtained spectra of 374 early-type field stars that are at least 28 pc from any other OB candidates. We also obtained spectra of an additional 23 field stars in the SMC bar identified from slightly different photometric criteria. Here, we present the observational catalog of stars in the RIOTS4 survey, including spectral classifications and radial velocities. For three multi-slit fields covering 8% of our sample, we carried out monitoring observations over 9-16 epochs to study binarity, finding a spectroscopic, massive binary frequency of at least $\sim$60% in this subsample. Classical Oe/Be stars represent a large fraction of RIOTS4 (42%), occurring at much higher frequency than in the Galaxy, consistent with expectation at low metallicity. RIOTS4 confirmed a steep upper IMF in the field, apparently caused by the inability of the most massive stars to form in the smallest clusters. Our survey also yields evidence for in-situ field OB star formation, and properties of field emission-line star populations, including sgB[e] stars and classical Oe/Be stars. We also discuss the radial velocity distribution and its relation to SMC kinematics and runaway stars. RIOTS4 presents a first quantitative characterization of field OB stars in an external galaxy, including the contributions of sparse, but normal, star formation; runaway stars; and candidate isolated star formation.

preprint2014arXiv

Discovery of the Earliest-Type Oe Stars in the Small Magellanic Cloud

No classical Oe/Be stars with spectral type earlier than O7.5e have been identified to date in the Milky Way (MW). This is consistent with the decretion disk model because strong stellar winds cause early-type O stars to lose angular momentum, thereby preventing them from rotating fast enough to spin out decretion disks. How- ever, metal-poor O stars have weaker stellar winds, allowing the stars to retain angular momentum. Therefore, low-metallicity environments should promote the formation of Oe stars, including those of earlier spectral types than observed in high-metallicity en- vironments. Using the RIOTS4 survey, a spatially complete spectroscopic survey of Small Magellanic Cloud (SMC) field OB stars taken with the IMACS multi-slit spec- trograph at the Magellan Baade Telescope, we identify 25-31 SMC field Oe stars, which account for 20-28% of SMC field O stars. This fraction is significantly higher than in the MW, where < 10-15% of O stars display the Be phenomenon. We also present 5-7 Oe stars of spectral type ranging from O5.5e to O7e, all earlier spectral types than the earliest MW Oe star. These early type Oe stars represent 20-23% of our SMC Oe stars, a dramatic increase compared to the MW, where no Oe stars have been identified with these early spectral types. Thus, the higher frequencies of Oe stars and their earlier spectral range in the metal-poor SMC are consistent with the decretion disk model.

preprint2014arXiv

Ionization by Massive Young Clusters as Revealed by Ionization-Parameter Mapping

Ionization-parameter mapping (IPM) is a powerful technique for tracing the optical depth of Lyman continuum radiation from massive stars. Using narrow-band line-ratio maps, we examine trends in radiative feedback from ordinary HII regions of the Magellanic Clouds and nearby starburst galaxies. We find that the aggregate escape fraction for the Lyman continuum is sufficient to ionize the diffuse, warm ionized medium in the Magellanic Clouds, and that more luminous nebulae are more likely to be optically thin. We apply ionization-parameter mapping to entire starburst galaxies, revealing ionization cones in two nearby starbursts. Within the limits of our small sample, we examine the conditions for the propagation of ionizing radiation beyond the host galaxies.

preprint2014arXiv

Linking Ly-alpha and Low-Ionization Transitions at Low Optical Depth

We suggest that low optical depth in the Lyman continuum (LyC) may relate the Ly-alpha emission, C II and Si II absorption, and C II* and Si II* emission seen in high-redshift galaxies. We base this analysis on Hubble Space Telescope COS spectra of four Green Pea (GP) galaxies, which may be analogs of z>2 Ly-alpha emitters (LAEs). In the two GPs with the strongest Ly-alpha emission, the Ly-alpha line profiles show reduced signs of resonant scattering. Instead, the Ly-alpha profiles resemble the H-alpha line profiles of evolved star ejecta, suggesting that the Ly-alpha emission originates from a low column density and similar outflow geometry. The weak C II absorption and presence of non-resonant C II* emission in these GPs support this interpretation and imply a low LyC optical depth along the line of sight. In two additional GPs, weak Ly-alpha emission and strong C II absorption suggest a higher optical depth. These two GPs differ in their Ly-alpha profile shapes and C II* emission strengths, however, indicating different inclinations of the outflows to our line of sight. With these four GPs as examples, we explain the observed trends linking Ly-alpha, C II, and C II* in stacked LAE spectra, in the context of optical depth and geometric effects. Specifically, in some galaxies with strong Ly-alpha emission, a low LyC optical depth may allow Ly-alpha to escape with reduced scattering. Furthermore, C II absorption, C II* emission, and Ly-alpha profile shape can reveal the optical depth, constrain the orientation of neutral outflows in LAEs, and identify candidate LyC emitters.

preprint2014arXiv

Non-LTE Abundances in OB stars: Preliminary Results for 5 Stars in the Outer Galactic Disk

The aim of this study is to analyse and determine elemental abundances for a large sample of distant B stars in the outer Galactic disk in order to constrain the chemical distribution of the Galactic disk and models of chemical evolution of the Galaxy. Here, we present preliminary results on a few stars along with the adopted methodology based on securing simultaneous O and Si ionization equilibria with consistent NLTE model atmospheres.

preprint2014arXiv

The Origin and Optical Depth of Ionizing Photons in the Green Pea Galaxies

Our understanding of radiative feedback and star formation in galaxies at high redshift is hindered by the rarity of similar systems at low redshift. However, the recently identified Green Pea (GP) galaxies are similar to high-redshift galaxies in their morphologies and star formation rates and are vital tools for probing the generation and transmission of ionizing photons. The GPs contain massive star clusters that emit copious amounts of high-energy radiation, as indicated by intense [OIII] 5007 emission and HeII 4686 emission. We focus on six GP galaxies with high ratios of [O III] 5007,4959/[O II] 3727 ~10 or more. Such high ratios indicate gas with a high ionization parameter or a low optical depth. The GP line ratios and ages point to chemically homogeneous massive stars, Wolf-Rayet stars, or shock ionization as the most likely sources of the He II emission. Models including shock ionization suggest that the GPs may have low optical depths, consistent with a scenario in which ionizing photons escape along passageways created by recent supernovae. The GPs and similar galaxies can shed new light on cosmic reionization by revealing how ionizing photons propagate from massive star clusters to the intergalactic medium.

preprint2013arXiv

A Sample of OB Stars That Formed in the Field

We present a sample of 14 OB stars in the Small Magellanic Cloud that meet strong criteria for having formed under extremely sparse star-forming conditions in the field. These stars are a minimum of 28 pc in projection from other OB stars, and they are centered within symmetric, round HII regions. They show no evidence of bow shocks, implying that the targets are not transverse runaway stars. Their radial velocities relative to local HI also indicate that they are not line-of-sight runaway stars. A friends-of-friends analysis shows that 9 of the objects present a few low-mass companion stars, with typical mass ratios for the two highest-mass stars of around 0.1. This further substantiates that these OB stars formed in place, and that they can and do form in extremely sparse conditions. This poses strong constraints on theories of star formation and challenges proposed relations between cluster mass and maximum stellar mass.

preprint2013arXiv

Exploring the Milky Way stellar disk. A detailed elemental abundance study of 714 F and G dwarf stars in the Solar neighbourhood

ABRIDGED: METHODS: We have conducted a high-resolution spectroscopic study of 714 F and G dwarf and subgiant stars in the Solar neighbourhood. The star sample has been kinematically selected to trace the Galactic thin and thick disks to their extremes...... The determination of stellar parameters and elemental abundances is based on a standard 1-D LTE analysis using equivalent width measurements in high-resolution (R=40000-110000) and high signal-to-noise (S/N=150-300) spectra obtained with..... RESULTS: .... Our data show that there is an old and alpha-enhanced disk population, and a younger and less alpha-enhanced disk population. While they overlap greatly in metallicity between -0.7<[Fe/H]<+0.1, they show a bimodal distribution in [alpha/Fe]. This bimodality becomes even clearer if stars where stellar parameters and abundances show larger uncertainties (Teff<5400 K) are discarded, showing that it is important to constrain the data set to a narrow range in the stellar parameters if small differences between stellar populations are to be revealed. We furthermore find that the alpha-enhanced population has orbital parameters placing the stellar birthplaces in the inner Galactic disk while the low-alpha stars mainly come from the outer Galactic disk........... We furthermore have discovered that a standard 1-D, LTE analysis, utilising ionisation and excitation balance of Fe I and Fe II lines produces a flat lower main sequence. As the exact cause for this effect is unclear we chose to apply an empirical correction. Turn-off, and more evolved, stars, appears to be un-affected.

preprint2013arXiv

New Constraints on the Escape of Ionizing Photons From Starburst Galaxies Using Ionization-Parameter Mapping

The fate of ionizing radiation in starburst galaxies is key to understanding cosmic reionization. However, the galactic parameters on which the escape fraction of ionizing radiation depend are not well understood. Ionization-parameter mapping provides a simple, yet effective, way to study the radiative transfer in starburst galaxies. We obtain emission-line ratio maps of [SIII]/[SII] for six, nearby, dwarf starbursts: NGC 178, NGC 1482, NGC 1705, NGC 3125, NGC 7126, and He 2-10. The narrow-band images are obtained with the Maryland-Magellan Tunable Filter at Las Campanas Observatory. Using these data, we previously reported the discovery of an optically thin ionization cone in NGC 5253, and here we also discover a similar ionization cone in NGC 3125. This latter cone has an opening angle of 40+/-5 degrees (0.4 ster), indicating that the passageways through which ionizing radiation may travel correspond to a small solid angle. Additionally, there are three sample galaxies that have winds and/or superbubble activity, which should be conducive to escaping radiation, yet they are optically thick. These results support the scenario that an orientation bias limits our ability to directly detect escaping Lyman continuum in many starburst galaxies. A comparison of the star-formation properties and histories of the optically thin and thick galaxies is consistent with the model that high escape fractions are limited to galaxies that are old enough (> 3 Myr) for mechanical feedback to have cleared optically thin passageways in the ISM, but young enough (< 5 Myr) that the ionizing stars are still present.

preprint2013arXiv

The Optical Depth of H II Regions in the Magellanic Clouds

We exploit ionization-parameter mapping as a powerful tool to measure the optical depth of star-forming HII regions. Our simulations using the photoionization code CLOUDY and our new, SURFBRIGHT surface brightness simulator demonstrate that this technique can directly diagnose most density-bounded, optically thin nebulae using spatially resolved emission line data. We apply this method to the Large and Small Magellanic Clouds, using the data from the Magellanic Clouds Emission Line Survey. We generate new HII region catalogs based on photoionization criteria set by the observed ionization structure in the [SII]/[OIII] ratio and Ha surface brightness. The luminosity functions from these catalogs generally agree with those from Ha-only surveys. We then use ionization-parameter mapping to crudely classify all the nebulae into optically thick vs optically thin categories, yielding fundamental new insights into Lyman continuum radiation transfer. We find that in both galaxies, the frequency of optically thin objects correlates with Ha luminosity, and that the numbers of these objects dominate above log L\geq37.0. The frequencies of optically thin objects are 40% and 33% in the LMC and SMC, respectively. Similarly, the frequency of optically thick regions correlates with H I column density, with optically thin objects dominating at the lowest N(HI). The integrated escape luminosity of ionizing radiation is dominated by the largest regions, and corresponds to luminosity-weighted, ionizing escape fractions from the H II region population of \geq0.42 and \geq0.40 in the LMC and SMC, respectively. These values correspond to global galactic escape fractions of 4% and 11%, respectively. This is sufficient to power the ionization rate of the observed diffuse ionized gas in both galaxies. Our results suggest the possibility of significant galactic escape fractions of Lyman continuum radiation.

preprint2013arXiv

The Origin and Optical Depth of Ionizing Radiation in the "Green Pea" Galaxies

Although Lyman continuum (LyC) radiation from star-forming galaxies likely drove the reionization of the Universe, observations of star-forming galaxies at low redshift generally indicate low LyC escape fractions. However, the extreme [O III]/[O II] ratios of the z=0.1-0.3 Green Pea galaxies may be due to high escape fractions of ionizing radiation. To analyze the LyC optical depths and ionizing sources of these rare, compact starbursts, we compare nebular photoionization and stellar population models with observed emission lines in the Peas' SDSS spectra. We focus on the six most extreme Green Peas, the galaxies with the highest [O III]/[O II] ratios and the best candidates for escaping ionizing radiation. The Balmer line equivalent widths and He I λ3819 emission in the extreme Peas support young ages of 3-5 Myr, and He II λ4686 emission in five extreme Peas signals the presence of hard ionizing sources. Ionization by active galactic nuclei or high-mass X-ray binaries is inconsistent with the Peas' line ratios and ages. Although stacked spectra reveal no Wolf-Rayet (WR) features, we tentatively detect WR features in the SDSS spectra of three extreme Peas. Based on the Peas' ages and line ratios, we find that WR stars, chemically homogeneous O stars, or shocks could produce the observed He II emission. If hot stars are responsible, the Peas' optical depths are ambiguous. However, accounting for emission from shocks lowers the inferred optical depth and suggests that the Peas may be optically thin. The Peas' ages likely optimize the escape of Lyman-continuum radiation; they are old enough for supernovae and stellar winds to reshape the interstellar medium, but young enough to possess large numbers of UV-luminous O or WR stars.

preprint2012arXiv

Abundance trends in the inner and outer Galactic disk

Based on high-resolution spectra obtained with the MIKE spectrograph on the Magellan telescopes we present detailed elemental abundances for 64 red giant stars in the inner and outer Galactic disk. For the inner disk sample (4-7 kpc from the Galactic centre) we find that stars with both thin and thick disk abundance patterns are present while for Galactocentric distances beyond 10 kpc, we only find chemical patterns associated with the local thin disk, even for stars far above the Galactic plane. Our results show that the relative densities of the thick and thin disks are dramatically different from the solar neighbourhood, and we therefore suggest that the radial scale length of the thick disk is much shorter than that of the thin disk. A thick disk scale-length of L_{thick}=2.0 kpc, and L_{thin}=3.8 kpc for the thin disk, better match the data.

preprint2012arXiv

Discovery of New, Dust-Poor B[e] Supergiants in the Small Magellanic Cloud

We present the discovery of three new B[e] supergiants (sgB[e] stars) in the Small Magellanic Cloud (SMC). All three stars (R15, R38, and R48) were identified in the course of our Runaways and Isolated O Type Star Spectroscopic Survey of the SMC (RIOTS4). The stars show optical spectra that closely resemble those of previously known B[e] stars, presenting numerous low-ionization forbidden and permitted emission lines such as [Fe II] and Fe II. Furthermore, our stars have luminosities of log(L/L_sun) > 4, demonstrating that they are supergiants. However, we find lower infrared excesses and weaker forbidden emission lines than for previously identified B[e] supergiants. Thus our stars appear to either have less material in their circumstellar disks than other sgB[e] stars, or the circumstellar material has lower dust content. We suggest that these may constitute a new subclass of dust-poor sgB[e] stars.

preprint2012arXiv

Projected Rotational Velocities and Stellar Characterization of 350 B Stars in the Nearby Galactic Disk

Projected rotational velocities (vsini) are presented for a sample of 350 early B-type main sequence stars in the nearby Galactic disk. The stars are located within ~1.5 kpc from the Sun, and the great majority within 700 pc. The analysis is based on high-resolution spectra obtained with the MIKE spectrograph on the Magellan Clay 6.5-m telescope at the Las Campanas Observatory in Chile.Spectral types were estimated based on relative intensities of some key line absorption ratios and comparisons to synthetic spectra. Effective temperatures were estimated from the reddening-free Q index, and projected rotational velocities were then determined via interpolation on a published grid that correlates the synthetic full width at half maximum of the He I lines at 4026, 4388 and 4471 A with vsini. As the sample has been selected solely on the basis of spectral types it contains an selection of B stars in the field, in clusters, and in OB associations. The vsini distribution obtained for the entire sample is found to be essentially flat for vsini values between 0-150 km/s, with only a modest peak at low projected rotational velocities. Considering subsamples of stars, there appears to be a gradation in the vsini distribution with the field stars presenting a larger fraction of the slow rotators and the cluster stars distribution showing an excess of stars with vsini between 70 and 130 km/s. Furthermore, for a subsample of potential runaway stars we find that the vsini distribution resembles the distribution seen in denser environments, which could suggest that these runaway stars have been subject to dynamical ejection mechanisms.

preprint2012arXiv

Single-Star HII Regions as a Probe of Massive Star Spectral Energy Distributions

The shape of the OB-star spectral energy distribution is a critical component in many diagnostics of the ISM and galaxy properties. We use single-star HII regions from the LMC to quantitatively examine the ionizing SEDs from widely available CoStar, TLUSTY, and WM-basic atmosphere grids. We evaluate the stellar atmosphere models by matching the emission-line spectra that they predict from CLOUDY photoionization simulations with those observed from the nebulae. The atmosphere models are able to reproduce the observed optical nebular line ratios, except at the highest energy transitions > 40 eV, assuming that the gas distribution is non-uniform. Overall we find that simulations using WM-basic produce the best agreement with the observed line ratios. The rate of ionizing photons produced by the model SEDs is consistent with the rate derived from the \Halpha\ luminosity for standard, log(g) = 4.0 models adopted from the atmosphere grids. However, there is a systematic offset between the rate of ionizing photons from different atmosphere models that is correlated with the relative hardness of the SEDs. In general WM-basic and TLUSTY atmosphere models predict similar effective temperatures, while CoStar predicts effective temperatures that are cooler by a few thousand degrees. We compare our effective temperatures, which depend on the nebular ionization balance, to conventional photospheric-based calibrations from the literature. We suggest that in the future, spectral type to effective temperature calibrations can be constructed from nebular data.

preprint2012arXiv

The Salpeter Slope of the IMF Explained

If we accept a paradigm that star formation is a self-similar, hierarchical process, then the Salpeter slope of the IMF for high-mass stars can be simply and elegantly explained as follows. If the instrinsic IMF at the smallest scales follows a simple -2 power-law slope, then the steepening to the -2.35 Salpeter value results when the most massive stars cannot form in the lowest-mass clumps of a cluster. It is stressed that this steepening MUST occur if clusters form hierarchically from clumps, and the lowest-mass clumps can form stars. This model is consistent with a variety of observations as well as theoretical simulations.

preprint2011arXiv

A First Constraint on the Thick Disk Scale Length: Differential Radial Abundances in K Giants at Galactocentric Radii 4, 8, and 12 kpc

Based on high-resolution spectra obtained with the MIKE spectrograph on the Magellan telescopes we present detailed elemental abundances for 20 red giant stars in the outer Galactic disk, located at Galactocentric distances between 9 and 13 kpc. The outer disk sample is complemented with samples of red giants from the inner Galactic disk and the solar neighbourhood, analysed using identical methods. For Galactocentric distances beyond 10 kpc, we only find chemical patterns associated with the local thin disk, even for stars far above the Galactic plane. Our results show that the relative densities of the thick and thin disks are dramatically different from the solar neighbourhood, and we therefore suggest that the radial scale length of the thick disk is much shorter than that of the thin disk. We make a first estimate of the thick disk scale-length of L_thick=2.0 kpc, assuming L_thin=3.8 kpc for the thin disk. We suggest that radial migration may explain the lack of radial age, metallicity, and abundance gradients in the thick disk, possibly also explaining the link between the thick disk and the metal-poor bulge.

preprint2011arXiv

An Ionization Cone in the Dwarf Starburst Galaxy NGC 5253

There are few observational constraints on how the escape of ionizing photons from starburst galaxies depends on galactic parameters. Here, we report on the first major detection of an ionization cone in NGC 5253, a nearby starburst galaxy. This high-excitation feature is identified by mapping the emission-line ratios in the galaxy using [S III] lambda 9069, [S II] lambda 6716, and H_alpha narrow-band images from the Maryland-Magellan Tunable Filter at Las Campanas Observatory. The ionization cone appears optically thin, which is suggestive of the escape of ionizing photons. The cone morphology is narrow with an estimated solid angle covering just 3% of 4pi steradians, and the young, massive clusters of the nuclear starburst can easily generate the radiation required to ionize the cone. Although less likely, we cannot rule out the possibility of an obscured AGN source. An echelle spectrum along the minor axis shows complex kinematics that are consistent with outflow activity. The narrow morphology of the ionization cone supports the scenario that an orientation bias contributes to the difficulty in detecting Lyman continuum emission from starbursts and Lyman break galaxies.

preprint2011arXiv

Completing the Massive Star Population: Striking Into the Field

As a population, field massive stars are relatively enigmatic, and this review attempts to illuminate this sector of the high-mass stellar population, which comprises 20 - 25% of the massive stars in star-forming galaxies. The statistical properties of the field population are vital diagnostics of star formation theory, cluster dynamical evolution, and stellar evolution. We present evidence that field massive stars originate both in situ and as runaways from clusters, based on the clustering law, IMF, rotation velocities, and individual observed in situ candidate field stars. We compare the known properties of field and cluster massive stars from studies in the Magellanic Clouds and the Galaxy, including our RIOTS4 complete spectroscopic survey of SMC OB stars. In addition to the origin of the field massive stars, we discuss additional properties including binarity, runaway mechanisms, and some evolved spectral types.

preprint2011arXiv

Dust Emission from Evolved and Unevolved HII Regions in the Large Magellanic Cloud

We present a study of the dust properties of 12 classical and superbubble HII regions in the Large Magellanic Cloud. We use infrared photometry from Spitzer (8, 24, 70, and 160 \mum bands), obtained as part of the Surveying the Agents of a Galaxy's Evolution (SAGE) program, along with archival spectroscopic classifications of the ionizing stars to examine the role of stellar sources on dust heating and processing. Our infrared observations show surprisingly little correlation between the emission properties of the dust and the effective temperatures or bolometric magnitudes of stars in the HII regions, suggesting that the HII region evolutionary timescale is not on the order of the dust processing timescale. We find that the infrared emission of superbubbles and classical HII regions shows little differentiation between the two classes, despite the significant differences in age and morphology. We do detect a correlation of the 24 \mum emission from hot dust with the ratio of 70 to 160 \mum flux. This correlation can be modeled as a trend in the temperature of a minority hot dust component, while a majority of the dust remains significantly cooler.

preprint2011arXiv

Observational Constraints on Superbubble X-ray Energy Budgets

The hot, X-ray-emitting gas in superbubbles imparts energy and enriched material to the interstellar medium (ISM) and generates the hot ionized medium, the ISM's high-temperature component. The evolution of superbubble energy budgets is not well understood, however, and the processes responsible for enhanced X-ray emission in superbubbles remain a matter of debate. We present Chandra ACIS-S observations of two X-ray-bright superbubbles in the Large Magellanic Cloud (LMC), DEM L50 (N186) and DEM L152 (N44), with an emphasis on disentangling the true superbubble X-ray emission from non-related diffuse emission and determining the spatial origin and spectral variation of the X-ray emission. An examination of the superbubble energy budgets shows that on the order of 50% of the X-ray emission comes from regions associated with supernova remnant (SNR) impacts. We find some evidence of mass-loading due to swept-up clouds and metallicity enrichment, but neither mechanism provides a significant contribution to the X-ray luminosities. We also find that one of the superbubbles, DEM L50, is likely not in collisional ionization equilibrium. We compare our observations to the predictions of the standard Weaver et al. model and to 1-D hydrodynamic simulations including cavity supernova impacts on the shell walls. Our observations show that mass-loading due to thermal evaporation from the shell walls and SNR impacts are the dominant source of enhanced X-ray luminosities in superbubbles. These two processes should affect most superbubbles, and their contribution to the X-ray luminosity must be considered when determining the energy available for transport to the ISM.

preprint2011arXiv

On the Origin of the Salpeter Slope for the Initial Mass Function

We suggest that the intrinsic, stellar initial mass function (IMF) follows a power-law slope gamma=2, inherited from hierarchical fragmentation of molecular clouds into clumps and clumps into stars. The well-known, logarithmic Salpeter slope GAMMA=1.35 in clusters is then the aggregate slope for all the star-forming clumps contributing to an individual cluster, and it is steeper than the intrinsic slope within individual clumps because the smallest star-forming clumps contributing to any given cluster are unable to form the highest-mass stars. Our Monte Carlo simulations demonstrate that the Salpeter power-law index is the limiting value obtained for the cluster IMF when the lower-mass limits for allowed stellar masses and star-forming clumps are effectively equal, m_lo = M_lo. This condition indeed is imposed for the high-mass IMF tail by the turn-over at the characteristic value m_c ~ 1 M_sun. IMF slopes of GAMMA ~ 2 are obtained if the stellar and clump upper-mass limits are also equal m_up = M_up ~ 100 M_sun, and so our model explains the observed range of IMF slopes between GAMMA ~ 1 to 2. Flatter slopes of GAMMA = 1 are expected when M_lo > m_up, which is a plausible condition in starbursts, where such slopes are suggested to occur. While this model is a simplistic parameterization of the star-formation process, it seems likely to capture the essential elements that generate the Salpeter tail of the IMF for massive stars. These principles also likely explain the IGIMF effect seen in low-density star-forming environments.

preprint2011arXiv

The IMF of Field OB Stars in the Small Magellanic Cloud

The population of field OB stars are an important component of a galaxy's stellar content, representing 20-30% of the massive stars. To study this population, we have undertaken the Runaways and Isolated O Type Star Spectroscopic Survey of the SMC (RIOTS4). RIOTS4 surveys a spatially complete sample of >350 field OB stars in the Small Magellanic Cloud and will serve as a key probe of runaways, binaries, and the stellar IMF in the field massive star population. Here, we focus on the field IMF, which provides an empirical probe of the star-forming process and is a fundamental property of a stellar population. Together with photometry from the OGLE survey, RIOTS4 will yield a definitive stellar IMF for the SMC field massive star population. We present preliminary results that suggests the field IMF is much steeper, Gamma = 2.9, than the canonical stellar IMF of Gamma = 1.35. Despite the steep slope, we see no evidence of a stellar upper mass limit, up to our highest mass star of 65 solar masses.

preprint2011arXiv

The Integrated Diffuse X-ray Emission of the Carina Nebula Compared to Other Massive Star-forming Regions

The Chandra Carina Complex Project (CCCP) has shown that the Carina Nebula displays bright, spatially-complex soft diffuse X-ray emission. Here we `sum up' the CCCP diffuse emission work by comparing the global morphology and spectrum of Carina's diffuse X-ray emission to other famous sites of massive star formation with pronounced diffuse X-ray emission: M17, NGC 3576, NGC 3603, and 30 Doradus. All spectral models require at least two diffuse thermal plasma components to achieve adequate spectral fits, a softer component with kT = 0.2--0.6 keV and a harder component with kT = 0.5--0.9 keV. In several cases these hot plasmas appear to be in a state of non-equilibrium ionization that may indicate recent and current strong shocks. A cavity north of the embedded giant HII region NGC 3576 is the only region studied here that exhibits hard diffuse X-ray emission; this emission appears to be nonthermal and is likely due to a recent cavity supernova, as evidenced by a previously-known pulsar and newly-discovered pulsar wind nebula also seen in this cavity. All of these targets exhibit X-ray emission lines that are not well-modeled by variable-abundance thermal plasmas and that might be attributed to charge exchange at the shock between the hot, tenuous, X-ray-emitting plasma and cold, dense molecular material; this is likely evidence for dust destruction at the many hot/cold interfaces that characterize massive star-forming regions.

preprint2010arXiv

A Multiwavelength Study on the Fate of Ionizing Radiation in Local Starbursts

The fate of ionizing radiation is vital for understanding cosmic ionization, energy budgets in the interstellar and intergalactic medium, and star formation rate indicators. The low observed escape fractions of ionizing radiation have not been adequately explained, and there is evidence that some starbursts have high escape fractions. We examine the spectral energy distributions of a sample of local star-forming galaxies, containing thirteen local starburst galaxies and ten of their ordinary star-forming counterparts, to determine if there exist significant differences in the fate of ionizing radiation in these galaxies. We find that the galaxy-to-galaxy variations in the SEDs is much larger than any systematic differences between starbursts and non-starbursts. For example, we find no significant differences in the total absorption of ionizing radiation by dust, traced by the 24um, 70um, and 160um MIPS bands of the Spitzer Space Telescope, although the dust in starburst galaxies appears to be hotter than that of non-starburst galaxies. We also observe no excess ultraviolet flux in the GALEX bands that could indicate a high escape fraction of ionizing photons in starburst galaxies. The small H-alpha fractions of the diffuse, warm ionized medium in starburst galaxies are apparently due to temporarily boosted H-alpha luminosity within the star-forming regions themselves, with an independent, constant WIM luminosity. This independence of the WIM and starburst luminosities contrasts with WIM behavior in non-starburst galaxies and underscores our poor understanding of radiation transfer in both ordinary and starburst galaxies.

preprint2010arXiv

Gemini GMOS spectroscopy of HeII nebulae in M33

We have carried out a narrow-band survey of the Local Group galaxy, M33, in the HeII4686 emission line, to identify HeII nebulae in this galaxy. With spectroscopic follow-up observations, we confirm three of seven candidate objects, including identification of two new HeII nebulae, BCLMP651, HBW673. We also obtain spectra of associated ionizing stars for all the HII regions, identifying two new WN stars. We demonstrate that the ionizing source for the known HeII nebula, MA 1, is consistent with being the early-type WN star MC8 (M33-WR14), by carrying out a combined stellar and nebular analysis of MC8 and MA1. We were unable to identify the helium ionizing sources for HBW 673 and BCLMP 651, which do not appear to be Wolf-Rayet stars. According to the [OIII]5007/Hbeta vs [NII]6584/Halpha diagnostic diagram, excitation mechanisms apart from hot stellar continuum are needed to account for the nebular emission in HBW 673, which appears to have no stellar source at all.

preprint2010arXiv

OB Stars in Stochastic Regimes

The highest-mass stars have the lowest frequency in the stellar IMF, and they are also the most easily observed stars. Thus, the counting statistics for OB stars provide important tests for the fundamental nature and quantitative parameters of the IMF. We first examine some local statistics for the stellar upper-mass limit itself. Then, we examine the parameter space and statistics for extremely sparse clusters that contain OB stars, in the SMC. We find that thus far, these locally observed counting statistics are consistent with a constant stellar upper-mass limit. The sparse OB star clusters easily fall within the parameter space of Monte Carlo simulations of cluster populations. If the observed objects are representative of their cluster birth masses, their existence implies that the maximum stellar mass is largely independent of the parent cluster mass.

preprint2010arXiv

Spitzer Analysis of HII Region Complexes in the Magellanic Clouds: Determining a Suitable Monochromatic Obscured Star Formation Indicator

HII regions are the birth places of stars, and as such they provide the best measure of current star formation rates (SFRs) in galaxies. The close proximity of the Magellanic Clouds allows us to probe the nature of these star forming regions at small spatial scales. We aim to determine the monochromatic IR band that most accurately traces the bolometric IR flux (TIR), which can then be used to estimate an obscured SFR. We present the spatial analysis, via aperture/annulus photometry, of 16 LMC and 16 SMC HII region complexes using the Spitzer IRAC and MIPS bands. UV rocket data and SHASSA H-alpha data are also included. We find that nearly all of the LMC and SMC HII region SEDs peak around 70um, from ~10 to ~400 pc from the central sources. As a result, the sizes of HII regions as probed by 70um is approximately equal to the sizes as probed by TIR (about 70 pc in radius); the radial profile of the 70um flux, normalized by TIR, is constant at all radii (70um ~ 0.45 TIR); the 1-sigma standard deviation of the 70um fluxes, normalized by TIR, is a lower fraction of the mean (0.05 to 0.12 out to ~220 pc) than the normalized 8, 24, and 160um normalized fluxes (0.12 to 0.52); and these results are invariant between the LMC and SMC. From these results, we argue that 70um is the most suitable IR band to use as a monochromatic obscured star formation indicator because it most accurately reproduces the TIR of HII regions in the LMC and SMC and over large spatial scales. We also explore the general trends of the 8, 24, 70, and 160um bands in the LMC and SMC HII region SEDs, radial surface brightness profiles, sizes, and normalized (by TIR) radial flux profiles. We derive an obscured SFR equation that is modified from the literature to use 70um luminosity, SFR(Mo/yr) = 9.7(0.7)x10^{-44} L(70)(ergs/s), which is applicable from 10 to 300 pc distance from the center of an HII region.

preprint2010arXiv

The first chemical abundance analysis of K giants in the inner Galactic disc

The elemental abundance structure of the Galactic disc has been extensively studied in the solar neighbourhood using long-lived stars such as F and G dwarfs or K and M giants. These are stars whose atmospheres preserve the chemical composition of their natal gas clouds, and are hence excellent tracers of the chemical evolution of the Galaxy. As far as we are aware, there are no such studies of the inner Galactic disc, which hampers our ability to constrain and trace the origin and evolution of the Milky Way. Therefore, we aim in this study to establish the elemental abundance trend(s) of the disc(s) in the inner regions of the Galaxy. Based on equivalent width measurements in high-resolution spectra obtained with the MIKE spectrograph on the Magellan II telescope on Las Campanas in Chile, we determine elemental abundances for 44 K-type red giant stars in the inner Galactic disc, located at Galactocentric distances of 4-7\,kpc. The analysis method is identical to the one recently used on red giant stars in the Galactic bulge and in the nearby thin and thick discs, enabling us to perform a truly differential comparison of the different stellar populations. We present the first detailed elemental abundance study of a significant number of red giant stars in the inner Galactic disc. We find that these inner disc stars show the same type of chemical and kinematical dichotomy as the thin and thick discs show in the solar neighbourhood. The abundance trends of the inner disc agree very well with those of the nearby thick disc, and also to those of the Bulge. The chemical similarities between the Bulge and the Galactic thick disc stellar populations indicate that they have similar chemical histories, and any model trying to understand the formation and evolution of either of the two should preferably incorporate both of them.

preprint2010arXiv

The Herschel revolution: unveiling the morphology of the high mass star formation sites N44 and N63 in the LMC

We study the structure of the medium surrounding sites of high-mass star formation to determine the interrelation between the HII regions and the environment from which they were formed. The density distribution of the surroundings is key in determining how the radiation of the newly formed stars interacts with the surrounds in a way that allows it to be used as a star formation tracer. We present new Herschel/SPIRE 250, 350 and 500 mum data of LHA 120-N44 and LHA 120-N63 in the LMC. We construct average spectral energy distributions (SEDs) for annuli centered on the IR bright part of the star formation sites. The annuli cover ~10-~100 pc. We use a phenomenological dust model to fit these SEDs to derive the dust column densities, characterise the incident radiation field and the abundance of polycyclic aromatic hydrocarbon molecules. We see a factor 5 decrease in the radiation field energy density as a function of radial distance around N63. N44 does not show a systematic trend. We construct a simple geometrical model to derive the 3-D density profile of the surroundings of these two regions. Herschel/SPIRE data have proven very efficient in deriving the dust mass distribution. We find that the radiation field in the two sources behaves very differently. N63 is more or less spherically symmetric and the average radiation field drops with distance. N44 shows no systematic decrease of the radiation intensity which is probably due to the inhomogeneity of the surrounding molecular material and to the complex distribution of several star forming clusters in the region.

preprint2010arXiv

The Sparsest Clusters With O Stars

There is much debate on how high-mass star formation varies with environment, and whether the sparsest star-forming environments are capable of forming massive stars. To address this issue, we have observed eight apparently isolated OB stars in the SMC using HST's Advanced Camera for Surveys. Five of these objects appear as isolated stars, two of which are confirmed to be runaways. The remaining three objects are found to exist in sparse clusters, with <10 companion stars revealed, having masses of 1-4 solar mass. Stochastic effects dominate in these sparse clusters, so we perform Monte Carlo simulations to explore how our observations fit within the framework of empirical, galactic cluster properties. We generate clusters using a simplistic -2 power-law distribution for either the number of stars per cluster (N_*) or cluster mass (M_cl). These clusters are then populated with stars randomly chosen from a Kroupa IMF. We find that simulations with cluster lower-mass limits of M_cl,lo >20 solar mass and N_*,lo >40 match best with observations of SMC and Galactic OB star populations. We examine the mass ratio of the second-most massive and most massive stars (m_max,2/m_max), finding that our observations all exist below the 20th percentile of our simulated clusters. However, all of our observed clusters lie within the parameter space spanned by the simulated clusters, although some are in the lowest 5th percentile frequency. These results suggest that clusters are built stochastically by randomly sampling stars from a universal IMF with a fixed stellar upper-mass limit. In particular, we see no evidence to suggest a m_max - M_cl relation. Our results may be more consistent with core accretion models of star formation than with competitive accretion models, and they are inconsistent with the proposed steepening of the integrated galaxy IMF (IGIMF).

preprint2009arXiv

Outlying HII Regions in HI-Selected Galaxies

We present results from the first systematic search for outlying HII regions, as part of a sample of 96 emission-line point sources (referred to as ELdots - emission-line dots) derived from the NOAO Survey for Ionization in Neutral Gas Galaxies (SINGG). Our automated ELdot-finder searches SINGG narrow-band and continuum images for high equivalent width point sources outside the optical radius of the target galaxy (> 2 X r25 in the R-band). Follow-up longslit spectroscopy and deep GALEX images (exposure time > 1000 s) distinguish outlying HII regions from background galaxies whose strong emission lines ([OIII], Hbeta or [OII]) have been redshifted into the SINGG bandpass. We find that these deep GALEX images can serve as a substitute for spectroscopic follow-up because outlying HII regions separate cleanly from background galaxies in color-color space. We identify seven SINGG systems with outlying massive star formation that span a large range in Halpha luminosities corresponding to a few O stars in the most nearby cases, and unresolved dwarf satellite companion galaxies in the most distant cases. Six of these seven systems feature galaxies with nearby companions or interacting galaxies. Furthermore, our results indicate that some outlying HII regions are linked to the extended-UV disks discovered by GALEX, representing emission from the most massive O stars among a more abundant population of lower mass (or older) star clusters. The overall frequency of outlying HII regions in this sample of gas-rich galaxies is 8 - 11% when we correct for background emission-line galaxy contamination (~75% of ELdots).

preprint2007arXiv

The Survey for Ionization in Neutral-Gas Galaxies: III. Diffuse, Warm Ionized Medium and Escape of Ionizing Radiation

We use the first data release from the SINGG H-alpha survey of HI-selected galaxies to study the quantitative behavior of the diffuse, warm ionized medium (WIM) across the range of properties represented by these 109 galaxies. The mean fraction f_WIM of diffuse ionized gas in this sample is 0.59+/- 0.19, slightly higher than found in previous samples. Since lower surface-brightness galaxies tend to have higher f_WIM, we believe that most of this difference is due to selection effects favoring large, optically-bright, nearby galaxies with high star-formation rates. As found in previous studies, there is no appreciable correlation with Hubble type or total star-formation rate. However, we find that starburst galaxies, defined here by an H-alpha surface brightness > 2.5x 10^39 erg s^-1 kpc^-2 within the H-alpha half-light radius, do show much lower fractions of diffuse H-alpha emission. The cause apparently is not dominated by a lower fraction of field OB stars. However, it is qualitatively consistent with an expected escape of ionizing radiation above a threshold star-formation rate, predicted from our model in which the ISM is shredded by pressure-driven supernova feedback. The HI gas fractions in the starburst galaxies are also lower, suggesting that the starbursts are consuming and ionizing all the gas, and thus promoting regions of density-bounded ionization. If true, these effects imply that some amount of Lyman continuum radiation is escaping from most starburst galaxies, and that WIM properties and outflows from mechanical feedback are likely to be pressure-driven. However, in view of previous studies showing that the escape fraction of ionizing radiation is generally low, it is likely that other factors also drive the low fractions of diffuse ionized gas in starbursts.

preprint2007arXiv

Tracing the Galactic thick disk to Solar metallicities

We show that the Galactic thick disk reaches at least solar metallicities, and that it experienced strong chemical enrichment during a period of ~3 Gyr, ending around 8-9 Gyr ago. This finding puts further constraints on the relation and interface between the thin and thick disks, and their formation processes. Our results are based on a detailed elemental abundance analysis of 261 kinematically selected F and G dwarf stars in the solar neighborhood: 194 likely members of the thick disk and 67 likely members of the thin disk, in the range -1.3<[Fe/H]<+0.4.

preprint2006arXiv

Spitzer Survey of the Large Magellanic Cloud, Surveying the Agents of a Galaxy's Evolution (SAGE) I: Overview and Initial Results

We are performing a uniform and unbiased, ~7x7 degrees imaging survey of the Large Magellanic Cloud (LMC), using the IRAC and MIPS instruments on board the Spitzer Space Telescope in order to survey the agents of a galaxy's evolution (SAGE), the interstellar medium (ISM) and stars in the LMC. The detection of diffuse ISM with column densities >1.2x10^21 H cm^-2 permits detailed studies of dust processes in the ISM. SAGE's point source sensitivity enables a complete census of newly formed stars with masses >3 solar masses that will determine the current star formation rate in the LMC. SAGE's detection of evolved stars with mass loss rates >1x10^-8 solar masses per year will quantify the rate at which evolved stars inject mass into the ISM of the LMC. The observing strategy includes two epochs in 2005, separated by three months, that both mitigate instrumental artifacts and constrain source variability. The SAGE data are non-proprietary. The data processing includes IRAC and MIPS pipelines and a database for mining the point source catalogs, which will be released to the community in support of Spitzer proposal cycles 4 and 5. We present initial results on the epoch 1 data with a special focus on the N79 and N83 region. The SAGE epoch 1 point source catalog has ~4 million sources. The point source counts are highest for the IRAC 3.6 microns band and decrease dramatically towards longer wavelengths consistent with the fact that stars dominate the point source catalogs and that the dusty objects, e.g. young stellar objects and dusty evolved stars that detected at the longer wavelengths, are rare in comparison. We outline a strategy for identifying foreground MW stars, that may comprise as much as 18% of the source list, and background galaxies, that may comprise ~12% of the source list.

preprint2006arXiv

The Survey for Ionization in Neutral Gas Galaxies- II. The Star Formation Rate Density of the Local Universe

We derive observed Halpha and R band luminosity densities of an HI-selected sample of nearby galaxies using the SINGG sample to be l_Halpha' = (9.4 +/- 1.8)e38 h_70 erg s^-1 Mpc^-3 for Halpha and l_R' = (4.4 +/- 0.7)e37 h_70 erg s^-1 A^-1 Mpc^-3 in the R band. This R band luminosity density is approximately 70% of that found by the Sloan Digital Sky Survey. This leads to a local star formation rate density of log(SFRD) = -1.80 +0.13/-0.07(random) +/- 0.03(systematic) + log(h_70) after applying a mean internal extinction correction of 0.82 magnitudes. The gas cycling time of this sample is found to be t_gas = 7.5 +1.3/-2.1 Gyr, and the volume-averaged equivalent width of the SINGG galaxies is EW(Halpha) = 28.8 +7.2/-4.7 A (21.2 +4.2/-3.5 A without internal dust correction). As with similar surveys, these results imply that SFRD(z) decreases drastically from z ~ 1.5 to the present. A comparison of the dynamical masses of the SINGG galaxies evaluated at their optical limits with their stellar and HI masses shows significant evidence of downsizing: the most massive galaxies have a larger fraction of their mass locked up in stars compared with HI, while the opposite is true for less massive galaxies. We show that the application of the Kennicutt star formation law to a galaxy having the median orbital time at the optical limit of this sample results in a star formation rate decay with cosmic time similar to that given by the SFRD(z) evolution. This implies that the SFRD(z) evolution is primarily due to the secular evolution of galaxies, rather than interactions or mergers. This is consistent with the morphologies predominantly seen in the SINGG sample.

preprint2006arXiv

The Survey for Ionization in Neutral Gas Galaxies: I. Description and Initial Results

We introduce the Survey for Ionization in Neutral Gas Galaxies (SINGG), a census of star formation in HI-selected galaxies. The survey consists of H-alpha and R-band imaging of a sample of 468 galaxies selected from the HI Parkes All Sky Survey (HIPASS). The sample spans three decades in HI mass and is free of many of the biases that affect other star forming galaxy samples. We present the criteria for sample selection, list the entire sample, discuss our observational techniques, and describe the data reduction and calibration methods. This paper focuses on 93 SINGG targets whose observations have been fully reduced and analyzed to date. The majority of these show a single Emission Line Galaxy (ELG). We see multiple ELGs in 13 fields, with up to four ELGs in a single field. All of the targets in this sample are detected in H-alpha indicating that dormant (non-star forming) galaxies with M(HI) > ~3e7 M_sun are very rare. A database of the measured global properties of the ELGs is presented. The ELG sample spans four orders of magnitude in luminosity (H-alpha and R-band), and H-alpha surface brightness, nearly three orders of magnitude in R surface brightness and nearly two orders of magnitude in H-alpha equivalent width (EW). The surface brightness distribution of our sample is broader than that of the Sloan Digital Sky Survey spectroscopic sample, the (EW) distribution is broader than prism-selected samples, and the morphologies found include all common types of star forming galaxies (e.g. irregular, spiral, blue compact dwarf, starbursts, merging and colliding systems, and even residual star formation in S0 and Sa spirals). (abridged)

preprint2005arXiv

Statistical Confirmation of a Stellar Upper Mass Limit

We derive the expectation value for the maximum stellar mass (m_max) in an ensemble of N stars, as a function of the IMF upper-mass cutoff (m_up) and N. We statistically demonstrate that the upper IMF of the local massive star census observed thus far in the Milky Way and Magellanic Clouds clearly exhibits a universal upper mass cutoff around 120 - 200 M_sun for a Salpeter IMF, although the result is more ambiguous for a steeper IMF.

preprint2003arXiv

Massive Field Stars and the Stellar Clustering Law

The distribution of N*, the number of OB stars per association or cluster, appears to follow a universal power-law form $N*^{-2}$ in the local Universe. We evaluate the distribution of N* in the Small Magellanic Cloud using recent broadband optical and space-ultraviolet data, with special attention to the lowest values of N*. We find that the power-law distribution in N* continues smoothly down to N*=1. This strongly suggests that the formation of field massive stars is a continuous process with those in associations, and that the field stars do not originate from a different star formation mode. Our results are consistent with the model that field massive stars represent the most massive members in groups of smaller stars, as expected if the clustering law applies to much lower masses as is expected from the stellar initial mass function (IMF). These results are consistent with the simultaneous existence of a universal IMF and a universal clustering law. Jointly, these laws imply that the fraction of field OB stars typically ranges from about 35% to 7% for most astrophysical situations, with an inverse logarithmic dependence on the most populous cluster, and hence, on galaxy size and/or star formation rate. There are important consequences for global feedback effects in galaxies: field stars should therefore contribute proportionately to the volume of the warm ionized medium, and equal relative contributions by superbubbles of all sizes to the interstellar porosity are expected.

preprint1997arXiv

Comparison of H II region luminosities with observed stellar ionizing sources in the Large Magellanic Cloud

We estimate the total predicted Lyc emission rates of OB associations for which the complete census of O star spectral types exists. The results are compared to the observed H-alpha luminosities of the host H II regions. We find evidence for substantial leakage of ionizing photons from some H II regions, while others appear to be radiation bounded. We estimate that overall for the LMC, 0-51% of the ionizing radiation escapes the local nebulae, and would be available to ionize the diffuse, warm, ionized medium (WIM) in that galaxy. This range of values is consistent with the observed 35% fraction of \Ha\ luminosity emitted by the WIM in the LMC, as well as the corresponding fractions observed in other nearby galaxies. It is therefore possible that photoionization by O stars is indeed the dominant ionization mechanism for the WIM.

preprint1997arXiv

The Superbubble Size Distribution in the Interstellar Medium of Galaxies

We use the standard, adiabatic shell evolution to predict the size distribution N(R) for populations of OB superbubbles in a uniform ISM. We derive N(R) for simple cases of superbubble creation rate and mechanical luminosity function (MLF). For R < the characteristic radius R_e, N(R) is dominated by stalled objects, while for R>R_e it is dominated by growing objects. We also briefly investigate N(R) resulting from momentum-conserving shell evolution. We predict a peak in N(R) corresponding to individual SNRs. To estimate the MLF, we also examine evolutionary effects on the HII region luminosity function (HII LF), finding that for nebular luminosity fading as a power law in time, there is a minimum observed slope for the HII LFs. Comparison with the largely complete HI hole catalog for the SMC shows surprising agreement in the predicted and observed slope of N(R), suggesting that no other fundamental process is needed to explain the size distribution of shells in the SMC. Further comparison with largely incomplete HI data for M31, M33, and Holmberg II is also encouraging. We present expressions for the ISM porosity parameters, and estimate that they are substantially <1 for all of the galaxies except Holmberg II. Most of these galaxies therefore may not be strongly dominated by a hot interstellar component. However, porosity results for the Galaxy remain inconclusive.