Source author record

Takanobu Amano

Takanobu Amano appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

22works
7topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

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

Published work

22 published item(s)

preprint2022arXiv

Theory of Electron Injection at Oblique Shock of Finite Thickness

A theory of electron injection into diffusive shock acceleration (DSA) for the generation of cosmic-ray electrons at collisionless shocks is presented. We consider a recently proposed particle acceleration mechanism called stochastic shock drift acceleration (SSDA). We find that SSDA may be understood as a diffusive particle acceleration mechanism at an oblique shock of finite thickness. More specifically, it is described by a solution to the diffusion-convection equation for particles with the characteristic diffusion length comparable to the shock thickness. On the other hand, the same equation yields the standard DSA if the diffusion length is much longer than the thickness. Although SSDA predicts, in general, a spectral index steeper than DSA, it is much more efficient for low-energy electron acceleration and is favorable for injection. The injection threshold energy corresponds to the transition energy between the two different regimes. It is of the order of $0.1\text{-}1$ MeV in typical interstellar and interplanetary conditions if the dissipation scale of turbulence around the shock is determined by the ion inertial length. The electron injection is more efficient at high $M_\textrm{A} / \cos θ_{Bn}$ where $M_\textrm{A}$ and $θ_{Bn}$ are the Alfvén Mach number and the shock obliquity. The theory suggests that efficient acceleration of electrons to ultra-relativistic energies will be more easily realized at high-Mach-number young supernova remnant shocks, but not at weak or moderate shocks in the heliosphere unless the upstream magnetic field is nearly perpendicular to the shock normal.

preprint2021arXiv

Critical comparison of collisionless fluid models: Nonlinear simulations of parallel firehose instability

Two different fluid models for collisionless plasmas are compared. One is based on the classical Chew-Goldberger-Low (CGL) model that includes a finite Larmor radius (FLR) correction and the Landau closure for the longitudinal mode. Another one takes into account the effect of cyclotron resonance in addition to Landau resonance, which is referred to as the cyclotron resonance closure (CRC) model. While the linear property of the parallel firehose instability is better described by the CGL model, the electromagnetic ion cyclotron instability driven unstable by the cyclotron resonance is reproduced only by the CRC model. Nonlinear simulation results for the parallel firehose instability performed with the two models are also discussed. Although the linear and quasilinear isotropization phases are consistent with theory in both models, long-term behaviors may be substantially different. The final state obtained by the CRC model may be reasonably understood in terms of the marginal stability condition. In contrast, the lack of cyclotron damping in the CGL model makes it rather difficult to predict the long-term behavior with a simple physical argument. This suggests that incorporating the collisionless damping both for longitudinal and transverse modes is crucial for a nonlinear fluid simulation model of collisionless plasmas.

preprint2021arXiv

Mildly relativistic magnetized shocks in electron-ion plasmas -- II. Particle acceleration and heating

Particle acceleration and heating at mildly relativistic magnetized shocks in electron-ion plasma are investigated with unprecedentedly high-resolution two-dimensional particle-in-cell simulations that include ion-scale shock rippling. Electrons are super-adiabatically heated at the shock, and most of the energy transfer from protons to electrons takes place at or downstream of the shock. We are the first to demonstrate that shock rippling is crucial for the energization of electrons at the shock. They remain well below equipartition with the protons. The downstream electron spectra are approximately thermal with a limited supra-thermal power-law component. Our results are discussed in the context of wakefield acceleration and the modelling of electromagnetic radiation from blazar cores.

preprint2021arXiv

Particle Acceleration by Pickup Process Upstream of Relativistic Shocks

Particle acceleration at magnetized purely perpendicular relativistic shocks in electron-ion plasmas are studied by means of two-dimensional particle-in-cell simulations. Magnetized shocks with the upstream bulk Lorentz factor $γ_1 \gg 1$ are known to emit intense electromagnetic waves from the shock front, which induce electrostatic plasma waves (wakefield) and transverse filamentary structures in the upstream region via the stimulated/induced Raman scattering and the filamentation instability, respectively. The wakefield and filaments inject a fraction of incoming particles into a particle acceleration process, in which particles are once decoupled from the upstream bulk flow by the wakefield, and are piked up again by the flow. The picked-up particles are accelerated by the motional electric field. The maximum attainable Lorentz factor is estimated as $γ_{max,e} \sim αγ_1^3$ for electrons and $γ_{max,i} \sim (1+m_eγ_1/m_i)γ_1^2$ for ions, where $α\sim 10$ is determined from our simulation results. $α$ can increase up to $γ_1$ for weakly magnetized shock if $γ_1$ is sufficiently large. This result indicates that highly relativistic astrophysical shocks such as external shocks of gamma-ray bursts can be an efficient particle accelerator.

preprint2020arXiv

Kinetic simulations of nonrelativistic perpendicular shocks of young supernova remnants. III. Magnetic reconnection

Fully kinetic two-dimensional particle-in-cell simulations are used to study electron acceleration at high-Mach-number nonrelativistic perpendicular shocks. SNR shocks are mediated by the Weibel instability which is excited because of an interaction between shock-reflected and upstream ions. Nonlinear evolution of the Weibel instability leads to the formation of current sheets. At the turbulent shock ramp the current sheets decay through magnetic reconnection. The number of reconnection sites strongly depends on the ion-to-electron mass ratio and the Alfvénic Mach number of the simulated shock. Electron acceleration is observed at locations where magnetic reconnection operates. For the highest mass ratios almost all electrons are involved in magnetic reconnection, which makes the magnetic reconnection the dominant acceleration process for electrons at these shocks. We discuss the relevance of our results for 3D systems with realistic ion-to-electron mass ratio.

preprint2020arXiv

Mildly relativistic magnetized shocks in electron-ion plasmas I. Electromagnetic shock structure

Mildly relativistic shocks in magnetized electron-ion plasmas are investigated with 2D kinetic particle-in-cell simulations of unprecedentedly high resolution and large scale for conditions that may be found at internal shocks in blazar cores. Ion-scale effects cause corrugations along the shock surface whose properties somewhat depend on the configuration of the mean perpendicular magnetic field, that is either in or out of the simulation plane. We show that the synchrotron maser instability persists to operate in mildly relativistic shocks in agreement with theoretical predictions and produces coherent emission of upstream-propagating electromagnetic waves. Shock front ripples are excited in both mean-field configurations and they engender effective wave amplification. The interaction of these waves with upstream plasma generates electrostatic wakefields.

preprint2016arXiv

A High-Order Relativistic Two-Fluid Electrodynamic Scheme with Consistent Reconstruction of Electromagnetic Fields and a Multidimensional Riemann Solver for Electromagnetism

In various astrophysics settings it is common to have a two-fluid relativistic plasma that interacts with the electromagnetic field. While it is common to ignore the displacement current in the ideal, classical magnetohydrodynamic limit, when the flows become relativistic this approximation is less than absolutely well-justified. In such a situation, it is more natural to consider a positively charged fluid made up of positrons or protons interacting with a negatively charged fluid made up of electrons. The two fluids interact collectively with the full set of Maxwell's equations. As a result, a solution strategy for that coupled system of equations is sought and found here. Our strategy extends to higher orders, providing increasing accuracy. Three important innovations are reported here. In our first innovation, the magnetic field within each zone is reconstructed in a divergence-free fashion while the electric field within each zone is reconstructed in a form that is consistent with Gauss' law. In our second innovation, a multidimensionally upwinded strategy is presented which ensures that the magnetic field can be updated via a discrete interpretation of Faraday's law and the electric field can be updated via a discrete interpretation of the generalized Ampere's law. Our third innovation consists of an efficient design for several popular IMEX schemes so that they provide strong coupling between the finite-volume-based fluid solver and the electromagnetic fields at high order. Several accuracy analyses are presented showing that our method meets its design accuracy in the MHD limit as well as in the limit of electromagnetic wave propagation. Several stringent test problems are also presented. We also present a relativistic version of the GEM problem, which shows that our algorithm can successfully adapt to challenging problems in high energy astrophysics.

preprint2016arXiv

A new framework for magnetohydrodynamic simulations with anisotropic pressure

We describe a new theoretical and numerical framework of the magnetohydrodynamic simulation incorporated with an anisotropic pressure tensor, which can play an important role in a collisionless plasma. A classical approach to handle the anisotropy is based on the double adiabatic approximation assuming that a pressure tensor is well described only by the components parallel and perpendicular to the local magnetic field. This gyrotropic assumption, however, fails around a magnetically neutral region, where the cyclotron period may get comparable to or even longer than a dynamical time in a system, and causes a singularity in the mathematical expression. In this paper, we demonstrate that this singularity can be completely removed away by the combination of direct use of the 2nd-moment of the Vlasov equation and an ingenious gyrotropization model. Numerical tests also verify that the present model properly reduces to the standard MHD or the double adiabatic formulation in an asymptotic manner under an appropriate limit.

preprint2016arXiv

A Second-order Divergence-constrained Multidimensional Numerical Scheme for Relativistic Two-Fluid Electrodynamics

A new multidimensional simulation code for relativistic two-fluid electrodynamics (RTFED) is described. The basic equations consist of the full set of Maxwell's equations coupled with relativistic hydrodynamic equations for separate two charged fluids, representing the dynamics of either an electron-positron or an electron-proton plasma. It can be recognized as an extension of conventional relativistic magnetohydrodynamics (RMHD). Finite resistivity may be introduced as a friction between the two species, which reduces to resistive RMHD in the long wavelength limit without suffering from a singularity at infinite conductivity. A numerical scheme based on HLL (Harten-Lax-Van Leer) Riemann solver is proposed that exactly preserves the two divergence constraints for Maxwell's equations simultaneously. Several benchmark problems demonstrate that it is capable of describing RMHD shocks/discontinuities at long wavelength limit, as well as dispersive characteristics due to the two-fluid effect appearing at small scales. This shows that the RTFED model is a promising tool for high energy astrophysics application.

preprint2015arXiv

Divergence-free Approximate Riemann Solver for the Quasi-neutral Two-fluid Plasma Model

A numerical method for the quasi-neutral two-fluid (QNTF) plasma model is described. The basic equations are ion and electron fluid equations and the Maxwell equations without displacement current. The neglect of displacement current is consistent with the assumption of charge neutrality. It thus reduces to the ideal magnetohydrodynamic (MHD) equations in the long wavelength limit, but the two-fluid effect appearing at ion and electron inertial scales is fully taken into account. It is shown that the basic equations may be rewritten in a form that has formally the same structure as the MHD equations. The total mass, momentum, and energy are all written in the conservative form. A new three-dimensional numerical simulation code has been developed for the QNTF equations. The HLL (Harten-Lax-van Leer) approximate Riemann solver combined with the upwind constrained transport (UCT) scheme is applied. The method was originally developed for MHD (Londrillo & Del Zanna, 2004), but works quite well for the present model as well. The simulation code is able to capture sharp multidimensional discontinuities as well as dispersive waves arising from the two-fluid effect at small scales without producing $\nabla \cdot \mathbf{B}$ errors. It is well known that conventional Hall-MHD codes often suffer a numerical stability issue associated with short wavelength whistler waves. On the other hand, since finite electron inertia introduces an upper bound to the phase speed of whistler waves in the present model, our code is free from the issue even without explicit dissipation terms or implicit time integration. Numerical experiments have confirmed that there is no need to resolve characteristic time scales such as plasma frequency or cyclotron frequency for numerical stability. Consequently, the QNTF model offers a better alternative to the Hall-MHD or fully electromagnetic two-fluid models.

preprint2014arXiv

A Robust Method for Handling Low Density Regions in Hybrid Simulations for Collisionless Plasmas

A robust method to handle vacuum and near vacuum regions in hybrid simulations for space and astrophysical plasmas is presented. The conventional hybrid simulation model dealing with kinetic ions and a massless charge-neutralizing electron fluid is known to be susceptible to numerical instability due to divergence of the whistler-mode wave dispersion, as well as division-by-density operation in regions of low density. Consequently, a pure vacuum region is not allowed to exist in the simulation domain unless some ad hoc technique is used. To resolve this difficulty, an alternative way to introduce finite electron inertia effect is proposed. Contrary to the conventional method, the proposed one introduces a correction to the electric field rather than the magnetic field. It is shown that the generalized Ohm's law correctly reduces to Laplace's equation in a vacuum which therefore does not involve any numerical problems. In addition, a variable ion-to-electron mass ratio is introduced to reduce the phase velocity of high frequency whistler waves at low density regions so that the stability condition is always satisfied. It is demonstrated that the proposed model is able to handle near vacuum regions generated as a result of nonlinear self-consistent development of the system, as well as pure vacuum regions set up at the initial condition, without losing the advantages of the standard hybrid code.

preprint2014arXiv

First-principles simulations of electrostatic interactions between dust grains

We investigated the electrostatic interaction between two identical dust grains of an infinite mass immersed in homogeneous plasma by employing first-principles N-body simulations combined with the Ewald method. We specifically tested the possibility of an attractive force due to overlapping Debye spheres (ODSs), as was suggested by Resendes et al. (1998). Our simulation results demonstrate that the electrostatic interaction is repulsive and even stronger than the standard Yukawa potential. We showed that the measured electric field acting on the grain is highly consistent with a model electrostatic potential around a single isolated grain that takes into account a correction due to the orbital motion limited theory. Our result is qualitatively consistent with the counterargument suggested by Markes and Williams (2000), indicating the absence of the ODS attractive force.

preprint2013arXiv

Stability of Cosmic Ray Modified Shocks: Two-Fluid Approach

The stability of the cosmic ray modified shock (CRMS) is studied by means of numerical simulations. Owing to the nonlinear feedback of cosmic-ray (CR) acceleration, a downstream state of the modified shock can no longer be uniquely determined for given upstream parameters. It is known that up to three distinct solutions exist, which are characterized by the CR production efficiency as "efficient", "intermediate" and "inefficient" branches. The stability of these solutions is investigated by performing direct time-dependent simulations of a two-fluid model. It is found that both the efficient and inefficient branches are stable even against a large-amplitude perturbation, while the intermediate one is always unstable and evolves into the inefficient state as a result of nonlinear time development. This bistable feature is robust in a wide range of parameters and does not depend on the injection model. Fully nonlinear time evolution of a hydrodynamic shock with injection results in the least efficient state in terms of the CR production, consistent with the bistable feature. This suggests that the CR production efficiency at supernova remnant shocks may be lower than previously discussed in the framework of the nonlinear shock acceleration theory considering the efficient solution of the CRMS.

preprint2013arXiv

The Role of Superluminal Electromagnetic Waves in Pulsar Wind Termination Shocks

The dynamics of a standing shock front in a Poynting-flux dominated relativistic flow is investigated by using a one-dimensional, relativistic, two-fluid simulation. An upstream flow containing a circularly polarized, sinusoidal magnetic shear wave is considered, mimicking a wave driven by an obliquely rotating pulsar. It is demonstrated that this wave is converted into large amplitude electromagnetic waves with superluminal phase speeds by interacting with the shock when the shock-frame frequency of the wave exceeds the proper plasma frequency. The superluminal waves propagate in the upstream, modify the shock structure substantially, and form a well-developed precursor region ahead of a subshock. Dissipation of Poynting flux occurs in the precursor as well as in the downstream region through a parametric instability driven by the superluminal waves. The Poynting flux remaining in the downstream region is carried entirely by the superluminal waves. The downstream plasma is therefore an essentially unmagnetized, relativistically hot plasma with a non-relativistic flow speed, as suggested by observations of pulsar wind nebulae.

preprint2012arXiv

Electron Accelerations at High Mach Number Shocks: Two-Dimensional Particle-In-Cell Simulations in Various Parameter Regimes

Electron accelerations at high Mach number collision-less shocks are investigated by means of two-dimensional electromagnetic Particle-in-Cell simulations with various Alfven Mach numbers, ion-to-electron mass ratios, and the upstream electron beta_e (the ratio of the thermal pressure to the magnetic pressure). We found electrons are effectively accelerated at a super-high Mach number shock (MA~30) with a mass ratio of M/m=100 and beta_e=0.5. The electron shock surfing acceleration is an effective mechanism for accelerating the particles toward the relativistic regime even in two dimensions with the large mass ratio. Buneman instability excited at the leading edge of the foot in the super-high Mach number shock results in a coherent electrostatic potential structure. While multi-dimensionality allows the electrons to escape from the trapping region, they can interact with the strong electrostatic field several times. Simulation runs in various parameter regimes indicate that the electron shock surfing acceleration is an effective mechanism for producing relativistic particles in extremely-high Mach number shocks in supernova remnants, provided that the upstream electron temperature is reasonably low.

preprint2012arXiv

Multi-Moment Advection scheme for Vlasov simulations

We present a new numerical scheme for solving the advection equation and its application to the Vlasov simulation. The scheme treats not only point values of a profile but also its zeroth to second order piecewise moments as dependent variables, and advances them on the basis of their governing equations. We have developed one- and two-dimensional schemes and show that they provide quite accurate solutions compared to other existing schemes with the same memory usage. The two-dimensional scheme can solve the solid body rotation problem of a gaussian profile with little numerical diffusion. This is a very important property for Vlasov simulations of magnetized plasma. The application of the scheme to the electromagnetic Vlasov simulation of collisionless shock waves is presented as a benchmark test.

preprint2012arXiv

Multi-moment advection scheme in three dimension for Vlasov simulations of magnetized plasma

We present an extension of the multi-moment advection scheme (Minoshima et al., 2011, J. Comput. Phys.) to the three-dimensional case, for full electromagnetic Vlasov simulations of magnetized plasma. The scheme treats not only point values of a profile but also its zeroth to second order piecewise moments as dependent variables, and advances them on the basis of their governing equations. Similar to the two-dimensional scheme, the three-dimensional scheme can accurately solve the solid body rotation problem of a gaussian profile with little numerical dispersion or diffusion. This is a very important property for Vlasov simulations of magnetized plasma. We apply the scheme to electromagnetic Vlasov simulations. Propagation of linear waves and nonlinear evolution of the electron temperature anisotropy instability are successfully simulated with a good accuracy of the energy conservation.

preprint2011arXiv

Molecular and Atomic Gas toward HESS J1745-303 in the Galactic Center: Further Support for the Hadronic Scenario

We have compared the TeV gamma-rays with the new 12CO J=2-1 data toward HESS J1745-303 in the Galactic center and confirmed that the molecular gas MG358.9-0.5 toward (l, b)=(358.9, -0.5 at VLSR=-100-0 km s-1 shows a reasonable positional agreement with the primary peak (northern part) of the gamma-ray source. For the southern part of HESS J1745-303, we see no CO counterpart, whereas the HI gas in the Parkes 21 cm HI dataset shows a possible counterpart to the gamma-ray source. This HI gas may be optically thick as supported by the HI line shape similar to the optically thick 12CO. We estimate the total mass of interstellar protons including both the molecular and atomic gas to be 2x10^6 Mo and the cosmic-ray proton energy to be 6x10^{48} ergs in the hadronic scenario. We discuss possible origins of the cosmic-ray protons including the nearby SNR G359.1-0.5. The SNR may be able to explain the northern gamma-ray source but the southern source seems to be too far to be energized by the SNR. As an alternative, we argue that the second-order Fermi acceleration in the inter-clump space surrounded by randomly moving high-velocity clumps may offer a possible mechanism to accelerate protons. The large turbulent motion with velocity dispersion of ~15 km s-1 has energy density two orders of magnitude higher than in the solar vicinity and is viable as the energy source.

preprint2011arXiv

Multi-Moment Advection scheme for Vlasov simulations

We present a new numerical scheme for solving the advection equation and its application to Vlasov simulations. The scheme treats not only point values of a profile but also its zeroth to second order piecewise moments as dependent variables, for better conservation of the information entropy. We have developed one- and two-dimensional schemes and show that they provide quite accurate solutions within reasonable usage of computational resources compared to other existing schemes. The two-dimensional scheme can accurately solve the solid body rotation problem of a gaussian profile for more than hundred rotation periods with little numerical diffusion. This is crucially important for Vlasov simulations of magnetized plasmas. Applications of the one- and two-dimensional schemes to electrostatic and electromagnetic Vlasov simulations are presented with some benchmark tests.

preprint2011arXiv

Stochastic Acceleration of Cosmic Rays in the Central Molecular Zone of the Galaxy

Particle acceleration in the inner ~ 200 pc of the Galaxy is discussed, where diffuse TeV gamma-rays have been detected by the High Energy Stereoscopic System (HESS) observation. The diffuse gamma-ray emission has a strong correlation with molecular clouds with large velocity dispersion, indicating the presence of turbulence. It is argued that the turbulence may contribute to the acceleration of cosmic rays via stochastic acceleration. The stochastic acceleration may energize cosmic-ray protons up to ~ 100 TeV and electrons to ~ 1 TeV in a relatively tenuous medium. The difference in the efficiency between protons and electrons supports the hadronic scenario of the diffuse TeV gamma-ray emission.

preprint2010arXiv

A critical Mach number for electron injection in collisionless shocks

Electron acceleration in collisionless shocks with arbitrary magnetic field orientations is discussed. It is shown that the injection of thermal electrons into diffusive shock acceleration process is achieved by an electron beam with a loss-cone in velocity space that is reflected back upstream from the shock through shock drift acceleration mechanism. The electron beam is able to excite whistler waves which can scatter the energetic electrons themselves when the Alfven Mach number of the shock is sufficiently high. A critical Mach number for the electron injection is obtained as a function of upstream parameters. The application to supernova remnant shocks is discussed.

preprint2010arXiv

Structure of a strong supernova shock wave and rapid electron acceleration confined in its transition region

A new rapid energization process within a supernova shock transition region (STR) is reported by utilizing numerical simulation. Although the scale of a STR as a main dissipation region is only several hundreds of thousands km, several interesting structures are found relating to generation of a root of the energetic particles. The nonlinear evolution of plasma instabilities lead to a dynamical change in the ion phase space distribution which associates with change of the field properties. As a result, different types of large-amplitude field structures appear. One is the leading wave packet and another is a series of magnetic solitary humps. Each field structure has a microscopic scale (~ the ion inertia length). Through the multiple nonlinear scattering between these large-amplitude field structures, electrons are accelerated directly. Within a STR, quick thermalization realizes energy equipartition between the ion and electron, hot electrons play an important role in keeping these large-amplitude field structures on the ion-acoustic mode. The hot electron shows non-Maxwellian distribution and could be the seed of further non-thermal population. The "shock system", where fresh incoming and reflected ions are supplied constantly, play an essential role in our result. With a perpendicular shock geometry, the maximum energy of the electron is estimated by equating a width of the STR to a length of the Larmor radius of the energetic electron. Under some realistic condition of M_A = 170 and ω_{pe}/Ω_{ce} = 120, maximum energy is estimated to ~ 10 MeV at an instant only within the STR.