Search

The afterglows of gamma-ray bursts are non-thermal electron synchrotron emissions from relativistic shocks. The origin of strong magnetic field in the emission region remains elusive, and two field amplification mechanisms via the plasma kinetic and magnetohydrodynamic instabilities have been discussed. The polarimetric observations are a powerful probe to distinguish these two mechanisms. So far, most theoretical works have focused on the former mechanism and constructed afterglow polarization models with microscopic-scale turbulence whose coherence length is much smaller than the thickness of the blast wave. In this work, focusing on the latter mechanism, we utilize our semi-analytic model of the synchrotron polarization with large-scale turbulence whose coherence length is comparable to the thickness of the blast wave to investigate the effect of magnetic field anisotropy and the observer viewing angle. We find that the polarization in our large-scale turbulence model can exhibit both behaviors characteristic of the microscopic-scale turbulence model and those not seen in the microscopic-scale model. Then we find that the large-scale model could explain all the polarimetric observational data to date that seem to be forward shock emission. We also examine the effect of ordered-field component, and find that polarization degree and polarization angle constant in time are realized only when the energy density ratio of the ordered and fluctuated components is $\gtrsim 50$. In this case, however, the polarization degree is much higher than the observed values., Comment: Submitted to ApJ, comments welcome

The Fermi and eROSITA bubbles, large diffuse structures in our Galaxy, can be the by-products of the steady star formation activity. To simultaneously explain the star formation history of the Milky Way and the metallicity of $\sim$ Z$_\odot$ at the Galactic disk, a steady Galactic wind driven by cosmic-rays is required. For tenuous gases with a density of $\lesssim$10$^{-3}$ cm$^{-3}$, the cosmic-ray heating dominates over radiative cooling, and the gas can maintain the virial temperature of $\sim$0.3 keV ideal for escape from the Galactic system as the wind. A part of the wind falls back onto the disk like a galactic fountain flow. We model the wind dynamics according to the Galactic evolution scenario and find that the scale height and surface brightness of the X-ray and the hadronic gamma-ray emissions from such fountain flow region can be consistent with the observed properties of the Fermi and eROSITA bubbles. This implies that the bubbles are persistent structures of the Milky Way existing over (at least) the last $\sim$1 Gyr, rather than evanescent structures formed by non-trivial, $\sim$10 Myr past Galactic Center transient activities., Comment: 13 pages, 9 figures, accepted for publication in ApJ

Magnetars, which are neutron stars with strong magnetic fields, exhibit occasional bursting activities. The shape of a magnetar is not perfectly spherical due to the Lorentz force exerted by its strong magnetic fields and is described as a triaxial body. We study the unstable free precession in a triaxial magnetar; one of the principal axes undergoes an upside-down flip. This flip is known as the Dzhanibekov effect. We find that during the flip, the Euler force can suddenly disturb the force balance on the surface layer of the magnetar, potentially leading to plastic flow of the layer. This, in turn, may trigger different forms of magnetar activity, such as the emission of the bursts and/or of gravitational waves., Comment: 8 pages, 4 figures, accepted for publication in ApJ

We study the long-term evolution of the Milky Way (MW) over cosmic time by modeling the star formation, cosmic rays, metallicity, stellar dynamics, outflows and inflows of the galactic system to obtain various insights into the galactic evolution. The mass accretion is modeled by the results of cosmological N-body simulations for the cold dark matter. We find that the star formation rate is about half the mass accretion rate of the disk, given the consistency between observed Galactic Diffuse X-ray Emissions (GDXEs) and possible conditions driving the Galactic wind. Our model simultaneously reproduces the quantities of star formation rate, cosmic rays, metals, and the rotation curve of the current MW. The most important predictions of the model are that there is an unidentified accretion flow with a possible number density of $\sim10^{-2}$ cm$^{-3}$ and the part of the GDXEs originates from a hot, diffuse plasma which is formed by consuming about 10 % of supernova explosion energy. The latter is the science case for future X-ray missions; XRISM, Athena, and so on. We also discuss further implications of our results for the planet formation and observations of externalgalaxies in terms of the multimessenger astronomy., Comment: 18 pages, 11 figures, accepted for publication in PASJ

We investigate the role of magnetic field on the gas dynamics in a galactic bulge region by three dimensional simulations with radiative cooling and heating. While high-temperature corona with $T>10^6\ {\rm K}$ is formed in the halo regions, the temperature near the mid-plane is $\lesssim 10^4\ {\rm K}$ following the thermal equilibrium curve determined by the radiative cooling and heating. Although the thermal energy of the interstellar gas is lost by radiative cooling, the saturation level of the magnetic field strength does not significantly depend on the radiative cooling and heating. The magnetic field strength is amplified to $10\ {\rm \mu G}$ on average, and reaches several hundred ${\rm \mu G}$ locally. We find the formation of magnetically dominated regions at mid-latitudes in the case with the radiative cooling and heating, which is not seen in the case without radiative effect. The vertical thickness of the mid-latitude regions is $50-150\ {\rm pc}$ at the radial location of $0.4-0.8\ {\rm kpc}$ from the galactic center, which is comparable to the observed vertical distribution of neutral atomic gas. When we take the average of different components of energy density integrated over the galactic bulge region, the magnetic energy is comparable to the thermal energy. We conclude that the magnetic field plays a substantial role in controlling the dynamical and thermal properties of the galactic bulge region., Comment: Accepted to ApJ; 24 pages, 20 figures 3 tables

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., Comment: Author's version. 41 pages, 9 figures, 4 tables

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.

A novel collisionless shock jump condition is suggested by modeling the entropy production at the shock transition region. We also calculate downstream developments of the atomic ionization balance and the ion temperature relaxation in supernova remnants (SNRs). The injection process and subsequent acceleration of cosmic-rays (CRs) in the SNR shocks are closely related to the formation process of the collisionless shocks. The formation of the shock is caused by wave-particle interactions. Since the wave-particle interactions result in energy exchanges between electromagnetic fields and charged particles, the randomization of particles associated with the shock transition may occur with the rate given by the scalar product of the electric field and current. We find that order-of-magnitude estimates of the randomization with reasonable strength of the electromagnetic fields in the SNR constrain the amount of the CR nuclei and ion temperatures. The constrained amount of the CR nuclei can be sufficient to explain the Galactic CRs. The ion temperature becomes significantly lower than in the case of no CRs. To distinguish the case without CRs, we perform synthetic observations of atomic line emissions from the downstream region of the SNR RCW~86. Future observations by {\it XRISM} and {\it Athena} can distinguish whether the SNR shock accelerates the CRs or not from the ion temperatures., Comment: 18 pages, 14 figures, 3 tables, accepted for publication in PASJ

The effects of cosmic-ray diffusion and radiative cooling on the structure of the Galactic wind are studied using a steady state approximation. It is known that realistic cooling processes suppress the wind from launching. The effects of cosmic-ray diffusion are also supposed to be unfavorable for launching the wind. Both of these effects have not been studied simultaneously in a steady-state approximation of the wind. We find 327,254 solutions of the steady-state Galactic wind and confirm that: the effect of cosmic-ray pressure depends on the Alfv{\'e}n Mach number, the mass flux carried by the wind does not depend on the cosmic-ray pressure directly (but depends on the thermal pressure), and the typical conditions found in the Galaxy may correspond to the wind solution that provides metal polluted matters at a height of $\sim300$~kpc from the disk., Comment: 18 pages, 13 figures, accepted for publication in ApJ

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., Comment: 16 pages, 12 figures, revised version of 2, accepted by ApJ in publication

A novel diagnostic of cosmic-ray modified shocks by polarimetry of H $\alpha$ emissions is suggested. In a cosmic-ray modified shock, the pressure of cosmic rays is sufficiently high compared to the upstream ram pressure to force the background plasma to decelerate (measured in the shock rest frame). Simultaneously, a fraction of the hydrogen atoms co-existing in the upstream plasma collide with the decelerated protons and undergo charge-exchange reactions. As a result, hydrogen atoms with the same bulk velocity of the decelerated protons are generated. We show that when the shock is observed from edge-on, the H $\alpha$ radiated by these upstream hydrogen atoms is linearly polarized with a sizable degree of a few per cent as a result of resonant scattering of Ly $\beta$. The polarization direction depends strongly on the velocity modification; the direction is parallel to the shock surface for the case of no modification, while the direction is parallel to the shock velocity for the case of a modified shock., Comment: 9 pages, 8 figures, accepted for publication in MNRAS

Radiative transfer in hydrogen lines in supernova remnant (SNR) shock waves is studied taking into account the population of the hydrogen atom 2s-state. Measurements of Balmer line emission, especially of H~$\alpha$, are often relied upon to derive physical conditions in the SNR shock. On the other hand, Lyman series photons, especially Ly~$\beta$, are mostly absorbed by upstream hydrogen atoms. As a result, atoms are excited to the 3p state, and then emit H~$\alpha$ by the spontaneous transition from 3p to 2s. Thus, the nature of H~$\alpha$ depends on how many Ly~$\beta$ photons are converted to H~$\alpha$ photons. Moreover, the Balmer lines can be scattered by the 2s-state hydrogen atoms, which are excited not only by collisional excitation but also by the Lyman-Balmer conversion. It is shown for example that the H~$\alpha$ photons are scattered if the shock propagates into an H~$_{\rm I}$ cloud with a density of $\sim30~{\rm cm^{-3}}$ and a size of $\sim 1$~pc. We find that the line profile of H~$\alpha$ becomes asymmetric resulting from the difference between line centre frequencies among the transitions from 3s to 2p, from 3p to 2s and from 3d to 2p. We also find that the broad-to-narrow ratio of H~$\alpha$, which is often used to estimate the ion-electron temperature equilibrium, varies at most $\simeq 10$ per cent depending on the ionization degree of the upstream medium because of incomplete conversion of Lyman lines to Balmer lines., Comment: 16 pages, accepted for publication in MNRAS

The spectral slope of the magnetic energy in supernova remnants (SNRs) can be obtained by analysis of spatial two-point correlation functions of synchrotron intensities. This method has been originally developed for the analysis of magnetic field structure in diffuse interstellar medium and applied when the geometry of the emission region is simple and known. In this paper, applying this correlation analysis to Tycho's SNR, for which the synchrotron emission region is known to be a spherical shell, we find that the magnetic energy spectrum shows the Kolmogorov-like scaling. Our results can be explained by turbulence developed downstream of the shock front via Richtmyer-Meshkov instability or an amplification of upstream magnetic field induced by cosmic rays. They could be discriminated by future observations with a sub arcsecond resolution such as Square Kilometer Array., Comment: 6 pages, 3 figures, accepted by MNRAS

Linearly polarized Balmer line emissions from supernova remnant shocks are studied taking into account the energy loss of the shock owing to the production of nonthermal particles. The polarization degree depends on the downstream temperature and the velocity difference between upstream and downstream regions. The former is derived once the line width of the broad component of the H$\alpha$ emission is observed. Then, the observation of the polarization degree tells us the latter. At the same time, the estimated value of the velocity difference independently predicts adiabatic downstream temperature that is derived from Rankine-Hugoniot relations for adiabatic shocks. If the actually observed downstream temperature is lower than the adiabatic temperature, there is a missing thermal energy which is consumed for particle acceleration. It is shown that a larger energy loss rate leads to more highly polarized H$\alpha$ emission. Furthermore, we find that polarized intensity ratio of H$\beta$ to H$\alpha$ also depends on the energy loss rate and that it is independent of uncertain quantities such as electron temperature, the effect of Lyman line trapping and our line of sight., Comment: 15 pages, accepted by MNRAS

We present high-resolution long-slit spectroscopy of a Balmer-dominated shock in the northeastern limb of the Cygnus Loop with the Subaru high dispersion spectrograph. By setting the slit angle along the shock normal, we investigate variations of the flux and profile of the H-alpha line from preshock to postshock regions with a spatial resolution of about 4 times 10^{15} cm. The H-alpha line profile can be represented by a narrow (28.9+/-0.7 km/s) Gaussian in a diffuse region ahead of the shock, i.e., a photoionization precursor, and narrow (33.1+/-0.2 km/s) plus broad (130-230 km/s) Gaussians at the shock itself. We find that the width of the narrow component abruptly increases up to 33.1+/-0.2 km/s, or 38.8+/-0.4 km/s if we eliminate projected emission originating from the photoionization precursor, in an unresolved thin layer (< 4 times 10^{15} cm at a distance of 540 pc) at the shock. We show that the sudden broadening can be best explained by heating via damping of Alfven waves in a thin cosmic-ray precursor, although other possibilities are not fully ruled out. The thickness of the cosmic-ray precursor in the Cygnus Loop (a soft gamma-ray emitter) is an order of magnitude thinner than that in Tycho's Knot g (a hard gamma-ray emitter), which may be caused by different energy distribution of accelerated particles between the two sources. In this context, systematic studies might reveal a positive correlation between the thickness of the cosmic-ray precursor and the hardness of the cosmic-ray energy distribution., Comment: 16 pages, 5 figures, accepted for publication in The Astrophysical Journal Letters

Using three-dimensional magnetohydrodynamics simulations, we show that the efficiency of cosmic-ray (CR) production at supernova remnants (SNRs) is over-predicted if it could be estimated based on proper motion measurements of H$\alpha$ filaments in combination with shock-jump conditions. Density fluctuations of upstream medium make shock waves rippled and oblique almost everywhere. The kinetic energy of the shock wave is transferred into that of downstream turbulence as well as thermal energy which is related to the shock velocity component normal to the shock surface. Our synthetic observation shows that the CR acceleration efficiency as estimated from a lower downstream plasma temperature, is overestimated by 10-40%, because rippled shock does not immediately dissipate all upstream kinetic energy., Comment: 7 pages, 3 figures, 1 table, accepted for publication in ApJ; the paper with full resolution images is http://www.phys.aoyama.ac.jp/~ryo/papers/shimoda2015.pdf

It has been suggested by radio observations of polarized synchrotron emissions that downstream magnetic field in some young supernova remnants are oriented radially. We study magnetic field distribution of turbulent supernova remnant driven by the Richtmyer-Meshkov instability -- in other words, the effect of rippled shock -- by using three-dimensional magnetohydrodynamics simulations. We find that the induced turbulence has radially biased anisotropic velocity dispersion that leads to a selective amplification of the radial component of the magnetic field. The Richtmyer-Meshkov instability is induced by the interaction between the shock and upstream density fluctuations. Future high-resolution polarization observations can distinguish the following candidates responsible for the upstream density fluctuations: (i) inhomogeneity caused by the cascade of large-scale turbulence in the ISM so-called the big-power-law-in-the-sky, (ii) structures generated by the Drury instability in the cosmic-ray modified shock, and (iii) fluctuations induced by the non-linear feedback of the cosmic-ray streaming instability., Comment: 5 pages, 3 figures, accepted by ApJ Letters