Source author record

M. Rafelski

M. Rafelski 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

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

15 published item(s)

preprint2022arXiv

A Cautionary Tale of LyC Escape Fraction Estimates from High Redshift Galaxies

Measuring the escape fraction, $f_{\rm esc}$, of ionizing, Lyman Continuum (LyC) radiation is key to our understanding of the process of cosmic reionization. In this paper we provide a methodology for recovering the posterior probability distribution of the LyC escape fraction, $f_{\rm esc}^{\rm PDF}$, considering both the observational uncertainties and ensembles of simulated transmission functions through the intergalactic medium (IGM). We present an example of this method applied to a VUDS galaxy at $z=3.64$ and find $f_{\rm esc}^{\rm PDF}$ = 0.51$^{+0.33}_{-0.34}$ and compare this to the values computed assuming averaged IGM transmission with and without consideration of detection bias along average sightlines yielding $f_{\rm esc}^{\langle T \rangle}$ = 1.40$^{+0.80}_{-0.42}$, and $f_{\rm esc}^{\rm bias}$ = 0.82$^{+0.33}_{-0.16}$. Our results highlight the limitations of methods assuming average, smooth transmission functions. We also present MOSFIRE data for a sample of seven LyC candidates selected based on photometric redshifts at $z > 3.4$, but find that all seven have overestimated photometric redshifts by $Δz \sim 0.2$ making them unsuitable for LyC measurements. This results likely due to a bias induced by our selection criteria.

preprint2022arXiv

METAL: The Metal Evolution, Transport, and Abundance in the Large Magellanic Cloud Hubble program. III. Interstellar Depletions, Dust-to-Metal, and Dust-to-Gas Ratios Versus Metallicity

The metallicity and gas density dependence of interstellar depletions, the dust-to-gas (D/G), and dust-to-metal (D/M) ratios have important implications for how accurately we can trace the chemical enrichment of the universe; either by using FIR dust emission as a tracer of the ISM; or by using spectroscopy of damped Lyman-alpha systems (DLAs) to measure chemical abundances over a wide range of redshifts. We collect and compare large samples of depletion measurements in the Milky Way (MW), LMC (Z=0.5 Zsun), and SMC (Z=0.2 Zsun). The relation between the depletions of different elements do not strongly vary between the three galaxies, implying that abundance ratios should trace depletions accurately down to 20% solar metallicity. From the depletions, we derive D/G and D/M. The D/G increases with density, consistent with the more efficient accretion of gas-phase metals onto dust grains in the denser ISM. For log N(H) > 21 cm^-2, the depletion of metallicity tracers (S, Zn) exceeds -0.5 dex, even at 20% solar metallicity. The gas fraction of metals increases from the MW to the LMC (factor 3) and SMC (factor 6), compensating the reduction in total heavy element abundances and resulting in those three galaxies having the same neutral gas-phase metallicities. The D/G derived from depletions are a factor of 2 (LMC) and 5 (SMC) higher than the D/G derived from FIR, 21 cm, and CO emission, likely due to the combined uncertainties on the dust FIR opacity and on the depletion of carbon and oxygen.

preprint2022arXiv

The average dust attenuation curve at z~1.3 based on HST grism surveys

We present the first characterisation of the average dust attenuation curve at $z\sim1.3$ by combining rest-frame ultraviolet through near-IR photometry with Balmer decrement ($\mathrm{H}α$/$\mathrm{H}β$) constraints for $\sim$900 galaxies with $8\lesssim\log (M_\star /M_\odot)<10.2$ at $0.75<z<1.5$ in the HST WFC3 IR Spectroscopic Parallel (WISP) and 3D-HST grism surveys. Using galaxies in SDSS, we establish that the ($\mathrm{H}α$+[NII])/[OIII] line ratio and stellar mass are good proxies for the Balmer decrement in low-spectral resolution grism data when only upper-limits on $\mathrm{H}β$ are available and/or $\mathrm{H}α$ is blended with [NII]. The slope of the $z\sim1.3$ attenuation curve ($A(0.15μm)/A(V)=3.15$) and its normalization ($R_V=3.26$) lie in-between the values found for $z=0$ and $z\sim2$ dust attenuation curves derived with similar methods. These provide supporting evidence that the average dust attenuation curve of star forming galaxies evolves continuously with redshift. The $z\sim1.3$ curve has a mild 2175Å feature (bump amplitude, $E_b=0.83$; $\sim$25% that of the MW extinction curve), which is comparable to several other studies at $0<z\lesssim3$, and suggests that the average strength of this feature may not evolve significantly with redshift. The methods we develop to constrain dust attenuation from HST grism data can be applied to future grism surveys with JWST, Euclid, and RST. These new facilities will detect millions of emission line galaxies and offer the opportunity to significantly improve our understanding of how and why dust attenuation curves evolve.

preprint2021arXiv

IGM Transmission Bias for $z$ $\geq$ 2.9 Lyman Continuum Detected Galaxies

Understanding the relationship between the underlying escape fraction of Lyman continuum (LyC) photons ($f_{\rm esc}$) emitted by galaxies and measuring the distribution of observed $f_{\rm esc}$ values at high redshift is fundamental to the interpretation of the reionization process. In this paper we perform a statistical exploration of the attenuation of LyC photons by neutral hydrogen in the intergalactic medium using ensembles of simulated transmission functions. We show that LyC detected galaxies are more likely to be found in sightlines with higher-than-average transmission of LyC photons. This means that adopting a mean transmission at a given redshift leads to an overestimate of the true $f_{\rm esc}$ for LyC detected galaxies. We note, however, that mean values are appropriate for $f_{\rm esc}$ estimates of larger parent samples that include LyC non-detected galaxies. We quantify this IGM transmission bias for LyC detections in photometric and spectroscopic surveys in the literature and show that the bias is stronger for both shallower observations and for fainter parent samples (i.e. Lyman $α$ emitters versus Lyman break galaxies). We also explore the effects of varying the underlying probability distribution function (PDF) of $f_{\rm esc}$ on recovered values, showing that the underlying $f_{\rm esc}$ PDF may depend on sample selection by comparing with observational surveys. This work represents a first step in improved interpretation of LyC detections in the context of understanding $f_{\rm esc}$ from high redshift galaxies.

preprint2020arXiv

Identification of single spectral lines through supervised machine learning in a large HST survey (WISP): a pilot study for Euclid and WFIRST

Future surveys focusing on understanding the nature of dark energy (e.g., Euclid and WFIRST) will cover large fractions of the extragalactic sky in near-IR slitless spectroscopy. These surveys will detect a large number of galaxies that will have only one emission line in the covered spectral range. In order to maximize the scientific return of these missions, it is imperative that single emission lines are correctly identified. Using a supervised machine-learning approach, we classified a sample of single emission lines extracted from the WFC3 IR Spectroscopic Parallel survey (WISP), one of the closest existing analogs to future slitless surveys. Our automatic software integrates a SED fitting strategy with additional independent sources of information. We calibrated it and tested it on a "gold" sample of securely identified objects with multiple lines detected. The algorithm correctly classifies real emission lines with an accuracy of 82.6%, whereas the accuracy of the SED fitting technique alone is low (~50%) due to the limited amount of photometric data available (<=6 bands). While not specifically designed for the Euclid and WFIRST surveys, the algorithm represents an important precursor of similar algorithms to be used in these future missions.

preprint2020arXiv

Outside the Lyman-break box: detecting Lyman continuum emitters at $3.5<z<5.1$ with CLAUDS

Identifying non-contaminated sample of high-redshift galaxies with escaping Lyman continuum (LyC) flux is important for understanding the sources and evolution of cosmic reionization. We present CLAUDS $u$-band photometry of the COSMOS field to probe LyC radiation from spectroscopically confirmed galaxies at $z\geq3.5$ and outside the standard Lyman-break galaxy colour selection expectations. Complementary to the CLAUDS data, we use Subaru multi-filter photometry, \textit {Hubble Space Telescope} ($HST$) multi-filter imaging, and the spectroscopic surveys D10K, VUDS and 3D-HST. We present a sample of Lyman continuum galaxy (LCG) candidates in the redshift range $3.5\lesssim z\lesssim5.1$. Here, we introduce 5 LCG candidates, where two are flagged quality 1 and three quality 2. The estimated $f_{\rm esc}^{\rm abs}$ for quality 1 candidates are in the range $\sim5\% - 73\%$ and $\sim30\% - 93\%$. These estimates are based on our derived parameters from individual galaxies as inputs to a range of BPASS models as well as mean intergalactic medium (IGM) and maximal intergalactic and circumgalactic media (IGM+CGM) transmission. We conclude that our search for LCGs is most likely biased to lines of sight with low HI densities or free from Lyman limit systems. Our two best LCG candidates have EW (Ly$α)\leq50$Å and we find no correlation or anti-correlation between EW (Ly$α$), $f_{\rm esc}^{\rm abs}$, and $R_{\rm obs}$, the ratio of ionizing to non-ionizing observed flux in the measured passbands. Stacking candidates without solid LyC detections ($S/N<3$) results in an estimated $f_{\rm esc}^{\rm abs}$ from galaxies not greater than $1\%$.

preprint2016arXiv

Are LGRBs biased tracers of star formation? Clues from the host galaxies of the Swift/BAT6 complete sample of bright LGRBs. II: star formation rates and metallicities at z < 1

Long gamma-ray bursts (LGRBs) are associated with the deaths of massive stars and could thus be a potentially powerful tool to trace cosmic star formation. However, especially at low redshifts (z < 1.5) LGRBs seem to prefer particular types of environment. Our aim is to study the host galaxies of a complete sample of bright LGRBs to investigate the impact of the environment on GRB formation. We study host galaxy spectra of the Swift/BAT6 complete sample of 14 z < 1 bright LGRBs. We use the detected nebular emission lines to measure the dust extinction, star formation rate (SFR) and nebular metallicity (Z) of the hosts and supplement the data set with previously measured stellar masses M$_{\star}$. The distributions of the obtained properties and their interrelations (e.g. mass-metallicity and SFR-M$_{\star}$ relations) are compared to samples of field star-forming galaxies.We find that LGRB hosts at z < 1 have on average lower SFRs than if they were direct star-formation tracers. By directly comparing metallicity distributions of LGRB hosts and star-forming galaxies, we find a good match between the two populations up to Z $\sim 8.4-8.5$, after which the paucity of metal-rich LGRB hosts becomes apparent. The LGRB host galaxies of our complete sample are not inconsistent with the mass-metallicity relation at similar mean redshift and stellar masses. The cutoff against high metallicities (and high masses) can explain the low SFR values of LGRB hosts. We find a hint of increased incidence of starburst galaxies in the Swift/BAT6 z < 1 sample with respect to that of a field star-forming population. Given that the SFRs are low on average, the latter is ascribed to low stellar masses. Nevertheless the limits on the completeness and metallicity availability of current surveys, coupled with the limited number of LGRB host galaxies, prevent us from reaching more quantitative conclusions.

preprint2016arXiv

The role of quenching time in the evolution of the mass-size relation of passive galaxies from the WISP survey

We analyze how passive galaxies at z $\sim$ 1.5 populate the mass-size plane as a function of their stellar age, to understand if the observed size growth with time can be explained with the appearance of larger quenched galaxies at lower redshift. We use a sample of 32 passive galaxies extracted from the Wide Field Camera 3 Infrared Spectroscopic Parallel (WISP) survey with spectroscopic redshift 1.3 $\lesssim$ z $\lesssim$ 2.05, specific star-formation rates lower than 0.01 Gyr$^{-1}$, and stellar masses above 4.5 $\times$ 10$^{10}$ M$_\odot$. All galaxies have spectrally determined stellar ages from fitting of their rest-frame optical spectra and photometry with stellar population models. When dividing our sample into young (age $\leq$ 2.1 Gyr) and old (age $>$ 2.1 Gyr) galaxies we do not find a significant trend in the distributions of the difference between the observed radius and the one predicted by the mass-size relation. This result indicates that the relation between the galaxy age and its distance from the mass-size relation, if it exists, is rather shallow, with a slope alpha $\gtrsim$ -0.6. At face value, this finding suggests that multiple dry and/or wet minor mergers, rather than the appearance of newly quenched galaxies, are mainly responsible for the observed time evolution of the mass-size relation in passive galaxies.

preprint2015arXiv

Unveiling the Secrets of Metallicity and Massive Star Formation Using DLAs along Gamma-ray Bursts

We present the largest, publicly available, sample of Damped Lyman-$α$ systems (DLAs) along Gamma-ray Bursts (GRB) line of sights in order to investigate the environmental properties of long GRBs in the $z=1.8-6$ redshift range. Compared with the most recent quasar DLAs sample (QSO-DLA), our analysis shows that GRB-DLAs probe a more metal enriched environment at $z\gtrsim3$, up to $[X/H]\sim-0.5$. In the $z=2-3$ redshift range, despite the large number of lower limits, there are hints that the two populations may be more similar (only at 90\% significance level). Also at \hiz, the GRB-DLA average metallicity seems to decline at a shallower rate than the QSO-DLAs: GRB-DLA hosts may be polluted with metals at least as far as $\sim 2$kpc from the GRB explosion site, probably due to previous star-formation episodes and/or supernovae explosions. This shallow metallicity trend, extended now up to $z\sim5$, confirms previous results that GRB hosts are star-forming and have, on average, higher metallicity than the general QSO-DLA population. Finally, our metallicity measurements are broadly consistent with the hypothesis of two channels of GRB progenitors, one of which is mildly affected by a metallicity bias. The metallicity evolution of modeled GRB hosts agrees reasonably well with our data up to intermediate redshift, while more data are needed to constrain the models at $z\gtrsim 4$.

preprint2014arXiv

The host galaxies of X-ray selected Active Galactic Nuclei to z=2.5: Structure, star-formation and their relationships from CANDELS and Herschel/PACS

We study the relationship between the structure and star-formation rate (SFR) of X-ray selected low and moderate luminosity active galactic nuclei (AGNs) in the two Chandra Deep Fields, using Hubble Space Telescope imaging from the Cosmic Assembly Near Infrared Extragalactic Legacy Survey (CANDELS) and deep far-infrared maps from the PEP+GOODS-Herschel survey. We derive detailed distributions of structural parameters and FIR luminosities from carefully constructed control samples of galaxies, which we then compare to those of the AGNs. At z~1, AGNs show slightly diskier light profiles than massive inactive (non-AGN) galaxies, as well as modestly higher levels of gross galaxy disturbance (as measured by visual signatures of interactions and clumpy structure). In contrast, at z~2, AGNs show similar levels of galaxy disturbance as inactive galaxies, but display a red central light enhancement, which may arise due to a more pronounced bulge in AGN hosts or due to extinguished nuclear light. We undertake a number of tests of these alternatives, but our results do not strongly favour one interpretation over the other. The mean SFR and its distribution among AGNs and inactive galaxies are similar at z>1.5. At z<1, however, clear and significant enhancements are seen in the SFRs of AGNs with bulge-dominated light profiles. These trends suggest an evolution in the relation between nuclear activity and host properties with redshift, towards a minor role for mergers and interactions at z>1.5.

preprint2013arXiv

A Rapidly Star-forming Galaxy 700 Million Years After the Big Bang at z=7.51

Out of several dozen z > 7 candidate galaxies observed spectroscopically, only five have been confirmed via Lyman-alpha emission, at z=7.008, 7.045, 7.109, 7.213 and 7.215. The small fraction of confirmed galaxies may indicate that the neutral fraction in the intergalactic medium (IGM) rises quickly at z > 6.5, as Lyman-alpha is resonantly scattered by neutral gas. However, the small samples and limited depth of previous observations makes these conclusions tentative. Here we report the results of a deep near-infrared spectroscopic survey of 43 z > 6.5 galaxies. We detect only a single galaxy, confirming that some process is making Lyman-alpha difficult to detect. The detected emission line at 1.0343 um is likely to be Lyman-alpha emission, placing this galaxy at a redshift z = 7.51, an epoch 700 million years after the Big Bang. This galaxy's colors are consistent with significant metal content, implying that galaxies become enriched rapidly. We measure a surprisingly high star formation rate of 330 Msol/yr, more than a factor of 100 greater than seen in the Milky Way. Such a galaxy is unexpected in a survey of our size, suggesting that the early universe may harbor more intense sites of star-formation than expected.

preprint2013arXiv

HST/WFC3 Near-Infrared spectroscopy of quenched galaxies at z~1.5 from the WISP Survey: Stellar population properties

We combine Hubble Space Telescope (HST) G102 & G141 NIR grism spectroscopy with HST/WFC3-UVIS, HST/WFC3-IR and Spitzer/IRAC[3.6μm] photometry to assemble a sample of massive (log(M_star/M_sun) ~ 11) and quenched galaxies at z~1.5. Our sample of 41 galaxies is the largest with G102+G141 NIR spectroscopy for quenched sources at these redshifts. In contrast to the local Universe, z~1.5 quenched galaxies in the high-mass range have a wide range of stellar population properties. We find their SEDs are well fitted with exponentially decreasing SFHs, and short star-formation time-scales (τ<100Myr). Quenched galaxies also show a wide distribution in ages, between 1-4Gyr. In the (u-r)_0-versus-mass space quenched galaxies have a large spread in rest-frame color at a given mass. Most quenched galaxies populate the z~1.5 red-sequence (RS), but an important fraction of them (32%) have substantially bluer colors. Although with a large spread, we find that the quenched galaxies ON the RS have older median ages (3.1Gyr) than the quenched galaxies OFF the RS (1.5Gyr). We also show that a rejuvenated SED cannot reproduce the observed stacked spectra of (the bluer) quenched galaxies OFF the RS. We derive the upper limit on the fraction of massive galaxies ON the RS at z~1.5 to be <43%. We speculate that the young quenched galaxies OFF the RS are in a transition phase between vigorous star formation at z>2 and the z~1.5 RS. According to their estimated ages, the time required for quenched galaxies OFF the RS to join their counterparts ON the z~1.5 RS is of the order of ~1Gyr.

preprint2012arXiv

Dust extinction from Balmer decrements of star-forming galaxies at 0.75<z<1.5 with HST/WFC3 spectroscopy from the WISP survey

Spectroscopic observations of Halpha and Hbeta emission lines of 128 star-forming galaxies in the redshift range 0.75<z<1.5 are presented. These data were taken with slitless spectroscopy using the G102 and G141 grisms of the Wide-Field-Camera 3 (WFC3) on board the Hubble Space Telescope as part of the WFC3 Infrared Spectroscopic Parallel (WISP) survey. Interstellar dust extinction is measured from stacked spectra that cover the Balmer decrement (Halpha/Hbeta). We present dust extinction as a function of Halpha luminosity (down to 3 x 10^{41} erg/s), galaxy stellar mass (reaching 4 x 10^{8} Msun), and rest-frame Halpha equivalent width. The faintest galaxies are two times fainter in Halpha luminosity than galaxies previously studied at z~1.5. An evolution is observed where galaxies of the same Halpha luminosity have lower extinction at higher redshifts, whereas no evolution is found within our error bars with stellar mass. The lower Halpha luminosity galaxies in our sample are found to be consistent with no dust extinction. We find an anti-correlation of the [OIII]5007/Halpha flux ratio as a function of luminosity where galaxies with L_{Halpha}<5 x 10^{41} erg/s are brighter in [OIII]5007 than Halpha. This trend is evident even after extinction correction, suggesting that the increased [OIII]5007/Halpha ratio in low luminosity galaxies is likely due to lower metallicity and/or higher ionization parameters.

preprint2011arXiv

The Transitional Stripped-Envelope SN 2008ax: Spectral Evolution and Evidence for Large Asphericity

Supernova (SN) 2008ax in NGC 4490 was discovered within hours after shock breakout, presenting the rare opportunity to study a core-collapse SN beginning with the initial envelope-cooling phase immediately following shock breakout. We present an extensive sequence of optical and near-infrared spectra, as well as three epochs of optical spectropolarimetry. Our initial spectra, taken two days after shock breakout, are dominated by hydrogen Balmer lines at high velocity. However, by maximum light, He I lines dominated the optical and near-infrared spectra, which closely resembled those of normal Type Ib supernovae (SNe Ib) such as SN 1999ex. This spectroscopic transition defines Type IIb supernovae, but the strong similarity of SN 2008ax to normal SNe Ib beginning near maximum light, including an absorption feature near 6270A due to H-alpha at high velocities, suggests that many objects classified as SNe Ib in the literature may have ejected similar amounts of hydrogen as SN 2008ax, roughly a few x 0.01 M_sun. Early-time spectropolarimetry (6 and 9 days after shock breakout) revealed strong line polarization modulations of 3.4% across H-alpha, indicating the presence of large asphericities in the outer ejecta. The continuum shares a common polarization angle with the hydrogen, helium, and oxygen lines, while the calcium and iron absorptions are oriented at different angles. This is clear evidence of deviations from axisymmetry even in the outer ejecta. Intrinsic continuum polarization of 0.64% only nine days after shock breakout shows that the outer layers of the ejecta were quite aspherical. A single epoch of late-time spectropolarimetry, as well as the shapes of the nebular line profiles, demonstrate that asphericities extended from the outermost layers all the way down to the center of this SN. [Abridged]

preprint2009arXiv

The Quasar SDSS J1536+0441: An Unusual Double-Peaked Emitter

The quasar SDSS J153636.22+044127.0, exhibiting peculiar broad emission-line profiles with multiple components, was proposed as a candidate sub-parsec binary supermassive black hole system. More recently, imaging revealed two spatially distinct sources, leading some to suggest the system to be a quasar pair separated by ~5 kpc. We present Palomar and Keck optical spectra of this system from which we identify a third velocity component to the emission lines. We argue that the system is more likely an unusual member of the class of active galactic nuclei (AGNs) known as "double-peaked emitters" than a sub-parsec black hole binary or quasar pair. We find no significant velocity evolution of the two main peaks over the course of 0.95 yr, with a 3-sigma upper limit on any secular change of 70 km/s/yr. We also find that the three velocity components of the emission lines are spatially coincident to within 0.015" along the slit, apparently ruling out the double-quasar hypothesis.