Search

This Letter presents the detection of a source at the position of the Type Ib/c supernova (SN) 2013ge more than four years after the radioactive component is expected to have faded. This source could mark the first post-SN direct detection of a surviving companion to a stripped-envelope Type Ib/c explosion. We test this hypothesis and find the shape of the source's spectral energy distribution is most consistent with that of a B5 I supergiant. While binary models tend to predict OB-type stars for stripped-envelope companions, the location of the source on a color-magnitude diagram (CMD) places it redward of its more likely position on the main sequence (MS). The source may be temporarily out of thermal equilibrium, or a cool and inflated non-MS companion, which is similar to the suggested companion of Type Ib SN 2019yvr that was constrained from pre-SN imaging. We also consider other possible physical scenarios for the source, including a fading SN, circumstellar shock interaction, line of site coincidence, and an unresolved host star cluster, all of which will require future observations to more definitively rule out. Ultimately, the fraction of surviving companions ("binary fraction") will provide necessary constraints on binary evolution models and the underlying physics., Accepted to ApJL. 9 pages, 4 figures, 1 table

We present proper motion measurements of the oxygen-rich ejecta of the LMC supernova remnant N132D using two epochs of Hubble Space Telescope Advanced Camera for Surveys data spanning 16 years. The proper motions of 120 individual knots of oxygen-rich gas were measured and used to calculate a center of expansion (CoE) of α = 5h25m01.ˢ71 and δ = −69°38′41.″64 (J2000) with a 1σ uncertainty of 2.″90. This new CoE measurement is 9.″2 and 10.″8 from two previous CoE estimates based on the geometry of the optically emitting ejecta. We also derive an explosion age of 2770 ± 500 yr, which is consistent with recent age estimates of ≈2500 yr made from 3D ejecta reconstructions. We verified our estimates of the CoE and age using a new automated procedure that detected and tracked the proper motions of 137 knots, with 73 knots that overlap with the visually identified knots. We find that the proper motions of the ejecta are still ballistic, despite the remnant’s age, and are consistent with the notion that the ejecta are expanding into an interstellar medium cavity. Evidence for explosion asymmetry from the parent supernova is also observed. Using the visually measured proper motion measurements and corresponding CoE and age, we compare N132D to other supernova remnants with proper motion ejecta studies.

We analyze a sample of 45 Type II supernovae from the Zwicky Transient Facility (ZTF) public survey using a grid of hydrodynamical models in order to assess whether theoretically-driven forecasts can intelligently guide follow up observations supporting all-sky survey alert streams. We estimate several progenitor properties and explosion physics parameters including zero-age-main-sequence (ZAMS) mass, mass-loss rate, kinetic energy, 56Ni mass synthesized, host extinction, and the time of explosion. Using complete light curves we obtain confident characterizations for 34 events in our sample, with the inferences of the remaining 11 events limited either by poorly constraining data or the boundaries of our model grid. We also simulate real-time characterization of alert stream data by comparing our model grid to various stages of incomplete light curves (t less than 25 days, t less than 50 days, all data), and find that some parameters are more reliable indicators of true values at early epochs than others. Specifically, ZAMS mass, time of explosion, steepness parameter beta, and host extinction are reasonably constrained with incomplete light curve data, whereas mass-loss rate, kinetic energy and 56Ni mass estimates generally require complete light curves spanning greater than 100 days. We conclude that real-time modeling of transients, supported by multi-band synthetic light curves tailored to survey passbands, can be used as a powerful tool to identify critical epochs of follow up observations. Our findings are relevant to identify, prioritize, and coordinate efficient follow up of transients discovered by Vera C. Rubin Observatory., Comment: 27 pages, 14 figures, Accepted to The Astrophysical Journal

Type IIb supernovae (SNe IIb) present a unique opportunity for investigating the evolutionary channels and mechanisms governing the evolution of stripped-envelope SN progenitors due to a variety of observational constraints available. Comparison of these constraints with the full distribution of theoretical properties not only help ascertain the prevalence of observed properties in nature, but can also reveal currently unobserved populations. In this follow-up paper, we use the large grid of models presented in Sravan et al. 2019 to derive distributions of single and binary SNe IIb progenitor properties and compare them to constraints from three independent observational probes: multi-band SN light-curves, direct progenitor detections, and X-ray/radio observations. Consistent with previous work, we find that while current observations exclude single stars as SN IIb progenitors, SN IIb progenitors in binaries can account for them. We also find that the distributions indicate the existence of an unobserved dominant population of binary SNe IIb at low metallicity that arise due to mass transfer initiated on the Hertzsprung Gap. In particular, our models indicate the existence of a group of highly stripped (envelope mass ~0.1-0.2 M_sun) progenitors that are compact (, Comment: Resubmitted to the Astrophysical Journal after incorporating suggestions from the referee

We present new proper motion measurements of optically emitting oxygen-rich knots of supernova remnant 1E 0102.2-7219 (E0102), which are used to estimate the remnant's center of expansion and age. Four epochs of high resolution Hubble Space Telescope images spanning 19 yr were retrieved and analyzed. We found a robust center of expansion of alpha=1:04:02.48 and delta=-72:01:53.92 (J2000) with 1-sigma uncertainty of 1.77 arcseconds using 45 knots from images obtained with the Advanced Camera for Surveys using the F475W filter in 2003 and 2013 having the highest signal-to-noise ratio. We also estimate an upper limit explosion age of 1738 +/- 175 yr by selecting knots with the highest proper motions, that are assumed to be the least decelerated. We find evidence of an asymmetry in the proper motions of the knots as a function of position angle. We conclude that these asymmetries were most likely caused by interaction between E0102's original supernova blast wave and an inhomogeneous surrounding environment, as opposed to intrinsic explosion asymmetry. The observed non-homologous expansion suggests that the use of a free expansion model inaccurately offsets the center of expansion and leads to an overestimated explosion age. We discuss our findings as they compare to previous age and center of expansion estimates of E0102 and their relevance to a recently identified candidate central compact object., 14 pages, 11 figures, revised according to referee comments and resubmitted to ApJ

We present a three-dimensional kinematic reconstruction of the optically-emitting, oxygen-rich ejecta of supernova remnant N132D in the Large Magellanic Cloud. Data were obtained with the 6.5 m Magellan telescope in combination with the IMACS+GISMO instrument and survey [O III] $\lambda\lambda$4959,5007 line emission in a ${\sim}$3$^{\prime}~\times$ 3$^{\prime}$ region centered on N132D. The spatial and spectral resolution of our data enable detailed examination of the optical ejecta structure. The majority of N132D's optically bright oxygen ejecta are arranged in a torus-like geometry tilted approximately 28$^{\circ}$ with respect to the plane of the sky. The torus has a radius of 4.4 pc ($D_{\rm LMC}$/50 kpc), exhibits a blue-shifted radial velocity asymmetry of $-3000$ to $+2300$ km s$^{-1}$, and has a conspicuous break in its circumference. Assuming homologous expansion from the geometric center of O-rich filaments, the average expansion velocity of 1745 km s$^{-1}$ translates to an age since explosion of 2450 $\pm$ 195 yr. A faint, spatially-separated "runaway knot" (RK) with total space velocity of 3650 km s$^{-1}$ is nearly perpendicular to the torus plane and coincident with X-ray emission that is substantially enhanced in Si relative to the LMC and N132D's bulk ejecta. These kinematic and chemical signatures suggest that the RK may have had its origin deep within the progenitor star. Overall, the main shell morphology and high-velocity, Si-enriched components of N132D have remarkable similarity with that of Cassiopeia A, which was the result of a Type IIb supernova explosion. Our results underscore the need for further observations and simulations that can robustly reconcile whether the observed morphology is dominated by explosion dynamics or shaped by interaction with the environment., Comment: 18 pages, 12 figures, accepted for publication in ApJ; 3D animation of N132D can found at: http://www.physics.purdue.edu/kaboom/n132d-anim.mp4

Type IIb supernovae (SNe) are important candidates to understand mechanisms that drive the stripping of stripped-envelope (SE) supernova (SN) progenitors. While binary interactions and their high incidence are generally cited to favor them as Type IIb SN progenitors, this idea has not been tested using models covering a broad parameter space. In this paper, we use non-rotating single- and binary-star models at solar and low metallicities spanning a wide parameter space in primary mass, mass ratio, orbital period, and mass transfer efficiencies. We find that our single- and binary-star models contribute to roughly equal, however small, numbers of Type IIb SNe at solar metallicity. Binaries only dominate as progenitors at low metallicity. We also find that our models can account for less than half the observationally inferred rate for Type IIb SNe at solar metallicity, with computed rates ~, Accepted to The Astrophysical Journal

Type IIb supernovae (SNe) present a unique opportunity for understanding the progenitors of stripped-envelope (SE) SNe as the stellar progenitor of several Type IIb SNe have been identified in pre-explosion images. In this paper, we use Bayesian inference and a large grid of non-rotating solar-metallicity single and binary stellar models to derive the associated probability distributions of single and binary progenitors of the Type IIb SN 2016gkg using existing observational constraints. We find that potential binary star progenitors have smaller pre-SN hydrogen-envelope and helium-core masses than potential single-star progenitors typically by 0.1 Msun and 2 Msun, respectively. We find that, a binary companion, if present, is a main-sequence or red-giant star. Apart from this, we do not find strong constraints on the nature of the companion star. We demonstrate that the range of progenitor helium-core mass inferred from observations could help improve constraints on the progenitor. We find that the probability that the progenitor of SN 2016gkg was a binary is 22% when we use constraints only on the progenitor luminosity and effective temperature. Imposing the range of pre-SN progenitor hydrogen-envelope mass and radius inferred from SN light-curves the probability the progenitor is a binary increases to 44%. However, there is no clear preference for a binary progenitor. This is in contrast to binaries being the currently favored formation channel for Type IIb SNe. Our analysis demonstrates the importance of statistical inference methods to constrain progenitor channels., Comment: Accepted to The Astrophysical Journal Letters

We use cosmological simulations from the Feedback In Realistic Environments (FIRE) project, which implement a comprehensive set of stellar feedback processes, to study ultra-violet (UV) metal line emission from the circum-galactic medium of high-redshift (z=2-4) galaxies. Our simulations cover the halo mass range Mh ~ 2x10^11 - 8.5x10^12 Msun at z=2, representative of Lyman break galaxies. Of the transitions we analyze, the low-ionization C III (977 A) and Si III (1207 A) emission lines are the most luminous, with C IV (1548 A) and Si IV (1394 A) also showing interesting spatially-extended structures. The more massive halos are on average more UV-luminous. The UV metal line emission from galactic halos in our simulations arises primarily from collisionally ionized gas and is strongly time variable, with peak-to-trough variations of up to ~2 dex. The peaks of UV metal line luminosity correspond closely to massive and energetic mass outflow events, which follow bursts of star formation and inject sufficient energy into galactic halos to power the metal line emission. The strong time variability implies that even some relatively low-mass halos may be detectable. Conversely, flux-limited samples will be biased toward halos whose central galaxy has recently experienced a strong burst of star formation. Spatially-extended UV metal line emission around high-redshift galaxies should be detectable by current and upcoming integral field spectrographs such as the Multi Unit Spectroscopic Explorer (MUSE) on the Very Large Telescope and Keck Cosmic Web Imager (KCWI)., 16 pages, 8 figures, accepted for publication in MNRAS

We perform a Bayesian analysis of the mass distribution of stellar-mass black holes using the observed masses of 15 low-mass X-ray binary systems undergoing Roche lobe overflow and five high-mass, wind-fed X-ray binary systems. Using Markov Chain Monte Carlo calculations, we model the mass distribution both parametrically---as a power law, exponential, gaussian, combination of two gaussians, or log-normal distribution---and non-parametrically---as histograms with varying numbers of bins. We provide confidence bounds on the shape of the mass distribution in the context of each model and compare the models with each other by calculating their relative Bayesian evidence as supported by the measurements, taking into account the number of degrees of freedom of each model. The mass distribution of the low-mass systems is best fit by a power-law, while the distribution of the combined sample is best fit by the exponential model. We examine the existence of a "gap" between the most massive neutron stars and the least massive black holes by considering the value, M_1%, of the 1% quantile from each black hole mass distribution as the lower bound of black hole masses. The best model (the power law) fitted to the low-mass systems has a distribution of lower-bounds with M_1% > 4.3 Msun with 90% confidence, while the best model (the exponential) fitted to all 20 systems has M_1% > 4.5 Msun with 90% confidence. We conclude that our sample of black hole masses provides strong evidence of a gap between the maximum neutron star mass and the lower bound on black hole masses. Our results on the low-mass sample are in qualitative agreement with those of Ozel, et al (2010)., 56 pages, 22 figures, 9 tables, as accepted by ApJ

Although not nearly as numerous as binaries with two white dwarfs, eccentric neutron star-white dwarf (NS-WD) binaries are important gravitational-wave (GW) sources for the next generation of space-based detectors sensitive to low frequency waves. Here we investigate periastron precession in these sources as a result of general relativistic, tidal, and rotational effects; such precession is expected to be detectable for at least some of the detected binaries of this type. Currently, two eccentric NS-WD binaries are known in the galactic field, PSR J1141-6545 and PSR B2303+46, both of which have orbits too wide to be relevant in their current state to GW observations. However, population synthesis studies predict the existence of a significant Galactic population of such systems. Though small in most of these systems, we find that tidally induced periastron precession becomes important when tides contribute to more than 3% of the total precession rate. For these systems, accounting for tides when analyzing periastron precession rate measurements can improve estimates of the inferred WD component mass and, in some cases, will prevent us from misclassifying the object. However, such systems are rare, due to rapid orbital decay. To aid the inclusion of tidal effects when using periastron precession as a mass measurement tool, we derive a function that relates the WD radius and periastron precession constant to the WD mass., Facultad de Ciencias Astronómicas y Geofísicas