Showing 1–50 of 86 results for author: Toma, K

Multi-wavelength Emission from Jets and Magnetically Arrested Disks in Nearby Radio Galaxies: Application to M87

Authors: Riku Kuze, Shigeo S. Kimura, Kenji Toma

Abstract: Radio galaxies are a subclass of active galactic nuclei that drive relativistic jets from their center and are observed in radio to very-high-energy gamma rays. The emission mechanisms and regions are still unknown. High-energy gamma rays can be explained by the emission from the magnetically arrested disks (MADs) around the central supermassive black hole, for which the magnetic flux threading th… ▽ More Radio galaxies are a subclass of active galactic nuclei that drive relativistic jets from their center and are observed in radio to very-high-energy gamma rays. The emission mechanisms and regions are still unknown. High-energy gamma rays can be explained by the emission from the magnetically arrested disks (MADs) around the central supermassive black hole, for which the magnetic flux threading the black hole is in a saturation level, although the emission from the MADs does not explain the optical and X-ray data. We construct a multi-wavelength emission model in which the optical and X-ray emission is emitted by jets and the gamma rays by MADs. Our model takes into account the particle injection by the magnetic reconnection at the jet base close to the black hole and particle entrainment from the ambient gas at the jet emission zone. We apply our model to M87 and find that our model can explain the simultaneous multi-wavelength data. In our model, the emission from the jets is the synchrotron radiation of the nonthermal electrons accelerated by magnetic reconnection, and the emission from the MADs is the synchrotron radiation mainly of the nonthermal protons accelerated by turbulence. We also find that the strong plasma entrainment is necessary to explain the multi-wavelength data. Our model will be tested by variability analysis among the multi-wavelength data. △ Less

Submitted 4 June, 2024; originally announced June 2024.

Comments: 17 pages, 8 figures, 1 table

Varying linear polarisation in the dust-free GRB 210610B

Authors: J. F. Agüí Fernández, A. de Ugarte Postigo, C. C. Thöne, S. Kobayashi, A. Rossi, K. Toma, M. Jelínek, D. A. Kann, S. Covino, K. Wiersema, D. Hartmann, P. Jakobsson, A. Martin-Carrillo, A. Melandri, M. De Pasquale, G. Pugliese, S. Savaglio, R. L. C. Starling, J. Štrobl, M. Della Valle, S. de Wet, T. Zafar

Abstract: Long gamma ray bursts (GRBs) are produced by the collapse of some very massive stars, which emit ultra-relativistic jets. When the jets collide with the interstellar medium they decelerate and generate the so-called afterglow emission, which has been observed to be polarised. In this work we study the polarimetric evolution of GRB 210610B afterglow, at $z = 1.1341$. This allows to evaluate the rol… ▽ More Long gamma ray bursts (GRBs) are produced by the collapse of some very massive stars, which emit ultra-relativistic jets. When the jets collide with the interstellar medium they decelerate and generate the so-called afterglow emission, which has been observed to be polarised. In this work we study the polarimetric evolution of GRB 210610B afterglow, at $z = 1.1341$. This allows to evaluate the role of geometric and/or magnetic mechanisms in the GRB afterglow polarisation. We observed GRB 210610B using imaging polarimetry with CAFOS on the 2.2 m Calar Alto Telescope and FORS2 on the 4 $\times$ 8.1 m Very Large Telescope. Complementary optical spectroscopy was obtained with OSIRIS on the 10.4 m Gran Telescopio Canarias. We study the GRB light-curve from X-rays to optical bands and the Spectral Energy Distribution (SED). This allows us to strongly constrain the line-of-sight extinction. Finally, we study the GRB host galaxy using optical/NIR data to fit the SED and derive its integrated properties. GRB 210610B had a bright afterglow with a negligible line-of-sight extinction. Polarimetry was obtained at three epochs: during an early plateau phase, at the time when the light curve breaks, and after the light curve steepened. We observe an initial polarisation of $\sim 4\%$ that goes to zero at the time of the break, and then increases again to $\sim 2\%$ with a change of the position angle of $54 \pm 9$ deg. The spectrum show features with very low equivalent widths, indicating a small amount of material in the line-of-sight within the host. The lack of dust and the low amount of material on the line-of-sight to GRB 210610B allow us to study the intrinsic polarisation of the GRB optical afterglow. We find the GRB polarisation signals are consistent with ordered magnetic fields in refreshed shock or/and hydrodynamics-scale turbulent fields in the forward shock. △ Less

Submitted 4 March, 2024; originally announced March 2024.

Comments: 17 pages, 7 figures, 5 tables. Accepted for publication in Astronomy and Astrophysics 23 Feb 2024

1D GRPIC Simulations of Stellar-Mass Black Hole Magnetospheres: Semi-Analytic Model of Gamma-Rays from Gaps

Authors: Koki Kin, Shota Kisaka, Kenji Toma, Shigeo S. Kimura, Amir Levinson

Abstract: In the absence of a sufficient amount of plasma injection into the black hole (BH) magnetosphere, the force-free state of the magnetosphere cannot be maintained, leading to the emergence of strong, time-dependent, longitudinal electric field (spark gap). Recent studies of supermassive BH magnetospheres by using analytical methods and particle-in-cell (PIC) simulations propose the possibility of th… ▽ More In the absence of a sufficient amount of plasma injection into the black hole (BH) magnetosphere, the force-free state of the magnetosphere cannot be maintained, leading to the emergence of strong, time-dependent, longitudinal electric field (spark gap). Recent studies of supermassive BH magnetospheres by using analytical methods and particle-in-cell (PIC) simulations propose the possibility of the efficient particle acceleration and consequent gamma-ray emissions in the spark gap. In this work, we perform one-dimensional general relativistic PIC simulations to examine the gamma-ray emission from stellar-mass BH magnetospheres. We find that intermittent spark gaps emerge and particles are efficiently accelerated, in a similar manner to the supermassive BH case. We build a semi-analytic model of the plasma dynamics and radiative processes which reproduces the maximum electron energies and peak gamma-ray luminosities in the simulation results. Based on this model, we show that gamma-ray signals from stellar-mass BHs wandering through the interstellar medium could be detected by gamma-ray telescopes such as the Fermi Large Area Telescope, or the Cherenkov Telescope Array. △ Less

Submitted 19 January, 2024; v1 submitted 19 October, 2023; originally announced October 2023.

Comments: 17 pages, 9 figures, 1 table, accepted in ApJ

Gamma rays from a reverse shock with turbulent magnetic fields in GRB 180720B

Authors: Makoto Arimoto, Katsuaki Asano, Koji S. Kawabata, Kenji Toma, Ramandeep Gill, Jonathan Granot, Masanori Ohno, Shuta Takahashi, Naoki Ogino, Hatsune Goto, Kengo Nakamura, Tatsuya Nakaoka, Kengo Takagi, Miho Kawabata, Masayuki Yamanaka, Mahito Sasada, Soebur Razzaque

Abstract: Gamma-ray bursts (GRBs) are the most electromagnetically luminous cosmic explosions. They are powered by collimated streams of plasma (jets) ejected by a newborn stellar-mass black hole or neutron star at relativistic velocities (near the speed of light). Their short-lived (typically tens of seconds) prompt $γ$-ray emission from within the ejecta is followed by long-lived multi-wavelength afterglo… ▽ More Gamma-ray bursts (GRBs) are the most electromagnetically luminous cosmic explosions. They are powered by collimated streams of plasma (jets) ejected by a newborn stellar-mass black hole or neutron star at relativistic velocities (near the speed of light). Their short-lived (typically tens of seconds) prompt $γ$-ray emission from within the ejecta is followed by long-lived multi-wavelength afterglow emission from the ultra-relativistic forward shock. This shock is driven into the circumburst medium by the GRB ejecta that are in turn decelerated by a mildly-relativistic reverse shock. Forward shock emission was recently detected up to teraelectronvolt-energy $γ$-rays, and such very-high-energy emission was also predicted from the reverse shock. Here we report the detection of optical and gigaelectronvolt-energy $γ$-ray emission from GRB 180720B during the first few hundred seconds, which is explained by synchrotron and inverse-Compton emission from the reverse shock propagating into the ejecta, implying a low-magnetization ejecta. Our optical measurements show a clear transition from the reverse shock to the forward shock driven into the circumburst medium, accompanied by a 90-degree change in the mean polarization angle and fluctuations in the polarization degree and angle. This indicates turbulence with large-scale toroidal and radially-stretched magnetic field structures in the reverse and forward shocks, respectively, which tightly couple to the physics of relativistic shocks and GRB jets -- launching, composition, dissipation and particle acceleration. △ Less

Submitted 6 October, 2023; originally announced October 2023.

Comments: 5 pages, 4 figures (main) plus Methods and Supplementary Methods, accepted for publication

High-energy neutrino emission associated with gravitational-wave signals: effects of cocoon photons and constraints on late-time emission

Authors: Riki Matsui, Shigeo S. Kimura, Kenji Toma, Kohta Murase

Abstract: We investigate prospects for the detection of high-energy neutrinos produced in the prolonged jets of short gamma-ray bursts (sGRBs). The X-ray lightcurves of sGRBs show extended emission components lasting for 100-1000 seconds, which are considered to be produced by prolonged engine activity. Jets by prolonged engine activity should interact with photons in the cocoon formed by the jet propagatio… ▽ More We investigate prospects for the detection of high-energy neutrinos produced in the prolonged jets of short gamma-ray bursts (sGRBs). The X-ray lightcurves of sGRBs show extended emission components lasting for 100-1000 seconds, which are considered to be produced by prolonged engine activity. Jets by prolonged engine activity should interact with photons in the cocoon formed by the jet propagation inside the ejecta of neutron star mergers. We calculate neutrino emission from jets by prolonged engine activity, taking account of the interaction between photons provided from the cocoon and cosmic rays accelerated in the jets. We find that IceCube-Gen2, a future neutrino telescope, with the second-generation gravitational wave detectors will probably be able to observe neutrino signals associated with gravitational waves with around 10 years of operation, regardless of the assumed value of the Lorentz factor of the jets. Neutrino observations may enable us to constrain the dissipation region of the jets. We apply this model to GRB 211211A, a peculiar long GRB whose origin may be a binary neutron-star merger. Our model predicts that IceCube is unlikely to detect any associated neutrino, but a few similar events will be able to put a meaningful constraint on the physical quantities of the prolonged engine activities. △ Less

Submitted 6 July, 2023; v1 submitted 8 February, 2023; originally announced February 2023.

Journal ref: Astrophys.J. 950 (2023) 190

Simultaneous Radio and Optical Polarimetry of GRB 191221B Afterglow

Authors: Yuji Urata, Kenji Toma, Stefano Covino, Klaas Wiersema, Kuiyun Huang, Jiro Shimoda, Asuka Kuwata, Sota Nagao, Keiichi Asada, Hiroshi Nagai, Satoko Takahashi, Chao-En Chung, Glen Petitpas, Kazutaka Yamaoka, Luca Izzo, Johan Fynbo, Antonio de Ugarte Postigo, Maryam Arabsalmani, Makoto Tashiro

Abstract: Gamma-ray bursts (GRBs) are the most luminous transients in the universe and are utilized as probes of early stars, gravitational wave counterparts, and collisionless shock physics. In spite of studies on polarimetry of GRBs in individual wavelengths that characterized intriguing properties of prompt emission and afterglow, no coordinated multi-wavelength measurements have yet been performed. Here… ▽ More Gamma-ray bursts (GRBs) are the most luminous transients in the universe and are utilized as probes of early stars, gravitational wave counterparts, and collisionless shock physics. In spite of studies on polarimetry of GRBs in individual wavelengths that characterized intriguing properties of prompt emission and afterglow, no coordinated multi-wavelength measurements have yet been performed. Here, we report the first coordinated simultaneous polarimetry in the optical and radio bands for the afterglow associated with the typical long GRB 191221B. Our observations successfully caught the radio emission, which is not affected by synchrotron self-absorption, and show that the emission is depolarized in the radio band compared to the optical one. Our simultaneous polarization angle measurement and temporal polarization monitoring indicate the existence of cool electrons that increase the estimate of jet kinetic energy by a factor of $>$ 4 for this GRB afterglow. Further coordinated multi-wavelength polarimetric campaigns would improve our understanding of the total jet energies and magnetic field configurations in the emission regions of various types of GRBs, which are required to comprehend the mass scales of their progenitor systems and the physics of collisionless shocks. △ Less

Submitted 9 December, 2022; originally announced December 2022.

Comments: Author's version. 41 pages, 9 figures, 4 tables

Synchrotron Polarization of Gamma-Ray Burst Afterglow Shocks with Hydrodynamic-scale Turbulent Magnetic Field

Authors: Asuka Kuwata, Kenji Toma, Shigeo S. Kimura, Sara Tomita, Jiro Shimoda

Abstract: Afterglows of gamma-ray bursts (GRBs) are emitted from expanding forward shocks, which are expected to have magnetic field much stronger than the interstellar field, although the origin of the field is a long-standing problem. Two field amplification mechanisms, plasma kinetic instabilities and magnetohydrodynamic instabilities, have been discussed so far. The coherence length scales of the fields… ▽ More Afterglows of gamma-ray bursts (GRBs) are emitted from expanding forward shocks, which are expected to have magnetic field much stronger than the interstellar field, although the origin of the field is a long-standing problem. Two field amplification mechanisms, plasma kinetic instabilities and magnetohydrodynamic instabilities, have been discussed so far. The coherence length scales of the fields amplified by these two processes are different by $7-10$ orders of magnitudes, and the polarimetric observations may distinguish them. We construct a semi-analytic model of the forward shock afterglow polarization under the assumption of hydrodynamic-scale turbulent magnetic field. We perform numerical calculations of synchrotron polarization for the isotropic turbulence and the zero viewing angle. We find that the polarization degrees are $ \sim1~\%$ when the field coherence length scale in the fluid comoving frame is of the order of the thickness of the shocked regions. This range of polarization degree is comparable to that of the observed late-phase optical afterglows. Our model also shows that the radio polarization degrees are comparable to the optical ones on average but can be higher than the optical ones at some time intervals. The polarization angles are shown to vary randomly and continuously. These polarimetric properties are clearly different from the case of plasma kinetic instability. Simultaneous polarimetric observations of GRB afterglows at the radio and optical bands have recently started, which will help us constrain the magnetic field amplification mechanism. △ Less

Submitted 16 December, 2022; v1 submitted 19 August, 2022; originally announced August 2022.

The Photon Ring in M87*

Authors: Avery E. Broderick, Dominic W. Pesce, Paul Tiede, Hung-Yi Pu, Roman Gold, Richard Anantua, Silke Britzen, Chiara Ceccobello, Koushik Chatterjee, Yongjun Chen, Nicholas S. Conroy, Geoffrey B. Crew, Alejandro Cruz-Osorio, Yuzhu Cui, Sheperd S. Doeleman, Razieh Emami, Joseph Farah, Christian M. Fromm, Peter Galison, Boris Georgiev, Luis C. Ho, David J. James, Britton Jeter, Alejandra Jimenez-Rosales, Jun Yi Koay , et al. (26 additional authors not shown)

Abstract: We report measurements of the gravitationally lensed secondary image -- the first in an infinite series of so-called "photon rings" -- around the supermassive black hole M87* via simultaneous modeling and imaging of the 2017 Event Horizon Telescope (EHT) observations. The inferred ring size remains constant across the seven days of the 2017 EHT observing campaign and is consistent with theoretical… ▽ More We report measurements of the gravitationally lensed secondary image -- the first in an infinite series of so-called "photon rings" -- around the supermassive black hole M87* via simultaneous modeling and imaging of the 2017 Event Horizon Telescope (EHT) observations. The inferred ring size remains constant across the seven days of the 2017 EHT observing campaign and is consistent with theoretical expectations, providing clear evidence that such measurements probe spacetime and a striking confirmation of the models underlying the first set of EHT results. The residual diffuse emission evolves on timescales comparable to one week. We are able to detect with high significance a southwestern extension consistent with that expected from the base of a jet that is rapidly rotating in the clockwise direction. This result adds further support to the identification of the jet in M87* with a black hole spin-driven outflow, launched via the Blandford-Znajek process. We present three revised estimates for the mass of M87* based on identifying the modeled thin ring component with the bright ringlike features seen in simulated images, one of which is only weakly sensitive to the astrophysics of the emission region. All three estimates agree with each other and previously reported values. Our strongest mass constraint combines information from both the ring and the diffuse emission region, which together imply a mass-to-distance ratio of $4.20^{+0.12}_{-0.06}~μ{\rm as}$ and a corresponding black hole mass of $(7.13\pm0.39)\times10^9M_\odot$, where the error on the latter is now dominated by the systematic uncertainty arising from the uncertain distance to M87*. △ Less

Submitted 18 August, 2022; originally announced August 2022.

Journal ref: ApJ, 935, 61 (2022)

Magnetic Reconnection in Black-Hole Magnetospheres: Lepton Loading into Jets, Superluminal Radio Blobs, and Multi-wavelength Flares

Authors: Shigeo S. Kimura, Kenji Toma, Hirofumi Noda, Kazuhiro Hada

Abstract: Supermassive black holes in active galactic nuclei launch relativistic jets, as indicated by observed superluminal radio blobs. The energy source of these jets is widely discussed in the theoretical framework of Blandford-Znajek process, the electromagnetic energy extraction from rotating black holes (BHs), while formation mechanism of the radio blobs in the electromagnetically-dominated jets has… ▽ More Supermassive black holes in active galactic nuclei launch relativistic jets, as indicated by observed superluminal radio blobs. The energy source of these jets is widely discussed in the theoretical framework of Blandford-Znajek process, the electromagnetic energy extraction from rotating black holes (BHs), while formation mechanism of the radio blobs in the electromagnetically-dominated jets has been a long-standing problem. Recent high-resolution magnetohydrodynamic simulations of magnetically arrested disks exhibited magnetic reconnection in a transient magnetically-dominated part of the equatorial disk near the BH horizon, which led to a promising scenario of efficient MeV gamma-ray production and subsequent electron-positron pair loading into BH magnetosphere. We develop this scenario to build a theoretical framework on energetics, timescales and particle number density of the superluminal radio blobs and discuss observable signatures in other wavebands. We analytically show that the non-thermal electrons emit broadband photons from optical to multi-MeV bands. The electron-positron pairs produced in the magnetosphere are optically thick for synchrotron-self absorption, so that the injected energy is stored in the plasma. The stored energy is enough to power the superluminal radio blobs observed in M87. This scenario predicts rather dim radio blobs around Sgr A*, which are consistent with no clear detection by current facilities. In addition, this scenario inevitably produces strong X-ray flares in a short timescale, which will be detectable by future X-ray satellites. △ Less

Submitted 31 August, 2022; v1 submitted 3 August, 2022; originally announced August 2022.

Comments: 14 pages, 5 figures, 1 table, accepted for publication in ApJL

High-energy Gamma-rays from Magnetically Arrested Disks in Nearby Radio Galaxies

Authors: Riku Kuze, Shigeo S. Kimura, Kenji Toma

Abstract: The origins of the GeV gamma-rays from nearby radio galaxies are unknown. Hadronic emission from magnetically arrested disks (MADs) around central black holes (BHs) is proposed as a possible scenario. Particles are accelerated in the MAD by magnetic reconnection and stochastic turbulence acceleration. We pick out the fifteen brightest radio galaxies in the GeV band from the Fermi 4LAC-DR2 catalog… ▽ More The origins of the GeV gamma-rays from nearby radio galaxies are unknown. Hadronic emission from magnetically arrested disks (MADs) around central black holes (BHs) is proposed as a possible scenario. Particles are accelerated in the MAD by magnetic reconnection and stochastic turbulence acceleration. We pick out the fifteen brightest radio galaxies in the GeV band from the Fermi 4LAC-DR2 catalog and apply the MAD model. We find that we can explain the data in the GeV bands by the MAD model if the accretion rate is lower than 0.1% of the Eddington rate. For a higher accretion rate, GeV gamma-rays are absorbed by two-photon interaction due to copious low-energy photons. If we assume another proposed prescription of the electron heating rate by magnetic reconnection, the MAD model fails to reproduce the GeV data for the majority of our sample. This indicates that the electron heating rate is crucial. We also apply the MAD model to Sgr A* and find that GeV gamma-rays observed at the Galactic center do not come from the MAD of Sgr A*. We estimate the cosmic ray intensity from Sgr A*, but it is too low to explain the high-energy cosmic ray intensity on Earth. △ Less

Submitted 21 July, 2022; v1 submitted 19 May, 2022; originally announced May 2022.

Comments: 15 pages, 10 figures, accepted for publication in ApJ

Interaction of a Relativistic Magnetized Collisionless Shock with a Dense Clump

Authors: Sara Tomita, Yutaka Ohira, Shigeo S. Kimura, Kengo Tomida, Kenji Toma

Abstract: The interactions between a relativistic magnetized collisionless shock and dense clumps have been expected to play a crucial role on the magnetic field amplification and cosmic-ray acceleration. We investigate this process by two-dimensional Particle-In-Cell (PIC) simulations for the first time, where the clump size is much larger than the gyroradius of downstream particles. We also perform relati… ▽ More The interactions between a relativistic magnetized collisionless shock and dense clumps have been expected to play a crucial role on the magnetic field amplification and cosmic-ray acceleration. We investigate this process by two-dimensional Particle-In-Cell (PIC) simulations for the first time, where the clump size is much larger than the gyroradius of downstream particles. We also perform relativistic magnetohydrodynamic (MHD) simulations for the same condition to see the kinetic effects. We find that particles escape from the shocked clump along magnetic field lines in the PIC simulations, so that the vorticity is lower than that in the MHD simulations. Moreover, in both the PIC and MHD simulations, the shocked clump quickly decelerates because of relativistic effects. Owing to the escape and the deceleration, the shocked clump cannot amplify the downstream magnetic field in relativistic collisionless shocks. This large-scale PIC simulation opens a new window to understand large-scale behaviors in collisionless plasma systems. △ Less

Submitted 15 August, 2022; v1 submitted 18 April, 2022; originally announced April 2022.

Comments: Accepted for publication in The Astrophysical Journal Letters

Polarimetric studies of GRBs, AGN jets, and axion dark matter

Authors: Kenji Toma

Abstract: Relativistic jets are collimated outflows with speeds close to light speed, which are associated with gamma-ray bursts (GRBs), active galactic nuclei (AGNs), and so on. This article mainly overviews recent developments of polarimetric studies of GRBs and their afterglows in the gamma-ray and optical wavebands as well as the first detections of their radio polarization. Polarimetric observations an… ▽ More Relativistic jets are collimated outflows with speeds close to light speed, which are associated with gamma-ray bursts (GRBs), active galactic nuclei (AGNs), and so on. This article mainly overviews recent developments of polarimetric studies of GRBs and their afterglows in the gamma-ray and optical wavebands as well as the first detections of their radio polarization. Polarimetric observations and theoretical modelings can address the emission mechanism, magnetic field structure, and energetics of GRB jets and related collisionless plasma physics. Some of the discussed key physics are common with AGN jets. Furthermore, we mention that polarimetry of AGN jets and protoplanetary disks may be a novel approach to search for ultra-light axion dark matter. △ Less

Submitted 27 December, 2021; originally announced December 2021.

Comments: 8 pages, 5 figures, to appear in proceedings of IAU symposium 360 "Astronomical Polarimetry -- New Era of Multi-Wavelength Polarimetry", held at Mar 22-26, 2021

The Variability of the Black-Hole Image in M87 at the Dynamical Time Scale

Authors: Kaushik Satapathy, Dimitrios Psaltis, Feryal Ozel, Lia Medeiros, Sean T. Dougall, Chi-kwan Chan, Maciek Wielgus, Ben S. Prather, George N. Wong, Charles F. Gammie, Kazunori Akiyama, Antxon Alberdi, Walter Alef, Juan Carlos Algaba, Richard Anantua, Keiichi Asada, Rebecca Azulay, Anne-Kathrin Baczko, David R. Ball, Mislav Baloković, John Barrett, Bradford A. Benson, Dan Bintley, Lindy Blackburn, Raymond Blundell , et al. (213 additional authors not shown)

Abstract: The black-hole images obtained with the Event Horizon Telescope (EHT) are expected to be variable at the dynamical timescale near their horizons. For the black hole at the center of the M87 galaxy, this timescale (5-61 days) is comparable to the 6-day extent of the 2017 EHT observations. Closure phases along baseline triangles are robust interferometric observables that are sensitive to the expect… ▽ More The black-hole images obtained with the Event Horizon Telescope (EHT) are expected to be variable at the dynamical timescale near their horizons. For the black hole at the center of the M87 galaxy, this timescale (5-61 days) is comparable to the 6-day extent of the 2017 EHT observations. Closure phases along baseline triangles are robust interferometric observables that are sensitive to the expected structural changes of the images but are free of station-based atmospheric and instrumental errors. We explored the day-to-day variability in closure phase measurements on all six linearly independent non-trivial baseline triangles that can be formed from the 2017 observations. We showed that three triangles exhibit very low day-to-day variability, with a dispersion of $\sim3-5^\circ$. The only triangles that exhibit substantially higher variability ($\sim90-180^\circ$) are the ones with baselines that cross visibility amplitude minima on the $u-v$ plane, as expected from theoretical modeling. We used two sets of General Relativistic magnetohydrodynamic simulations to explore the dependence of the predicted variability on various black-hole and accretion-flow parameters. We found that changing the magnetic field configuration, electron temperature model, or black-hole spin has a marginal effect on the model consistency with the observed level of variability. On the other hand, the most discriminating image characteristic of models is the fractional width of the bright ring of emission. Models that best reproduce the observed small level of variability are characterized by thin ring-like images with structures dominated by gravitational lensing effects and thus least affected by turbulence in the accreting plasmas. △ Less

Submitted 1 November, 2021; originally announced November 2021.

Comments: Accepted for Publication in ApJ

The response of black hole spark gaps to external changes: A production mechanism of rapid TeV flares?

Authors: Shota Kisaka, Amir Levinson, Kenji Toma, Idan Niv

Abstract: We study the response of a starved Kerr black hole magnetosphere to abrupt changes in the intensity of disk emission and in the global magnetospheric current, by means of 1D general relativistic particle-in-cell simulations. Such changes likely arise from the intermittency of the accretion process. We find that in cases where the pair production opacity contributed by the soft disk photons is mode… ▽ More We study the response of a starved Kerr black hole magnetosphere to abrupt changes in the intensity of disk emission and in the global magnetospheric current, by means of 1D general relativistic particle-in-cell simulations. Such changes likely arise from the intermittency of the accretion process. We find that in cases where the pair production opacity contributed by the soft disk photons is modest, as in, e.g., M87, such changes can give rise to delayed, strong TeV flares, dominated by curvature emission of particles accelerated in the gap. The flare rise time, and the delay between the external variation and the onset of the flare emitted from the outer gap boundary, are of the order of the light crossing time of the gap. The rapid, large amplitude TeV flares observed in M87 and, perhaps, other AGNs may be produced by such a mechanism. △ Less

Submitted 2 November, 2021; v1 submitted 6 August, 2021; originally announced August 2021.

Comments: 7 pages, 4 figures, 1 table. Accepted for publication in Astrophysical Journal

Backreaction of Mass and Angular Momentum Accretion on Black Holes: General Formulation of the Metric Perturbations and Application to the Blandford-Znajek Process

Authors: Masashi Kimura, Tomohiro Harada, Atsushi Naruko, Kenji Toma

Abstract: We study the metric backreaction of mass and angular momentum accretion on black holes. We first develop the formalism of monopole and dipole linear gravitational perturbations around the Schwarzschild black holes in the Eddington-Finkelstein coordinates against the generic time-dependent matters. We derive the relation between the time dependence of the mass and angular momentum of the black hole… ▽ More We study the metric backreaction of mass and angular momentum accretion on black holes. We first develop the formalism of monopole and dipole linear gravitational perturbations around the Schwarzschild black holes in the Eddington-Finkelstein coordinates against the generic time-dependent matters. We derive the relation between the time dependence of the mass and angular momentum of the black hole and the energy-momentum tensors of accreting matters. As a concrete example, we apply our formalism to the Blandford-Znajek process around the slowly rotating black holes. We find that the time dependence of the monopole and dipole perturbations can be interpreted as the slowly rotating Kerr metric with time-dependent mass and spin parameters, which are determined from the energy and angular momentum extraction rates of the Blandford-Znajek process. We also show that the Komar angular momentum and the area of the apparent horizon are decreasing and increasing in time, respectively, while they are consistent with the Blandford-Znajek argument of energy extraction in term of black hole mechanics if we regard the time-dependent mass parameter as the energy of the black hole. △ Less

Submitted 31 July, 2021; v1 submitted 12 May, 2021; originally announced May 2021.

Comments: 35 pages, v2: minor revisions, v3: minor revisions, accepted for publication in PTEP

Report number: RUP-21-7, YITP-21-43

Journal ref: PTEP 2021 (2021) 9, 093E03

The Polarized Image of a Synchrotron Emitting Ring of Gas Orbiting a Black Hole

Authors: Ramesh Narayan, Daniel C. M. Palumbo, Michael D. Johnson, Zachary Gelles, Elizabeth Himwich, Dominic O. Chang, Angelo Ricarte, Jason Dexter, Charles F. Gammie, Andrew A. Chael, The Event Horizon Telescope Collaboration, :, Kazunori Akiyama, Antxon Alberdi, Walter Alef, Juan Carlos Algaba, Richard Anantua, Keiichi Asada, Rebecca Azulay, Anne-Kathrin Baczko, David Ball, Mislav Balokovic, John Barrett, Bradford A. Benson, Dan Bintley , et al. (215 additional authors not shown)

Abstract: Synchrotron radiation from hot gas near a black hole results in a polarized image. The image polarization is determined by effects including the orientation of the magnetic field in the emitting region, relativistic motion of the gas, strong gravitational lensing by the black hole, and parallel transport in the curved spacetime. We explore these effects using a simple model of an axisymmetric, equ… ▽ More Synchrotron radiation from hot gas near a black hole results in a polarized image. The image polarization is determined by effects including the orientation of the magnetic field in the emitting region, relativistic motion of the gas, strong gravitational lensing by the black hole, and parallel transport in the curved spacetime. We explore these effects using a simple model of an axisymmetric, equatorial accretion disk around a Schwarzschild black hole. By using an approximate expression for the null geodesics derived by Beloborodov (2002) and conservation of the Walker-Penrose constant, we provide analytic estimates for the image polarization. We test this model using currently favored general relativistic magnetohydrodynamic simulations of M87*, using ring parameters given by the simulations. For a subset of these with modest Faraday effects, we show that the ring model broadly reproduces the polarimetric image morphology. Our model also predicts the polarization evolution for compact flaring regions, such as those observed from Sgr A* with GRAVITY. With suitably chosen parameters, our simple model can reproduce the EVPA pattern and relative polarized intensity in Event Horizon Telescope images of M87*. Under the physically motivated assumption that the magnetic field trails the fluid velocity, this comparison is consistent with the clockwise rotation inferred from total intensity images. △ Less

Submitted 13 May, 2021; v1 submitted 4 May, 2021; originally announced May 2021.

Comments: 29 pages, 11 figures, published in ApJ on May 3

Journal ref: ApJ 912 35 (2021)

Matter density distribution of general relativistic highly magnetized jets driven by black holes

Authors: Taiki Ogihara, Takumi Ogawa, Kenji Toma

Abstract: High-resolution very long baseline interferometry (VLBI) radio observations have resolved the detailed emission structures of active galactic nucleus jets. General relativistic magnetohydrodynamic (GRMHD) simulations have improved the understanding of jet production physics, although theoretical studies still have difficulties in constraining the origin and distribution of jetted matter. We constr… ▽ More High-resolution very long baseline interferometry (VLBI) radio observations have resolved the detailed emission structures of active galactic nucleus jets. General relativistic magnetohydrodynamic (GRMHD) simulations have improved the understanding of jet production physics, although theoretical studies still have difficulties in constraining the origin and distribution of jetted matter. We construct a new steady, axisymmetric GRMHD jet model to obtain approximate solutions of black hole (BH) magnetospheres, and examine the matter density distribution of jets. By assuming fixed poloidal magnetic field shapes that mimic force-free analytic solutions and GRMHD simulation results and assuming constant poloidal velocity at the separation surface, which divides the inflow and outflow, we numerically solve the force-balance between the field lines at the separation surface and analytically solve the distributions of matter velocity and density along the field lines. We find that the densities at the separation surface in our parabolic field models roughly follow $\propto r_{ss}^{-2}$ in the far zone from the BH, where $r_{ss}$ is the radius of the separation surface. When the BH spin is larger or the velocity at the separation surface is smaller, the density at the separation surface becomes concentrated more near the jet edge. Our semi-analytic model, combined with radiative transfer calculations, may help interpret the high-resolution VLBI observations and understand the origin of jetted matter. △ Less

Submitted 16 February, 2021; originally announced February 2021.

Comments: 15 pages, 10 figures, 1 table, Accepted for publication in ApJ

A Jet-Bases Emission Model of the EHT 2017 Image of M87*

Authors: Tomohisa Kawashima, Kenji Toma, Motoki Kino, Kazunori Akiyama, Masanori Nakamura, Kotaro Moriyama

Abstract: We carry out general relativistic ray-tracing radiative-transfer calculations to study whether a localized emission from $e^{\pm}$ plasma rings created at the stagnation surface in the jet funnel, to which we refer as stagnation rings, can explain the ring image of M87 observed by Event Horizon Telescope (EHT) 2017. The resultant images consist of the direct image of the stagnation rings and the r… ▽ More We carry out general relativistic ray-tracing radiative-transfer calculations to study whether a localized emission from $e^{\pm}$ plasma rings created at the stagnation surface in the jet funnel, to which we refer as stagnation rings, can explain the ring image of M87 observed by Event Horizon Telescope (EHT) 2017. The resultant images consist of the direct image of the stagnation rings and the ring images formed via the strong deflection by the black-hole (BH) gravity, to which we refer as "quasi photon-ring". For the model with the BH spin $a_* = 0.99$, the direct image of the counter-jet and quasi photon-ring are almost coincident to the photon ring with diameter $\sim 40 μ{\rm as}$, while the approaching jet shows the small ring-image inside them. The synthetic observation image assuming the EHT 2017 array is consistent with that observed in M87, because the array is a bit sparse to detect the inner ring image. This indicates that the ring image in M87 might contain the important feature of the jet bases in addition to the photon ring. We find that forthcoming EHT observations can resolve the stagnation-ring image and may enable us to explore the plasma-injection mechanism into the jet funnel. △ Less

Submitted 21 December, 2020; v1 submitted 18 September, 2020; originally announced September 2020.

Comments: 8 pages, 3 figures, accepted for publication in ApJ

Probing Particle Acceleration through Broadband Early Afterglow Emission of MAGIC Gamma-Ray Burst GRB 190114C

Authors: Katsuaki Asano, Kohta Murase, Kenji Toma

Abstract: Major Atmospheric Gamma Imaging Cherenkov Telescopes (MAGIC) detected the gamma-ray afterglow of GRB 190114C, which can constrain microscopic parameters of the shock-heated plasma emitting non-thermal emission. Focusing on the early afterglow of this event, we numerically simulate the spectrum and multi-wavelength light curves with constant and wind-like circumstellar medium using a time-dependent… ▽ More Major Atmospheric Gamma Imaging Cherenkov Telescopes (MAGIC) detected the gamma-ray afterglow of GRB 190114C, which can constrain microscopic parameters of the shock-heated plasma emitting non-thermal emission. Focusing on the early afterglow of this event, we numerically simulate the spectrum and multi-wavelength light curves with constant and wind-like circumstellar medium using a time-dependent code. Our results show that the electron acceleration timescale at the highest energies is likely shorter than 20 times the gyroperiod to reproduce the GeV gamma-ray flux and its spectral index reported by {\it Fermi}. This gives an interesting constraint on the acceleration efficiency for Weibel-mediated shocks. We also constrain the number fraction of non-thermal electrons $f_{\rm e}$, and the temperature of the thermal electrons. The early optical emission can be explained by the thermal synchrotron emission with $f_{\rm e} \lesssim 0.01$. On the other hand, the X-ray light curves restrict efficient energy transfer from protons to the thermal electrons, and $f_{\rm e}\sim1$ is required if the energy fraction of the thermal electrons is larger than $\sim10$\%. The parameter constraints obtained in this work give important clues to probing plasma physics with relativistic shocks. △ Less

Submitted 24 November, 2020; v1 submitted 13 July, 2020; originally announced July 2020.

Comments: 9 pages, 10figures, ApJ in press

Journal ref: Astrophys.J. 905 (2020) 105

Comprehensive analysis of magnetospheric gaps around Kerr black holes using 1D GRPIC simulations

Authors: Shota Kisaka, Amir Levinson, Kenji Toma

Abstract: Spark gaps are likely the source of plasma in active black hole (BH) magnetospheres. In this paper, we present results of 1D general-relativistic particle-in-cell simulations of a starved BH magnetosphere with a realistic treatment of inverse Compton scattering and pair production, for a broad range of conditions, run times longer than in previous studies, and different setups. We find that follow… ▽ More Spark gaps are likely the source of plasma in active black hole (BH) magnetospheres. In this paper, we present results of 1D general-relativistic particle-in-cell simulations of a starved BH magnetosphere with a realistic treatment of inverse Compton scattering and pair production, for a broad range of conditions, run times longer than in previous studies, and different setups. We find that following the initial discharge, the system undergoes gradual evolution over prolonged time until either, restoring the vacuum state or reaching a state of quasi-periodic oscillations, depending on the spectral shape and luminosity of the ambient radiation field. The oscillations occur near the null charge surface in cases where the global magnetospheric current is in the direction defined by the product of the asymptotic Goldreich-Julian charge density and the radial velocity, while they occur near the boundary of the simulation box when it is the opposite direction (return current). Their amplitude and the resultant luminosity of TeV photons emitted from the gap depend sensitively on the conditions; for the cases studied here the ratio of TeV luminosity to the Blandford-Znajek power ranges from $10^{-5}$ to $10^{-2}$, suggesting that strong flares may be generated by moderate changes in disk emission. We also examined the dependence of the solution on the initial number of particles per cell (PPC) and found convergence for PPC of about 50 for the cases examined. At lower PPC values the pair multiplicity is found to be artificially high, affecting the solution considerably. △ Less

Submitted 31 August, 2020; v1 submitted 6 July, 2020; originally announced July 2020.

Comments: 12 pages, 9 figures, 1 table. Accepted for publication in Astrophysical Journal

Multi-Wave band Synchrotron Polarization of Gamma-Ray Burst Afterglows

Authors: Jiro Shimoda, Kenji Toma

Abstract: Multi-wave band synchrotron linear polarization of gamma-ray burst (GRB) afterglows is studied under the assumption of an anisotropic turbulent magnetic field with a coherence length of the plasma skin-depth scale in the downstream of forward shocks. We find that for typical GRBs, in comparison to the optical polarization, the degree of radio polarization shows a similar temporal evolution but a s… ▽ More Multi-wave band synchrotron linear polarization of gamma-ray burst (GRB) afterglows is studied under the assumption of an anisotropic turbulent magnetic field with a coherence length of the plasma skin-depth scale in the downstream of forward shocks. We find that for typical GRBs, in comparison to the optical polarization, the degree of radio polarization shows a similar temporal evolution but a significantly smaller peak value. This results from differences in observed intensity image shapes between the radio and optical bands. We also show that the degree of the polarization spectrum undergoes a gradual variation from the low- to the high-polarization regime above the intensity of the spectral peak frequency, and that the difference in polarization angles in the two regimes is zero or 90 degrees. Thus, simultaneous multi-wave band polarimetric observations of GRB afterglows would be a new determinative test of the plasma-scale magnetic field model. We also discuss theoretical implications from the recent detection of radio linear polarization in GRB 171205A with ALMA and other models of magnetic field configuration. △ Less

Submitted 27 May, 2021; v1 submitted 7 May, 2020; originally announced May 2020.

Comments: 16 pages, 12 figures, revised version of 2, accepted by ApJ in publication

SYMBA: An end-to-end VLBI synthetic data generation pipeline

Authors: F. Roelofs, M. Janssen, I. Natarajan, R. Deane, J. Davelaar, H. Olivares, O. Porth, S. N. Paine, K. L. Bouman, R. P. J. Tilanus, I. M. van Bemmel, H. Falcke, K. Akiyama, A. Alberdi, W. Alef, K. Asada, R. Azulay, A. Baczko, D. Ball, M. Baloković, J. Barrett, D. Bintley, L. Blackburn, W. Boland, G. C. Bower , et al. (183 additional authors not shown)

Abstract: Realistic synthetic observations of theoretical source models are essential for our understanding of real observational data. In using synthetic data, one can verify the extent to which source parameters can be recovered and evaluate how various data corruption effects can be calibrated. These studies are important when proposing observations of new sources, in the characterization of the capabili… ▽ More Realistic synthetic observations of theoretical source models are essential for our understanding of real observational data. In using synthetic data, one can verify the extent to which source parameters can be recovered and evaluate how various data corruption effects can be calibrated. These studies are important when proposing observations of new sources, in the characterization of the capabilities of new or upgraded instruments, and when verifying model-based theoretical predictions in a comparison with observational data. We present the SYnthetic Measurement creator for long Baseline Arrays (SYMBA), a novel synthetic data generation pipeline for Very Long Baseline Interferometry (VLBI) observations. SYMBA takes into account several realistic atmospheric, instrumental, and calibration effects. We used SYMBA to create synthetic observations for the Event Horizon Telescope (EHT), a mm VLBI array, which has recently captured the first image of a black hole shadow. After testing SYMBA with simple source and corruption models, we study the importance of including all corruption and calibration effects. Based on two example general relativistic magnetohydrodynamics (GRMHD) model images of M87, we performed case studies to assess the attainable image quality with the current and future EHT array for different weather conditions. The results show that the effects of atmospheric and instrumental corruptions on the measured visibilities are significant. Despite these effects, we demonstrate how the overall structure of the input models can be recovered robustly after performing calibration steps. With the planned addition of new stations to the EHT array, images could be reconstructed with higher angular resolution and dynamic range. In our case study, these improvements allowed for a distinction between a thermal and a non-thermal GRMHD model based on salient features in reconstructed images. △ Less

Submitted 2 April, 2020; originally announced April 2020.

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

Hadronic High-energy Emission from Magnetically Arrested Disks in Radio Galaxies

Authors: Shigeo S. Kimura, Kenji Toma

Abstract: We propose a novel interpretation that gamma-rays from nearby radio galaxies are hadronic emission from magnetically arrested disks (MADs) around central black holes (BHs). The magnetic energy in MADs is higher than the thermal energy of the accreting plasma, where the magnetic reconnection or turbulence may efficiently accelerate non-thermal protons. They emit gamma-rays via hadronic processes, w… ▽ More We propose a novel interpretation that gamma-rays from nearby radio galaxies are hadronic emission from magnetically arrested disks (MADs) around central black holes (BHs). The magnetic energy in MADs is higher than the thermal energy of the accreting plasma, where the magnetic reconnection or turbulence may efficiently accelerate non-thermal protons. They emit gamma-rays via hadronic processes, which can account for the observed gamma-rays for M87 and NGC 315. Non-thermal electrons are also accelerated with protons and produce MeV gamma-rays, which is useful to test our model by proposed MeV satellites. The hadronic emission from the MADs may significantly contribute to the GeV gamma-ray background and produce the multi-PeV neutrino background detectable by IceCube-Gen2. In addition, gamma-rays from MADs provide electron-positron pairs through two-photon pair production at the BH magnetosphere. These pairs can screen the vacuum gap, which affects high-energy emission and jet-launching mechanisms in radio galaxies. △ Less

Submitted 6 January, 2021; v1 submitted 29 March, 2020; originally announced March 2020.

Comments: 14 pages, 5 figures, 1 table, published in ApJ. Fixed typos

Journal ref: Astrophys.J. 905 (2020) 178

Efficient acceleration of cylindrical jets: effects of radiative cooling and tangled magnetic field

Authors: Shuta J. Tanaka, Kenji Toma

Abstract: Diverging supersonic flows are accelerating, as in the case of a de Laval nozzle, and the same concept has been applied for acceleration of magnetohydrodynamic flows in the universe. Here, we study the dynamics of "non-diverging" cylindrical supersonic flows and show that they can be accelerated by effects of radiative cooling and the tangled magnetic field. In addition to radiative cooling of the… ▽ More Diverging supersonic flows are accelerating, as in the case of a de Laval nozzle, and the same concept has been applied for acceleration of magnetohydrodynamic flows in the universe. Here, we study the dynamics of "non-diverging" cylindrical supersonic flows and show that they can be accelerated by effects of radiative cooling and the tangled magnetic field. In addition to radiative cooling of the jet materials (cooling effect), conversion of the ordered magnetic field into the turbulent one (conversion effect) and dissipation of the turbulent magnetic field (dissipation effect) are formulated according to our study on pulsar wind nebulae. Although each of the cooling and conversion effects is an ineffective acceleration process, the terminal velocity of magnetized cylindrical jets attains about half of the maximum possible value when the cooling, conversion and dissipation effects work simultaneously. The radiation efficiency is also about half of the total luminosity of the jet in the case of maximal acceleration. The concept for flow acceleration by the non-ideal MHD effects may be useful for studying relativistic jets in active galactic nuclei, in which the region near the jet axis is expected to be cylindrical and kink unstable. △ Less

Submitted 13 March, 2020; originally announced March 2020.

Comments: Accepted for publication in MNRAS, 12 pages, 4 figures, 1 table

ALMA Polarimetry of AT2018cow

Authors: K. Y. Huang, J. Shimoda, Y. Urata, K. Toma, K. Yamaoka, K. Asada, H. Nagai, S. Takahashi, G. Petitpas, M. Tashiro

Abstract: We present the first radio polarimetric observations of a fast-rising blue optical transient, AT2018cow. Two epochs of polarimetry with additional coincident photometry were performed with the Atacama Large Millimeter/submillimeter Array (ALMA). The overall photometric results based on simultaneous observations in the 100 and 230 GHz bands are consistent with the non-thermal radiation model report… ▽ More We present the first radio polarimetric observations of a fast-rising blue optical transient, AT2018cow. Two epochs of polarimetry with additional coincident photometry were performed with the Atacama Large Millimeter/submillimeter Array (ALMA). The overall photometric results based on simultaneous observations in the 100 and 230 GHz bands are consistent with the non-thermal radiation model reported by Ho et al. (2019) and indicate that the spectral peaks ($\sim110$ GHz at the first epoch and $\sim67$ GHz at the second epoch) represent the synchrotron self-absorption frequency. The non-detection of linear polarization with $<$0.15% in the 230 GHz band at the phase when the effect of synchrotron self-absorption was quite small in the band may be explained by internal Faraday depolarization with high circumburst density and strong magnetic field. This result supports the stellar explosion scenario rather than the tidal disruption model. The maximum energy of accelerating particles at the shocks of AT2018cow-like objects is also discussed. △ Less

Submitted 23 May, 2019; originally announced May 2019.

Comments: 7 pages, 3 figures, 1 table, accepted for publication in ApJL

First Detection of Radio Linear Polarization in a Gamma Ray Burst Afterglow

Authors: Y. Urata, K. Toma, K. Huang, K. Asada, H. Nagai, S. Takahashi, G. Petitpas, M. Tashiro, K. Yamaoka

Abstract: We report the first detection of radio polarization of a GRB afterglow with the first intensive combined use of telescopes in the millimeter and submillimeter ranges for GRB171205A. The linear polarization degree in the millimeter band at the sub-percent level ($0.27 \pm 0.04\%$) is lower than those observed in late-time optical afterglows (weighted average of $\sim 1\%$). The Faraday depolarizati… ▽ More We report the first detection of radio polarization of a GRB afterglow with the first intensive combined use of telescopes in the millimeter and submillimeter ranges for GRB171205A. The linear polarization degree in the millimeter band at the sub-percent level ($0.27 \pm 0.04\%$) is lower than those observed in late-time optical afterglows (weighted average of $\sim 1\%$). The Faraday depolarization by non-accelerated, cool electrons in the shocked region is one of possible mechanisms for the low value. In this scenario, larger total energy by a factor of $\sim 10$ than ordinary estimate without considering non-accelerated electrons is required. The polarization position angle varies by at least 20 degrees across the millimeter band, which is not inconsistent with this scenario. This result indicates that polarimetry in the millimeter and submillimeter ranges is a unique tool for investigating GRB energetics, and coincident observations with multiple frequencies or bands would provide more accurate measurements of the non-accelerated electron fraction. △ Less

Submitted 29 September, 2019; v1 submitted 17 April, 2019; originally announced April 2019.

Comments: Accepted to ApJ Letters; 8 pages, 4 figures, 1 table,

A Mechanism for Triple-Ridge Emission Structure of AGN Jets

Authors: Taiki Ogihara, Kazuya Takahashi, Kenji Toma

Abstract: Recent radio VLBI observations of the relativistic jet in M87 radio galaxy have shown a triple-ridge structure that consists of the conventional limb-brightened feature and a central narrow ridge. Motivated by these observations, we examine a steady axisymmetric force-free model of a jet driven by the central black hole (BH) with its electromagnetic structure being consistent with general relativi… ▽ More Recent radio VLBI observations of the relativistic jet in M87 radio galaxy have shown a triple-ridge structure that consists of the conventional limb-brightened feature and a central narrow ridge. Motivated by these observations, we examine a steady axisymmetric force-free model of a jet driven by the central black hole (BH) with its electromagnetic structure being consistent with general relativistic magnetohydrodynamic simulations, and find that it can produce triple-ridge images even if we assume a simple Gaussian distribution of emitting electrons at the base of the jet. We show that the fluid velocity field associated with the electromagnetic field produces the central ridge component due to the relativistic beaming effect, while the limb-brightened feature arises due to strong magnetic field around the jet edge which also induces the electrons to be dense there. We also show that the computed image strongly depends on the electromagnetic field structure, viewing angle, and parameters related to the electrons' spatial distribution at the jet base. This study will help constraining the non-thermal electron injection mechanism of BH jets and be complementary to theoretical analyses of the upcoming data of Event Horizon Telescope. △ Less

Submitted 15 April, 2019; originally announced April 2019.

Comments: 15 Pages, 8 figures, accepted for publication in the Astrophysical Journal

Hunting Axion Dark Matter with Protoplanetary Disk Polarimetry

Authors: Tomohiro Fujita, Ryo Tazaki, Kenji Toma

Abstract: We find that the polarimetric observations of protoplanetary disks are useful to search for ultra-light axion dark matter. Axion dark matter predicts the rotation of the linear polarization plane of propagating light, and protoplanetary disks are ideal targets to observe it. We show that a recent observation puts the tightest constraint on the axion-photon coupling constant for axion mass… ▽ More We find that the polarimetric observations of protoplanetary disks are useful to search for ultra-light axion dark matter. Axion dark matter predicts the rotation of the linear polarization plane of propagating light, and protoplanetary disks are ideal targets to observe it. We show that a recent observation puts the tightest constraint on the axion-photon coupling constant for axion mass $m\lesssim10^{-21}$eV. △ Less

Submitted 1 May, 2019; v1 submitted 8 November, 2018; originally announced November 2018.

Comments: 5 pages, 3 figures, Updated to match the version accepted by PRL

Journal ref: Phys. Rev. Lett. 122, 191101 (2019)

Parabolic Jets from the Spinning Black Hole in M87

Authors: Masanori Nakamura, Keiichi Asada, Kazuhiro Hada, Hung-Yi Pu, Scott Noble, Chihyin Tseng, Kenji Toma, Motoki Kino, Hiroshi Nagai, Kazuya Takahashi, Juan-Carlos Algaba, Monica Orienti, Kazunori Akiyama, Akihiro Doi, Gabriele Giovannini, Marcello Giroletti, Mareki Honma, Shoko Koyama, Rocco Lico, Kotaro Niinuma, Fumie Tazaki

Abstract: The M87 jet is extensively examined by utilizing general relativistic magnetohydrodynamic (GRMHD) simulations as well as the steady axisymmetric force-free electrodynamic (FFE) solution. Quasi-steady funnel jets are obtained in GRMHD simulations up to the scale of $\sim 100$ gravitational radius ($r_{\rm g}$) for various black hole (BH) spins. As is known, the funnel edge is approximately determin… ▽ More The M87 jet is extensively examined by utilizing general relativistic magnetohydrodynamic (GRMHD) simulations as well as the steady axisymmetric force-free electrodynamic (FFE) solution. Quasi-steady funnel jets are obtained in GRMHD simulations up to the scale of $\sim 100$ gravitational radius ($r_{\rm g}$) for various black hole (BH) spins. As is known, the funnel edge is approximately determined by the following equipartitions; i) the magnetic and rest-mass energy densities and ii) the gas and magnetic pressures. Our numerical results give an additional factor that they follow the outermost parabolic streamline of the FFE solution, which is anchored to the event horizon on the equatorial plane. We also identify the matter dominated, non-relativistic corona/wind play a dynamical role in shaping the funnel jet into the parabolic geometry. We confirm a quantitative overlap between the outermost parabolic streamline of the FFE jet and the edge of jet sheath in VLBI observations at $\sim 10^{1}$-$10^{5} \, r_{\rm g}$, suggesting that the M87 jet is likely powered by the spinning BH. Our GRMHD simulations also indicate a lateral stratification of the bulk acceleration (i.e., the spine-sheath structure) as well as an emergence of knotty superluminal features. The spin characterizes the location of the jet stagnation surface inside the funnel. We suggest that the limb-brightened feature could be associated with the nature of the BH-driven jet, if the Doppler beaming is a dominant factor. Our findings can be examined with (sub-)mm VLBI observations, giving a clue for the origin of the M87 jet. △ Less

Submitted 23 October, 2018; originally announced October 2018.

Comments: 29 pages, 23 figures, accepted for publication in ApJ

The origin of polarization in kilonovae and the case of the gravitational-wave counterpart AT 2017gfo

Authors: M. Bulla, S. Covino, K. Kyutoku, M. Tanaka, J. R. Maund, F. Patat, K. Toma, K. Wiersema, J. Bruten, Z. P. Jin, V. Testa

Abstract: The Gravitational Wave (GW) event GW 170817 was generated by the coalescence of two neutron stars (NS) and produced an electromagnetic transient, labelled AT 2017gfo, that was target of a massive observational campaign. Polarimetry, a powerful diagnostic tool for probing the geometry and emission processes of unresolved sources, was obtained for this event. The observed linear polarization was con… ▽ More The Gravitational Wave (GW) event GW 170817 was generated by the coalescence of two neutron stars (NS) and produced an electromagnetic transient, labelled AT 2017gfo, that was target of a massive observational campaign. Polarimetry, a powerful diagnostic tool for probing the geometry and emission processes of unresolved sources, was obtained for this event. The observed linear polarization was consistent with being mostly induced by intervening dust, suggesting that the intrinsic emission was weakly polarized ($P < 0.4-0.5$ %). In this paper, we present and discuss a detailed analysis of the linear polarization expected from a merging NS binary system by means of 3D Monte Carlo radiative transfer simulations assuming a range of possible configurations, wavelengths, epochs and viewing angles. We find that polarization originates from the non-homogeneous opacity distribution within the ejecta and can reach levels of $P\sim1$ % at early times (1$-$2 days after the merger) and in the optical R band. Smaller polarization signals are expected at later epochs and/or different wavelengths. From the viewing-angle dependence of the polarimetric signal, we constrain the observer orientation of AT 2017gfo within $\sim$65$^\circ$ from the polar direction. The detection of non-zero polarization in future events will unambiguously reveal the presence of a lanthanide-free ejecta component and unveil its spatial and angular distribution. △ Less

Submitted 20 November, 2018; v1 submitted 11 September, 2018; originally announced September 2018.

Comments: published in Nature Astronomy, www.nature.com/articles/s41550-018-0593-y/

Confinement of the Crab Nebula with tangled magnetic field by its supernova remnant

Authors: Shuta J. Tanaka, Kenji Toma, Nozomu Tominaga

Abstract: A pulsar wind is a relativistic outflow dominated by Poynting energy at its base. Based on the standard ideal magnetohydrodynamic (MHD) model of pulsar wind nebulae (PWNe) with the ordered magnetic field, the observed slow expansion $v_{\rm PWN} \ll c$ requires the wind to be dominated by kinetic energy at the upstream of its termination shock, which conflicts with the pulsar wind theory ($σ$-prob… ▽ More A pulsar wind is a relativistic outflow dominated by Poynting energy at its base. Based on the standard ideal magnetohydrodynamic (MHD) model of pulsar wind nebulae (PWNe) with the ordered magnetic field, the observed slow expansion $v_{\rm PWN} \ll c$ requires the wind to be dominated by kinetic energy at the upstream of its termination shock, which conflicts with the pulsar wind theory ($σ$-problem). In this paper, we extend the standard model of PWNe by phenomenologically taking into account conversion of the ordered to turbulent magnetic field and dissipation of the turbulent magnetic field. Disordering of the magnetic structure is inferred from the recent three-dimensional relativistic ideal MHD simulations, while magnetic dissipation is a non-ideal MHD effect requiring a finite resistivity. We apply this model to the Crab Nebula and find that the conversion effect is important for the flow deceleration, while the dissipation effect is not. Even for Poynting-dominated pulsar wind, we obtain the Crab Nebula's $v_{\rm PWN}$ by adopting a finite conversion time-scale of $\sim 0.3$ yr. Magnetic dissipation primarily affects the synchrotron radiation properties. Any values of the pulsar wind magnetization $σ_{\rm w}$ are allowed within the present model of the PWN dynamics alone, and even a small termination shock radius of $\ll 0.1$ pc reproduces the observed dynamical features of the Crab Nebula. In order to establish a high-$σ_{\rm w}$ model of PWNe, it is important to extend the present model by taking into account the broadband spectrum and its spacial profiles. △ Less

Submitted 21 May, 2018; originally announced May 2018.

Comments: Accepted for publication in MNRAS, 13 pages, 4 figures, 1 table

Fast-spinning black holes inferred from symmetrically limb-brightened radio jets

Authors: Kazuya Takahashi, Kenji Toma, Motoki Kino, Masanori Nakamura, Kazuhiro Hada

Abstract: This paper theoretically investigates the relations between the structure of relativistic jets and produced synchrotron images, by using a steady axisymmetric force-free jet model. We especially focus on the limb-brightened jets that are largely symmetric to the jet axes and observed in some active galactic nuclei such as M87, Mrk~501, Cyg~A, and 3C84. We find that symmetrically limb-brightened im… ▽ More This paper theoretically investigates the relations between the structure of relativistic jets and produced synchrotron images, by using a steady axisymmetric force-free jet model. We especially focus on the limb-brightened jets that are largely symmetric to the jet axes and observed in some active galactic nuclei such as M87, Mrk~501, Cyg~A, and 3C84. We find that symmetrically limb-brightened images can be produced when magnetic field lines of the jet penetrate a fast-spinning black hole as motivated by the Blandford-Znajek mechanism. On the other hand, jets with magnetic field lines that pass through a slowly spinning black hole or the Keplerian accretion disk produce highly asymmetric radio images. In addition, the edge of a counter jet tends to be luminous in the accretion-disk model even for rather small viewing angles, which may be problematic for some observed jets. We also suggest that the site of particle accelerations in relativistic jets can be constrained by fitting the radio images to observations. This kind of study focusing on the jet images far away from the central engine is complimentary to those concentrating directly on the innermost region with upcoming data of the Event Horizon Telescope. △ Less

Submitted 10 October, 2018; v1 submitted 1 February, 2018; originally announced February 2018.

Comments: modified after the second referee report, accepted for publication in ApJ

The unpolarized macronova associated with the gravitational wave event GW170817

Authors: S. Covino, K. Wiersema, Y. Z. Fan, K. Toma, A. B. Higgins, A. Melandri, P. D'Avanzo, C. G. Mundell, E. Palazzi, N. R. Tanvir, M. G. Bernardini, M. Branchesi, E. Brocato, S. Campana, S. di Serego Alighieri, D. Gotz, J. P. U. Fynbo, W. Gao, A. Gomboc, B. Gompertz, J. Greiner, J. Hjorth, Z. P. Jin, L. Kaper, S. Klose , et al. (15 additional authors not shown)

Abstract: The merger of two dense stellar remnants including at least one neutron star (NS) is predicted to produce gravitational waves (GWs) and short duration gamma ray bursts (GRBs). In the process, neutron-rich material is ejected from the system and heavy elements are synthesized by r-process nucleosynthesis. The radioactive decay of these heavy elements produces additional transient radiation termed "… ▽ More The merger of two dense stellar remnants including at least one neutron star (NS) is predicted to produce gravitational waves (GWs) and short duration gamma ray bursts (GRBs). In the process, neutron-rich material is ejected from the system and heavy elements are synthesized by r-process nucleosynthesis. The radioactive decay of these heavy elements produces additional transient radiation termed "kilonova" or "macronova". We report the detection of linear optical polarization P = (0.50 +/- 0.07)% at 1.46 days after detection of the GWs from GW170817, a double neutron star merger associated with an optical macronova counterpart and a short GRB. The optical emission from a macronova is expected to be characterized by a blue, rapidly decaying, component and a red, more slowly evolving, component due to material rich of heavy elements, the lanthanides. The polarization measurement was made when the macronova was still in its blue phase, during which there is an important contribution from a lanthanide-free outflow. The low degree of polarization is consistent with intrinsically unpolarized emission scattered by Galactic dust, suggesting a symmetric geometry of the emitting region and low inclination of the merger system. Stringent upper limits to the polarization degree from 2.45 - 9.48 days post-burst are consistent with the lanthanides-rich macronova interpretation. △ Less

Submitted 16 October, 2017; originally announced October 2017.

Comments: 18 pages, 1 figure, 2 tables, Nature Astronomy, in press

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.

Comments: Index of Cherenkov Telescope Array conference proceedings at the ICRC2017, Busan, Korea

Rayleigh-Taylor Instability in Two-Component Relativistic Jets

Authors: Kenji Toma, Serguei S. Komissarov, Oliver Porth

Abstract: Relativistic jets associated with active galactic nuclei and gamma-ray bursts propagate over huge distances without significant loss of momentum. At the same time they are bright emitters, which is indicative of strong energy dissipation. This points towards a mechanism of internal dissipation which does not result in a global disruption of the flow. One possibility is internal shocks and another… ▽ More Relativistic jets associated with active galactic nuclei and gamma-ray bursts propagate over huge distances without significant loss of momentum. At the same time they are bright emitters, which is indicative of strong energy dissipation. This points towards a mechanism of internal dissipation which does not result in a global disruption of the flow. One possibility is internal shocks and another one is turbulence driven by local instabilities. Such instabilities can be triggered when a freely expanding jet is reconfined by either the cocoon or external gas pressure. In this paper we study the dynamics of two-component spine-sheath hydrodynamic jets coming into pressure equilibrium with external gas using 2D computer simulations. We find that the jet oscillations lead to a rapid onset of Rayleigh-Taylor-type instabilities, which results in additional internal dissipation and mixing of the jet components. Although slightly different in details, this outcome holds both for the heavy-spine-light-sheath and light-spine-heavy-sheath configurations. The results may provide an explanation to the spatial flaring observed in some AGN jets on kpc-scales. △ Less

Submitted 13 July, 2017; v1 submitted 29 May, 2017; originally announced May 2017.

Comments: 6 pages, 4 figures; accepted for publication in MNRAS; texts revised from v1

Where is the electric current driven in the Blandford-Znajek process?

Authors: Kenji Toma, Fumio Takahara

Abstract: The Blandford-Znajek process, the steady electromagnetic energy extraction from a rotating black hole, is widely believed to work for driving relativistic jets, although it is still under debate where the electric current is driven. We address this issue analytically by investigating the time-dependent state in the Boyer-Lindquist and Kerr-Schild coordinate systems. This analysis suggests that a n… ▽ More The Blandford-Znajek process, the steady electromagnetic energy extraction from a rotating black hole, is widely believed to work for driving relativistic jets, although it is still under debate where the electric current is driven. We address this issue analytically by investigating the time-dependent state in the Boyer-Lindquist and Kerr-Schild coordinate systems. This analysis suggests that a non-ideal magnetohydrodynamic region is required in the time-dependent state, while not in the steady state. △ Less

Submitted 28 November, 2016; originally announced November 2016.

Comments: 4 pages, 2 figures, contribution to proceedings IAU Symposium No. 324, "New Frontiers in Black Hole Astrophysics", C. Mundell & A. Gomboc, eds

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.

Comments: Index of CTA conference proceedings for the Gamma 2016, Heidelberg, Germany

Evolution of an Accretion Disk in Binary Black Hole Systems

Authors: Shigeo S. Kimura, Sanemichi Z. Takahashi, Kenji Toma

Abstract: We investigate evolution of an accretion disc in binary black hole (BBH) systems and possible electromagnetic counterparts of the gravitational waves from mergers of BBHs. Perna et al. (2016) proposed a novel evolutionary scenario of an accretion disc in BBHs in which a disc eventually becomes "dead", i.e., the magnetorotational instability (MRI) becomes inactive. In their scenario, the dead disc… ▽ More We investigate evolution of an accretion disc in binary black hole (BBH) systems and possible electromagnetic counterparts of the gravitational waves from mergers of BBHs. Perna et al. (2016) proposed a novel evolutionary scenario of an accretion disc in BBHs in which a disc eventually becomes "dead", i.e., the magnetorotational instability (MRI) becomes inactive. In their scenario, the dead disc survives until {\it a few seconds before} the merger event. We improve the dead disc model and propose another scenario, taking account of effects of the tidal torque from the companion and the critical ionization degree for MRI activation more carefully. We find that the mass of the dead disc is much lower than that in the Perna's scenario. When the binary separation sufficiently becomes small, the mass inflow induced by the tidal torque reactivates MRI, restarting mass accretion onto the black hole. We also find that this disc "revival" happens {\it more than thousands of years before} the merger. The mass accretion induced by the tidal torque increases as the separation decreases, and a relativistic jet could be launched before the merger. The emissions from these jets are too faint compared to GRBs, but detectable if the merger events happen within $\lesssim 10$ Mpc or if the masses of the black holes are as massive as $\sim 10^5 M_{\odot}$. △ Less

Submitted 21 November, 2016; v1 submitted 7 July, 2016; originally announced July 2016.

Comments: 8 pages, 3 figures, accepted for publication in MNRAS

Journal ref: Mon.Not.Roy.Astron.Soc. 465 (2017) 4406

Causal production of the electromagnetic energy flux and role of the negative energies in Blandford-Znajek process

Authors: Kenji Toma, Fumio Takahara

Abstract: Blandford-Znajek process, the steady electromagnetic energy extraction from a rotating black hole (BH), is widely believed to work for driving relativistic jets in active galactic nuclei, gamma-ray bursts and Galactic microquasars, although it is still under debate how the Poynting flux is causally produced and how the rotational energy of the BH is reduced. We generically discuss the Kerr BH magn… ▽ More Blandford-Znajek process, the steady electromagnetic energy extraction from a rotating black hole (BH), is widely believed to work for driving relativistic jets in active galactic nuclei, gamma-ray bursts and Galactic microquasars, although it is still under debate how the Poynting flux is causally produced and how the rotational energy of the BH is reduced. We generically discuss the Kerr BH magnetosphere filled with a collisionless plasma screening the electric field along the magnetic field, extending the arguments of Komissarov and our previous paper, and propose a new picture for resolving the issues. For the magnetic field lines threading the equatorial plane in the ergosphere, we find that the inflow of particles with negative energy as measured in the coordinate basis is generated near that plane as a feedback from the Poynting flux production, which appears to be a similar process to the mechanical Penrose process. For the field lines threading the event horizon, we first show that the concept of the steady inflow of negative electromagnetic energy is not physically essential, partly because the sign of the electromagnetic energy density depends on the coordinates. Then we build an analytical toy model of a time-dependent process both in the Boyer-Lindquist and Kerr-Schild coordinate systems in which the force-free plasma injected continuously is filling a vacuum, and suggest that the structure of the steady outward Poynting flux is causally constructed by the displacement current and the cross-field current at the in-going boundary between the plasma and the vacuum. In the steady state, the Poynting flux is maintained without any electromagnetic source. △ Less

Submitted 11 May, 2016; originally announced May 2016.

Comments: 28 pages, 6 figures; accepted for publication in PTEP

Journal ref: PTEP, 2016, 3E01 (2016)

Gamma-ray bursts and Population III stars

Authors: Kenji Toma, Sung-Chul Yoon, Volker Bromm

Abstract: Gamma-ray bursts (GRBs) are ideal probes of the epoch of the first stars and galaxies. We review the recent theoretical understanding of the formation and evolution of the first (so-called Population III) stars, in light of their viability of providing GRB progenitors. We proceed to discuss possible unique observational signatures of such bursts, based on the current formation scenario of long GRB… ▽ More Gamma-ray bursts (GRBs) are ideal probes of the epoch of the first stars and galaxies. We review the recent theoretical understanding of the formation and evolution of the first (so-called Population III) stars, in light of their viability of providing GRB progenitors. We proceed to discuss possible unique observational signatures of such bursts, based on the current formation scenario of long GRBs. These include signatures related to the prompt emission mechanism, as well as to the afterglow radiation, where the surrounding intergalactic medium might imprint a telltale absorption spectrum. We emphasize important remaining uncertainties in our emerging theoretical framework. △ Less

Submitted 15 March, 2016; originally announced March 2016.

Comments: 27 pages, 5 figures, an output from ISSI-BJ Workshop on Gamma-Ray Bursts: a tool to explore the young Universe, accepted for publication in Space Science Reviews

Stochastic Particle Acceleration in Turbulence Generated by the Magnetorotational Instability

Authors: Shigeo S. Kimura, Kenji Toma, Takeru K. Suzuki, Shu-ichiro Inutsuka

Abstract: We investigate stochastic particle acceleration in accretion flows. It is believed that the magnetorotational instability (MRI) generates turbulence inside accretion flows and that cosmic rays (CRs) are accelerated by the turbulence. We calculate equations of motion for CRs in the turbulent fields generated by MRI with the shearing box approximation without back reaction to the field. The results… ▽ More We investigate stochastic particle acceleration in accretion flows. It is believed that the magnetorotational instability (MRI) generates turbulence inside accretion flows and that cosmic rays (CRs) are accelerated by the turbulence. We calculate equations of motion for CRs in the turbulent fields generated by MRI with the shearing box approximation without back reaction to the field. The results show that the CRs randomly gain or lose their energies through the interaction with the turbulent fields. The CRs diffuse in the configuration space anisotropically: The diffusion coefficient in direction of the unperturbed flow is about twenty times higher than the Bohm coefficient, while those in the other directions are only a few times higher than the Bohm. The momentum distribution is isotropic, and its evolution can be described by the diffusion equation in momentum space where the diffusion coefficient is a power-law function of the CR momentum. We show that the shear acceleration efficiently works for energetic particles. We also cautiously note that in the shearing box approximation, particles that cross the simulation box many times along the radial direction suffer unphysical runaway acceleration by the Lorentz transformation, which needs to be taken with special care. △ Less

Submitted 24 February, 2016; originally announced February 2016.

Comments: 11 pages, 9 figures, re-submitted to ApJ with minor revision

Journal ref: Astrophys. J, 822, 88, 2016

X-Ray Polarimetry with the Polarization Spectroscopic Telescope Array (PolSTAR)

Authors: Henric S. Krawczynski, Daniel Stern, Fiona A. Harrison, Fabian F. Kislat, Anna Zajczyk, Matthias Beilicke, Janie Hoormann, Qingzhen Guo, Ryan Endsley, Adam R. Ingram, Hiromasa Miyasaka, Kristin K. Madsen, Kim M. Aaron, Rashied Aminia, Matthew G. Baring, Banafsheh Beheshtipour, Arash Bodaghee, Jeffrey Booth, Chester Borden, Markus Boettcher, Finn E. Christensen, Paolo S. Coppi, Ramanath Cowsik, Shane Davis, Jason Dexter , et al. (27 additional authors not shown)

Abstract: This paper describes the Polarization Spectroscopic Telescope Array (PolSTAR), a mission proposed to NASA's 2014 Small Explorer (SMEX) announcement of opportunity. PolSTAR measures the linear polarization of 3-50 keV (requirement; goal: 2.5-70 keV) X-rays probing the behavior of matter, radiation and the very fabric of spacetime under the extreme conditions close to the event horizons of black hol… ▽ More This paper describes the Polarization Spectroscopic Telescope Array (PolSTAR), a mission proposed to NASA's 2014 Small Explorer (SMEX) announcement of opportunity. PolSTAR measures the linear polarization of 3-50 keV (requirement; goal: 2.5-70 keV) X-rays probing the behavior of matter, radiation and the very fabric of spacetime under the extreme conditions close to the event horizons of black holes, as well as in and around magnetars and neutron stars. The PolSTAR design is based on the technology developed for the Nuclear Spectroscopic Telescope Array (NuSTAR) mission launched in June 2012. In particular, it uses the same X-ray optics, extendable telescope boom, optical bench, and CdZnTe detectors as NuSTAR. The mission has the sensitivity to measure ~1% linear polarization fractions for X-ray sources with fluxes down to ~5 mCrab. This paper describes the PolSTAR design as well as the science drivers and the potential science return. △ Less

Submitted 28 October, 2015; originally announced October 2015.

Comments: (Astroparticle Physics in press, 34 pages, 23 figures, 6 tables)

CTA Contributions to the 34th International Cosmic Ray Conference (ICRC2015)

Authors: The CTA Consortium, :, A. Abchiche, U. Abeysekara, Ó. Abril, F. Acero, B. S. Acharya, M. Actis, G. Agnetta, J. A. Aguilar, F. Aharonian, A. Akhperjanian, A. Albert, M. Alcubierre, R. Alfaro, E. Aliu, A. J. Allafort, D. Allan, I. Allekotte, R. Aloisio, J. -P. Amans, E. Amato, L. Ambrogi, G. Ambrosi, M. Ambrosio , et al. (1290 additional authors not shown)

Abstract: List of contributions from the CTA Consortium presented at the 34th International Cosmic Ray Conference, 30 July - 6 August 2015, The Hague, The Netherlands. List of contributions from the CTA Consortium presented at the 34th International Cosmic Ray Conference, 30 July - 6 August 2015, The Hague, The Netherlands. △ Less

Submitted 11 September, 2015; v1 submitted 24 August, 2015; originally announced August 2015.

Comments: Index of CTA conference proceedings at the ICRC2015, The Hague (The Netherlands). v1: placeholder with no arXiv links yet, to be replaced once individual contributions have been all submitted; v2: final with arXiv links to all CTA contributions and full author list

Synchrotron Self-Compton Emission by Relativistic Electrons under Stochastic Acceleration: Application to Mrk 421 and Mrk 501

Authors: Jun Kakuwa, Kenji Toma, Katsuaki Asano, Masaaki Kusunose, Fumio Takahara

Abstract: We examine the applicability of the stochastic electron acceleration to two high synchrotron peaked blazars, Mrk 421 and Mrk 501, assuming synchrotron self-Compton emission of gamma-rays. Our model considers an emitting region moving at relativistic speed, where non-thermal electrons are accelerated and attain a steady-state energy spectrum together with the photons they emit. The kinetic equation… ▽ More We examine the applicability of the stochastic electron acceleration to two high synchrotron peaked blazars, Mrk 421 and Mrk 501, assuming synchrotron self-Compton emission of gamma-rays. Our model considers an emitting region moving at relativistic speed, where non-thermal electrons are accelerated and attain a steady-state energy spectrum together with the photons they emit. The kinetic equations of the electrons and photons are solved numerically, given a stationary wave number spectrum of the magnetohydrodynamic (MHD) disturbances, which are responsible for the electron acceleration and escape. Our simple formulation appears to reproduce the two well-sampled, long-term averaged photon spectra. In order to fit the model to the emission component from the radio to the X-ray bands, we need both a steeper wave spectral index than the Kolmogorov spectrum and efficient particle escape. Although the model provides a natural explanation for the high-energy cutoff of the electron energy distribution, the derived physical parameters raise a problem with an energy budget if the MHD waves with the Alfv{é}n velocity are assumed to be the acceleration agent. △ Less

Submitted 11 February, 2015; originally announced February 2015.

Comments: 9 pages, 4 figures, accepted for publication in MNRAS

Neutrino and Cosmic-Ray Emission and Cumulative Background from Radiatively Inefficient Accretion Flows in Low-Luminosity Active Galactic Nuclei

Authors: Shigeo S. Kimura, Kohta Murase, Kenji Toma

Abstract: We study high-energy neutrino and cosmic-ray (CR) emission from the cores of low-luminosity active galactic nuclei (LLAGN). In LLAGN, the thermalization of particles is expected to be incomplete in radiatively inefficient accretion flows (RIAFs), allowing the existence of non-thermal particles. In this work, assuming stochastic particle acceleration due to turbulence in RIAFs, we solve the Fokker-… ▽ More We study high-energy neutrino and cosmic-ray (CR) emission from the cores of low-luminosity active galactic nuclei (LLAGN). In LLAGN, the thermalization of particles is expected to be incomplete in radiatively inefficient accretion flows (RIAFs), allowing the existence of non-thermal particles. In this work, assuming stochastic particle acceleration due to turbulence in RIAFs, we solve the Fokker-Planck equation and calculate spectra of escaping neutrinos and CRs. The RIAF in LLAGN can emit CR protons with $\gtrsim10$ PeV energies and TeV-PeV neutrinos generated via $pp$ and/or $pγ$ reactions. We find that, if $\sim1$\% of the accretion luminosity is carried away by non-thermal ions, the diffuse neutrino intensity from the cores of LLAGN may be as high as $E_ν^2Φ_ν\sim3\times{10}^{-8} {\rm GeV} {\rm cm}^{-2} {\rm s}^{-1}{\rm sr}^{-1}$, which can be compatible with the observed IceCube data. This result does not contradict either of the diffuse gamma-ray background observed by {\it Fermi} or observed diffuse cosmic-ray flux. Our model suggests that, although very-high-energy gamma rays may not escape, radio-quiet AGN with RIAFs can emit GeV gamma-rays, which could be used for testing the model. We also calculate the neutron luminosity from RIAFs of LLAGN, and discuss a strong constraint on the model of jet mass loading mediated by neutrons from the diffuse neutrino observation. △ Less

Submitted 29 June, 2015; v1 submitted 13 November, 2014; originally announced November 2014.

Comments: 16 pages, 8 figures, Published in the Astrophysical Journal

Journal ref: Astrophys. J. 806 (2015) 159

Circular polarization in the optical afterglow of GRB 121024A

Authors: K. Wiersema, S. Covino, K. Toma, A. J. van der Horst, K. Varela, M. Min, J. Greiner, R. L. C. Starling, N. R. Tanvir, R. A. M. J. Wijers, S. Campana, P. A. Curran, Y. Fan, J. P. U. Fynbo, J. Gorosabel, A. Gomboc, D. Gotz, J. Hjorth, Z. P. Jin, S. Kobayashi, C. Kouveliotou, C. Mundell, P. T. O'Brien, E. Pian, A. Rowlinson , et al. (15 additional authors not shown)

Abstract: Gamma-ray bursts (GRBs) are most probably powered by collimated relativistic outflows (jets) from accreting black holes at cosmological distances. Bright afterglows are produced when the outflow collides with the ambient medium. Afterglow polarization directly probes the magnetic properties of the jet, when measured minutes after the burst, and the geometric properties of the jet and the ambient m… ▽ More Gamma-ray bursts (GRBs) are most probably powered by collimated relativistic outflows (jets) from accreting black holes at cosmological distances. Bright afterglows are produced when the outflow collides with the ambient medium. Afterglow polarization directly probes the magnetic properties of the jet, when measured minutes after the burst, and the geometric properties of the jet and the ambient medium when measured hours to days after the burst. High values of optical polarization detected minutes after burst in GRB 120308A indicate the presence of large-scale ordered magnetic fields originating from the central engine (the power source of the GRB). Theoretical models predict low degrees of linear polarization and negligable circular polarization at late times, when the energy in the original ejecta is quickly transferred to the ambient medium and propagates farther into the medium as a blastwave. Here we report the detection of circularly polarized optical light in the afterglow of GRB 121024A, measured 0.15 days after the burst. We show that the circular polarization is intrinsic to the afterglow and unlikely to be produced by dust scattering or plasma propagation effects. A possible explanation is to invoke anisotropic (rather than the commonly assumed isotropic) electron pitch angle distributions, and we suggest that new models are required to produce the complex microphysics of realistic shocks in relativistic jets. △ Less

Submitted 2 October, 2014; originally announced October 2014.

Comments: Published as Wiersema et al. 2014, Nature 509, 201-204. This is the version prior to final editing; please see official published version for the final version and higher quality images

Effects of High-Energy Particles on Accretion Flows onto a Supermassive Black Hole

Authors: Shigeo S. Kimura, Kenji Toma, Fumio Takahara

Abstract: We study effects of high-energy particles on the accretion flows onto a supermassive black hole and luminosities of escaping particles such as protons, neutrons, gamma-rays, and neutrinos. We formulate a one-dimensional model of the two-component accretion flow consisting of thermal particles and high-energy particles, supposing that some fraction of the released energy is converted to the acceler… ▽ More We study effects of high-energy particles on the accretion flows onto a supermassive black hole and luminosities of escaping particles such as protons, neutrons, gamma-rays, and neutrinos. We formulate a one-dimensional model of the two-component accretion flow consisting of thermal particles and high-energy particles, supposing that some fraction of the released energy is converted to the acceleration of the high-energy particles. The thermal component is governed by fluid dynamics while the high-energy particles obey the moment equations of the diffusion-convection equation. By solving the time evolution of these equations, we obtain advection dominated flows as the steady state solutions. Effects of the high-energy particles on the flow structures turn out to be small even if the pressure of the high-energy particles dominates over the thermal pressure. For a model in which the escaping protons take away almost all the released energy, the high-energy particles have large influence enough to make the flow have the Keplerian angular velocity at the inner region. We calculate the luminosities of the escaping particles for these steady solutions. The escaping particles can extract the energy from about $10^{-4}\dot M c^2$ to $10^{-2}\dot M c^2$, where $\dot M$ is the mass accretion rates. The luminosities of the escaping particles depend on the parameters such as the injection Lorentz factors, the mass accretion rates, and the diffusion coefficients. We also discuss some implications on the relativistic jet production by the escaping particles. △ Less

Submitted 1 July, 2014; originally announced July 2014.

Comments: 43 pages, 7 figures, 3 tables, accepted for publication in ApJ

Journal ref: Astrophys.J. 791 (2014) 100

Electromotive Force in the Blandford-Znajek Process

Authors: Kenji Toma, Fumio Takahara

Abstract: One of the mechanisms widely considered for driving relativistic jets in active galactic nuclei, galactic microquasars, and gamma-ray bursts is the electromagnetic extraction of the rotational energy of a central black hole, i.e., the Blandford-Znajek process, although the origin of the electromotive force in this process is still under debate. We study this process as the steady unipolar inductio… ▽ More One of the mechanisms widely considered for driving relativistic jets in active galactic nuclei, galactic microquasars, and gamma-ray bursts is the electromagnetic extraction of the rotational energy of a central black hole, i.e., the Blandford-Znajek process, although the origin of the electromotive force in this process is still under debate. We study this process as the steady unipolar induction in the Kerr black hole magnetosphere filled with a collisionless plasma screening the electric field (the D field) along the magnetic field (the B field), i.e., D dot B = 0. We extend the formulations and arguments made by Komissarov, and generally show that the origin of the electromotive force is ascribed to the ergosphere. It is explicitly shown that open magnetic field lines penetrating the ergosphere have a region where the D field is stronger than the B field in the ergosphere, and it keeps driving the poloidal currents and generating the electromotive force and the outward Poynting flux. The range of the possible value of the so-called angular velocity of the magnetic field line Omega_F is deduced for the field lines threading the equatorial plane in the ergosphere. We briefly discuss the relation between our conclusion and the ideal magnetohydrodynamic condition. △ Less

Submitted 28 May, 2014; originally announced May 2014.

Comments: 13 pages, 8 figures; Accepted for publication in MNRAS

Time-Dependent Models for Blazar Emissions with the Second-Order Fermi Acceleration

Authors: Katsuaki Asano, Fumio Takahara, Masaaki Kusunose, Kenji Toma, Jun Kakuwa

Abstract: The second-order Fermi acceleration (Fermi-II) driven by turbulence may be responsible for the electron acceleration in blazar jets. We test this model with time-dependent simulations. The hard electron spectrum predicted by the Fermi-II process agrees with the hard photon spectrum of 1ES 1101-232. For other blazars that show softer spectra, the Fermi-II model requires radial evolution of the elec… ▽ More The second-order Fermi acceleration (Fermi-II) driven by turbulence may be responsible for the electron acceleration in blazar jets. We test this model with time-dependent simulations. The hard electron spectrum predicted by the Fermi-II process agrees with the hard photon spectrum of 1ES 1101-232. For other blazars that show softer spectra, the Fermi-II model requires radial evolution of the electron injection rate and/or diffusion coefficient in the outflow. Such evolutions can yield a curved electron spectrum, which can reproduce the synchrotron spectrum of Mrk 421 from the radio to the X-ray regime.The photon spectrum in the GeV energy range of Mrk 421 is hard to fit with a synchrotron self-Compton model. However, if we introduce an external radio photon field with a luminosity of $4.9 \times 10^{38}~\mbox{erg}~\mbox{s}^{-1}$, GeV photons are successfully produced via inverse Compton scattering. The temporal variability of the diffusion coefficient or injection rate causes flare emission. The observed synchronicity of X-ray and TeV flares implies a decrease of the magnetic field in the flaring source region. △ Less

Submitted 22 November, 2013; v1 submitted 4 November, 2013; originally announced November 2013.

Comments: 34 pages, 15 figures, accepted for publication in the Astrophysical Journal

Polarization of GRB Prompt Emission

Authors: Kenji Toma

Abstract: We review the recent observational results of the gamma-ray linear polarization of Gamma-Ray Bursts (GRBs), and discuss some theoretical implications for the prompt emission mechanism and the magnetic composition of GRB jets. We also report a strict observational verification of CPT invariance in the photon sector as a result of the GRB polarization measurements. We review the recent observational results of the gamma-ray linear polarization of Gamma-Ray Bursts (GRBs), and discuss some theoretical implications for the prompt emission mechanism and the magnetic composition of GRB jets. We also report a strict observational verification of CPT invariance in the photon sector as a result of the GRB polarization measurements. △ Less

Submitted 26 August, 2013; originally announced August 2013.

Comments: 13 pages, 5 figures, 7th Huntsville Gamma-Ray Burst Symposium, GRB 2013: paper 45 in eConf Proceedings C1304143