Researcher profile

Kunihito Ioka

Kunihito Ioka contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
15works
0followers
7topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

15 published item(s)

preprint2026arXiv

One-dimensional PIC Simulation of Induced Compton Scattering in Magnetized Electron-Positron Pair Plasma

We investigate induced Compton scattering of a circularly polarized Alfvén wave propagating in a magnetized electron-positron pair plasma using one-dimensional Particle-in-Cell (PIC) simulations. In this system, two distinct modes of density fluctuations, referred to as the charged mode and the neutral mode, are theoretically expected to arise through parametric instabilities. Our simulations confirm these predictions: in the charged mode, the electron and positron densities fluctuate oppositely (Langmuir-like), while in the neutral mode, the charge is Debye-screened and both species fluctuate in phase (acoustic-like). The linear growth rates obtained from the simulations are in good agreement with analytical estimates for both modes. We also find that, in some cases, the linear growth saturates before full scattering occurs, allowing the incident wave to propagate without significant attenuation. Our results allow us to determine whether induced Compton scattering grows linearly in magnetized pair plasmas, offering a foundation for studies of fast radio bursts and laser-plasma experiments.

preprint2023arXiv

Expanding Fireball in Magnetar Bursts and Fast Radio Bursts

A fireball of radiation plasma created near the surface of a neutron star (NS) expands under its own pressure along magnetic field lines, and produces photon emission and relativistic matter outflow. We comprehensively classify the expanding fireball evolution into five cases and obtain the photospheric luminosity and the kinetic energy of the outflow, taking into account key processes; lateral diffusion of photons escaping from a magnetic flux tube, effects of strong magnetic field, baryon loading from the NS surface, and radiative acceleration via cyclotron resonant scattering, some of which have not been considered in the context of gamma-ray bursts. Applying our model to magnetar bursts with fast radio bursts (FRBs), in particular the X-ray short bursts from SGR 1935+2154 associated with the Galactic FRB 20200428A, we show that the burst radiation can accelerate the outflow to high Lorentz factor with sufficient energy to power FRBs.

preprint2022arXiv

Cocoon breakout and escape from the ejecta of neutron star mergers

The cocoon is an inevitable product of a jet propagating through ambient matter, and takes a fair fraction of the jet energy. In short gamma-ray bursts, the ambient matter is the ejecta from the merger of neutron stars, expanding with a high velocity $\sim 0.2 c$, in contrast to the static stellar envelope in collapsars. Using 2D relativistic hydrodynamic simulations with the ejecta density profile as $ρ\propto r^{-2}$, we find that the expansion makes a big difference; only 0.5--5\% of the cocoon mass escapes from (faster than) the ejecta, with an opening angle $20^{\circ}$--$30^{\circ}$, while it is $\sim 100\%$ and spherical in collapsars. We also analytically obtain the shares of mass and energies for the escaped and trapped cocoons. Because the mass of the escaped cocoon is small and the trapped cocoon is concealed by the ejecta and the escaped cocoon, we conclude that it is unlikely that the cocoon emission was observed as a counterpart to the gravitational wave event GW170817.

preprint2022arXiv

Probing cosmology and gastrophysics with fast radio bursts: Cross-correlations of dark matter haloes and cosmic dispersion measures

For future surveys of fast radio bursts (FRBs), we clarify information available from cosmic dispersion measures (DMs) through cross-correlation analyses of foreground dark matter haloes (hosting galaxies and galaxy clusters) with their known redshifts. With a halo-model approach, we predict that the cross-correlation with cluster-sized haloes is less affected by the details of gastrophysics, providing robust cosmological information. For less massive haloes, the cross-correlation at angular scales of $<10\, \mathrm{arcmin}$ is sensitive to gas expelled from the halo centre due to galactic feedback. Assuming $20000$ FRBs over $20000 \, {\rm deg}^2$ with a localisation error being 3 arcmin, we expect that the cross-correlation signal at halo masses of $10^{12}-10^{14}\, M_\odot$ can be measured with a level of $\sim 1\%$ precision in a redshift range of $0<z<1$. Such precise measurements enable one to put a 1.5\% level constraint on $σ_8\, (Ω_\mathrm{M}/0.3)^{0.5}$ and a 3\% level constraint on $(Ω_\mathrm{b}/0.049)(h/0.67)(f_\mathrm{e}/0.95)$ ($σ_8$, $Ω_\mathrm{M}$, $Ω_\mathrm{b}$, $h$ and $f_\mathrm{e}$ are the linear mass variance smoothed at $8\, h^{-1}\mathrm{Mpc}$, mean mass density, mean baryon density, the present-day Hubble parameter and fraction of free electrons in cosmic baryons today), whereas the gas-to-halo mass relation in galaxies and clusters can be constrained with a level of $10\%-20\%$. Furthermore the cross-correlation analyses can break the degeneracy among $Ω_\mathrm{b}$, $h$ and $f_\mathrm{e}$, inherent in the DM-redshift relation. Our proposal opens new possibilities for FRB cosmology, while it requires extensive galaxy redshift catalogues and further improvement of the halo model.

preprint2021arXiv

Binary comb models for FRB 121102

The first repeating fast radio burst source, FRB 121102, is observed to emit bursts periodically. We show that FRB 121102 can be interpreted as an interacting neutron star binary system with an orbital period of 159 days. We develop a binary comb model by introducing an eccentricity in the orbit. Besides the original funnel mode of the binary comb model, which was applied to FRB 180916.J0158+65 by Ioka and Zhang 2020, we also identify two new modes of the binary comb model, i.e. the tau-crossing mode and the inverse funnel mode, and apply them to interpret FRB 121102. These new developments expand the applicable parameter space, allowing the companion star to be a massive star, a massive black hole, or a supermassive black hole, with the latter two having larger parameter spaces. These models are also consistent with other observations, such as the persistent bright radio counterpart associated with the source. We also argue that the observed frequency dependence of the active window does not disfavor the binary comb model, in contrast to recent claims, and propose two possible scenarios to interpret the data.

preprint2021arXiv

Diverse Jet Structures Consistent with the Off-axis Afterglow of GRB 170817A

The jet structure of short gamma-ray bursts (GRBs) has been controversial after the detection of GRB 170817A as the electromagnetic counterparts to the gravitational wave event GW170817. Different authors use different jet structures for calculating the afterglow light curves. We formulated a method to inversely reconstruct a jet structure from a given off-axis GRB afterglow, without assuming any functional form of the structure. By systematically applying our inversion method, we find that more diverse jet structures are consistent with the observed afterglow of GRB 170817A within errors: such as hollow-cone, spindle, Gaussian, and power-law jet structures. In addition, the total energy of the reconstructed jet is arbitrary, proportional to the ambient density $n_0$, with keeping the same jet shape if the parameters satisfy the degeneracy combination $n_0 \varepsilon_\mathrm{B}^{(p+1)/(p+5)} \varepsilon_\mathrm{e}^{4(p-1)/(p+5)} = \mathrm{const.}$. Observational accuracy less than $\sim 6$ per cent is necessary to distinguish the different shapes, while the degeneracy of the energy scaling would be broken by observing the spectral breaks and viewing angle. Future events in denser environment with brighter afterglows and observable spectral breaks are ideal for our inversion method to pin down the jet structure, providing the key to the jet formation and propagation.

preprint2021arXiv

Statistical modelling of the cosmological dispersion measure

We have investigated the basic statistics of the cosmological dispersion measure (DM) -- such as its mean, variance, probability distribution, angular power spectrum and correlation function -- using the state-of-the-art hydrodynamic simulations, IllustrisTNG300, for the fast radio burst (FRB) cosmology. To model the DM statistics, we first measured the free-electron abundance and the power spectrum of its spatial fluctuations. The free-electron power spectrum turns out to be consistent with the dark matter power spectrum at large scales, but it is strongly damped at small scales ($\lesssim 1$Mpc) owing to the stellar and active galactic nucleus feedback. The free-electron power spectrum is well modelled using a scale-dependent bias factor (the ratio of its fluctuation amplitude to that of the dark matter). We provide analytical fitting functions for the free-electron abundance and its bias factor. We next constructed mock sky maps of the DM by performing standard ray-tracing simulations with the TNG300 data. The DM statistics are calculated analytically from the fitting functions of the free-electron distribution, which agree well with the simulation results measured from the mock maps. We have also obtained the probability distribution of source redshift for a given DM, which helps in identifying the host galaxies of FRBs from the measured DMs. The angular two-point correlation function of the DM is described by a simple power law, $ξ(θ) \approx 2400 (θ/{\rm deg})^{-1} \, {\rm pc}^2 \, {\rm cm}^{-6}$, which we anticipate will be confirmed by future observations when thousands of FRBs are available.

preprint2020arXiv

A Binary Comb Model for Periodic Fast Radio Bursts

We show that the periodic FRB 180916.J0158+65 can be interpreted by invoking an interacting neutron star binary system with an orbital period of $\sim 16$ days. The FRBs are produced by a highly magnetized pulsar, whose magnetic field is ``combed&#39;&#39; by the strong wind from a companion star, either a massive star or a millisecond pulsar. The FRB pulsar wind retains a clear funnel in the companion&#39;s wind that is otherwise opaque to induced Compton or Raman scatterings for repeating FRB emission. The 4 day active window corresponds to the time when the funnel points toward Earth. The interaction also perturbs the magnetosphere of the FRB pulsar and may trigger emission of FRBs. We derive the physical constraints on the comb and the FRB pulsar from the observations and estimate the event rate of FRBs. In this scenario, a lower limit on the period of observable FRBs is predicted. We speculate that both the intrinsic factors (strong magnetic field and young age) and the extrinsic factor (interaction) may be needed to generate FRBs in neutron star binary systems.

preprint2020arXiv

Analytic properties of electromagnetic field of binary compact stars and electromagnetic precursors to gravitational waves

We analytically study the properties of the electromagnetic field in vacuum around close binary compact stars containing at least one neutron star. We show that the orbital motion of the neutron star induces high multipole modes of the electromagnetic field just before the merger. These modes are superimposed to form a spiral arm configuration and its edge is found to be a likely site for magnetic reconnection. These modes also enhance the total Poynting flux from neutron star binaries by a factor of 2--4. We also indicate that the electric field induced by the orbital motion lead to a magnetosphere around binaries and estimate its plasma density, which has a different parameter dependence than Goldreich-Julian density. With these properties, we discuss possible electromagnetic counterparts to gravitational wave events, and identify radio precursors, such as fast radio bursts, as the most promising observational targets.

preprint2020arXiv

Current status of space gravitational wave antenna DECIGO and B-DECIGO

Deci-hertz Interferometer Gravitational Wave Observatory (DECIGO) is the future Japanese space mission with a frequency band of 0.1 Hz to 10 Hz. DECIGO aims at the detection of primordial gravitational waves, which could be produced during the inflationary period right after the birth of the universe. There are many other scientific objectives of DECIGO, including the direct measurement of the acceleration of the expansion of the universe, and reliable and accurate predictions of the timing and locations of neutron star/black hole binary coalescences. DECIGO consists of four clusters of observatories placed in the heliocentric orbit. Each cluster consists of three spacecraft, which form three Fabry-Perot Michelson interferometers with an arm length of 1,000 km. Three clusters of DECIGO will be placed far from each other, and the fourth cluster will be placed in the same position as one of the three clusters to obtain the correlation signals for the detection of the primordial gravitational waves. We plan to launch B-DECIGO, which is a scientific pathfinder of DECIGO, before DECIGO in the 2030s to demonstrate the technologies required for DECIGO, as well as to obtain fruitful scientific results to further expand the multi-messenger astronomy.

preprint2020arXiv

Inverse reconstruction of jet structure from off-axis gamma-ray burst afterglows

The gravitational wave event GW170817 and the slowly-rising afterglows of short gamma-ray burst GRB 170817A clearly suggest that the GRB jet has an angular structure. However the actual jet structure remains unclear as different authors give different structures. We formulate a novel method to inversely reconstruct the jet structure from off-axis GRB afterglows, without assuming any functional form of the structure in contrast to the previous studies. The jet structure is uniquely determined from the rising part of a light curve for a given parameter set by integrating an ordinary differential equation, which is derived from the standard theory of GRB afterglows. Applying to GRB 170817A, we discover that a non-trivial hollow-cone jet is consistent with the observed afterglows, as well as Gaussian and power-law jets within errors, which implies the Blandford-Znajek mechanism or an ejecta-jet interaction. The current observations only constrain the jet core, not in principle the outer jet structure around the line of sight. More precise and high-cadence observations with our inversion method will fix the jet structure, providing a clue to the jet formation and propagation.

preprint2019arXiv

Jet Propagation in Neutron Star Mergers and GW170817

The gravitational wave event from the binary neutron star (BNS) merger GW170817 and the following multi-messenger observations present strong evidence for i) merger ejecta expanding with substantial velocities and ii) a relativistic jet which had to propagate through the merger ejecta. The ejecta&#39;s expansion velocity is not negligible for the jet head motion, which is a fundamental difference from the other systems like collapsars and active galactic nuclei. Here we present an analytic model of the jet propagation in an expanding medium. In particular, we notice a new term in the expression of the breakout time and velocity. In parallel, we perform a series of over a hundred 2D numerical simulations of jet propagation. The BNS merger ejecta is prepared based on numerical relativity simulations of a BNS merger with the highest-resolution to date. We show that our analytic results agree with numerical simulations over a wide parameter space. Then we apply our analytic model to GW170817, and obtain two solid constraints on: i) the central engine luminosity as $L_{iso,0} \sim 3\times10^{49}-2.5\times10^{52}$ erg s$^{-1}$, and on ii) the delay time between the merger and engine activation $t_0-t_m < 1.3$ s. The engine power implies that the apparently-faint \textit{short} gamma-ray burst (\textit{s}GRB) \textit{s}GRB 170817A is similar to typical \textit{s}GRBs if observed on-axis.

preprint2007arXiv

Low-Luminosity GRB 060218: A Collapsar Jet from a Neutron Star, Leaving a Magnetar as a Remnant?

The gamma-ray burst (GRB) 060218 has ~ 10^5 times lower luminosity than typical long GRBs, and is associated with a supernova (SN). The radio afterglow displays no jet break, so that this burst might arise from a mildly-relativistic spherical outflow produced by the SN shock sweeping the stellar surface. Since this model is energetically difficult, we propose that the radio afterglow is produced by a non-relativistic phase of an initially collimated outflow (jet). Our jet model is supported by the detection of optical linear polarization in the SN component. We also show analytically that the jet can penetrate a progenitor star. We analyzed the observational data of the prompt emission of this burst and obtained a smooth power-law light curve which might last longer than 10^6s. This behavior contrasts with the long intermittent activities with the X-ray flares of typical GRBs, implying that the central engine of this burst is different from those of typical GRBs. This argument is consistent with the analysis of the SN component of this burst, which suggests that the progenitor star was less massive and collapsed to a neutron star instead of a black hole. The collimation-corrected event rate of such low-luminosity GRBs is estimated to be ~ 10 times higher than that of typical long GRBs, and they might form a different GRB population: low-luminosity GRBs are produced by mildly-relativistic jets from neutron stars at the collapses of massive stars, while typical long GRBs by highly-relativistic jets from black holes. We suggest that the central engine of GRB 060218 is a pulsar (or a magnetar) with the initial rotation period P_0 \~ 10ms and the magnetic field B ~ 10^16G. A giant flare from the magnetar might be observed in future.

preprint2006arXiv

High Energy Neutrinos and Cosmic-Rays from Low-Luminosity Gamma-Ray Bursts?

The recently discovered gamma-ray burst (GRB) GRB 060218/SN 2006aj is classified as an X-ray Flash with very long duration driven possibly by a neutron star. Since GRB 060218 is very near about 140 Mpc and very dim, one-year observation by Swift suggests that the rate of GRB 060218-like events might be very high so that such low luminosity GRBs (LL-GRBs) might form a different population from the cosmological high luminosity GRBs (HL-GRBs). We found that the high energy neutrino background from LL-GRBs could be comparable with that from HL-GRBs. If each neutrino event is detected by IceCube, later optical-infrared follow-up observations such as by Subaru and HST have possibilities to identify a Type Ibc supernova associated with LL-GRBs, even if gamma- and X-rays are not observed by Swift. This is in a sense a new window from neutrino astronomy, which might enable us to confirm the existence of LL-GRBs and to obtain information about their rate and origin. We also argue LL-GRBs as high energy gamma-ray and cosmic-ray sources.

preprint2004arXiv

A Unified Model of Short and Long Gamma-Ray Bursts, X-Ray Rich Gamma-Ray Bursts, and X-Ray Flashes

We propose a possible unified model of short and long gamma-ray bursts (GRBs), X-ray rich GRBs, and X-ray flashes. The jet of a GRB is assumed to consist of multiple sub-jets or sub-shells (i.e., an inhomogeneous jet model). The multiplicity of the sub-jets along a line of sight n_s is an important parameter. If n_s is large (>> 1) the event looks like a long GRB, while if n_s=1, the event looks like a short GRB. Finally, when n_s=0, the event looks like an X-ray flash or an X-ray rich GRB. Furthermore, our model may also explain the bimodal distributions of T_{90} duration of BATSE-GRBs.