Source author record

V. V. Dwarkadas

V. V. Dwarkadas appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

20works
4topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

20 published item(s)

preprint2022arXiv

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

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.

preprint2021arXiv

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

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.

preprint2020arXiv

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

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.

preprint2020arXiv

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

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.

preprint2015arXiv

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

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.

preprint2015arXiv

Science Highlights from VERITAS

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.

preprint2014arXiv

Investigating the TeV Morphology of MGRO J1908+06 with VERITAS

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.

preprint2013arXiv

Acceleration of cosmic rays by young core-collapse supernova remnants

Context. Supernova remnants (SNRs) are thought to be the primary candidates for the sources of Galactic cosmic rays. According to the diffusive shock acceleration theory, SNR shocks produce a power-law spectrum with an index of s = 2, perhaps nonlinearly modified to harder spectra at high energy. Observations of SNRs often indicate particle spectra that are softer than that and show features not expected from classical theory. Known drawbacks of the standard approach are the assumption that SNRs evolve in a uniform environment, and that the reverse shock does not accelerate particles. Relaxing these assumptions increases the complexity of the problem, because one needs reliable hydrodynamical data for the plasma flow as well as good estimates for the magnetic field (MF) at the reverse shock. Aims. We show that these two factors are especially important when modeling young core-collapse SNRs that evolve in a complicated circumstellar medium shaped by the winds of progenitor stars. Methods. We used high-resolution numerical simulations for the hydrodynamical evolution of the SNR. Instead of parametrizations of the MF profiles inside the SNR, we followed the advection of the frozen-in MF inside the SNR, and thus obtained the B-field value at all locations, in particular at the reverse shock. To model cosmic-ray acceleration we solved the cosmic-ray transport equation in test-particle approximation. Results. We find that the complex plasma-flow profiles of core-collapse SNRs significantly modify the particle spectra. Additionally, the reverse shock strongly affects the emission spectra and the surface brightness.

preprint2013arXiv

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

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.

preprint2013arXiv

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

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.

preprint2012arXiv

Discovery of TeV Gamma-ray Emission from CTA 1 by VERITAS

We report the discovery of TeV gamma-ray emission coincident with the shell-type radio supernova remnant (SNR) CTA 1 using the VERITAS gamma-ray observatory. The source, VER J0006+729, was detected as a 6.5 standard deviation excess over background and shows an extended morphology, approximated by a two-dimensional Gaussian of semi-major (semi-minor) axis 0.30 degree (0.24 degree) and a centroid 5' from the Fermi gamma-ray pulsar PSR J0007+7303 and its X-ray pulsar wind nebula (PWN). The photon spectrum is well described by a power-law dN/dE = N_0 (E/3 TeV)^(-Γ), with a differential spectral index of Γ= 2.2 +- 0.2_stat +- 0.3_sys, and normalization N_0 = (9.1 +- 1.3_stat +- 1.7_sys) x 10^(-14) cm^(-2) s^(-1) TeV^(-1). The integral flux, F_γ= 4.0 x 10^(-12) erg cm^(-2) s^(-1) above 1 TeV, corresponds to 0.2% of the pulsar spin-down power at 1.4 kpc. The energetics, co-location with the SNR, and the relatively small extent of the TeV emission strongly argue for the PWN origin of the TeV photons. We consider the origin of the TeV emission in CTA 1.

preprint2012arXiv

Evolution and Hydrodynamics of the Very-Broad X-ray Line Emission in SN 1987A

Observations of SN 1987A by the Chandra High Energy Transmission Grating (HETG) in 1999 and the XMM-Newton Reflection Grating Spectrometer (RGS) in 2003 show very broad (v-b) lines with a full-width at half-maximum (FWHM) of order 10^4 kms; at these times the blast wave was primarily interacting with the HII region around the progenitor. Since then, the X-ray emission has been increasingly dominated by narrower components as the blast wave encounters dense equatorial ring (ER) material. Even so, continuing v-b emission is seen in the grating spectra suggesting that interaction with HII region material is on-going. Based on the deep HETG 2007 and 2011 data sets, and confirmed by RGS and other HETG observations, the v-b component has a width of 9300 +/-2000 kms FWHM and contributes of order 20% of the current 0.5--2 keV flux. Guided by this result, SN 1987A's X-ray spectra are modeled as the weighted sum of the non-equilibrium-ionization (NEI) emission from two simple 1D hydrodynamic simulations, this "2x1D" model reproduces the observed radii, light curves, and spectra with a minimum of free parameters. The interaction with the HII region (rho_init \sim 130 amu/cc, +/- 15 degrees opening angle) produces the very-broad emission lines and most of the 3-10 keV flux. Our ER hydrodynamics, admittedly a crude approximation to the multi-D reality, gives ER densities of order 10^4 amu/cc, requires dense clumps (x5.5 density enhancement in \sim 30% of the volume), and it predicts that the 0.5-2 keV flux will drop at a rate of \sim 17% per year once no new dense ER material is being shocked.

preprint2012arXiv

On the Evolution of the Maximum Energy of Accelerated Particles in Young Supernova Remnants

It has generally been assumed in the literature that while young supernova remnants (SNRs) accelerate particles even in the early stages, the particles do not escape until the start of the Sedov-Taylor or adiabatic stage, when the maximum energy of accelerated particles is reached. These calculations however do not take into account the detailed hydrodynamical expansion in the ejecta-dominated stage, and the approach to the Sedov stage. Using analytic approximations, we explore different environments in which the SNR may evolve, and investigate how the maximum energy to which particles are accelerated, and its time evolution, depends on various parameters. We take into account the ambient magnetic field and its amplification by resonant or non-resonant modes. Our studies reveal that the maximum energy to which particles are accelerated is generally reached in the ejecta-dominated stage, much before the start of Sedov stage. For SNe evolving within the winds of their massive stars, the maximum energy is reached very early in the evolution. We briefly explore the consequences for supernova remnants expanding in surroundings such as wind-bubbles or superbubbles.

preprint2012arXiv

Time-Dependent Escape of Cosmic Rays from Supernova Remnants, and their Interaction with Dense Media

Context. Supernova remnants (SNRs) are thought to be the main source of Galactic cosmic rays (CRs) up to the "knee" in CR spectrum. During the evolution of a SNR, the bulk of the CRs are confined inside the SNR shell. The highest-energy particles leave the system continuously, while the remaining adiabatically cooled particles are released when the SNR has expanded sufficiently and decelerated so that the magnetic field at the shock is no longer able to confine them. Particles escaping from the parent system may interact with nearby molecular clouds, producing $γ$-rays in the process via pion decay. The soft gamma-ray spectra observed for a number of SNRs interacting with molecular clouds, however, challenge current theories of non-linear particle acceleration that predict harder spectra. Aims. We study how the spectrum of escaped particles depends on the time-dependent acceleration history in both Type Ia and core-collapse SNRs, as well as on different assumptions about the diffusion coefficient in the vicinity of the SNR. Methods. We solve the CR transport equation in a test-particle approach combined with numerical simulations of SNR evolution. Results. We extend our method for calculating the cosmic-ray acceleration in SNRs to trace the escaped particles in a large volume around SNRs. We calculate the evolution of the spectra of CRs that have escaped from a SNR into a molecular cloud or dense shell for two diffusion models. We find a strong confinement of CRs in a close region around the SNR, and a strong dilution effect for CRs that were able to propagate out as far as a few SNR radii.

preprint2012arXiv

VERITAS Observations of the Nova in V407 Cygni

We report on very high energy (E > 100 GeV) gamma-ray observations of V407 Cygni, a symbiotic binary that underwent a nova outburst producing 0.1-10 GeV gamma rays during 2010 March 10-26. Observations were made with the Very Energetic Radiation Imaging Telescope Array System during 2010 March 19-26 at relatively large zenith angles, due to the position of V407 Cyg. An improved reconstruction technique for large zenith angle observations is presented and used to analyze the data. We do not detect V407 Cygni and place a differential upper limit on the flux at 1.6 TeV of 2.3 \times 10^(-12) erg cm^(-2) s^(-1) (at the 95% confidence level). When considered jointly with data from Fermi-LAT, this result places limits on the acceleration of very high energy particles in the nova.

preprint2011arXiv

Radio Imaging of SN 1993J: The Story Continues

We present the most recent VLBI images of SN 1993J, taken at 1.7 GHz on 2010 March 5-6, along with a discussion of its evolution with time. The new image is the latest in a sequence covering almost the entire lifetime of the supernova. For these latest observations we used an "in beam calibrator" technique, and obtained a background rms brightness of 3.7 micro-Jy/beam. The supernova shell remains quite circular in outline. Modulations in brightness are seen around the rim which evolve relatively slowly, having remained generally similar over the last several years of observation. We determine the outer radius of the supernova using visibility-plane model-fitting. The supernova has slowed down to around 30% of its original expansion velocity, and continues to expand with radius approximately proportional to t^0.8, however, deviations from a strict power-law evolution are seen. We do not find any clear-cut evidence for systematically frequency-dependent evolution, suggesting that the radii as determined from visibility-plane model-fitting continue to provide reasonable estimates of the physical outer shock-front radius.

preprint2009arXiv

Observation of Extended VHE Emission from the Supernova Remnant IC 443 with VERITAS

We present evidence that the very-high-energy (VHE, E > 100 GeV) gamma-ray emission coincident with the supernova remnant IC 443 is extended. IC 443 contains one of the best-studied sites of supernova remnant/molecular cloud interaction and the pulsar wind nebula CXOU J061705.3+222127, both of which are important targets for VHE observations. VERITAS observed IC 443 for 37.9 hours during 2007 and detected emission above 300 GeV with an excess of 247 events, resulting in a significance of 8.3 standard deviations (sigma) before trials and 7.5 sigma after trials in a point-source search. The emission is centered at 06 16 51 +22 30 11 (J2000) +- 0.03_stat +- 0.08_sys degrees, with an intrinsic extension of 0.16 +- 0.03_stat +- 0.04_sys degrees. The VHE spectrum is well fit by a power law (dN/dE = N_0 * (E/TeV)^-Gamma) with a photon index of 2.99 +- 0.38_stat +- 0.3_sys and an integral flux above 300 GeV of (4.63 +- 0.90_stat +- 0.93_sys) * 10^-12 cm^-2 s^-1. These results are discussed in the context of existing models for gamma-ray production in IC 443.