Source author record

Maria R. Drout

Maria R. Drout 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

26works
11topics
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

26 published item(s)

preprint2022arXiv

A low-cost ultraviolet-to-infrared absolute quantum efficiency characterization system of detectors

We present a low-cost ultraviolet to infrared absolute quantum efficiency detector characterization system developed using commercial off-the-shelf components. The key components of the experiment include a light source,a regulated power supply, a monochromator, an integrating sphere, and a calibrated photodiode. We provide a step-by-step procedure to construct the photon and quantum efficiency transfer curves of imaging sensors. We present results for the GSENSE 2020 BSI CMOS sensor and the Sony IMX 455 BSI CMOS sensor. As a reference for similar characterizations, we provide a list of parts and associated costs along with images of our setup.

preprint2022arXiv

Discovery of a Short-Period and Unusually Helium-Deficient Dwarf Nova KSP-OT-201701a by the KMTNet Supernova Program

We present the first ever discovery of a short-period and unusually helium-deficient dwarf nova KSP-OT-201701a by the Korea Microlensing Telescope Network Supernova Program. The source shows three superoutbursts, each led by a precursor outburst, and several normal outbursts in BVI during the span of ~2.6 years with supercycle and normal cycle lengths of about 360 and 76 days, respectively. Spectroscopic observations near the end of a superoutburst reveal the presence of strong double-peaked HI emission lines together with weak HeI emission lines. The helium-to-hydrogen intensity ratios measured by HeIλ5876 and Hα lines are 0.10 {\pm} 0.01 at a quiescent phase and 0.26 {\pm} 0.04 at an outburst phase, similar to the ratios found in long-period dwarf novae while significantly lower than those in helium cataclysmic variables (He CVs). Its orbital period of 51.91 {\pm} 2.50 minutes, which is estimated based on time series spectroscopy, is a bit shorter than the superhump period of 56.52 {\pm} 0.19 minutes, as expected from the gravitational interaction between the eccentric disk and the secondary star. We measure its mass ratio to be 0.37^{+0.32}_{-0.21} using the superhump period excess of 0.089 {\pm} 0.053. The short orbital period, which is under the period minimum, the unusual helium deficiency, and the large mass ratio suggest that KSP-OT-201701a is a transition object evolving to a He CV from a long-period dwarf nova with an evolved secondary star.

preprint2022arXiv

Infant-phase reddening by surface Fe-peak elements in a normal Type Ia Supernova

Type Ia Supernovae are thermonuclear explosions of white dwarf stars. They play a central role in the chemical evolution of the Universe and are an important measure of cosmological distances. However, outstanding questions remain about their origins. Despite extensive efforts to obtain natal information from their earliest signals, observations have thus far failed to identify how the majority of them explode. Here, we present infant-phase detections of SN 2018aoz from a brightness of -10.5 absolute AB magnitudes -- the lowest luminosity early Type Ia signals ever detected -- revealing a hitherto unseen plateau in the $B$-band that results in a rapid redward color evolution between 1.0 and 12.4 hours after the estimated epoch of first light. The missing $B$-band flux is best-explained by line-blanket absorption from Fe-peak elements in the outer 1% of the ejected mass. The observed $B-V$ color evolution of the SN also matches the prediction from an over-density of Fe-peak elements in the same outer 1% of the ejected mass, whereas bluer colors are expected from a purely monotonic distribution of Fe-peak elements. The presence of excess nucleosynthetic material in the extreme outer layers of the ejecta points to enhanced surface nuclear burning or extended sub-sonic mixing processes in some normal Type Ia Supernova explosions.

preprint2021arXiv

An Early-Time Optical and Ultraviolet Excess in the type-Ic SN 2020oi

We present photometric and spectroscopic observations of Supernova 2020oi (SN 2020oi), a nearby ($\sim$17 Mpc) type-Ic supernova (SN Ic) within the grand-design spiral M100. We undertake a comprehensive analysis to characterize the evolution of SN 2020oi and constrain its progenitor system. We detect flux in excess of the fireball rise model $δt \approx 2.5$ days from the date of explosion in multi-band optical and UV photometry from the Las Cumbres Observatory and the Neil Gehrels Swift Observatory, respectively. The derived SN bolometric luminosity is consistent with an explosion with $M_{\rm ej} = 0.81 \pm 0.03 M_{\odot}$, $E_{k}= 0.79 \pm 0.09 \times 10^{51} \rm{erg} \rm{s}^{-1}$, and $M_{\rm Ni56} = 0.08 \pm 0.02 M_{\odot}$. Inspection of the event's decline reveals the highest $Δm_{15,\rm{bol}}$ reported for a stripped-envelope event to date. Modeling of optical spectra near event peak indicates a partially mixed ejecta comparable in composition to the ejecta observed in SN 1994I, while the earliest spectrum shows signatures of a possible interaction with material of a distinct composition surrounding the SN progenitor. Further, Hubble Space Telescope (HST) pre-explosion imaging reveals a stellar cluster coincident with the event. From the cluster photometry, we derive the mass and age of the SN progenitor using stellar evolution models implemented in the BPASS library. Our results indicate that SN 2020oi occurred in a binary system from a progenitor of mass $M_{\rm ZAMS} \approx 9.5 \pm 1.0 M_{\odot}$, corresponding to an age of $27 \pm 7$ Myr. SN 2020oi is the dimmest SN Ic event to date for which an early-time flux excess has been observed, and the first in which an early excess is unlikely to be associated with shock-cooling.

preprint2021arXiv

Final Moments I: Precursor Emission, Envelope Inflation, and Enhanced Mass loss Preceding the Luminous Type II Supernova 2020tlf

We present panchromatic observations and modeling of supernova (SN) 2020tlf, the first normal type II-P/L SN with confirmed precursor emission, as detected by the Young Supernova Experiment transient survey with the Pan-STARRS1 telescope. Pre-explosion emission was detected in $riz-$bands at 130 days prior to SN 2020tlf and persisted at relatively constant flux until first light. Soon after discovery, "flash" spectroscopy of SN 2020tlf revealed prominent narrow symmetric emission lines ($v_w < 300$ km s$^{-1}$) that resulted from the photo-ionization of unshocked circumstellar material (CSM) shedded in progenitor mass loss episodes in the final weeks to months before explosion. Surprisingly, this novel display of pre-SN emission and associated mass loss occurred in a RSG progenitor with ZAMS mass of only 10-12 M$_{\odot}$, as inferred from nebular spectra. Modeling of the light curve and multi-epoch spectra with the non-LTE radiative transfer code CMFGEN and radiation-hydrodynamical (RHD) code HERACLES suggests a dense CSM limited to $r \approx 10^{15}$ cm, and mass loss rate of $10^{-2}$ M$_{\odot}$ yr$^{-1}$. The subsequent luminous light-curve plateau and persistent blue excess indicates an extended progenitor, compatible with a RSG model with $R_{\star} = 1100$ R$_{\odot}$. Inferences from the limits on the shock-powered X-ray and radio luminosity are consistent with these conclusions and suggest a CSM density of $ρ< 2 \times 10^{-16}$ g cm$^{-3}$ for distances of $r \approx 5 \times 10^{15}$ cm, as well as a mass loss rate of $\dot M<1.3 \times 10^{-5}\,\rm{M_{\odot}\,yr^{-1}}$ at larger distances. A promising power source for the observed precursor emission is the ejection of stellar material following energy disposition into the stellar envelope as a result of gravity waves emitted during either neon/oxygen burning or a nuclear flash from silicon combustion.

preprint2021arXiv

FIRST J153350.8+272729: the radio afterglow of a decades-old tidal disruption event

We present the discovery of the fading radio transient FIRST J153350.8+272729. The source had a maximum observed 5-GHz radio luminosity of $8\times10^{39}$ erg s$^{-1}$ in 1986, but by 2019 had faded by a factor of nearly 400. It is located 0.15 arcsec from the center of a galaxy (SDSS J153350.89+272729) at 147 Mpc, which shows weak Type II Seyfert activity. We show that a tidal disruption event (TDE) is the preferred scenario for FIRST J153350.8+272729, although it could plausibly be interpreted as the afterglow of a long-duration gamma-ray burst. This is only the second TDE candidate to be first discovered at radio wavelengths. Its luminosity fills a gap between the radio afterglows of sub-relativistic TDEs in the local universe, and relativistic TDEs at high redshifts. The unusual properties of FIRST J153350.8+272729 (ongoing nuclear activity in the host galaxy, high radio luminosity) motivate more extensive TDE searches in untargeted radio surveys.

preprint2021arXiv

Gaia EDR3 confirms that Westerlund 1 is closer and older than previously thought

Using {\it Gaia} Early Data Release 3 (EDR3) parallaxes and Bayesian inference, we infer a parallax of the Westerlund 1 (Wd1) cluster. We find a parallax of $0.34\pm{0.05}$ mas corresponding to a distance of $2.8^{+0.7}_{-0.6}$ kpc. The new {\it Gaia} EDR3 distance is consistent with our previous result using {\it Gaia} DR2 parallaxes. This confirms that Wd1 is less massive and older than previously assumed. Compared to DR2, the EDR3 individual parallax uncertainties for each star decreased by 30\%. However, the aggregate parallax uncertainty for the cluster remained the same. This suggests that the uncertainty is dominated by systematics, which is possibly due to crowding, motions within the cluster, or motions due to binary orbits.

preprint2021arXiv

Revised Stellar Parameters for V471 Tau, A Post-common Envelope Binary in the Hyades

V471 Tau is a post-common-envelope binary consisting of an eclipsing DA white dwarf and a K-type main-sequence star in the Hyades star cluster. We analyzed publicly available photometry and spectroscopy of V471 Tau to revise the stellar and orbital parameters of the system. We used archival K2 photometry, archival Hubble Space Telescope spectroscopy, and published radial-velocity measurements of the K-type star. Employing Gaussian processes to fit for rotational modulation of the system flux by the main-sequence star, we recovered the transits of the white dwarf in front of the main-sequence star for the first time. The transits are shallower than would be expected from purely geometric occultations owing to gravitational microlensing during transit, which places an additional constraint on the white-dwarf mass. Our revised mass and radius for the main-sequence star is consistent with single-star evolutionary models given the age and metallicity of the Hyades. However, as noted previously in the literature, the white dwarf is too massive and too hot to be the result of single-star evolution given the age of the Hyades, and may be the product of a merger scenario. We independently estimate the conditions of the system at the time of common envelope that would result in the measured orbital parameters today.

preprint2021arXiv

Searching for Anomalies in the ZTF Catalog of Periodic Variable Stars

Periodic variables illuminate the physical processes of stars throughout their lifetime. Wide-field surveys continue to increase our discovery rates of periodic variable stars. Automated approaches are essential to identify interesting periodic variable stars for multi-wavelength and spectroscopic follow-up. Here, we present a novel unsupervised machine learning approach to hunt for anomalous periodic variables using phase-folded light curves presented in the Zwicky Transient Facility Catalogue of Periodic Variable Stars by \citet{Chen_2020}. We use a convolutional variational autoencoder to learn a low dimensional latent representation, and we search for anomalies within this latent dimension via an isolation forest. We identify anomalies with irregular variability. Most of the top anomalies are likely highly variable Red Giants or Asymptotic Giant Branch stars concentrated in the Milky Way galactic disk; a fraction of the identified anomalies are more consistent with Young Stellar Objects. Detailed spectroscopic follow-up observations are encouraged to reveal the nature of these anomalies.

preprint2020arXiv

A Deep CFHT Optical Search for a Counterpart to the Possible Neutron Star -- Black Hole Merger GW190814

We present a wide-field optical imaging search for electromagnetic counterparts to the likely neutron star - black hole (NS-BH) merger GW190814/S190814bv. This compact binary merger was detected through gravitational waves by the LIGO/Virgo interferometers, with masses suggestive of a NS-BH merger. We imaged the LIGO/Virgo localization region using the MegaCam instrument on the Canada-France-Hawaii Telescope. We describe our hybrid observing strategy of both tiling and galaxy-targeted observations, as well as our image differencing and transient detection pipeline. Our observing campaign produced some of the deepest multi-band images of the region between 1.7 and 8.7 days post-merger, reaching a 5sigma depth of g > 22.8 (AB mag) at 1.7 days and i > 23.1 and i > 23.9 at 3.7 and 8.7 days, respectively. These observations cover a mean total integrated probability of 67.0% of the localization region. We find no compelling candidate transient counterparts to this merger in our images, which suggests that either the lighter object was tidally disrupted inside of the BH's innermost stable circular orbit, the transient lies outside of the observed sky footprint, or the lighter object is a low-mass BH. We use 5sigma source detection upper limits from our images in the NS-BH interpretation of this merger to constrain the mass of the kilonova ejecta to be Mej < 0.015Msun for a 'blue' (kappa = 0.5 cm^2 g^-1) kilonova, and Mej < 0.04Msun for a 'red' (kappa = 5-10 cm^2 g^-1) kilonova. Our observations emphasize the key role of large-aperture telescopes and wide-field imagers such as CFHT MegaCam in enabling deep searches for electromagnetic counterparts to gravitational wave events.

preprint2020arXiv

A Demonstration of Extremely Low Latency $γ$-ray, X-Ray & UV Follow-Up of a Millisecond Radio Transient

We report results of a novel high-energy follow-up observation of a potential Fast Radio Burst. The radio burst was detected by VLA/realfast and followed-up by the Neil Gehrels Swift Observatory in very low latency utilizing new operational capabilities of Swift (arXiv:2005.01751), with pointed soft X-ray and UV observations beginning at T0+32 minutes, and hard X-ray/gamma-ray event data saved around T0. These observations are $>10$x faster than previous X-ray/UV follow-up of any radio transient to date. No emission is seen coincident with the FRB candidate at T0, with a 0.2s fluence $5σ$ upper limit of $1.35\times10^{-8}$ erg cm$^{-2}$ (14-195 keV) for a SGR 1935+2154-like flare, nor at T0+32 minutes down to $3σ$ upper limits of 22.18 AB mag in UVOT u band, and $3.33\times10^{-13}$ erg cm$^{-2}$ s$^{-1}$ from 0.3-10 keV for the 2 ks observation. The candidate FRB alone is not significant enough to be considered astrophysical, so this note serves as a technical demonstration. These new Swift operational capabilities will allow future FRB detections to be followed up with Swift at even lower latencies than demonstrated here: 15-20 minutes should be regularly achievable, and 5-10 minutes occasionally achievable. We encourage FRB detecting facilities to release alerts in low latency to enable this science.

preprint2020arXiv

Inferring the parallax of Westerlund 1 from Gaia DR2

Westerlund 1 (Wd1) is potentially the largest star cluster in the Galaxy. That designation critically depends upon the distance to the cluster, yet the cluster is highly obscured, making luminosity-based distance estimates difficult. Using {\it Gaia} Data Release 2 (DR2) parallaxes and Bayesian inference, we infer a parallax of $0.35^{+0.07}_{-0.06}$ mas corresponding to a distance of $2.6^{+0.6}_{-0.4}$ kpc. To leverage the combined statistics of all stars in the direction of Wd1, we derive the Bayesian model for a cluster of stars hidden among Galactic field stars; this model includes the parallax zero-point. Previous estimates for the distance to Wd1 ranged from 1.0 to 5.5 kpc, although values around 5 kpc have usually been adopted. The {\it Gaia} DR2 parallaxes reduce the uncertainty from a factor of 3 to 18\% and rules out the most often quoted value of 5 kpc with 99\% confidence. This new distance allows for more accurate mass and age determinations for the stars in Wd1. For example, the previously inferred initial mass at the main-sequence turn-off was around 40 M$_{\odot}$; the new {\it Gaia} DR2 distance shifts this down to about 22 M$_{\odot}$. This has important implications for our understanding of the late stages of stellar evolution, including the initial mass of the magnetar and the LBV in Wd1. Similarly, the new distance suggests that the total cluster mass is about four times lower than previously calculated.

preprint2020arXiv

The $R$-Process Alliance: Fourth Data Release from the Search for $r$-Process-Enhanced Stars in the Galactic Halo

This compilation is the fourth data release from the $R$-Process Alliance (RPA) search for $r$-process-enhanced stars, and the second release based on "snapshot" high-resolution ($R \sim 30,000$) spectra collected with the du Pont 2.5m Telescope. In this data release, we propose a new delineation between the $r$-I and $r$-II stellar classes at $\mathrm{[Eu/Fe]} = +0.7$, instead of the empirically chosen $\mathrm{[Eu/Fe]} = +1.0$ level previously in use, based on statistical tests of the complete set of RPA data released to date. We also statistically justify the minimum level of [Eu/Fe] for definition of the $r$-I stars, [Eu/Fe] $> +0.3$. Redefining the separation between $r$-I and $r$-II stars will aid in analysis of the possible progenitors of these two classes of stars and whether these signatures arise from separate astrophysical sources at all. Applying this redefinition to previous RPA data, the number of identified $r$-II and $r$-I stars changes to 51 and 121, respectively, from the initial set of data releases published thus far. In this data release, we identify 21 new $r$-II, 111 new $r$-I (plus three re-identified), and 7 new (plus one re-identified) limited-$r$ stars out of a total of 232 target stars, resulting in a total sample of 72 new $r$-II stars, 232 new $r$-I stars, and 42 new limited-$r$ stars identified by the RPA to date.

preprint2020arXiv

The Red Supergiant Binary Fraction of the Large Magellanic Cloud

The binary fraction of unevolved massive stars is thought to be 70-100% but there are few observational constraints on the binary fraction of the evolved version of a subset of these stars, the red supergiants (RSGs). Here we identify a complete sample of RSGs in the Large Magellanic Cloud (LMC) using new spectroscopic observations and archival UV, IR and broadband optical photometry. We find 4090 RSGs with log L/Lo > 3.5 with 1820 of them having log L/Lo > 4, which we believe is our completeness limit. We additionally spectroscopically confirmed 38 new RSG+B star binaries in the LMC, bringing the total known up to 55. We then estimated the binary fraction using a k-nearest neighbors algorithm that classifies stars as single or binary based on photometry with a spectroscopic sample as a training set. We take into account observational biases such as line-of-sight stars and binaries in eclipse while also calculating model-dependent corrections for RSGs with companions that our observations were not designed to detect. Based on our data, we find an initial result of 13.5 +7.56/-6.67% for RSGs with O or B-type companions. Using the Binary Population and Spectral Synthesis (BPASS) models to correct for unobserved systems, this corresponds to a total RSG binary fraction of 19.5 +7.6/-6.7%. This number is in broad agreement with what we would expect given an initial OB binary distribution of 70%, a predicted merger fraction of 20-30% and a binary interaction fraction of 40-50%.

preprint2020arXiv

The Red Supergiant Content of M31 and M33

We identify red supergiants (RSGs) in our spiral neighbors M31 and M33 using near-IR (NIR) photometry complete to a luminosity limit of log L/Lo=4.0. Our archival survey data cover 5 deg^2 of M31, and 3 deg^2 for M33, and are likely spatially complete for these massive stars. Gaia is used to remove foreground stars, after which the RSGs can be separated from asymptotic giant branch (AGB) stars in the color-magnitude diagram. The photometry is used to derive effective temperatures and bolometric luminosities via MARCS stellar atmosphere models. The resulting H-R diagrams show superb agreement with the evolutionary tracks of the Geneva evolutionary group. Our census includes 6400 RSGs in M31 and 2850 RSGs in M33 within their Holmberg radii; by contrast, only a few hundred RSGs are known so far in the Milky Way. Our catalog serves as the basis for a study of the RSG binary frequency being published separately, as well as future studies relating to the evolution of massive stars. Here we use the matches between the NIR-selected RSGs and their optical counterparts to show that the apparent similarity in the reddening of OB stars in M31 and M33 is the result of Malmquist bias; the average extinction in M31 is likely higher than that of M33. As expected, the distribution of RSGs follows that of the spiral arms, while the much older AGB population is more uniformly spread across each galaxy's disk.

preprint2020arXiv

Thirty Years of Radio Observations of Type Ia SN 1972E and SN 1895B: Constraints on Circumstellar Shells

We have imaged over 35 years of archival Very Large Array (VLA) observations of the nearby (d$_{\rm{L}}$ $=$ 3.15 Mpc) Type Ia supernovae SN\,1972E and SN\,1895B between 9 and 121 years post-explosion. No radio emission is detected, constraining the 8.5 GHz luminosities of SN\,1972E and SN\,1895B to be L$_{ν,8.5\rm{GHz}}$ $<$ 6.0 $\times$ 10$^{23}$ erg s$^{-1}$ Hz$^{-1}$ 45 years post-explosion and L$_{ν,8.5\rm{GHz}}$ $<$ 8.9 $\times$ 10$^{23}$ erg s$^{-1}$ Hz$^{-1}$ 121 years post-explosion, respectively. These limits imply a clean circumstellar medium (CSM), with $n$ $<$ 0.9 cm$^{-3}$ out to radii of a few $\times$ 10$^{18}$ cm, if the SN blastwave is expanding into uniform density material. Due to the extensive time coverage of our observations, we also constrain the presence of CSM shells surrounding the progenitor of SN\,1972E. We rule out essentially all medium and thick shells with masses of 0.05$-$0.3 M$_\odot$ at radii between $\sim$10$^{17}$ and 10$^{18}$ cm, and thin shells at specific radii with masses down to $\lesssim$0.01 M$_\odot$. These constraints rule out swaths of parameter space for a range of single and double degenerate progenitor scenarios, including recurrent nova, core-degenerate objects, ultra-prompt explosions and white dwarf (WD) mergers with delays of a few hundred years between the onset of merger and explosion. Allowed progenitors include WD-WD systems with a significant ($>$ 10$^{4}$ years) delay from the last episode of common envelope evolution and single degenerate systems undergoing recurrent nova, provided that the recurrence timescale i short and the system has been in the nova phase for $\gtrsim$10$^{4}$ yr, such that a large ($>$ 10$^{18}$ cm) cavity has been evacuated. Future multi-epoch observations of additional intermediate-aged Type Ia SNe will provide a comprehensive view of the large-scale CSM environments around these explosions.

preprint2020arXiv

Variability of massive stars in M31 from the Palomar Transient Factory

Using data from the (intermediate) Palomar Transient Factory (iPTF), we characterize the time variability of ~500 massive stars in M31. Our sample is those stars which are spectrally typed by Massey and collaborators, including Luminous Blue Variables, Wolf-Rayets, and warm and cool supergiants. We use the high-cadence, long-baseline (~5 years) data from the iPTF survey, coupled with data-processing tools that model complex features in the light curves. We find widespread photometric (R-band) variability in the upper Hertzsprung Russell diagram (or CMD) with an increasing prevalence of variability with later spectral types. Red stars (V-I>1.5) exhibit larger amplitude fluctuations than their bluer counterparts. We extract a characteristic variability timescale, tch, via wavelet transformations that are sensitive to both continuous and localized fluctuations. Cool supergiants are characterized by longer timescales (>100 days) than the hotter stars. The latter have typical timescales of tens of days but cover a wider range, from our resolution limit of a few days to longer than 100 days timescales. Using a 60-night block of data straddling two nights with a cadence of around 2 minutes, we extracted tch in the range 0.1--10 days with amplitudes of a few percent for 13 stars. Though there is broad agreement between the observed variability characteristics in the different parts of the upper CMD with theoretical predictions, detailed comparison requires models with a more comprehensive treatment of the various physical processes operating in these stars such as pulsation, subsurface convection, and the effect of binary companions.

preprint2019arXiv

LRP2020: The cosmic origin and evolution of the elements

The origin of many elements of the periodic table remains an unsolved problem. While many nucleosynthetic channels are broadly understood, significant uncertainties remain regarding certain groups of elements such as the intermediate and rapid neutron-capture processes, the p-process, or the origin of odd-Z elements in the most metal-poor stars. Canada has a long tradition of leadership in nuclear astrophysics, dating back to the work of Alastair Cameron in the 1950s. Recent faculty hires have further boosted activity in the field, including transient observation and theory, survey science on galactic nucleosynthesis, and nuclear experiments. This white paper contains a brief overview of recent activity in the community, highlighting strengths in each sub-field, and provides recommendations to improve interdisciplinary collaboration. Sustaining Canadian leadership in the next decade will require, on the observational side, access to transient and non-transient surveys like LSST, SKA, or MSE, support for target-of-opportunity observing in current and future Canadian telescopes, and participation in next-generation X-ray telescopes such as ATHENA. State-of-the-art theoretical predictions will require an ambitious succession plan for the Niagara supercomputer to support large parallel jobs. We propose a funding instrument for postdoctoral training that reflects the interdisciplinary nature of nuclear astrophysics research, and the creation of a national collaborative funding program that allows for joint projects and workshop organization.

preprint2019arXiv

The Luminosity Function of Red Supergiants in M31

The mass-loss rates of red supergiant stars (RSGs) are poorly constrained by direct measurements, and yet the subsequent evolution of these stars depends critically on how much mass is lost during the RSG phase. In 2012 the Geneva evolutionary group updated their mass-loss prescription for RSGs with the result that a 20 solar mass star now loses 10x more mass during the RSG phase than in the older models. Thus, higher mass RSGs evolve back through a second yellow supergiant phase rather than exploding as Type II-P supernovae, in accord with recent observations (the so-called "RSG Problem"). Still, even much larger mass-loss rates during the RSG phase cannot be ruled out by direct measurements of their current dust-production rates, as such mass-loss is episodic. Here we test the models by deriving a luminosity function for RSGs in the nearby spiral galaxy M31 which is sensitive to the total mass loss during the RSG phase. We carefully separate RSGs from asymptotic giant branch stars in the color-magnitude diagram following the recent method exploited by Yang and collaborators in their Small Magellanic Cloud studies. Comparing our resulting luminosity function to that predicted by the evolutionary models shows that the new prescription for RSG mass-loss does an excellent job of matching the observations, and we can readily rule out significantly larger values.

preprint2013arXiv

GRB 130606A as a Probe of the Intergalactic Medium and the Interstellar Medium in a Star-forming Galaxy in the First Gyr After the Big Bang

We present high signal-to-noise ratio Gemini and MMT spectroscopy of the optical afterglow of the gamma-ray burst (GRB) 130606A at redshift z=5.913, discovered by Swift. This is the first high-redshift GRB afterglow to have spectra of comparable quality to those of z~6 quasars. The data exhibit a smooth continuum at near-infrared wavelengths that is sharply cut off blueward of 8410 Angs due to absorption from Ly-alpha at redshift z~5.91, with some flux transmitted through the Ly-alpha forest between 7000-7800 Angs. We use column densities inferred from metal absorption lines to constrain the metallicity of the host galaxy between a lower limit of [Si/H]>-1.7 and an upper limit of [S/H]<-0.5 set by the non-detection of S II absorption. We demonstrate consistency between the dramatic evolution in the transmission fraction of Ly-alpha seen in this spectrum over the redshift range z=4.9 to 5.85 with that previously measured from observations of high-redshift quasars. There is an extended redshift interval of Delta-z=0.12 in the Ly-alpha forest at z=5.77 with no detected transmission, leading to a 3-sigma upper limit on the mean Ly-alpha transmission fraction of <0.2% (or tau_eff(Ly-alpha) > 6.4). This is comparable to the lowest-redshift Gunn-Peterson troughs found in quasar spectra. We set a 2-sigma upper limit of 0.11 on the neutral fraction of the IGM at the redshift of the GRB from the lack of a Ly-alpha red damping wing, assuming a model with a constant neutral density. Some Ly-beta and Ly-gamma transmission is detected in this redshift window, indicating that it is not completely opaque, and hence that the IGM is nonetheless mostly ionized at these redshifts. GRB 130606A thus for the first time realizes the promise of GRBs as probes of the first galaxies and cosmic reionization.

preprint2012arXiv

A Spectroscopic Study of Type Ibc Supernova Host Galaxies from Untargeted Surveys

We present the largest spectroscopic study of the host environments of Type Ibc supernovae (SN Ibc) discovered exclusively by untargeted SN searches. Past studies of SN Ibc host environments have been biased towards high-mass, high-metallicity galaxies by focusing on SNe discovered in galaxy-targeted SN searches. Our new observations more than double the total number of spectroscopic stellar population age and metallicity measurements published for untargeted SN Ibc host environments, and extend to a median redshift about twice as large as previous statistical studies (z = 0.04). For the 12 SNe Ib and 21 SNe Ic in our metallicity sample, we find median metallicities of log(O/H)+12 = 8.48 and 8.61, respectively, but determine that the discrepancy in the full distribution of metallicities is not statistically significant. This median difference would correspond to only a small difference in the mass loss via metal-line driven winds (<30%), suggesting this does not play the dominant role in distinguishing SN Ib and Ic progenitors. However, the median metallicity of the 7 broad-lined SN Ic (SN Ic-BL) in our sample is significantly lower, log(O/H)+12 = 8.34. The age of the young stellar population of SN Ic-BL host environments also seems to be lower than for SN Ib and Ic, but our age sample is small. A synthesis of SN Ibc host environment spectroscopy to date does not reveal a significant difference in SN Ib and Ic metallicities, but reinforces the significance of the lower metallicities for SN Ic-BL. This combined sample demonstrates that galaxy-targeted SN searches introduce a significant bias for studies seeking to infer the metallicity distribution of SN progenitors, and we identify and discuss other systematic effects that play smaller roles. We discuss the path forward for making progress on SN Ibc progenitor studies in the LSST era.

preprint2012arXiv

The Yellow and Red Supergiants of M33

Yellow and red supergiants are evolved massive stars whose numbers and locations on the HR diagram can provide a stringent test for models of massive star evolution. Previous studies have found large discrepancies between the relative number of yellow supergiants observed as a function of mass and those predicted by evolutionary models, while a disagreement between the predicted and observed locations of red supergiants on the HR diagram was only recently resolved. Here we extend these studies by examining the yellow and red supergiant populations of M33. Unfortunately, identifying these stars is difficult as this portion of the color-magnitude diagram is heavily contaminated by foreground dwarfs. We identify the red supergiants through a combination of radial velocities and a two-color surface gravity discriminant and, after re-characterizing the rotation curve of M33 with our newly selected red supergiants, we identify the yellow supergiants through a combination of radial velocities and the strength of the OI $λ$7774 triplet. We examine ~1300 spectra in total and identify 121 yellow supergiants (a sample which is unbiased in luminosity above log(L/L\odot) ~ 4.8) and 189 red supergiants. After placing these objects on the HR diagram, we find that the latest generation of Geneva evolutionary tracks show excellent agreement with the observed locations of our red and yellow supergiants, the observed relative number of yellow supergiants with mass and the observed red supergiant upper mass limit. These models therefore represent a drastic improvement over previous generations.

preprint2011arXiv

The First Systematic Study of Type Ibc Supernova Multi-band Light Curves

We present detailed optical photometry for 25 Type Ibc supernovae within d\approx150 Mpc obtained with the robotic Palomar 60-inch telescope in 2004-2007. This study represents the first uniform, systematic, and statistical sample of multi-band SNe Ibc light curves available to date. We correct the light curves for host galaxy extinction using a new technique based on the photometric color evolution, namely, we show that the (V-R) color of extinction-corrected SNe Ibc at t\approx10 days after V-band maximum is tightly distributed, (V-R)=0.26+-0.06 mag. Using this technique, we find that SNe Ibc typically suffer from significant host galaxy extinction, E(B-V)\approx0.4 mag. A comparison of the extinction-corrected light curves for SNe Ib and Ic reveals that they are statistically indistinguishable, both in luminosity and decline rate. We report peak absolute magnitudes of M_R=-17.9+-0.9 mag and M_R=-18.3+-0.6 mag for SNe Ib and Ic, respectively. Focusing on the broad-lined SNe Ic, we find that they are more luminous than the normal SNe Ibc sample, M_R=-19.0+-1.1 mag, with a probability of only 1.6% that they are drawn from the same population of explosions. By comparing the peak absolute magnitudes of SNe Ic-BL with those inferred for local engine-driven explosions (GRB-SN 1998bw, XRF-SN 2006aj, and SN2009bb) we find a 25% probability that they are drawn from the SNe Ic-BL population. Finally, we fit analytic models to the light-curves to derive typical Ni-56 masses of M_Ni \approx0.2 and 0.5 M_sun for SNe Ibc and SNe Ic-BL, respectively. With reasonable assumptions for the photospheric velocities, we extract kinetic energy and ejecta mass values of M_ej \approx 2 M_sun and E_K\approx1e+51 erg for SNe Ibc, while for SNe Ic-BL we find higher values, M_ej\approx5 M_sun and E_K\approx1e+52 erg. We discuss the implications for the progenitors of SNe Ibc and their relation to engine-driven explosions [ABRIDGED].

preprint2010arXiv

Yellow Supergiants in the Small Magellanic Cloud (SMC): Putting Current Evolutionary Theory to the Test

The yellow supergiant content of nearby galaxies provides a critical test of massive star evolutionary theory. While these stars are the brightest in a galaxy, they are difficult to identify because a large number of foreground Milky Way stars have similar colors and magnitudes. We previously conducted a census of yellow supergiants within M31 and found that the evolutionary tracks predict a yellow supergiant duration an order of magnitude longer than we observed. Here we turn our attention to the SMC, where the metallicity is 10x lower than that of M31, which is important as metallicity strongly affects massive star evolution. The SMC's large radial velocity (~160 km/s) allows us to separate members from foreground stars. Observations of ~500 candidates yielded 176 near-certain SMC supergiants, 16 possible SMC supergiants, along with 306 foreground stars and provide good relative numbers of yellow supergiants down to 12Mo. Of the 176 near-certain SMC supergiants, the kinematics predicted by the Besancon model of the Milky Way suggest a foreground contamination of >4%. After placing the SMC supergiants on the H-R diagram and comparing our results to the Geneva evolutionary tracks, we find results similar to those of the M31 study: while the locations of the stars on the H-R diagram match the locations of evolutionary tracks well, the models over-predict the yellow supergiant lifetime by a factor of ten. Uncertainties about the mass-loss rates on the main-sequence thus cannot be the primary problem with the models.

preprint2009arXiv

Yellow Supergiants in the Andromeda Galaxy (M31)

The yellow supergiant content of nearby galaxies can provide a critical test of stellar evolution theory, bridging the gap between the hot, massive stars and the cool red supergiants. But, this region of the color-magnitude diagram is dominated by foreground contamination, requiring membership to somehow be determined. Fortunately, the large negative systemic velocity of M31, coupled to its high rotation rate, provides the means for separating the contaminating foreground dwarfs from the bona fide yellow supergiants within M31. Using the MMT, we obtained spectra of about 2900 stars, selected using the color and magnitude range to be yellow supergiants. Comparing the velocities to that of M31's rotation curve, we identified 54 certain, and 66 probable yellow supergiants from among the sea of foreground dwarfs. We find excellent agreement between the location of yellow supergiants in the H-R diagram and that predicted by the latest Geneva evolutionary tracks which include rotation. However, the relative number of yellow supergiants seen as a function of mass varies from that predicted by the models by a factor of more than 10, in the sense that more high mass yellow supergiants are predicted than are actually observed. Comparing the total number of yellow supergiants with masses above 20Mo with the estimated number of unevolved O stars indicates that the duration of the yellow supergiant phase is about 3000 years. This is consistent with what the 12Mo and 15Mo evolutionary tracks predict, but disagrees with the 20,000-80,000 year time scales predicted by the models for higher masses.

preprint2008arXiv

The magnetic environment in the central region of nearby galaxies

The central regions of galaxies harbor some of the most extreme physical phenomena, including dense stellar clusters, non-circular motions of molecular clouds and strong and pervasive magnetic field structures. In particular, radio observations have shown that the central few hundred parsecs of our Galaxy has a striking magnetic field configuration. It is not yet clear whether these magnetic structures are unique to our Milky Way or a common feature of all similar galaxies. Therefore, we report on (a) a new radio polarimetric survey of the central 200 pc of the Galaxy to better characterize the magnetic field structure and (b) a search for large-scale and organized magnetized structure in the nuclear regions of nearby galaxies using data from the Very Large Array (VLA) archive. The high angular resolution of the VLA allows us to study the central 1 kpc of the nearest galaxies to search for magnetized nuclear features similar to what is detected in our own Galactic center. Such magnetic features play a important role in the nuclear regions of galaxies in terms of gas transport and the physical conditions of the interstellar medium in this unusual region of galaxies.