Kishalay De, Morgan MacLeod, Jacob E. Jencson, Elizabeth Lovegrove, Andrea Antoni, Erin Kara, Mansi M. Kasliwal, Ryan M. Lau, Abraham Loeb, Megan Masterson, Aaron M. Meisner, Christos Panagiotou, Eliot Quataert, Robert Simcoe Stellar mass black holes are formed from the terminal collapse of massive stars if the ensuing neutrino shock is unable to eject the stellar envelope. Direct observations of black hole formation remain inconclusive. We report observations of M31-2014-DS1, a massive, hydrogen-depleted supergiant in the Andromeda galaxy identified via a mid-infrared brightening in 2014. Its total luminosity remained nearly constant for the subsequent thousand days, before fading dramatically over the next thousand days by $\gtrsim 10\times$ and $\gtrsim 10^4\times$ in total and visible light, respectively. Together with the lack of a detected optical outburst, the observations are explained by the fallback of the stellar envelope into a newly formed black hole, moderated by the injection of a $\sim 10^{48}$ erg shock. Unifying these observations with a candidate in NGC 6946, we present a concordant picture for the birth of stellar mass black holes from stripped massive stars.
Melissa Shahbandeh, Ori D. Fox, Tea Temim, Eli Dwek, Arkaprabha Sarangi, Nathan Smith, Luc Dessart, Bryony Nickson, Michael Engesser, Alexei V. Filippenko, Thomas G. Brink, Weikang Zheng, Tamás Szalai, Joel Johansson, Armin Rest, Schuyler D. Van Dyk, Jennifer Andrews, Chris Ashall, Geoffrey C. Clayton, Ilse De Looze, et al (25) Dust from core-collapse supernovae (CCSNe), specifically Type IIP SNe, has been suggested to be a significant source of the dust observed in high-redshift galaxies. CCSNe eject large amounts of newly formed heavy elements, which can condense into dust grains in the cooling ejecta. However, infrared (IR) observations of typical CCSNe generally measure dust masses that are too small to account for the dust production needed at high redshifts. Type IIn SNe, classified by their dense circumstellar medium (CSM), are also known to exhibit strong IR emission from warm dust, but the dust origin and heating mechanism have generally remained unconstrained because of limited observational capabilities in the mid-IR. Here, we present a JWST/MIRI Medium Resolution Spectrograph (MRS) spectrum of the Type IIn SN 2005ip nearly 17 years post-explosion. The Type IIn SN 2005ip is one of the longest-lasting and most well-studied SNe observed to date. Combined with a Spitzer mid-IR spectrum of SN 2005ip obtained in 2008, this data set provides a rare 15-year baseline, allowing for a unique investigation of the evolution of dust. The JWST spectrum shows a new high-mass dust component ($\gtrsim0.08$ M$_{\odot}$) that is not present in the earlier Spitzer spectrum. Our analysis shows dust likely formed over the past 15 years in the cold, dense shell (CDS), between the forward and reverse shocks. There is also a smaller mass of carbonaceous dust ($\gtrsim0.005$ M$_{\odot}$) in the ejecta. These observations provide new insights into the role of SN dust production, particularly within the CDS, and its potential contribution to the rapid dust enrichment of the early Universe.
Noel D. Richardson, Gail H. Schaefer, Jan J. Eldridge, Rebecca Spejcher, Amanda Holdsworth, Ryan M. Lau, John D. Monnier, Anthony F. J. Moffat, Gerd Weigelt, Peredur M. Williams, Stefan Kraus, Jean-Baptiste Le Bouquin, Narsireddy Anugu, Sorabh Chhabra, Isabelle Codron, Jacob Ennis, Tyler Gardner, Mayra Gutierrez, Noura Ibrahim, Aaron Labdon, et al (2) Classical Wolf-Rayet stars are the descendants of massive OB stars that have lost their hydrogen envelopes and are burning helium in their cores prior to exploding as type Ib/c supernovae. The mechanisms for losing their hydrogen envelopes are either through binary interactions or through strong stellar winds potentially coupled with episodic mass-loss. Amongst the bright classical WR stars, the binary system WR\u2009137 (HD\u2009192641; WC7d + O9e) is the subject of this paper. This binary is known to have a 13-year period and produces dust near periastron. Here we report on interferometry with the CHARA Array collected over a decade of time and providing the first visual orbit for the system. We combine these astrometric measurements with archival radial velocities to measure masses of the stars of $M_{\rm WR} = 9.5\pm3.4 M_\odot$ and $M_{\rm O} = 17.3\pm 1.9 M_\odot$ when we use the most recent \textitGaia distance. These results are then compared to predicted dust distribution using these orbital elements, which match the observed imaging from \textitJWST as discussed recently by Lau et al. Furthermore, we compare the system to the BPASS models, finding that the WR star likely formed through stellar winds and not through binary interactions. However, the companion O star did likely accrete some material from the WR's mass-loss to provide the rotation seen today that drives its status as an Oe star.
Masafumi Niwano, Michael M. Fausnaugh, Ryan M. Lau, Kishalay De, Roberto Soria, George R. Ricker, Roland Vanderspek, Michael C. B. Ashley, Nicholas Earley, Matthew J. Hankins, Mansi M. Kasliwal, Anna M. Moore, Jamie Soon, Tony Travouillon, Mahito Sasada, Ichiro Takahashi, Yoichi Yatsu, Nobuyuki Kawai The mechanism of X-ray outbursts in Be X-ray binaries remains a mystery, and understanding their circumstellar disks is crucial for a solution of the mass-transfer problem. In particular, it is important to identify the Be star activities (e.g., pulsations) that cause mass ejection and, hence, disk formation. Therefore, we investigated the relationship between optical flux oscillations and the infrared (IR) excess in a sample of five Be X-ray binaries. Applying the Lomb-Scargle technique to high-cadence optical light curves from the Transiting Exoplanet Survey Satellite (TESS), we detected several significant oscillation modes in the 3 to 24 hour period range for each source. We also measured the IR excess (a proxy for disk growth) of those five sources, using J-band light curves from Palomar Gattini-IR. In four of the five sources, we found anti-correlations between the IR excess and the amplitude of the main flux oscillation modes. This result is inconsistent with the conventional idea that non-radial pulsations drive mass ejections. We propose an alternative scenario where internal temperature variations in the Be star cause transitions between pulsation-active and mass-ejection-active states.
We present the substellar mass function of star-forming clusters ($\simeq$0.1 Myr old) in a low-metallicity environment ($\simeq$$-$0.7 dex). We performed deep JWST/NIRCam and MIRI imaging of two star-forming clusters in Digel Cloud 2, a star-forming region in the Outer Galaxy ($R_G \gtrsim 15$ kpc). The very high sensitivity and spatial resolution of JWST enable us to resolve cluster members clearly down to a mass detection limit of 0.02 $M_\odot$, enabling the first detection of brown dwarfs in low-metallicity clusters. Fifty-two and ninety-one sources were extracted in mass-$A_V$-limited samples in the two clusters, from which Initial mass functions (IMFs) were derived by model-fitting the F200W band luminosity function, resulting in IMF peak masses (hereafter $M_C$) $\log M_C / M_\odot \simeq -1.5 \pm 0.5$ for both clusters. Although the uncertainties are rather large, the obtained $M_C$ values are lower than those in any previous study ($\log M_C / M_\odot \sim -0.5$). Comparison with the local open clusters with similar ages to the target clusters ($\sim$$10^6$-$10^7$ yr) suggests a metallicity dependence of $M_C$, with lower $M_C$ at lower metallicities, while the comparison with globular clusters, similarly low metallicities but considerably older ($\sim$$10^{10}$ yr), suggests that the target clusters have not yet experienced significant dynamical evolution and remain in their initial physical condition. The lower $M_C$ is also consistent with the theoretical expectation of the lower Jeans mass due to the higher gas density under such low metallicity. The $M_C$ values derived from observations in such an environment would place significant constraints on the understanding of star formation.
The James Webb Space Telescope (JWST) has opened up a new window to study highly reddened explosive transients. We present results from late-time (1421 days post-explosion) JWST follow-up spectroscopic observations with NIRSpec and MIRI LRS of the intermediate luminosity red transient (ILRT) AT 2019abn located in the nearby Messier 51 galaxy (8.6 Mpc). ILRTs represent a mysterious class of transients which exhibit peak luminosities between those of classical novae and supernovae and which are known to be highly dust obscured. Similar to the prototypical examples of this class of objects, NGC 300 2008-OT and SN 2008S, AT 2019abn has an extremely red and dusty progenitor detected only in pre-explosion Spitzer/IRAC imaging at 3.6 and 4.5 micron and not in deep optical or near-infrared HST images. We find that late time observations of AT 2019abn from NEOWISE and JWST are consistent with the late time evolution of SN 2008S. In part because they are so obscured by dust, it is unknown what produces an ILRT with hypotheses ranging from high mass stellar merger events, non-terminal stellar outbursts, or terminal supernovae explosions through electron-capture in super-AGB stars. Our JWST observations show strong mid-IR Class C PAH features at 6.3 and 8.25 micron typical of carbon-rich post-AGB sources. These features suggest the dust around AT 2019abn, either pre-existing or newly formed in the ejecta, is composed of carbonaceous grains which are not typically observed around red supergiants. However, depending on the strength and temperature of hot bottom burning, SAGBs may be expected to exhibit a carbon-rich chemistry. Thus our JWST observations are consistent with AT 2019abn having an SAGB progenitor.
Angela S. Cotera, Matthew J. Hankins, John Bally, Ashley T. Barnes, Cara D. Battersby, H Perry Hatchfield, Terry L. Herter, Ryan M. Lau, Steven N. Longmore, Elisabeth A. C. Mills, Mark R. Morris, James T. Radomski, Janet P. Simpson, Zachary Stephens, Daniel L. Walker The central regions of the Milky Way constitute a unique laboratory for a wide swath of astrophysical studies, consequently the inner $\sim$400 pc has been the target of numerous large surveys at all accessible wavelengths. In this paper we present a catalog of sources at 25 and 37 $\mu$m located within all of the regions observed with the SOFIA/FORCAST instrument in the inner $\sim$200 pc of the Galaxy. The majority of the observations were obtained as part of the SOFIA Cycle 7 Galactic Center Legacy program survey, which was designed to complement the Spitzer/MIPS 24 $\mu$m catalog in regions saturated in the MIPS observations. Due to the wide variety of source types captured by our observations at 25 and 37 $\mu$m, we do not limit the FORCAST source catalog to unresolved point sources, or treat all sources as if they are point-like sources. The catalog includes all detectable sources in the regions, resulting in a catalog of 950 sources, including point sources, compact sources, and extended sources. We also provide the user with metrics to discriminate between the source types.
Viraj R. Karambelkar, Mansi M. Kasliwal, Patrick Tisserand, Shreya Anand, Michael C. B. Ashley, Lars Bildsten, Geoffrey C. Clayton, Courtney C. Crawford, Kishalay De, Nicholas Earley, Matthew J. Hankins, Xander Hall, Astrid Lamberts, Ryan M. Lau, Dan McKenna, Anna Moore, Eran O. Ofek, Roger M. Smith, Roberto Soria, Jamie Soon, et al (1) We present results from a systematic infrared (IR) census of R Coronae Borealis (RCB) stars in the Milky Way, using data from the Palomar Gattini IR (PGIR) survey. R Coronae Borealis stars are dusty, erratic variable stars presumably formed from the merger of a He-core and a CO-core white dwarf (WD). PGIR is a 30 cm $J$-band telescope with a 25 deg$^{2}$ camera that surveys 18000 deg$^{2}$ of the northern sky ($\delta>-28^{o}$) at a cadence of 2 days. Using PGIR J-band lightcurves for $\sim$60 million stars together with mid-IR colors from WISE, we selected a sample of 530 candidate RCB stars. We obtained near-IR spectra for these candidates and identified 53 RCB stars in our sample. Accounting for our selection criteria, we find that there are a total of $\approx350^{+150}_{-100}$ RCB stars in the Milky Way. Assuming typical RCB lifetimes, this corresponds to an RCB formation rate of 0.8 - 5 $\times$ 10$^{-3}$ yr$^{-1}$, consistent with observational and theoretical estimates of the He-CO WD merger rate. We searched for quasi-periodic pulsations in the PGIR lightcurves of RCB stars and present pulsation periods for 16 RCB stars. We also examined high-cadenced TESS lightcurves for RCB and the chemically similar, but dustless hydrogen-deficient carbon (dLHdC) stars. We find that dLHdC stars show variations on timescales shorter than RCB stars, suggesting that they may have lower masses than RCB stars. Finally, we identified 3 new spectroscopically confirmed and 12 candidate Galactic DY Per type stars - believed to be colder cousins of RCB stars - doubling the sample of Galactic DY Per type stars.
The extreme outer Galaxy (EOG), which we define as the region of the Milky Way with a galactocentric radius of more than 18 kpc, provides an excellent opportunity to study star formation in an environment significantly different from that in the solar neighborhood because of its lower metallicity and lower gas density. We carried out near- and mid-infrared (NIR and MIR) imaging observations toward two star-forming clusters located in the EOG using JWST NIRCam and MIRI with nine filters: F115W, F150W, F200W, F350W, F405N, F444W, F770W, F1280W, and F2100W. In this paper, we present an overview of the observations, data reduction, and initial results. The NIR sensitivity is approximately 10--80 times better than our previous observation with the Subaru 8.2 m telescope. Accordingly, the mass detection limit reaches to about 0.01--0.05 $M_\odot$, which is about 10 times better than the previous observations. At MIR wavelengths, the high sensitivity and resolution data enable us to resolve individual young stellar objects in such a distant region for the first time. The mass detection limit at MIR F770W filter reaches about 0.1--0.3 $M_\odot$. With these new observations, we have identified components of the clusters that previous surveys did not detect, including class 0 candidates, outflow/jet components, and distinctive nebular structures. These data will enable us to investigate the properties of star formation in the EOG at the same depth of detail as previous observations of star formation in the solar neighborhood.
Jayadev K. Rajagopal, Ryan M. Lau, Isack Padilla, Stephen T. Ridgway, Chaohan Cui, Brittany McClinton, Aqil Sajjad, Stuartt Corder, Mark Rawlings, Fredrik Rantakyro, J. Gabriel Richardson, Amit Ashok, Saikat Guha Microarcsecond resolutions afforded by an optical-NIR array with kilometer-baselines would enable breakthrough science. However significant technology barriers exist in transporting weakly coherent photon states over these distances: primarily photon loss and phase errors. Quantum telescopy, using entangled states to link spatially separated apertures, offers a possible solution to the loss of photons. We report on an initiative launched by NSF NOIRLab in collaboration with the Center for Quantum Networks and Arizona Quantum Initiative at the University of Arizona, Tucson, to explore these concepts further. A brief description of the quantum concepts and a possible technology roadmap towards a quantum-enhanced very long baseline optical-NIR interferometric array is presented. An on-sky demonstration of measuring spatial coherence of photons with apertures linked through the simplest Gottesman protocol over short baselines and with limited phase fluctuations is envisaged as the first step.
Palomar Gattini-IR is a wide-field, synoptic infrared time domain survey covering $\approx 15000$ sq. deg. of the sky at $\approx 1-3$ night cadence to a depth of $J\approx 13.0$ and $\approx 14.9$ Vega mag in and outside the Galactic plane, respectively. Here, we present the first data release of $J$-band light curves of 2MASS sources within the survey footprint covering approximately the first four years of operations. We describe the construction of the source catalog based on the 2MASS point sources, followed by exposure filtering criteria and forced PSF photometry. The catalog contains light curves of $\approx 286$ million unique sources with 2MASS magnitudes of $J < 15.5$ mag, with a total of $\approx 50$ billion photometric measurements and $\approx 20$ billion individual source detections at signal-to-noise-ratio $> 3$. We demonstrate the photometric fidelity of the catalog by i) quantifying the magnitude-dependent accuracy and uncertainty of the photometry with respect to 2MASS and ii) comparing against forced PGIR aperture photometry for known variable sources. We present simple filtering criteria for selecting reliable photometric measurements as well as example Python notebooks for users. This catalog is the largest compilation of nightly cadence, synoptic infrared light curves to date, comparable to those in the largest optical surveys, providing a stepping stone to upcoming infrared surveys in the coming decade.
Charlotte Myers, Kishalay De, Lin Yan, Jacob E. Jencson, Nicholas Earley, Christoffer Fremling, Daichi Hiramatsu, Mansi M. Kasliwal, Ryan M. Lau, Morgan MacLeod, Megan Masterson, Christos Panagiotou, Robert Simcoe, Samaporn Tinyanont While core-collapse supernovae (SNe) often show early and consistent signs of circumstellar (CSM) interaction, some exhibit delayed signatures due to interaction with distant material around the progenitor star. Here we present the discovery in NEOWISE data of WTP19aalnxx, a luminous mid-infrared (IR) transient in the outskirts of the galaxy KUG 0022-007 at $\approx 190$ Mpc. First detected in 2018, WTP19aalnxx reaches a peak absolute (Vega) magnitude of $\approx-22$ at $4.6 \, \mu$m in $\approx3$ yr, comparable to the most luminous interacting SNe. Archival data reveal a $\gtrsim 5\times$ fainter optical counterpart detected since 2015, while follow-up near-IR observations in 2022 reveal an extremely red ($Ks-W2 \approx 3.7$ mag) active transient. Deep optical spectroscopy confirm strong CSM interaction signatures via intermediate-width Balmer emission lines and coronal metal lines. Modeling the broadband spectral energy distribution, we estimate the presence of $\gtrsim 10^{-2}$ M$_\odot$ of warm dust, likely formed in the shock interaction region. Together with the lack of nebular Fe emission, we suggest that WTP19aalnxx is a missed, low (optical) luminosity SN in an emerging family of core-collapse SNe distinguished by their CSM-interaction-powered mid-IR emission that outshines the optical bands. Investigating the Zwicky Transient Facility sample of SNe in NEOWISE data, we find $17$ core-collapse SNe ($\gtrsim 3$% in a volume-limited sample) without early signs of CSM interaction that exhibit delayed IR brightening, suggestive of dense CSM shells at $\lesssim 10^{17}$cm. We suggest that synoptic IR surveys offer a new route to revealing late-time CSM interaction and the prevalence of intense terminal mass loss in massive stars.
Noel D. Richardson, Andrea R. Daly, Peredur M. Williams, Grant M. Hill, Victor I. Shenavrin, Izumi Endo, André-Nicolas Chené, Nicole Karnath, Ryan M. Lau, Anthony F. J. Moffat, Gerd Weigelt Several long-period binaries with a carbon-rich Wolf-Rayet star and an O star produce dust in their wind collisions. In eccentric binaries, this is seen most strongly near periastron passage. The exact conditions leading to dust creation require orbital properties to be determined, which is difficult owing to their long periods. Recently, the binary system WR 125 (WC7+O9III) began a dust creation episode seen through an infrared outburst first detected by NEOWISE-R, which was the first outburst detected since 1991. We present new near- and mid-infrared photometry, which we use to show consistency between the two outbursts and derive an orbital period of 28.12$^{+0.10}_{-0.05}$ yr. We use a long time-series of optical spectra to place the first constraints on its orbital elements, on the assumption that this system will produce dust near periastron. The orbit has a mild eccentricity of 0.29$\pm$0.12 and is only derived for the Wolf-Rayet component, as the O star's radial velocities have noise that is likely larger than the expected semi-amplitude of the orbit. We also present SOFIA/FORCAST grism spectroscopy to examine the infrared spectral energy distribution (SED) of the dust during this outburst, comparing its properties to other WCd binaries, deriving a dust temperature of 580 K in 2021. This collection of observations will allow us to plan future observations of this system and place the system in the context of dust-creating Wolf-Rayet binaries.
C. Fransson, M. J. Barlow, P. J. Kavanagh, J. Larsson, O. C. Jones, B. Sargent, M. Meixner, P. Bouchet, T. Temim, G. S. Wright, J. A. D. L. Blommaert, N. Habel, A. S. Hirschauer, J. Hjorth, L. Lenkić, T. Tikkanen, R. Wesson, A. Coulais, O. D. Fox, R. Gastaud, et al (14) The nearby Supernova 1987A was accompanied by a burst of neutrino emission, which indicates that a compact object (a neutron star or black hole) was formed in the explosion. There has been no direct observation of this compact object. In this work, we observe the supernova remnant with JWST spectroscopy finding narrow infrared emission lines of argon and sulphur. The line emission is spatially unresolved and blueshifted in velocity relative to the supernova rest frame. We interpret the lines as gas illuminated by a source of ionizing photons located close to the center of the expanding ejecta. Photoionization models show that the line ratios are consistent with ionization by a cooling neutron star or pulsar wind nebula. The velocity shift could be evidence for a neutron star natal kick.
P. Bouchet, R. Gastaud, A. Coulais, M. J. Barlow, C. Fransson, P. J. Kavanagh, J. Larsson, T. Temim, O. C. Jones, A. S. Hirschauer, T. Tikkanen, J. A. D. L. Blommaert, O. D. Fox, A. Glasse, N. Habel, J. Hjorth, J. Jaspers, O. Krause, R. M. Lau, L. Lenkić, et al (14) There exist very few mid-infrared (IR) observations of supernovae (SNe) in general. Therefore, SN 1987A, the closest visible SN in 400 years, gives us the opportunity to explore the mid-IR properties of SNe, the dust in their ejecta and surrounding medium, and to witness the birth of a SN remnant (SNR). The James Webb Space Telescope (JWST), with its high spatial resolution and extreme sensitivity, gives a new view on these issues. We report on the first imaging observations obtained with the Mid-InfraRed Instrument (MIRI). We build temperature maps and discuss the morphology of the nascent SNR. Our results show that the temperatures in the equatorial ring (ER) are quite non-uniform. This could be due to dust destruction in some parts of the ring, as had been assumed in some previous works. We show that the IR emission extends beyond the ER, illustrating the fact that the shock wave has now passed through this ring to affect the circumstellar medium on a larger scale. Finally, while sub-mm Atacama Large Millimeter Array (ALMA) observations have hinted at the location of the compact remnant of SN 1987A, we note that our MIRI data have found no such evidence.
2023 Windows on, Universe Workshop White Paper Working Group, T. Ahumada, J. E. Andrews, S. Antier, E. Blaufuss, P. R. Brady, A. M. Brazier, E. Burns, S. B. Cenko, P. Chandra, D. Chatterjee, A. Corsi, M. W. Coughlin, D. A. Coulter, S. Fu, A. Goldstein, L. P. Guy, E. J. Hooper, S. B. Howell, et al (27) In this White Paper, we present recommendations for the scientific community and funding agencies to foster the infrastructure for a collaborative multi-messenger and time-domain astronomy (MMA/TDA) ecosystem. MMA/TDA is poised for breakthrough discoveries in the coming decade. In much the same way that expanding beyond the optical bandpass revealed entirely new and unexpected discoveries, cosmic messengers beyond light (i.e., gravitational waves, neutrinos, and cosmic rays) open entirely new windows to answer some of the most fundamental questions in (astro)physics: heavy element synthesis, equation of state of dense matter, particle acceleration, etc. This field was prioritized as a frontier scientific pursuit in the 2020 Decadal Survey on Astronomy and Astrophysics via its "New Windows on the Dynamic Universe" theme. MMA/TDA science presents technical challenges distinct from those experienced in other disciplines. Successful observations require coordination across myriad boundaries -- different cosmic messengers, ground vs. space, international borders, etc. -- all for sources that may not be well localized, and whose brightness may be changing rapidly with time. Add that all of this work is undertaken by real human beings, with distinct backgrounds, experiences, cultures, and expectations, that often conflict. To address these challenges and help MMA/TDA realize its full scientific potential in the coming decade (and beyond), the second in a series of community workshops sponsored by the U.S. National Science Foundation (NSF) and NASA titled "Windows on the Universe: Establishing the Infrastructure for a Collaborative Multi-Messenger Ecosystem" was held on October 16-18, 2023 in Tucson, AZ. Here we present the primary recommendations from this workshop focused on three key topics -- hardware, software, and people and policy. [abridged]
Ryan M. Lau, Matthew J. Hankins, Joel Sanchez-Bermudez, Deepashri Thatte, Anthony Soulain, Rachel A. Cooper, Anand Sivaramakrishnan, Michael F. Corcoran, Alexandra Z. Greenbaum, Theodore R. Gull, Yinuo Han, Olivia C. Jones, Thomas Madura, Anthony F. J. Moffat, Mark R. Morris, Takashi Onaka, Christopher M. P. Russell, Noel D. Richardson, Nathan Smith, Peter Tuthill, et al (3) We present infrared aperture masking interferometry (AMI) observations of newly formed dust from the colliding winds of the massive binary system Wolf-Rayet (WR) 137 with JWST using the Near Infrared Imager and Slitless Spectrograph (NIRISS). NIRISS AMI observations of WR 137 and a point-spread-function calibrator star, HD~228337, were taken using the F380M and F480M filters in 2022 July and August as part of the Director's Discretionary Early Release Science (DD-ERS) program 1349. Interferometric observables (squared visibilities and closure phases) from the WR 137 "interferogram" were extracted and calibrated using three independent software tools: ImPlaneIA, AMICAL, and SAMpip. The analysis of the calibrated observables yielded consistent values except for slightly discrepant closure phases measured by ImPlaneIA. Based on all three sets of calibrated observables, images were reconstructed using three independent software tools: BSMEM, IRBis, and SQUEEZE. All reconstructed image combinations generated consistent images in both F380M and F480M filters. The reconstructed images of WR 137 reveal a bright central core with a $\sim300$ mas linear filament extending to the northwest. A geometric colliding-wind model with dust production constrained to the orbital plane of the binary system and enhanced as the system approaches periapsis provided a general agreement with the interferometric observables and reconstructed images. Based on a colliding-wind dust condensation analysis, we suggest that dust formation within the orbital plane of WR 137 is induced by enhanced equatorial mass-loss from the rapidly rotating O9 companion star, whose axis of rotation is aligned with that of the orbit.
Rainer Schoedel, Steve Longmore, Jonny Henshaw, Adam Ginsburg, John Bally, Anja Feldmeier, Matt Hosek, Francisco Nogueras Lara, Anna Ciurlo, Mélanie Chevance, J. M. Diederik Kruijssen, Ralf Klessen, Gabriele Ponti, Pau Amaro-Seoane, Konstantina Anastasopoulou, Jay Anderson, Maria Arias, Ashley T. Barnes, Cara Battersby, Giuseppe Bono, et al (90) The inner hundred parsecs of the Milky Way hosts the nearest supermassive black hole, largest reservoir of dense gas, greatest stellar density, hundreds of massive main and post main sequence stars, and the highest volume density of supernovae in the Galaxy. As the nearest environment in which it is possible to simultaneously observe many of the extreme processes shaping the Universe, it is one of the most well-studied regions in astrophysics. Due to its proximity, we can study the center of our Galaxy on scales down to a few hundred AU, a hundred times better than in similar Local Group galaxies and thousands of times better than in the nearest active galaxies. The Galactic Center (GC) is therefore of outstanding astrophysical interest. However, in spite of intense observational work over the past decades, there are still fundamental things unknown about the GC. JWST has the unique capability to provide us with the necessary, game-changing data. In this White Paper, we advocate for a JWST NIRCam survey that aims at solving central questions, that we have identified as a community: i) the 3D structure and kinematics of gas and stars; ii) ancient star formation and its relation with the overall history of the Milky Way, as well as recent star formation and its implications for the overall energetics of our galaxy's nucleus; and iii) the (non-)universality of star formation and the stellar initial mass function. We advocate for a large-area, multi-epoch, multi-wavelength NIRCam survey of the inner 100\u2009pc of the Galaxy in the form of a Treasury GO JWST Large Program that is open to the community. We describe how this survey will derive the physical and kinematic properties of ~10,000,000 stars, how this will solve the key unknowns and provide a valuable resource for the community with long-lasting legacy value.
Szanna Zsíros, Tamás Szalai, Ilse De Looze, Arkaprabha Sarangi, Melissa Shahbandeh, Ori D. Fox, Tea Temim, Dan Milisavljevic, Schuyler D. Van Dyk, Nathan Smith, Alexei V. Filippenko, Thomas G. Brink, WeiKang Zheng, Luc Dessart, Jacob Jencson, Joel Johansson, Justin Pierel, Armin Rest, Samaporn Tinyanont, Maria Niculescu-Duvaz, et al (22) We present mid-infrared (mid-IR) imaging of the Type IIL supernova (SN) 1980K with the James Webb Space Telescope (JWST) more than 40 yr post-explosion. SN 1980K, located in the nearby ($D\approx7$ Mpc) "SN factory" galaxy NGC 6946, was serendipitously captured in JWST/MIRI images taken of the field of SN 2004et in the same galaxy. SN 1980K serves as a promising candidate for studying the transitional phase between young SNe and older SN remnants and also provides a great opportunity to investigate its the close environment. SN 1980K can be identified as a clear and bright point source in all eight MIRI filters from F560W up to F2550W. We fit analytical dust models to the mid-IR spectral energy distribution that reveal a large amount ($M_d \approx 0.002 {M}_{\odot}$) of Si-dominated dust at $T_{dust}\approx 150$ K (accompanied by a hotter dust/gas component), and also computed numerical SED dust models. Radiative transfer modeling of a late-time optical spectrum obtained recently with Keck discloses that an even larger ($\sim 0.24-0.58~{M}_{\odot}$) amount of dust is needed in order for selective extinction to explain the asymmetric line profile shapes observed in SN 1980K. As a conclusion, with JWST, we may see i) pre-existing circumstellar dust heated collisionally (or, partly radiatively), analogous to the equatorial ring of SN 1987A, or ii) the mid-IR component of the presumed newly-formed dust, accompanied by much more colder dust present in the ejecta (as suggested by the late-time the optical spectra).
WR 137 (HD 192641) is a binary system consisting of a carbon-rich Wolf-Rayet star and an Oe companion star in a 13-year orbit. Near periastron, the winds of the two stars collide and form carbonaceous dust. We obtained three mid-infrared grism spectra of the system with SOFIA and FORCAST during the last year of SOFIA's operations in July 2021, February 2021, and May 2022 (Cycle 9). Within these spectra, we have identified several wind lines from He I, He II, C III, and C IV that are emitted from the Wolf-Rayet wind as well as a weak emission feature around 6.3-6.4 $\mu$m that may have shifted its peak flux from 6.29 to 6.41$\mu$m through this time period. The weak feature grew as the continuum dust emission grew while the WR emission appeared to decline due to lower contrast with the continuum. Furthermore, we observe that the peak of the feature shifts to redder wavelengths during the observations. We compare this feature to the UIR feature and other emission lines identified in dusty WC binaries. For WR 137, we speculate that mixing of the winds in the system with the Oe star's disk is important for starting the dust formation and that it is less important as dust formation continues. Previous infrared photometry shows "mini-eruptions" of dust production which could then be explained with variations of the Oe star disk.
O. C. Jones, P. J. Kavanagh, M. J. Barlow, T. Temim, C. Fransson, J. Larsson, J. A. D. L. Blommaert, M. Meixner, R. M. Lau, B. Sargent, P. Bouchet, J. Hjorth, G. S. Wright, A. Coulais, O. D. Fox, R. Gastaud, A. Glasse, N. Habel, A. S. Hirschauer, J. Jaspers, et al (14) Supernova (SN) 1987A is the nearest supernova in $\sim$400 years. Using the \em JWST MIRI Medium Resolution Spectrograph, we spatially resolved the ejecta, equatorial ring (ER) and outer rings in the mid-infrared 12,927 days after the explosion. The spectra are rich in line and dust continuum emission, both in the ejecta and the ring. Broad emission lines (280-380~km~s$^{-1}$ FWHM) seen from all singly-ionized species originate from the expanding ER, with properties consistent with dense post-shock cooling gas. Narrower emission lines (100-170~km~s$^{-1}$ FWHM) are seen from species originating from a more extended lower-density component whose high ionization may have been produced by shocks progressing through the ER, or by the UV radiation pulse associated with the original supernova event. The asymmetric east-west dust emission in the ER has continued to fade, with constant temperature, signifying a reduction in dust mass. Small grains in the ER are preferentially destroyed, with larger grains from the progenitor surviving the transition from SN into SNR. The ER is fit with a single set of optical constants, eliminating the need for a secondary featureless hot dust component. We find several broad ejecta emission lines from [Ne~\sc ii], [Ar~\sc ii], [Fe~\sc ii], and [Ni~\sc ii]. With the exception of [Fe~\sc ii]~25.99$\mu$m, these all originate from the ejecta close to the ring and are likely being excited by X-rays from the interaction. The [Fe~\sc ii]~5.34$\mu$m to 25.99$\mu$m line ratio indicates a temperature of only a few hundred K in the inner core, consistent with being powered by ${}^{44}$Ti decay.
R. Jain, E. F. Brown, H. Schatz, A. V. Afanasjev, M. Beard, L. R. Gasques, S. S. Gupta, G. W. Hitt, W. R. Hix, R. Lau, P. Moller, W. J. Ong, M. Wiescher, Y. Xu The observation of X-rays during quiescence from transiently accreting neutron stars provides unique clues about the nature of dense matter. This, however, requires extensive modeling of the crusts and matching the results to observations. The pycnonuclear fusion reaction rates implemented in these models are theoretically calculated by extending phenomenological expressions and have large uncertainties spanning many orders of magnitude. We present the first sensitivity studies of these pycnonuclear fusion reactions in realistic network calculations. We also couple the reaction network with the thermal evolution code dStar to further study their impact on the neutron star cooling curves in quiescence. Varying the pycnonuclear fusion reaction rates alters the depth at which nuclear heat is deposited although the total heating remains constant. The enhancement of the pycnonuclear fusion reaction rates leads to an overall shallower deposition of nuclear heat. The impurity factors are also altered depending on the type of ashes deposited on the crust. These total changes correspond to a variation of up to 9 eV in the modeled cooling curves. While this is not sufficient to explain the shallow heat source, it is comparable to the observational uncertainties and can still be important for modeling the neutron star crust.
Jacob E. Jencson, Jeniveve Pearson, Emma R. Beasor, Ryan M. Lau, Jennifer E. Andrews, K. Azalee Bostroem, Yize Dong, Michael Engesser, Sebastian Gomez, Muryel Guolo, Emily Hoang, Griffin Hosseinzadeh, Saurabh W. Jha, Viraj Karambelkar, Mansi M. Kasliwal, Michael Lundquist, Nicolas E. Meza Retamal, Armin Rest, David J. Sand, Melissa Shahbandeh, et al (6) We analyze pre-explosion near- and mid-infrared (IR) imaging of the site of SN 2023ixf in the nearby spiral galaxy M101 and characterize the candidate progenitor star. The star displays compelling evidence of variability with a possible period of $\approx$1000 days and an amplitude of $\Delta m \approx 0.6$ mag in extensive monitoring with the Spitzer Space Telescope since 2004, likely indicative of radial pulsations. Variability consistent with this period is also seen in the near-IR $J$ and $K_{s}$ bands between 2010 and 2023, up to just 10 days before the explosion. Beyond the periodic variability, we do not find evidence for any IR-bright pre-supernova outbursts in this time period. The IR brightness ($M_{K_s} = -10.7$ mag) and color ($J-K_{s} = 1.6$ mag) of the star suggest a luminous and dusty red supergiant. Modeling of the phase-averaged spectral energy distribution (SED) yields constraints on the stellar temperature ($T_{\mathrm{eff}} = 3500_{-1400}^{+800}$ K) and luminosity ($\log L/L_{\odot} = 5.1\pm0.2$). This places the candidate among the most luminous Type II supernova progenitors with direct imaging constraints, with the caveat that many of these rely only on optical measurements. Comparison with stellar evolution models gives an initial mass of $M_{\mathrm{init}} = 17\pm4 M_{\odot}$. We estimate the pre-supernova mass-loss rate of the star between 3 and 19 yr before explosion from the SED modeling at $\dot M \approx 3\times10^{-5}$ to $3\times10^{-4} M_{\odot}$ yr$^{-1}$ for an assumed wind velocity of $v_w = 10$ km s$^{-1}$, perhaps pointing to enhanced mass loss in a pulsation-driven wind.
Ryan M. Lau, Jason Wang, Matthew J. Hankins, Thayne Currie, Vincent Deo, Izumi Endo, Olivier Guyon, Yinuo Han, Anthony P. Jones, Nemanja Jovanovic, Julien Lozi, Anthony F. J. Moffat, Takashi Onaka, Garreth Ruane, Andreas A. C. Sander, Samaporn Tinyanont, Peter G. Tuthill, Gerd Weigelt, Peredur M. Williams, Sebastien Vievard Wolf-Rayet (WR) 140 is the archetypal periodic dust-forming colliding-wind binary that hosts a carbon-rich WR (WC) star and an O-star companion with an orbital period of 7.93 years and an orbital eccentricity of 0.9. Throughout the past several decades, multiple dust-formation episodes from WR 140 have been observed that are linked to the binary orbit and occur near the time of periastron passage. Given its predictable dust-formation episodes, WR 140 presents an ideal astrophysical laboratory for investigating the formation and evolution of dust in the hostile environment around a massive binary system. In this paper, we present near- and mid-infrared (IR) spectroscopic and imaging observations of WR 140 with Subaru/SCExAO+CHARIS, Keck/NIRC2+PyWFS, and Subaru/COMICS taken between 2020 June and Sept that resolve the circumstellar dust emission linked to its most recent dust-formation episode in 2016 Dec. Our spectral energy distribution (SED) analysis of WR 140's resolved circumstellar dust emission reveals the presence of a hot ($T_\mathrm{d}\sim1000$ K) near-IR dust component that is co-spatial with the previously known and cooler ($T_\mathrm{d}\sim500$ K) mid-IR dust component composed of $300-500$ Å-sized dust grains. We attribute the hot near-IR dust emission to the presence of nano-sized ("nanodust") grains and suggest they were formed from grain-grain collisions or the rotational disruption of the larger grain size population by radiative torques in the strong radiation field from the central binary. Lastly, we speculate on the astrophysical implications of nanodust formation around colliding-wind WC binaries, which may present an early source of carbonaceous nanodust in the interstellar medium.
J. Larsson, C. Fransson, B. Sargent, O. C. Jones, M. J. Barlow, P. Bouchet, M. Meixner, J. A. D. L. Blommaert, A. Coulais, O. D. Fox, R. Gastaud, A. Glasse, N. Habel, A. S. Hirschauer, J. Hjorth, J. Jaspers, P. J. Kavanagh, O. Krause, R. M. Lau, L. Lenkic, et al (6) We present initial results from JWST NIRSpec integral field unit observations of the nearby Supernova (SN) 1987A. The observations provide the first spatially-resolved spectroscopy of the ejecta and equatorial ring (ER) over the 1-5 \mu m range. We construct 3D emissivity maps of the [Fe I] 1.443 \mu m line from the inner ejecta and the He I 1.083 \mu m line from the reverse shock (RS), where the former probes the explosion geometry and the latter traces the structure of the circumstellar medium. We also present a model for the integrated spectrum of the ejecta. The [Fe I] 3D map reveals a highly-asymmetric morphology resembling a broken dipole, dominated by two large clumps with velocities of ~2300 km/s. We also find evidence that the Fe-rich inner ejecta have started to interact with the RS. The RS surface traced by the He I line extends from just inside the ER to higher latitudes on both sides of the ER with a half-opening angle ~45 degrees, forming a bubble-like structure. The spectral model for the ejecta allows us to identify the many emission lines, including numerous H_2 lines. We find that the H_2 is most likely excited by far-UV emission, while the metal lines ratios are consistent with a combination of collisional excitation and recombination in the low-temperature ejecta. We also find several high-ionization coronal lines from the ER, requiring a temperature > 2 \times 10^6 K.
Melissa Shahbandeh, Arkaprabha Sarangi, Tea Temim, Tamas Szalai, Ori D. Fox, Samaporn Tinyanont, Eli Dwek, Luc Dessart, Alexei V. Filippenko, Thomas G. Brink, Ryan J. Foley, Jacob Jencson, Justin Pierel, Szanna Zsiros, Armin Rest, WeiKang Zheng, Jennifer Andrews, Geoffrey C. Clayton, Kishalay De, Michael Engesser, et al (19) Supernova (SN) explosions have been sought for decades as a possible source of dust in the Universe, providing the seeds of galaxies, stars, and planetary systems. SN 1987A offers one of the most promising examples of significant SN dust formation, but until the James Webb Space Telescope (JWST), instruments have traditionally lacked the sensitivity at both late times (>1 yr post-explosion) and longer wavelengths (i.e., >10 um) to detect analogous dust reservoirs. Here we present JWST/MIRI observations of two historic Type IIP SNe, 2004et and SN 2017eaw, at nearly 18 and 5 yr post-explosion, respectively. We fit the spectral energy distributions as functions of dust mass and temperature, from which we are able to constrain the dust geometry, origin, and heating mechanism. We place a 90% confidence lower limit on the dust masses for SNe 2004et and 2017eaw of >0.014 and >4e-4 M_sun, respectively. More dust may exist at even colder temperatures or may be obscured by high optical depths. We conclude dust formation in the ejecta to be the most plausible and consistent scenario. The observed dust is radiatively heated to ~100-150 K by ongoing shock interaction with the circumstellar medium. Regardless of the best fit or heating mechanism adopted, the inferred dust mass for SN 2004et is the second highest (next to SN 1987A) inferred dust mass in extragalactic SNe thus far, promoting the prospect of SNe as potential significant sources of dust in the Universe.
Mojgan Aghakhanloo, Nathan Smith, Peter Milne, Jennifer E. Andrews, Schuyler D. Van Dyk, Alexei V. Filippenko, Jacob E. Jencson, Ryan M. Lau, David J. Sand, Samuel Wyatt, WeiKang Zheng We present the first photometric analysis of the supernova (SN) impostor AT 2016blu in NGC 4559. This transient was discovered by the Lick Observatory Supernova Search in 2012 and has continued its outbursts since then. Optical and infrared photometry of AT 2016blu reveals at least 19 outbursts in 2012-2022. Similar photometry from 1999-2009 shows no outbursts, indicating that the star was relatively stable in the decade before discovery. Archival \it Hubble Space Telescope observations suggest that the progenitor had a minimum initial mass of $M >= 33$ M$_{\odot}$ and a luminosity of $L >= 10^{5.7}$ L$_{\odot}$. AT 2016blu's outbursts show irregular variability with multiple closely spaced peaks having typical amplitudes of 1-2 mag and durations of 1-4 weeks. While individual outbursts have irregular light curves, concentrations of these peaks recur with a period of $\sim 113 \pm 2$ d. Based on this period, we predict times for upcoming outbursts in 2023 and 2024. AT 2016blu shares similarities with SN 2000ch in NGC 3432, where outbursts may arise from periastron encounters in an eccentric binary containing a luminous blue variable (LBV). We propose that AT 2016blu's outbursts are also driven by interactions that intensify around periastron in an eccentric system. Intrinsic variability of the LBV-like primary star may cause different intensity and duration of binary interaction at each periastron passage. AT 2016blu also resembles the periastron encounters of $\eta$ Carinae prior to its Great Eruption and the erratic pre-SN eruptions of SN 2009ip. This similarity and the onset of eruptions in the past decade hint that AT 2016blu may also be headed for a catastrophe, making it a target of great interest.
The Wolf-Rayet (WR) binary system WR140 is a close (0.9-16.7 mas) binary star consisting of an O5 primary and WC7 companion and is known as the archetype of episodic dust-producing WRs. Dust in WR binaries is known to form in a confined stream originating from the collision of the two stellar winds, with orbital motion of the binary sculpting the large-scale dust structure into arcs as dust is swept radially outwards. It is understood that sensitive conditions required for dust production in WR140 are only met around periastron when the two stars are sufficiently close. Here we present multiepoch imagery of the circumstellar dust shell of WR140. We constructed geometric models that closely trace the expansion of the intricately structured dust plume, showing that complex effects induced by orbital modulation may result in a 'Goldilocks zone' for dust production. We find that the expansion of the dust plume cannot be reproduced under the assumption of a simple uniform-speed outflow, finding instead the dust to be accelerating. This constitutes a direct kinematic record of dust motion under acceleration by radiation pressure and further highlights the complexity of the physical conditions in colliding-wind binaries.
Ryan M. Lau, Matthew J. Hankins, Yinuo Han, Ioannis Argyriou, Michael F. Corcoran, Jan J. Eldridge, Izumi Endo, Ori D. Fox, Macarena Garcia Marin, Theodore R. Gull, Olivia C. Jones, Kenji Hamaguchi, Astrid Lamberts, David R. Law, Thomas Madura, Sergey V. Marchenko, Hideo Matsuhara, Anthony F. J. Moffat, Mark R. Morris, Patrick W. Morris, et al (12) Massive colliding-wind binaries that host a Wolf-Rayet (WR) star present a potentially important source of dust and chemical enrichment in the interstellar medium (ISM). However, the chemical composition and survival of dust formed from such systems is not well understood. The carbon-rich WR (WC) binary WR~140 presents an ideal astrophysical laboratory for investigating these questions given its well-defined orbital period and predictable dust-formation episodes every 7.93 years around periastron passage. We present observations from our Early Release Science program (ERS1349) with the James Webb Space Telescope (JWST) Mid-Infrared Instrument (MIRI) Medium-Resolution Spectrometer (MRS) and Imager that reveal the spectral and spatial signatures of nested circumstellar dust shells around WR~140. MIRI MRS spectroscopy of the second dust shell and Imager detections of over 17 shells formed throughout the past $\gtrsim130$ years confirm the survival of carbonaceous dust grains from WR~140 that are likely carriers of "unidentified infrared" (UIR)-band features at 6.4 and 7.7 $\mu$m. The observations indicate that dust-forming WC binaries can enrich the ISM with organic compounds and carbonaceous dust.
A handful of binary Wolf-Rayet stars are known to harbour spectacular spiral structures spanning a few hundred AU. These systems host some of the highest dust production rates in the Universe and are therefore interesting candidates to address the origin of the enigmatic dust excess observed across galactic evolution. The substantial interaction between the winds of the Wolf-Rayet star and its companion constitutes a unique laboratory to study the mechanisms of dust nucleation in a hostile environment. Using the grid-based $\texttt{RAMSES}$ code, we investigate this problem by performing a 3D hydrodynamic simulation of the inner region of the prototypical spiral nebula around WR104. We then process the $\texttt{RAMSES}$ results using the radiative transfer code $\texttt{RADMC3d}$ to generate a candidate observable scene. This allows us to estimate the geometrical parameters of the shocked region. We link those quantities to the specific chemical pathway for dust nucleation, where the hydrogen-rich companion's wind catalyses dust formation. The scaling laws we derive constitute a unique tool that can be directly compared to observations. Depending on the dust nucleation locus, the velocity field reveals a differential wind speed. Thus, the initial dust speed could be more balanced between the speeds of the two stellar winds ($\sim$1600 km/s). With $\texttt{RADMC3d}$, we provide constraints on the dust nucleation radius for different combinations of dust-to-gas ratio, hydrogen enrichment and dust grain properties. Finally, our models reveal that dust may escape beyond the boundaries of the spiral due to hydrodynamical instabilities in the wind collision zone.
Kishalay De, Ilya Mereminskiy, Roberto Soria, Charlie Conroy, Erin Kara, Shreya Anand, Michael C. B. Ashley, Martha L. Boyer, Deepto Chakrabarty, Brian Grefenstette, Matthew J. Hankins, Lynne A. Hillenbrand, Jacob E. Jencson, Viraj Karambelkar, Mansi M. Kasliwal, Ryan M. Lau, Alexander Lutovinov, Anna M. Moore, Mason Ng, Christos Panagiotou, et al (7) We present the discovery and multi-wavelength characterization of SRGA J181414.6-225604, a Galactic hard X-ray transient discovered during the ongoing SRG/ART-XC sky survey. Using data from the Palomar Gattini-IR survey, we identify a spatially and temporally coincident variable infrared (IR) source, IRAS 18111-2257, and classify it as a very late-type (M7-M8), long period ($1502 \pm 24$ days) and luminous ($M_K\approx -9.9 \pm 0.2$) O-rich Mira donor star located at a distance of $\approx 14.6^{+2.9}_{-2.3}$ kpc. Combining multi-color photometric data over the last $\approx 25$ years, we show that the IR counterpart underwent a recent (starting $\approx 800$ days before the X-ray flare) enhanced mass loss (reaching $\approx 2.1 \times 10^{-5}$ M$_\odot$ yr$^{-1}$) episode resulting in an expanding dust shell obscuring the underlying star. Multi-epoch follow-up from Swift, NICER and NuSTAR reveal a $\approx 200$ day long X-ray outburst reaching a peak luminosity of $L_X \approx 2.5 \times 10^{36}$ erg s$^{-1}$, characterized by a heavily absorbed ($N_{\rm H} \approx 6\times 10^{22}$ cm$^{-2}$) X-ray spectrum consistent with an optically thick Comptonized plasma. The X-ray spectral and timing behavior suggest the presence of clumpy wind accretion together with a dense ionized nebula overabundant in silicate material surrounding the compact object. Together, we show that SRGA J181414.6-225604 is a new symbiotic X-ray binary in outburst, triggered by an intense dust formation episode of a highly evolved donor. Our results offer the first direct confirmation for the speculated connection between enhanced late-stage donor mass loss and active lifetimes of the symbiotic X-ray binaries.
We present the detection of a broad 8 $\mu$m feature in newly formed dust around the carbon-rich Wolf-Rayet (WC) binary WR 125 from N-band low-resolution (NL; R$\sim$250) spectroscopy between 7.3-13.6 $\mu$m and N-band (11.7 $\mu$m) and Q-band (18.8 $\mu$m) imaging with Subaru/COMICS in 2019 October. WR 125 is a colliding wind binary (${\rm WC7+O9}$) that exhibited renewed dust formation starting in 2018, $\sim$28 years after its first dust formation episode had been observed. We also compare our infrared photometry with historical observations and revise the dust-formation period of WR 125 to 28.1 years. Archival infrared spectra of five dusty WC stars, WR 48a, WR 98a, WR 104, WR 112 and WR 118, obtained with ISO/SWS are reanalyzed and compared with the WR 125 spectrum to search for a similar feature. We analyze the dusty WC spectra using two different extinction curves to investigate the impact of interstellar extinction correction on the presence and/or properties of the 8 $\mu$m feature. All of the dusty WC spectra dereddened with the two different extinction curves show a broad feature around 8 $\mu$m (FWHM$\sim$1-2 $\mu$m). We suggest that these 8 $\mu$m features seen in the dusty WC spectra are related to the Class C unidentified infrared (UIR) features.
Kishalay De, Deepto Chakrabarty, Roberto Soria, Michael C. B. Ashley, Charlie Conroy, Matthew J. Hankins, Mansi M. Kasliwal, Ryan M. Lau, Anna M. Moore, Robert Simcoe, Jamie Soon, Tony Travouillon The symbiotic X-ray binary Sct X-1 was suggested as the first known neutron star accreting from a red supergiant companion. Although known for nearly 50 years, detailed characterization of the donor remains lacking, particularly due to the extremely high reddening towards the source ($A_V\gtrsim25$ mag). Here, we present i) improved localization of the counterpart using Gaia and Chandra observations, ii) the first broadband infrared spectrum ($\approx1-5\,\mu$m; $R\approx 2000$) obtained with SpeX on the NASA Infrared Telescope Facility and iii) $J$-band light curve from the Palomar Gattini-IR survey. The infrared spectrum is characterized by i) deep water absorption features (H$_2$O index $\approx 40$%), ii) strong TiO, VO and CO features, and iii) weak/absent CN lines. We show that these features are inconsistent with known red supergiants, but suggest a M8-9 III type O-rich Mira donor star. We report the discovery of large amplitude ($\Delta J\approx3.5$ mag) periodic photometric variability suggesting a pulsation period of $621\pm36\,{\rm(systematic)}\pm8\,{\rm(statistical)}$ days, which we use to constrain the donor to be a relatively luminous Mira ($M_K=-8.6\pm0.3$ mag) at a distance of $3.6^{+0.8}_{-0.7}$ kpc. Comparing these characteristics to recent models, we find the donor to be consistent with a $\approx 3-5$ M$_\odot$ star at an age of $\approx 0.1-0.3$ Gyr. Together, we show that Sct X-1 was previously mis-classified as an evolved High Mass X-ray Binary; instead it is an intermediate mass system with the first confirmed Mira donor in an X-ray binary. We discuss the implications of Mira donors in symbiotic X-ray binaries, and highlight the potential of wide-field infrared time domain surveys and broadband infrared spectroscopy to unveil their demographics.
H. Schatz, Z. Meisel, E. F. Brown, S. S. Gupta, G. W. Hitt, W. R. Hix, R. Jain, R. Lau, P. Möller, W.-J. Ong, P. S. Shternin, Y. Xu, M. Wiescher Nuclear reactions heat and cool the crust of accreting neutron stars and need to be understood to interpret observations of X-ray bursts and of long-term cooling in transiently accreting systems. It was recently suggested that previously neglected neutron transfer reactions may play a significant role in the nuclear processes. We present results from full nuclear network calculations that now include these reactions and determine their impact on crust composition, crust impurity, heating, and cooling. We find that a large number of neutron transfer reactions indeed occur and impact crust models. In particular, we identify a new type of reaction cycle that brings a pair of nuclei across the nuclear chart into equilibrium via alternating neutron capture and neutron release, interspersed with a neutron transfer. While neutron transfer reactions lead to changes in crust model predictions, and need to be considered in future studies, previous conclusions concerning heating, cooling, and compositional evolution are remarkably robust.
Kylie Y. Hansen, Kishalay De, Michael C. B. Ashley, Mansi M. Kasliwal, Alexander Delacroix, Tim Greffe, David Hale, Matthew J. Hankins, Ryan Lau, Chengkui Li, Daniel McKenna, Anna M. Moore, Eran O. Ofek, Roger M. Smith, Jamie Soon, Roberto Soria, Gokul P. Srinivasaragavan, Tony Travouillon We report second-timescale infrared photometry of the nova V1674 Her using Palomar Gattini-IR. These observations constitute the first infrared and highest temporal resolution data (resolution of ~ 0.84 s) of the nova reported to date. PGIR observed in this fast readout mode for more than an hour on three nights between 3 and 6 days after discovery. We searched for periodic variability using a Lomb-Scargle periodogram and did not detect anything down to a three sigma upper limit of 0.093 mag. This suggests that the periodic variability detected in the nova by Patterson et al. 2021 was lower by at least a factor of about 1.65 in the first week of the eruption.
G. P. Srinivasaragavan, I. Sfaradi, J. Jencson, K. De, A. Horesh, M. M. Kasliwal, S. Tinyanont, M. Hankins, S. Schulze, M. C. B. Ashley, M. J. Graham, V. Karambelkar, R. Lau, A. A. Mahabal, A. M. Moore, E. O. Ofek, Y. Sharma, J. Sollerman, J. Soon, R. Soria, et al (2) We present a detailed analysis of SN 2020qmp, a nearby type IIP core-collapse supernova (CCSN), discovered by the Palomar Gattini-IR (PGIR) survey in the galaxy UGC07125. We illustrate how the multiwavelength study of this event helps our general understanding of stellar progenitors and circumstellar medium (CSM) interactions in CCSNe. We also highlight the importance of near-infrared (NIR) surveys for early detections of SNe in dusty environments. SN 2020qmp displays characteristic hydrogen lines in its optical spectra, as well as a plateau in its optical LC, hallmarks of a type IIP SN. We do not detect linear polarization during the plateau phase, with a 3 sigma upper limit of 0.78%. Through hydrodynamical LC modeling and an analysis of its nebular spectra, we estimate a progenitor mass of around 11 solar masses, and an explosion energy of around 0.8e51 erg. We find that the spectral energy distribution cannot be explained by a simple CSM interaction model, assuming a constant shock velocity and steady mass-loss rate, and the excess X-ray luminosity compared with the synchrotron radio luminosity suggests deviations from equipartition. Finally, we demonstrate the advantages of NIR surveys over optical surveys for the detection of dust-obscured CCSNe in the local universe. Specifically, our simulations show that the Wide-Field Infrared Transient Explorer will detect about 14 more CCSNe out of 75 expected in its footprint within 40 Mpc, over five years than an optical survey equivalent to the Zwicky Transient Facility would detect. We have determined or constrained the main properties of SN 2020qmp and of its progenitor, highlighting the value of multiwavelength follow-up observations of nearby CCSNe. We have also shown that forthcoming NIR surveys will finally enable us to do a nearly complete census of CCSNe in the local universe.
Ryan M. Lau, Samaporn Tinyanont, Matthew J. Hankins, Michael C. B. Ashley, Kishalay De, Alexei V. Filippenko, Lynne A. Hillenbrand, Mansi M. Kasliwal, Jon C. Mauerhan, Anthony F. J. Moffat, Anna M. Moore, Nathan Smith, Jamie Soon, Roberto Soria, Tony Travouillon, Karel A. van der Hucht, Peredur M. Williams, WeiKang Zheng We present optical and infrared (IR) light curves of NaSt1, also known as Wolf-Rayet (WR) 122, with observations from Palomar Gattini-IR (PGIR), the Zwicky Transient Facility (ZTF), the Katzman Automatic Imaging Telescope (KAIT), the Asteroid Terrestrial-impact Last Alert System (ATLAS), and the All-Sky Automated Survey for Supernovae (ASAS-SN). We identify a $P=309.7\pm0.7$ d photometric period from the optical and IR light curves that reveal periodic, sinusoidal variability between 2014 July and 2021 July. We also present historical IR light curves taken between 1983 July and 1989 May, which show variability consistent with the period of the present-day light curves. In the past, NaSt1 was brighter in the $J$ band with larger variability amplitudes than the present-day PGIR values, suggesting that NaSt1 exhibits variability on longer ($\gtrsim$ decade) timescales. Sinusoidal fits to the recent optical and IR light curves show that the amplitude of NaSt1's variability differs at various wavelengths and also reveal significant phase offsets of $17.0\pm2.5$ d between the ZTF $r$ and PGIR $J$ light curves. We interpret the $310$ d photometric period from NaSt1 as the orbital period of an enshrouded massive binary. We suggest that the photometric variability of NaSt1 may arise from variations in the line-of-sight optical depth toward circumstellar optical/IR emitting regions throughout its orbit due to colliding-wind dust formation. We speculate that past mass transfer in NaSt1 may have been triggered by Roche-lobe overflow (RLOF) during an eruptive phase of an Ofpe/WN9 star. Lastly, we argue that NaSt1 is no longer undergoing RLOF mass transfer.
Jacob E. Jencson, Jennifer E. Andrews, Howard E. Bond, Viraj Karambelkar, David J. Sand, Schuyler D. van Dyk, Nadejda Blagorodnova, Martha L. Boyer, Mansi M. Kasliwal, Ryan M. Lau, Shazrene Mohamed, Robert Williams, Patricia A. Whitelock, Rachael C. Amaro, K. Azalee Bostroem, Yize Dong, Michael J. Lundquist, Stefano Valenti, Samuel D. Wyatt, Jamie Burke, et al (22) Nova eruptions, thermonuclear explosions on the surfaces of white dwarfs (WDs), are now recognized to be among the most common shock-powered astrophysical transients. We present the early discovery and rapid ultraviolet (UV), optical, and infrared (IR) temporal development of AT 2019qyl, a recent nova in the nearby Sculptor Group galaxy NGC 300. The light curve shows a rapid rise lasting $\lesssim 1$ day, reaching a peak absolute magnitude of $M_V = -9.2$ mag, and a very fast decline, fading by 2 mag over 3.5 days. A steep dropoff in the light curves after 71 days and the rapid decline timescale suggest a low-mass ejection from a massive WD with $M_{\rm WD} \gtrsim 1.2~M_{\odot}$. We present an unprecedented view of the early spectroscopic evolution of such an event. Three spectra prior to the peak reveal a complex, multicomponent outflow giving rise to internal collisions and shocks in the ejecta of an He/N-class nova. We identify a coincident IR-variable counterpart in the extensive preeruption coverage of the transient location and infer the presence of a symbiotic progenitor system with an O-rich asymptotic-giant-branch donor star, as well as evidence for an earlier UV-bright outburst in 2014. We suggest that AT 2019qyl is analogous to the subset of Galactic recurrent novae with red-giant companions such as RS Oph and other embedded nova systems like V407 Cyg. Our observations provide new evidence that internal shocks between multiple, distinct outflow components likely contribute to the generation of the shock-powered emission from such systems.
Nadejda Blagorodnova, Jakub Klencki, Ondrej Pejcha, Paul M. Vreeswijk, Howard E. Bond, Kevin B. Burdge, Kishalay De, Christoffer Fremling, Robert D. Gehrz, Jacob E. Jencson, Mansi M. Kasliwal, Thomas Kupfer, Ryan M. Lau, Frank J. Masci, R. Michael Rich Luminous Red Novae (LRNe) are astrophysical transients associated with the partial ejection of a binary system's common envelope (CE) shortly before its merger. Here we present the results of our photometric and spectroscopic follow-up campaign of AT2018bwo (DLT18x), a LRN discovered in NGC45, and investigate its progenitor system using binary stellar-evolution models. The transient reached a peak magnitude of $M_r=-10.97\pm0.11$ and maintained this brightness during its optical plateau of $t_p = 41\pm5$days. During this phase, it showed a rather stable photospheric temperature of ~3300K and a luminosity of ~$10^{40}$erg/s. The photosphere of AT2018bwo at early times appeared larger and cooler than other similar LRNe, likely due to an extended mass-loss episode before the merger. Towards the end of the plateau, optical spectra showed a reddened continuum with strong molecular absorption bands. The reprocessed emission by the cooling dust was also detected in the mid-infrared bands ~1.5 years after the outburst. Archival Spitzer and Hubble Space Telescope data taken 10-14 years before the transient event suggest a progenitor star with $T_{prog}\sim 6500$K, $R_{prog}\sim 100R_{\odot}$ and $L_{prog}\sim 2\times10^4L_{\odot}$, and an upper limit for optically thin warm (1000 K) dust mass of $M_d<10^{-6}M_{\odot}$. Using stellar binary-evolution models, we determined the properties of binary systems consistent with the progenitor parameter space. For AT2018bwo, we infer a primary mass of 12-16 $M_{\odot}$, which is 9-45% larger than the ~11$M_{\odot}$ obtained using single-star evolution models. The system, consistent with a yellow-supergiant primary, was likely in a stable mass-transfer regime with -2.4<log ($\dot{M}/M_{\odot}$/yr)<-1.2 a decade before the main instability occurred. During the dynamical merger, the system would have ejected 0.15-0.5$M_{\odot}$ with a velocity of ~500 km/s.
Lynne A. Hillenbrand, Kishalay De, Matthew Hankins, Mansi M. Kasliwal, Luisa M. Rebull, Ryan M. Lau, Roc M. Cutri, Michael C.B. Ashley, Viraj R. Karambelkar, Anna M. Moore, T. Travouillon, A.K. Mainzer We report the discovery of a likely outbursting Class I young stellar object, associated with the star-forming region NGC 281-W (distance $\sim 2.8$ kpc). The source is currently seen only at infrared wavelengths, appearing in both the Palomar Gattini InfraRed ($1.2~\mu$m) and the Near Earth Object Widefield Infrared Survey Explorer ($3.4$ and $4.6~\mu$m) photometric time-domain surveys. Recent near-infrared imaging reveals a new, extended scattered light nebula. Recent near-infrared spectroscopy confirms the similarity of PGIR 20dci to FU Ori type sources, based on strong molecular absorption in CO, H$_2$O, and OH, weak absorption in several atomic lines, and a warm wind/outflow as indicated by a P Cygni profile in the HeI 10830 A line. This is a rare case of an FU Ori star with a well-measured long term photometric rise before a sharper outburst, and the second instance of an FU Ori star with a documented two-step brightening in the mid-infrared.
Kishalay De, Mansi M. Kasliwal, Matthew J. Hankins, Jennifer L. Sokoloski, Scott M. Adams, Michael C. B. Ashley, Aliya-Nur Babul, Ashot Bagdasaryan, Alexandre Delacroix, Richard Dekany, Timothee Greffe, David Hale, Jacob E. Jencson, Viraj R. Karambelkar, Ryan M. Lau, Ashish Mahabal, Daniel McKenna, Anna M. Moore, Eran O. Ofek, Manasi Sharma, et al (8) The nova rate in the Milky Way remains largely uncertain, despite its vital importance in constraining models of Galactic chemical evolution as well as understanding progenitor channels for Type Ia supernovae. The rate has been previously estimated in the range of $\approx10-300$ yr$^{-1}$, either based on extrapolations from a handful of very bright optical novae or the nova rates in nearby galaxies; both methods are subject to debatable assumptions. The total discovery rate of optical novae remains much smaller ($\approx5-10$ yr$^{-1}$) than these estimates, even with the advent of all-sky optical time domain surveys. Here, we present a systematic sample of 12 spectroscopically confirmed Galactic novae detected in the first 17 months of Palomar Gattini-IR (PGIR), a wide-field near-infrared time domain survey. Operating in $J$-band ($\approx1.2$ $\mu$m) that is relatively immune to dust extinction, the extinction distribution of the PGIR sample is highly skewed to large extinction values ($> 50$% of events obscured by $A_V\gtrsim5$ mag). Using recent estimates for the distribution of mass and dust in the Galaxy, we show that the observed extinction distribution of the PGIR sample is commensurate with that expected from dust models. The PGIR extinction distribution is inconsistent with that reported in previous optical searches (null hypothesis probability $< 0.01$%), suggesting that a large population of highly obscured novae have been systematically missed in previous optical searches. We perform the first quantitative simulation of a $3\pi$ time domain survey to estimate the Galactic nova rate using PGIR, and derive a rate of $\approx 46.0^{+12.5}_{-12.4}$ yr$^{-1}$. Our results suggest that all-sky near-infrared time-domain surveys are well poised to uncover the Galactic nova population.
Viraj R. Karambelkar, Mansi M. Kasliwal, Patrick Tisserand, Kishalay De, Shreya Anand, Michael C. B. Ashley, Alex Delacroix, Matthew Hankins, Jacob E. Jencson, Ryan M. Lau, Dan McKenna, Anna Moore, Eran O. Ofek, Roger M. Smith, Roberto Soria, Jamie Soon, Samaporn Tinyanont, Tony Travouillon, Yuhan Yao We are undertaking the first systematic infrared (IR) census of R Coronae Borealis (RCB) stars in the Milky Way, beginning with IR light curves from the Palomar Gattini IR (PGIR) survey. PGIR is a 30 cm $J$-band telescope with a 25 deg$^{2}$ camera that is surveying 18000 deg$^{2}$ of the northern sky ($\delta>-28^{o}$) at a cadence of 2 days. We present PGIR light curves for 922 RCB candidates selected from a mid-IR color-based catalog (Tisserand et al. 2020). Of these 922, 149 are promising RCB candidates as they show pulsations or declines similar to RCB stars. Majority of the candidates that are not RCB stars are either long period variables (LPVs) or RV-Tauri stars. We identify IR color-based criteria to better distinguish between RCB stars and LPVs. As part of a pilot spectroscopic run, we obtained NIR spectra for 26 out of the 149 promising candidates and spectroscopically confirm 11 new RCB stars. We detect strong He I $\lambda 10830$ features in spectra of all RCB stars, likely originating within high velocity (200-400 km-s$^{-1}$) winds in their atmospheres. 9 of these RCB stars show $^{12}$C$^{16}$O and $^{12}$C$^{18}$O molecular absorption features, suggesting that they are formed through a white dwarf merger. We detect quasi-periodic pulsations in the light curves of 5 RCB stars. The periods range between 30-125 days and likely originate from the strange-mode instability in these stars. Our pilot run results motivate a dedicated IR spectroscopic campaign to classify all RCB candidates.
We present Spitzer/IRAC observations of dust formation from six extragalactic carbon-rich Wolf-Rayet (WC) binary candidates in low-metallicity (Z $\lesssim0.65$ Z$_\odot$) environments using multi-epoch mid-infrared (IR) imaging data from the SPitzer InfraRed Intensive Transients Survey (SPIRITS). Optical follow-up spectroscopy of SPIRITS~16ln, 19q, 16df, 18hb, and 14apu reveals emission features from C IV $\lambda5801\text{-}12$~and/or the C III-IV $\lambda4650$ He II $\lambda4686$~blend that are consistent with early-type WC stars. We identify SPIRITS~16ln as the variable mid-IR counterpart of the recently discovered colliding-wind WC4+O binary candidate, N604-WRXc, located in the sub-solar metallicity NGC 604 H II~region in M33. We interpret the mid-IR variability from SPIRITS~16ln as a dust-formation episode in an eccentric colliding-wind WC binary. SPIRITS~19q, 16df, 14apu, and 18hb exhibit absolute [3.6] magnitudes exceeding one of most IR-luminous dust-forming WC systems known, WR~104 (M$_\mathrm{[3.6]}\lesssim-12.3$). An analysis of dust formation in the mid-IR outburst from SPIRITS~19q reveals a high dust production rate of $\dot{M}_d\gtrsim2\times10^{-6}$ M$_\odot$ yr$^{-1}$, which may therefore exceed that of the most efficient dust-forming WC systems known. We demonstrate that efficient dust-formation is feasible from early-type WC binaries in the theoretical framework of colliding-wind binary dust formation if the systems host an O-type companion with a high mass-loss rate ($\dot{M}\gtrsim1.6\times10^{-6}$ M$_\odot$ yr$^{-1}$). This efficient dust-formation from early-type WC binaries highlights their potential role as significant sources of dust in low-metallicity environments.
The recent discovery of a spectacular dust plume in the system 2XMM J160050.7-514245 (referred to as "Apep") suggested a physical origin in a colliding-wind binary by way of the "Pinwheel" mechanism. Observational data pointed to a hierarchical triple-star system, however several extreme and unexpected physical properties seem to defy the established physics of such objects. Most notably, a stark discrepancy was found in the observed outflow speed of the gas as measured spectroscopically in the line-of-sight direction compared to the proper motion expansion of the dust in the sky plane. This enigmatic behaviour arises at the wind base within the central Wolf-Rayet binary: a system that has so far remained spatially unresolved. Here we present an updated proper motion study deriving the expansion speed of Apep's dust plume over a two-year baseline that is four times slower than the spectroscopic wind speed, confirming and strengthening the previous finding. We also present the results from high-angular-resolution near-infrared imaging studies of the heart of the system, revealing a close binary with properties matching a Wolf-Rayet colliding-wind system. Based on these new observational constraints, an improved geometric model is presented yielding a close match to the data, constraining the orbital parameters of the Wolf-Rayet binary and lending further support to the anisotropic wind model.
Ryan M. Lau, Matthew J. Hankins, Yinuo Han, Izumi Endo, Anthony F. J. Moffat, Michael E. Ressler, Itsuki Sakon, Joel Sanchez-Bermudez, Anthony Soulain, Ian R. Stevens, Peter G. Tuthill, Peredur M. Williams WR 112 is a dust-forming carbon-rich Wolf-Rayet (WC) binary with a dusty circumstellar nebula that exhibits a complex asymmetric morphology, which traces the orbital motion and dust formation in the colliding winds of the central binary. Unraveling the complicated circumstellar dust emission around WR 112 therefore provides an opportunity to understand the dust formation process in colliding-wind WC binaries. In this work, we present a multi-epoch analysis of the circumstellar dust around WR 112 using seven high spatial resolution (FWHM $\sim0.3-0.4''$) N-band ($\lambda \sim12$ $\mu$m) imaging observations spanning almost 20 years and includes newly obtained images from Subaru/COMICS in Oct 2019. In contrast to previous interpretations of a face-on spiral morphology, we observe clear evidence of proper motion of the circumstellar dust around WR 112 consistent with a nearly edge-on spiral with a $\theta_s=55^\circ$ half-opening angle and a $\sim20$-yr period. The revised near edge-on geometry of WR 112 reconciles previous observations of highly variable non-thermal radio emission that was inconsistent with a face-on geometry. We estimate a revised distance to WR 112 of $d = 3.39^{+0.89}_{-0.84}$ kpc based on the observed dust expansion rate and a spectroscopically derived WC terminal wind velocity of $v_\infty= 1230\pm260$ km s$^{-1}$. With the newly derived WR 112 parameters we fit optically-thin dust spectral energy distribution models and determine a dust production rate of $\dot{M}_d=2.7^{+1.0}_{-1.3}\times10^{-6}$ M$_\odot$ yr$^{-1}$, which demonstrates that WR 112 is one of the most prolific dust-making WC systems known.
We present a dust spectral energy distribution (SED) and binary stellar population analysis revisiting the dust production rates (DPRs) in the winds of carbon-rich Wolf-Rayet (WC) binaries and their impact on galactic dust budgets. DustEM SED models of 19 Galactic WC ``dustars" reveal DPRs of $\dot{M}_d\sim10^{-10}-10^{-6}$ M$_\odot$ yr$^{-1}$ and carbon dust condensation fractions, $\chi_C$, between $0.002 - 40\%$. A large ($0.1 - 1.0$ $\mu$m) dust grain size composition is favored for efficient dustars where $\chi_C\gtrsim1\%$. Results for dustars with known orbital periods verify a power-law relation between $\chi_C$, orbital period, WC mass-loss rate, and wind velocity consistent with predictions from theoretical models of dust formation in colliding-wind binaries. We incorporated dust production into Binary Population and Spectral Synthesis (BPASS) models to analyze dust production rates from WC dustars, asymptotic giant branch stars (AGBs), red supergiants (RSGs), and core-collapse supernovae (SNe). BPASS models assuming constant star formation (SF) and a co-eval $10^6$ M$_\odot$ stellar population were performed at low, Large Magellanic Cloud (LMC)-like, and solar metallicities (Z = 0.001, 0.008, and 0.020). Both constant SF and co-eval models show that SNe are net dust destroyers at all metallicities. Constant SF models at LMC-like metallicities show that AGB stars slightly outproduce WC binaries and RSGs by factors of $2-3$, whereas at solar metallicites WC binaries are the dominant source of dust for $\sim60$ Myr until the onset of AGBs, which match the dust input of WC binaries. Co-eval population models show that for "bursty" SF, AGB stars dominate dust production at late times ($t\gtrsim 70$ Myr).
In 2010 May, an intermediate luminosity optical transient was discovered in the nearby galaxy NGC 300 by a South African amateur astronomer. In the decade since its discovery, multi-wavelength observations of the misnamed ``SN 2010da'' have continually re-shaped our understanding of this high mass X-ray binary system. In this review, we present an overview of the multi-wavelength observations and attempts to understand the 2010 transient event and, later, the re-classification of this system as NGC~300 ULX-1: a red supergiant + neutron star ultraluminous X-ray source.
Matthew J. Hankins, Ryan M. Lau, James T. Radomski, Angela S. Cotera, Mark R. Morris, Elisabeth A. C. Mills, Daniel L. Walker, Ashley T. Barnes, Janet P. Simpson, Terry L. Herter, Steven N. Longmore, John Bally, Mansi M. Kasliwal, Nadeen B. Sabha, Macarena Garcia-Marin The Galactic Center contains some of the most extreme conditions for star formation in our Galaxy as well as many other phenomena that are unique to this region. Given our relative proximity to the Galactic Center, we are able to study details of physical processes to a level that is simply not yet possible for more distant galaxies, yielding an otherwise inaccessible view of the nuclear region of a galaxy. We recently carried out a targeted imaging survey of mid-infrared bright portions of the Galactic Center at 25 and 37 $\mu$m using the FORCAST instrument on SOFIA. This survey was one of the inaugural Legacy Programs from SOFIA cycle 7, observing a total area of 403 arcmin$^2$ (2180 pc$^2$), including the Sgr A, B, and C complexes. Here we present an overview of the survey strategy, observations, and data reduction as an accompaniment to the initial public release of the survey data. We discuss interesting regions and features within the data including extended features near the circumnuclear disk, structures in the Arched Filaments and Sickle H II regions, and signs of embedded star formation in Sgr B2 and Sgr C. We also feature a handful of less well studied mid-infrared sources located between Sgr A and Sgr C that could be sites of relatively isolated star formation activity. Last, we discuss plans for subsequent publications and future data releases from the survey.
Kishalay De, Matthew J. Hankins, Mansi M. Kasliwal, Anna M. Moore, Eran O. Ofek, Scott M. Adams, Michael C. B. Ashley, Aliya-Nur Babul, Ashot Bagdasaryan, Kevin B. Burdge, Jill Burnham, Richard G. Dekany, Alexander Declacroix, Antony Galla, Tim Greffe, David Hale, Jacob E. Jencson, Ryan M. Lau, Ashish Mahabal, Daniel McKenna, et al (7) (Abridged) Palomar Gattini-IR is a new wide-field, near-infrared robotic time domain survey operating at Palomar Observatory. Using a 30 cm telescope mounted with a H2RG detector, Gattini-IR achieves a field of view of 25 sq. deg. with a pixel scale of 8.7" in J-band. Here, we describe the system design, survey operations, data processing system and on-sky performance of Palomar Gattini-IR. As a part of the nominal survey, Gattini-IR scans $\approx 7500$ square degrees of the sky every night to a median 5$\sigma$ depth of $15.7$ AB mag outside the Galactic plane. The survey covers $\approx 15000$ square degrees of the sky visible from Palomar with a median cadence of 2 days. A real-time data processing system produces stacked science images from dithered raw images taken on sky, together with PSF-fit source catalogs and transient candidates identified from subtractions within a median delay of $\approx 4$ hours from the time of observation. The calibrated data products achieve an astrometric accuracy (RMS) of $\approx 0.7$" with respect to Gaia DR2 for sources with S/N $> 10$, and better than $\approx 0.35$" for sources brighter than $\approx 12$ Vega mag. The photometric accuracy (RMS) achieved in the PSF-fit source catalogs is better than $\approx 3$% for sources brighter than $\approx 12$ Vega mag, as calibrated against the 2MASS catalog. With a field of view $\approx 40\times$ larger than any other existing near infrared imaging instrument, Gattini-IR is probing the reddest and dustiest transients in the local universe such as dust obscured supernovae in nearby galaxies, novae behind large columns of extinction within the galaxy, reddened micro-lensing events in the Galactic plane and variability from cool and dust obscured stars. We present results from transients and variables identified since the start of the commissioning period.
Samaporn Tinyanont, Ryan M Lau, Mansi M Kasliwal, Keiichi Maeda, Nathan Smith, Ori D Fox, Robert D Gehrz, Kishalay De, Jacob Jencson, John Bally, Frank Masci Supernova (SN) 2014C is a unique explosion where a seemingly typical hydrogen-poor stripped envelope SN started to interact with a dense, hydrogen-rich circumstellar medium (CSM) a few months after the explosion. The delayed interaction suggests a detached CSM shell, unlike in a typical SN IIn where the CSM is much closer and the interaction commences earlier post-explosion; indicating a different mass loss history. We present near- to mid-infrared observations of SN 2014C from 1-5 years after the explosion, including uncommon 9.7 $\mu$m imaging with COMICS on the Subaru telescope. Spectroscopy shows that the interaction is still ongoing, with the intermediate-width He I 1.083 $\mu$m emission present out to our latest epoch 1639 days post-explosion. The last Spitzer/IRAC photometry at 1920 days post-explosion further confirms ongoing CSM interaction. The 1-10 $\mu$m spectral energy distributions (SEDs) can be explained by a dust model with a mixture of 69% carbonaceous and 31% silicate dust, pointing to a chemically inhomogeneous CSM. The inference of silicate dust is the first among interacting SNe. An SED model with purely carbonaceous CSM dust is possible, but would require more than 0.22 $M_{\odot}$ of dust, which is an order of magnitude larger than what observed in any other SNe, measured in the same way, at this epoch. The light curve beyond 500 days is well fit by an interaction model with a wind-driven CSM and a mass loss rate of $\sim 10^{-3} \, M_{\odot}\,\rm yr^{-1}$, which presents an additional CSM density component exterior to the constant density shell reported previously in the literature. SN 2014C could originate in a binary system, similar to RY Scuti, which would explain the observed chemical and density profile inhomogeneity in the CSM.
SN2010da/NGC 300 ULX-1 was first detected as a supernova impostor in May 2010 and was recently discovered to be a pulsating ultraluminous X-ray source. In this letter, we present VLT/X-shooter spectra of this source obtained in October 2018, covering the wavelength range 350-2300 nm. The $J$- and $H$-bands clearly show the presence of a red supergiant donor star that is best matched by a MARCS stellar atmosphere with $T_{\rm eff} = 3650 - 3900$ K and $\log(L_{\rm bol}/L_{\odot}) = 4.25\pm0.10$, which yields a stellar radius $R = 310 \pm 70 R_{\odot}$. To fit the full spectrum, two additional components are required: a blue excess that can be fitted either by a hot blackbody (T $\gtrsim 20,000$ K) or a power law (spectral index $\alpha \approx 4$) and is likely due to X-ray emission reprocessed in the outer accretion disk or the donor star; and a red excess that is well fitted by a blackbody with a temperature of $\sim 1100$ K, and is likely due to warm dust in the vicinity of SN2010da. The presence of a red supergiant in this system implies an orbital period of at least 0.8-2.1 years, assuming Roche lobe overflow. Given the large donor-to-compact object mass ratio, orbital modulations of the radial velocity of the red supergiant are likely undetectable. However, the radial velocity amplitude of the neutron star is large enough (up to 40-60 km s$^{-1}$) to potentially be measured in the future, unless the system is viewed at a very unfavorable inclination.
Ryan M. Lau, Marianne Heida, Dominic J. Walton, Mansi M. Kasliwal, Scott M. Adams, Ann Marie Cody, Kishalay De, Robert D. Gehrz, Felix Furst, Jacob E. Jencson, jamie A. Kennea, Frank Masci We present a mid-infrared (IR) sample study of nearby ultraluminous X-ray sources (ULXs) using multi-epoch observations with the Infrared Array Camera (IRAC) on the Spitzer Space Telescope. Spitzer/IRAC observations taken after 2014 were obtained as part of the Spitzer Infrared Intensive Transients Survey (SPIRITS). Our sample includes 96 ULXs located within 10 Mpc. Of the 96~ULXs, 12 have candidate counterparts consistent with absolute mid-IR magnitudes of supergiants, and 16 counterparts exceeded the mid-IR brightness of single supergiants and are thus more consistent with star clusters or non-ULX background active galactic nuclei (AGN). The supergiant candidate counterparts exhibit a bi-modal color distribution in a Spitzer/IRAC color-magnitude diagram, where "red" and "blue" ULXs fall in IRAC colors $[3.6] - [4.5]\sim0.7$ and $[3.6] - [4.5]\sim0.0$, respectively. The mid-IR colors and absolute magnitudes of 4 "red" and 5 "blue" ULXs are consistent with that of supergiant B[e] (sgB[e]) and red supergiant (RSG) stars, respectively. While "blue", RSG-like mid-IR ULX counterparts likely host RSG mass donors, we propose the "red" counterparts are ULXs exhibiting the "B[e] phenomenon'' rather than hosts of sgB[e] mass donors. We show that the mid-IR excess from the "red" ULXs is likely due to thermal emission from circumstellar or circumbinary dust. Using dust as a probe for total mass, we estimate mass-loss rates of $\dot{M}\sim1\times10^{-4}$ M$_\odot$ yr$^{-1}$ in dust-forming outflows of red ULXs. Based on the transient mid-IR behavior and its relatively flat spectral index, $\alpha=-0.19\pm0.1$, we suggest that the mid-IR emission from Holmberg IX X-1 originates from a variable jet.
Jacob E. Jencson, Scott M. Adams, Howard E. Bond, Schuyler D. van Dyk, Mansi M. Kasliwal, John Bally, Nadejda Blagorodnova, Kishalay De, Christoffer Fremling, Yuhan Yao, Andrew Fruchter, David Rubin, Cristina Barbarino, Jesper Sollerman, Adam A. Miller, Erin K. S. Hicks, Matthew A. Malkan, Igor Andreoni, Eric C. Bellm, Robert Buchheim, et al (32) We present the discovery of an optical transient (OT) in Messier 51, designated M51 OT2019-1 (also ZTF19aadyppr, AT 2019abn, ATLAS19bzl), by the Zwicky Transient Facility (ZTF). The OT rose over 15 days to an observed luminosity of $M_r=-13$ (${\nu}L_{\nu}=9\times10^6~L_{\odot}$), in the luminosity gap between novae and typical supernovae (SNe). Spectra during the outburst show a red continuum, Balmer emission with a velocity width of $\approx400$ km s$^{-1}$, Ca II and [Ca II] emission, and absorption features characteristic of an F-type supergiant. The spectra and multiband light curves are similar to the so-called "SN impostors" and intermediate-luminosity red transients (ILRTs). We directly identify the likely progenitor in archival Spitzer Space Telescope imaging with a $4.5~\mu$m luminosity of $M_{[4.5]}\approx-12.2$ and a $[3.6]-[4.5]$ color redder than 0.74 mag, similar to those of the prototype ILRTs SN 2008S and NGC 300 OT2008-1. Intensive monitoring of M51 with Spitzer further reveals evidence for variability of the progenitor candidate at [4.5] in the years before the OT. The progenitor is not detected in pre-outburst Hubble Space Telescope optical and near-IR images. The optical colors during outburst combined with spectroscopic temperature constraints imply a higher reddening of $E(B-V)\approx0.7$ mag and higher intrinsic luminosity of $M_r\approx-14.9$ (${\nu}L_{\nu}=5.3\times10^7~L_{\odot}$) near peak than seen in previous ILRT candidates. Moreover, the extinction estimate is higher on the rise than on the plateau, suggestive of an extended phase of circumstellar dust destruction. These results, enabled by the early discovery of M51 OT2019-1 and extensive pre-outburst archival coverage, offer new clues about the debated origins of ILRTs and may challenge the hypothesis that they arise from the electron-capture induced collapse of extreme asymptotic giant branch stars.
Mansi M. Kasliwal, Scott Adams, Igor Andreoni, Michael Ashley, Nadia Blagorodnova, Kishalay De, Danielle Frostig, Gabor Furesz, Jacob Jencson, Matthew Hankins, George Helou, Ryan Lau, Anna Moore, Eran Ofek, Rob Simcoe, Jennifer Sokoloski, Jamie Soon, Samaporn Tinyanont, Tony Travouillon Opening up the dynamic infrared sky for systematic time-domain exploration would yield many scientific advances. Multi-messenger pursuits such as localizing gravitational waves from neutron star mergers and quantifying the nucleosynthetic yields require the infrared. Another multi-messenger endeavor that needs infrared surveyors is the study of the much-awaited supernova in our own Milky Way. Understanding shocks in novae, true rates of supernovae and stellar mergers are some other examples of stellar evolution and high energy physics wherein the answers are buried in the infrared. We discuss some of the challenges in the infrared and pathfinders to overcome them. We conclude with recommendations on both infrared discovery engines and infrared follow-up machines that would enable this field to flourish in the next decade.
V. R. Karambelkar, S. M. Adams, P. A. Whitelock, M. M. Kasliwal, J. E. Jencson, M. L. Boyer, S. R. Goldman, F. Masci, A. M. Cody, J. Bally, H. E. Bond, R. D. Gehrz, M. Parthasarathy, R. M. Lau We present a catalog of 417 luminous infrared variable stars with periods exceeding 250 days. These were identified in 20 nearby galaxies by the ongoing SPIRITS survey with the Spitzer Space Telescope. Of these, 359 variables have $M_{[4.5]}$ (phase-weighted mean magnitudes) fainter than $-12$ and periods and luminosities consistent with previously reported variables in the Large Magellanic Cloud. However, 58 variables are more luminous than $M_{[4.5]} = -12$, including 11 that are brighter than $M_{[4.5]} = -13$ with the brightest having $M_{[4.5]} = -15.51$. Most of these bright variable sources have quasi-periods longer than 1000 days, including four over 2000 days. We suggest that the fundamental period-luminosity relationship, previously measured for the Large Magellanic Cloud, extends to much higher luminosities and longer periods in this large galaxy sample. We posit that these variables include massive AGB stars (possibly super-AGB stars), red supergiants experiencing exceptionally high mass-loss rates, and interacting binaries. We also present 3.6, 4.5, 5.8 and 8.0 $\mu$m photometric catalogs for all sources in these 20 galaxies.
Samaporn Tinyanont, Mansi M Kasliwal, Kelsie Krafton, Ryan Lau, Jeonghee Rho, Douglas C Leonard, Kishalay De, Jacob Jencson, Dimitri Mawet, Maxwell Millar-Blanchaer, Ricky Nilsson, Lin Yan, Robert D Gehrz, George Helou, Schuyler D Van Dyk, Eugene Serabyn, Ori D Fox, Geoffrey Clayton We present infrared (IR) photometry and spectroscopy of the Type II-P SN 2017eaw and its progenitor in the nearby galaxy NGC 6946. Progenitor observations in the Ks band in 4 epochs from 1 year to 1 day before the explosion reveal no significant variability in the progenitor star greater than 6% that last longer than 200 days. SN 2017eaw is a typical SN II-P with near-IR and mid-IR photometric evolution similar to those of SNe 2002hh and 2004et, other normal SNe II-P in the same galaxy. Spectroscopic monitoring between 389 and 480 days post explosion reveals strong CO first overtone emission at 389 d, with a line profile matching that of SN 1987A from the same epoch, indicating $\sim 10^{-3} \, M_{\odot}$ of CO at 1,800 K. From the 389 d epoch until the most recent observation at 566 d, the first overtone feature fades while the 4.5 $\mu$m excess, likely from the CO fundamental band, remains. This behavior indicates that the CO has not been destroyed, but that the gas has cooled enough that the levels responsible for first overtone emissions are no longer populated. Finally, the evolution of Spitzer 3.6 $\mu$m photometry shows evidence for dust formation in SN 2017eaw, with a dust mass of $10^{-6}$ or $10^{-4}\,M_{\odot}$ assuming carbonaceous or silicate grains respectively.
Jacob E. Jencson, Mansi M. Kasliwal, Scott M. Adams, Howard E. Bond, Kishalay De, Joel Johansson, Viraj Karambelkar, Ryan M. Lau, Samaporn Tinyanont, Stuart D. Ryder, Ann Marie Cody, Frank J. Masci, John Bally, Nadia Blagorodnova, Sergio Castellón, Christoffer Fremling, Robert D. Gehrz, George Helou, Charles D. Kilpatrick, Peter A. Milne, et al (6) We present a systematic study of the most luminous ($M_{\mathrm{IR}}$ [Vega magnitudes] brighter than $-14$) infrared (IR) transients discovered by the SPitzer InfraRed Intensive Transients Survey (SPIRITS) between 2014 and 2018 in nearby galaxies ($D < 35$ Mpc). The sample consists of nine events that span peak IR luminosities of $M_{[4.5],\mathrm{peak}}$ between $-14$ and $-18.2$, show IR colors between $0.2 < ([3.6]{-}[4.5]) < 3.0$, and fade on timescales between $55$ days $< t_{\mathrm{fade}} < 480$ days. The two reddest events ($A_V > 12$) show multiple, luminous IR outbursts over several years and have directly detected, massive progenitors in archival imaging. With analyses of extensive, multiwavelength follow-up, we suggest the following possible classifications: five obscured core-collapse supernovae (CCSNe), two erupting massive stars, one luminous red nova, and one intermediate-luminosity red transient. We define a control sample of all optically discovered transients recovered in SPIRITS galaxies and satisfying the same selection criteria. The control sample consists of eight CCSNe and one Type Iax SN. We find that 7 of the 13 CCSNe in the SPIRITS sample have lower bounds on their extinction of $2 < A_V < 8$. We estimate a nominal fraction of CCSNe in nearby galaxies that are missed by optical surveys as high as $38.5^{+26.0}_{-21.9}$% (90% confidence). This study suggests that a significant fraction of CCSNe may be heavily obscured by dust and therefore undercounted in the census of nearby CCSNe from optical searches.
Mansi M. Kasliwal, Daniel Kasen, Ryan M. Lau, Daniel A. Perley, Stephan Rosswog, Eran O. Ofek, Kenta Hotokezaka, Ranga-Ram Chary, Jesper Sollerman, Ariel Goobar, David L. Kaplan We report our Spitzer Space Telescope observations and detections of the binary neutron star merger GW170817. At 4.5um, GW170817 is detected at 21.9 mag AB at +43 days and 23.9 mag AB at +74 days after merger. At 3.6um, GW170817 is not detected to a limit of 23.2 mag AB at +43 days and 23.1 mag AB at +74 days. Our detections constitute the latest and reddest constraints on the kilonova/macronova emission and composition of heavy elements. The 4.5um luminosity at this late phase cannot be explained by elements exclusively from the first abundance peak of the r-process. Moreover, the steep decline in the Spitzer band, with a power-law index of 3.4 +/- 0.2, can be explained by a few of the heaviest isotopes in the third abundance peak with half-life around 14 days dominating the luminosity (e.g. 140Ba, 143Pr, 147Nd, 156Eu, 191Os, 223Ra, 225Ra, 233Pa, 234Th) or a model with lower deposition efficiency. This data offers evidence that the heaviest elements in the second and third r-process abundance peak were indeed synthesized. Our conclusion is verified by both analytics and network simulations and robust despite intricacies and uncertainties in the nuclear physics. Future observations with Spitzer and James Webb Space Telescope will further illuminate the relative abundance of the synthesized heavy elements.
We present the first spatially resolved mid-infrared (37.1 $\mu$m) image of the Fomalhaut debris disk. We use PSF fitting and subtraction to distinctly measure the flux from the unresolved component and the debris disk. We measure an infrared excess in the point source of $0.9 \pm 0.2$ Jy, consistent with emission from warm dust in an inner disk structure (Su et al. 2016), and inconsistent with a stellar wind origin. We cannot confirm or rule out the presence of a pileup ring (Su et al. 2016) near the star. In the cold region, the 37 $\mu$m imaging is sensitive to emission from small, blowout grains, which is an excellent probe of the dust production rate from planetesimal collisions. Under the assumptions that the dust grains are icy aggregates and the debris disk is in steady state, this result is consistent with the dust production rates predicted by Kenyon & Bromley (2008) from theoretical models of icy planet formation. We find a dust luminosity of $(7.9 \pm 0.8) \times 10^{-4}$ L$_\odot$ and a dust mass of 8 -- 16 lunar masses, depending on grain porosity, with $\sim 1$ lunar mass in grains with radius 1 $\mu$m -- 1 mm. If the grains are icy and highly porous, meter-sized objects must be invoked to explain the far-IR, submm, and mm emission. If the grains are composed of astronomical silicates, there is a dearth of blowout grains (Pawellek et al. 2014) and the mass loss rate is well below the predicted dust production values.
Christoph Baranec, Mark Chun, Donald Hall, Michael Connelley, Klaus Hodapp, Daniel Huber, Michael Liu, Eugene Magnier, Karen Meech, Marianne Takamiya, Richard Griffiths, Reed Riddle, Richard Dekany, Mansi Kasliwal, Ryan Lau, Nicholas M. Law, Olivier Guyon, Imke de Pater, Mike Wong, Eran Ofek, et al (6) We are building a next-generation laser adaptive optics system, Robo-AO-2, for the UH 2.2-m telescope that will deliver robotic, diffraction-limited observations at visible and near-infrared wavelengths in unprecedented numbers. The superior Maunakea observing site, expanded spectral range and rapid response to high-priority events represent a significant advance over the prototype. Robo-AO-2 will include a new reconfigurable natural guide star sensor for exquisite wavefront correction on bright targets and the demonstration of potentially transformative hybrid AO techniques that promise to extend the faintness limit on current and future exoplanet adaptive optics systems.
D. J. Walton, M. Bachetti, F. Fuerst, D. Barret, M. Brightman, A. C. Fabian, B. W. Grefenstette, F. A. Harrison, M. Heida, J. Kennea, P. Kosec, R. M. Lau, K. K. Madsen, M. J. Middleton, C. Pinto, J. F. Steiner, N. Webb Based on phase-resolved broadband spectroscopy using $XMM$-$Newton$ and $NuSTAR$, we report on a potential cyclotron resonant scattering feature at $E \sim 13$ keV in the pulsed spectrum of the recently discoverd ULX pulsar NGC 300 ULX1. If this interpretation is correct, the implied magnetic field of the central neutron star is $B \sim 10^{12}$ G (assuming scattering off electrons), similar to that estimated from the observed spin-up of the star, and also similar to known Galactic X-ray pulsars. We discuss the implications of this result for the connection between NGC 300 ULX1 and the other known ULX pulsars, particularly in light of the recent discovery of a likely proton Cyclotron line in another ULX, M51 ULX-8.
R. Lau, M. Beard, S. S. Gupta, H. Schatz, A. V. Afanasjev, E. F. Brown, A. Deibel, L. R. Gasques, G. W. Hitt, W. R. Hix, L. Keek, P. Möller, P. S. Shternin, A. Steiner, M. Wiescher, Y. Xu X-ray observations of transiently accreting neutron stars during quiescence provide information about the structure of neutron star crusts and the properties of dense matter. Interpretation of the observational data requires an understanding of the nuclear reactions that heat and cool the crust during accretion, and define its nonequilibrium composition. We identify here in detail the typical nuclear reaction sequences down to a depth in the inner crust where the mass density is 2E12 g/cm^3 using a full nuclear reaction network for a range of initial compositions. The reaction sequences differ substantially from previous work. We find a robust reduction of crust impurity at the transition to the inner crust regardless of initial composition, though shell effects can delay the formation of a pure crust somewhat to densities beyond 2E12 g/cm^3. This naturally explains the small inner crust impurity inferred from observations of a broad range of systems. The exception are initial compositions with A >= 102 nuclei, where the inner crust remains impure with an impurity parameter of Qimp~20 due to the N = 82 shell closure. In agreement with previous work we find that nuclear heating is relatively robust and independent of initial composition, while cooling via nuclear Urca cycles in the outer crust depends strongly on initial composition. This work forms a basis for future studies of the sensitivity of crust models to nuclear physics and provides profiles of composition for realistic crust models.
Jacob E. Jencson, Mansi M. Kasliwal, Scott M. Adams, Howard E. Bond, Ryan M. Lau, Joel Johansson, Assaf Horesh, Kunal P. Mooley, Robert Fender, Kishalay De, Dónal O'Sullivan, Frank J. Masci, Ann Marie Cody, Nadia Blagorodnova, Ori D. Fox, Robert D. Gehrz, Peter A. Milne, Daniel A. Perley, Nathan Smith, Schuyler D. Van Dyk We present the discovery by the SPitzer InfraRed Intensive Transients Survey (SPIRITS) of a likely supernova (SN) in NGC 3556 at only 8.8 Mpc, which, despite its proximity, was not detected by optical searches. A luminous infrared (IR) transient at $M_{[4.5]} = -16.7$ mag (Vega), SPIRITS 16tn is coincident with a dust lane in the inclined, star-forming disk of the host. Using IR, optical, and radio observations, we attempt to determine the nature of this event. We estimate $A_V \approx$ 8 - 9 mag of extinction, placing it among the three most highly obscured IR-discovered SNe to date. The [4.5] light curve declined at a rate of 0.013 mag day$^{-1}$, and the $[3.6] - [4.5]$ color grew redder from 0.7 to $\gtrsim$ 1.0 mag by 184.7 days post discovery. Optical/IR spectroscopy shows a red continuum, but no clearly discernible features, preventing a definitive spectroscopic classification. Deep radio observations constrain the radio luminosity of SPIRITS 16tn to $L_{\nu} \lesssim 10^{24}$ erg s$^{-1}$ Hz$^{-1}$ between 3 - 15 GHz, excluding many varieties of radio core-collapse SNe. A type Ia SN is ruled out by the observed red IR color, and lack of features normally attributed to Fe-peak elements in the optical and IR spectra. SPIRITS 16tn was fainter at [4.5] than typical stripped-envelope SNe by $\approx$ 1 mag. Comparison of the spectral energy distribution to SNe II suggests SPIRITS 16tn was both highly obscured, and intrinsically dim, possibly akin to the low-luminosity SN 2005cs. We infer the presence of an IR dust echo powered by a peak luminosity of the transient of $5 \times 10^{40}$ erg s$^{-1} < L_{\mathrm{peak}} < 4\times10^{43}$ erg s$^{-1}$, consistent with the observed range for SNe II. This discovery illustrates the power of IR surveys to overcome the compounding effects of visible extinction and optically sub-luminous events in completing the inventory of nearby SNe.
We present an optical and infrared (IR) study of IC 10 X-2, a high-mass X-ray binary in the galaxy IC 10. Previous optical and X-ray studies suggest X-2 is a Supergiant Fast X-ray Transient: a large-amplitude (factor of $\sim$ 100), short-duration (hours to weeks) X-ray outburst on 2010 May 21. We analyze R- and g-band light curves of X-2 from the intermediate Palomar Transient Factory taken between 2013 July 15 and 2017 Feb 14 show high-amplitude ($\gtrsim$ 1 mag), short-duration ($\lesssim8$ d) flares and dips ($\gtrsim$ 0.5 mag). Near-IR spectroscopy of X-2 from Palomar/TripleSpec show He I, Paschen-$\gamma$, and Paschen-$\beta$ emission lines with similar shapes and amplitudes as those of luminous blue variables (LBVs) and LBV candidates (LBVc). Mid-IR colors and magnitudes from Spitzer/IRAC photometry of X-2 resemble those of known LBV/LBVcs. We suggest that the stellar companion in X-2 is an LBV/LBVc and discuss possible origins of the optical flares. Dips in the optical light curve are indicative of eclipses from optically thick clumps formed in the winds of the stellar counterpart. Given the constraints on the flare duration ($0.02 - 0.8$ d) and the time between flares ($15.1\pm7.8$ d), we estimate the clump volume filling factor in the stellar winds, $f_V$, to be $0.01 < f_V < 0.71$, which overlaps with values measured from massive star winds. In X-2, we interpret the origin of the optical flares as the accretion of clumps formed in the winds of an LBV/LBVc onto the compact object.
The asymptotic giant branch (AGB) star R Sculptoris (R Scl) is one of the most extensively studied stars on the AGB. R Scl is a carbon star with a massive circumstellar shell ($M_{shell}\sim 7.3\times10^{-3}~M_{\odot}$) which is thought to have been produced during a thermal pulse event $\sim2200$ years ago. To study the thermal dust emission associated with its circumstellar material, observations were taken with the Faint Object InfraRed CAMera for the SOFIA Telescope (FORCAST) at 19.7, 25.2, 31.5, 34.8, and 37.1 $\mu$m. Maps of the infrared emission at these wavelengths were used to study the morphology and temperature structure of the spatially extended dust emission. Using the radiative transfer code DUSTY and fitting the spatial profile of the emission, we find that a geometrically thin dust shell cannot reproduce the observed spatially resolved emission. Instead, a second dust component in addition to the shell is needed to reproduce the observed emission. This component, which lies interior to the dust shell, traces the circumstellar envelope of R Scl. It is best fit by a density profile with $n \propto r^{\alpha}$ where $\alpha=0.75^{+0.45}_{-0.25}$ and dust mass of $M_d=9.0^{+2.3}_{-4.1}\times10^{-6}~M_{\odot}$. The strong departure from an $r^{-2}$ law indicates that the mass-loss rate of R Scl has not been constant. This result is consistent with a slow decline in the post-pulse mass-loss which has been inferred from observations of the molecular gas.
M. M. Kasliwal, E. Nakar, L. P. Singer, D. L. Kaplan, D. O. Cook, A. Van Sistine, R. M. Lau, C. Fremling, O. Gottlieb, J. E. Jencson, S.M. Adams, U. Feindt, K. Hotokezaka, S. Ghosh, D. A. Perley, P.-C. Yu, T. Piran, J. R. Allison, G. C. Anupama, A.Balasubramanian, et al (61) Merging neutron stars offer an exquisite laboratory for simultaneously studying strong-field gravity and matter in extreme environments. We establish the physical association of an electromagnetic counterpart EM170817 to gravitational waves (GW170817) detected from merging neutron stars. By synthesizing a panchromatic dataset, we demonstrate that merging neutron stars are a long-sought production site forging heavy elements by r-process nucleosynthesis. The weak gamma-rays seen in EM170817 are dissimilar to classical short gamma-ray bursts with ultra-relativistic jets. Instead, we suggest that breakout of a wide-angle, mildly-relativistic cocoon engulfing the jet elegantly explains the low-luminosity gamma-rays, the high-luminosity ultraviolet-optical-infrared and the delayed radio/X-ray emission. We posit that all merging neutron stars may lead to a wide-angle cocoon breakout; sometimes accompanied by a successful jet and sometimes a choked jet.
I. Andreoni, K. Ackley, J. Cooke, A. Acharyya, J. R. Allison, G. E. Anderson, M. C. B. Ashley, D. Baade, M. Bailes, K. Bannister, A. Beardsley, M. S. Bessell, F. Bian, P. A. Bland, M. Boer, T. Booler, A. Brandeker, I. S. Brown, D. Buckley, S.-W. Chang, et al (105) The discovery of the first electromagnetic counterpart to a gravitational wave signal has generated follow-up observations by over 50 facilities world-wide, ushering in the new era of multi-messenger astronomy. In this paper, we present follow-up observations of the gravitational wave event GW170817 and its electromagnetic counterpart SSS17a/DLT17ck (IAU label AT2017gfo) by 14 Australian telescopes and partner observatories as part of Australian-based and Australian-led research programs. We report early- to late-time multi-wavelength observations, including optical imaging and spectroscopy, mid-infrared imaging, radio imaging, and searches for fast radio bursts. Our optical spectra reveal that the transient source afterglow cooled from approximately 6400K to 2100K over a 7-day period and produced no significant optical emission lines. The spectral profiles, cooling rate, and photometric light curves are consistent with the expected outburst and subsequent processes of a binary neutron star merger. Star formation in the host galaxy probably ceased at least a Gyr ago, although there is evidence for a galaxy merger. Binary pulsars with short (100 Myr) decay times are therefore unlikely progenitors, but pulsars like PSR B1534+12 with its 2.7 Gyr coalescence time could produce such a merger. The displacement (about 2.2 kpc) of the binary star system from the centre of the main galaxy is not unusual for stars in the host galaxy or stars originating in the merging galaxy, and therefore any constraints on the kick velocity imparted to the progenitor are poor.
We present constraints on Ks-band emission from one of the nearest short hard gamma-ray bursts, GRB 160821B, at z=0.16, at three epochs. We detect a reddened relativistic afterglow from the jetted emission in the first epoch but do not detect any excess kilonova emission in the second two epochs. We compare upper limits obtained with Keck I/MOSFIRE to multi-dimensional radiative transfer models of kilonovae, that employ composition-dependent nuclear heating and LTE opacities of heavy elements. We discuss eight models that combine toroidal dynamical ejecta and two types of wind and one model with dynamical ejecta only. We also discuss simple, empirical scaling laws of predicted emission as a function of ejecta mass and ejecta velocity. Our limits for GRB 160821B constrain the ejecta mass to be lower than 0.03 Msun for velocities greater than 0.1c. At the distance sensitivity range of advanced LIGO, similar ground-based observations would be sufficiently sensitive to the full range of predicted model emission including models with only dynamical ejecta. The color evolution of these models shows that I-K color spans 7--16 mag, which suggests that even relatively shallow infrared searches for kilonovae could be as constraining as optical searches.
We present mid-infrared (IR) light curves of the Ultraluminous X-ray Source (ULX) Holmberg II X-1 from observations taken between 2014 January 13 and 2017 January 5 with the \textitSpitzer Space Telescope at 3.6 and 4.5 $\mu$m in the \textitSpitzer Infrared Intensive Transients Survey (SPIRITS). The mid-IR light curves, which reveal the first detection of mid-IR variability from a ULX, is determined to arise primarily from dust emission rather than from a jet or an accretion disk outflow. We derived the evolution of the dust temperature ($T_\mathrm{d}\sim600 - 800$ K), IR luminosity ($L_\mathrm{IR}\sim3\times10^4$ $\mathrm{L}_\odot$), mass ($M_\mathrm{d}\sim1-3\times10^{-6}$ $\mathrm{M}_\odot$), and equilibrium temperature radius ($R_\mathrm{eq}\sim10-20$ AU). A comparison of X-1 with a sample spectroscopically identified massive stars in the Large Magellanic Cloud on a mid-IR color-magnitude diagram suggests that the mass donor in X-1 is a supergiant (sg) B[e]-star. The sgB[e]-interpretation is consistent with the derived dust properties and the presence of the [Fe II] ($\lambda=1.644$ $\mu$m) emission line revealed from previous near-IR studies of X-1. We attribute the mid-IR variability of X-1 to increased heating of dust located in a circumbinary torus. It is unclear what physical processes are responsible for the increased dust heating; however, it does not appear to be associated with the X-ray flux from the ULX given the constant X-ray luminosities provided by serendipitous, near-contemporaneous X-ray observations around the first mid-IR variability event in 2014. Our results highlight the importance of mid-IR observations of luminous X-ray sources traditionally studied at X-ray and radio wavelengths.
Massive stellar clusters provide radiation ($\mathrm{\sim 10^7-10^8~L_{\odot}}$) and winds ($\mathrm{\sim 1000~km/s}$) that act to heat dust and shape their surrounding environment. In this paper, the Arched Filaments in the Galactic center were studied to better understand the influence of the Arches cluster on its nearby interstellar medium (ISM). The Arched Filaments were observed with the Faint Object InfraRed CAMera for the SOFIA Telescope (FORCAST) at 19.7, 25.2, 31.5, and 37.1 $\mu$m. Color-temperature maps of the region created with the 25.2 and 37.1 $\mu$m data reveal relatively uniform dust temperatures (70-100 K) over the extent of the filaments ($\sim 25$ pc). Distances between the cluster and the filaments were calculated assuming equilibrium heating of standard size ISM dust grains ($\sim$0.1 $\mu$m). The distances inferred by this method are in conflict with the projected distance between the filaments and the cluster, although this inconsistency can be explained if the characteristic grain size in the filaments is smaller ($\sim$0.01 $\mu$m) than typical values. DustEM models of selected locations within the filaments show evidence of depleted abundances of polycyclic aromatic hydrocarbons (PAHs) by factors of $\sim$1.6-10 by mass compared to the diffuse ISM. The evidence for both PAH depletion and a smaller characteristic grain size points to processing of the ISM within the filaments. We argue that the eroding of dust grains within the filaments is not likely attributable to the radiation or winds from the Arches cluster, but may be related to the physical conditions in the Galactic center.