-
Unveiling the white dwarf in J191213.72-441045.1 through ultraviolet observations
Authors:
Ingrid Pelisoli,
Snehalata Sahu,
Maxim Lyutikov,
Maxim Barkov,
Boris T. Gaensicke,
Jaco Brink,
David A. H. Buckley,
Stephen B. Potter,
Axel Schwope,
S. H. Ramirez
Abstract:
J191213.72-441045.1 is a binary system composed of a white dwarf and an M-dwarf in a 4.03-hour orbit. It shows emission in radio, optical, and X-ray, all modulated at the white dwarf spin period of 5.3 min, as well as various orbital sideband frequencies. Like in the prototype of the class of radio-pulsing white dwarfs, AR Scorpii, the observed pulsed emission seems to be driven by the binary inte…
▽ More
J191213.72-441045.1 is a binary system composed of a white dwarf and an M-dwarf in a 4.03-hour orbit. It shows emission in radio, optical, and X-ray, all modulated at the white dwarf spin period of 5.3 min, as well as various orbital sideband frequencies. Like in the prototype of the class of radio-pulsing white dwarfs, AR Scorpii, the observed pulsed emission seems to be driven by the binary interaction. In this work, we present an analysis of far-ultraviolet spectra obtained with the Cosmic Origins Spectrograph at the Hubble Space Telescope, in which we directly detect the white dwarf in J191213.72-441045.1. We find that the white dwarf has an effective temperature of 11485+/-90 K and mass of 0.59+/-0.05 solar masses. We place a tentative upper limit on the magnetic field of ~50 MG. If the white dwarf is in thermal equilibrium, its physical parameters would imply that crystallisation has not started in the core of the white dwarf. Alternatively, the effective temperature could have been affected by compressional heating, indicating a past phase of accretion. The relatively low upper limit to the magnetic field and potential lack of crystallisation that could generate a strong field pose challenges to pulsar-like models for the system and give preference to propeller models with a low magnetic field. We also develop a geometric model of the binary interaction which explains many salient features of the system.
△ Less
Submitted 9 November, 2023;
originally announced November 2023.
-
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.
-
Relativistic coronal mass ejections from magnetars
Authors:
Praveen Sharma,
Maxim Barkov,
Maxim Lyutikov
Abstract:
We study dynamics of relativistic Coronal Mass Ejections (CMEs), from launching by shearing of foot-points (either slowly - the ``Solar flare'' paradigm, or suddenly - the ``star quake" paradigm), to propagation in the preceding magnetar wind. For slow shear, most of the energy injected into the CME is first spent on the work done on breaking through the over-laying magnetic field. At later stages…
▽ More
We study dynamics of relativistic Coronal Mass Ejections (CMEs), from launching by shearing of foot-points (either slowly - the ``Solar flare'' paradigm, or suddenly - the ``star quake" paradigm), to propagation in the preceding magnetar wind. For slow shear, most of the energy injected into the CME is first spent on the work done on breaking through the over-laying magnetic field. At later stages, sufficiently powerful CMEs may experience ``detonation" and lead to opening of the magnetosphere beyond some equipartition radius $r_{eq}$, where the energy of the CME becomes larger than the decreasing external magnetospheric energy. Post-CME magnetosphere relaxes via formation of a plasmoid-mediated current sheet, initially at $\sim r_{eq}$ and slowly reaching the light cylinder (this transient stage has much higher spindown rate and may produce an ``anti-glitch''). Both the location of the foot-point shear and the global magnetospheric configuration affect the frequent-and-weak versus rare-and-powerful CME dichotomy - to produce powerful flares the slow shear should be limited to field lines that close near the star. After the creation of a topologically disconnected flux tube, the tube quickly (at $\sim$ the light cylinder) comes into force-balance with the preceding wind, and is passively advected/frozen in the wind afterward. For fast shear (a local rotational glitch), the resulting large amplitude Alfven waves lead to opening of the magnetosphere (which later recovers similarly to the slow shear case). At distances much larger than the light cylinder, the resulting shear Alfven waves propagate through the wind non-dissipatively. Implications to Fast Radio Bursts are discussed.
△ Less
Submitted 17 February, 2023;
originally announced February 2023.
-
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.
-
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.
-
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.
-
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.
-
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.
-
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.
-
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.
-
Radio afterglow of magnetars' giant flares
Authors:
Riddhi Mehta,
Maxim Barkov,
Maxim Lyutikov
Abstract:
We develop a model for the radio afterglow of the giant flare of SGR 1806-20 arising due to the interaction of magnetically-dominated cloud, an analogue of Solar Coronal Mass Ejections (CMEs), with the interstellar medium (ISM). The CME is modeled as a spheromak-like configuration. The CME is first advected with the magnetar's wind and later interacts with the ISM, creating a strong forward shock…
▽ More
We develop a model for the radio afterglow of the giant flare of SGR 1806-20 arising due to the interaction of magnetically-dominated cloud, an analogue of Solar Coronal Mass Ejections (CMEs), with the interstellar medium (ISM). The CME is modeled as a spheromak-like configuration. The CME is first advected with the magnetar's wind and later interacts with the ISM, creating a strong forward shock and complicated backwards exhaust flow. Using three-dimensional magnetohydrodynamic simulations, we study various relative configurations of the magnetic field of the CME with respect to the ISM's magnetic field. We show that the dynamics of the forward shock mostly follows the Sedov-Taylor blastwave, while the internal structure of the shocked medium is considerably modified by the back flow, creating a multiple shock configuration. We calculate synthetic synchrotron emissivity maps and light curves using two assumptions: (i) magnetic field compression; (ii) amplification of the magnetic field at the shock.We find that models with magnetic field amplification account better for the observed radio emission.
△ Less
Submitted 8 August, 2021; v1 submitted 5 June, 2021;
originally announced June 2021.
-
Angled pulsar magnetospheres
Authors:
Maxim Lyutikov,
Praveen Sharma,
Maxim Barkov
Abstract:
We consider magnetospheric structure of rotating neutron stars with internally twisted axisymmetric magnetic fields. The twist-induced and rotation-induced toroidal magnetic fields align/counter-align in different hemispheres. Using analytical and numerical calculations (with PHAEDRA code) we show that as a result the North-South symmetry is broken: the magnetosphere and the wind become "angled",…
▽ More
We consider magnetospheric structure of rotating neutron stars with internally twisted axisymmetric magnetic fields. The twist-induced and rotation-induced toroidal magnetic fields align/counter-align in different hemispheres. Using analytical and numerical calculations (with PHAEDRA code) we show that as a result the North-South symmetry is broken: the magnetosphere and the wind become "angled", of conical shape. Angling of the magnetosphere affects the spindown (making it smaller for mild twists), makes the return current split unequally at the Y-point, produces anisotropic wind and linear acceleration that may dominate over gravitational acceleration in the Galactic potential and give a total kick up to $\sim 100$ km/s. We also consider analytically the structure of the Y-point in the twisted magnetosphere, and provide estimate of the internal twist beyond which no stable solutions exist: over-twisted magnetospheres must produce plasma ejection events.
△ Less
Submitted 16 April, 2021; v1 submitted 31 March, 2021;
originally announced March 2021.
-
Peeking Between the Pulses: The Far-UV Spectrum of the Previously Unseen White Dwarf in AR Scorpii
Authors:
Peter Garnavich,
Colin Littlefield,
Maxim Lyutikov,
Maxim Barkov
Abstract:
The compact object in the interacting binary AR Sco has widely been presumed to be a rapidly rotating, magnetized white dwarf (WD), but it has never been detected directly. Isolating its spectrum has proven difficult because the spin-down of the WD generates pulsed synchrotron radiation that far outshines the WD's photosphere. As a result, a previous study of AR Sco was unable to detect the WD in…
▽ More
The compact object in the interacting binary AR Sco has widely been presumed to be a rapidly rotating, magnetized white dwarf (WD), but it has never been detected directly. Isolating its spectrum has proven difficult because the spin-down of the WD generates pulsed synchrotron radiation that far outshines the WD's photosphere. As a result, a previous study of AR Sco was unable to detect the WD in the averaged far-ultraviolet spectrum from a Hubble Space Telescope (HST) observation. In an effort to unveil the WD's spectrum, we reanalyze these HST observations by calculating the average spectrum in the troughs between synchrotron pulses. We identify weak spectral features from the previously unseen WD and estimate its surface temperature to be 11500$\pm$500K. Additionally, during the synchrotron pulses, we detect broad Lyman-$α$ absorption consistent with hot WD spectral models. We infer the presence of a pair of hotspots, with temperatures between 23000K and 28000K, near the magnetic poles of the WD. As the WD is not expected to be accreting from its companion, we describe two possible mechanisms for heating the magnetic poles. The Lyman-$α$ absorption of the hotspots appears relatively undistorted by Zeeman splitting, constraining the WD's field strength to be 100 MG, but the data are insufficient to search for the subtle Zeeman splits expected at lower field strengths.
△ Less
Submitted 17 December, 2020;
originally announced December 2020.
-
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.
-
Sensitivity of the Cherenkov Telescope Array for probing cosmology and fundamental physics with gamma-ray propagation
Authors:
The Cherenkov Telescope Array Consortium,
:,
H. Abdalla,
H. Abe,
F. Acero,
A. Acharyya,
R. Adam,
I. Agudo,
A. Aguirre-Santaella,
R. Alfaro,
J. Alfaro,
C. Alispach,
R. Aloisio,
R. Alves B,
L. Amati,
E. Amato,
G. Ambrosi,
E. O. Angüner,
A. Araudo,
T. Armstrong,
F. Arqueros,
L. Arrabito,
K. Asano,
Y. Ascasíbar,
M. Ashley
, et al. (474 additional authors not shown)
Abstract:
The Cherenkov Telescope Array (CTA), the new-generation ground-based observatory for $γ$-ray astronomy, provides unique capabilities to address significant open questions in astrophysics, cosmology, and fundamental physics. We study some of the salient areas of $γ$-ray cosmology that can be explored as part of the Key Science Projects of CTA, through simulated observations of active galactic nucle…
▽ More
The Cherenkov Telescope Array (CTA), the new-generation ground-based observatory for $γ$-ray astronomy, provides unique capabilities to address significant open questions in astrophysics, cosmology, and fundamental physics. We study some of the salient areas of $γ$-ray cosmology that can be explored as part of the Key Science Projects of CTA, through simulated observations of active galactic nuclei (AGN) and of their relativistic jets. Observations of AGN with CTA will enable a measurement of $γ$-ray absorption on the extragalactic background light with a statistical uncertainty below 15% up to a redshift $z=2$ and to constrain or detect $γ$-ray halos up to intergalactic-magnetic-field strengths of at least 0.3pG. Extragalactic observations with CTA also show promising potential to probe physics beyond the Standard Model. The best limits on Lorentz invariance violation from $γ$-ray astronomy will be improved by a factor of at least two to three. CTA will also probe the parameter space in which axion-like particles could constitute a significant fraction, if not all, of dark matter. We conclude on the synergies between CTA and other upcoming facilities that will foster the growth of $γ$-ray cosmology.
△ Less
Submitted 26 February, 2021; v1 submitted 3 October, 2020;
originally announced October 2020.
-
Sensitivity of the Cherenkov Telescope Array to a dark matter signal from the Galactic centre
Authors:
The Cherenkov Telescope Array Consortium,
:,
A. Acharyya,
R. Adam,
C. Adams,
I. Agudo,
A. Aguirre-Santaella,
R. Alfaro,
J. Alfaro,
C. Alispach,
R. Aloisio,
R. Alves Batista,
L. Amati,
G. Ambrosi,
E. O. Angüner,
L. A. Antonelli,
C. Aramo,
A. Araudo,
T. Armstrong,
F. Arqueros,
K. Asano,
Y. Ascasíbar,
M. Ashley,
C. Balazs,
O. Ballester
, et al. (427 additional authors not shown)
Abstract:
We provide an updated assessment of the power of the Cherenkov Telescope Array (CTA) to search for thermally produced dark matter at the TeV scale, via the associated gamma-ray signal from pair-annihilating dark matter particles in the region around the Galactic centre. We find that CTA will open a new window of discovery potential, significantly extending the range of robustly testable models giv…
▽ More
We provide an updated assessment of the power of the Cherenkov Telescope Array (CTA) to search for thermally produced dark matter at the TeV scale, via the associated gamma-ray signal from pair-annihilating dark matter particles in the region around the Galactic centre. We find that CTA will open a new window of discovery potential, significantly extending the range of robustly testable models given a standard cuspy profile of the dark matter density distribution. Importantly, even for a cored profile, the projected sensitivity of CTA will be sufficient to probe various well-motivated models of thermally produced dark matter at the TeV scale. This is due to CTA's unprecedented sensitivity, angular and energy resolutions, and the planned observational strategy. The survey of the inner Galaxy will cover a much larger region than corresponding previous observational campaigns with imaging atmospheric Cherenkov telescopes. CTA will map with unprecedented precision the large-scale diffuse emission in high-energy gamma rays, constituting a background for dark matter searches for which we adopt state-of-the-art models based on current data. Throughout our analysis, we use up-to-date event reconstruction Monte Carlo tools developed by the CTA consortium, and pay special attention to quantifying the level of instrumental systematic uncertainties, as well as background template systematic errors, required to probe thermally produced dark matter at these energies.
"Full likelihood tables complementing our analysis are provided here [ https://doi.org/10.5281/zenodo.4057987 ]"
△ Less
Submitted 30 January, 2021; v1 submitted 31 July, 2020;
originally announced July 2020.
-
Tilting Instability of Magnetically Confined Spheromaks
Authors:
Riddhi Mehta,
Maxim Barkov,
Lorenzo Sironi,
Maxim Lyutikov
Abstract:
We consider the tilting instability of a magnetically confined spheromak using 3D MHD and relativistic PIC calculations with an application to astrophysical plasmas, specifically those occurring in magnetar magnetospheres. The instability is driven by the counter alignment of the spheromak's intrinsic magnetic dipole with the external magnetic field. Initially the spheromak rotates - tilts - tryin…
▽ More
We consider the tilting instability of a magnetically confined spheromak using 3D MHD and relativistic PIC calculations with an application to astrophysical plasmas, specifically those occurring in magnetar magnetospheres. The instability is driven by the counter alignment of the spheromak's intrinsic magnetic dipole with the external magnetic field. Initially the spheromak rotates - tilts - trying to lower its magnetic potential energy. As a result a current sheet forms between the internal magnetic field of a spheromak and the confining field. Magnetic reconnection sets in; this leads to the annihilation of the newly counter-aligned magnetic flux of the spheromak. This occurs on few Alfvén time scales. In the case of higher order (second order) spheromak, the internal core is first pushed out of the envelope, resulting in formation of two nearly independent tilting spheromaks. Thus, the magnetically twisted outer shell cannot stabilize the inner core. During dissipation, helicity of the initial spheromak is carried away by torsional Alfvén waves, violating the assumptions of the Taylor relaxation theorem. In applications to magnetars' giant flares, fast development of tilting instabilities, and no stabilization of the higher order spheromaks, make it unlikely that trapped spheromaks are responsible for the tail emission lasting hundreds of seconds.
△ Less
Submitted 9 July, 2020; v1 submitted 25 June, 2020;
originally announced June 2020.
-
Dynamics and emission of wind-powered afterglows of gamma-ray bursts: flares, plateaus and steep decays
Authors:
Maxim Barkov,
Yonggang Luo,
Maxim Lyutikov
Abstract:
We develop a model of early X-ray afterglows of gamma-ray bursts originating from the reverse shock (RS) propagating through ultra-relativistic, highly magnetized pulsar-like winds produced by long-lasting central engines. We first perform fluid and MHD numerical simulations of relativistic double explosions. We demonstrate that even for constant properties of the wind a variety of temporal behavi…
▽ More
We develop a model of early X-ray afterglows of gamma-ray bursts originating from the reverse shock (RS) propagating through ultra-relativistic, highly magnetized pulsar-like winds produced by long-lasting central engines. We first perform fluid and MHD numerical simulations of relativistic double explosions. We demonstrate that even for constant properties of the wind a variety of temporal behaviors can be produced, depending on the energy of the initial explosion and the wind power, the delay time for the switch-on of the wind, and magnetization of the wind. X-ray emission of the highly magnetized RS occurs in the fast cooling regime - this ensures high radiative efficiency and allows fast intensity variations. We demonstrate that: (i) RS emission naturally produces light curves showing power-law temporal evolution with various temporal indices; (ii) mild wind power, of the order of $\sim 10^{46}$ erg s$^{-1}$ (equivalent isotropic), can reproduce the afterglows' plateau phase; (iii) termination of the wind can produce sudden steep decays; (iv) short-duration afterglow flares are due to mild variations in the wind luminosity, with small total injected energy.
△ Less
Submitted 16 February, 2021; v1 submitted 28 April, 2020;
originally announced April 2020.
-
Magnetospheric interaction in white dwarf binaries AR Sco and AE Aqr
Authors:
Maxim Lyutikov,
Maxim Barkov,
Matthew Route,
Dinshaw Balsara,
Peter Garnavich,
Colin Littlefield
Abstract:
We develop a model of the white dwarf (WD) - red dwarf (RD) binaries AR Sco and AE Aqr as systems in a transient propeller stage of highly asynchronous intermediate polars. The WDs are relatively weakly magnetized with magnetic field of $\sim 10^6$ G. We explain the salient observed features of the systems due to the magnetospheric interaction of two stars. Currently, the WD's spin-down is determi…
▽ More
We develop a model of the white dwarf (WD) - red dwarf (RD) binaries AR Sco and AE Aqr as systems in a transient propeller stage of highly asynchronous intermediate polars. The WDs are relatively weakly magnetized with magnetic field of $\sim 10^6$ G. We explain the salient observed features of the systems due to the magnetospheric interaction of two stars. Currently, the WD's spin-down is determined by the mass loading of the WD's magnetosphere from the RD's at a mild rate of $\dot{M}_{WD} \sim 10^{-11} M_\odot $/yr. Typical loading distance is determined by the ionization of the RD's wind by the WD's UV flux. The WD was previously spun up by a period of high accretion rate from the RD via Roch lobe overflow with $\dot{M} \sim 10^{-9} M_\odot $/yr, acting for as short a period as tens of thousands of years. The non-thermal X-ray and optical synchrotron emitting particles originate in reconnection events in the magnetosphere of the WD due to the interaction with the flow from the RD. In the case of AR Sco, the reconnection events produce signals at the WD's rotation and beat periods - this modulation is due to the changing relative orientation of the companions' magnetic moments and resulting variable reconnection conditions. Radio emission is produced in the magnetosphere of the RD, we hypothesize, in a way that it is physically similar to the Io-induced Jovian decametric radiation.
△ Less
Submitted 23 April, 2020;
originally announced April 2020.
-
Fast moving pulsars as probes of interstellar medium
Authors:
Maxim Barkov,
Maxim Lyutikov,
Dmitry Khangulyan
Abstract:
Pulsars moving through ISM produce bow shocks detected in hydrogen H$α$ line emission. The morphology of the bow shock nebulae allows one to probe the properties of ISM on scales $\sim 0.01$ pc and smaller. We performed 2D RMHD modeling of the pulsar bow shock and simulated the corresponding H$α$ emission morphology. We find that even a mild spatial inhomogeneity of ISM density, $δρ/ρ\sim 1$, lead…
▽ More
Pulsars moving through ISM produce bow shocks detected in hydrogen H$α$ line emission. The morphology of the bow shock nebulae allows one to probe the properties of ISM on scales $\sim 0.01$ pc and smaller. We performed 2D RMHD modeling of the pulsar bow shock and simulated the corresponding H$α$ emission morphology. We find that even a mild spatial inhomogeneity of ISM density, $δρ/ρ\sim 1$, leads to significant variations of the shape of the shock seen in H$α$ line emission. We successfully reproduce the morphology of the Guitar Nebula. We infer quasi-periodic density variations in the warm component of ISM with a characteristic length of $\sim0.1$~pc. Structures of this scale might be also responsible for the formation of the fine features seen at the forward shock of Tycho SNR in X-rays. Formation of such short periodic density structures in the warm component of ISM is puzzling, and bow-shock nebulae provide unique probes to study this phenomenon.
△ Less
Submitted 3 June, 2020; v1 submitted 27 February, 2020;
originally announced February 2020.
-
FRB-periodicity: mild pulsars in tight O/B-star binaries
Authors:
Maxim Lyutikov,
Maxim Barkov,
Dimitrios Giannios
Abstract:
Periodicities observed in two Fast Radio Burst (FRB) sources (16 days in FRB 180916.J0158+65 and 160 days in FRB 121102) are consistent with that of tight, stellar mass binary systems. In the case of FRB 180916.J0158+65 the primary is an early OB-type star with mass loss rate $\dot{M} \sim 10^{-8}- 10^{-7} M_\odot$ yr$^{-1}$, and the secondary a neutron star. The observed periodicity is not intrin…
▽ More
Periodicities observed in two Fast Radio Burst (FRB) sources (16 days in FRB 180916.J0158+65 and 160 days in FRB 121102) are consistent with that of tight, stellar mass binary systems. In the case of FRB 180916.J0158+65 the primary is an early OB-type star with mass loss rate $\dot{M} \sim 10^{-8}- 10^{-7} M_\odot$ yr$^{-1}$, and the secondary a neutron star. The observed periodicity is not intrinsic to the FRB's source, but is due to the orbital phase-dependent modulation of the absorption conditions in the massive star's wind. The observed relatively narrow FRB activity window implies that the primary's wind dynamically dominates that of the pulsar, $η= L_{sd}/(\dot{M} v_w c) \leq 1$, where $L_{sd} $ is pulsar spin-down, $\dot{M}$ is the primary's wind mass loss rate and $v_w$ is its velocity. The condition $η\leq 1$ requires mildly powerful pulsar with $L_{sd} \lesssim 10^{37}$ erg $s^{-1}$. The observations are consistent with magnetically-powered radio emission originating in the magnetospheres of strongly magnetized neutron stars, the classical magnetars.
△ Less
Submitted 7 August, 2022; v1 submitted 5 February, 2020;
originally announced February 2020.
-
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.
-
Monte Carlo studies for the optimisation of the Cherenkov Telescope Array layout
Authors:
A. Acharyya,
I. Agudo,
E. O. Angüner,
R. Alfaro,
J. Alfaro,
C. Alispach,
R. Aloisio,
R. Alves Batista,
J. -P. Amans,
L. Amati,
E. Amato,
G. Ambrosi,
L. A. Antonelli,
C. Aramo,
T. Armstrong,
F. Arqueros,
L. Arrabito,
K. Asano,
H. Ashkar,
C. Balazs,
M. Balbo,
B. Balmaverde,
P. Barai,
A. Barbano,
M. Barkov
, et al. (445 additional authors not shown)
Abstract:
The Cherenkov Telescope Array (CTA) is the major next-generation observatory for ground-based very-high-energy gamma-ray astronomy. It will improve the sensitivity of current ground-based instruments by a factor of five to twenty, depending on the energy, greatly improving both their angular and energy resolutions over four decades in energy (from 20 GeV to 300 TeV). This achievement will be possi…
▽ More
The Cherenkov Telescope Array (CTA) is the major next-generation observatory for ground-based very-high-energy gamma-ray astronomy. It will improve the sensitivity of current ground-based instruments by a factor of five to twenty, depending on the energy, greatly improving both their angular and energy resolutions over four decades in energy (from 20 GeV to 300 TeV). This achievement will be possible by using tens of imaging Cherenkov telescopes of three successive sizes. They will be arranged into two arrays, one per hemisphere, located on the La Palma island (Spain) and in Paranal (Chile). We present here the optimised and final telescope arrays for both CTA sites, as well as their foreseen performance, resulting from the analysis of three different large-scale Monte Carlo productions.
△ Less
Submitted 2 April, 2019;
originally announced April 2019.
-
Nuclear $γ$-ray emission from very hot accretion flows
Authors:
Ervin Kafexhiu,
Felix Aharonian,
Maxim Barkov
Abstract:
Optically thin accretion plasmas can reach ion temperatures $T_{\rm i} \geq 10^{10}$K and thus trigger nuclear reactions. Using a large nuclear interactions network, we studied the radial evolution of the chemical composition of the accretion flow toward the black hole and computed the emissivity in nuclear $γ$-ray lines. In the advection dominated accretion flow (ADAF) regime, CNO and heavier nuc…
▽ More
Optically thin accretion plasmas can reach ion temperatures $T_{\rm i} \geq 10^{10}$K and thus trigger nuclear reactions. Using a large nuclear interactions network, we studied the radial evolution of the chemical composition of the accretion flow toward the black hole and computed the emissivity in nuclear $γ$-ray lines. In the advection dominated accretion flow (ADAF) regime, CNO and heavier nuclei are destroyed before reaching the last stable orbit. The overall luminosity in the de-excitation lines for a solar composition of plasma can be as high as few times $10^{-5}$ the accretion luminosity ($\dot{M}c^2$) and can be increased for heavier compositions up to $10^{-3}$. The efficiency of transformation of the kinetic energy of the outflow into high energy ($\geq 100$~MeV) $γ$-rays through the production and decay of $π^0$-mesons can be higher, up to $10^{-2}$ of the accretion luminosity. We show that in the ADAF model up to 15 percent of the mass of accretion matter can `evaporate' in the form of neutrons.
△ Less
Submitted 27 March, 2019; v1 submitted 25 July, 2018;
originally announced July 2018.
-
Gamma-ray emission of hot astrophysical plasmas
Authors:
Ervin Kafexhiu,
Felix Aharonian,
Maxim Barkov
Abstract:
Very hot plasmas with ion temperature exceeding $10^{10}$ K can be formed in certain astrophysical environments. The distinct radiation signature of such plasmas is the $γ$-ray emission dominated by the prompt de-excitation nuclear lines and $π^0$-decay $γ$-rays. Using a large nuclear reaction network, we compute the time evolution of the chemical composition of such hot plasmas and their $γ$-ray…
▽ More
Very hot plasmas with ion temperature exceeding $10^{10}$ K can be formed in certain astrophysical environments. The distinct radiation signature of such plasmas is the $γ$-ray emission dominated by the prompt de-excitation nuclear lines and $π^0$-decay $γ$-rays. Using a large nuclear reaction network, we compute the time evolution of the chemical composition of such hot plasmas and their $γ$-ray line emissivity. At higher energies, we provide simple but accurate analytical presentations for the $π^0$-meson production rate and the corresponding $π^0\to2γ$ emissivity derived for the Maxwellian distribution of protons. We discuss the impact of the possible deviation of the high energy tail of the particle distribution function from the "nominal" Maxwellian distribution on the plasma $γ$-ray emissivity.
△ Less
Submitted 27 March, 2019; v1 submitted 16 July, 2018;
originally announced July 2018.
-
A hydrodynamics-informed, radiation model for HESS J0632$+$057 from radio to gamma rays
Authors:
Maxim Barkov,
Valenti Bosch-Ramon
Abstract:
Relativistic hydrodynamical simulations of the eccentric gamma-ray binary HESS J0632$+$057, show that the energy of a putative pulsar wind should accumulate in the binary surroundings between periastron and apastron, being released by fast advection close to apastron. To assess whether this could lead to a maximum of the non-thermal emission before apastron, we derive simple prescriptions for the…
▽ More
Relativistic hydrodynamical simulations of the eccentric gamma-ray binary HESS J0632$+$057, show that the energy of a putative pulsar wind should accumulate in the binary surroundings between periastron and apastron, being released by fast advection close to apastron. To assess whether this could lead to a maximum of the non-thermal emission before apastron, we derive simple prescriptions for the non-thermal energy content, the radiation efficiency, and the impact of energy losses on non-thermal particles, in the simulated hydrodynamical flow. These prescriptions are used to estimate the non-thermal emission in radio, X-rays, GeV, and TeV, from the shocked pulsar wind in a binary system simulated using a simplified 3-dimensional scheme for several orbital cycles. Lightcurves at different wavelengths are derived, together with synthetic radio images for different orbital phases. The dominant peak in the computed lightcurves is broad and appears close to, but before, apastron. This peak is followed by a quasi-plateau shape, and a minor peak only in gamma rays right after periastron. The radio maps show ejection of radio blobs before apastron in the periastron-apastron direction. The results show that a scenario with a highly eccentric high-mass binary hosting a young pulsar can explain the general phenomenology of HESS J0632$+$057: despite its simplicity, the adopted approach yields predictions that are robust at a semi-quantitative level and consistent with multiwavelength observations.
△ Less
Submitted 14 June, 2018;
originally announced June 2018.
-
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.
-
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.
-
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.
-
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.
-
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.
-
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.
-
Science with the Cherenkov Telescope Array
Authors:
The Cherenkov Telescope Array Consortium,
:,
B. S. Acharya,
I. Agudo,
I. Al Samarai,
R. Alfaro,
J. Alfaro,
C. Alispach,
R. Alves Batista,
J. -P. Amans,
E. Amato,
G. Ambrosi,
E. Antolini,
L. A. Antonelli,
C. Aramo,
M. Araya,
T. Armstrong,
F. Arqueros,
L. Arrabito,
K. Asano,
M. Ashley,
M. Backes,
C. Balazs,
M. Balbo,
O. Ballester
, et al. (558 additional authors not shown)
Abstract:
The Cherenkov Telescope Array, CTA, will be the major global observatory for very high energy gamma-ray astronomy over the next decade and beyond. The scientific potential of CTA is extremely broad: from understanding the role of relativistic cosmic particles to the search for dark matter. CTA is an explorer of the extreme universe, probing environments from the immediate neighbourhood of black ho…
▽ More
The Cherenkov Telescope Array, CTA, will be the major global observatory for very high energy gamma-ray astronomy over the next decade and beyond. The scientific potential of CTA is extremely broad: from understanding the role of relativistic cosmic particles to the search for dark matter. CTA is an explorer of the extreme universe, probing environments from the immediate neighbourhood of black holes to cosmic voids on the largest scales. Covering a huge range in photon energy from 20 GeV to 300 TeV, CTA will improve on all aspects of performance with respect to current instruments.
The observatory will operate arrays on sites in both hemispheres to provide full sky coverage and will hence maximize the potential for the rarest phenomena such as very nearby supernovae, gamma-ray bursts or gravitational wave transients. With 99 telescopes on the southern site and 19 telescopes on the northern site, flexible operation will be possible, with sub-arrays available for specific tasks. CTA will have important synergies with many of the new generation of major astronomical and astroparticle observatories. Multi-wavelength and multi-messenger approaches combining CTA data with those from other instruments will lead to a deeper understanding of the broad-band non-thermal properties of target sources.
The CTA Observatory will be operated as an open, proposal-driven observatory, with all data available on a public archive after a pre-defined proprietary period. Scientists from institutions worldwide have combined together to form the CTA Consortium. This Consortium has prepared a proposal for a Core Programme of highly motivated observations. The programme, encompassing approximately 40% of the available observing time over the first ten years of CTA operation, is made up of individual Key Science Projects (KSPs), which are presented in this document.
△ Less
Submitted 21 January, 2018; v1 submitted 22 September, 2017;
originally announced September 2017.
-
Cherenkov Telescope Array Contributions to the 35th International Cosmic Ray Conference (ICRC2017)
Authors:
F. Acero,
B. S. Acharya,
V. Acín Portella,
C. Adams,
I. Agudo,
F. Aharonian,
I. Al Samarai,
A. Alberdi,
M. Alcubierre,
R. Alfaro,
J. Alfaro,
C. Alispach,
R. Aloisio,
R. Alves Batista,
J. -P. Amans,
E. Amato,
L. Ambrogi,
G. Ambrosi,
M. Ambrosio,
J. Anderson,
M. Anduze,
E. O. Angüner,
E. Antolini,
L. A. Antonelli,
V. Antonuccio
, et al. (1117 additional authors not shown)
Abstract:
List of contributions from the Cherenkov Telescope Array Consortium presented at the 35th International Cosmic Ray Conference, July 12-20 2017, Busan, Korea.
List of contributions from the Cherenkov Telescope Array Consortium presented at the 35th International Cosmic Ray Conference, July 12-20 2017, Busan, Korea.
△ Less
Submitted 24 October, 2017; v1 submitted 11 September, 2017;
originally announced September 2017.
-
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.
-
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.
-
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.
-
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.
-
Prospects for CTA observations of the young SNR RX J1713.7-3946
Authors:
The CTA Consortium,
:,
F. Acero,
R. Aloisio,
J. Amans,
E. Amato,
L. A. Antonelli,
C. Aramo,
T. Armstrong,
F. Arqueros,
K. Asano,
M. Ashley,
M. Backes,
C. Balazs,
A. Balzer,
A. Bamba,
M. Barkov,
J. A. Barrio,
W. Benbow,
K. Bernlöhr,
V. Beshley,
C. Bigongiari,
A. Biland,
A. Bilinsky,
E. Bissaldi
, et al. (359 additional authors not shown)
Abstract:
We perform simulations for future Cherenkov Telescope Array (CTA) observations of RX~J1713.7$-$3946, a young supernova remnant (SNR) and one of the brightest sources ever discovered in very-high-energy (VHE) gamma rays. Special attention is paid to explore possible spatial (anti-)correlations of gamma rays with emission at other wavelengths, in particular X-rays and CO/H{\sc i} emission. We presen…
▽ More
We perform simulations for future Cherenkov Telescope Array (CTA) observations of RX~J1713.7$-$3946, a young supernova remnant (SNR) and one of the brightest sources ever discovered in very-high-energy (VHE) gamma rays. Special attention is paid to explore possible spatial (anti-)correlations of gamma rays with emission at other wavelengths, in particular X-rays and CO/H{\sc i} emission. We present a series of simulated images of RX J1713.7$-$3946 for CTA based on a set of observationally motivated models for the gamma-ray emission. In these models, VHE gamma rays produced by high-energy electrons are assumed to trace the non-thermal X-ray emission observed by {\it XMM-Newton}, whereas those originating from relativistic protons delineate the local gas distributions. The local atomic and molecular gas distributions are deduced by the NANTEN team from CO and H{\sc i} observations. Our primary goal is to show how one can distinguish the emission mechanism(s) of the gamma rays (i.e., hadronic vs leptonic, or a mixture of the two) through information provided by their spatial distribution, spectra, and time variation. This work is the first attempt to quantitatively evaluate the capabilities of CTA to achieve various proposed scientific goals by observing this important cosmic particle accelerator.
△ Less
Submitted 13 April, 2017;
originally announced April 2017.
-
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.
-
Contributions of the Cherenkov Telescope Array (CTA) to the 6th International Symposium on High-Energy Gamma-Ray Astronomy (Gamma 2016)
Authors:
The CTA Consortium,
:,
A. Abchiche,
U. Abeysekara,
Ó. Abril,
F. Acero,
B. S. Acharya,
C. Adams,
G. Agnetta,
F. Aharonian,
A. Akhperjanian,
A. Albert,
M. Alcubierre,
J. Alfaro,
R. Alfaro,
A. J. Allafort,
R. Aloisio,
J. -P. Amans,
E. Amato,
L. Ambrogi,
G. Ambrosi,
M. Ambrosio,
J. Anderson,
M. Anduze,
E. O. Angüner
, et al. (1387 additional authors not shown)
Abstract:
List of contributions from the Cherenkov Telescope Array (CTA) Consortium presented at the 6th International Symposium on High-Energy Gamma-Ray Astronomy (Gamma 2016), July 11-15, 2016, in Heidelberg, Germany.
List of contributions from the Cherenkov Telescope Array (CTA) Consortium presented at the 6th International Symposium on High-Energy Gamma-Ray Astronomy (Gamma 2016), July 11-15, 2016, in Heidelberg, Germany.
△ Less
Submitted 17 October, 2016;
originally announced October 2016.
-
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.
-
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.
-
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.
-
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.
-
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.
-
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.
-
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.
-
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.
-
A multi-dimensional numerical scheme for two-fluid Relativistic MHD
Authors:
Maxim Barkov,
Serguei S. Komissarov,
Vitaly Korolev,
Andrey Zankovich
Abstract:
The paper describes an explicit multi-dimensional numerical scheme for Special Relativistic Two-Fluid Magnetohydrodynamics of electron-positron plasma and a suit of test problems. The scheme utilizes Cartesian grid and the third order WENO interpolation. The time integration is carried out using the third order TVD method of Runge-Kutta type, thus ensuring overall third order accuracy on smooth so…
▽ More
The paper describes an explicit multi-dimensional numerical scheme for Special Relativistic Two-Fluid Magnetohydrodynamics of electron-positron plasma and a suit of test problems. The scheme utilizes Cartesian grid and the third order WENO interpolation. The time integration is carried out using the third order TVD method of Runge-Kutta type, thus ensuring overall third order accuracy on smooth solutions. The magnetic field is kept near divergence-free by means of the method of generalized Lagrange multiplier. The test simulations, which include linear and non-linear continuous plasma waves, shock waves, strong explosions and the tearing instability, show that the scheme is sufficiently robust and confirm its accuracy.
△ Less
Submitted 16 April, 2014; v1 submitted 20 September, 2013;
originally announced September 2013.