Showing 1–41 of 41 results for author: Vurm, I

Superbubbles as the source of dynamical friction: gas migration, stellar and dark matter contributions

Authors: Rain Kipper, Indrek Vurm, Aikaterini Niovi Triantafyllaki, Peeter Tenjes, Elmo Tempel

Abstract: The gas distribution in galaxies is smooth on large scales but is usually inhomogeneous as well as time-dependent on smaller scales. The time-dependence originates from processes such as cloud formation, their collisions and supernovae (SNe) explosions. The inhomogeneities in the matter distribution give rise to variations in the local galactic gravitational potential, which can contribute towards… ▽ More The gas distribution in galaxies is smooth on large scales but is usually inhomogeneous as well as time-dependent on smaller scales. The time-dependence originates from processes such as cloud formation, their collisions and supernovae (SNe) explosions. The inhomogeneities in the matter distribution give rise to variations in the local galactic gravitational potential, which can contribute towards dynamically coupling the gas disk to the stellar and the dark matter (DM) components of the galaxy. Specifically, multiple supernovae occurring in young stellar clusters give rise to superbubbles (SB), which modify the local acceleration field and alter the energy and momentum of stars of DM particles traversing them, in broad analogy to the dynamical friction caused by a massive object. We aim to quantify how the acceleration field from SBs causes dynamical friction and contributes to the secular evolution of galaxies. In order to assess this, we construct the density modifications to the gas distribution that mimics a SB. By evaluating the acceleration field from these density modifications, we see how the momentum or angular momentum of the gas hosting the SBs changes when stars pass through the SB. Combining the effects of all the stars and SBs we, construct an empirical approximation formula for the momentum loss in homogeneous and isotropic cases. We find that the rate at which the gas disc loses its specific angular momentum via the above process is up to 4% per Gyr, which translates to under half of its original value over the lifetime of the disc. Finally, we studied how the dynamical coupling of the gas disk with the DM halo depends on assumptions on the halo kinematics (e.g. rotation) and found a ~0.3% variation in the gas disc secular evolution between different DM kinematic models. △ Less

Submitted 28 April, 2024; originally announced April 2024.

Comments: 14 pages, 12 figures

TESS photometry of the nova eruption in V606 Vul: asymmetric photosphere and multiple ejections?

Authors: Kirill V. Sokolovsky, Elias Aydi, Konstantin Malanchev, Colin J. Burke, Koji Mukai, J. L. Sokoloski, Brian D. Metzger, Kirill E. Atapin, Aleksandre A. Belinski, Yu-Ching Chen, Laura Chomiuk, Pavol A. Dubovsky, Claude-Andre Faucher-Giguere, Rebekah A. Hounsell, Natalia P. Ikonnikova, Vsevolod Yu. Lander, Junyao Li, Justin D. Linford, Amy J. Mioduszewski, Isabella Molina, Ulisse Munari, Sergey A. Potanin, Robert M. Quimby, Michael P. Rupen, Simone Scaringi , et al. (46 additional authors not shown)

Abstract: Lightcurves of many classical novae deviate from the canonical "fast rise - smooth decline" pattern and display complex variability behavior. We present the first TESS-space-photometry-based investigation of this phenomenon. We use TESS Sector 41 full-frame images to extract a lightcurve of the slow Galactic nova V606 Vul that erupted nine days prior to the start of the TESS observations. The ligh… ▽ More Lightcurves of many classical novae deviate from the canonical "fast rise - smooth decline" pattern and display complex variability behavior. We present the first TESS-space-photometry-based investigation of this phenomenon. We use TESS Sector 41 full-frame images to extract a lightcurve of the slow Galactic nova V606 Vul that erupted nine days prior to the start of the TESS observations. The lightcurve covers the first of two major peaks of V606 Vul that was reached 19 days after the start of the eruption. The nova reached its brightest visual magnitude V=9.9 in its second peak 64 days after the eruption onset, following the completion of Sector 41 observations. To increase the confidence level of the extracted lightcurve, we performed the analysis using four different codes implementing the aperture photometry (Lightkurve, VaST) and image subtraction (TESSreduce, tequila_shots) and find good agreement between them. We performed ground-based photometric and spectroscopic monitoring to complement the TESS data. The TESS lightcurve reveals two features: periodic variations (0.12771d, 0.01mag average peak-to-peak amplitude) that disappeared when the source was within 1mag of peak optical brightness and a series of isolated mini-flares (with peak-to-peak amplitudes of up to 0.5mag) appearing at seemingly random times. We interpret the periodic variations as the result of azimuthal asymmetry of the photosphere engulfing the nova-hosting binary that was distorted by and rotating with the binary. Whereas we use spectra to associate the two major peaks in the nova lightcurve with distinct episodes of mass ejection, the origin of mini-flares remains elusive. △ Less

Submitted 8 November, 2023; originally announced November 2023.

Comments: 31 pages, 10 figures, submitted to ApJ; comments welcome

Back to the present: A general treatment for the tidal field from the wake of dynamical friction

Authors: Rain Kipper, Peeter Tenjes, María Benito, Punyakoti Ganeshaiah Veena, Aikaterini Niovi Triantafyllaki, Indrek Vurm, Moorits Mihkel Muru, Maret Einasto, Elmo Tempel

Abstract: Dynamical friction can be a valuable tool for inferring dark matter properties that are difficult to constrain by other methods. Most applications of dynamical friction calculations are concerned with the long-term angular momentum loss and orbital decay of the perturber within its host. This, however, assumes knowledge of the unknown initial conditions of the system. We advance an alternative met… ▽ More Dynamical friction can be a valuable tool for inferring dark matter properties that are difficult to constrain by other methods. Most applications of dynamical friction calculations are concerned with the long-term angular momentum loss and orbital decay of the perturber within its host. This, however, assumes knowledge of the unknown initial conditions of the system. We advance an alternative methodology to infer the host properties from the perturber's shape distortions induced by the tides of the wake of dynamical friction, which we refer to as the tidal dynamical friction. As the shape distortions rely on the tidal field that has a predominantly local origin, we present a strategy to find the local wake by integrating the stellar orbits back in time along with the perturber, then removing the perturber's potential and re-integrating them back to the present. This provides perturbed and unperturbed coordinates and hence a change in coordinates, density, and acceleration fields, which yields the back-reaction experienced by the perturber. The method successfully recovers the tidal field of the wake based on a comparison with N-body simulations. We show that similar to the tidal field itself, the noise and randomness of the dynamical friction force due to the finite number of stars is also dominated by regions close to the perturber. Stars near the perturber influence it more but are smaller in number, causing a high variance in the acceleration field. These fluctuations are intrinsic to dynamical friction. We show that a stellar density of $0.0014 {\rm M_\odot\, kpc^{-3}}$ yields an inherent variance of 10% to the dynamical friction. The current method extends the family of dynamical friction methods that allow for the inference of host properties from tidal forces of the wake. It can be applied to specific galaxies, such as Magellanic Clouds, with Gaia data. △ Less

Submitted 15 November, 2023; v1 submitted 7 November, 2023; originally announced November 2023.

Comments: Accepted to Astronomy and Astrophysics

Cosmic gas highways in C-EAGLE simulations

Authors: I. Vurm, J. Nevalainen, S. E. Hong, Y. M. Bahé, C. Dalla Vecchia, P. Heinämäki

Abstract: A substantial fraction of the cosmic baryons is expected to hide in the form of diffuse warm-hot intergalactic medium (WHIM), the majority of which resides in the filaments of the Cosmic Web and has proven very difficult to detect due to its low density. Close to galaxy clusters, the filament gas is affected by the cluster's gravitational potential and attains substantial infall velocities, eventu… ▽ More A substantial fraction of the cosmic baryons is expected to hide in the form of diffuse warm-hot intergalactic medium (WHIM), the majority of which resides in the filaments of the Cosmic Web and has proven very difficult to detect due to its low density. Close to galaxy clusters, the filament gas is affected by the cluster's gravitational potential and attains substantial infall velocities, eventually undergoing a termination shock which may boost its X-ray signal. We aim to identify the optimal locations of the enhanced X-ray emission and absorption arising from cluster-filament interactions, as well as improve our understanding of the various physical processes affecting the WHIM as it approaches the cluster. We applied the DisPerSE filament finder to the galaxy distribution in the surroundings of a Coma-like ($M_{200} \sim 10^{15.4}~{\rm M}_{\odot}$) simulated C-EAGLE galaxy cluster. We characterised the thermodynamic properties of the gas in such filaments as well as their dependence on distance from the cluster, and provided a physical interpretation for the results. The identified filaments account for $\sim 50$% of the hot WHIM ($T > 10^{5.5}$ K) in the cluster vicinity. The filament gas is in approximate free-fall all the way down to $\sim 2 \ r_{200}$ from the cluster, at which stage it begins to slow down due to the increasing pressure of the ambient gas. The deceleration is accompanied by the conversion of gas bulk kinetic energy into heat and the increase of density and temperature from the general Cosmic Web level of $ρ\sim 10ρ_{\rm av}$ and $T = 10^5-10^6$ K towards $ρ\sim 100ρ_{\rm av}$ and $T = 10^7-10^8$ K near the cluster boundary. We conclude that the detection of the cosmic filaments of galaxies around clusters may provide a practical observational avenue for locating the densest and hottest phase of the missing baryons. △ Less

Submitted 6 March, 2023; originally announced March 2023.

Comments: 20 pages, 16 figures, accepted for publication in A&A

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

The multi-wavelength view of shocks in the fastest nova V1674 Her

Authors: K. V. Sokolovsky, T. J. Johnson, S. Buson, P. Jean, C. C. Cheung, K. Mukai, L. Chomiuk, E. Aydi, B. Molina, A. Kawash, J. D. Linford, A. J. Mioduszewski, M. P. Rupen, J. L. Sokoloski, M. N. Williams, E. Steinberg, I. Vurm, B. D. Metzger, K. L. Page, M. Orio, R. M. Quimby, A. W. Shafter, H. Corbett, S. Bolzoni, J. DeYoung , et al. (19 additional authors not shown)

Abstract: Classical novae are shock-powered multi-wavelength transients triggered by a thermonuclear runaway on an accreting white dwarf. V1674 Her is the fastest nova ever recorded (time to declined by two magnitudes is t_2=1.1 d) that challenges our understanding of shock formation in novae. We investigate the physical mechanisms behind nova emission from GeV gamma-rays to cm-band radio using coordinated… ▽ More Classical novae are shock-powered multi-wavelength transients triggered by a thermonuclear runaway on an accreting white dwarf. V1674 Her is the fastest nova ever recorded (time to declined by two magnitudes is t_2=1.1 d) that challenges our understanding of shock formation in novae. We investigate the physical mechanisms behind nova emission from GeV gamma-rays to cm-band radio using coordinated Fermi-LAT, NuSTAR, Swift and VLA observations supported by optical photometry. Fermi-LAT detected short-lived (18 h) 0.1-100 GeV emission from V1674 Her that appeared 6 h after the eruption began; this was at a level of (1.6 +/- 0.4)x10^-6 photons cm^-2 s^-1. Eleven days later, simultaneous NuSTAR and Swift X-ray observations revealed optically thin thermal plasma shock-heated to kT_shock = 4 keV. The lack of a detectable 6.7 keV Fe K_alpha emission suggests super-solar CNO abundances. The radio emission from V1674 Her was consistent with thermal emission at early times and synchrotron at late times. The radio spectrum steeply rising with frequency may be a result of either free-free absorption of synchrotron and thermal emission by unshocked outer regions of the nova shell or the Razin-Tsytovich effect attenuating synchrotron emission in dense plasma. The development of the shock inside the ejecta is unaffected by the extraordinarily rapid evolution and the intermediate polar host of this nova. △ Less

Submitted 21 March, 2023; v1 submitted 6 February, 2023; originally announced February 2023.

Comments: 20 pages, 9 figures, 3 tables. Accepted to MNRAS

Evidence for multiple shocks from the $γ$-ray emission of RS Ophiuchi

Authors: Rebecca Diesing, Brian D. Metzger, Elias Aydi, Laura Chomiuk, Indrek Vurm, Siddhartha Gupta, Damiano Caprioli

Abstract: In August of 2021, Fermi-LAT, H.E.S.S., and MAGIC detected GeV and TeV $γ$-ray emission from an outburst of recurrent nova RS Ophiuchi. This detection represents the first very high energy $γ$-rays observed from a nova, and opens a new window to study particle acceleration. Both H.E.S.S. and MAGIC described the observed $γ$-rays as arising from a single, external shock. In this paper, we perform d… ▽ More In August of 2021, Fermi-LAT, H.E.S.S., and MAGIC detected GeV and TeV $γ$-ray emission from an outburst of recurrent nova RS Ophiuchi. This detection represents the first very high energy $γ$-rays observed from a nova, and opens a new window to study particle acceleration. Both H.E.S.S. and MAGIC described the observed $γ$-rays as arising from a single, external shock. In this paper, we perform detailed, multi-zone modeling of RS Ophiuchi's 2021 outburst including a self-consistent prescription for particle acceleration and magnetic field amplification. We demonstrate that, contrary to previous work, a single shock cannot simultaneously explain RS Ophiuchi's GeV and TeV emission, particularly the spectral shape and distinct light curve peaks. Instead, we put forward a model involving multiple shocks that reproduces the observed $γ$-ray spectrum and temporal evolution. The simultaneous appearance of multiple distinct velocity components in the nova optical spectrum over the first several days of the outburst supports the presence of distinct shocks, which may arise either from the strong latitudinal dependence of the density of the external circumbinary medium (e.g., in the binary equatorial plane versus the poles) or due to internal collisions within the white dwarf ejecta (as powers the $γ$-ray emission in classical novae). △ Less

Submitted 9 May, 2023; v1 submitted 3 November, 2022; originally announced November 2022.

Comments: 19 pages, 10 figures, published in ApJ

$γ$-ray Emission from Classical Nova V392 Per: Measurements from Fermi and HAWC

Authors: A. Albert, R. Alfaro, C. Alvarez, J. C. Arteaga-Velázquez, D. Avila Rojas, H. A. Ayala Solares, R. Babu, E. Belmont-Moreno, C. Blochwitz, K. S. Caballero-Mora, T. Capistrán, A. Carramiñana, S. Casanova, O. Chaparro-Amaro, U. Cotti, J. Cotzomi, E. De la Fuente. C. de León. S. Coutiño de León, R. Diaz Hernandez, B. L. Dingus, M. A. DuVernois, M. Durocher, J. C. Díaz-Vélez, K. Engel, C. Espinoza, K. L. Fan , et al. (62 additional authors not shown)

Abstract: This paper reports on the $γ$-ray properties of the 2018 Galactic nova V392 Per, spanning photon energies $\sim$0.1 GeV to 100 TeV by combining observations from the Fermi Gamma-ray Space Telescope and the HAWC Observatory. In one of the most rapidly evolving $γ$-ray signals yet observed for a nova, GeV $γ$ rays with a power law spectrum with index $Γ= 2.0 \pm 0.1$ were detected over eight days fo… ▽ More This paper reports on the $γ$-ray properties of the 2018 Galactic nova V392 Per, spanning photon energies $\sim$0.1 GeV to 100 TeV by combining observations from the Fermi Gamma-ray Space Telescope and the HAWC Observatory. In one of the most rapidly evolving $γ$-ray signals yet observed for a nova, GeV $γ$ rays with a power law spectrum with index $Γ= 2.0 \pm 0.1$ were detected over eight days following V392 Per's optical maximum. HAWC observations constrain the TeV $γ$-ray signal during this time and also before and after. We observe no statistically significant evidence of TeV $γ$-ray emission from V392 Per, but present flux limits. Tests of the extension of the Fermi/LAT spectrum to energies above 5 TeV are disfavored by 2 standard deviations (95\%) or more. We fit V392 Per's GeV $γ$ rays with hadronic acceleration models, incorporating optical observations, and compare the calculations with HAWC limits. △ Less

Submitted 16 December, 2022; v1 submitted 25 January, 2022; originally announced January 2022.

Comments: 18 pages, 14 figures; revised for referee comments; section 5 revised

Journal ref: Astroph. J. 940(2) (2022) 141

The first nova eruption in a novalike variable: YZ Ret as seen in X-rays and gamma-rays

Authors: Kirill V. Sokolovsky, Kwan-Lok Li, Raimundo Lopes de Oliveira, Jan-Uwe Ness, Koji Mukai, Laura Chomiuk, Elias Aydi, Elad Steinberg, Indrek Vurm, Brian D. Metzger, Aliya-Nur Babul, Adam Kawash, Justin D. Linford, Thomas Nelson, Kim L. Page, Michael P. Rupen, Jennifer L. Sokoloski, Jay Strader, David Kilkenny

Abstract: Peaking at 3.7 mag on 2020 July 11, YZ Ret was the second-brightest nova of the decade. The nova's moderate proximity (2.7 kpc from Gaia) provided an opportunity to explore its multi-wavelength properties in great detail. Here we report on YZ Ret as part of a long-term project to identify the physical mechanisms responsible for high-energy emission in classical novae. We use simultaneous Fermi/LAT… ▽ More Peaking at 3.7 mag on 2020 July 11, YZ Ret was the second-brightest nova of the decade. The nova's moderate proximity (2.7 kpc from Gaia) provided an opportunity to explore its multi-wavelength properties in great detail. Here we report on YZ Ret as part of a long-term project to identify the physical mechanisms responsible for high-energy emission in classical novae. We use simultaneous Fermi/LAT and NuSTAR observations complemented by XMM-Newton X-ray grating spectroscopy to probe the physical parameters of the shocked ejecta and the nova-hosting white dwarf. The XMM-Newton observations revealed a super-soft X-ray emission which is dominated by emission lines of CV, CVI, NVI, NVII, and OVIII rather than a blackbody-like continuum, suggesting CO-composition of the white dwarf in a high-inclination binary system. Fermi/LAT detected YZ Ret for 15 days with the gamma-ray spectrum best described by a power law with an exponential cut-off at 1.9 +/-0.6 GeV. In stark contrast with theoretical predictions and in keeping with previous NuSTAR observations of Fermi-detected classical novae (V5855 Sgr and V906 Car), the 3.5-78 keV X-ray emission is found to be two orders of magnitude fainter than the GeV emission. The X-ray emission observed by NuSTAR is consistent with a single-temperature thermal plasma. We detect no non-thermal tail of the GeV emission expected to extend down to the NuSTAR band. NuSTAR observations continue to challenge theories of high-energy emission from shocks in novae. △ Less

Submitted 21 May, 2022; v1 submitted 6 August, 2021; originally announced August 2021.

Comments: 22 pages, 5 figures, 4 tables; accepted to MNRAS

Gamma-ray Thermalization and Leakage from Millisecond Magnetar Nebulae: Towards a Self-Consistent Model for Superluminous Supernovae

Authors: Indrek Vurm, Brian D. Metzger

Abstract: Superluminous supernovae (SLSNe) are massive star explosions too luminous to be powered by traditional energy sources, such as radioactive 56Ni. These transients may instead be powered by a central engine, such as a millisecond pulsar or magnetar, whose relativistic wind inflates a nebula of high energy particles and radiation behind the expanding ejecta. We present 3D Monte Carlo radiative transf… ▽ More Superluminous supernovae (SLSNe) are massive star explosions too luminous to be powered by traditional energy sources, such as radioactive 56Ni. These transients may instead be powered by a central engine, such as a millisecond pulsar or magnetar, whose relativistic wind inflates a nebula of high energy particles and radiation behind the expanding ejecta. We present 3D Monte Carlo radiative transfer calculations which follow the production and thermalization of high energy radiation from the nebula into optical radiation and, conversely, determine the gamma-ray emission that escapes the ejecta without thermalizing. We track the evolution of photons and matter in a coupled two-zone ("wind/nebula" and "ejecta") model, accounting for the range of radiative processes. We identify a novel mechanism by which gamma-gamma pair creation in the upstream pulsar wind regulates the mean energy of particles entering the nebula over the first several years after the explosion, rendering our results on this timescale insensitive to the (uncertain) intrinsic wind pair multiplicity. To explain the observed late-time steepening of SLSNe optical light curves as being the result of gamma-ray leakage, the nebular magnetization must be very low, epsB <~ 1e-6-1e-4. For higher epsB, synchrotron emission quickly comes to dominate the thermalized nebula radiation, and being readily absorbed because of its lower photon energies, results in the SN optical light curve tracking the spin-down power even to late times >~ 1 yr, inconsistent with observations. For magnetars to remain viable contenders for powering SLSNe, we conclude that either magnetic dissipation in the wind/nebula is extremely efficient, or that the spin-down luminosity decays significantly faster than the canonical dipole rate ~1/t^2 in a way that coincidentally mimicks gamma-ray escape. △ Less

Submitted 13 January, 2021; originally announced January 2021.

Comments: 25 pages plus appendices, 13 figures, submitted to ApJ

High-Energy Neutrinos and Gamma-Rays from Non-Relativistic Shock-Powered Transients

Authors: Ke Fang, Brian D. Metzger, Indrek Vurm, Elias Aydi, Laura Chomiuk

Abstract: Shock interaction has been argued to play a role in powering a range of optical transients, including supernovae (particularly the superluminous class), classical novae, stellar mergers, tidal disruption events, and fast blue optical transients. These same shocks can accelerate relativistic ions, generating high-energy neutrino and gamma-ray emission via hadronic pion production. The recent discov… ▽ More Shock interaction has been argued to play a role in powering a range of optical transients, including supernovae (particularly the superluminous class), classical novae, stellar mergers, tidal disruption events, and fast blue optical transients. These same shocks can accelerate relativistic ions, generating high-energy neutrino and gamma-ray emission via hadronic pion production. The recent discovery of time-correlated optical and gamma-ray emission in classical novae has revealed the important role of radiative shocks in powering these events, enabling an unprecedented view of the properties of ion acceleration, including its efficiency and energy spectrum, under similar physical conditions to shocks in extragalactic transients. Here we introduce a model for connecting the radiated optical fluence of non-relativistic transients to their maximal neutrino and gamma-ray fluence. We apply this technique to a wide range of extragalactic transient classes in order to place limits on their contributions to the cosmological high-energy gamma-ray and neutrino backgrounds. Based on a simple model for diffusive shock acceleration at radiative shocks, calibrated to novae, we demonstrate that several of the most luminous transients can accelerate protons up to energies $E_{\rm max} \gtrsim 10^{16}$ eV, sufficient to contribute to the IceCube astrophysical background. Furthermore, several of the considered sources$-$particularly hydrogen-poor supernovae$-$may serve as "hidden" gamma-ray sources due to the high gamma-ray opacity of their ejecta, evading constraints imposed by the non-blazar Fermi-LAT background. However, adopting an ion acceleration efficiency $\sim$ 0.3-1$\%$ motivated by nova observations, we find that currently known classes of non-relativistic, potentially shock-powered transients contribute at most a few percent of the total IceCube background. △ Less

Submitted 1 January, 2021; v1 submitted 30 July, 2020; originally announced July 2020.

Comments: 24 pages, 7 figures, 2 tables

Journal ref: ApJ 904, 4, 2020

X-ray Spectroscopy of the gamma-ray Brightest Nova V906 Car (ASASSN-18fv)

Authors: Kirill V. Sokolovsky, Koji Mukai, Laura Chomiuk, Raimundo Lopes de Oliveira, Elias Aydi, Kwan-Lok Li, Elad Steinberg, Indrek Vurm, Brian D. Metzger, Adam Kawash, Justin D. Linford, Amy J. Mioduszewski, Thomas Nelson, Jan-Uwe Ness, Kim L. Page, Michael P. Rupen, Jennifer L. Sokoloski, Jay Strader

Abstract: Shocks in gamma-ray emitting classical novae are expected to produce bright thermal and non-thermal X-rays. We test this prediction with simultaneous NuSTAR and Fermi/LAT observations of nova V906 Car, which exhibited the brightest GeV gamma-ray emission to date. The nova is detected in hard X-rays while it is still gamma-ray bright, but contrary to simple theoretical expectations, the detected 3.… ▽ More Shocks in gamma-ray emitting classical novae are expected to produce bright thermal and non-thermal X-rays. We test this prediction with simultaneous NuSTAR and Fermi/LAT observations of nova V906 Car, which exhibited the brightest GeV gamma-ray emission to date. The nova is detected in hard X-rays while it is still gamma-ray bright, but contrary to simple theoretical expectations, the detected 3.5-78 keV emission of V906 Car is much weaker than the simultaneously observed >100 MeV emission. No non-thermal X-ray emission is detected, and our deep limits imply that the gamma-rays are likely hadronic. After correcting for substantial absorption (N_H ~ 2 x 10^23 cm^-2), the thermal X-ray luminosity (from a 9 keV optically-thin plasma) is just ~2% of the gamma-ray luminosity. We consider possible explanations for the low thermal X-ray luminosity, including the X-rays being suppressed by corrugated, radiative shock fronts or the X-rays from the gamma-ray producing shock are hidden behind an even larger absorbing column (N_H >10^25 cm^-2). Adding XMM-Newton and Swift/XRT observations to our analysis, we find that the evolution of the intrinsic X-ray absorption requires the nova shell to be expelled 24 days after the outburst onset. The X-ray spectra show that the ejecta are enhanced in nitrogen and oxygen, and the nova occurred on the surface of a CO-type white dwarf. We see no indication of a distinct super-soft phase in the X-ray lightcurve, which, after considering the absorption effects, may point to a low mass of the white dwarf hosting the nova. △ Less

Submitted 15 July, 2020; originally announced July 2020.

Comments: 19 pages, 6 figures, 7 tables; accepted to MNRAS

Direct evidence for shock-powered optical emission in a nova

Authors: Elias Aydi, Kirill V. Sokolovsky, Laura Chomiuk, Elad Steinberg, Kwan Lok Li, Indrek Vurm, Brian D. Metzger, Jay Strader, Koji Mukai, Ondřej Pejcha, Ken J. Shen, Gregg A. Wade, Rainer Kuschnig, Anthony F. J. Moffat, Herbert Pablo, Andrzej Pigulski, Adam Popowicz, Werner Weiss, Konstanze Zwintz, Luca Izzo, Karen R. Pollard, Gerald Handler, Stuart D. Ryder, Miroslav D. Filipović, Rami Z. E. Alsaberi , et al. (17 additional authors not shown)

Abstract: Classical novae are thermonuclear explosions that occur on the surfaces of white dwarf stars in interacting binary systems (Bode & Evans 2008). It has long been thought that the luminosity of classical novae is powered by continued nuclear burning on the surface of the white dwarf after the initial runaway (Gallaher & Starrfield 1978). However, recent observations of GeV $γ$-rays from classical no… ▽ More Classical novae are thermonuclear explosions that occur on the surfaces of white dwarf stars in interacting binary systems (Bode & Evans 2008). It has long been thought that the luminosity of classical novae is powered by continued nuclear burning on the surface of the white dwarf after the initial runaway (Gallaher & Starrfield 1978). However, recent observations of GeV $γ$-rays from classical novae have hinted that shocks internal to the nova ejecta may dominate the nova emission. Shocks have also been suggested to power the luminosity of events as diverse as stellar mergers (Metzger & Pejcha 2017), supernovae (Moriya et al. 2018), and tidal disruption events (Roth et al. 2016), but observational confirmation has been lacking. Here we report simultaneous space-based optical and $γ$-ray observations of the 2018 nova V906 Carinae (ASASSN-18fv), revealing a remarkable series of distinct correlated flares in both bands. The optical and $γ$-ray flares occur simultaneously, implying a common origin in shocks. During the flares, the nova luminosity doubles, implying that the bulk of the luminosity is shock-powered. Furthermore, we detect concurrent but weak X-ray emission from deeply embedded shocks, confirming that the shock power does not appear in the X-ray band and supporting its emergence at longer wavelengths. Our data, spanning the spectrum from radio to $γ$-ray, provide direct evidence that shocks can power substantial luminosity in classical novae and other optical transients. △ Less

Submitted 12 April, 2020; originally announced April 2020.

Comments: 49 pages including supplementary information. Accepted and published in Nature Astronomy

Journal ref: Nature Astronomy 2020

Radio Time-Domain Signatures of Magnetar Birth

Authors: Casey J. Law, Ben Margalit, Nipuni T. Palliyaguru, Brian D. Metzger, Lorenzo Sironi, Yong Zheng, Edo Berger, Raffaella Margutti, Andrei Beloborodov, Matt Nicholl, Tarraneh Eftekhari, Indrek Vurm, Peter K. G. Williams

Abstract: The last decade has seen the rapid, concurrent development of new classes of energetic astrophysical transients, including Fast Radio Bursts, Superluminous Supernovae, and ultralong gamma-ray bursts. We consider the implications of an emerging consensus on the origin of these diverse phenomena as different signatures of magnetar birth. We describe how magnetar-powered radio transients from millise… ▽ More The last decade has seen the rapid, concurrent development of new classes of energetic astrophysical transients, including Fast Radio Bursts, Superluminous Supernovae, and ultralong gamma-ray bursts. We consider the implications of an emerging consensus on the origin of these diverse phenomena as different signatures of magnetar birth. We describe how magnetar-powered radio transients from milliseconds to gigaseconds will be used to address diverse physical questions related to the stellar evolution and death, the circum- and intergalactic medium, and fundamental physics on scales ranging from the quantum to cosmological. We argue that high-sensitivity radio interferometer concepts in development will be critical to accessing this science. There is also a growing need for scientists with expertise in new statistical tools and software techniques. △ Less

Submitted 11 March, 2019; originally announced March 2019.

Comments: Science white paper for Astro2020

NuSTAR Detection of X-rays Concurrent with Gamma Rays in the Nova V5855 Sgr

Authors: Thomas Nelson, Koji Mukai, Kwan-Lok Li, Indrek Vurm, Brian D. Metzger, Laura Chomiuk, J. L. Sokoloski, Justin D. Linford, Terry Bohlsen, Paul Luckas

Abstract: We report the first detection of hard (>10 keV) X-ray emission simultaneous with gamma rays in a nova eruption. Observations of the nova V5855 Sgr carried out with the NuSTAR satellite on Day 12 of the eruption revealed faint, highly absorbed thermal X-rays. The extreme equivalent hydrogen column density towards the X-ray emitting region (~3 x 10$^{24}$ cm$^{-2}$) indicates that the shock producin… ▽ More We report the first detection of hard (>10 keV) X-ray emission simultaneous with gamma rays in a nova eruption. Observations of the nova V5855 Sgr carried out with the NuSTAR satellite on Day 12 of the eruption revealed faint, highly absorbed thermal X-rays. The extreme equivalent hydrogen column density towards the X-ray emitting region (~3 x 10$^{24}$ cm$^{-2}$) indicates that the shock producing the X-rays was deeply embedded within the nova ejecta. The slope of the X-ray spectrum favors a thermal origin for the bulk of the emission, and the constraints of the temperature in the shocked region suggest a shock velocity compatible with the ejecta velocities inferred from optical spectroscopy. While we do not claim the detection of non-thermal X-rays, the data do not allow us to rule out an additional, fainter component dominating at energy above 20 keV, for which we obtained upper limits. The inferred luminosity of the thermal X-rays is too low to be consistent with the gamma-ray luminosities if both are powered by the same shock under standard assumptions regarding the efficiency of non-thermal particle acceleration and the temperature distribution of the shocked gas. △ Less

Submitted 31 December, 2018; originally announced January 2019.

Comments: Accepted for publication in the Astrophysical Journal - 10 pages, 5 figures

An embedded X-ray source shines through the aspherical AT2018cow: revealing the inner workings of the most luminous fast-evolving optical transients

Authors: Raffaella Margutti, B. D. Metzger, R. Chornock, I. Vurm, N. Roth, B. W. Grefenstette, V. Savchenko, R. Cartier, J. F. Steiner, G. Terreran, G. Migliori, D. Milisavljevic, K. D. Alexander, M. Bietenholz, P. K. Blanchard, E. Bozzo, D. Brethauer, I. V. Chilingarian, D. L. Coppejans, L. Ducci, C. Ferrigno, W. Fong, D. GÖtz, C. Guidorzi, A. Hajela , et al. (19 additional authors not shown)

Abstract: We present the first extensive radio to gamma-ray observations of a fast-rising blue optical transient (FBOT), AT2018cow, over its first ~100 days. AT2018cow rose over a few days to a peak luminosity $L_{pk}\sim4\times 10^{44}$ erg/s exceeding those of superluminous supernovae (SNe), before declining as $\propto t^{-2}$. Initial spectra at $\lesssim 15$ days were mostly featureless and indicated l… ▽ More We present the first extensive radio to gamma-ray observations of a fast-rising blue optical transient (FBOT), AT2018cow, over its first ~100 days. AT2018cow rose over a few days to a peak luminosity $L_{pk}\sim4\times 10^{44}$ erg/s exceeding those of superluminous supernovae (SNe), before declining as $\propto t^{-2}$. Initial spectra at $\lesssim 15$ days were mostly featureless and indicated large expansion velocities v~0.1c and temperatures reaching 30000 K. Later spectra revealed a persistent optically-thick photosphere and the emergence of H and He emission features with v~sim 4000 km/s with no evidence for ejecta cooling. Our broad-band monitoring revealed a hard X-ray spectral component at $E\ge 10$ keV, in addition to luminous and highly variable soft X-rays, with properties unprecedented among astronomical transients. An abrupt change in the X-ray decay rate and variability appears to accompany the change in optical spectral properties. AT2018cow showed bright radio emission consistent with the interaction of a blastwave with $v_{sh}$~0.1c with a dense environment ($\dot M\sim10^{-3}-10^{-4}\,M_{\odot}yr^{-1}$ for $v_w=1000$ km\s). While these properties exclude Ni-powered transients, our multi-wavelength analysis instead indicates that AT2018cow harbored a "central engine", either a compact object (magnetar or black hole) or an embedded internal shock produced by interaction with a compact, dense circumstellar medium. The engine released $\sim10^{50}-10^{51.5}$ erg over $\sim10^3-10^5$ s and resides within low-mass fast-moving material with equatorial-polar density asymmetry ($M_{ej,fast}\lesssim0.3\,\rm{M_{\odot}}$). Successful SNe from low-mass H-rich stars (like electron-capture SNe) or failed explosions from blue supergiants satisfy these constraints. Intermediate-mass black-holes are disfavored by the large environmental density probed by the radio observations. △ Less

Submitted 25 October, 2018; originally announced October 2018.

Comments: 26 pages + tables. Submitted

A strong limit on the very-high-energy emission from GRB 150323A

Authors: A. U. Abeysekara, A. Archer, W. Benbow, R. Bird, R. Brose, M. Buchovecky, V. Bugaev, M. P. Connolly, W. Cui, M. Errando, A. Falcone, Q. Feng, J. P. Finley, A. Flinders, L. Fortson, A. Furniss, G. H. Gillanders, M. Hütten, D. Hanna, O. Hervet, J. Holder, G. Hughes, T. B. Humensky, C. A. Johnson, P. Kaaret , et al. (41 additional authors not shown)

Abstract: On 2015 March 23, VERITAS responded to a $Swift$-BAT detection of a gamma-ray burst, with observations beginning 270 seconds after the onset of BAT emission, and only 135 seconds after the main BAT emission peak. No statistically significant signal is detected above 140 GeV. The VERITAS upper limit on the fluence in a 40 minute integration corresponds to about 1% of the prompt fluence. Our limit i… ▽ More On 2015 March 23, VERITAS responded to a $Swift$-BAT detection of a gamma-ray burst, with observations beginning 270 seconds after the onset of BAT emission, and only 135 seconds after the main BAT emission peak. No statistically significant signal is detected above 140 GeV. The VERITAS upper limit on the fluence in a 40 minute integration corresponds to about 1% of the prompt fluence. Our limit is particularly significant since the very-high-energy (VHE) observation started only $\sim$2 minutes after the prompt emission peaked, and $Fermi$-LAT observations of numerous other bursts have revealed that the high-energy emission is typically delayed relative to the prompt radiation and lasts significantly longer. Also, the proximity of GRB~150323A ($z=0.593$) limits the attenuation by the extragalactic background light to $\sim 50$ % at 100-200 GeV. We conclude that GRB 150323A had an intrinsically very weak high-energy afterglow, or that the GeV spectrum had a turnover below $\sim100$ GeV. If the GRB exploded into the stellar wind of a massive progenitor, the VHE non-detection constrains the wind density parameter to be $A\gtrsim 3\times 10^{11}$ g cm$^{-1}$, consistent with a standard Wolf-Rayet progenitor. Alternatively, the VHE emission from the blast wave would be weak in a very tenuous medium such as the ISM, which therefore cannot be ruled out as the environment of GRB 150323A. △ Less

Submitted 3 March, 2018; originally announced March 2018.

Comments: 6 pages, 1 figure. Accepted for publication in The Astrophysical Journal

A Nova Outburst Powered by Shocks

Authors: Kwan-Lok Li, Brian D. Metzger, Laura Chomiuk, Indrek Vurm, Jay Strader, Thomas Finzell, Andrei M. Beloborodov, Thomas Nelson, Benjamin J. Shappee, Christopher S. Kochanek, Jose L. Prieto, Stella Kafka, Thomas W. -S. Holoien, Todd A. Thompson, Paul J. Luckas, Hiroshi Itoh

Abstract: Classical novae are runaway thermonuclear burning events on the surfaces of accreting white dwarfs in close binary star systems, sometimes appearing as new naked-eye sources in the night sky. The standard model of novae predicts that their optical luminosity derives from energy released near the hot white dwarf which is reprocessed through the ejected material. Recent studies with the Fermi Large… ▽ More Classical novae are runaway thermonuclear burning events on the surfaces of accreting white dwarfs in close binary star systems, sometimes appearing as new naked-eye sources in the night sky. The standard model of novae predicts that their optical luminosity derives from energy released near the hot white dwarf which is reprocessed through the ejected material. Recent studies with the Fermi Large Area Telescope have shown that many classical novae are accompanied by gigaelectronvolt gamma-ray emission. This emission likely originates from strong shocks, providing new insights into the properties of nova outflows and allowing them to be used as laboratories to study the unknown efficiency of particle acceleration in shocks. Here we report gamma-ray and optical observations of the Milky Way nova ASASSN-16ma, which is among the brightest novae ever detected in gamma-rays. The gamma-ray and optical light curves show a remarkable correlation, implying that the majority of the optical light comes from reprocessed emission from shocks rather than the white dwarf. The ratio of gamma-ray to optical flux in ASASSN-16ma directly constrains the acceleration efficiency of non-thermal particles to be ~0.005, favouring hadronic models for the gamma-ray emission. The need to accelerate particles up to energies exceeding 100 gigaelectronvolts provides compelling evidence for magnetic field amplification in the shocks. △ Less

Submitted 3 September, 2017; originally announced September 2017.

Comments: Published in Nature Astronomy. This is the authors' version with 55 pages, 8 figures, and 3 tables

Journal ref: Nature Astronomy, 1, 0222 (2017)

Radiation mediated shocks in gamma-ray bursts: Pair creation

Authors: Christoffer Lundman, Andrei Beloborodov, Indrek Vurm

Abstract: Sub-photospheric shock dissipation is one of the main proposed mechanisms for producing the prompt gamma-ray burst (GRB) emission. Such shocks are mediated by scattering of radiation. We introduce a time dependent, special relativistic code which dynamically couples Monte Carlo radiative transfer to the flow hydrodynamics. The code also self-consistently implements electron-positron pair productio… ▽ More Sub-photospheric shock dissipation is one of the main proposed mechanisms for producing the prompt gamma-ray burst (GRB) emission. Such shocks are mediated by scattering of radiation. We introduce a time dependent, special relativistic code which dynamically couples Monte Carlo radiative transfer to the flow hydrodynamics. The code also self-consistently implements electron-positron pair production and annihilation. We simulate shocks with properties relevant for GRBs and study the steady-state solutions, which are accurate deep below the jet photosphere. The shock generates a power-law photon spectrum through the first-order Fermi mechanism, extending upwards from the typical upstream photon energy. Strong shocks (for which the downstream pressure is much larger than the upstream pressure) have rising $νF_ν$ shock spectra. The spectrum extends up to $ε_{max} \equiv E_{max}/m_e c^2 \sim v^2$ for non-relativistic shocks, where $m_e$ is the electron rest mass and $v$ is the relative speed between the upstream and downstream in units of the speed of light $c$. For mildly relativistic shocks the power law softens at $ε\gtrsim 10^{-1}$ due to Klein-Nishina effects, and shocks with $vγ\gtrsim 1$, where $γ\equiv (1-v^2)^{-1/2}$, produce electron-positron pairs. As an example, a strong shock with $vγ= 3$ and a photon-to-proton ratio of $n_γ/n_p = 2 \times 10^5$ has a peak pair-to-proton ratio of $Z_\pm \approx 225$. The main effect of pairs in a steady-state shock is to decrease its spatial width by a factor of $\sim Z_\pm$. The post-shock spectrum thermalizes in the downstream. In absence of emission and absorption processes, kinetic equilibrium at temperature $θ_d \equiv kT_d/m_e c^2 \approx ε_d/3$ is reached at an optical depth of $τ\gg θ_d^{-1}$ behind the shock, where $ε_d$ is the average downstream photon energy. △ Less

Submitted 8 August, 2017; originally announced August 2017.

Comments: 14 pages, 8 figures, submitted to ApJ

On the prospects of gamma-ray burst detection in the TeV band

Authors: Indrek Vurm, Andrei M. Beloborodov

Abstract: A gamma-ray burst (GRB) jet running into an external medium is expected to generate luminous GeV-TeV emission lasting from minutes to several hours. The high-energy emission results from inverse Compton upscattering of prompt and afterglow photons by shock-heated {\it thermal} plasma. At its peak the high-energy radiation carries a significant fraction of the power dissipated at the forward shock.… ▽ More A gamma-ray burst (GRB) jet running into an external medium is expected to generate luminous GeV-TeV emission lasting from minutes to several hours. The high-energy emission results from inverse Compton upscattering of prompt and afterglow photons by shock-heated {\it thermal} plasma. At its peak the high-energy radiation carries a significant fraction of the power dissipated at the forward shock. We discuss in detail the expected TeV luminosity, using a robust "minimal" emission model. Then, using the statistical properties of the GRB population (luminosity function, redshift distribution, afterglow energy) we simulate the expected detection rates of GRBs by current and upcoming atmospheric Cherenkov instruments. We find that GRBs exploding into a low-density interstellar medium must produce TeV emission that would have already been detected by the currently operating Cherenkov telescopes. The absence of detections is consistent with explosions into a dense wind of the GRB progenitor. If the typical environment of long GRBs is a Wolf-Rayet progenitor wind, as suggested by the recent analysis of Fermi LAT data, then the upcoming Cherenkov Telescope Array (CTA) should be able to detect a fraction of GRBs that trigger the space-borne detectors. Since absorption by the extragalactic background light limits the detectability above 0.1 TeV for all but the most nearby bursts ($z\lesssim 1$), the reduced energy threshold of CTA is the key improvement over current instruments, which should boost the detection rate by at least a factor of 3. △ Less

Submitted 15 November, 2016; originally announced November 2016.

Comments: 15 pages, 10 figures, submitted to ApJ

High-energy emission from non-relativistic radiative shocks: application to gamma-ray novae

Authors: Indrek Vurm, Brian D. Metzger

Abstract: Multiwavelength radiation from relativistic particles accelerated at shocks in novae and other astrophysical sources carries a wealth of information about the outflow properties and the microphysical processes at work near the shocks. The observation of GeV gamma-rays from novae by Fermi/LAT demonstrates that the shocks in these systems can accelerate particles to energies of at least $\sim 10$ Ge… ▽ More Multiwavelength radiation from relativistic particles accelerated at shocks in novae and other astrophysical sources carries a wealth of information about the outflow properties and the microphysical processes at work near the shocks. The observation of GeV gamma-rays from novae by Fermi/LAT demonstrates that the shocks in these systems can accelerate particles to energies of at least $\sim 10$ GeV. The low-energy extension of the same non-thermal particle distribution inevitably gives rise to emission extending into the X-ray band. Above $\gtrsim 10$ keV this radiation can escape the system without significant absorption/attenuation, and can potentially be detected by NuSTAR. We present theoretical models for hard X-ray and gamma-ray emission from radiative shocks in both leptonic and hadronic scenarios, accounting for the rapid evolution of the downstream properties due to the fast cooling of thermal plasma. Due to strong Coulomb cooling of the mildly relativistic electrons nominally responsible for producing hard X-ray emission, only a fraction of $10^{-4} - 10^{-3}$ of the gamma-ray luminosity is radiated in the NuSTAR band; nevertheless, this emission could be detectable simultaneous with the LAT emission in bright gamma-ray novae with a $\sim 50$ ks exposure. Our work demonstrates how combined hard X-ray and gamma-ray observations can be used to constrain properties of the nova outflow (velocity, density and mass outflow rate) and particle acceleration at the shock. A very low X-ray to gamma-ray luminosity ratio ($L_{\rm X}/L_γ \lesssim 5 \times 10^{-4}$) would disfavor leptonic models for the gamma-ray emission. Our model can also be applied to other astrophysical environments with radiative shocks, including Type IIn supernovae and colliding winds in massive star binaries. △ Less

Submitted 14 November, 2016; originally announced November 2016.

Comments: 17 pages, 12 figures, submitted to ApJ

Polarization of gamma-ray bursts in the dissipative photosphere model

Authors: Christoffer Lundman, Indrek Vurm, Andrei M. Beloborodov

Abstract: The MeV spectral peak of gamma-ray bursts (GRBs) is best explained as photospheric emission from a dissipative relativistic jet. The observed non-blackbody spectrum shows that sub-photospheric dissipation involves both thermal plasma heating and injection of nonthermal particles, which quickly cool through inverse Compton scattering and emission of synchrotron radiation. Synchrotron photons emitte… ▽ More The MeV spectral peak of gamma-ray bursts (GRBs) is best explained as photospheric emission from a dissipative relativistic jet. The observed non-blackbody spectrum shows that sub-photospheric dissipation involves both thermal plasma heating and injection of nonthermal particles, which quickly cool through inverse Compton scattering and emission of synchrotron radiation. Synchrotron photons emitted around and above the photosphere are predicted to dominate the low-energy part of the GRB spectrum, starting from roughly a decade in energy below the MeV peak. We show that this leads to a unique polarization signature: a rise in GRB polarization toward lower energies. We compute the polarization degree of GRB radiation as a function of photon energy for a generic jet model, and show the predictions for GRBs 990123, 090902B and 110721A. The expected polarization is significant in the X-ray band, in particular for bursts similar to GRB 090902B. Radiation in the MeV peak (and at higher energies) is unpolarized as long as the jet is approximately uniform on angular scales $δθ\gtrsim Γ^{-1}$ where $Γ$ is the bulk Lorentz factor of the jet. △ Less

Submitted 4 November, 2016; originally announced November 2016.

Comments: 10 pages, 3, figures, submitted to ApJ

Shocks in nova outflows. II. Synchrotron radio emission

Authors: Andrey Vlasov, Indrek Vurm, Brian Metzger

Abstract: The discovery of GeV gamma-rays from classical novae indicates that shocks and relativistic particle acceleration are energetically key in these events. Further evidence for shocks comes from thermal keV X-ray emission and an early peak in the radio light curve on a timescale of months with a brightness temperature which is too high to result from freely expanding photo-ionized gas. Paper I develo… ▽ More The discovery of GeV gamma-rays from classical novae indicates that shocks and relativistic particle acceleration are energetically key in these events. Further evidence for shocks comes from thermal keV X-ray emission and an early peak in the radio light curve on a timescale of months with a brightness temperature which is too high to result from freely expanding photo-ionized gas. Paper I developed a one dimensional model for the thermal emission from nova shocks. This work concluded that the shock-powered radio peak cannot be thermal if line cooling operates in the post-shock gas at the rate determined by collisional ionization equilibrium. Here we extend this calculation to include non-thermal synchrotron emission. Applying our model to three classical novae, we constrain the amplification of the magnetic field $ε_B$ and the efficiency $ε_e$ of accelerating relativistic electrons of characteristic Lorentz factor $γ\sim 100$. If the shocks are radiative (low velocity $v_{\rm sh} \lesssim 1000$ km s$^{-1}$) and cover a large solid angle of the nova outflow, as likely characterize those producing gamma-rays, then values of $ε_e \sim 0.01-0.1$ are required to achieve the peak radio brightness for $ε_B = 10^{-2}$. Such high efficiencies exclude secondary pairs from pion decay as the source of the radio-emitting particles, instead favoring the direct acceleration of electrons at the shock. If the radio-emitting shocks are instead adiabatic (high velocity), as likely characterize those responsible for the thermal X-rays, then much higher brightness temperatures are possible, allowing the radio-emitting shocks to cover a smaller outflow solid angle. △ Less

Submitted 25 August, 2016; v1 submitted 16 March, 2016; originally announced March 2016.

Comments: v2 abridged to match MNRAS version

Novae as Tevatrons: Prospects for CTA and IceCube

Authors: Brian D. Metzger, Damiano Caprioli, Indrek Vurm, Andrei Beloborodov, Imre Bartos, Andrey Vlasov

Abstract: The discovery of novae as sources of ~GeV gamma-rays highlights the key role of shocks and relativistic particle acceleration in these transient systems. Although there is evidence for a spectral cut-off above energies ~1-100 GeV at particular epochs in some novae, the maximum particle energy achieved in these accelerators has remained an open question. The high densities of the nova ejecta (~10 o… ▽ More The discovery of novae as sources of ~GeV gamma-rays highlights the key role of shocks and relativistic particle acceleration in these transient systems. Although there is evidence for a spectral cut-off above energies ~1-100 GeV at particular epochs in some novae, the maximum particle energy achieved in these accelerators has remained an open question. The high densities of the nova ejecta (~10 orders of magnitude larger than in supernova remnants) render the gas far upstream of the shock neutral and shielded from ionizing radiation. The amplification of the magnetic field needed for diffusive shock acceleration requires ionized gas, thus confining the acceleration process to a narrow photo-ionized layer immediately ahead of the shock. Based on the growth rate in this layer of the hybrid non-resonant cosmic ray current-driven instability (considering also ion-neutral damping), we quantify the maximum particle energy, Emax, across the range of shock velocities and upstream densities of interest. We find values of Emax ~ 10 GeV - 10 TeV, which are broadly consistent with the inferred spectral cut-offs, but which could also in principle lead to emission extending to higher energies >100 GeV accessible to atmosphere Cherenkov telescopes, such as the planned Cherenkov Telescope Array (CTA). Detecting TeV neutrinos with IceCube in hadronic scenarios appears to be more challenging, although the prospects are improved if the shock power during the earliest, densest phases of the nova outburst is higher than is implied by the observed GeV light curves, due to downscattering of the gamma-rays by electrons within the ejecta. Novae provide ideal nearby laboratories to study magnetic field amplification and the onset of cosmic ray acceleration, because other time-dependent sources (e.g. radio supernovae) typically occur too distant to detect as gamma-ray sources. △ Less

Submitted 26 October, 2015; originally announced October 2015.

Comments: 9 pages, 4 figures, submitted to MNRAS

Radiative Transfer Models for Gamma-Ray Bursts

Authors: Indrek Vurm, Andrei M. Beloborodov

Abstract: We present global radiative transfer models for heated relativistic jets. The simulations include all relevant radiative processes, starting deep in the opaque zone and following the evolution of radiation to and beyond the photosphere of the jet. The transfer models are compared with three gamma-ray bursts GRB 990123, GRB 090902B, and GRB 130427A, which have well-measured and different spectra. T… ▽ More We present global radiative transfer models for heated relativistic jets. The simulations include all relevant radiative processes, starting deep in the opaque zone and following the evolution of radiation to and beyond the photosphere of the jet. The transfer models are compared with three gamma-ray bursts GRB 990123, GRB 090902B, and GRB 130427A, which have well-measured and different spectra. The models provide good fits to the observed spectra in all three cases, and we obtain estimates for the jet magnetization parameter $\varepsilon_{\rm B}$ and the Lorentz factor $Γ$. In the small sample of three bursts, $\varepsilon_{\rm B}$ varies between 0.01 and 0.05, and $Γ$ varies between 400 and 1200. △ Less

Submitted 26 July, 2016; v1 submitted 2 June, 2015; originally announced June 2015.

Comments: 18 pages, 8 figures, accepted to ApJ

Measuring Ambient Densities and Lorentz Factors of Gamma-Ray Bursts from GeV and Optical Observations

Authors: Romain Hascoët, Indrek Vurm, Andrei M. Beloborodov

Abstract: Fermi satellite discovered that cosmological gamma-ray bursts (GRBs) are accompanied by long GeV flashes. In two GRBs, an optical counterpart of the GeV flash has been detected. Recent work suggests that the GeV+optical flash is emitted by the external blast wave from the explosion in a medium loaded with copious $e^\pm$ pairs. The full light curve of the flash is predicted by a first-principle ra… ▽ More Fermi satellite discovered that cosmological gamma-ray bursts (GRBs) are accompanied by long GeV flashes. In two GRBs, an optical counterpart of the GeV flash has been detected. Recent work suggests that the GeV+optical flash is emitted by the external blast wave from the explosion in a medium loaded with copious $e^\pm$ pairs. The full light curve of the flash is predicted by a first-principle radiative transfer simulation and can be tested against observations. Here we examine a sample of 7 bursts with best GeV+optical data and test the model. We find that the observed light curves are in agreement with the theoretical predictions and allow us to measure three parameters for each burst: the Lorentz factor of the explosion, its isotropic kinetic energy, and the external density. With one possible exception of GRB 090510 (which is the only short burst in the sample) the ambient medium is consistent with a wind from a Wolf-Rayet progenitor. The wind density parameter $A=ρr^2$ varies in the sample around $10^{11}$g/cm. The initial Lorentz factor of the blast wave varies from 200 to 540 and correlates with the burst luminosity. Radiative efficiency of the prompt emission in the sample is between 0.1 and 0.8. For the two bursts with detected optical flash, GRB 120711A and GRB 130427A, we also estimate the magnetization of the external blast wave. Remarkably, the model reproduces the entire optical light curve of GRB 120711A (with its sharp peak, fast decay, plateau, and break) as well as the GeV data. The spectrum of GeV flashes is predicted to extend above 0.1 TeV, where they can be detected by ground-based Cherenkov telescopes. △ Less

Submitted 23 April, 2015; originally announced April 2015.

Comments: 21 pages, 9 figures, submitted to ApJ

Gamma-ray novae as probes of relativistic particle acceleration at non-relativistic shocks

Authors: Brian D. Metzger, Thomas Finzell, Indrek Vurm, Romain Hascoet, Andrei M. Beloborodov, Laura Chomiuk

Abstract: The Fermi LAT discovery that classical novae produce >100 MeV gamma-rays establishes that shocks and relativistic particle acceleration are key features of these events. These shocks are likely to be radiative due to the high densities of the nova ejecta at early times coincident with the gamma-ray emission. Thermal X-rays radiated behind the shock are absorbed by neutral gas and reprocessed into… ▽ More The Fermi LAT discovery that classical novae produce >100 MeV gamma-rays establishes that shocks and relativistic particle acceleration are key features of these events. These shocks are likely to be radiative due to the high densities of the nova ejecta at early times coincident with the gamma-ray emission. Thermal X-rays radiated behind the shock are absorbed by neutral gas and reprocessed into optical emission, similar to Type IIn (interacting) supernovae. Gamma-rays are produced by collisions between relativistic protons with the nova ejecta (hadronic scenario) or Inverse Compton/bremsstrahlung emission from relativistic electrons (leptonic scenario), where in both scenarios the efficiency for converting relativistic particle energy into LAT gamma-rays is at most a few tens of per cent. The ratio of gamma-ray and optical luminosities, L_gam/L_opt, thus sets a lower limit on the fraction of the shock power used to accelerate relativistic particles, e_nth. The measured values of L_gam/L_opt for two classical novae, V1324 Sco and V339 Del, constrains e_nth > 1e-2 and > 1e-3, respectively. Inverse Compton models for the gamma-ray emission are disfavored given the low electron acceleration efficiency, e_nth ~ 1e-4-1e-3, inferred from observations of Galactic cosmic rays and particle-in-cell (PIC) numerical simulations. A fraction > 100(0.01/e_nth) and > 10(0.01/e_nth) per cent of the optical luminosity is powered by shocks in V1324 Sco and V339 Del, respectively. Such high fractions challenge standard models that instead attribute all nova optical emission to the direct outwards transport of thermal energy released near the white dwarf surface. △ Less

Submitted 21 April, 2015; v1 submitted 21 January, 2015; originally announced January 2015.

Comments: 9 pages, 3 figures, 1 table, accepted to MNRAS after minor edits

Simulations of gamma-ray burst afterglows with a relativistic kinetic code

Authors: T. Pennanen, I. Vurm, J. Poutanen

Abstract: This paper introduces a kinetic code that simulates gamma-ray burst (GRB) afterglow emission from the external forward shock and presents examples of some of its applications. One interesting research topic discussed in the paper is the high-energy radiation produced by Compton scattering of the prompt GRB photons against the shock-accelerated electrons. The difference between the forward shock em… ▽ More This paper introduces a kinetic code that simulates gamma-ray burst (GRB) afterglow emission from the external forward shock and presents examples of some of its applications. One interesting research topic discussed in the paper is the high-energy radiation produced by Compton scattering of the prompt GRB photons against the shock-accelerated electrons. The difference between the forward shock emission in a wind-type and a constant-density medium is also studied, and the emission due to Maxwellian electron injection is compared to the case with pure power-law electrons. The code calculates the time-evolving photon and electron distributions in the emission region by solving the relativistic kinetic equations for each particle species. For the first time, the full relativistic equations for synchrotron emission/absorption, Compton scattering, and pair production/annihilation were applied to model the forward shock emission. The synchrotron self-absorption thermalization mechanism, which shapes the low-energy end of the electron distribution, was also included in the electron equation. The simulation results indicate that inverse Compton scattering of the prompt GRB photons can produce a luminous TeV emission component, even when pair production in the emission region is taken into account. This very high-energy radiation may be observable in low-redshift GRBs. The test simulations also show that the low-energy end of a pure power-law distribution of electrons can thermalize owing to synchrotron self-absorption in a wind-type environment, but without an observable impact on the radiation spectrum. Moreover, a flattening in the forward shock X-ray light curve may be expected when the electron injection function is assumed to be purely Maxwellian instead of a power law. △ Less

Submitted 21 March, 2014; originally announced March 2014.

Comments: 16 pages, 11 figures, accepted for publication in A&A

Journal ref: A&A, 564 (2014), A77

Shocks in nova outflows. I. Thermal emission

Authors: Brian D. Metzger, Romain Hascoet, Indrek Vurm, Andrei M. Beloborodov, Laura Chomiuk, J. L. Sokoloski, Thomas Nelson

Abstract: Evidence for shocks in nova outflows include (1) multiple velocity components in the optical spectra; (2) keV X-ray emission weeks to months after the outburst; (3) early radio flare on timescales of months, in excess of that predicted from the freely expanding photo-ionized gas; and (4) ~ GeV gamma-rays. We present a 1D model for the shock interaction between the fast nova outflow and a dense ext… ▽ More Evidence for shocks in nova outflows include (1) multiple velocity components in the optical spectra; (2) keV X-ray emission weeks to months after the outburst; (3) early radio flare on timescales of months, in excess of that predicted from the freely expanding photo-ionized gas; and (4) ~ GeV gamma-rays. We present a 1D model for the shock interaction between the fast nova outflow and a dense external shell (DES) and its associated thermal X-ray, optical, and radio emission. The forward shock is radiative initially when the density of shocked gas is highest, at which times radio emission originates from the dense cooling layer immediately downstream of the shock. The radio light curve is characterized by sharper rises to maximum and later peak times at progressively lower frequencies, with a peak brightness temperature that is approximately independent of frequency. We apply our model to the recent gamma-ray classical nova V1324 Sco, obtaining an adequate fit to the early radio maximum for reasonable assumptions about the fast nova outflow and assuming the DES possesses a velocity ~1e3 km/s and mass ~ 2e-4 M_sun; the former is consistent with the velocities of narrow line absorption systems observed previously in nova spectra, while the total ejecta mass of the DES and fast outflow is consistent with that inferred independently by modeling the late radio peak. Rapid evolution of the early radio light curves require the DES possess a steep outer density profile, which may indicate that the onset of mass loss from the white dwarf was rapid, providing indirect evidence that the DES was expelled by the thermonuclear runaway event. Reprocessed X-rays from the shock absorbed by the DES at early times may contribute significantly to the optical/UV emission, which we speculate is responsible for the previously unexplained `plateaus' and secondary maxima in nova optical light curves. △ Less

Submitted 6 March, 2014; originally announced March 2014.

Comments: 17 pages, 10 figures, submitted to MNRAS

Pair-dominated GeV-optical flash in GRB 130427A

Authors: Indrek Vurm, Romain Hascoet, Andrei M. Beloborodov

Abstract: We show that the light curve of the double GeV+optical flash in GRB 130427A is consistent with radiation from the blast wave in a wind-type medium with density parameter $A=ρr^2\sim 5\times 10^{10}$ g cm$^{-1}$. The peak of the flash is emitted by copious $e^\pm$ pairs created and heated in the blast wave; our first-principle calculation determines the pair-loading factor and temperature of the sh… ▽ More We show that the light curve of the double GeV+optical flash in GRB 130427A is consistent with radiation from the blast wave in a wind-type medium with density parameter $A=ρr^2\sim 5\times 10^{10}$ g cm$^{-1}$. The peak of the flash is emitted by copious $e^\pm$ pairs created and heated in the blast wave; our first-principle calculation determines the pair-loading factor and temperature of the shocked plasma. Using detailed radiative transfer simulations we reconstruct the observed double flash. The optical flash is dominated by synchrotron emission from the thermal plasma behind the forward shock, and the GeV flash is produced via inverse Compton (IC) scattering by the same plasma. The seed photons for IC scattering are dominated by the prompt MeV radiation during the first tens of seconds, and by the optical to X-ray afterglow thereafter. IC cooling of the thermal plasma behind the forward shock reproduces all GeV data from a few seconds to $\sim 1$ day. We find that the blast wave Lorentz factor at the peak of the flash is $Γ\approx 200$, and the forward shock magnetization is $ε_B\sim 2\times 10^{-4}$. An additional source is required by the data in the optical and X-ray bands at times $>10^2$ s; we speculate that this additional source may be a long-lived reverse shock in the explosion ejecta. △ Less

Submitted 9 May, 2014; v1 submitted 11 February, 2014; originally announced February 2014.

Comments: 5 pages, 3 figures, accepted to ApJL

Ionization Break-Out from Millisecond Pulsar Wind Nebulae: an X-ray Probe of the Origin of Superluminous Supernovae

Authors: Brian D. Metzger, Indrek Vurm, Romain Hascoet, Andrei M. Beloborodov

Abstract: Magnetic spin-down of a millisecond neutron star has been proposed as the power source of hydrogen-poor "superluminous" supernovae (SLSNe-I). However, producing an unambiguous test that can distinguish this model from alternatives, such as circumstellar interaction, has proven challenging. After the supernova explosion, the pulsar wind inflates a hot cavity behind the expanding stellar ejecta: the… ▽ More Magnetic spin-down of a millisecond neutron star has been proposed as the power source of hydrogen-poor "superluminous" supernovae (SLSNe-I). However, producing an unambiguous test that can distinguish this model from alternatives, such as circumstellar interaction, has proven challenging. After the supernova explosion, the pulsar wind inflates a hot cavity behind the expanding stellar ejecta: the nascent millisecond pulsar wind nebula. Electron/positron pairs injected by the wind cool through inverse Compton scattering and synchrotron emission, producing a pair cascade and hard X-ray spectrum inside the nebula. These X-rays ionize the inner exposed side of the ejecta, driving an ionization front that propagates outwards with time. Under some conditions this front can breach the ejecta surface within months after the optical supernova peak, allowing ~0.1-1 keV photons to escape the nebula unattenuated with a characteristic luminosity L_X ~ 1e43-1e45 erg/s. This "ionization break-out" may explain the luminous X-ray emission observed from the transient SCP 06F, providing direct evidence that this SLSN was indeed engine-powered. Luminous break-out requires a low ejecta mass and that the spin-down time of the pulsar be comparable to the photon diffusion timescale at optical maximum, the latter condition similar to that required for a supernova with a high optical fluence. These relatively special requirements may explain why most SLSNe-I are not accompanied by detectable X-ray emission. Global asymmetry of the supernova ejecta increases the likelihood of an early break-out along the direction of lowest density. Atomic states with lower threshold energies are more readily ionized at earlier times near optical maximum, allowing UV break-out across a wider range of pulsar and ejecta properties than X-ray break-out, possibly contributing to the blue/UV colors of SLSNe-I. △ Less

Submitted 30 July, 2013; originally announced July 2013.

Comments: 15 pages, 12 figures, 2 appendices, submitted to MNRAS

On the origin of GeV emission in gamma-ray bursts

Authors: Andrei M. Beloborodov, Romain Hascoet, Indrek Vurm

Abstract: The most common progenitors of gamma-ray bursts (GRBs) are massive stars with strong stellar winds. We show that the GRB blast wave in the wind should emit a bright GeV flash. It is produced by inverse Compton cooling of the thermal plasma behind the forward shock wave. The main part of the flash is shaped by scattering of the prompt MeV radiation (emitted at smaller radii) which streams through t… ▽ More The most common progenitors of gamma-ray bursts (GRBs) are massive stars with strong stellar winds. We show that the GRB blast wave in the wind should emit a bright GeV flash. It is produced by inverse Compton cooling of the thermal plasma behind the forward shock wave. The main part of the flash is shaped by scattering of the prompt MeV radiation (emitted at smaller radii) which streams through the external blast wave. The inverse-Compton flash is bright due to the huge e+- enrichment of the external medium. At late times, the blast wave switches to normal synchrotron-self-Compton cooling. The mechanism is demonstrated by a detailed transfer simulation. The observed prompt MeV radiation is taken as an input of the simulation; we use GRB 080916C as an example. The result reproduces the GeV flash observed by the Fermi telescope. It explains the delayed onset, the steep rise, the peak flux, the time of the peak, the long smooth decline, and the spectral slope of GeV emission. The wind density required to reproduce all these features is typical of Wolf-Rayet stars. Our simulation predicts strong TeV emission 1 min after the burst trigger; then a cutoff in the observed high-energy spectrum is expected from absorption by extragalactic background light. In addition, a bright optical counterpart of the GeV flash is predicted for plausible values of the magnetic field; such a double (optical+GeV) flash has been observed in GRB 130427A. △ Less

Submitted 28 January, 2014; v1 submitted 9 July, 2013; originally announced July 2013.

Comments: 20 pages, 9 figures, accepted to ApJ

Hot accretion flow in black hole binaries: a link connecting X-rays to the infrared

Authors: Alexandra Veledina, Juri Poutanen, Indrek Vurm

Abstract: Multiwavelength observations of Galactic black hole transients have opened a new path to understanding the physics of the innermost parts of the accretion flows. While the processes giving rise to their X-ray continuum have been studied extensively, the emission in the optical and infrared (OIR) energy bands was less investigated and remains poorly understood. The standard accretion disc, which ma… ▽ More Multiwavelength observations of Galactic black hole transients have opened a new path to understanding the physics of the innermost parts of the accretion flows. While the processes giving rise to their X-ray continuum have been studied extensively, the emission in the optical and infrared (OIR) energy bands was less investigated and remains poorly understood. The standard accretion disc, which may contribute to the flux at these wavelengths, is not capable of explaining a number of observables: the infrared excesses, fast OIR variability and a complicated correlation with the X-rays. It was suggested that these energy bands are dominated by the jet emission, however, this scenario does not work in a number of cases. We propose here an alternative, namely that most of the OIR emission is produced by the extended hot accretion flow. In this scenario, the OIR bands are dominated by the synchrotron radiation from the non-thermal electrons. An additional contribution is expected from the outer irradiated part of the accretion disc heated by the X-rays. We discuss properties of the model and compare them to the data. We show that the hot-flow scenario is consistent with many of the observed spectral data, at the same time naturally explaining X-ray timing properties, fast OIR variability and its correlation with the X-rays. △ Less

Submitted 5 February, 2013; v1 submitted 30 September, 2012; originally announced October 2012.

Comments: 17 pages, 12 figures; an expanded version, MNRAS, in press

Journal ref: MNRAS, 430 (2013), 3196-3212

On thermalization in gamma-ray burst jets and the peak energies of photospheric spectra

Authors: Indrek Vurm, Yuri Lyubarsky, Tsvi Piran

Abstract: The low energy spectral slopes of the prompt emission of most gamma-ray bursts (GRBs) are difficult to reconcile with radiatively efficient optically thin emission models irrespective of the radiation mechanism. An alternative is to ascribe the radiation around the spectral peak to a thermalization process occurring well inside the Thomson photosphere. This quasi-thermal spectrum can evolve into t… ▽ More The low energy spectral slopes of the prompt emission of most gamma-ray bursts (GRBs) are difficult to reconcile with radiatively efficient optically thin emission models irrespective of the radiation mechanism. An alternative is to ascribe the radiation around the spectral peak to a thermalization process occurring well inside the Thomson photosphere. This quasi-thermal spectrum can evolve into the observed non-thermal shape by additional energy release at moderate to small Thomson optical depths, which can readily give rise to the hard spectral tail. The position of the spectral peak is determined by the temperature and Lorentz factor of the flow in the termalization zone, where the total number of photons carried by the jet is established. To reach thermalization, dissipation alone is not sufficient and photon generation requires an efficient emission/absorption process in addition to scattering. We perform a systematic study of all relevant photon production mechanisms searching for possible conditions in which thermalization can take place. We find that a significant fraction of the available energy should be dissipated at intermediate radii, $\sim 10^{10}$-- a few$\times 10^{11}$ cm and the flow there should be relatively slow: the bulk Lorentz factor could not exceed a few tens for all but the most luminous bursts with the highest $\Epk$-s. The least restrictive constraint for successful thermalization, $Γ\lesssim 20$, is obtained if synchrotron emission acts as the photon source. This requires, however, a non-thermal acceleration deep below the Thomson photosphere transferring a significant fraction of the flow energy to relativistic electrons with Lorentz factors between 10 and 100. Other processes require bulk flow Lorentz factors of order of a few for typical bursts. We examine the implications of these results to different GRB photospheric emission models. △ Less

Submitted 21 December, 2012; v1 submitted 4 September, 2012; originally announced September 2012.

Comments: 18 pages, 10 figures, accepted to ApJ

A synchrotron self-Compton -- disk reprocessing model for optical/X-ray correlation in black hole X-ray binaries

Authors: Alexandra Veledina, Juri Poutanen, Indrek Vurm

Abstract: Physical picture of the emission mechanisms operating in the X-ray binaries was put under question by the simultaneous optical/X-ray observations with high time resolution. The light curves of the two energy bands appeared to be connected and the cross-correlation functions observed in three black hole binaries exhibited a complicated shape. They show a dip of the optical emission a few seconds be… ▽ More Physical picture of the emission mechanisms operating in the X-ray binaries was put under question by the simultaneous optical/X-ray observations with high time resolution. The light curves of the two energy bands appeared to be connected and the cross-correlation functions observed in three black hole binaries exhibited a complicated shape. They show a dip of the optical emission a few seconds before the X-ray peak and the optical flare just after the X-ray peak. This behavior could not be explained in terms of standard optical emission candidates (e.g., emission from the cold accretion disk or a jet). We propose a novel model, which explains the broadband optical to the X-ray spectra and the variability properties. We suggest that the optical emission consists of two components: synchrotron radiation from the non-thermal electrons in the hot accretion flow and the emission produced by reprocessing of the X-rays in the outer part of the accretion disk. The first component is anti-correlated with the X-rays, while the second one is correlated, but delayed and smeared relative to the X-rays. The interplay of the components explains the complex shape of the cross-correlation function, the features in the optical power spectral density as well as the time lags. △ Less

Submitted 27 June, 2011; v1 submitted 13 May, 2011; originally announced May 2011.

Comments: 5 pages, 3 figures, ApJ Letters in press

Journal ref: Astroph, J. Lett,, 737:L17, 2011

Gamma-ray bursts from magnetized collisionally-heated jets

Authors: Indrek Vurm, Andrei M. Beloborodov, Juri Poutanen

Abstract: Jets producing gamma-ray bursts (GRBs) are likely to carry a neutron component that drifts with respect to the proton component. The neutron-proton collisions strongly heat the jet and generate electron-positron pairs. We investigate radiation produced by this heating using a new numerical code. Our results confirm the recent claim that collisional heating generates the observed Band-type spectrum… ▽ More Jets producing gamma-ray bursts (GRBs) are likely to carry a neutron component that drifts with respect to the proton component. The neutron-proton collisions strongly heat the jet and generate electron-positron pairs. We investigate radiation produced by this heating using a new numerical code. Our results confirm the recent claim that collisional heating generates the observed Band-type spectrum of GRBs. We extend the model to study the effects of magnetic fields on the emitted spectrum. We find that the spectrum peak remains near 1 MeV for the entire range of the magnetization parameter $0<ε_B<2$ that is explored in our simulations. The low-energy part of the spectrum softens with increasing $ε_B$, and a visible soft excess appears in the keV band. The high-energy part of the spectrum extends well above the GeV range and can contribute to the prompt emission observed by Fermi/LAT. Overall, the radiation spectrum created by the collisional mechanism appears to agree with observations, with no fine-tuning of parameters. △ Less

Submitted 26 June, 2011; v1 submitted 3 April, 2011; originally announced April 2011.

Comments: 13 pages, 6 figures, accepted to ApJ

Journal ref: Astrophys. J., 738:77, 2011

A self-consistent hybrid Comptonization model for broad-band spectra of accreting supermassive black holes

Authors: Alexandra Veledina, Indrek Vurm, Juri Poutanen

Abstract: The nature of the broad-band spectra of supermassive accreting black holes in active galactic nuclei (AGNs) is still unknown. The hard X-ray spectra of Seyferts as well as of Galactic stellar-mass black holes (GBHs) are well represented by thermal Comptonization, but the origin of the seed photons is less certain. The MeV tails observed in GBHs provide evidence in favour of non-thermal electron ta… ▽ More The nature of the broad-band spectra of supermassive accreting black holes in active galactic nuclei (AGNs) is still unknown. The hard X-ray spectra of Seyferts as well as of Galactic stellar-mass black holes (GBHs) are well represented by thermal Comptonization, but the origin of the seed photons is less certain. The MeV tails observed in GBHs provide evidence in favour of non-thermal electron tails and it is possible that such electrons are also present in the X-ray emitting regions of AGNs. Using simulations with the kinetic code that self-consistently models electron and photon distributions, we find that the power-law-like X-ray spectra in AGNs can be explained in terms of the synchrotron self-Compton radiation of hybrid thermal/non-thermal electrons, similarly to the hard/low state of GBHs. Under a very broad range of parameters the model predicts a rather narrow distribution of photon spectral slopes consistent with that observed from LINERs and Seyferts at luminosities less than 3 per cent of the Eddington luminosity. The entire infrared to X-ray spectrum of these objects can be described in terms of our model, suggesting a tight correlation between the two energy bands. We show that the recently found correlation between slope and the Eddington ratio at higher luminosities can be described by the increasing fraction of disc photons in the emitting region, which may be associated with the decreasing inner radius of the optically thick accretion disc. The increasing flux of soft photons is also responsible for the transformation of the electron distribution from nearly thermal to almost completely non-thermal. The softer X-ray spectra observed in narrow-line Seyfert galaxies may correspond to non-thermal Comptonization of the disc photons, predicting that no cutoff should be observed up to MeV energies in these sources, similarly to the soft-state GBHs. △ Less

Submitted 9 September, 2011; v1 submitted 2 December, 2010; originally announced December 2010.

Comments: MNRAS published

Journal ref: MNRAS, 414: 3330-3343, 2011

Theory of Compton scattering by anisotropic electrons

Authors: Juri Poutanen, Indrek Vurm

Abstract: Compton scattering plays an important role in various astrophysical objects such as accreting black holes and neutron stars, pulsars, and relativistic jets, clusters of galaxies as well as the early Universe. In most of the calculations it is assumed that the electrons have isotropic angular distribution in some frame. However, there are situations where the anisotropy may be significant due to th… ▽ More Compton scattering plays an important role in various astrophysical objects such as accreting black holes and neutron stars, pulsars, and relativistic jets, clusters of galaxies as well as the early Universe. In most of the calculations it is assumed that the electrons have isotropic angular distribution in some frame. However, there are situations where the anisotropy may be significant due to the bulk motions, or anisotropic cooling by synchrotron radiation, or anisotropic source of seed soft photons. We develop here an analytical theory of Compton scattering by anisotropic distribution of electrons that can simplify significantly the calculations. Assuming that the electron angular distribution can be represented by a second order polynomial over cosine of some angle (dipole and quadrupole anisotropy), we integrate the exact Klein-Nishina cross-section over the angles. Exact analytical and approximate formulae valid for any photon and electron energies are derived for the redistribution functions describing Compton scattering of photons with arbitrary angular distribution by anisotropic electrons. The analytical expressions for the corresponding photon scattering cross-section on such electrons as well as the mean energy of scattered photons, its dispersion and radiation pressure force are also derived. We applied the developed formalism to the accurate calculations of the thermal and kinematic Sunyaev-Zeldovich effects for arbitrary electron distributions. △ Less

Submitted 11 June, 2010; originally announced June 2010.

Comments: 23 pages, 12 figures, ApJ Supplement Series, in press

Journal ref: Astrophys. J. Suppl. Ser.,189:286-308, 2010

On the origin of spectral states in accreting black holes

Authors: Juri Poutanen, Indrek Vurm

Abstract: The origin of dramatically different electron distributions responsible for Comptonization in black hole X-ray binaries (BHBs) in their various states is discussed. We solve the coupled kinetic equations for photons and electrons without approximations on the relevant cross-sections accounting for Compton scattering, synchrotron radiation, and Coulomb collisions. In the absence of external soft… ▽ More The origin of dramatically different electron distributions responsible for Comptonization in black hole X-ray binaries (BHBs) in their various states is discussed. We solve the coupled kinetic equations for photons and electrons without approximations on the relevant cross-sections accounting for Compton scattering, synchrotron radiation, and Coulomb collisions. In the absence of external soft photons, the electrons are efficiently thermalized by synchrotron self-absorption and Coulomb scattering even for pure nonthermal electron injection. The resulting quasi-thermal synchrotron self-Compton spectra have very stable slopes and electron temperatures similar to the hard states of BHBs. The hard spectral slopes observed in the X-rays, the cutoff at 100 keV, and the MeV tail together require low magnetic fields ruling out the magnetic dissipation mechanism. The motion of the accretion disk toward the black hole results in larger Compton cooling and lower equilibrium electron temperature. Our self-consistent simulations show that in this case both electron and photon distributions attain a power-law-dominated shape similar to what is observed in the soft state. The electron distribution in the Cyg X-1 soft state might require a strong magnetic field, being consistent with the magnetically dominated corona. △ Less

Submitted 15 December, 2008; v1 submitted 19 July, 2008; originally announced July 2008.

Comments: 4 pages, 2 figures; version published in ApJ Letters

Journal ref: Astrophys.J.Lett.690:L97-L100,2009

Time-dependent modeling of radiative processes in hot magnetized plasmas

Authors: Indrek Vurm, Juri Poutanen

Abstract: Numerical simulations of radiative processes in magnetized compact sources such as hot accretion disks around black holes, relativistic jets in active galaxies and gamma-ray bursts are complicated because the particle and photon distributions span many orders of magnitude in energy, they also strongly depend on each other, the radiative processes behave significantly differently depending on the… ▽ More Numerical simulations of radiative processes in magnetized compact sources such as hot accretion disks around black holes, relativistic jets in active galaxies and gamma-ray bursts are complicated because the particle and photon distributions span many orders of magnitude in energy, they also strongly depend on each other, the radiative processes behave significantly differently depending on the energy regime, and finally due to the enormous difference in the time-scales of the processes. We have developed a novel computer code for the time-dependent simulations that overcomes these problems. The processes taken into account are Compton scattering, electron-positron pair production and annihilation, Coulomb scattering as well as synchrotron emission and absorption. No approximation has been made on the corresponding rates. For the first time, we solve coupled integro-differential kinetic equations for photons and electrons/positrons without any limitations on the photon and lepton energies. A numerical scheme is proposed to guarantee energy conservation when dealing with synchrotron processes in electron and photon equations. We apply the code to model non-thermal pair cascades in the blackbody radiation field, to study the synchrotron self-absorption as particle thermalization mechanism, and to simulate time evolution of stochastically heated pairs and corresponding synchrotron self-Compton photon spectra which might be responsible for the prompt emission of gamma-ray bursts. Good agreement with previous works is found in the parameter regimes where comparison is feasible, with the differences attributable to our improved treatment of the microphysics. △ Less

Submitted 3 March, 2009; v1 submitted 16 July, 2008; originally announced July 2008.

Comments: 23 pages, 4 figures, ApJ, in press

Journal ref: Astrophys.J.698:293-316,2009

Electron thermalization and photon emission from magnetized compact sources

Authors: Indrek Vurm, Juri Poutanen

Abstract: We present detailed calculations of the electron thermalization by synchrotron self-absorption accounting for cooling by Compton scattering. For the first time, we solve coupled kinetic equations for electrons and photons without any approximations on the relevant cross-sections and compute self-consistently the resulting electron and photon distributions. The results presented in this contribut… ▽ More We present detailed calculations of the electron thermalization by synchrotron self-absorption accounting for cooling by Compton scattering. For the first time, we solve coupled kinetic equations for electrons and photons without any approximations on the relevant cross-sections and compute self-consistently the resulting electron and photon distributions. The results presented in this contribution may be applied to the magnetized coronae around Galactic black hole accretion disks. △ Less

Submitted 25 February, 2008; originally announced February 2008.

Comments: 4 pages, 8 figures; Proceedings of the workshop 'High Energy Phenomena in Relativistic Outflows' (HEPRO), Dublin, Ireland

Journal ref: Int.J.Mod.Phys.D17:1629-1634,2008