Showing 1–50 of 84 results for author: Dwarkadas, V V

X-ray Observations of Two Type Ia Supernovae with an H$α$ Line in their Optical Spectrum

Authors: Vikram V. Dwarkadas

Abstract: We report on Chandra X-ray observations of SN 2016jae and SN 2018cqj, both low luminosity Type Ia supernova that showed the presence of a H line in their early optical spectrum. No X-ray emission is detected at the location of either SN. Upper limits to the luminosity of up to 2 $\times 10^{40}\,$ erg s$^{-1}$ are calculated for each SN, depending on the assumed spectral model, temperature and col… ▽ More We report on Chandra X-ray observations of SN 2016jae and SN 2018cqj, both low luminosity Type Ia supernova that showed the presence of a H line in their early optical spectrum. No X-ray emission is detected at the location of either SN. Upper limits to the luminosity of up to 2 $\times 10^{40}\,$ erg s$^{-1}$ are calculated for each SN, depending on the assumed spectral model, temperature and column density. This luminosity is comparable to that of another low-luminosity Type Ia SN, SN 2018fhw, that was observed with Chandra. It is generally lower than upper limits calculated for Type Ia-CSM SNe observed in X-rays, and also below that of SN 2012ca, the only Type Ia-CSM SN to have been detected in X-rays. Comparisons are made to other Type Ia SN with a H line observed in X-rays. The observations suggest that while the density into which the SN is expanding may have been high at the time the H$α$ line was detected, it had decreased considerably by the time of X-ray observations. △ Less

Submitted 28 June, 2024; originally announced July 2024.

Comments: 9 pages, 4 figures and 3 Tables. Accepted to MNRAS. arXiv admin note: text overlap with arXiv:2211.09274

Dark Matter Line Searches with the Cherenkov Telescope Array

Authors: S. Abe, J. Abhir, A. Abhishek, F. Acero, A. Acharyya, R. Adam, A. Aguasca-Cabot, I. Agudo, A. Aguirre-Santaella, J. Alfaro, R. Alfaro, N. Alvarez-Crespo, R. Alves Batista, J. -P. Amans, E. Amato, G. Ambrosi, L. Angel, C. Aramo, C. Arcaro, T. T. H. Arnesen, L. Arrabito, K. Asano, Y. Ascasibar, J. Aschersleben, H. Ashkar , et al. (540 additional authors not shown)

Abstract: Monochromatic gamma-ray signals constitute a potential smoking gun signature for annihilating or decaying dark matter particles that could relatively easily be distinguished from astrophysical or instrumental backgrounds. We provide an updated assessment of the sensitivity of the Cherenkov Telescope Array (CTA) to such signals, based on observations of the Galactic centre region as well as of sele… ▽ More Monochromatic gamma-ray signals constitute a potential smoking gun signature for annihilating or decaying dark matter particles that could relatively easily be distinguished from astrophysical or instrumental backgrounds. We provide an updated assessment of the sensitivity of the Cherenkov Telescope Array (CTA) to such signals, based on observations of the Galactic centre region as well as of selected dwarf spheroidal galaxies. We find that current limits and detection prospects for dark matter masses above 300 GeV will be significantly improved, by up to an order of magnitude in the multi-TeV range. This demonstrates that CTA will set a new standard for gamma-ray astronomy also in this respect, as the world's largest and most sensitive high-energy gamma-ray observatory, in particular due to its exquisite energy resolution at TeV energies and the adopted observational strategy focussing on regions with large dark matter densities. Throughout our analysis, we use up-to-date instrument response functions, and we thoroughly model the effect of instrumental systematic uncertainties in our statistical treatment. We further present results for other potential signatures with sharp spectral features, e.g.~box-shaped spectra, that would likewise very clearly point to a particle dark matter origin. △ Less

Submitted 7 March, 2024; originally announced March 2024.

Comments: 43 pages JCAP style (excluding author list and references), 19 figures

Chasing Gravitational Waves with the Cherenkov Telescope Array

Authors: Jarred Gershon Green, Alessandro Carosi, Lara Nava, Barbara Patricelli, Fabian Schüssler, Monica Seglar-Arroyo, Cta Consortium, :, Kazuki Abe, Shotaro Abe, Atreya Acharyya, Remi Adam, Arnau Aguasca-Cabot, Ivan Agudo, Jorge Alfaro, Nuria Alvarez-Crespo, Rafael Alves Batista, Jean-Philippe Amans, Elena Amato, Filippo Ambrosino, Ekrem Oguzhan Angüner, Lucio Angelo Antonelli, Carla Aramo, Cornelia Arcaro, Luisa Arrabito , et al. (545 additional authors not shown)

Abstract: The detection of gravitational waves from a binary neutron star merger by Advanced LIGO and Advanced Virgo (GW170817), along with the discovery of the electromagnetic counterparts of this gravitational wave event, ushered in a new era of multimessenger astronomy, providing the first direct evidence that BNS mergers are progenitors of short gamma-ray bursts (GRBs). Such events may also produce very… ▽ More The detection of gravitational waves from a binary neutron star merger by Advanced LIGO and Advanced Virgo (GW170817), along with the discovery of the electromagnetic counterparts of this gravitational wave event, ushered in a new era of multimessenger astronomy, providing the first direct evidence that BNS mergers are progenitors of short gamma-ray bursts (GRBs). Such events may also produce very-high-energy (VHE, > 100GeV) photons which have yet to be detected in coincidence with a gravitational wave signal. The Cherenkov Telescope Array (CTA) is a next-generation VHE observatory which aims to be indispensable in this search, with an unparalleled sensitivity and ability to slew anywhere on the sky within a few tens of seconds. New observing modes and follow-up strategies are being developed for CTA to rapidly cover localization areas of gravitational wave events that are typically larger than the CTA field of view. This work will evaluate and provide estimations on the expected number of of gravitational wave events that will be observable with CTA, considering both on- and off-axis emission. In addition, we will present and discuss the prospects of potential follow-up strategies with CTA. △ Less

Submitted 5 February, 2024; v1 submitted 11 October, 2023; originally announced October 2023.

Comments: Presented at the 38th International Cosmic Ray Conference (ICRC 2023), 2023 (arXiv:2309.08219)

Report number: CTA-ICRC/2023/30

Prospects for $γ$-ray observations of the Perseus galaxy cluster with the Cherenkov Telescope Array

Authors: The Cherenkov Telescope Array Consortium, :, K. Abe, S. Abe, F. Acero, A. Acharyya, R. Adam, A. Aguasca-Cabot, I. Agudo, A. Aguirre-Santaella, J. Alfaro, R. Alfaro, N. Alvarez-Crespo, R. Alves Batista, J. -P. Amans, E. Amato, E. O. Angüner, L. A. Antonelli, C. Aramo, M. Araya, C. Arcaro, L. Arrabito, K. Asano, Y. Ascasíbar, J. Aschersleben , et al. (542 additional authors not shown)

Abstract: Galaxy clusters are expected to be dark matter (DM) reservoirs and storage rooms for the cosmic-ray protons (CRp) that accumulate along the cluster's formation history. Accordingly, they are excellent targets to search for signals of DM annihilation and decay at gamma-ray energies and are predicted to be sources of large-scale gamma-ray emission due to hadronic interactions in the intracluster med… ▽ More Galaxy clusters are expected to be dark matter (DM) reservoirs and storage rooms for the cosmic-ray protons (CRp) that accumulate along the cluster's formation history. Accordingly, they are excellent targets to search for signals of DM annihilation and decay at gamma-ray energies and are predicted to be sources of large-scale gamma-ray emission due to hadronic interactions in the intracluster medium. We estimate the sensitivity of the Cherenkov Telescope Array (CTA) to detect diffuse gamma-ray emission from the Perseus galaxy cluster. We perform a detailed spatial and spectral modelling of the expected signal for the DM and the CRp components. For each, we compute the expected CTA sensitivity. The observing strategy of Perseus is also discussed. In the absence of a diffuse signal (non-detection), CTA should constrain the CRp to thermal energy ratio within the radius $R_{500}$ down to about $X_{500}<3\times 10^{-3}$, for a spatial CRp distribution that follows the thermal gas and a CRp spectral index $α_{\rm CRp}=2.3$. Under the optimistic assumption of a pure hadronic origin of the Perseus radio mini-halo and depending on the assumed magnetic field profile, CTA should measure $α_{\rm CRp}$ down to about $Δα_{\rm CRp}\simeq 0.1$ and the CRp spatial distribution with 10% precision. Regarding DM, CTA should improve the current ground-based gamma-ray DM limits from clusters observations on the velocity-averaged annihilation cross-section by a factor of up to $\sim 5$, depending on the modelling of DM halo substructure. In the case of decay of DM particles, CTA will explore a new region of the parameter space, reaching models with $τ_χ>10^{27}$s for DM masses above 1 TeV. These constraints will provide unprecedented sensitivity to the physics of both CRp acceleration and transport at cluster scale and to TeV DM particle models, especially in the decay scenario. △ Less

Submitted 7 September, 2023; originally announced September 2023.

Comments: 93 pages (including author list, appendix and references), 143 figures. Submitted to JCAP

X-Ray Observations of Type Ia Supernova 2018cqj

Authors: Vikram V. Dwarkadas

Abstract: I report on Chandra X-ray observations of SN 2018cqj, a low luminosity Type Ia supernova that showed an H$α$ line in its optical spectrum. No X-ray emission was detected at the location of the SN, with an upper limit to the X-ray luminosity of 2 $\times 10^{39}$ erg s$^{-1}$. I report on Chandra X-ray observations of SN 2018cqj, a low luminosity Type Ia supernova that showed an H$α$ line in its optical spectrum. No X-ray emission was detected at the location of the SN, with an upper limit to the X-ray luminosity of 2 $\times 10^{39}$ erg s$^{-1}$. △ Less

Submitted 21 June, 2023; originally announced June 2023.

Comments: 3 pages, 1 table. Published in RNAAS

Journal ref: RNAAS, 2023, 7, 129

On the Evolution of, and Hot Gas in, Wind-Blown Bubbles around Massive Stars -- Wind Bubbles Are Not Energy-Conserving

Authors: Vikram V. Dwarkadas

Abstract: The structure and evolution of wind-blown bubbles (WBBs) around massive stars has primarily been investigated using an energy-conserving model of wind-blown bubbles. While this model is useful in explaining the general properties of the evolution, several problems remain, including inconsistencies between observed wind luminosities and those derived using this formulation. Major difficulties inclu… ▽ More The structure and evolution of wind-blown bubbles (WBBs) around massive stars has primarily been investigated using an energy-conserving model of wind-blown bubbles. While this model is useful in explaining the general properties of the evolution, several problems remain, including inconsistencies between observed wind luminosities and those derived using this formulation. Major difficulties include the low X-ray temperature and X-ray luminosity, compared to the model. In this paper, we re-examine the evolution, dynamics, and kinematics of WBBs around massive stars, using published ionization gasdynamic simulations of wind-blown bubbles. We show that WBBs can cool efficiently due to the presence of various instabilities and turbulence within the bubble. The expansion of WBBs is more consistent with a momentum-conserving solution, rather than an energy-conserving solution. This compares well with the dynamics and kinematics of observed wind bubbles. Despite the cooling of the bubble, the shocked wind temperature is not reduced to the observed values. We argue that the X-ray emission arise mainly from clumps and filaments within the hot shocked wind region, with temperatures just above 10$^6$ K. The remainder of the plasma can contribute to a lesser extent. △ Less

Submitted 19 June, 2023; originally announced June 2023.

Comments: 22 pages, 7 figures. Published in Galaxies. Only low-resolution figures are included due to size requirements; original high-resolution figures, especially Figure 1, can be obtained from the published paper in the journal

Journal ref: Galaxies 2023, 11(3), 78

Sensitivity of the Cherenkov Telescope Array to spectral signatures of hadronic PeVatrons with application to Galactic Supernova Remnants

Authors: The Cherenkov Telescope Array Consortium, F. Acero, A. Acharyya, R. Adam, A. Aguasca-Cabot, I. Agudo, A. Aguirre-Santaella, J. Alfaro, R. Aloisio, N. Álvarez Crespo, R. Alves Batista, L. Amati, E. Amato, G. Ambrosi, E. O. Angüner, C. Aramo, C. Arcaro, T. Armstrong, K. Asano, Y. Ascasibar, J. Aschersleben, M. Backes, A. Baktash, C. Balazs, M. Balbo , et al. (334 additional authors not shown)

Abstract: The local Cosmic Ray (CR) energy spectrum exhibits a spectral softening at energies around 3~PeV. Sources which are capable of accelerating hadrons to such energies are called hadronic PeVatrons. However, hadronic PeVatrons have not yet been firmly identified within the Galaxy. Several source classes, including Galactic Supernova Remnants (SNRs), have been proposed as PeVatron candidates. The pote… ▽ More The local Cosmic Ray (CR) energy spectrum exhibits a spectral softening at energies around 3~PeV. Sources which are capable of accelerating hadrons to such energies are called hadronic PeVatrons. However, hadronic PeVatrons have not yet been firmly identified within the Galaxy. Several source classes, including Galactic Supernova Remnants (SNRs), have been proposed as PeVatron candidates. The potential to search for hadronic PeVatrons with the Cherenkov Telescope Array (CTA) is assessed. The focus is on the usage of very high energy $γ$-ray spectral signatures for the identification of PeVatrons. Assuming that SNRs can accelerate CRs up to knee energies, the number of Galactic SNRs which can be identified as PeVatrons with CTA is estimated within a model for the evolution of SNRs. Additionally, the potential of a follow-up observation strategy under moonlight conditions for PeVatron searches is investigated. Statistical methods for the identification of PeVatrons are introduced, and realistic Monte--Carlo simulations of the response of the CTA observatory to the emission spectra from hadronic PeVatrons are performed. Based on simulations of a simplified model for the evolution for SNRs, the detection of a $γ$-ray signal from in average 9 Galactic PeVatron SNRs is expected to result from the scan of the Galactic plane with CTA after 10 hours of exposure. CTA is also shown to have excellent potential to confirm these sources as PeVatrons in deep observations with $\mathcal{O}(100)$ hours of exposure per source. △ Less

Submitted 27 March, 2023; originally announced March 2023.

Comments: 34 pages, 16 figures, Accepted for publication in Astroparticle Physics

VERITAS and Fermi-LAT constraints on the Gamma-ray Emission from Superluminous Supernovae SN2015bn and SN2017egm

Authors: A. Acharyya, C. B. Adams, P. Bangale, W. Benbow, J. H. Buckley, M. Capasso, V. V. Dwarkadas, M. Errando, A. Falcone, Q. Feng, J. P. Finley, G. M. Foote, L. Fortson, A. Furniss, G. Gallagher, A. Gent, W. F Hanlon, O. Hervet, J. Holder, T. B. Humensky, W. Jin, P. Kaaret, M. Kertzman, M. Kherlakian, D. Kieda , et al. (34 additional authors not shown)

Abstract: Superluminous supernovae (SLSNe) are a rare class of stellar explosions with luminosities ~10-100 times greater than ordinary core-collapse supernovae. One popular model to explain the enhanced optical output of hydrogen-poor (Type I) SLSNe invokes energy injection from a rapidly spinning magnetar. A prediction in this case is that high-energy gamma rays, generated in the wind nebula of the magnet… ▽ More Superluminous supernovae (SLSNe) are a rare class of stellar explosions with luminosities ~10-100 times greater than ordinary core-collapse supernovae. One popular model to explain the enhanced optical output of hydrogen-poor (Type I) SLSNe invokes energy injection from a rapidly spinning magnetar. A prediction in this case is that high-energy gamma rays, generated in the wind nebula of the magnetar, could escape through the expanding supernova ejecta at late times (months or more after optical peak). This paper presents a search for gamma-ray emission in the broad energy band from 100 MeV to 30 TeV from two Type I SLSNe, SN2015bn, and SN2017egm, using observations from Fermi-LAT and VERITAS. Although no gamma-ray emission was detected from either source, the derived upper limits approach the putative magnetar's spin-down luminosity. Prospects are explored for detecting very-high-energy (VHE; 100 GeV - 100 TeV) emission from SLSNe-I with existing and planned facilities such as VERITAS and CTA. △ Less

Submitted 13 February, 2023; originally announced February 2023.

Comments: 20 pages, 7 figures, 2 tables

Chandra Fails to Detect X-ray Emission from Type Ia SN 2018fhw/ASASSN-18tb

Authors: Vikram V. Dwarkadas

Abstract: We report on Chandra X-ray observations of ASASSN-18tb/SN 2018fhw, a low luminosity Type Ia supernova that showed a H line in its optical spectrum. No X-ray emission was detected at the location of the SN. Upper limits to the luminosity of up to 3 $\times 10^{39}$ erg s$^{-1}$ are calculated, depending on the assumed spectral model, temperature and column density. These are compared to two Type Ia… ▽ More We report on Chandra X-ray observations of ASASSN-18tb/SN 2018fhw, a low luminosity Type Ia supernova that showed a H line in its optical spectrum. No X-ray emission was detected at the location of the SN. Upper limits to the luminosity of up to 3 $\times 10^{39}$ erg s$^{-1}$ are calculated, depending on the assumed spectral model, temperature and column density. These are compared to two Type Ia-CSM SNe, SN 2005gj and SN 2002ic, that have been observed with Chandra in the past. The upper limits are lower than the X-ray luminosity found for the Type Ia-CSM SN 2012ca, the only Type Ia SN to have been detected in X-rays. Consideration of various scenarios for the H$α$ line suggests that the density of the surrounding medium at the time of H$α$ line detection could have been as high as 10$^8$ cm$^{-3}$, but must have decreased below 5 $\times\, 10^6$ cm$^{-3}$ at the time of X-ray observation. Continual X-ray observations of SNe which show a H line in their spectrum are necessary in order to establish Type Ia SNe as an X-ray emitting class. △ Less

Submitted 16 November, 2022; originally announced November 2022.

Comments: 9 pages, 1 figure, 3 tables. Accepted to MNRAS

The Bright Supernova 1996cr in the Circinus Galaxy Imaged with VLBI: Shell Structure with Complex Evolution

Authors: Michael F. Bietenholz, Norbert Bartel, Vikram V. Dwarkadas, Leon Mtshweni, Carlos Orquera-Rojas, Simon Ellingsen, Shinji Horiuchi, Anastasios Tzioumis

Abstract: We present broadband radio flux-density measurements supernova (SN) 1996cr, made with MeerKAT, ATCA and ALMA, and images made from very long baseline interferometry (VLBI) observations with the Australian Long Baseline Array. The spectral energy distribution of SN 1996cr in 2020, at age, $t \sim$8700 d, is a power-law, with flux density, $S \propto ν^{-0.588 \pm 0.011}$ between 1 and 34 GHz, but m… ▽ More We present broadband radio flux-density measurements supernova (SN) 1996cr, made with MeerKAT, ATCA and ALMA, and images made from very long baseline interferometry (VLBI) observations with the Australian Long Baseline Array. The spectral energy distribution of SN 1996cr in 2020, at age, $t \sim$8700 d, is a power-law, with flux density, $S \propto ν^{-0.588 \pm 0.011}$ between 1 and 34 GHz, but may steepen at $>35$ GHz. The spectrum has flattened since $t = 5370$ d (2010). Also since $t = 5370$ d, the flux density has declined rapidly, with $S_{\rm 9 \, GHz} \propto t^{-2.9}$. The VLBI image at $t = 8859$ d shows an approximately circular structure, with a central minimum reminiscent of an optically-thin spherical shell of emission. For a distance of 3.7 Mpc, the average outer radius of the radio emission at $t = 8859$ d was $(5.1 \pm 0.3) \times 10^{17}$ cm, and SN 1996cr has been expanding with a velocity of $4650 \pm 1060$ km s$^{-1}$ between $t=4307$ and 8859 d. It must have undergone considerable deceleration before $t = 4307$ d. Deviations from a circular shell structure in the image suggest a range of velocities up to $\sim$7000 km s$^{-1}$, and hint at the presence of a ring- or equatorial-belt-like structure rather than a complete spherical shell. △ Less

Submitted 27 February, 2023; v1 submitted 30 October, 2022; originally announced October 2022.

Comments: 10 Pages, 5 figures; Accepted for publication in MNRAS; this version (3) differs from version 2 by further minor revision

An Exploration of X-ray Supernova Remnants in the Milky Way and Nearby Galaxies

Authors: Chris Albert, Vikram V Dwarkadas

Abstract: We probe the environmental properties of X-ray supernova remnants (SNRs) at various points along their evolutionary journey, especially the S-T phase, and their conformance with theoretically derived models of SNR evolution. The remnant size is used as a proxy for the age of the remnant. Our data set includes 34 Milky Way, 59 Large Magellanic Cloud (LMC), and 5 Small Magellanic Cloud (SMC) SNRs. W… ▽ More We probe the environmental properties of X-ray supernova remnants (SNRs) at various points along their evolutionary journey, especially the S-T phase, and their conformance with theoretically derived models of SNR evolution. The remnant size is used as a proxy for the age of the remnant. Our data set includes 34 Milky Way, 59 Large Magellanic Cloud (LMC), and 5 Small Magellanic Cloud (SMC) SNRs. We select remnants that have been definitively typed as either core-collapse (CC) or Type Ia supernovae, with well-defined size estimates, and a thermal X-ray flux measured over the entire remnant. A catalog of SNR size and X-ray luminosity is presented and plotted, with ambient density and age estimates from the literature. Model remnants with a given density, in the Sedov-Taylor (S-T) phase, are overplotted on the diameter-vs-luminosity plot, allowing the evolutionary state and physical properties of SNRs to be compared to each other, and to theoretical models. We find that small, young remnants are predominantly Type Ia remnants or high luminosity CCs, suggesting that many CC SNRs are not detected until after they have emerged from the progenitor's wind-blown bubble. An examination of the distribution of SNR diameters in the Milky Way and LMC reveals that LMC SNRs must be evolving in an ambient medium which is 30% as dense as that in the Milky Way. This is consistent with ambient density estimates for the Galaxy and LMC. △ Less

Submitted 10 May, 2022; originally announced May 2022.

Comments: 22 pages, 6 figures. Accepted to MNRAS. Online version of Fig 2 with all SNR names is available at https://voices.uchicago.edu/vikramvd/snrspaper/

Seven Years of SN 2014C: a Multi-Wavelength Synthesis of an Extraordinary Supernova

Authors: Benjamin P. Thomas, J. Craig Wheeler, Vikram V. Dwarkadas, Christopher Stockdale, Jozsef Vinko, David Pooley, Yerong Xu, Greg Zeimann, Phillip MacQueen

Abstract: SN 2014C was originally classified as a Type Ib supernova, but at phase φ = 127 d post-explosion strong Hα emission was observed. SN 2014C has since been observed in radio, infrared, optical and X-ray bands. Here we present new optical spectroscopic and photometric data spanning φ = 947 - 2494 d post-explosion. We address the evolution of the broadened Hα emission line, as well as broad [O III] em… ▽ More SN 2014C was originally classified as a Type Ib supernova, but at phase φ = 127 d post-explosion strong Hα emission was observed. SN 2014C has since been observed in radio, infrared, optical and X-ray bands. Here we present new optical spectroscopic and photometric data spanning φ = 947 - 2494 d post-explosion. We address the evolution of the broadened Hα emission line, as well as broad [O III] emission and other lines. We also conduct a parallel analysis of all publicly available multi-wavelength data. From our spectra, we find a nearly constant Hα FWHM velocity width of {\sim}2000 km/s that is significantly lower than that of other broadened atomic transitions ({\sim}3000 - 7000 km/s) present in our spectra ([O I] λ6300; [O III] λλ4959,5007; He I λ7065; [Ca II] λλ7291,7324). The late radio data demand a fast forward shock ({\sim}10,000 km/s at φ = 1700 d) in rarified matter that contrasts with the modest velocity of the Hα. We propose that the infrared flux originates from a toroidal-like structure of hydrogen surrounding the progenitor system, while later emission at other wavelengths (radio, X-ray) likely originates predominantly from the reverse shock in the ejecta and the forward shock in the quasi-spherical progenitor He wind. We propose that the Hα emission arises in the boundary layer between the ejecta and torus. We also consider the possible roles of a pulsar and a binary companion. △ Less

Submitted 23 March, 2022; originally announced March 2022.

Comments: 33 pages, 12 figures

Ionization-Gasdynamic Simulations of Wind-Blown Nebulae around Massive Stars

Authors: Vikram V. Dwarkadas

Abstract: Using a code that employs a self-consistent method for computing the effects of photo-ionization on circumstellar gas dynamics, we model the formation of wind-driven nebulae around massive stars. We take into account changes in stellar properties and mass-loss over the star's evolution. Our simulations show how various properties, such as the density and ionization fraction, change throughout the… ▽ More Using a code that employs a self-consistent method for computing the effects of photo-ionization on circumstellar gas dynamics, we model the formation of wind-driven nebulae around massive stars. We take into account changes in stellar properties and mass-loss over the star's evolution. Our simulations show how various properties, such as the density and ionization fraction, change throughout the evolution of the star. The multi-dimensional simulations reveal the presence of strong ionization front instabilities in the main-sequence phase, similar to those seen in galactic ionization fronts. Hydrodynamic instabilities at the interfaces lead to the formation of filaments and clumps that are continually being stripped off and mixed with the low density interior. Even though the winds start out as completely radial, the spherical symmetry is quickly destroyed, and the shocked wind region is manifestly asymmetrical. The simulations demonstrate that it is important to include the effects of the photoionizing photons from the star, and simulations that do not include this may fail to reproduce the observed density profile and ionization structure of wind-blown bubbles around massive stars. △ Less

Submitted 18 February, 2022; originally announced February 2022.

Comments: 3 figures, 1 table. 13 pages. Author version of published paper, with reduced size figures 2 and 3. For full resolution figures see the published paper at https://www.mdpi.com/2075-4434/10/1/37. Accepted to Galaxies. This article belongs to the special issue of the conference Asymmetrical Post Main-Sequence Nebulae 8

Journal ref: Galaxies 2022, 10, 37

The first days of type II-P core collapse supernovae in the gamma-ray range

Authors: P. Cristofari, A. Marcowith, M. Renaud, V. V. Dwarkadas, V. Tatischeff, G. Giacinti, E. Peretti, H. Sol

Abstract: Type II-P supernovæ~(SNe), the most common core-collapse SNe type, result from the explosions of red supergiant stars. Their detection in the radio domain testifies of the presence of relativistic electrons, and shows that they are potentially efficient energetic particle accelerators. If hadrons can also be accelerated, these energetic particles are expected to interact with the surrounding mediu… ▽ More Type II-P supernovæ~(SNe), the most common core-collapse SNe type, result from the explosions of red supergiant stars. Their detection in the radio domain testifies of the presence of relativistic electrons, and shows that they are potentially efficient energetic particle accelerators. If hadrons can also be accelerated, these energetic particles are expected to interact with the surrounding medium to produce a gamma-ray signal even in the multi--TeV range. The intensity of this signal depends on various factors, but an essential one is the density of the circumstellar medium. Such a signal should however be limited by electron-positron pair production arising from the interaction of the gamma-ray photons with optical photons emitted by the supernova photosphere, which can potentially degrade the gamma-ray signal by over ten orders of magnitude in the first days/weeks following the explosion. We calculate the gamma-gamma opacity from a detailed modelling of the time evolution of the forward shock and supernova photosphere, taking a full account of the non-isotropy of the photon interactions. We discuss the time-dependent gamma-ray TeV emission from type II-P SNe as a function of the stellar progenitor radius and mass-loss rate, as well as the explosion energy and mass of the ejected material. We evaluate the detectability of the SNe with the next generation of Cherenkov telescopes. We find that, while most extragalactic events may be undetectable, type II-P SNe exploding in our Galaxy or in the Magellanic Clouds should be detected by gamma-ray observatories such as the upcoming Cherenkov Telescope Array. △ Less

Submitted 24 January, 2022; originally announced January 2022.

Comments: 9 pages, 8 figures, Accepted for publication in MNRAS

A Long-term Study of Ultraluminous X-ray Sources in NGC 891

Authors: Nicholas M. Earley, Vikram V. Dwarkadas, Victoria Cirillo

Abstract: We perform empirical fits to the \emph{Chandra} and \emph{XMM-Newton} spectra of three ultraluminous X-ray sources (ULXs) in the edge-on spiral galaxy NGC 891, monitoring the region over a seventeen year time window. One of these sources has been visible since the early 1990s with \emph{ROSAT} and has been observed multiple times with \emph{Chandra} and \emph{XMM-Newton}. Another has been visible… ▽ More We perform empirical fits to the \emph{Chandra} and \emph{XMM-Newton} spectra of three ultraluminous X-ray sources (ULXs) in the edge-on spiral galaxy NGC 891, monitoring the region over a seventeen year time window. One of these sources has been visible since the early 1990s with \emph{ROSAT} and has been observed multiple times with \emph{Chandra} and \emph{XMM-Newton}. Another has been visible since 2011. We build upon prior analyses of these sources by analyzing all available data at all epochs. Where possible \emph{Chandra} data is used, since its superior spatial resolution allows for more effective isolation of the emission from each individual source, thus providing a better determination of their spectral properties. We also identify a new transient ULX, CXOU J022230.1+421937, which faded from view over the course of a two month period from Nov 2016 to Jan 2017. Modeling of each source at every epoch was conducted using six different models ranging from thermal bremsstrahlung to accretion disk models. Unfortunately, but as is common with many ULXs, no single model yielded a much better fit than the others. The two known sources had unabsorbed luminosities that remained fairly consistent over five or more years. Various possibilities for the new transient ULX are explored. △ Less

Submitted 22 December, 2021; originally announced December 2021.

Comments: 23 pages, 4 tables, 5 figures. Accepted to Universe journal

Journal ref: Universe 2022, 8, 18

Can the Fe K-alpha Line Reliably Predict Supernova Remnant Progenitors?

Authors: Jared Siegel, Vikram V. Dwarkadas, Kari A. Frank, David N. Burrows

Abstract: The centroid energy of the Fe K$α$ line has been used to identify the progenitors of supernova remnants (SNRs). These investigations generally considered the energy of the centroid derived from the spectrum of the entire remnant. Here we use {\it XMM-Newton} data to investigate the Fe K$α$ centroid in 6 SNRs: 3C~397, N132D, W49B, DEM L71, 1E 0102.2-7219, and Kes 73. In Kes 73 and 1E 0102.2-7219, w… ▽ More The centroid energy of the Fe K$α$ line has been used to identify the progenitors of supernova remnants (SNRs). These investigations generally considered the energy of the centroid derived from the spectrum of the entire remnant. Here we use {\it XMM-Newton} data to investigate the Fe K$α$ centroid in 6 SNRs: 3C~397, N132D, W49B, DEM L71, 1E 0102.2-7219, and Kes 73. In Kes 73 and 1E 0102.2-7219, we fail to detect any Fe K$α$ emission. We report a tentative first detection of Fe K$α$ emission in SNR DEM L71, with a centroid energy consistent with its Type Ia designation. In the remaining remnants, the spatial and spectral sensitivity is sufficient to investigate spatial variations of the Fe K$α$ centroid. We find in N132D and W49B that the centroids in different regions are consistent with that derived from the overall spectrum, although not necessarily with the remnant type identified via other means. However, in SNR 3C~397, we find statistically significant variation in the centroid of up to 100 eV, aligning with the variation in the density structure around the remnant. These variations span the intermediate space between centroid energies signifying core-collapse and Type Ia remnants. Shifting the dividing line downwards by 50 eV can place all the centroids in the CC region, but contradicts the remnant type obtained via other means. Our results show that caution must be used when employing the Fe K$α$ centroid of the entire remnant as the sole diagnostic for typing a remnant. △ Less

Submitted 2 September, 2021; originally announced September 2021.

Comments: 10 pages, 3 figures. Accepted to the Astrophysical Journal

X-ray Luminosity of Decades-Old Supernovae

Authors: Vandana Ramakrishnan, Vikram V. Dwarkadas

Abstract: The transition from supernovae (SNe) to supernova remnants (SNRs) remains poorly understood, given the age gap between well-studied examples of the two. In order to bridge this gap, we analysed archival Chandra data for some of the oldest supernovae detected in X-rays, in order to extend their light curves out to late times. We fitted the spectra with thermal models. All the SNe with multiple X-ra… ▽ More The transition from supernovae (SNe) to supernova remnants (SNRs) remains poorly understood, given the age gap between well-studied examples of the two. In order to bridge this gap, we analysed archival Chandra data for some of the oldest supernovae detected in X-rays, in order to extend their light curves out to late times. We fitted the spectra with thermal models. All the SNe with multiple X-ray data points were found to have similar X-ray luminosity, which was decreasing with time. The X-ray luminosity will likely continue to decrease while the SNe are evolving in a wind medium, but is anticipated to increase in the Sedov phase when the SNe are interacting with a constant density interstellar medium, bringing it in line with observed SNRs. △ Less

Submitted 21 August, 2021; originally announced August 2021.

Comments: Published in RNAAS.3 pages, 1 figure. Part of Focus on AAS 238 issue

Journal ref: Research Notes of the AAS, Volume 5, Number 8, 191 (2021)

Two Decades of Education and Public Outreach with Chicago Public Schools

Authors: Vikram V. Dwarkadas

Abstract: Over the past two decades, I have been actively involved in teaching astronomy and astrophysics to Chicago Public School (CPS) students and their teachers, in collaboration with various groups as well as by myself. Valuable resources that we have created for schools include the Multiwavelength Astronomy Website, with modules for infrared, optical, ultraviolet, X-ray and gamma-ray astronomy. The co… ▽ More Over the past two decades, I have been actively involved in teaching astronomy and astrophysics to Chicago Public School (CPS) students and their teachers, in collaboration with various groups as well as by myself. Valuable resources that we have created for schools include the Multiwavelength Astronomy Website, with modules for infrared, optical, ultraviolet, X-ray and gamma-ray astronomy. The content of each lesson is derived from interviews with scientists, archived oral histories, and/or memoirs. Lessons are evaluated by a science educator and at least one subject matter expert before being produced for the web. They are supplemented by NASA media, archival material from the University of Chicago Library and other archives, and participant contributed photographs, light curves, and spectra. Summer programs provided training to CPS teachers to use the resources in their classrooms. Currently, I lead the Chicago Area Research Mentoring (CHARM) initiative. In the past academic year I worked with a class of 17 diverse 11th grade honors students at the University of Chicago Charter School, Woodlawn. Through frequent lectures ($\sim$ every 4 weeks), these students were exposed to astrophysical topics and concepts not normally covered in a school curriculum. CHARM aims to develop the student's critical thinking, introduce them to astrophysical research methods and techniques, and prepare them for a career in science, technology, engineering and mathematics (STEM), particularly a research-oriented one. In this article, I highlight some projects, educational resources, results achieved, and lessons learned along the way. △ Less

Submitted 17 August, 2021; originally announced August 2021.

Comments: 8 pages, 4 figures. To appear in the proceedings of the International Cosmic Ray Conference 2021. To be published in the Proceedings of Science. All comments appreciated

Supernova X-Ray Database (SNaX) Updated to Ensure Long-term Stability

Authors: Alexandra Nisenoff, Vikram V. Dwarkadas, Mathias C. Ross

Abstract: The Supernova X-Ray Database (SNaX) was established a few years ago to make X-ray data on supernovae (SNe) publicly available via an elegant searchable web interface. The database has recently been updated to PhP7, had security updates done, and moved to a new server, ensuring its long-term stability. We urge astronomers to continue to download the data as needed for their work. Those with X-ray d… ▽ More The Supernova X-Ray Database (SNaX) was established a few years ago to make X-ray data on supernovae (SNe) publicly available via an elegant searchable web interface. The database has recently been updated to PhP7, had security updates done, and moved to a new server, ensuring its long-term stability. We urge astronomers to continue to download the data as needed for their work. Those with X-ray data on SNe are requested to upload it to the database via the easily fillable spreadsheet, making it accessible to everyone. △ Less

Submitted 3 November, 2020; originally announced November 2020.

Comments: 3 pages, 1 figure. Published in the Research Notes of the American Astronomical Society. Just a small reminder to astronomers to use the database for X-ray data on SNe, and to upload your data on X-ray SNe to the database

Journal ref: Res. Notes AAS, 4, 195 (2020)

Analysis of XMM-Newton Observations of Supernova Remnant W49B and Clues to the Progenitor

Authors: Jared Siegel, Vikram V. Dwarkadas, Kari A. Frank, David N. Burrows

Abstract: W49B is a supernova remnant (SNR) discovered over 60 years ago in early radio surveys. It has since been observed over the entire wavelength range, with the X-ray morphology resembling a centrally-filled SNR. The nature of its progenitor star is still debated. Applying Smoothed Particle Inference techniques to analyze the X-Ray emission from W49B, we characterize the morphology and abundance distr… ▽ More W49B is a supernova remnant (SNR) discovered over 60 years ago in early radio surveys. It has since been observed over the entire wavelength range, with the X-ray morphology resembling a centrally-filled SNR. The nature of its progenitor star is still debated. Applying Smoothed Particle Inference techniques to analyze the X-Ray emission from W49B, we characterize the morphology and abundance distribution over the entire remnant. We also infer the density structure and derive the mass of individual elements present in the plasma. The morphology is consistent with an interaction between the remnant and a dense medium along the eastern edge, and some obstruction towards the west. We find a total mass of 130 $(\pm 16)$ M$_{\odot}$ and an estimated ejecta mass of 1.2 $(\pm 0.2)$ M$_{\odot}$. Comparison of the inferred abundance values and individual element masses with a wide selection of SN models suggests that deflagration-to-detonation (DDT) Type Ia models are the most compatible, with Fe abundance being the major discriminating factor. The general agreement between our abundance measurements and those from previous studies suggests that disagreement between various authors is more likely due to the choice of models used for comparison, rather than the abundance values themselves. While our abundance results lean toward a Type Ia origin, ambiguities in the interpretation of various morphological and spectral characteristics of W49B do not allow us to provide a definitive classification. △ Less

Submitted 9 October, 2020; originally announced October 2020.

Comments: 14 pages, 8 figures, 5 tables. Accepted to the Astrophysical Journal

Sensitivity of the Cherenkov Telescope Array for probing cosmology and fundamental physics with gamma-ray propagation

Authors: The Cherenkov Telescope Array Consortium, :, H. Abdalla, H. Abe, F. Acero, A. Acharyya, R. Adam, I. Agudo, A. Aguirre-Santaella, R. Alfaro, J. Alfaro, C. Alispach, R. Aloisio, R. Alves B, L. Amati, E. Amato, G. Ambrosi, E. O. Angüner, A. Araudo, T. Armstrong, F. Arqueros, L. Arrabito, K. Asano, Y. Ascasíbar, M. Ashley , et al. (474 additional authors not shown)

Abstract: The Cherenkov Telescope Array (CTA), the new-generation ground-based observatory for $γ$-ray astronomy, provides unique capabilities to address significant open questions in astrophysics, cosmology, and fundamental physics. We study some of the salient areas of $γ$-ray cosmology that can be explored as part of the Key Science Projects of CTA, through simulated observations of active galactic nucle… ▽ More The Cherenkov Telescope Array (CTA), the new-generation ground-based observatory for $γ$-ray astronomy, provides unique capabilities to address significant open questions in astrophysics, cosmology, and fundamental physics. We study some of the salient areas of $γ$-ray cosmology that can be explored as part of the Key Science Projects of CTA, through simulated observations of active galactic nuclei (AGN) and of their relativistic jets. Observations of AGN with CTA will enable a measurement of $γ$-ray absorption on the extragalactic background light with a statistical uncertainty below 15% up to a redshift $z=2$ and to constrain or detect $γ$-ray halos up to intergalactic-magnetic-field strengths of at least 0.3pG. Extragalactic observations with CTA also show promising potential to probe physics beyond the Standard Model. The best limits on Lorentz invariance violation from $γ$-ray astronomy will be improved by a factor of at least two to three. CTA will also probe the parameter space in which axion-like particles could constitute a significant fraction, if not all, of dark matter. We conclude on the synergies between CTA and other upcoming facilities that will foster the growth of $γ$-ray cosmology. △ Less

Submitted 26 February, 2021; v1 submitted 3 October, 2020; originally announced October 2020.

Comments: 71 pages (including affiliations and references), 13 figures, 6 tables. Accepted in JCAP; matches published version. Corresponding authors: Jonathan Biteau, Julien Lefaucheur, Humberto Martinez-Huerta, Manuel Meyer, Santiago Pita, Ievgen Vovk

Journal ref: JCAP 02 (2021) 048

From Supernova to Remnant: Tracking the Evolution of the Oldest Known X-ray Supernovae

Authors: Vandana Ramakrishnan, Vikram V. Dwarkadas

Abstract: Core-collapse supernovae (SNe) expand into a medium created by winds from the pre-SN progenitor. The SN explosion and resulting shock wave(s) heat up the surrounding plasma, giving rise to thermal X-ray emission, which depends on the density of the emitting material. Tracking the variation of the X-ray luminosity over long periods of time thus allows for investigation of the kinematics of the SN s… ▽ More Core-collapse supernovae (SNe) expand into a medium created by winds from the pre-SN progenitor. The SN explosion and resulting shock wave(s) heat up the surrounding plasma, giving rise to thermal X-ray emission, which depends on the density of the emitting material. Tracking the variation of the X-ray luminosity over long periods of time thus allows for investigation of the kinematics of the SN shock waves, the structure of the surrounding medium, and the nature of the progenitor star. In this paper X-ray observations of five of the oldest known X-ray supernovae - SN 1970G, SN 1968D, SN 1959D, SN 1957D and SN 1941C - are analyzed, with the aim of reconstructing their light curves over several decades. For those supernovae for which we can extract multi-epoch data, the X-ray luminosity appears to decline with time, although with large error bars. No increase in the X-ray emission from SN 1970G is found at later epochs, contrary to previous reports. All five SNe show X-ray luminosities that are of comparable magnitude. We compare the late-time X-ray luminosities of these SNe to those of supernova remnants (SNRs) in the Galaxy which are a few hundred years old, and find that when the tentative decline is taken into account, the luminosity of the old SNe studied herein could fall below the luminosity of some of the younger SNRs within a few hundred years. However, the X-ray luminosity should begin to increase as the SNe expand in the Sedov phase, thus reaching that of the observed SNRs. △ Less

Submitted 20 August, 2020; originally announced August 2020.

Comments: Accepted to the Astrophysical Journal. 24 pages, 12 tables, 16 figures

Stochastic re-acceleration and magnetic-field damping in Tycho's supernova remnant

Authors: A. Wilhelm, I. Telezhinsky, V. V. Dwarkadas, M. Pohl

Abstract: A number of studies suggest that shock acceleration with particle feedback and very efficient magnetic-field amplification combined with Alfvénic drift are needed to explain the rather soft radio spectrum and the narrow rims observed for Tycho's SNR. We show that the broadband spectrum of Tycho's SNR can alternatively be well explained when accounting for stochastic acceleration as a secondary pro… ▽ More A number of studies suggest that shock acceleration with particle feedback and very efficient magnetic-field amplification combined with Alfvénic drift are needed to explain the rather soft radio spectrum and the narrow rims observed for Tycho's SNR. We show that the broadband spectrum of Tycho's SNR can alternatively be well explained when accounting for stochastic acceleration as a secondary process. The re-acceleration of particles in the turbulent region immediately downstream of the shock should be efficient enough to impact particle spectra over several decades in energy. The so-called Alfvénic drift and particle feedback on the shock structure are not required in this scenario. Additionally, we investigate whether synchrotron losses or magnetic-field damping play a more profound role in the formation of the non-thermal filaments. We solve the full particle transport equation in test-particle mode using hydrodynamic simulations of the SNR plasma flow. The background magnetic field is either computed from the induction equation or follows analytic profiles, depending on the model considered. Fast-mode waves in the downstream region provide the diffusion of particles in momentum space. We show that the broadband spectrum of Tycho can be well explained if magnetic-field damping and stochastic re-acceleration of particles are taken into account. Although not as efficient as standard DSA, stochastic acceleration leaves its imprint on the particle spectra, which is especially notable in the emission at radio wavelengths. We find a lower limit for the post-shock magnetic-field strength $\sim330\,\mathrm{μG}$, implying efficient amplification even for the magnetic-field damping scenario. For the formation of the filaments in the radio range magnetic-field damping is necessary, while the X-ray filaments are shaped by both the synchrotron losses and magnetic-field damping. △ Less

Submitted 8 June, 2020; originally announced June 2020.

Comments: 14 pages, 10 figures, accepted for publication in Astronomy and Astrophysics

Journal ref: A&A 639, A124 (2020)

Time-dependent high-energy gamma-ray signal from accelerated particles in core-collapse supernovae: the case of SN 1993J

Authors: Pierre Cristofari, Matthieu Renaud, Alexandre Marcowith, Vikram V. Dwarkadas, Vincent Tatischeff

Abstract: Some core-collapse supernovae are likely to be efficient cosmic-ray accelerators up to the PeV range, and therefore, to potentially play an important role in the overall Galactic cosmic-ray population. The TeV gamma-ray domain can be used to study particle acceleration in the multi-TeV and PeV range. This motivates the study of the detectability of such supernovae by current and future gamma-ray f… ▽ More Some core-collapse supernovae are likely to be efficient cosmic-ray accelerators up to the PeV range, and therefore, to potentially play an important role in the overall Galactic cosmic-ray population. The TeV gamma-ray domain can be used to study particle acceleration in the multi-TeV and PeV range. This motivates the study of the detectability of such supernovae by current and future gamma-ray facilities. The gamma-ray emission of core-collapse supernovae strongly depends on the level of the two-photon annihilation process: high-energy gamma-ray photons emitted at the expanding shock wave following the supernova explosion can interact with soft photons from the supernova photosphere through the pair production channel, thereby strongly suppressing the flux of gamma rays leaving the system. In the case of SN 1993J, whose photospheric and shock-related parameters are well measured, we calculate the temporal evolution of the expected gamma-ray attenuation by accounting for the temporal and geometrical effects. We find the attenuation to be of about $10$ orders of magnitude in the first few days after the SN explosion. The probability of detection of a supernova similar to SN 1993J with the Cherenkov Telescope Array is highest if observations are performed either earlier than 1 day, or later than 10 days after the explosion, when the gamma-ray attenuation decreases to about $2$ orders of magnitude. △ Less

Submitted 6 April, 2020; originally announced April 2020.

Comments: Accepted to MNRAS, 7 pages, 7 figures

Evidence for proton acceleration up to TeV energies based on VERITAS and Fermi-LAT observations of the Cas A SNR

Authors: A. U. Abeysekara, A. Archer, W. Benbow, R. Bird, R. Brose, M. Buchovecky, J. H. Buckley, A. J. Chromey, W. Cui, M. K. Daniel, S. Das, V. V. Dwarkadas, A. Falcone, Q. Feng, J. P. Finley, L. Fortson, A. Gent, G. H. Gillanders, C. Giuri, O. Gueta, D. Hanna, T. Hassan, O. Hervet, J. Holder, G. Hughes , et al. (38 additional authors not shown)

Abstract: We present a study of $γ$-ray emission from the core-collapse supernova remnant Cas~A in the energy range from 0.1GeV to 10TeV. We used 65 hours of VERITAS data to cover 200 GeV - 10 TeV, and 10.8 years of \textit{Fermi}-LAT data to cover 0.1-500 GeV. The spectral analysis of \textit{Fermi}-LAT data shows a significant spectral curvature around $1.3 \pm 0.4_{stat}$ GeV that is consistent with the… ▽ More We present a study of $γ$-ray emission from the core-collapse supernova remnant Cas~A in the energy range from 0.1GeV to 10TeV. We used 65 hours of VERITAS data to cover 200 GeV - 10 TeV, and 10.8 years of \textit{Fermi}-LAT data to cover 0.1-500 GeV. The spectral analysis of \textit{Fermi}-LAT data shows a significant spectral curvature around $1.3 \pm 0.4_{stat}$ GeV that is consistent with the expected spectrum from pion decay. Above this energy, the joint spectrum from \textit{Fermi}-LAT and VERITAS deviates significantly from a simple power-law, and is best described by a power-law with spectral index of $2.17\pm 0.02_{stat}$ with a cut-off energy of $2.3 \pm 0.5_{stat}$ TeV. These results, along with radio, X-ray and $γ$-ray data, are interpreted in the context of leptonic and hadronic models. Assuming a one-zone model, we exclude a purely leptonic scenario and conclude that proton acceleration up to at least 6 TeV is required to explain the observed $γ$-ray spectrum. From modeling of the entire multi-wavelength spectrum, a minimum magnetic field inside the remnant of $B_{\mathrm{min}}\approx150\,\mathrm{μG}$ is deduced. △ Less

Submitted 30 March, 2020; originally announced March 2020.

Comments: 33 pages, 9 Figures, 6 Tables

SPI Analysis and Abundance Calculations of DEM L71, and Comparison to SN explosion Models

Authors: Jared Siegel, Vikram V. Dwarkadas, Kari Frank, David N. Burrows, Aldo Panfichi

Abstract: We analyze the X-Ray emission from the supernova remnant DEM L71 using the Smoothed Particle Inference (SPI) technique. The high Fe abundance found appears to confirm the Type Ia origin. Our method allows us to separate the material ejected in the supernova explosion from the material swept-up by the supernova shock wave. We are able to calculate the total mass of this swept-up material to be abou… ▽ More We analyze the X-Ray emission from the supernova remnant DEM L71 using the Smoothed Particle Inference (SPI) technique. The high Fe abundance found appears to confirm the Type Ia origin. Our method allows us to separate the material ejected in the supernova explosion from the material swept-up by the supernova shock wave. We are able to calculate the total mass of this swept-up material to be about 228 $\pm$ 23 M$_{\odot}$. We plot the posterior distribution for the number density parameter, and create a map of the density structure within the remnant. While the observed density shows substantial variations, we find our results are generally consistent with a two-dimensional hydrodynamical model of the remnant that we have run. Assuming the ejected material arises from a Type Ia explosion, with no hydrogen present, we use the predicted yields from Type Ia models available in the literature to characterize the emitting gas. We find that the abundance of various elements match those predicted by deflagration to detonation transition (DDT) models. Our results, compatible with the Type Ia scenario, highlight the complexity of the remnant and the nature of the surrounding medium. △ Less

Submitted 27 January, 2020; originally announced January 2020.

Comments: 6 pages, 7 Figures. Accepted for publication in Astronomische Nachrichten. To appear as part of the Proceedings of the XMM-Newton 2019 Workshop on "Astrophysics of Hot Plasma in Extended X-Ray Sources"

Interaction of SN~Ib 2004dk with a Previously-Expelled Envelope

Authors: David Pooley, J. Craig Wheeler, Jozsef Vinkó, Vikram V. Dwarkadas, Tamas Szalai, Jeffrey M. Silverman, Madelaine Griesel, Molly McCullough, G. H. Marion, Phillip MacQueen

Abstract: The interaction between the expanding supernova (SN) ejecta with the circumstellar material (CSM) that was expelled from the progenitor prior to explosion is a long-sought phenomenon, yet observational evidence is scarce. Here we confirm a new example: SN 2004dk, originally a hydrogen-poor, helium-rich Type Ib SN that reappeared as a strong H$α$-emitting point-source on narrowband H$α$ images. We… ▽ More The interaction between the expanding supernova (SN) ejecta with the circumstellar material (CSM) that was expelled from the progenitor prior to explosion is a long-sought phenomenon, yet observational evidence is scarce. Here we confirm a new example: SN 2004dk, originally a hydrogen-poor, helium-rich Type Ib SN that reappeared as a strong H$α$-emitting point-source on narrowband H$α$ images. We present follow-up optical spectroscopy that reveals the presence of a broad H$α$ component with full width at half maximum of ~290 km/s in addition to the narrow H$α$ +[NII] emission features from the host galaxy. Such a broad component is a clear sign of an ejecta-CSM interaction. We also present observations with the XMM-Newton Observatory, the Swift satellite, and the Chandra X-ray Observatory that span 10 days to 15 years after discovery. The detection of strong radio, X-ray, and H$α$ emission years after explosion allows various constraints to be put on pre-SN mass-loss processes. We present a wind-bubble model in which the CSM is "pre-prepared" by a fast wind interacting with a slow wind. Much of the outer density profile into which the SN explodes corresponds to no steady-state mass-loss process. We estimate that the shell of compressed slow wind material was ejected ~1400 yr prior to explosion, perhaps during carbon burning, and that the SN shock had swept up about 0.04 M_sun of material. The region emitting the H$α$ has a density of order $10^{-20}$ g/cc. △ Less

Submitted 14 October, 2019; originally announced October 2019.

Comments: 16 pages

Journal ref: 2019, ApJ, 833, 120

The Exceptional X-ray Evolution of SN1996cr in High Resolution

Authors: J. Quirola-Vasquez, F. E. Bauer, V. V. Dwarkadas, C. Badenes, W. N. Brandt, T. Nymark, D. Walton

Abstract: We present X-ray spectra spanning 18 years of evolution for SN\,1996cr, one of the five nearest SNe detected in the modern era. {\it Chandra} HETG exposures in 2000, 2004, and 2009 allow us to resolve spectrally the velocity profiles of Ne, Mg, Si, S, and Fe emission lines and monitor their evolution as tracers of the ejecta-circumstellar medium (CSM) interaction. To explain the diversity of X-ray… ▽ More We present X-ray spectra spanning 18 years of evolution for SN\,1996cr, one of the five nearest SNe detected in the modern era. {\it Chandra} HETG exposures in 2000, 2004, and 2009 allow us to resolve spectrally the velocity profiles of Ne, Mg, Si, S, and Fe emission lines and monitor their evolution as tracers of the ejecta-circumstellar medium (CSM) interaction. To explain the diversity of X-ray line profiles, we explore several possible geometrical models. Based on the highest signal-to-noise 2009 epoch, we find that a polar geometry with two distinct opening angle configurations and internal obscuration can successfully reproduce all of the observed line profiles. The best fit model consists of two plasma components: (1) a mildly absorbed (2$\times$10$^{21}$\,cm$^{-2}$), cooler ($\approx$2\,keV) with high Ne, Mg, Si, and S abundances associated with a wide polar interaction region (half-opening angle $\approx$58$^{\circ}$); (2) a moderately absorbed (2$\times$10$^{22}$\,cm$^{-2}$), hotter (>~20\,keV) plasma with high Fe abundances and strong internal obscuration associated with a narrow polar interaction region (half-opening angle $\approx$20$^{\circ}$). We extend this model to seven further epochs with lower signal-to-noise ratio and/or lower spectral-resolution between 2000-2018, yielding several interesting trends in absorption, flux, geometry and expansion velocity. We argue that the hotter and colder components are associated with reflected and forward shocks, respectively, at least at later epochs. We discuss the physical implications of our results and plausible explosion scenarios to understand the X-ray data of SN\,1996cr. △ Less

Submitted 8 October, 2019; originally announced October 2019.

Comments: 27 pages, 16 figures (+2 figures in the appendix), accepted for publication MNRAS

Excavating the Explosion and Progenitor Properties of Type IIP Supernovae via Modelling of their Optical Lightcurves

Authors: Wilson Ricks, Vikram V. Dwarkadas

Abstract: The progenitors of Type IIP supernovae (SNe) are known to be red supergiants, but their properties are not well determined. We employ hydrodynamical modelling to investigate the explosion characteristics of eight Type IIP supernovae, and the properties of their progenitor stars. We create evolutionary models using the {\sc MESA} stellar evolution code, explode these models, and simulate the optica… ▽ More The progenitors of Type IIP supernovae (SNe) are known to be red supergiants, but their properties are not well determined. We employ hydrodynamical modelling to investigate the explosion characteristics of eight Type IIP supernovae, and the properties of their progenitor stars. We create evolutionary models using the {\sc MESA} stellar evolution code, explode these models, and simulate the optical lightcurves using the {\sc STELLA} code. We fit the optical lightcurves, Fe II 5169Å velocity, and photospheric velocity, to the observational data. Recent research has suggested that the progenitors of Type IIP SNe have a zero age main sequence (ZAMS) mass not exceeding $\sim18$ M$_{\odot}$. Our fits give a progenitor ZAMS mass $\leq18$ M$_{\odot}$ for seven of the supernovae. Where previous progenitor mass estimates exist, from various sources such as hydrodynamical modelling, multi-wavelength observations, or semi-analytic calculations, our modelling generally tends towards the lower mass values. This result is in contrast to results from previous hydrodynamical modelling, but is consistent with those obtained using general-relativistic radiation-hydrodynamical codes. We do find that one event, SN 2015ba, has a progenitor whose mass is closer to 24 M$_{\odot}$ , although we are unable to fit it well. We also derive the amount of $^{56}$Ni required to reproduce the tail of the lightcurve, and find values generally larger than previous estimates. Overall, we find that it is difficult to characterize the explosion by a single parameter, and that a range of parameters is needed. △ Less

Submitted 17 June, 2019; originally announced June 2019.

Comments: 21 pages, 16 figures. Accepted to the Astrophysical Journal

Near-Earth Supernova Explosions: Evidence, Implications, and Opportunities

Authors: Brian D. Fields, John R. Ellis, Walter R. Binns, Dieter Breitschwerdt, Georgia A. de Nolfo, Roland Diehl, Vikram V. Dwarkadas, Adrienne Ertel, Thomas Faestermann, Jenny Feige, Caroline Fitoussi, Priscilla Frisch, David Graham, Brian Haley, Alexander Heger, Wolfgang Hillebrandt, Martin H. Israel, Thomas Janka, Michael Kachelriess, Gunther Korschinek, Marco Limongi, Maria Lugaro, Franciole Marinho, Adrian Melott, Richard A. Mewaldt , et al. (14 additional authors not shown)

Abstract: There is now solid experimental evidence of at least one supernova explosion within 100 pc of Earth within the last few million years, from measurements of the short-lived isotope 60Fe in widespread deep-ocean samples, as well as in the lunar regolith and cosmic rays. This is the first established example of a specific dated astrophysical event outside the Solar System having a measurable impact o… ▽ More There is now solid experimental evidence of at least one supernova explosion within 100 pc of Earth within the last few million years, from measurements of the short-lived isotope 60Fe in widespread deep-ocean samples, as well as in the lunar regolith and cosmic rays. This is the first established example of a specific dated astrophysical event outside the Solar System having a measurable impact on the Earth, offering new probes of stellar evolution, nuclear astrophysics, the astrophysics of the solar neighborhood, cosmic-ray sources and acceleration, multi-messenger astronomy, and astrobiology. Interdisciplinary connections reach broadly to include heliophysics, geology, and evolutionary biology. Objectives for the future include pinning down the nature and location of the established near-Earth supernova explosions, seeking evidence for others, and searching for other short-lived isotopes such as 26Al and 244Pu. The unique information provided by geological and lunar detections of radioactive 60Fe to assess nearby supernova explosions make now a compelling time for the astronomy community to advocate for supporting multi-disciplinary, cross-cutting research programs. △ Less

Submitted 11 March, 2019; originally announced March 2019.

Comments: 11 pages, 2 figures. Astro2020 Science White Paper submitted to the 2020 Decadal Survey on Astronomy and Astrophysics

Effect of Supernovae on the Local Interstellar Material

Authors: Priscilla Frisch, Vikram V. Dwarkadas

Abstract: A range of astronomical data indicates that ancient supernovae created the galactic environment of the Sun and sculpted the physical properties of the interstellar medium near the heliosphere. In this paper we review the characteristics of the local interstellar medium that have been affected by supernovae. The kinematics, magnetic field, elemental abundances, and configuration of the nearest inte… ▽ More A range of astronomical data indicates that ancient supernovae created the galactic environment of the Sun and sculpted the physical properties of the interstellar medium near the heliosphere. In this paper we review the characteristics of the local interstellar medium that have been affected by supernovae. The kinematics, magnetic field, elemental abundances, and configuration of the nearest interstellar material support the view that the Sun is at the edge of the Loop I superbubble, which has merged into the low density Local Bubble. The energy source for the higher temperature X-ray emitting plasma pervading the Local Bubble is uncertain. Winds from massive stars and nearby supernovae, perhaps from the Sco-Cen Association, may have contributed radioisotopes found in the geologic record and galactic cosmic ray population. Nested supernova shells in the Orion and Sco-Cen regions suggest spatially distinct sites of episodic star formation. The heliosphere properties vary with the pressure of the surrounding interstellar cloud. A nearby supernova would modify this pressure equilibrium and thereby severely disrupt the heliosphere as well as the local interstellar medium. △ Less

Submitted 18 January, 2018; originally announced January 2018.

Comments: 30 pages, 7 figures. Author version, updated and modified (several updated and new paragraphs, one new subsection), of an article that was published in the Handbook of Supernovae, A.W. Alsabti, P. Murdin (eds.), Springer, 2016

Triggered Star Formation Inside the Shell of a Wolf-Rayet Bubble as the origin of the Solar System

Authors: Vikram V. Dwarkadas, Nicolas Dauphas, Bradley Meyer, Peter Boyajian, Michael Bojazi

Abstract: A critical constraint on solar system formation is the high $^{26}$Al/$^{27}$Al abundance ratio of 5 $\times 10^{-5}$ at the time of formation, which was about 17 times higher than the average Galactic ratio, while the $^{60}$Fe/$^{56}$Fe value was about $2 \times 10^{-8}$, lower than the Galactic value. This challenges the assumption that a nearby supernova was responsible for the injection of th… ▽ More A critical constraint on solar system formation is the high $^{26}$Al/$^{27}$Al abundance ratio of 5 $\times 10^{-5}$ at the time of formation, which was about 17 times higher than the average Galactic ratio, while the $^{60}$Fe/$^{56}$Fe value was about $2 \times 10^{-8}$, lower than the Galactic value. This challenges the assumption that a nearby supernova was responsible for the injection of these short-lived radionuclides into the early solar system. We show that this conundrum can be resolved if the Solar System was formed by triggered star formation at the edge of a Wolf-Rayet (W-R) bubble. Aluminium-26 is produced during the evolution of the massive star, released in the wind during the W-R phase, and condenses into dust grains that are seen around W-R stars. The dust grains survive passage through the reverse shock and the low density shocked wind, reach the dense shell swept-up by the bubble, detach from the decelerated wind and are injected into the shell. Some portions of this shell subsequently collapses to form the dense cores that give rise to solar-type systems. The subsequent aspherical supernova does not inject appreciable amounts of $^{60}$Fe into the proto-solar-system, thus accounting for the observed low abundance of $^{60}$Fe. We discuss the details of various processes within the model and conclude that it is a viable model that can explain the initial abundances of $^{26}$Al and $^{60}$Fe. We estimate that 1-16% of all Sun-like stars could have formed in such a setting of triggered star formation in the shell of a WR bubble. △ Less

Submitted 28 December, 2017; originally announced December 2017.

Comments: 19 pages, 6 figures. Accepted version. Final published version with proof corrections can be found on the Astrophysical Journal web page for ApJ, 2017, 851, 147

Journal ref: Astrophysical Journal, 2017, Vol 851, 147

Science with the Cherenkov Telescope Array

Authors: The Cherenkov Telescope Array Consortium, :, B. S. Acharya, I. Agudo, I. Al Samarai, R. Alfaro, J. Alfaro, C. Alispach, R. Alves Batista, J. -P. Amans, E. Amato, G. Ambrosi, E. Antolini, L. A. Antonelli, C. Aramo, M. Araya, T. Armstrong, F. Arqueros, L. Arrabito, K. Asano, M. Ashley, M. Backes, C. Balazs, M. Balbo, O. Ballester , et al. (558 additional authors not shown)

Abstract: The Cherenkov Telescope Array, CTA, will be the major global observatory for very high energy gamma-ray astronomy over the next decade and beyond. The scientific potential of CTA is extremely broad: from understanding the role of relativistic cosmic particles to the search for dark matter. CTA is an explorer of the extreme universe, probing environments from the immediate neighbourhood of black ho… ▽ More The Cherenkov Telescope Array, CTA, will be the major global observatory for very high energy gamma-ray astronomy over the next decade and beyond. The scientific potential of CTA is extremely broad: from understanding the role of relativistic cosmic particles to the search for dark matter. CTA is an explorer of the extreme universe, probing environments from the immediate neighbourhood of black holes to cosmic voids on the largest scales. Covering a huge range in photon energy from 20 GeV to 300 TeV, CTA will improve on all aspects of performance with respect to current instruments. The observatory will operate arrays on sites in both hemispheres to provide full sky coverage and will hence maximize the potential for the rarest phenomena such as very nearby supernovae, gamma-ray bursts or gravitational wave transients. With 99 telescopes on the southern site and 19 telescopes on the northern site, flexible operation will be possible, with sub-arrays available for specific tasks. CTA will have important synergies with many of the new generation of major astronomical and astroparticle observatories. Multi-wavelength and multi-messenger approaches combining CTA data with those from other instruments will lead to a deeper understanding of the broad-band non-thermal properties of target sources. The CTA Observatory will be operated as an open, proposal-driven observatory, with all data available on a public archive after a pre-defined proprietary period. Scientists from institutions worldwide have combined together to form the CTA Consortium. This Consortium has prepared a proposal for a Core Programme of highly motivated observations. The programme, encompassing approximately 40% of the available observing time over the first ten years of CTA operation, is made up of individual Key Science Projects (KSPs), which are presented in this document. △ Less

Submitted 21 January, 2018; v1 submitted 22 September, 2017; originally announced September 2017.

Comments: 213 pages, including references and glossary. Version 2: credits and references updated, some figures updated, and author list updated

X-ray Emission from SN 2012ca: A Type Ia-CSM Supernova Explosion in a Dense Surrounding Medium

Authors: C. D. Bochenek, Vikram. V. Dwarkadas, Jeffrey M. Silverman, Ori D. Fox, Roger A. Chevalier, Nathan Smith, Alexei V. Filippenko

Abstract: X-ray emission is one of the signposts of circumstellar interaction in supernovae (SNe), but until now, it has been observed only in core-collapse SNe. The level of thermal X-ray emission is a direct measure of the density of the circumstellar medium (CSM), and the absence of X-ray emission from Type Ia SNe has been interpreted as a sign of a very low density CSM. In this paper, we report late-tim… ▽ More X-ray emission is one of the signposts of circumstellar interaction in supernovae (SNe), but until now, it has been observed only in core-collapse SNe. The level of thermal X-ray emission is a direct measure of the density of the circumstellar medium (CSM), and the absence of X-ray emission from Type Ia SNe has been interpreted as a sign of a very low density CSM. In this paper, we report late-time (500--800 days after discovery) X-ray detections of SN 2012ca in {\it Chandra} data. The presence of hydrogen in the initial spectrum led to a classification of Type Ia-CSM, ostensibly making it the first SN~Ia detected with X-rays. Our analysis of the X-ray data favors an asymmetric medium, with a high-density component which supplies the X-ray emission. The data suggest a number density $> 10^8$ cm$^{-3}$ in the higher-density medium, which is consistent with the large observed Balmer decrement if it arises from collisional excitation. This is high compared to most core-collapse SNe, but it may be consistent with densities suggested for some Type IIn or superluminous SNe. If SN 2012ca is a thermonuclear SN, the large CSM density could imply clumps in the wind, or a dense torus or disk, consistent with the single-degenerate channel. A remote possibility for a core-degenerate channel involves a white dwarf merging with the degenerate core of an asymptotic giant branch star shortly before the explosion, leading to a common envelope around the SN. △ Less

Submitted 23 August, 2017; originally announced August 2017.

Comments: 11 pages, 4 figures. Accepted to MNRAS

SNaX - A Database of Supernova X-ray Lightcurves

Authors: Mathias Ross, Vikram V. Dwarkadas

Abstract: We present the Supernova X-ray Database (SNaX), a compilation of the X-ray data from young supernovae (SNe). The database includes the X-ray flux and luminosity of young SNe, days to years after outburst. The original goal and intent were to present a database of Type IIn SNe. After having accomplished this we are slowly expanding it to include all SNe for which published data are available. The i… ▽ More We present the Supernova X-ray Database (SNaX), a compilation of the X-ray data from young supernovae (SNe). The database includes the X-ray flux and luminosity of young SNe, days to years after outburst. The original goal and intent were to present a database of Type IIn SNe. After having accomplished this we are slowly expanding it to include all SNe for which published data are available. The interface allows one to search for SNe using various criteria, plot all or selected data-points, and download both the data and the plot. The plotting facility allows for significant customization. There is also a facility for the user to submit data that can be directly incorporated into the database. We include an option to fit the decay of any given SN lightcurve with a power-law. The database includes a conversion of most datapoints to a common 0.3-8 keV band so that SN lightcurves may be directly compared with each other. A mailing list has been set up to disseminate information about the database. We outline the structure and function of the database, describe its various features and outline the plans for future expansion. △ Less

Submitted 19 April, 2017; originally announced April 2017.

Comments: 9 pages, 4 figures. Accepted to the Astronomical Journal

Gamma-ray observations of Tycho's SNR with VERITAS and Fermi

Authors: VERITAS Collaboration, S. Archambault, A. Archer, W. Benbow, R. Bird, E. Bourbeau, M. Buchovecky, J. H. Buckley, V. Bugaev, M. Cerruti, M. P. Connolly, W. Cui, V. V. Dwarkadas, M. Errando, A. Falcone, Q. Feng, J. P. Finley, H. Fleischhack, L. Fortson, A. Furniss, S. Griffin, M. Hutten, D. Hanna, J. Holder, C. A. Johnson , et al. (42 additional authors not shown)

Abstract: High-energy gamma-ray emission from supernova remnants (SNRs) has provided a unique perspective for studies of Galactic cosmic-ray acceleration. Tycho's SNR is a particularly good target because it is a young, type Ia SNR that is well-studied over a wide range of energies and located in a relatively clean environment. Since the detection of gamma-ray emission from Tycho's SNR by VERITAS and Fermi-… ▽ More High-energy gamma-ray emission from supernova remnants (SNRs) has provided a unique perspective for studies of Galactic cosmic-ray acceleration. Tycho's SNR is a particularly good target because it is a young, type Ia SNR that is well-studied over a wide range of energies and located in a relatively clean environment. Since the detection of gamma-ray emission from Tycho's SNR by VERITAS and Fermi-LAT, there have been several theoretical models proposed to explain its broadband emission and high-energy morphology. We report on an update to the gamma-ray measurements of Tycho's SNR with 147 hours of VERITAS and 84 months of Fermi-LAT observations, which represents about a factor of two increase in exposure over previously published data. About half of the VERITAS data benefited from a camera upgrade, which has made it possible to extend the TeV measurements toward lower energies. The TeV spectral index measured by VERITAS is consistent with previous results, but the expanded energy range softens a straight power-law fit. At energies higher than 400 GeV, the power-law index is $2.92 \pm 0.42_{\mathrm{stat}} \pm 0.20_{\mathrm{sys}}$. It is also softer than the spectral index in the GeV energy range, $2.14 \pm 0.09_{\mathrm{stat}} \pm 0.02_{\mathrm{sys}}$, measured by this study using Fermi--LAT data. The centroid position of the gamma-ray emission is coincident with the center of the remnant, as well as with the centroid measurement of Fermi--LAT above 1 GeV. The results are consistent with an SNR shell origin of the emission, as many models assume. The updated spectrum points to a lower maximum particle energy than has been suggested previously. △ Less

Submitted 24 January, 2017; originally announced January 2017.

Comments: Accepted for publication in The Astrophysical Journal

X-Ray and Radio Emission from the Luminous Supernova 2005kd

Authors: Vikram V. Dwarkadas, C. Romero-Cañizales, R. Reddy, F. E. Bauer

Abstract: SN 2005kd is among the most luminous supernovae (SNe) to be discovered at X-ray wavelengths. We have re-analysed all good angular resolution (better than $20"$ FWHM PSF) archival X-ray data for SN 2005kd. The data reveal an X-ray light curve that decreases as t$^{-1.62 \pm 0.06}$. Our modelling of the data suggests that the early evolution is dominated by emission from the forward shock in a high-… ▽ More SN 2005kd is among the most luminous supernovae (SNe) to be discovered at X-ray wavelengths. We have re-analysed all good angular resolution (better than $20"$ FWHM PSF) archival X-ray data for SN 2005kd. The data reveal an X-ray light curve that decreases as t$^{-1.62 \pm 0.06}$. Our modelling of the data suggests that the early evolution is dominated by emission from the forward shock in a high-density medium. Emission from the radiative reverse shock is absorbed by the cold dense shell formed behind the reverse shock. Our results suggest a progenitor with a mass-loss rate towards the end of its evolution of $\ge$ 4.3 $\times$ 10$^{-4} M_{\odot} \,{\rm yr}^{-1}$, for a wind velocity of 10 km s$^{-1}$, at 4.0 $\times$ 10$^{16}$ cm. This mass-loss rate is too high for most known stars, except perhaps hypergiant stars. A higher wind velocity would lead to a correspondingly higher mass-loss rate. A Luminous Blue Variable star undergoing a giant eruption could potentially fulfill this requirement, but would need a high mass-loss rate lasting for several hundred years, and need to explain the plateau observed in the optical light curve. The latter could perhaps be due to the ejecta expanding in the dense circumstellar material at relatively small radii. These observations are consistent with the fact that Type IIn SNe appear to expand into high density and high mass-loss rate environments, and also suggest rapid variability in the wind mass-loss parameters within at least the last 5000 years of stellar evolution prior to core-collapse. △ Less

Submitted 20 July, 2016; originally announced July 2016.

Comments: 11 pages, 2 tables and 4 figures. Accepted to MNRAS

Probing Final Stages of Stellar Evolution with X-Ray Observations of SN 2013ej

Authors: Sayan Chakraborti, Alak Ray, Randall Smith, Raffaella Margutti, David Pooley, Subhash Bose, Firoza Sutaria, Poonam Chandra, Vikram V. Dwarkadas, Stuart Ryder, Keiichi Maeda

Abstract: Massive stars shape their surroundings with mass loss from winds during their lifetimes. Fast ejecta from supernovae, from these massive stars, shocks this circumstellar medium. Emission generated by this interaction provides a window into the final stages of stellar evolution, by probing the history of mass loss from the progenitor. Here we use Chandra and Swift x-ray observations of the type II-… ▽ More Massive stars shape their surroundings with mass loss from winds during their lifetimes. Fast ejecta from supernovae, from these massive stars, shocks this circumstellar medium. Emission generated by this interaction provides a window into the final stages of stellar evolution, by probing the history of mass loss from the progenitor. Here we use Chandra and Swift x-ray observations of the type II-P/L SN 2013ej to probe the history of mass loss from its progenitor. We model the observed x-rays as emission from both heated circumstellar matter and supernova ejecta. The circumstellar density profile probed by the supernova shock reveals a history of steady mass loss during the final 400 years. The inferred mass loss rate of $3 \times 10^{-6} {\rm \; M_\odot \; yr^{-1}}$ points back to a 14 $M_\odot$ progenitor. Soon after the explosion we find significant absorption of reverse shock emission by a cooling shell. The column depth of this shell observed in absorption provides an independent and consistent measurement of the circumstellar density seen in emission. We also determine the efficiency of cosmic ray acceleration from x-rays produced by Inverse Compton scattering of optical photons by relativistic electrons. Only about 1 percent of the thermal energy is used to accelerate electrons. Our x-ray observations and modeling provides stringent tests for models of massive stellar evolution and micro-physics of shocks. △ Less

Submitted 25 November, 2015; v1 submitted 20 October, 2015; originally announced October 2015.

Comments: 8 pages, 6 figures, 2 tables, accepted for publication in ApJ

VERITAS Collaboration Contributions to the 34th International Cosmic Ray Conference

Authors: The VERITAS Collaboration, A. U. Abeysekara, S. Archambault, A. Archer, A. Barnacka, W. Benbow, R. Bird, J. Biteau, M. Buchovecky, J. H. Buckley, V. Bugaev, K. Byrum, J. V Cardenzana, M. Cerruti, X. Chen, J. L. Christiansen, L. Ciupik, E. Collins-Hughes, M. P. Connolly, P. Coppi, W. Cui, V. V. Dwarkadas, J. D. Eisch, M. Errando, A. Falcone , et al. (68 additional authors not shown)

Abstract: Compilation of papers presented by the VERITAS Collaboration at the 34th International Cosmic Ray Conference (ICRC), held July 30 through August 6, 2015 in The Hague, The Netherlands. Compilation of papers presented by the VERITAS Collaboration at the 34th International Cosmic Ray Conference (ICRC), held July 30 through August 6, 2015 in The Hague, The Netherlands. △ Less

Submitted 6 October, 2015; originally announced October 2015.

Comments: 1 html page. 2015 ICRC, The Hague, The Netherlands. A full list of author affiliations can be found at this link: http://veritas.sao.arizona.edu/conferences/authors?icrc2015

Science Highlights from VERITAS

Authors: D. Staszak, A. U. Abeysekara, S. Archambault, A. Archer, A. Barnacka, W. Benbow, R. Bird, J. Biteau, M. Buchovecky, J. H. Buckley, V. Bugaev, K. Byrum, J. V Cardenzana, M. Cerruti, X. Chen, J. L. Christiansen, L. Ciupik, E. Collins-Hughes, M. P. Connolly, P. Coppi, W. Cui, V. V. Dwarkadas, J. D. Eisch, M. Errando, A. Falcone , et al. (66 additional authors not shown)

Abstract: The Very Energetic Radiation Imaging Telescope Array System (VERITAS) is a ground-based array located at the Fred Lawrence Whipple Observatory in southern Arizona and is one of the world's most sensitive gamma-ray instruments at energies of 85 GeV to $>$30 TeV. VERITAS has a wide scientific reach that includes the study of extragalactic and Galactic objects as well as the search for astrophysical… ▽ More The Very Energetic Radiation Imaging Telescope Array System (VERITAS) is a ground-based array located at the Fred Lawrence Whipple Observatory in southern Arizona and is one of the world's most sensitive gamma-ray instruments at energies of 85 GeV to $>$30 TeV. VERITAS has a wide scientific reach that includes the study of extragalactic and Galactic objects as well as the search for astrophysical signatures of dark matter and the measurement of cosmic rays. In this paper, we will summarize the current status of the VERITAS observatory and present some of the scientific highlights from the last two years, focusing in particular on those results shown at the 2015 ICRC in The Hague, Netherlands. △ Less

Submitted 5 October, 2015; originally announced October 2015.

Comments: To appear in the Proceedings of the 34th International Cosmic Ray Conference (ICRC2015), The Hague, The Netherlands

Probing Efficient Cosmic-Ray Acceleration in Young Supernovae

Authors: Vikram V. Dwarkadas, M. Renaud, A. Marcowith, V. Tatischeff

Abstract: The formation of a core collapse supernovae (SNe) results in a fast (but non- or mildly-relativistic) shock wave expanding outwards into the surrounding medium. The medium itself is likely modified due to the stellar mass-loss from the massive star progenitor, which may be Wolf-Rayet stars (for Type Ib/c SNe), red supergiant stars (for type IIP and perhaps IIb and IIL SNe), or some other stellar t… ▽ More The formation of a core collapse supernovae (SNe) results in a fast (but non- or mildly-relativistic) shock wave expanding outwards into the surrounding medium. The medium itself is likely modified due to the stellar mass-loss from the massive star progenitor, which may be Wolf-Rayet stars (for Type Ib/c SNe), red supergiant stars (for type IIP and perhaps IIb and IIL SNe), or some other stellar type. The wind mass-loss parameters determine the density structure of the surrounding medium. Combined with the velocity of the SN shock wave, this regulates the shock acceleration process. In this article we discuss the essential parameters that control the particle acceleration and gamma-ray emission in SNe, with particular reference to the Type IIb SN 1993J. The shock wave expanding into the high density medium leads to fast particle acceleration, giving rise to rapidly-growing plasma instabilities driven by the acceleration process itself. The instabilities grow over intraday timescales. This growth, combined with the interplay of non-linear processes, results in the amplification of the magnetic field at the shock front, which can adequately account for the magnetic field strengths deduced from radio monitoring of the source. The maximum particle energy can reach, and perhaps exceed, 1 PeV, depending on the dominant instability. The gamma-ray signal is found to be heavily absorbed by pair production process during the first week after the outburst. We derive the time dependent particle spectra and associated hadronic signatures of secondary particles (gamma-ray, leptons and neutrinos) arising from proton proton interactions. We find that the Cherenkov Telescope Array (CTA) should be able to detect objects like SN 1993J above 1 TeV. We predict a low neutrino flux above 10 TeV, implying a detectability horizon with current or planned neutrino telescopes of 1 Mpc. △ Less

Submitted 2 September, 2015; originally announced September 2015.

Comments: 8 pages, 1 figure. To be published in Proceedings of Science, as part of the proceedings of the 34th International Cosmic Ray Conference (ICRC), held in the Hague, Netherlands. arXiv admin note: text overlap with arXiv:1409.3670

X-ray Emission from Ionized Wind-Bubbles around Wolf-Rayet Stars

Authors: Vikram V. Dwarkadas, Duane Rosenberg

Abstract: Using a code that employs a self-consistent method for computing the effects of photoionization on circumstellar gas dynamics, we model the formation of wind-driven nebulae around massive Wolf-Rayet (W-R) stars. Our algorithm incorporates a simplified model of the photo-ionization source, computes the fractional ionization of hydrogen due to the photoionizing flux and recombination, and determines… ▽ More Using a code that employs a self-consistent method for computing the effects of photoionization on circumstellar gas dynamics, we model the formation of wind-driven nebulae around massive Wolf-Rayet (W-R) stars. Our algorithm incorporates a simplified model of the photo-ionization source, computes the fractional ionization of hydrogen due to the photoionizing flux and recombination, and determines self-consistently the energy balance due to ionization, photo-heating and radiative cooling. We take into account changes in stellar properties and mass-loss over the star's evolution. Our multi-dimensional simulations clearly reveal the presence of strong ionization front instabilities. Using various X-ray emission models, and abundances consistent with those derived for W-R nebulae, we compute the X-ray flux and spectra from our wind bubble models. We show the evolution of the X-ray spectral features with time over the evolution of the star, taking the absorption of the X-rays by the ionized bubble into account. Our simulated X-ray spectra compare reasonably well with observed spectra of Wolf-Rayet bubbles. They suggest that X-ray nebulae around massive stars may not be easily detectable, consistent with observations. △ Less

Submitted 19 August, 2015; originally announced August 2015.

Comments: 3 pages, 4 figures. To be published in the Proceedings of the International Workshop on Wolf-Rayet Stars, Potsdam, 1-5 June 2015. Editors W.-R. Hamann, A. Sander, and H. Todt

Analysis of GeV-band gamma-ray emission from SNR RX J1713.7-3946

Authors: S. Federici, M. Pohl, I. Telezhinsky, A. Wilhelm, V. V. Dwarkadas

Abstract: RX J1713.7-3946 is the brightest shell-type Supernova remnant (SNR) of the TeV gamma-ray sky. Earlier Fermi-LAT results on low-energy gamma-ray emission suggested that, despite large uncertainties in the background determination, the spectrum is inconsistent with a hadronic origin. We update the GeV-band spectra using improved estimates for the diffuse galactic gamma-ray emission and more than dou… ▽ More RX J1713.7-3946 is the brightest shell-type Supernova remnant (SNR) of the TeV gamma-ray sky. Earlier Fermi-LAT results on low-energy gamma-ray emission suggested that, despite large uncertainties in the background determination, the spectrum is inconsistent with a hadronic origin. We update the GeV-band spectra using improved estimates for the diffuse galactic gamma-ray emission and more than doubled data volume. We further investigate the viability of hadronic emission models for RX J1713.7-3946. We produced a high-resolution map of the diffuse Galactic gamma-ray background corrected for HI self-absorption and used it in the analysis of more than 5~years worth of Fermi-LAT data. We used hydrodynamic scaling relations and a kinetic transport equation to calculate the acceleration and propagation of cosmic-rays in SNR. We then determined spectra of hadronic gamma-ray emission from RX J1713.7-3946, separately for the SNR interior and the cosmic-ray precursor region of the forward shock, and computed flux variations that would allow to test the model with observations. We find that RX J1713.7-3946 is now detected by Fermi-LAT with very high statistical significance, and the source morphology is best described by that seen in the TeV band. The measured spectrum of RX J1713.7-3946 is hard with index gamma=1.53 +/- 0.07, and the integral flux above 500 MeV is F = (5.5 +/- 1.1)e-9 photons/cm^2/s. We demonstrate that scenarios based on hadronic emission from the cosmic-ray precursor region are acceptable for RX J1713.7-3946, and we predict a secular flux increase at a few hundred GeV at the level of around 15% over 10 years, which may be detectable with the upcoming CTA observatory. △ Less

Submitted 23 February, 2015; originally announced February 2015.

Comments: 9 pages, accepted for publication in Astronomy & Astrophysics

Journal ref: A&A 577, A12 (2015)

Investigating the TeV Morphology of MGRO J1908+06 with VERITAS

Authors: E. Aliu, S. Archambault, T. Aune, B. Behera, M. Beilicke, W. Benbow, K. Berger, R. Bird, J. H. Buckley, V. Bugaev, J. V Cardenzana, M. Cerruti, X. Chen, L. Ciupik, E. Collins-Hughes, M. P. Connolly, W. Cui, J. Dumm, V. V. Dwarkadas, M. Errando, A. Falcone, S. Federici, Q. Feng, J. P. Finley, H. Fleischhack , et al. (63 additional authors not shown)

Abstract: We report on deep observations of the extended TeV gamma-ray source MGRO J1908+06 made with the VERITAS very high energy (VHE) gamma-ray observatory. Previously, the TeV emission has been attributed to the pulsar wind nebula (PWN) of the Fermi-LAT pulsar PSR J1907+0602. We detect MGRO J1908+06 at a significance level of 14 standard deviations (14 sigma) and measure a photon index of 2.20 +/- 0.10_… ▽ More We report on deep observations of the extended TeV gamma-ray source MGRO J1908+06 made with the VERITAS very high energy (VHE) gamma-ray observatory. Previously, the TeV emission has been attributed to the pulsar wind nebula (PWN) of the Fermi-LAT pulsar PSR J1907+0602. We detect MGRO J1908+06 at a significance level of 14 standard deviations (14 sigma) and measure a photon index of 2.20 +/- 0.10_stat +/- 0.20_sys. The TeV emission is extended, covering the region near PSR J1907+0602 and also extending towards SNR G40.5--0.5. When fitted with a 2-dimensional Gaussian, the intrinsic extension has a standard deviation of sigma_src = 0.44 +/- 0.02 degrees. In contrast to other TeV PWNe of similar age in which the TeV spectrum softens with distance from the pulsar, the TeV spectrum measured near the pulsar location is consistent with that measured at a position near the rim of G40.5--0.5, 0.33 degrees away. △ Less

Submitted 28 April, 2014; originally announced April 2014.

Comments: To appear in ApJ, 8 pages

On the Lack of X-ray Bright Type IIP Supernovae

Authors: Vikram V. Dwarkadas

Abstract: Type IIP Supernovae (SNe) are expected to arise from Red Supergiant stars (RSGs). These stars have observed mass-loss rates that span more than two orders of magnitude, from $< 10^{-6}$ solar masses yr$^{-1}$ to almost $ 10^{-4} $ solar masses yr$^{-1}$. Thermal bremsstrahlung X-ray emission from at least some IIP's should reflect the larger end of the high mass-loss rates. Strangely, no IIP SNe a… ▽ More Type IIP Supernovae (SNe) are expected to arise from Red Supergiant stars (RSGs). These stars have observed mass-loss rates that span more than two orders of magnitude, from $< 10^{-6}$ solar masses yr$^{-1}$ to almost $ 10^{-4} $ solar masses yr$^{-1}$. Thermal bremsstrahlung X-ray emission from at least some IIP's should reflect the larger end of the high mass-loss rates. Strangely, no IIP SNe are seen where the X-ray luminosity is large enough to suggest mass-loss rates greater than about $ 10^{-5} $ solar masses yr$^{-1}$. We investigate if this could be due to absorption of the X-ray emission. After carefully studying all the various aspects, we conclude that absorption would not be large enough to prevent us from having detected X-ray emission from high mass-loss rate IIP's. This leads us to the conclusion that there may be an upper limit of $\sim 10^{-5} $ solar masses yr$^{-1}$ to the mass-loss rate of Type IIP progenitors, and therefore to the luminosity of RSGs that explode to form Type IIPs. This is turn suggests an upper limit of $\leq 19 $ solar masses for the progenitor mass of a Type IIP SN. This limit is close to that obtained by direct detection of IIP progenitors, as well as that suggested by recent stellar evolution calculations. Although the statistics need to be improved, many current indicators support the notion that RSGs above $\sim 19 $ solar masses do not explode to form Type IIP SNe. △ Less

Submitted 20 February, 2014; originally announced February 2014.

Comments: 9 pages, 1 figure, 1 table. Accepted to MNRAS

VERITAS contributions to the 33rd International Cosmic Ray Conference

Authors: VERITAS Collaboration, E. Aliu, S. Archambault, T. Aune, B. Behera, M. Beilicke, W. Benbow, K. Berger, R. Bird, A. Bouvier, J. H. Buckley, V. Bugaev, K. Byrum, J. V Cardenzana, M. Cerruti, L. Ciupik, M. P. Connolly, W. Cui, J. Dumm, V. V. Dwarkadas, M. Errando, A. Falcone, S. Federici, Q. Feng, J. P. Finley , et al. (64 additional authors not shown)

Abstract: Compilation of papers contributed by the VERITAS Collaboration to the 33rd International Cosmic Ray Conference, held 2-9 July, 2013, in Rio de Janeiro, Brazil. Compilation of papers contributed by the VERITAS Collaboration to the 33rd International Cosmic Ray Conference, held 2-9 July, 2013, in Rio de Janeiro, Brazil. △ Less

Submitted 28 August, 2013; originally announced August 2013.

Exploring the Gamma-Ray Emissivity of Young Supernova Remnants I: Hadronic Emission

Authors: Vikram V Dwarkadas

Abstract: Using a simplified model for the hadronic emission from young supernova remnants (SNRs), we derive an expression to calculate the hadronic luminosity with time, depending on the SN ejecta density profile and the density structure of the surrounding medium. Our analysis shows that the hadronic emission will decrease with time for core-collapse SNe expanding in the winds of their progenitor stars, b… ▽ More Using a simplified model for the hadronic emission from young supernova remnants (SNRs), we derive an expression to calculate the hadronic luminosity with time, depending on the SN ejecta density profile and the density structure of the surrounding medium. Our analysis shows that the hadronic emission will decrease with time for core-collapse SNe expanding in the winds of their progenitor stars, but increase with time for SNe expanding into a constant density medium, typical of Type Ia SNe. Using our expressions, we can compute the time-dependent hadronic flux from some well-known young SNe and SNRs with time, and where applicable reproduce previous results in the appropriate parameter regime. Using our calculations, we also emphasize the exciting possibility that SN 1987A may become a visible gamma-ray source in the next decade. △ Less

Submitted 16 July, 2013; originally announced July 2013.

Comments: Accepted to MNRAS. 11 pages, 0 figures

Discovery of TeV Gamma-ray Emission Toward Supernova Remnant SNR G78.2+2.1

Authors: E. Aliu, S. Archambault, T. Arlen, T. Aune, M. Beilicke, W. Benbow, R. Bird, A. Bouvier, S. M. Bradbury, J. H. Buckley, V. Bugaev, K. Byrum, A. Cannon, A. Cesarini, L. Ciupik, E. Collins-Hughes, M. P. Connolly, W. Cui, R. Dickherber, C. Duke, J. Dumm, V. V. Dwarkadas, M. Errando, A. Falcone, S. Federici , et al. (75 additional authors not shown)

Abstract: We report the discovery of an unidentified, extended source of very-high-energy (VHE) gamma-ray emission, VER J2019+407, within the radio shell of the supernova remnant SNR G78.2+2.1, using 21.4 hours of data taken by the VERITAS gamma-ray observatory in 2009. These data confirm the preliminary indications of gamma-ray emission previously seen in a two-year (2007-2009) blind survey of the Cygnus r… ▽ More We report the discovery of an unidentified, extended source of very-high-energy (VHE) gamma-ray emission, VER J2019+407, within the radio shell of the supernova remnant SNR G78.2+2.1, using 21.4 hours of data taken by the VERITAS gamma-ray observatory in 2009. These data confirm the preliminary indications of gamma-ray emission previously seen in a two-year (2007-2009) blind survey of the Cygnus region by VERITAS. VER J2019+407, which is detected at a post-trials significance of 7.5 standard deviations in the 2009 data, is localized to the northwestern rim of the remnant in a region of enhanced radio and X-ray emission. It has an intrinsic extent of 0.23^{\circ} \pm 0.03^{\circ} (stat)+0.04^{\circ}_{-0.02}^{\circ}(sys) and its spectrum is well-characterized by a differential power law (dN/dE = N_0 \times (E/TeV)^{-Γ}) with a photon index of Γ = 2.37 \pm 0.14 (stat) \pm 0.20 (sys) and a flux normalization of N0 = 1.5 \pm 0.2 (stat) \pm 0.4(sys) \times 10^-12 ph TeV^{-1} cm^{-2} s^{-1}. This yields an integral flux of 5.2 \pm 0.8 (stat) \pm 1.4 (sys) \times 10^-12 ph cm^{-2} s^{-1} above 320 GeV, corresponding to 3.7% of the Crab Nebula flux. We consider the relationship of the TeV gamma-ray emission with the GeV gamma-ray emission seen from SNR G78.2+2.1 as well as that seen from a nearby cocoon of freshly accelerated cosmic rays. Multiple scenarios are considered as possible origins for the TeV gamma-ray emission, including hadronic particle acceleration at the supernova remnant shock. △ Less

Submitted 28 May, 2013; originally announced May 2013.

Comments: 18 pages, 5 figures, accepted for publication in ApJ

The progenitor of SN 2011ja: Clues from circumstellar interaction

Authors: Sayan Chakraborti, Alak Ray, Randall Smith, Stuart Ryder, Naveen Yadav, Firoza Sutaria, Vikram V. Dwarkadas, Poonam Chandra, David Pooley, Rupak Roy

Abstract: Massive stars, possibly red supergiants, which retain extended hydrogen envelopes until core collapse, produce Type II Plateau (IIP) supernovae. The ejecta from these explosions shock the circumstellar matter originating from the mass loss of the progenitor during the final phases of its life. This interaction accelerates particles to relativistic energies which then lose energy via synchrotron ra… ▽ More Massive stars, possibly red supergiants, which retain extended hydrogen envelopes until core collapse, produce Type II Plateau (IIP) supernovae. The ejecta from these explosions shock the circumstellar matter originating from the mass loss of the progenitor during the final phases of its life. This interaction accelerates particles to relativistic energies which then lose energy via synchrotron radiation in the shock-amplified magnetic fields and inverse Compton scattering against optical photons from the supernova. These processes produce different signatures in the radio and X-ray part of the electromagnetic spectrum. Observed together, they allow us to break the degeneracy between shock acceleration and magnetic field amplification. In this work we use X-rays observations from the Chandra and radio observations from the ATCA to study the relative importance of processes which accelerate particles and those which amplify magnetic fields in producing the non-thermal radiation from SN 2011ja. We use radio observations to constrain the explosion date. Multiple Chandra observations allow us to probe the history of variable mass loss from the progenitor. The ejecta expands into a low density bubble followed by interaction with a higher density wind from a red supergiant consistent with ZAMS mass greater than 16 solar masses. Our results suggest that a fraction of type IIP supernovae may interact with circumstellar media set up by non-steady winds. △ Less

Submitted 11 March, 2013; v1 submitted 27 February, 2013; originally announced February 2013.

Comments: 9 pages, 5 figures, updated

Performing a stellar autopsy using the radio-bright remnant of SN 1996cr

Authors: C. Meunier, F. E. Bauer, V. V. Dwarkadas, B. Koribalski, B. Emonts, R. W. Hunstead, D. Campbell-Wilson, C. Stockdale, S. J. Tingay

Abstract: We present archival radio observations of SN1996cr from ATCA and MOST, and model its radio lightcurves using X-ray constrained hydrodynamical simulations of the interaction between the SN ejecta and the CSM. The early radio data show signs of free-free absorption (FFA), which decreases gradually and is minimal above 1.4 GHz after day ~3000. FFA optical depth constraints provide estimates of the CS… ▽ More We present archival radio observations of SN1996cr from ATCA and MOST, and model its radio lightcurves using X-ray constrained hydrodynamical simulations of the interaction between the SN ejecta and the CSM. The early radio data show signs of free-free absorption (FFA), which decreases gradually and is minimal above 1.4 GHz after day ~3000. FFA optical depth constraints provide estimates of the CSM free electron density, which allows insight into the ionisation of SN1996cr's CSM and offers a test on the density distribution adopted by the hydrodynamical simulation. The intrinsic spectral index of the radiation shows evidence for spectral flattening, which is characterised by alpha = 0.852 +/- 0.002 at day 3000 and a decay rate of d_alpha = -0.014 +/- 0.001 yr^-1. The similarity in the spectral flattening of SN1987A, SN1993J, and SN1996cr suggests this may be a relatively common feature of SNe/CSM shocks. We adopt this spectral index variation to model the synchrotron radio emission of the shock, and consider several scalings that relate the parameters of the hydrodynamical simulation to the magnetic field and electron distribution. The simulated light curves match the large-scale features of the observed light curves, but fail to match certain tightly constraining sections. This suggests that simple energy density scalings may not be able to account for the complexities of the true physical processes at work, or alternatively, that the parameters of the simulation require modification in order to accurately represent SN1996cr's CSM. △ Less

Submitted 21 February, 2013; originally announced February 2013.

Comments: 17 pages, 6 figures, 5 tables, accepted for publication in MNRAS