Showing 1–50 of 54 results for author: Barkov, M V

Numerical simulation of photospheric emission in long gamma-ray bursts: prompt correlations, spectral shapes, and polarizations

Authors: Hirotaka Ito, Jin Matsumoto, Shigehiro Nagataki, Donald C. Warren, Maxim V. Barkov, Daisuke Yonetoku

Abstract: We explore the properties of photospheric emission in the context of long gamma-ray bursts (LGRBs) using three numerical models that combine relativistic hydrodynamical simulations and Monte Carlo radiation transfer calculations in three dimensions. Our simulations confirm that the photospheric emission gives rise to correlations between the spectral peak energy and luminosity that agree with the… ▽ More We explore the properties of photospheric emission in the context of long gamma-ray bursts (LGRBs) using three numerical models that combine relativistic hydrodynamical simulations and Monte Carlo radiation transfer calculations in three dimensions. Our simulations confirm that the photospheric emission gives rise to correlations between the spectral peak energy and luminosity that agree with the observed Yonetoku, Amati, and Golenetskii correlations. It is also shown that the spectral peak energy and luminosity correlate with the bulk Lorentz factor, as indicated in the literature. On the other hand, synthetic spectral shapes tend to be narrower than those of the observations. The result indicates that an additional physical process that can provide non-thermal broadening is needed to reproduce the spectral features. Furthermore, the polarization analysis finds that, while the degree of polarization is low for the emission from the jet core ($Π< 4~\%$), it tends to increase with the viewing angle outside the core and can be as high as $Π\sim 20-40~\%$ in an extreme case. This suggests that the typical GRBs show systematically low polarization compared to softer, dimmer counterparts (X-ray-rich GRBs and X-ray flashes). Interestingly, our simulations indicate that photospheric emission exhibits large temporal variation in the polarization position angle ($Δψ\sim 90^{\circ}$), which may be compatible with those inferred in observations. A notable energy dependence of the polarization property is another characteristic feature found in the current study. Particularly, the difference in the position angle among different energy bands can be as large as $\sim 90^{\circ}$. △ Less

Submitted 19 July, 2023; originally announced July 2023.

Comments: 20 pages, 12 figures, accepted for publication in ApJ

Report number: RIKEN-iTHEMS-Report-23

3D Relativistic MHD simulations of the gamma-ray binaries

Authors: Maxim V. Barkov, Evgeniy Kalinin, Maxim Lyutikov

Abstract: In gamma-ray binaries neutron star is orbiting a companion that produces a strong stellar wind. We demonstrate that observed properties of "stellar wind"-"pulsar wind" interaction depend both on the overall wind thrust ratio, as well as more subtle geometrical factors: the relative direction of the pulsar's spin, the plane of the orbit, the direction of motion, and the instantaneous line of sight.… ▽ More In gamma-ray binaries neutron star is orbiting a companion that produces a strong stellar wind. We demonstrate that observed properties of "stellar wind"-"pulsar wind" interaction depend both on the overall wind thrust ratio, as well as more subtle geometrical factors: the relative direction of the pulsar's spin, the plane of the orbit, the direction of motion, and the instantaneous line of sight. Using fully 3D relativistic magnetohydrodynamical simulations we find that the resulting intrinsic morphologies can be significantly orbital phase-dependent: a given system may change from tailward-open to tailward-closed shapes. As a result, the region of unshocked pulsar wind can change by an order of magnitude over a quarter of the orbit. We calculate radiation maps and synthetic light curves for synchrotron (X-ray) and Inverse-Compton emission (GeV-TeV), taking into account $γ-γ$ absorption. Our modeled light curves are in agreement with the phase-dependent observed light curves of LS5039. △ Less

Submitted 25 June, 2024; v1 submitted 22 November, 2022; originally announced November 2022.

Comments: 22 pages, 20 figures

Relativistic magnetic explosions

Authors: Maxim V. Barkov, Praveen Sharma, Konstantinos N. Gourgouliatos, Maxim Lyutikov

Abstract: Many explosive astrophysical events, like magnetars' bursts and flares, are magnetically driven. We consider dynamics of such magnetic explosions - relativistic expansion of highly magnetized and highly magnetically over-pressurized clouds. The corresponding dynamics is qualitatively different from fluid explosions due to the topological constraint of the conservation of the magnetic flux. Using a… ▽ More Many explosive astrophysical events, like magnetars' bursts and flares, are magnetically driven. We consider dynamics of such magnetic explosions - relativistic expansion of highly magnetized and highly magnetically over-pressurized clouds. The corresponding dynamics is qualitatively different from fluid explosions due to the topological constraint of the conservation of the magnetic flux. Using analytical, relativistic MHD as well as force-free calculations, we find that the creation of a relativistically expanding, causally disconnected flow obeys a threshold condition: it requires sufficiently high initial over-pressure and sufficiently quick decrease of the pressure in the external medium (the pre-explosion wind). In the subcritical case the magnetic cloud just "puffs-up" and quietly expands with the pre-flare wind. We also find a compact analytical solution to the Prendergast's problem - expansion of force-free plasma into vacuum. △ Less

Submitted 2 May, 2022; originally announced May 2022.

Comments: 31 pages, 23 figures, submitted to ApJ

Formation of periodic FRB in binary systems with eccentricity

Authors: Maxim V. Barkov, Sergei B. Popov

Abstract: Long-term periodicity in the rate of flares is observed for two repeating sources of fast radio bursts (FRBs). In this paper We present a hydrodynamical modeling of a massive binary consisting of a magnetar and an early-type star. We model the interaction of the pulsar wind from the magnetar with an intense stellar wind. It is shown that only during a fraction of the orbital period radio emission… ▽ More Long-term periodicity in the rate of flares is observed for two repeating sources of fast radio bursts (FRBs). In this paper We present a hydrodynamical modeling of a massive binary consisting of a magnetar and an early-type star. We model the interaction of the pulsar wind from the magnetar with an intense stellar wind. It is shown that only during a fraction of the orbital period radio emission can escape the system. This explains the duty cycle of the two repeating FRB sources with periodic activity. The width of the transparency window depends on the eccentricity, stellar wind properties, and the viewing angle. To describe properties of the known sources it is necessary to assume large eccentricities $\gtrsim 0.5$. We apply the maser cyclotron mechanism of the radio emission generation to model spectral properties of the sources. The produced spectrum is not wide: $Δν/ν\sim 0.2$ and the typical frequency depends on the radius of the shock where the emission is generated. The shock radius changes along the orbit. This, together with changing parameters of the medium, allows us to explain the frequency drift during the phase of visibility. Frequency dependence of the degree of polarization at few GHz can be a consequence of a small scale turbulence in the shocked stellar wind. It is much more difficult to explain huge ($\sim 10^5$ [rad/m$^2$]) and variable value of the rotation measure observed for FRB 121102. We suggest that this can be explained if the supernova explosion which produced the magnetar happened near a dense interstellar cloud with $n \sim100$ cm$^{-3}$. △ Less

Submitted 6 June, 2022; v1 submitted 28 April, 2022; originally announced April 2022.

Comments: 12 pages, 9 figures, accepted by MNRAS

The major role of eccentricity in the evolution of colliding pulsar-stellar winds

Authors: Maxim V. Barkov, Valenti Bosch-Ramon

Abstract: Binary systems that host a massive star and a non-accreting pulsar can be powerful non-thermal emitters. The relativistic pulsar wind and the non-relativistic stellar outflows interact along the orbit, producing ultrarelativistic particles that radiate from radio to gamma rays. To properly characterize the physics of these sources, and better understand their emission and impact on the environment… ▽ More Binary systems that host a massive star and a non-accreting pulsar can be powerful non-thermal emitters. The relativistic pulsar wind and the non-relativistic stellar outflows interact along the orbit, producing ultrarelativistic particles that radiate from radio to gamma rays. To properly characterize the physics of these sources, and better understand their emission and impact on the environment, careful modelling of the outflow interactions, spanning a broad range of spatial and temporal scales, is needed. Full 3-dimensional approaches are very computationally expensive, but simpler approximate approaches, while still realistic at the semi-quantitative level, are available. We present here the results of calculations done with a quasi 3-dimensional scheme to compute the evolution of the interacting flows in a region spanning in size up to a thousand times the size of the binary. In particular, we analyze for the first time the role of different eccentricities in the large scale evolution of the shocked flows. We find that the higher the eccentricity, the closer the flows behave like a one-side outflow, which becomes rather collimated for eccentricity values $\gtrsim 0.75$. The simulations also unveil that the pulsar and the stellar winds become fully mixed within the grid for low eccentricity systems, presenting a more stochastic behavior at large scales than in the highly eccentric systems. △ Less

Submitted 15 April, 2022; originally announced April 2022.

Comments: 18 pages, 11 figures

Journal ref: July 2021, Universe, vol. 7, issue 8, p. 277

Relativistic hydrodynamical simulations of the effects of the stellar wind and the orbit on high-mass microquasar jets

Authors: Maxim V. Barkov, Valenti Bosch-Ramon

Abstract: High-mass microquasar jets, produced in an accreting compact object in orbit around a massive star, must cross a region filled with stellar wind. The combined effects of the wind and orbital motion can strongly affect the jet properties on binary scales and beyond. The study of such effects can shed light on how high-mass microquasar jets propagate and terminate in the interstellar medium. We stud… ▽ More High-mass microquasar jets, produced in an accreting compact object in orbit around a massive star, must cross a region filled with stellar wind. The combined effects of the wind and orbital motion can strongly affect the jet properties on binary scales and beyond. The study of such effects can shed light on how high-mass microquasar jets propagate and terminate in the interstellar medium. We study for the first time, using relativistic hydrodynamical simulations, the combined impact of the stellar wind and orbital motion on the properties of high-mass microquasar jets on binary scales and beyond. We have performed 3-dimensional relativistic hydrodynamic simulations, using the PLUTO code, of a microquasar scenario in which a strong weakly relativistic wind from a star interacts with a relativistic jet under the effect of the binary orbital motion. The parameters of the orbit are chosen such that the results can provide insight on the jet-wind interaction in compact systems like for instance Cyg~X-1 or Cyg~X-3. The wind and jet momentum rates are set to values that may be realistic for these sources and lead to moderate jet bending, which together with the close orbit and jet instabilities could trigger significant jet precession and disruption. For high-mass microquasars with orbit size $a\sim 0.1$~AU, and (relativistic) jet power $L_j\sim 10^{37}(\dot M_w/10^{-6}\,{\rm M}_\odot\,{\rm yr}^{-1})$~erg~s$^{-1}$, where $\dot M_w$ is the stellar wind mass rate, the combined effects of the stellar wind and orbital motion can induce relativistic jet disruption on scales $\sim 1$~AU. △ Less

Submitted 8 December, 2021; originally announced December 2021.

Comments: 16 pages, 13 figures, MNRAS accepted

A semi-analytic afterglow with thermal electrons and synchrotron self-Compton emission

Authors: Donald C. Warren, Maria Dainotti, Maxim V. Barkov, Bjorn Ahlgren, Hirotaka Ito, Shigehiro Nagataki

Abstract: We extend previous work on gamma-ray burst (GRB) afterglows involving hot thermal electrons at the base of a shock-accelerated tail. Using a physically-motivated electron distribution based on first-principles simulations, we compute broadband emission from radio to TeV gamma-rays. For the first time, we present the effects of a thermal distribution of electrons on synchrotron self-Compton (SSC) e… ▽ More We extend previous work on gamma-ray burst (GRB) afterglows involving hot thermal electrons at the base of a shock-accelerated tail. Using a physically-motivated electron distribution based on first-principles simulations, we compute broadband emission from radio to TeV gamma-rays. For the first time, we present the effects of a thermal distribution of electrons on synchrotron self-Compton (SSC) emission. The presence of thermal electrons causes temporal and spectral structure across the entire observable afterglow, which is substantively different from models that assume a pure power-law distribution for the electrons. We show that early-time TeV emission is enhanced by more than an order of magnitude for our fiducial parameters, with a time-varying spectral index that does not occur for a pure power law of electrons. We further show that the X-ray "closure relations" take a very different, also time-dependent, form when thermal electrons are present; the shape traced out by the X-ray afterglows is a qualitative match to observations of the traditional decay phase. △ Less

Submitted 15 September, 2021; originally announced September 2021.

Comments: 32 pages, 15 figures. Includes changes requested by referee

Report number: RIKEN-iTHEMS-Report-21

Fast radio bursts by high-frequency synchrotron maser emission generated at the reverse shock of a powerful magnetar flare

Authors: Dmitry Khangulyan, Maxim V. Barkov, Sergey B. Popov

Abstract: We consider a magnetar flare model for fast radio bursts (FRBs). We show that millisecond burst of sufficient power can be generated by synchrotron maser emission ignited at the reverse shock propagating through the weakly magnetized material that forms the magnetar flare. If the maser emission is generated in an anisotropic regime (due to the geometry of the production region or presence of an in… ▽ More We consider a magnetar flare model for fast radio bursts (FRBs). We show that millisecond burst of sufficient power can be generated by synchrotron maser emission ignited at the reverse shock propagating through the weakly magnetized material that forms the magnetar flare. If the maser emission is generated in an anisotropic regime (due to the geometry of the production region or presence of an intense external source of stimulating photons) the duration of the maser flashes is similar to the magnetar flare duration even if the shock front radius is large. Our scenario allows relaxing the requirements for several key parameters: the magnetic field strength at the production site, luminosity of the flare, and the production site bulk Lorentz factor. To check the feasibility of this model, we study the statistical relation between powerful magnetar flares and the rate of FRBs. The expected ratio is derived by convoluting the redshift-dependent magnetar density with their flare luminosity function above the energy limit determined by the FRB detection threshold. We obtain that only a small fraction, $\sim10^{-5}$, of powerful magnetar flares trigger FRBs. This ratio agrees surprisingly well with our estimates: we obtained that $10\%$ of magnetars should be in the evolutionary phase suitable for the production of FRBs, and only $10^{-4}$ of all flares are expected to be weakly magnetized, which is a necessary condition for the high-frequency maser emission. △ Less

Submitted 16 February, 2022; v1 submitted 17 June, 2021; originally announced June 2021.

Comments: 16 pages, 2 figure, ApJ in press

The maximum energy of shock-accelerated electrons in a microturbulent magnetic field

Authors: Donald C. Warren, Catherine A. A. Beauchemin, Maxim V. Barkov, Shigehiro Nagataki

Abstract: Relativistic shocks propagating into a medium with low magnetization are generated and sustained by small-scale but very strong magnetic field turbulence. This so-called "microturbulence" modifies the typical shock acceleration process, and in particular that of electrons. In this work we perform Monte Carlo (MC) simulations of electrons encountering shocks with microturbulent fields. The simulati… ▽ More Relativistic shocks propagating into a medium with low magnetization are generated and sustained by small-scale but very strong magnetic field turbulence. This so-called "microturbulence" modifies the typical shock acceleration process, and in particular that of electrons. In this work we perform Monte Carlo (MC) simulations of electrons encountering shocks with microturbulent fields. The simulations cover a three-dimensional parameter space in shock speed, acceleration efficiency, and peak magnetic field strength. From these, a Markov Chain Monte Carlo (MCMC) method was employed to estimate the maximum electron momentum from the MC-simulated electron spectra. Having estimated this quantity at many points well-distributed over an astrophysically relevant parameter space, an MCMC method was again used to estimate the parameters of an empirical formula that computes the maximum momentum of a Fermi-accelerated electron population anywhere in this parameter space. The maximum energy is well-approximated as a broken power-law in shock speed, with the break occurring when the shock decelerates to the point where electrons can begin to escape upstream from the shock. △ Less

Submitted 30 October, 2020; v1 submitted 13 October, 2020; originally announced October 2020.

Comments: 18 pages, 2 figures, accepted in ApJ; arXiv version includes changes requested by referee

Report number: RIKEN-iTHEMS-Report-20

On the nature of radio filaments near the Galactic Center

Authors: Maxim V. Barkov, Maxim Lyutikov

Abstract: We suggest that narrow, long radio filaments near the Galactic Center arise as kinetic jets - streams of high energy particles escaping from ram-pressure confined pulsar wind nebulae (PWNe). The reconnection between the PWN and interstellar magnetic field allows pulsar wind particles to escape, creating long narrow features. They are the low frequency analogs of kinetic jets seen around some fast-… ▽ More We suggest that narrow, long radio filaments near the Galactic Center arise as kinetic jets - streams of high energy particles escaping from ram-pressure confined pulsar wind nebulae (PWNe). The reconnection between the PWN and interstellar magnetic field allows pulsar wind particles to escape, creating long narrow features. They are the low frequency analogs of kinetic jets seen around some fast-moving pulsars, such as The Guitar and The Lighthouse PWNe. The radio filaments trace a population of pulsars also responsible for the Fermi GeV excess produced by the Inverse Compton scattering by the pulsar wind particles. The magnetic flux tubes are stretched radially by the large scale Galactic winds. In addition to PWNe accelerated particles can be injected at supernovae remnants. The model predicts variations of the structure of the largest filaments on scales of $\sim$ dozens of years - smaller variations can occur on shorter time scales. We also encourage targeted observations of the brightest sections of the filaments and of the related unresolved point sources in search of the powering PWNe and pulsars. △ Less

Submitted 30 May, 2019; originally announced May 2019.

Comments: 5 pages, 2 figures

The photospheric origin of the Yonetoku relation in gamma-ray bursts

Authors: Hirotaka Ito, Jin Matsumoto, Shigehiro Nagataki, Donald C. Warren, Maxim V. Barkov, Daisuke Yonetoku

Abstract: Long duration gamma-ray bursts (GRBs), the brightest events since the Big Bang itself, are believed to originate in an ultra-relativistic jet breaking out from a massive stellar envelope. Despite decades of study, there is still no consensus on their emission mechanism. One unresolved question is the origin of the tight correlation between the spectral peak energy Ep and peak luminosity Lp discove… ▽ More Long duration gamma-ray bursts (GRBs), the brightest events since the Big Bang itself, are believed to originate in an ultra-relativistic jet breaking out from a massive stellar envelope. Despite decades of study, there is still no consensus on their emission mechanism. One unresolved question is the origin of the tight correlation between the spectral peak energy Ep and peak luminosity Lp discovered in observations. This Yonetoku relation is the tightest correlation found in the properties of the prompt phase of GRB emission, providing the best diagnostic for the radiation mechanism. Here we present 3D hydrodynamical simulations, and post-process radiation transfer calculations, of photospheric emission from a relativistic jet. Our simulations reproduce the Yonetoku relation as a natural consequence of viewing angle. Although jet dynamics depend sensitively on luminosity, the Ep-Lp correlation holds regardless. This result strongly suggests that photospheric emission is the dominant component in the prompt phase of GRBs. △ Less

Submitted 3 April, 2019; v1 submitted 2 June, 2018; originally announced June 2018.

Comments: 18 pages, 5 figures, The original submitted version of an article published in Nature Communications. The final authenticated version is available online at: https://www.nature.com/articles/s41467-019-09281-z

Report number: RIKEN-iTHEMS-Report-18

Journal ref: Nature Communications volume 10, Article number: 1504 (2019)

Prediction of the second peak in the afterglow of GW170817

Authors: Maxim V. Barkov, Adithan Kathirgamaraju, Yonggang Luo, Maxim Lyutikov, Dimitrios Giannios

Abstract: We performed calculations of the late radio and X-ray afterglow of GRB/GW170817 in the cocoon-jet paradigm, predicting the appearance of a second peak in the afterglow light curve ~ one-three years after the explosion. The model assumes that the prompt emission and early afterglows originate from a cocoon generated during break-out of the delayed magnetically powered jet. As the jet breaks out fro… ▽ More We performed calculations of the late radio and X-ray afterglow of GRB/GW170817 in the cocoon-jet paradigm, predicting the appearance of a second peak in the afterglow light curve ~ one-three years after the explosion. The model assumes that the prompt emission and early afterglows originate from a cocoon generated during break-out of the delayed magnetically powered jet. As the jet breaks out from the torus-generated wind, a nearly isotropic mildly relativistic outflow is generated; at the same time the primary jet accelerates to high Lorentz factors and avoids detection. As the fast jet slows down, it should become visible to the off-axis observer. Thus, the model has a clear prediction: the X-ray and radio afterglows should first experience a decay, as the cocoon slows down, followed by a rebrightening when the primary jet starts emitting toward an observer. △ Less

Submitted 21 May, 2018; originally announced May 2018.

Comments: 6 pages, 4 figures, submitted to MNRAS

Kinetic "jets" from fast moving pulsars

Authors: Maxim V. Barkov, Maxim Lyutikov, Noel Klingler, Pol Bordas

Abstract: Some fast-moving pulsars, like the Guitar and the Lighthouse, exhibit asymmetric non-thermal emission features that extend well beyond their ram pressure confined pulsar wind nebulae (PWNe). Using 3D relativistic simulations we explain these features as kinetically streaming pulsar wind particles that escaped into the interstellar medium (ISM) due to reconnection between the PWN and ISM magnetic f… ▽ More Some fast-moving pulsars, like the Guitar and the Lighthouse, exhibit asymmetric non-thermal emission features that extend well beyond their ram pressure confined pulsar wind nebulae (PWNe). Using 3D relativistic simulations we explain these features as kinetically streaming pulsar wind particles that escaped into the interstellar medium (ISM) due to reconnection between the PWN and ISM magnetic fields. The structure of the reconnecting magnetic fields at the incoming and outgoing regions produce highly asymmetric magnetic bottles, and result in asymmetric extended features. For the features to become visible, the ISM magnetic field should be sufficiently high, $B_{\rm ISM}>10$~$μ$G. We also discuss archival observations of PWNe displaying evidence of kinetic jets: the Dragonfly PWN (PSR J2021+3651), G327.1--1.1, and MSH 11--62, the latter two of which exhibit "snail eyes" morphologies. We suggest that in those cases the pulsar is moving along the ambient magnetic field in a frisbee-type configuration. △ Less

Submitted 23 April, 2018; v1 submitted 19 April, 2018; originally announced April 2018.

Comments: 13 pages, 12 figures, submitted to MNRAS

3D dynamics and morphology of bow-shock Pulsar Wind Nebulae

Authors: Maxim V. Barkov, Maxim Lyutikov, Dmitry Khangulyan

Abstract: Bow-shock pulsar wind nebulae (PWNe) show a variety of morphological shapes. We attribute this diversity to the geometrical factors: relative orientations of the pulsar rotation axis, proper velocity, and the line of sight (magnetic inclination angle may also have a certain influence on the morphology). We identify three basic types of bow-shock nebulae: (i) a "Rifle Bullet" (pulsar spin axis and… ▽ More Bow-shock pulsar wind nebulae (PWNe) show a variety of morphological shapes. We attribute this diversity to the geometrical factors: relative orientations of the pulsar rotation axis, proper velocity, and the line of sight (magnetic inclination angle may also have a certain influence on the morphology). We identify three basic types of bow-shock nebulae: (i) a "Rifle Bullet" (pulsar spin axis and proper velocity are aligned); (ii) a "Frisbee" (pulsar spin axis and proper velocity are orthogonal with the spin axis lying in the plane of the sky), and (iii) a Cart Wheel" (like frisbee but the spin axis is perpendicular to the plane of the sky). Using 3D RMHD simulations, as well as analytical calculations, we reproduce the key morphological features of the bow-shock PWNe, as well as variations, are seen across different systems. magnetic stresses within the shocked pulsar wind affect the overall structure strongly, producing "whiskers", "tails", "filled-in" and "mushroom" shapes, as well as non-symmetric morphologies. On the other hand, the interstellar medium inhomogeneities and the anisotropy of the energy flux in the pulsar wind have only a mild impact of the PWN morphology. In a few cases, when we clearly identify specific morphological structures, our results do not favor alignment of the pulsar spin axis and proper velocity. Our calculations of the underlying emission processes explain the low synchrotron X-ray efficiency (in terms of the spin-down luminosity) and imply an energetical subdominant contribution of the inverse Compton process. △ Less

Submitted 22 January, 2019; v1 submitted 19 April, 2018; originally announced April 2018.

Comments: 27 pages, 27 figures, accepted by MNRAS

Synchrotron self-absorption in GRB afterglows: the effects of a thermal electron population

Authors: Donald C. Warren, Maxim V. Barkov, Hirotaka Ito, Shigehiro Nagataki, Tanmoy Laskar

Abstract: In the standard synchrotron afterglow model, a power law of electrons is responsible for all aspects of photon production and absorption. Recent numerical work has shown that the vast majority of particles in the downstream medium are actually "thermal" particles, which were shock-heated but did not enter the Fermi acceleration process (the name stands in contrast to the nonthermal high-energy tai… ▽ More In the standard synchrotron afterglow model, a power law of electrons is responsible for all aspects of photon production and absorption. Recent numerical work has shown that the vast majority of particles in the downstream medium are actually "thermal" particles, which were shock-heated but did not enter the Fermi acceleration process (the name stands in contrast to the nonthermal high-energy tail, rather than connoting a Maxwellian distribution). There are substantial differences at optical and higher energies when these thermal electrons participate in the afterglow, but early work along these lines ignored the radio end of the electromagnetic spectrum. We report here on an extension of previous Monte Carlo simulations of gamma-ray burst afterglows. The model now includes the synchrotron self-absorption (SSA) process and so can simulate afterglows across the entire EM spectrum, and several orders of magnitude in time. In keeping with earlier work, inclusion of the thermal electrons increases the SSA frequency by a factor of 30, and the radio intensity by a factor of 100. Furthermore, these changes happen with no modification to the late optical or X-ray afterglow. Our results provide very strong evidence that thermal electrons must be considered in any multiwavelength model for afterglows. △ Less

Submitted 16 April, 2018; originally announced April 2018.

Comments: Submitted to MNRAS; 11 pages, 4 figures

GRB170817A associated with GW170817: multifrequency observations and modeling of prompt gamma-ray emission

Authors: Alexei Pozanenko, Maxim V. Barkov, Pavel Yu. Minaev, Alina A. Volnova, Elena D. Mazaeva, Alexander S. Moskvitin, Maxim A. Krugov, Vladimir A. Samodurov, Vladimir M. Loznikov, Maxim Lyutikov

Abstract: We present our observations of electromagnetic transients associated with GW170817/GRB 170817A using optical telescopes of Chilescope observatory and Big Scanning Antenna (BSA) of Pushchino Radio Astronomy Observatory at 110MHz. The Chilescope observatory detected an optical transient of $\sim19^m$ on the third day in the outskirts of the galaxy NGC 4993; we continued observations following its ra… ▽ More We present our observations of electromagnetic transients associated with GW170817/GRB 170817A using optical telescopes of Chilescope observatory and Big Scanning Antenna (BSA) of Pushchino Radio Astronomy Observatory at 110MHz. The Chilescope observatory detected an optical transient of $\sim19^m$ on the third day in the outskirts of the galaxy NGC 4993; we continued observations following its rapid decrease. We put an upper limit of $1.5\times10^{4}$ Jy on any radio source with a duration of 10-60 s which may be associated with GW170817/GRB 170817A. The prompt gamma-ray emission consists of two distinctive components - a hard short pulse delayed by $\sim2$ seconds with respect to the LIGO signal and softer thermal pulse with $T\sim 10 $ keV lasting for another $\sim2$ seconds. The appearance of a thermal component at the end of the burst is unusual for short GRBs. Both the hard and the soft components do not satisfy the Amati relation, making GRB 170817A distinctively different from other short GRBs. Based on gamma-ray and optical observations, we develop a model for the prompt high-energy emission associated with GRB 170817A. The merger of two neutron stars creates an accretion torus of $\sim10^{-2} M_\odot$, which supplies the black hole with magnetic flux and confines the Blandford-Znajek-powered jet. We associate the hard prompt spike with the quasispherical breakout of the jet from the disk wind. As the jet plows through the wind with subrelativistic velocity, it creates a radiation-dominated shock that heats the wind material to tens of kiloelectron volts, producing the soft thermal component. △ Less

Submitted 22 January, 2018; v1 submitted 16 October, 2017; originally announced October 2017.

Comments: 20 pages, 12 figures

Journal ref: The Astrophysical Journal Letters, Volume 852, Issue 2, article id. L30, 18 pp. (2018)

HESS J0632+057: hydrodynamics and nonthermal emission

Authors: Valentí Bosch-Ramon, Maxim V. Barkov, Andrea Mignone, Pol Bordas

Abstract: HESS J0632+057 is an eccentric gamma-ray Be binary that produces non-thermal radio, X-rays, GeV, and very high-energy gamma rays. The non-thermal emission of HESS J0632+057 is modulated with the orbital period, with a dominant maximum before apastron passage. The nature of the compact object in HESS J0632+057 is not known, although it has been proposed to be a young pulsar as in PSR B1259-63, the… ▽ More HESS J0632+057 is an eccentric gamma-ray Be binary that produces non-thermal radio, X-rays, GeV, and very high-energy gamma rays. The non-thermal emission of HESS J0632+057 is modulated with the orbital period, with a dominant maximum before apastron passage. The nature of the compact object in HESS J0632+057 is not known, although it has been proposed to be a young pulsar as in PSR B1259-63, the only gamma-ray emitting high-mass binary known to host a non-accreting pulsar. In this Letter, we present hydrodynamical simulations of HESS J0632+057 in the context of a pulsar and a stellar wind interacting in an eccentric binary, and propose a scenario for the non-thermal phenomenology of the source. In this scenario, the non-thermal activity before and around apastron is linked to the accumulation of non-thermal particles in the vicinity of the binary, and the sudden drop of the emission before apastron is produced by the disruption of the two-wind interaction structure, allowing these particles to efficiently escape. In addition to providing a framework to explain the non-thermal phenomenology of the source, this scenario predicts extended, moving X-ray emitting structures similar to those observed in PSR B1259-63. △ Less

Submitted 31 July, 2017; originally announced August 2017.

Comments: 6 pages, 5 figures, accepted by MNRAS Letters

The impact of red giant/AGB winds on AGN jet propagation

Authors: Manel Perucho, Valentí Bosch-Ramon, Maxim V. Barkov

Abstract: Dense stellar winds may mass-load the jets of active galactic nuclei, although it is unclear what are the time and spatial scales in which the mixing takes place. We study the first steps of the interaction between jets and stellar winds, and also the scales at which the stellar wind may mix with the jet and mass-load it. We present a detailed two-dimensional simulation, including thermal cooling,… ▽ More Dense stellar winds may mass-load the jets of active galactic nuclei, although it is unclear what are the time and spatial scales in which the mixing takes place. We study the first steps of the interaction between jets and stellar winds, and also the scales at which the stellar wind may mix with the jet and mass-load it. We present a detailed two-dimensional simulation, including thermal cooling, of a bubble formed by the wind of a star. We also study the first interaction of the wind bubble with the jet using a three-dimensional simulation in which the star enters the jet. Stability analysis is carried out for the shocked wind structure, to evaluate the distances over which the jet-dragged wind, which forms a tail, can propagate without mixing with the jet flow. The two-dimensional simulations point at quick wind bubble expansion and fragmentation after about one bubble shock crossing time. Three-dimensional simulations and stability analysis point at local mixing in the case of strong perturbations and relatively small density ratios between the jet and the jet dragged-wind, and to a possibly more stable shocked wind structure at the phase of maximum tail mass flux. Analytical estimates also indicate that very early stages of the star jet-penetration time may be also relevant for mass loading. The combination of these and previous results from the literature suggest highly unstable interaction structures and efficient wind-jet flow mixing on the scale of the jet interaction height, possibly producing strong inhomogeneities within the jet. In addition, the initial wind bubble shocked by the jet leads to a transient, large interaction surface. The interaction structure can be a source of significant non-thermal emission. △ Less

Submitted 20 June, 2017; originally announced June 2017.

Comments: Accepted for publication in Astronomy & Astrophysics

Why does Einasto profile index $n\sim 6$ occur so frequently?

Authors: Anton N. Baushev, Maxim V. Barkov

Abstract: We consider the behavior of spherically symmetric Einasto halos composed of gravitating particles in the Fokker-Planck approximation. This approach allows us to consider the undesirable influence of close encounters in the N-body simulations more adequately than the generally accepted criteria. The Einasto profile with index $n \approx 6$ is a stationary solution of the Fokker-Planck equation in t… ▽ More We consider the behavior of spherically symmetric Einasto halos composed of gravitating particles in the Fokker-Planck approximation. This approach allows us to consider the undesirable influence of close encounters in the N-body simulations more adequately than the generally accepted criteria. The Einasto profile with index $n \approx 6$ is a stationary solution of the Fokker-Planck equation in the halo center. There are some reasons to believe that the solution is an attractor. Then the Fokker-Planck diffusion tends to transform a density profile to the equilibrium one with the Einasto index $n \approx 6$. We suggest this effect as a possible reason why the Einasto index $n \approx 6$ occurs so frequently in the interpretation of N-body simulation results. The results obtained cast doubt on generally accepted criteria of N-body simulation convergence. △ Less

Submitted 22 February, 2018; v1 submitted 15 May, 2017; originally announced May 2017.

Comments: 7 pages, 2 figures, Accepted to JCAP

Scenarios for ultrafast gamma-ray variability in AGN

Authors: F. A. Aharonian, M. V. Barkov, D. Khangulyan

Abstract: We analyze three scenarios to address the challenge of ultrafast gamma-ray variability reported from active galactic nuclei. We focus on the energy requirements imposed by these scenarios: (i) external cloud in the jet, (ii) relativistic blob propagating through the jet material, and (iii) production of high-energy gamma rays in the magnetosphere gaps. We show that while the first two scenarios ar… ▽ More We analyze three scenarios to address the challenge of ultrafast gamma-ray variability reported from active galactic nuclei. We focus on the energy requirements imposed by these scenarios: (i) external cloud in the jet, (ii) relativistic blob propagating through the jet material, and (iii) production of high-energy gamma rays in the magnetosphere gaps. We show that while the first two scenarios are not constrained by the flare luminosity, there is a robust upper limit on the luminosity of flares generated in the black hole magnetosphere. This limit depends weakly on the mass of the central black hole and is determined by the accretion disk magnetization, viewing angle, and the pair multiplicity. For the most favorable values of these parameters, the luminosity for 5-minute flares is limited by $2\times10^{43}\rm\,erg\,s^{-1}$, which excludes a black hole magnetosphere origin of the flare detected from IC310. In the scopes of scenarios (i) and (ii), the jet power, which is required to explain the IC310 flare, exceeds the jet power estimated based on the radio data. To resolve this discrepancy in the framework of the scenario (ii), it is sufficient to assume that the relativistic blobs are not distributed isotropically in the jet reference frame. A realization of scenario (i) demands that the jet power during the flare exceeds by a factor $10^2$ the power of the radio jet relevant to a timescale of $10^8$ years. △ Less

Submitted 16 May, 2017; v1 submitted 26 April, 2017; originally announced April 2017.

Comments: 15 pages, accepted by ApJ

Nonlinear particle acceleration and thermal particles in GRB afterglows

Authors: Donald C. Warren, Donald C. Ellison, Maxim V. Barkov, Shigehiro Nagataki

Abstract: The standard model for GRB afterglow emission treats the accelerated electron population as a simple power law, $N(E) \propto E^{-p}$ for $p \gtrsim 2$. However, in standard Fermi shock acceleration a substantial fraction of the swept-up particles do not enter the acceleration process at all. Additionally, if acceleration is efficient then the nonlinear backreaction of accelerated particles on the… ▽ More The standard model for GRB afterglow emission treats the accelerated electron population as a simple power law, $N(E) \propto E^{-p}$ for $p \gtrsim 2$. However, in standard Fermi shock acceleration a substantial fraction of the swept-up particles do not enter the acceleration process at all. Additionally, if acceleration is efficient then the nonlinear backreaction of accelerated particles on the shock structure modifies the shape of the non-thermal tail of the particle spectra. Both of these modifications to the standard synchrotron afterglow impact the luminosity, spectra, and temporal variation of the afterglow. To examine the effects of including thermal particles and nonlinear particle acceleration on afterglow emission, we follow a hydrodynamical model for an afterglow jet and simulate acceleration at numerous points during the evolution. When thermal particles are included, we find that the electron population is at no time well-fitted by a single power law, though the highest-energy electrons are; if the acceleration is efficient, then the power law region is even smaller. Our model predicts hard-soft-hard spectral evolution at X-ray energies, as well as an uncoupled X-ray and optical light curve. Additionally, we show that including emission from thermal particles has drastic effects (factors of 100 and 30, respectively) on the observed flux at optical and GeV energies. This enhancement of GeV emission makes afterglow detections by future $γ$-ray observatories, such as CTA, very likely. △ Less

Submitted 16 January, 2017; originally announced January 2017.

Comments: 18 pages, 9 figures

The effects of the stellar wind and orbital motion on the jets of high-mass microquasars

Authors: V. Bosch-Ramon, M. V. Barkov

Abstract: High-mass microquasar jets propagate under the effect of the wind from the companion star, and the orbital motion of the binary system. The stellar wind and the orbit may be dominant factors determining the jet properties beyond the binary scales.} An analytical study is performed to characterize the effects of the stellar wind and the orbital motion on the jet properties. Accounting for the wind… ▽ More High-mass microquasar jets propagate under the effect of the wind from the companion star, and the orbital motion of the binary system. The stellar wind and the orbit may be dominant factors determining the jet properties beyond the binary scales.} An analytical study is performed to characterize the effects of the stellar wind and the orbital motion on the jet properties. Accounting for the wind thrust transferred to the jet, we derive analytical estimates to characterize the jet evolution under the impact of the stellar wind. We include the Coriolis force effect, induced by orbital motion and enhanced by the wind presence. Large-scale evolution of the jet is sketched accounting for wind-to-jet thrust transfer, total energy conservation, and wind-jet flow mixing. If the angle of the wind-induced jet bending is larger than its half-opening angle, the following is expected: (i) a strong recollimation shock; (ii) bending against orbital motion, caused by Coriolis forces and enhanced by the wind presence; and (iii) non-ballistic helical propagation further away. Even if disrupted, the jet can reaccelerate due to ambient pressure gradients, but wind entrainment can weaken this acceleration. On large scales, the opening angle of the helical structure is determined by the wind-jet thrust relation, and the wind-loaded jet flow can be rather slow. The impact of stellar winds on high-mass microquasar jets can yield non-ballistic helical jet trajectories, jet partial disruption and wind mixing, shocks, and possibly non-thermal emission. Among several observational diagnostics at different bands, the radio morphology on milliarcsecond scales can be particularly insightful regarding the wind-jet interaction. △ Less

Submitted 21 April, 2016; originally announced April 2016.

Comments: 7 pages, 3 figures, accepted for publication in Astronomy & Astrophysics

Journal ref: A&A 590, A119 (2016)

Relativistic Tearing and Drift-kink Instabilities in Two-fluid Simulations

Authors: Maxim V. Barkov, Serguei S. Komissarov

Abstract: The stability of current sheets in collisionless relativistic pair plasma was studied via two-dimensional two-fluid relativistic magnetohydrodynamic simulations with vanishing internal friction between fluids. In particular, we investigated the linear growth of the tearing and drift-kink modes in the current sheets both with and without the guide field and obtained the growth rates which are very… ▽ More The stability of current sheets in collisionless relativistic pair plasma was studied via two-dimensional two-fluid relativistic magnetohydrodynamic simulations with vanishing internal friction between fluids. In particular, we investigated the linear growth of the tearing and drift-kink modes in the current sheets both with and without the guide field and obtained the growth rates which are very similar to what has been found in the corresponding PIC simulations. This suggests that the two-fluid simulations can be useful in studying the large-scale dynamics of astrophysical relativistic plasmas in problems involving magnetic reconnection. △ Less

Submitted 8 February, 2016; originally announced February 2016.

Comments: 11 pages, 7 figures, submitted to MNRAS

Photospheric Emission from Collapsar Jets in 3D Relativistic Hydrodynamics

Authors: Hirotaka Ito, Jin Matsumoto, Shigehiro Nagataki, Donald C. Warren, Maxim V. Barkov

Abstract: We explore the photospheric emission from a relativistic jet breaking out from a massive stellar envelope based on relativistic hydrodynamical simulations and post-process radiation transfer calculations in three dimensions. To investigate the impact of three-dimensional (3D) dynamics on the emission, two models of injection conditions are considered for the jet at the center of the progenitor sta… ▽ More We explore the photospheric emission from a relativistic jet breaking out from a massive stellar envelope based on relativistic hydrodynamical simulations and post-process radiation transfer calculations in three dimensions. To investigate the impact of three-dimensional (3D) dynamics on the emission, two models of injection conditions are considered for the jet at the center of the progenitor star: one with periodic precession and another without precession. We show that structures developed within the jet due to the interaction with the stellar envelope, as well as due to the precession, have a significant imprint on the resulting emission. Particularly, we find that the signature of precession activity by the central engine is not smeared out and can be directly observed in the light curve as a periodic signal. We also show non-thermal features, which can account for observations of gamma-ray bursts, are produced in the resulting spectra, even though only thermal photons are injected initially and the effect of non-thermal particles is not considered. △ Less

Submitted 27 November, 2015; v1 submitted 11 November, 2015; originally announced November 2015.

Comments: 7 pages, 4 figures, accepted for publication in ApJ Letters

The origin of the X-ray-emitting object moving away from PSR B1259-63

Authors: Maxim V. Barkov, Valenti Bosch-Ramon

Abstract: A mysterious X-ray-emitting object has been detected moving away from the high-mass gamma-ray binary PSR B1259-63, which contains a non-accreting pulsar and a Be star whose winds collide forming a complex interaction structure. Given the strong eccentricity of this binary, the interaction structure should be strongly anisotropic, which together with the complex evolution of the shocked winds, coul… ▽ More A mysterious X-ray-emitting object has been detected moving away from the high-mass gamma-ray binary PSR B1259-63, which contains a non-accreting pulsar and a Be star whose winds collide forming a complex interaction structure. Given the strong eccentricity of this binary, the interaction structure should be strongly anisotropic, which together with the complex evolution of the shocked winds, could explain the origin of the observed moving X-ray feature. We propose here that a fast outflow made of a pulsar-stellar wind mixture is always present moving away from the binary in the apastron direction, with the injection of stellar wind occurring at orbital phases close to periastron passage. This outflow periodically loaded with stellar wind would move with a high speed, and likely host non-thermal activity due to shocks, on scales similar to those of the observed moving X-ray object. Such an outflow is thus a very good candidate to explain this X-ray feature. This, if confirmed, would imply pulsar-to-stellar wind thrust ratios of $\sim 0.1$, and the presence of a jet-like structure on the larger scales, up to its termination in the ISM. △ Less

Submitted 26 October, 2015; originally announced October 2015.

Comments: 5 pages, 5 figures

Star-Jet Interactions and Gamma-Ray Outbursts from 3c454.3

Authors: M. V. Barkov, D. V. Khangulyan, V. Bosch-Ramon, F. A. Aharonian, A. V. Dorodnitsyn

Abstract: We propose a model to explain the ultra-bright GeV gamma-ray flares observed from the blazar \c454. The model is based on the concept of a relativistic jet interacting with compact gas condensations produced when a star (red giant) crosses the jet close to the central black hole. The study includes an analytical treatment of the evolution of the envelop lost by the star within the jet, and calcula… ▽ More We propose a model to explain the ultra-bright GeV gamma-ray flares observed from the blazar \c454. The model is based on the concept of a relativistic jet interacting with compact gas condensations produced when a star (red giant) crosses the jet close to the central black hole. The study includes an analytical treatment of the evolution of the envelop lost by the star within the jet, and calculations of the related high-energy radiation. The model readily explains the day-long, variable on timescales of hours, GeV gamma-ray flare from \c454, observed during November 2010 on top of a weeks-long plateau. In the proposed scenario, the plateau state is caused by a strong wind generated by the heating of the star atmosphere by nonthermal particles accelerated at the jet-star interaction region. The flare itself could be produced by a few clouds of matter lost by the red giant after the initial impact of the jet. In the framework of the proposed scenario, the observations constrain the key model parameters of the source, including the mass of the central black hole: $M_{\rm BH}\simeq 10^9 M_{\odot}$, the total jet power: $L_{\rm j}\simeq 10^{48}\,\rm erg\,s^{-1}$, and the Doppler factor of the gamma-ray emitting clouds, $δ\simeq 20$. Whereas we do not specify the particle acceleration mechanisms, the potential gamma-ray production processes are discussed and compared in the context of the proposed model. We argue that synchrotron radiation of protons has certain advantages compared to other radiation channels of directly accelerated electrons. △ Less

Submitted 9 February, 2015; originally announced February 2015.

Comments: 2014 Fermi Symposium proceedings - eConf C14102.1, 7 pages, 3 figures

Orbital evolution of colliding star and pulsar winds in 2D and 3D; effects of: dimensionality, EoS, resolution, and grid size

Authors: V. Bosch-Ramon, M. V. Barkov, M. Perucho

Abstract: (abridged)The structure formed by the shocked winds of a massive star and a non-accreting pulsar in a binary suffers periodic and random variations of orbital and non-linear dynamical origin. For the 1st time, we simulate in 3 D the interaction of isotropic stellar and relativistic pulsar winds along 1 full orbit, on scales well beyond the binary size. We also investigate the impact of grid resolu… ▽ More (abridged)The structure formed by the shocked winds of a massive star and a non-accreting pulsar in a binary suffers periodic and random variations of orbital and non-linear dynamical origin. For the 1st time, we simulate in 3 D the interaction of isotropic stellar and relativistic pulsar winds along 1 full orbit, on scales well beyond the binary size. We also investigate the impact of grid resolution and size, and of different EoOs: a gamma-constant ideal gas, and an ideal gas with gamma dependent on temperature. We carry out, with the code PLUTO, relativistic HD simulations in 2 and 3 D of the interaction of a slow wind and a relativistic wind with Gamma=2 along 1 full orbit up to ~100 x the binary size. The different 2-D simulations are carried out with equal and larger grid resolution and size, and 1 of them is done with a more realistic equation of state, than in 3 D. The simulations in 3 D confirm previous results in 2 D. The shocked flows are subject to a faster instabilities growth in 3 D, which enhances the presence of shocks, mixing, and large-scale disruption. In 2 D, higher resolution simulations confirm lower resolution results, simulations with larger grid sizes strengthen the case for the loss of global coherence of the shocked-wind structure, and simulations with 2 different EoOs yield very similar results. In addition to the KHI, we find that the Richtmyer-Meshkov and the RTI are likely acting together in the shocked flow evolution. Simulations in 3 D confirm that the interaction of stellar and pulsar winds yields structures that evolve non-linearly and get strongly entangled. The evolution is accompanied by strong kinetic energy dissipation, rapid changes in flow orientation and speed, and turbulent motion. The results strengthen the case for the loss of global coherence of the shocked structure on large scales, although higher pulsar wind speed simulations are needed. △ Less

Submitted 6 March, 2015; v1 submitted 28 November, 2014; originally announced November 2014.

Comments: 13 pages, 12 figures, accepted for publication in Astronomy and Astrophysics

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

Signatures of very massive stars: supercollapsars and their cosmological rate

Authors: Umberto Maio, Maxim V. Barkov

Abstract: We compute the rate of supercollapsars by using cosmological, N-body, hydro, chemistry simulations of structure formation, following detailed stellar evolution according to proper yields (for He, C, N, O, Si, S, Fe, Mg, Ca, Ne, etc.) and lifetimes for stars having different masses and metallicities, and for different stellar populations (population III and population II-I). We find that supercolla… ▽ More We compute the rate of supercollapsars by using cosmological, N-body, hydro, chemistry simulations of structure formation, following detailed stellar evolution according to proper yields (for He, C, N, O, Si, S, Fe, Mg, Ca, Ne, etc.) and lifetimes for stars having different masses and metallicities, and for different stellar populations (population III and population II-I). We find that supercollapsars are usually associated to dense, collapsing gas with little metal pollution and with abundances dominated by oxygen. The resulting supercollapsar rate is about $10^{-2}\,\rm yr^{-1} sr^{-1}$ at redshift $z=0$, and their contribution to the total rate is $ < 0.1 $ per cent, which explains why they have never been detected so far. Expected rates at redshift $z\simeq 6$ are of the order of $\sim 10^{-3}\,\rm yr^{-1} sr^{-1}$ and decrease further at higher $z$. Because of the strong metal enrichment by massive, short-lived stars, only $\sim 1$ supercollapsar generation is possible in the same star forming region. Given their sensitivity to the high-mass end of the primordial stellar mass function, they are suitable candidates to probe pristine population III star formation and stellar evolution at low metallicities. △ Less

Submitted 28 January, 2014; originally announced January 2014.

Comments: 6 pages; accepted, MNRAS. "Apri la mente a quel ch'io ti paleso" (Par. V, 40)

Formation of large-scale magnetic structures associated with the Fermi bubbles

Authors: Maxim V. Barkov, Valenti Bosch-Ramon

Abstract: The Fermi bubbles are part of a complex region of the Milky Way. This region presents broadband extended non-thermal radiation, apparently coming from a physical structure rooted in the Galactic Centre and with a partly-ordered magnetic field threading it. We explore the possibility of an explosive origin for the Fermi bubble region to explain its morphology, in particular that of the large-scale… ▽ More The Fermi bubbles are part of a complex region of the Milky Way. This region presents broadband extended non-thermal radiation, apparently coming from a physical structure rooted in the Galactic Centre and with a partly-ordered magnetic field threading it. We explore the possibility of an explosive origin for the Fermi bubble region to explain its morphology, in particular that of the large-scale magnetic fields, and provide context for the broadband non-thermal radiation. We perform 3D magnetohydrodynamical simulations of an explosion from a few million years ago that pushed and sheared a surrounding magnetic loop, anchored in the molecular torus around the Galactic Centre. Our results can explain the formation of the large-scale magnetic structure in the Fermi bubble region. Consecutive explosive events may match better the morphology of the region. Faster velocities at the top of the shocks than at their sides may explain the hardening with distance from the Galactic Plane found in the GeV emission. In the framework of our scenario, we estimate the lifetime of the Fermi bubbles as $2\times10^6$ yr, with a total energy injected in the explosion(s) $> 10^{55}$ ergs. The broadband non-thermal radiation from the region may be explained by leptonic emission, more extended in radio and X-rays, and confined to the Fermi bubbles in gamma rays. △ Less

Submitted 26 November, 2013; originally announced November 2013.

Comments: 5 pages, 4 figures, accepted for A&A

Star-Jet Interactions and Gamma-Ray Outbursts from 3C454.3

Authors: D. V. Khangulyan, M. V. Barkov, V. Bosch-Ramon, F. A. Aharonian, A. V. Dorodnitsyn

Abstract: We propose a model to explain the ultra-bright GeV gamma-ray flares observed from the blazar 3C454.3. The model is based on the concept of a relativistic jet interacting with compact gas condensations produced when a star (red giant) crosses the jet close to the central black hole. The study includes an analytical treatment of the evolution of the envelop lost by the star within the jet, and calcu… ▽ More We propose a model to explain the ultra-bright GeV gamma-ray flares observed from the blazar 3C454.3. The model is based on the concept of a relativistic jet interacting with compact gas condensations produced when a star (red giant) crosses the jet close to the central black hole. The study includes an analytical treatment of the evolution of the envelop lost by the star within the jet, and calculations of the related high-energy radiation. The model readily explains the day-long, variable on timescales of hours, GeV gamma-ray flare from 3C454.3, observed during November 2010 on top of a weeks-long plateau. In the proposed scenario, the plateau state is caused by a strong wind generated by the heating of the star atmosphere by nonthermal particles accelerated at the jet-star interaction region. The flare itself could be produced by a few clouds of matter lost by the red giant after the initial impact of the jet. In the framework of the proposed scenario, the observations constrain the key model parameters of the source, including the mass of the central black hole: $M_{\rm BH}\simeq 10^9 M_{\odot}$, the total jet power: $L_{\rm j}\simeq 10^{48}\,\rm erg\,s^{-1}$, and the Doppler factor of the gamma-ray emitting clouds, $δ\simeq 20$. Whereas we do not specify the particle acceleration mechanisms, the potential gamma-ray production processes are discussed and compared in the context of the proposed model. We argue that synchrotron radiation of protons has certain advantages compared to other radiation channels of directly accelerated electrons. △ Less

Submitted 27 May, 2013; v1 submitted 22 May, 2013; originally announced May 2013.

Comments: 16 pages, 5 figures, submitted to ApJ

Jets and gamma-ray emission from isolated accreting black holes

Authors: M. V. Barkov, D. V. Khangulyan, S. B. Popov

Abstract: The large number of isolated black holes (IBHs) in the Galaxy, estimated to be 10^8, implies a very high density of 10^-4 pc^-3 and an average distance between IBHs of 10 pc. Our study shows that the magnetic flux, accumulated on the horizon of an IBH because of accretion of interstellar matter, allows the Blandford-Znajeck mechanism to be activated. Thus, electron-positron jets can be launched. W… ▽ More The large number of isolated black holes (IBHs) in the Galaxy, estimated to be 10^8, implies a very high density of 10^-4 pc^-3 and an average distance between IBHs of 10 pc. Our study shows that the magnetic flux, accumulated on the horizon of an IBH because of accretion of interstellar matter, allows the Blandford-Znajeck mechanism to be activated. Thus, electron-positron jets can be launched. We have performed 2D numerical modelling which allowed the jet power to be estimated. Their inferred properties make such jets a feasible electron accelerator which, in molecular clouds, allows electron energy to be boosted up to 1 PeV. For the conditions expected in molecular clouds the radiative cooling time should be comparable to the escape time. Thus these sources can contribute both to the population of unidentified point-like sources and to the local cosmic ray (CR) electron spectrum. The impact of the generated electron CRs depends on the diffusion rate inside molecular clouds (MCs). If the diffusion regime in a MC is similar to Galactic diffusion, the produced electrons should rapidly escape the cloud and contribute to the Galactic CR population at very high energies >100 TeV. However, due to the modest jet luminosity (at the level of 10^35 erg s^-1) and low filling factor of MC, these sources cannot make a significant contribution to the spectrum of cosmic ray electrons at lower energies. On the other hand, if the diffusion within MCs operates at a rate close to the Bohm limit, the CR electrons escaping from the source should be confined in the cloud, significantly contributing to the local density of CRs. The IC emission of these locally-generated CRs may explain the variety of gamma ray spectra detected from nearby MCs. △ Less

Submitted 22 November, 2012; v1 submitted 3 September, 2012; originally announced September 2012.

Comments: 6 pages, accepted by MNRAS

Journal ref: Monthly Notices of the Royal Astronomical Society, Volume 427, Issue 1, pp. 589-594, (2012)

Simulations of stellar/pulsar wind interaction along one full orbit

Authors: V. Bosch-Ramon, M. V. Barkov, D. Khangulyan, M. Perucho

Abstract: The winds from a non-accreting pulsar and a massive star in a binary system collide forming a bow-shaped shock structure. The Coriolis force induced by orbital motion deflects the shocked flows, strongly affecting their dynamics. We study the evolution of the shocked stellar and pulsar winds on scales in which the orbital motion is important. Potential sites of non-thermal activity are investigate… ▽ More The winds from a non-accreting pulsar and a massive star in a binary system collide forming a bow-shaped shock structure. The Coriolis force induced by orbital motion deflects the shocked flows, strongly affecting their dynamics. We study the evolution of the shocked stellar and pulsar winds on scales in which the orbital motion is important. Potential sites of non-thermal activity are investigated. Relativistic hydrodynamical simulations in two dimensions, performed with the code PLUTO and using the adaptive mesh refinement technique, are used to model interacting stellar and pulsar winds on scales ~80 times the distance between the stars. The hydrodynamical results suggest the suitable locations of sites for particle acceleration and non-thermal emission. In addition to the shock formed towards the star, the shocked and unshocked components of the pulsar wind flowing away from the star terminate by means of additional strong shocks produced by the orbital motion. Strong instabilities lead to the development of turbulence and an effective two-wind mixing in both the leading and trailing sides of the interaction structure, which starts to merge with itself after one orbit. The adopted moderate pulsar-wind Lorentz factor already provides a good qualitative description of the phenomena involved in high-mass binaries with pulsars, and can capture important physical effects that would not appear in non-relativistic treatments. Simulations show that shocks, instabilities, and mass-loading yield efficient mass, momentum, and energy exchanges between the pulsar and the stellar winds. This renders a rapid increase in the entropy of the shocked structure, which will likely be disrupted on scales beyond the simulated ones. Several sites of particle acceleration and low- and high-energy emission can be identified. Doppler boosting will have significant and complex effects on radiation. △ Less

Submitted 17 July, 2012; v1 submitted 25 March, 2012; originally announced March 2012.

Comments: 8 pages, 11 figures, Astronomy and Astrophysics, in press, minor changes after acceptance

Interpretation of the flares of M87 at TeV energies in the cloud-jet interaction scenario

Authors: Maxim V. Barkov, Valenti Bosch-Ramon, Felix A. Aharonian

Abstract: Active galactic nuclei with misaligned jets have been recently established as a class of high-energy gamma-ray sources. M87, a nearby representative of this class, shows fast TeV variability on timescales less than one day. We present calculations performed in the framework of the scenario in which gamma-ray flares in non-blazar active galactic nuclei are produced by a red giant or a gas cloud int… ▽ More Active galactic nuclei with misaligned jets have been recently established as a class of high-energy gamma-ray sources. M87, a nearby representative of this class, shows fast TeV variability on timescales less than one day. We present calculations performed in the framework of the scenario in which gamma-ray flares in non-blazar active galactic nuclei are produced by a red giant or a gas cloud interacting with the jet. We show that both the light curve and energy spectrum of the spectacular April 2010 flare can be reproduced by this model, assuming that a relatively massive cloud of approx 1.e29 g penetrates into the jet at few tens of Schwarzschild radii from the super-massive black hole. △ Less

Submitted 27 June, 2012; v1 submitted 27 February, 2012; originally announced February 2012.

Comments: 8 pages, 8 figures, accepted by ApJ

Clouds and red giants interacting with the base of AGN jets

Authors: Valenti Bosch-Ramon, Manel Perucho, Maxim V. Barkov

Abstract: Extragalactic jets are formed close to supermassive black-holes in the center of galaxies. Large amounts of gas, dust, and stars cluster in the galaxy nucleus, and interactions between this ambient material and the jet base should be frequent, having dynamical as well as radiative consequences. This work studies the dynamical interaction of an obstacle, a clump of matter or the atmosphere of an ev… ▽ More Extragalactic jets are formed close to supermassive black-holes in the center of galaxies. Large amounts of gas, dust, and stars cluster in the galaxy nucleus, and interactions between this ambient material and the jet base should be frequent, having dynamical as well as radiative consequences. This work studies the dynamical interaction of an obstacle, a clump of matter or the atmosphere of an evolved star, with the innermost region of an extragalactic jet. Jet mass-loading and the high-energy outcome of this interaction are briefly discussed. Relativistic hydrodynamical simulations with axial symmetry have been carried out for homogeneous and inhomogeneous obstacles inside a relativistic jet. These obstacles may represent a medium inhomogeneity or the disrupted atmosphere of a red giant star. Once inside the jet, an homogeneous obstacle expands and gets disrupted after few dynamical timescales, whereas in the inhomogeneous case, a solid core can smoothen the process, with the obstacle mass-loss dominated by a dense and narrow tail pointing in the direction of the jet. In either case, matter is expected to accelerate and eventually get incorporated to the jet. Particles can be accelerated in the interaction region, and produce variable gamma-rays in the ambient matter, magnetic and photon fields. The presence of matter clumps or red giants into the base of an extragalactic jet likely implies significant jet mass-loading and slowing down. Fast flare-like gamma-ray events, and some level of persistent emission, are expected due to these interactions. △ Less

Submitted 25 January, 2012; originally announced January 2012.

Comments: 13 pages, 15 Figures, accepted for publication in Astronomy and Astrophysics

Direct Wind Accretion and Jet Launch in Binary Systems

Authors: Maxim V. Barkov, Dmitry V. Khangulyan

Abstract: In this paper we study the wind accretion onto a rotating black hole in the close binary system harboring a young massive star. It is shown that the angular momentum of the accreted stellar wind material is not sufficient for the formation of an accretion disk. On the other hand, in the considered conditions the Blanford-Znajek mechanism can be activated, thus powerful jets can be launched in the… ▽ More In this paper we study the wind accretion onto a rotating black hole in the close binary system harboring a young massive star. It is shown that the angular momentum of the accreted stellar wind material is not sufficient for the formation of an accretion disk. On the other hand, in the considered conditions the Blanford-Znajek mechanism can be activated, thus powerful jets can be launched in the direction of the rotation axis of the black hole. Importantly, no observational signatures of accretion, as typically seen from the thermal X-ray emission from the accretion disks, are expected in the suggested scenario. Here, properties of the generated jet are studied numerically in the framework of a 2D general relativity magnetohydrodynamical approach. Due to the accumulation of the magnetic flux at the black hole horizon, the jet power is expected to be modulated on a sub-second time-scale. Although the intervals between jet active phases depend on the magnetic flux escape from the black hole horizon (which can be modeled self-consistently only using a 3D code), a general estimate of the averaged jet power is obtained. It is expected that for the black hole rotation, expected in stellar binary system (the dimensionless rotation parameter a=0.5), approximately 10% of the accreted rest energy can be channeled into the jets. In the specific case of the gamma-ray binary system LS 5039, the obtained jet luminosity can be responsible for the observed GeV radiation if one invokes Doppler boosting, which can enhanced the apparent flux from the system. △ Less

Submitted 27 September, 2011; originally announced September 2011.

Comments: 8 pages, 4 figures

Large-scale flow dynamics and radiation in pulsar gamma-ray binaries

Authors: V. Bosch-Ramon, M. V. Barkov

Abstract: Several gamma-ray binaries show extended X-ray emission that may be associated to interactions of an outflow with the medium. Some of these systems are, or may be, high-mass binaries harboring young nonaccreting pulsars, in which the stellar and the pulsar winds collide, generating a powerful outflow that should terminate at some point in the ambient medium. This work studies the evolution and ter… ▽ More Several gamma-ray binaries show extended X-ray emission that may be associated to interactions of an outflow with the medium. Some of these systems are, or may be, high-mass binaries harboring young nonaccreting pulsars, in which the stellar and the pulsar winds collide, generating a powerful outflow that should terminate at some point in the ambient medium. This work studies the evolution and termination, as well as the related radiation, of the shocked-wind flow generated in high-mass binaries hosting powerful pulsars. A characterization, based on previous numerical work, is given for the stellar/pulsar wind interaction. Then, an analytical study of the further evolution of the shocked flow and its dynamical impact on the surrounding medium is carried out. Finally, the expected nonthermal emission from the flow termination shock, likely the dominant emitting region, is calculated. The shocked wind structure, initially strongly asymmetric, becomes a quasi-spherical, supersonically expanding bubble, with its energy coming from the pulsar and mass from the stellar wind. This bubble eventually interacts with the environment on ~pc scales, producing a reverse and, sometimes, a forward shock. Nonthermal leptonic radiation can be efficient in the reverse shock. Radio emission is expected to be faint, whereas X-rays can easily reach detectable fluxes. Under very low magnetic fields and large nonthermal luminosities, gamma rays may also be significant. We conclude that the complexity of the stellar/pulsar wind interaction is likely to be smoothed out outside the binary system, where the wind-mixed flow accelerates and eventually terminates in a strong reverse shock. This shock may be behind the extended X-rays observed in some binary systems. For very powerful pulsars, part of the unshocked pulsar wind may directly interact with the large-scale environment. △ Less

Submitted 14 October, 2011; v1 submitted 31 May, 2011; originally announced May 2011.

Comments: 11 pages, 6 Figures, 1 Table, Astronomy and Astrophysics, in press (small corrections after proofs)

Model of extended emission of short Gamma-ray Bursts

Authors: Maxim V. Barkov, Alexei S. Pozanenko

Abstract: Until now the existence of extended emission is an intriguing property of short bursts. It is not clear what is the nature of the extended emission. It might be a rising x-ray afterglow, or it could be a manifestation of the prolonged activity of a central engine. We consider short duration gamma-ray bursts, emphasizing the common properties of short bursts and short burst with extended emission.… ▽ More Until now the existence of extended emission is an intriguing property of short bursts. It is not clear what is the nature of the extended emission. It might be a rising x-ray afterglow, or it could be a manifestation of the prolonged activity of a central engine. We consider short duration gamma-ray bursts, emphasizing the common properties of short bursts and short burst with extended emission. Assuming that the extended emission with broad dynamic range is a common property of short bursts, we propose a two jet model which can describe both short main episode of hard spectra emission, specific for short bursts, and softer spectra extended emission by different off axis position of observer. The toy model involves a short duration jet powered by heating due to $ν\tildeν$ annihilation and long-lived Blandford-Znajek jet with significantly narrow opening angle. Our proposed model is a plausible mechanism for short duration burst energization, and can explain both short burst with and without extended emission within a single class of progenitor. △ Less

Submitted 6 June, 2011; v1 submitted 22 March, 2011; originally announced March 2011.

Comments: 6 pages, 4 figures

Recycling of Neutron Stars in Common Envelopes and Hypernova Explosions

Authors: Maxim V. Barkov, Serguei S. Komissarov

Abstract: In this paper we propose a new plausible mechanism of supernova explosions specific to close binary systems. The starting point is the common envelope phase in the evolution of a binary consisting of a red super giant and a neutron star. As the neutron star spirals towards the center of its companion it spins up via disk accretion. Depending on the specific angular momentum of gas captured by the… ▽ More In this paper we propose a new plausible mechanism of supernova explosions specific to close binary systems. The starting point is the common envelope phase in the evolution of a binary consisting of a red super giant and a neutron star. As the neutron star spirals towards the center of its companion it spins up via disk accretion. Depending on the specific angular momentum of gas captured by the neutron star via the Bondi-Hoyle mechanism, it may reach millisecond periods either when it is still inside the common envelope or after it has merged with the companion core. The high accretion rate may result in strong differential rotation of the neutron star and generation of the magnetar-strength magnetic field. The magnetar wind can blow away the common envelope if its magnetic field is as strong as $10^{15}\,$G, and can destroy the entire companion if it is as strong as $10^{16}\,$G. The total explosion energy can be comparable to the rotational energy of a millisecond pulsar and reach $10^{52}\,$erg. However, only a small amount of $^{56}$Ni is expected to be produced this way. The result is an unusual type-II supernova with very high luminosity during the plateau phase, followed by a sharp drop in brightness and a steep light-curve tail. The remnant is either a solitary magnetar or a close binary involving a Wolf-Rayet star and a magnetar. When this Wolf-Rayet star explodes this will be a third supernovae explosion in the same binary. △ Less

Submitted 29 March, 2011; v1 submitted 21 December, 2010; originally announced December 2010.

Comments: 16 pages, 4 figures

Rapid TeV variability in Blazars as result of Jet-Star Interaction

Authors: Maxim V. Barkov, Felix A. Aharonian, Sergey V. Bogovalov, Stanislav R. Kelner, Dmitriy V. Khangulyan

Abstract: We propose a new model for the description of ultra-short flares from TeV blazars by compact magnetized condensations (blobs), produced when red giant stars cross the jet close to the central black hole. Our study includes a simple dynamical model for the evolution of the envelope lost by the star in the jet, and its high energy nonthermal emission through different leptonic and hadronic radiation… ▽ More We propose a new model for the description of ultra-short flares from TeV blazars by compact magnetized condensations (blobs), produced when red giant stars cross the jet close to the central black hole. Our study includes a simple dynamical model for the evolution of the envelope lost by the star in the jet, and its high energy nonthermal emission through different leptonic and hadronic radiation mechanisms. We show that the fragmented envelope of the star can be accelerated to Lorentz factors up to 100 and radiate effectively the available energy in gamma-rays predominantly through proton synchrotron radiation or external inverse Compton scattering of electrons. The model can readily explain the minute-scale TeV flares on top of longer (typical time-scales of days) gamma-ray variability as observed from the blazar PKS 2155-304. In the framework of the proposed scenario, the key parameters of the source are robustly constrained. In the case of proton synchrotron origin of the emission a mass of the central black hole of $M_{\rm BH}\approx 10^8 M_{\odot}$, a total jet power of $L_{\rm j} \approx 2\times 10^{47} \, \rm erg\,s^{-1}$ and a Doppler factor, of the gamma-ray emitting blobs, of $δ\geq 40$ are required. Whilst for the external inverse Compton model, parameters of $M_{\rm BH}\approx 10^8 M_{\odot}$, $L_{\rm j} \approx 10^{46} \, \rm erg\,s^{-1}$ and the $δ\geq 150$ are required. △ Less

Submitted 23 December, 2011; v1 submitted 8 December, 2010; originally announced December 2010.

Comments: 25 pages, 11 figures, Submitted to ApJ

Gamma-ray flares from red giant/jet interactions in AGN

Authors: Maxim V. Barkov, Felix A. Aharonian, Valenti Bosch-Ramon

Abstract: Non-blazar AGN have been recently established as a class of gamma-ray sources. M87, a nearby representative of this class, show fast TeV variability on timescales of a few days. We suggest a scenario of flare gamma-ray emission in non-blazar AGN based on a red giant interacting with the jet at the base. We solve the hydrodynamical equations that describe the evolution of the envelope of a red gian… ▽ More Non-blazar AGN have been recently established as a class of gamma-ray sources. M87, a nearby representative of this class, show fast TeV variability on timescales of a few days. We suggest a scenario of flare gamma-ray emission in non-blazar AGN based on a red giant interacting with the jet at the base. We solve the hydrodynamical equations that describe the evolution of the envelope of a red giant blown by the impact of the jet. If the red giant is at least slightly tidally disrupted by the supermassive black hole, enough stellar material will be blown by the jet, expanding quickly until a significant part of the jet is shocked. This process can render suitable conditions for energy dissipation and proton acceleration, which could explain the detected day-scale TeV flares from M87 via proton-proton collisions. Since the produced radiation would be unbeamed, such an events should be mostly detected from non-blazar AGN. They may be frequent phenomena, detectable in the GeV-TeV range even up to distances of $\sim 1$ Gpc for the most powerful jets. The counterparts at lower energies are expected to be not too bright.} {M87, and nearby non-blazar AGN in general, can be fast variable sources of gamma-rays through red giant/jet interactions. △ Less

Submitted 3 September, 2010; v1 submitted 28 May, 2010; originally announced May 2010.

Comments: 8 pages, 4 figures

Journal ref: Astrophys.J.724:1517-1523,2010

Hard X-ray Bursts from Collapse of the Super Massive Stars

Authors: Maxim V. Barkov

Abstract: The very first stars in the Universe can be very massive, up to $10^6M_\odot$. They would leave behind massive black holes that could act as seeds for growing super massive black holes of active galactic nuclei. Given the anticipated fast rotation such stars would end their live as super massive collapsars and drive powerful magnetically-dominated jets. In this paper we investigate the possibility… ▽ More The very first stars in the Universe can be very massive, up to $10^6M_\odot$. They would leave behind massive black holes that could act as seeds for growing super massive black holes of active galactic nuclei. Given the anticipated fast rotation such stars would end their live as super massive collapsars and drive powerful magnetically-dominated jets. In this paper we investigate the possibility of observing the bursts of high-energy emission similar to the Long Gamma Ray Bursts associated with normal collapsars. We show that during the collapse of supercollapsars, the Blandford-Znajek mechanism can produce jets as powerful as few$\times10^{51}$erg/s and release up to $10^{56}$erg of the black hole rotational energy. Due to the higher intrinsic time scale and higher redshift the initial bright phase of the burst can last for about $10^5$ seconds whereas the central engine would remain active for about 10 days. Due to the high redshift the burst spectrum is expected to be soft, with the spectral energy distribution peaking at around 60keV. The peak total flux density is relatively low, few$\times 10^{-7}erg\, cm^{-2} s^{-1}$, but not prohibitive. The such events should be rear 0.03 year$^{-1}$, the observations needs long term program and could be done in future. △ Less

Submitted 23 March, 2010; originally announced March 2010.

Comments: 4 pages, 2 figures. proceedings of workshop "Many faces of GRB phenomena - optics vs high energy", SAO Russia, October 12-16, 2009,

Supercollapsars and their X-ray Bursts

Authors: S. S. Komissarov, M. V. Barkov

Abstract: The very first stars in the Universe can be very massive, up to $10^3M_\odot$. If born in large numbers, such massive stars can have strong impact on the subsequent star formation, producing strong ionising radiation and contaminating the primordial gas with heavy elements. They would leave behind massive black holes that could act as seeds for growing supermassive black holes of active galactic… ▽ More The very first stars in the Universe can be very massive, up to $10^3M_\odot$. If born in large numbers, such massive stars can have strong impact on the subsequent star formation, producing strong ionising radiation and contaminating the primordial gas with heavy elements. They would leave behind massive black holes that could act as seeds for growing supermassive black holes of active galactic nuclei. Given the anticipated fast rotation, such stars would end their live as supermassive collapsars and drive powerful magnetically-dominated jets. In this letter we investigate the possibility of observing the bursts of high-energy emission similar to the Long Gamma Ray Bursts associated with normal collapsars. We show that during the collapse of supercollapsars, the Blandford-Znajek mechanism can produce jets as powerful as few$\times10^{52}$erg/s and release up to $10^{54}$erg of the black hole rotational energy. Due to the higher intrinsic time scale and higher redshift the initial bright phase of the burst can last for about $10^4$ seconds, whereas the central engine would remain active for about one day. Due to the high redshift the burst spectrum is expected to be soft, with the spectral energy distribution peaking at around 20-30keV. The peak total flux density is relatively low, $10^{-7}{erg} {cm}^{-2} {s}^{-1}$, but not prohibitive. If one supercollapsar is produced per every mini-halo of dark matter arising from the 3-$σ$ cosmological fluctuations then the whole sky frequency of such bursts could reach several tens per year. △ Less

Submitted 25 November, 2009; v1 submitted 4 September, 2009; originally announced September 2009.

Comments: Accepted for publication in MNRAS, New calculations based on the Bethe's model of progenitor are given

Close Binary Progenitors of Long Gamma Ray Bursts

Authors: M. V. Barkov, S. S. Komissarov

Abstract: The strong dependence of the neutrino annihilation mechanism on the mass accretion rate makes it difficult to explain the LGRBs with duration in excess of 100 seconds as well as the precursors separated from the main gamma-ray pulse by few hundreds of seconds. Even more difficult is to explain the Swift observations of the shallow decay phase and X-ray flares, if they indeed indicate activity of… ▽ More The strong dependence of the neutrino annihilation mechanism on the mass accretion rate makes it difficult to explain the LGRBs with duration in excess of 100 seconds as well as the precursors separated from the main gamma-ray pulse by few hundreds of seconds. Even more difficult is to explain the Swift observations of the shallow decay phase and X-ray flares, if they indeed indicate activity of the central engine for as long as 10,000 seconds. These data suggest that some other, most likely magnetic mechanisms have to be considered. The magnetic models do not require the development of accretion disk within the first few seconds of the stellar collapse and hence do not require very rapidly rotating stellar cores at the pre-supernova state. This widens the range of potential LGRB progenitors. In this paper, we re-examine the close binary scenario allowing for the possibility of late development of accretion disks in the collapsar model and investigate the available range of mass accretion rates, black hole masses, and spins. A particularly interesting version of the binary progenitor involves merger of a WR star with an ultra-compact companion, neutron star or black hole. In this case we expect the formation of very long-lived accretion disks, that may explain the phase of shallow decay and X-ray flares observed by Swift. Similarly long-lived magnetic central engines are expected in the current single star models of LGRB progenitors due to their assumed exceptionally fast rotation. △ Less

Submitted 5 August, 2009; originally announced August 2009.

Comments: Submitted to MNRAS

Accretion of a massive magnetized torus on a rotating black hole

Authors: Maxim V. Barkov, Anton N. Baushev

Abstract: We present numerical simulations of the axisymmetric accretion of a massive magnetized plasma torus on a rotating black hole. We use a realistic equation of state, which takes into account neutrino cooling and energy loss due to nucleus dissociations. We simulated various magnetic field configurations and torus models, both optically thick and thin for neutrinos. It is shown that the neutrino co… ▽ More We present numerical simulations of the axisymmetric accretion of a massive magnetized plasma torus on a rotating black hole. We use a realistic equation of state, which takes into account neutrino cooling and energy loss due to nucleus dissociations. We simulated various magnetic field configurations and torus models, both optically thick and thin for neutrinos. It is shown that the neutrino cooling does not significantly change either the structure of the accretion flow or the total energy release of the system. The calculations evidence heating of the wind surrounding the collapsar by the shock waves generated at the jet-wind border. This mechanism can give rise to a hot corona around the binary system like SS433. Angular momentum of the accreting matter defines the time scale of the accretion. Due to the absence of the magnetic dynamo in our calculations, the initial strength and topology of the magnetic field determines magnetization of the black hole, jet formation properties and the total energy yield. We estimated the total energy transformed to jets as $1.3\times 10^{52}$ {ergs} which was sufficient to explain hypernova explosions like GRB 980425 or GRB 030329. △ Less

Submitted 27 May, 2009; originally announced May 2009.

Comments: 11 pages, 9 figures, submitted to MNRAS

Journal ref: New Astron.16:46-56,2011

Activation of the Blandford-Znajek mechanism in collapsing stars

Authors: S. S. Komissarov, M. V. Barkov

Abstract: Collapse of massive stars may result in formation of accreting black holes in their interior. The accreting stellar matter may advect substantial magnetic flux onto the black hole and promote release of its rotational energy via magnetic stresses (the Blandford-Znajek mechanism). In this paper we explore whether this process can explain the stellar explosions and relativistic jets associated wit… ▽ More Collapse of massive stars may result in formation of accreting black holes in their interior. The accreting stellar matter may advect substantial magnetic flux onto the black hole and promote release of its rotational energy via magnetic stresses (the Blandford-Znajek mechanism). In this paper we explore whether this process can explain the stellar explosions and relativistic jets associated with long Gamma-ray-bursts. In particularly, we show that the Blandford-Znajek mechanism is activated when the rest mass-energy density of matter drops below the energy density of magnetic field in the very vicinity of the black hole (within its ergosphere). We also discuss whether such a strong magnetic field is in conflict with the rapid rotation of stellar core required in the collapsar model and suggest that the conflict can be avoided if the progenitor star is a component of close binary. In this case the stellar rotation can be sustained via spin-orbital interaction. In an alternative scenario the magnetic field is generated in the accretion disk but in this case the magnetic flux through the black hole ergosphere is not expected to be sufficiently high to explain the energetics of hypernovae by the BZ mechanism alone. However, this energy deficit can be recovered via additional power provided by the disk. △ Less

Submitted 5 May, 2009; v1 submitted 17 February, 2009; originally announced February 2009.

Comments: submitted to MNRAS, new references and minor modifications in version 2

Central engines of Gamma Ray Bursts. Magnetic mechanism in the collapsar model

Authors: Maxim V. Barkov, Serguei S. Komissarov

Abstract: In this study we explore the magnetic mechanism of hypernovae and relativistic jets of long duration gamma ray bursts within the collapsar scenario. This is an extension of our earlier work [1]. We track the collapse of massive rotating stars onto a rotating central black hole using axisymmetric general relativistic magnetohydrodynamic code that utilizes a realistic equation of state and takes i… ▽ More In this study we explore the magnetic mechanism of hypernovae and relativistic jets of long duration gamma ray bursts within the collapsar scenario. This is an extension of our earlier work [1]. We track the collapse of massive rotating stars onto a rotating central black hole using axisymmetric general relativistic magnetohydrodynamic code that utilizes a realistic equation of state and takes into account the cooling associated with emission of neutrinos and the energy losses due to dissociation of nuclei. The neutrino heating is not included. We describe solutions with different black hole rotation, mass accretion rate, and strength of progenitor's magnetic field. Some of them exhibits strong explosions driven by Poynting-dominated jets with power up to $12\times10^{51} {erg s}^{-1}$. These jets originate from the black hole and powered via the Blandford-Znajek mechanism. A provisional criterion for explosion is derived. A number of simulation movies can be downloaded from http://www.maths.leeds.ac.uk/~serguei/research/movies/anim.html △ Less

Submitted 8 September, 2008; originally announced September 2008.

Comments: 4 pages, 3 figures, GAMMA08 Heidelberg Symposium

Hyper-accreting tori of Gamma Ray Bursters

Authors: M. V. Barkov

Abstract: We present numerical simulations of axisymmetric magnetised massive tori around rotating black holes taking into account the energy losses due to emission of neutrinos. A realistic equation of state is used which takes into account the energy losses due to dissociation of nuclei. The heating due to neutrino-antineutrino annihilation is not included. We study the cases of optically thick, semi-tr… ▽ More We present numerical simulations of axisymmetric magnetised massive tori around rotating black holes taking into account the energy losses due to emission of neutrinos. A realistic equation of state is used which takes into account the energy losses due to dissociation of nuclei. The heating due to neutrino-antineutrino annihilation is not included. We study the cases of optically thick, semi-transparent, and optically thin to neutrino discs. We show that neutrino cooling does not change significantly the structure of accretion flow and the total energy release. The time scale of accretion is set by the torus angular momentum. Due to the lack of magnetic dynamo in our calculations, it is the initial strength of magnetic field and its topology that determine the process of jet formation and its energetics. Extrapolation of our results gives the total energy released in the jet $\sim 10^{52}$ erg. This is sufficient to explain the hypernovae explosions associated with GRB 980425 and GRB 030329. △ Less

Submitted 4 May, 2008; originally announced May 2008.

Comments: 6 pages, 2 figures, conference "Cool discs, hot flows", Funasdalen, Sweden, 2008

Journal ref: AIP Conf.Proc.1054:79-85,2008

Magnetic acceleration of ultra-relativistic GRB and AGN jets

Authors: Maxim V. Barkov, Serguei S. Komissarov

Abstract: We present numerical simulations of cold, axisymmetric, magnetically driven relativistic outflows. The outflows are initially sub-Alfvénic and Poynting flux-dominated, with total--to--rest-mass energy flux ratio up to $μ\sim 620$. To study the magnetic acceleration of jets we simulate flows confined within a funnel with rigid wall of prescribed shape, which we take to be $z\propto r^a$ (in cylin… ▽ More We present numerical simulations of cold, axisymmetric, magnetically driven relativistic outflows. The outflows are initially sub-Alfvénic and Poynting flux-dominated, with total--to--rest-mass energy flux ratio up to $μ\sim 620$. To study the magnetic acceleration of jets we simulate flows confined within a funnel with rigid wall of prescribed shape, which we take to be $z\propto r^a$ (in cylindrical coordinates, with $a$ ranging from 1 to 2). This allows us to eliminate the numerical dissipative effects induced by a free boundary with an ambient medium. We find that in all cases they converge to a steady state characterized by a spatially extended acceleration region. For the jet solutions the acceleration process is very efficient - on the outermost scale of the simulation more than half of the Poynting flux has been converted into kinetic energy flux, and the terminal Lorentz factor approached its maximum possible value ($Γ_\infty \simeq μ$). The acceleration is accompanied by the collimation of magnetic field lines in excess of that dictated by the funnel shape. The numerical solutions are generally consistent with the semi-analytic self-similar jets solutions and the spatially extended acceleration observed in some astrophysical relativistic jets. In agreement with previous studies we also find that the acceleration is significantly less effective for wind solutions suggesting that pulsar winds may remain Poynting dominated when they reach the termination shock. △ Less

Submitted 16 April, 2008; v1 submitted 31 January, 2008; originally announced January 2008.

Comments: 4 pages, 3 figures, HEPRO-2007 Dublin

Stellar explosions powered by the Blandford-Znajek mechanism

Authors: M. V. Barkov, S. S. Komissarov

Abstract: In this letter we briefly describe the first results of our numerical study on the possibility of magnetic origin of relativistic jets of long duration gamma ray bursters within the collapsar scenario. We track the collapse of massive rotating stars onto a rotating central black hole using axisymmetric general relativistic magnetohydrodynamic code that utilizes a realistic equation of state of s… ▽ More In this letter we briefly describe the first results of our numerical study on the possibility of magnetic origin of relativistic jets of long duration gamma ray bursters within the collapsar scenario. We track the collapse of massive rotating stars onto a rotating central black hole using axisymmetric general relativistic magnetohydrodynamic code that utilizes a realistic equation of state of stellar matter, takes into account the cooling associated with emission of neutrinos, and the energy losses due to dissociation of nuclei. The neutrino heating is not included. We describe the solution for one particular model where the progenitor star has magnetic field $B=10^{10}$G. The solution exhibits strong explosion driven by the Poynting-dominated jets whose power exceeds $2\times10^{51} {erg/s}$. The jets originate mainly from the black hole and they are powered via the Blandford-Znajek mechanism. The full details of the simulations together with the results of parameter study will be presented elsewhere. A number of simulation movies can be downloaded from http://www.maths.leeds.ac.uk/~serguei/research/movies/anim.html △ Less

Submitted 9 January, 2008; v1 submitted 14 October, 2007; originally announced October 2007.

Comments: minor revision, accepted by MNRAS Letters, simulation movies can be downloaded from http://www.maths.leeds.ac.uk/~serguei/research/movies/anim.html

Magnetar-energized supernova explosions and GRB-jets

Authors: S. S. Komissarov, M. V. Barkov

Abstract: In this paper we report on the early evolution of core-collapse supernova explosion following the birth of a magnetar with the dipolar magnetic field of B=10^{15}G and the rotational period of 2ms, which was studied by means of axisymmetric general relativistic MHD simulations. The numerical models exhibit highly collimated magnetically-driven jets very early on. The jets are super-Alfvenic but… ▽ More In this paper we report on the early evolution of core-collapse supernova explosion following the birth of a magnetar with the dipolar magnetic field of B=10^{15}G and the rotational period of 2ms, which was studied by means of axisymmetric general relativistic MHD simulations. The numerical models exhibit highly collimated magnetically-driven jets very early on. The jets are super-Alfvenic but remain sub-fast until the end of the simulations (t=0.2s). The power released in the jets is about 3x10^{50}erg/s which implies the spin-down time of ~37s. The total rotational energy of the magnetar, E~10^{52}erg, is sufficient to drive hypernova but it is not clear as to how large a fraction of this energy can be transfered to the stellar envelope. Given the observed propagation speed of the jets, v_p~0.17c, they are expected to traverse the progenitor in few seconds and after this most of the released rotational energy would be simply carried away by these jets into the surrounding space. Our results provide the first more or less self-consistent numerical model of a central engine capable of producing, in the supernova setting and on a long-term basis, collimated jets with sufficient power to explain long duration GRBs and their afterglows. Although the flow speed of our jets is relatively low, v_j~0.5c$, the cooling of proto-neutron star will eventually result in much higher magnetization of its magnetospheres and ultra-relativistic asymptotic speeds of the jets. Given the relatively long cooling time-scale we still expect the jets to be only weakly relativistic by the time of break out. This leads to a model of GRB jets with systematic longitudinal variation of Lorentz factor which may have specific observational signatures both in the prompt and the afterglow emission. △ Less

Submitted 22 September, 2007; v1 submitted 2 July, 2007; originally announced July 2007.

Comments: accepted by MNRAS, additional simulations, expanded discussion