BZPEERReading, trust and collaboration for research
Menu

Discover

GraphFollow connections around a paper, topic or saved view.

Workspaces

HomePick up where your research flow left off.InboxHandle requests, replies and notifications from one place.CollectionsCurate reading lists, evidence sets and shareable dossiers.ResearchSee your private provenance graph, trail and imports.CollaborationMove from discovery into focused discussion rooms.PublishDraft, preview and publish full-text research articles.

Network

ExpertsFind specialists worth following or asking for help.CommunitiesJoin research circles organized around shared work.PartnersSee external organizations active around the same topics.

Opportunities

ListingsBrowse roles, calls and collaborations with clearer fit signals.

Account

Log inReturn to your reading and collaboration workspace.Create accountStart saving work, following people and joining teams.
Log inCreate account

Researcher profile

Ken-Ichi Nishikawa

Ken-Ichi Nishikawa contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate
astro-ph.HEgr-qcphysics.plasm-phphysics.comp-phastro-ph.SRphysics.optics

Trust snapshot

Quick read

Trust 21 - Emerging
14works
0followers
6topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

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

14 published item(s)

preprint2023arXiv

Femtosecond laser-induced sub-wavelength plasma inside dielectrics: III. Terahertz radiation emission

Electromagnetic radiation within the terahertz (THz) frequency range is of great interest for applications in remote sensing and time-domain spectroscopy. The laser-induced plasmas are promising mediums for generating THz radiation. It has been recently reported that focusing femtosecond Bessel pulses inside dielectrics induces a high aspect ratio over-critical plasmas. Here we show that the intense resonantly driven electrostatic fields at the so-called critical surface lead to THz radiation emission. Through three-dimensional particle-in-cell simulation and analytical derivation, we have investigated the emission of THz radiation. We show that the THz radiation is associated with a hot population of electrons trapped in ambipolar electric fields of the double layers.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Two-dimensional Particle-in-Cell simulations of axisymmetric black hole magnetospheres

We investigate the temporal evolution of an axisymmetric magnetosphere around a rapidly rotating, stellar-mass black hole, applying a two-dimensional particle-in-cell simulation scheme. Adopting a homogeneous pair production, and assuming that the mass accretion rate is much less than the Eddington limit, we find that the black hole's rotational energy is preferentially extracted from the middle latitudes, and that this outward energy flux exhibits an enhancement that lasts approximately 160 dynamical time scales. It is demonstrated that the magnetohydrodynamic approximations cannot be justified in such a magnetically-dominated magnetosphere, because the Ohm's law completely breaks down, and because the charge-separated electron-positron plasmas are highly non-neutral. An implication is given regarding the collimation of relativistic jets.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Collisionless electron-ion shocks in relativistic unmagnetized jet-ambient interactions: Non-thermal electron injection by double layer

The course of non-thermal electron ejection in relativistic unmagnetized electron-ion shocks is investigated by performing self-consistent particle-in-cell simulations. The shocks are excited through the injection of relativistic jet into ambient plasma, leading to two distinct shocks (named as the trailing shock and leading shock) and a contact discontinuity. The Weibel-like instabilities heat the electrons up to approximately half of ion kinetic energy. The double layers formed in the trailing and leading edges then accelerated the electrons by the ion kinetic energy. The electron distribution function in the leading edge shows a clear non-thermal power-law tail which contains $\sim1\%$ of electrons and $\sim8\%$ of electron energy. Its power-law index is -2.6. The acceleration efficiency is $\sim23\%$ by number and $\sim50\%$ by energy and the power-law index is -1.8 for electron distribution function in the trailing edge. The effect of the dimensionality is examined by comparing results of 3D simulation with 2D ones. It exhibits that the electron acceleration is more efficient in 2D.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Microscopic Processes in Global Relativistic Jets Containing Helical Magnetic Fields

In the study of relativistic jets one of the key open questions is their interaction with the environment on the microscopic level. Here, we study the initial evolution of both electron$-$proton ($e^{-}-p^{+}$) and electron$-$positron ($e^{\pm}$) relativistic jets containing helical magnetic fields, focusing on their interaction with an ambient plasma. We have performed simulations of "global" jets containing helical magnetic fields in order to examine how helical magnetic fields affect kinetic instabilities such as the Weibel instability, the kinetic Kelvin-Helmholtz instability (kKHI) and the Mushroom instability (MI). In our initial simulation study these kinetic instabilities are suppressed and new types of instabilities can grow. In the $e^{-}-p^{+}$ jet simulation a recollimation-like instability occurs and jet electrons are strongly perturbed. In the $e^{\pm}$ jet simulation a recollimation-like instability occurs at early times followed by a kinetic instability and the general structure is similar to a simulation without helical magnetic field. Simulations using much larger systems are required in order to thoroughly follow the evolution of global jets containing helical magnetic fields.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Non-relativistic perpendicular shocks modeling young supernova remnants: nonstationary dynamics and particle acceleration at forward and reverse shocks

For parameters that are applicable to the conditions at young supernova remnants, we present results of 2D3V particle-in-cell simulations of a non-relativistic plasma shock with a large-scale perpendicular magnetic field inclined at 45-deg angle to the simulation plane to approximate 3D physics. We developed an improved clean setup that uses the collision of two plasma slabs with different density and velocity, leading to the development of two distinctive shocks and a contact discontinuity. The shock formation is mediated by Weibel-type filamentation instabilities that generate magnetic turbulence. Cyclic reformation is observed in both shocks with similar period, for which we note global variations on account of shock rippling and local variations arising from turbulent current filaments. The shock rippling occurs on spatial and temporal scales given by gyro-motions of shock-reflected ions. The drift motion of electrons and ions is not a gradient drift, but commensurates with E x B drift. We observe a stable suprathermal tail in the ion spectra, but no electron acceleration because the amplitude of Buneman modes in the shock foot is insufficient for trapping relativistic electrons. We see no evidence of turbulent reconnection. A comparison with other 2D simulation results suggests that the plasma beta and the ion-to-electron mass ratio are not decisive for efficient electron acceleration, but pre-acceleration efficacy might be reduced with respect to the 2D results once three-dimensional effects are fully accounted for. Other microphysical factors may also be at play to limit the amplitude of Buneman waves or prevent return of electrons to the foot region.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Collisionless Weibel shocks and electron acceleration in gamma-ray bursts

A study of collisionless external shocks in gamma-ray bursts is presented. The shock structure, electromagnetic fields, and process of electron acceleration are assessed by performing a self-consistent 3D particle-in-cell (PIC) simulation. In accordance with hydrodynamic shock systems, the shock consists of a reverse shock (RS) and forward shock (FS) separated by a contact discontinuity (CD). The development and structure are controlled by the ion Weibel instability. The ion filaments are sources of strong transverse electromagnetic fields at both sides of the double shock structure over a length of 30 - 100 ion skin depths. Electrons are heated up to a maximum energy $ε_{\rm ele}\approx \sqrt{ε_{\rm b}}$, where $ε$ is the energy normalized to the total incoming energy. Jet electrons are trapped in the RS transition region due to the presence of an ambipolar electric field and reflection by the strong transverse magnetic fields in the shocked region. In a process similar to shock surfing acceleration (SSA) for ions, electrons experience drift motion and acceleration by ion filament transverse electric fields in the plane perpendicular to the shock propagation direction. Ultimately accelerated jet electrons are convected back into the upstream.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Recollimation Shocks in Magnetized Relativistic Jets

We have performed two-dimensional special-relativistic magnetohydrodynamic simulations of non-equilibrium over-pressured relativistic jets in cylindrical geometry. Multiple stationary recollimation shock and rarefaction structures are produced along the jet by the nonlinear interaction of shocks and rarefaction waves excited at the interface between the jet and the surrounding ambient medium. Although initially the jet is kinematically dominated, we have considered axial, toroidal and helical magnetic fields to investigate the effects of different magnetic-field topologies and strengths on the recollimation structures. We find that an axial field introduces a larger effective gas-pressure and leads to stronger recollimation shocks and rarefactions, resulting in larger flow variations. The jet boost grows quadratically with the initial magnetic field. On the other hand, a toroidal field leads to weaker recollimation shocks and rarefactions, modifying significantly the jet structure after the first recollimation rarefaction and shock. The jet boost decreases systematically. For a helical field, instead, the behaviour depends on the magnetic pitch, with a phenomenology that ranges between the one seen for axial and toroidal magnetic fields, respectively. In general, however, a helical magnetic field yields a more complex shock and rarefaction substructure close to the inlet that significantly modifies the jet structure. The differences in shock structure resulting from different field configurations and strengths may have observable consequences for disturbances propagating through a stationary recollimation shock.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Magnetic Field Amplification and Saturation in Turbulence Behind a Relativistic Shock

We have investigated via two-dimensional relativistic MHD simulations the long-term evolution of turbulence created by a relativistic shock propagating through an inhomogeneous medium. In the postshock region, magnetic field is strongly amplified by turbulent motions triggered by preshock density inhomogeneities. Using a long-simulation box we have followed the magnetic-field amplification until it is fully developed and saturated. The turbulent velocity is sub-relativistic even for a strong shock. Magnetic-field amplification is controled by the turbulent motion and saturation occurs when the magnetic energy is comparable to the turbulent kinetic energy. Magnetic-field amplification and saturation depend on the initial strength and direction of the magnetic field in the preshock medium, and on the shock strength. If the initial magnetic field is perpendicular to the shock normal, the magnetic field is first compressed at the shock and then can be amplified by turbulent motion in the postshock region. Saturation occurs when the magnetic energy becomes comparable to the turbulent kinetic energy in the postshock region. If the initial magnetic field in the preshock medium is strong, the postshock region becomes turbulent but significant field amplification does not occur. If the magnetic energy after shock compression is larger than the turbulent kinetic energy in the postshock region, significant field amplification does not occur. We discuss possible applications of our results to gamma-ray bursts and active galactic nuclei.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Spatial Growth of the Current-Driven Instability in Relativistic Jets

We have investigated the influence of velocity shear and a radial density profile on the spatial development of the current driven kink instability along helically magnetized relativistic jets via three-dimensional relativistic magnetohydrodynamic simulations. In this study, we use a non-periodic computational box, the jet flow is initially established across the computational grid, and a precessional perturbation at the inlet triggers growth of the kink instability. If the velocity shear radius is located inside the characteristic radius of the helical magnetic field, a static non-propagating current driven kink is excited as the perturbation propagates down the jet. Temporal growth disrupts the initial flow across the computational grid not too far from the inlet. On the other hand, if the velocity shear radius is outside the characteristic radius of the helical magnetic field, the kink is advected with the flow and grows spatially down the jet. In this case flow is maintained to much larger distances from the inlet. The effect of different radial density profiles is more subtle. When the density increases with radius, the kink appears to saturate by the end of the simulation without apparent disruption of the helical twist. This behavior suggests that relativistic jets consisting of a tenuous spine surrounded by a denser medium with a velocity shear radius outside the radius of maximum toroidal magnetic field have a relatively stable configuration.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Three-Dimensional Relativistic Magnetohydrodynamic Simulations of Current-Driven Instability. III. Rotating Relativistic Jets

We have investigated the influence of jet rotation and differential motion on the linear and nonlinear development of the current-driven (CD) kink instability of force-free helical magnetic equilibria via three-dimensional relativistic magnetohydrodynamic simulations. In this study, we follow the temporal development within a periodic computational box. Displacement of the initial helical magnetic field leads to the growth of the CD kink instability. We find that, in accord with linear stability theory, the development of the instability depends on the lateral distribution of the poloidal magnetic field. If the poloidal field significantly decreases outwards from the axis, the initial small perturbations grow strongly, and if multiple wavelengths are excited non-linear interaction eventually disrupts the initial cylindrical configuration. When the profile of the poloidal field is shallow, the instability develops slowly and eventually saturates. We briefly discuss implications of our findings for Poynting dominated jets.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Three-Dimensional Relativistic Magnetohydrodynamic Simulations of Current-Driven Instability with A Sub-Alfvenic Jet: Temporal Properties

We have investigated the influence of a velocity shear surface on the linear and non-linear development of the CD kink instability of force-free helical magnetic equilibria in 3D. In this study we follow the temporal development within a periodic computational box and concentrate on flows that are sub-Alfvenic on the cylindrical jet's axis. Displacement of the initial force-free helical magnetic field leads to the growth of CD kink instability. We find that helically distorted density structure propagates along the jet with speed and flow structure dependent on the radius of the velocity shear surface relative to the characteristic radius of the helically twisted force-free magnetic field. At small velocity shear surface radius the plasma flows through the kink with minimal kink propagation speed. The kink propagation speed increases as the velocity shear radius increases and the kink becomes more embedded in the plasma flow. A decreasing magnetic pitch profile and faster flow enhance the influence of velocity shear. Simulations show continuous transverse growth in the nonlinear phase of the instability. The growth rate of the CD kink instability and the nonlinear behavior also depend on the velocity shear surface radius and flow speed, and the magnetic pitch radial profile. Larger velocity shear radius leads to slower linear growth, makes a later transition to the nonlinear stage, and with larger maximum amplitude than occur for a static plasma column. However, when the velocity shear radius is much greater than the characteristic radius of the helical magnetic field, linear and non-linear development can be similar to the development of a static plasma column.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Magnetic Field Amplification by Turbulence in A Relativistic Shock Propagating through An Inhomogeneous Medium

We perform two-dimensional relativistic magnetohydrodynamic simulations of a mildly relativistic shock propagating through an inhomogeneous medium. We show that the postshock region becomes turbulent owing to preshock density inhomogeneity, and the magnetic field is strongly amplified due to the stretching and folding of field lines in the turbulent velocity field. The amplified magnetic field evolves into a filamentary structure in two-dimensional simulations. The magnetic energy spectrum is flatter than the Kolmogorov spectrum and indicates that a so-called small-scale dynamo is occurring in the postshock region. We also find that the amount of magnetic-field amplification depends on the direction of the mean preshock magnetic field, and the time scale of magnetic-field growth depends on the shock strength.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Three-Dimensional Relativistic Magnetohydrodynamic Simulations of Current-Driven Instability. II. Relaxation of Pulsar Wind Nebula

We have investigated the relaxation of a hydrostatic hot plasma column containing toroidal magnetic field by the Current-Driven (CD) kink instability as a model of pulsar wind nebulae. In our simulations the CD kink instability is excited by a small initial velocity perturbation and develops turbulent structure inside the hot plasma column. We demonstrate that, as envisioned by Begelman, the hoop stress declines and the initial gas pressure excess near the axis decreases. The magnetization parameter σ, the ratio of the Poynting to the kinetic energy flux, declines from an initial value of 0.3 to about 0.01 when the CD kink instability saturates. Our simulations demonstrate that axisymmetric models strongly overestimate the elongation of the pulsar wind nebulae. Therefore, the previous requirement for an extremely low pulsar wind magnetization can be abandoned. The observed structure of the pulsar wind nebulae do not contradict the natural assumption that the magnetic energy flux still remains a good fraction of the total energy flux after dissipation of alternating fields.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

A method of incorporating general relativity in electromagnetic particle-in-cell code

An algorithm is presented that incorporates the tensor form of Maxwell's equations in a general relativistic electromagnetic particle-in-cell code. The code simplifies to Schwartzschild space-time for a non-spinning central mass. The particle advance routine uses a fourth-order Runge-Kutta algorithm to integrate the four-velocity form of Lorentz force. The current density is calculated using the curved space-time of the metric.

0 citations0 reviewsNo code linkedOpen access