-
High-power laser experiment forming a supercritical collisionless shock in a magnetized uniform plasma at rest
Authors:
Ryo Yamazaki,
S. Matsukiyo,
T. Morita,
S. J. Tanaka,
T. Umeda,
K. Aihara,
M. Edamoto,
S. Egashira,
R. Hatsuyama,
T. Higuchi,
T. Hihara,
Y. Horie,
M. Hoshino,
A. Ishii,
N. Ishizaka,
Y. Itadani,
T. Izumi,
S. Kambayashi,
S. Kakuchi,
N. Katsuki,
R. Kawamura,
Y. Kawamura,
S. Kisaka,
T. Kojima,
A. Konuma
, et al. (29 additional authors not shown)
Abstract:
We present a new experimental method to generate quasi-perpendicular supercritical magnetized collisionless shocks. In our experiment, ambient nitrogen (N) plasma is at rest and well-magnetized, and it has uniform mass density. The plasma is pushed by laser-driven ablation aluminum (Al) plasma. Streaked optical pyrometry and spatially resolved laser collective Thomson scattering clarify structures…
▽ More
We present a new experimental method to generate quasi-perpendicular supercritical magnetized collisionless shocks. In our experiment, ambient nitrogen (N) plasma is at rest and well-magnetized, and it has uniform mass density. The plasma is pushed by laser-driven ablation aluminum (Al) plasma. Streaked optical pyrometry and spatially resolved laser collective Thomson scattering clarify structures of plasma density and temperatures, which are compared with one-dimensional particle-in-cell simulations. It is indicated that just after the laser irradiation, the Al plasma is magnetized by a self-generated Biermann battery field, and the plasma slaps the incident N plasma. The compressed external field in the N plasma reflects N ions, leading to counter-streaming magnetized N flows. Namely we identify the edge of the reflected N ions. Such interacting plasmas form a magnetized collisionless shock.
△ Less
Submitted 7 February, 2022; v1 submitted 19 January, 2022;
originally announced January 2022.
-
Mach number and plasma beta dependence of the ion temperature perpendicular to the external magnetic field in the transition region of perpendicular collisionless shocks
Authors:
Ryo Yamazaki,
Ayato Shinoda,
Takayuki Umeda,
Shuichi Matsukiyo
Abstract:
Ion temperature anisotropy is a common feature for (quasi-)perpendicular collisionless shocks. By using two-dimensional full particle simulations, it is shown, that the ion temperature component perpendicular to the shock magnetic field at the shock foot region is proportional to the square of the Alfven Mach number divided by the plasma beta. This result is also explained by a simple analytical a…
▽ More
Ion temperature anisotropy is a common feature for (quasi-)perpendicular collisionless shocks. By using two-dimensional full particle simulations, it is shown, that the ion temperature component perpendicular to the shock magnetic field at the shock foot region is proportional to the square of the Alfven Mach number divided by the plasma beta. This result is also explained by a simple analytical argument, in which the reflected ions get energy from upstream plasma flow. By comparing our analytic and numerical results, it is also confirmed that the fraction of the reflected ions hardly depends on the plasma beta and the Alfven Mach number when the square of the Alfven Mach number divided by the plasma beta is larger than about 20.
△ Less
Submitted 19 November, 2019;
originally announced November 2019.
-
Full particle-in-cell simulation of the interaction between two plasmas for laboratory experiments on the generation of magnetized collisionless shocks with high-power lasers
Authors:
T. Umeda,
R. Yamazaki,
Y. Ohira,
N. Ishizaka,
S. Kakuchi,
Y. Kuramitsu,
S. Matsukiyo,
I. Miyata,
T. Morita,
Y. Sakawa,
T. Sano,
S. Sei,
S. J. Tanaka,
H. Toda,
S. Tomita
Abstract:
A preliminary numerical experiment is conducted for laboratory experiments on the generation of magnetized collisionless shocks with high-power lasers by using one-dimensional particle-in-cell simulation. The present study deals with the interaction between a moving Aluminum plasma and a Nitrogen plasma at rest. In the numerical experiment, the Nitrogen plasma is unmagnetized or magnetized by a we…
▽ More
A preliminary numerical experiment is conducted for laboratory experiments on the generation of magnetized collisionless shocks with high-power lasers by using one-dimensional particle-in-cell simulation. The present study deals with the interaction between a moving Aluminum plasma and a Nitrogen plasma at rest. In the numerical experiment, the Nitrogen plasma is unmagnetized or magnetized by a weak external magnetic field. Since the previous study suggested the generation of spontaneous magnetic field in the piston (Aluminum) plasma due to the Biermann battery, the effect of the magnetic field is of interest. Sharp jumps of electron density and magnetic field are observed around the interface between the two plasmas as long as one of the two plasmas is magnetized, which indicates the formation of tangential electron-magneto-hydro-dynamic discontinuity. When the Aluminum plasma is magnetized, strong compression of both density and magnetic field takes place in the pure Aluminum plasma during the gyration of Nitrogen ions in the Aluminum plasma region. The formation of a shock downstream is indicated from the shock jump condition. The result suggests that the spontaneous magnetic field in the piston (Aluminum) plasma plays an essential role in the formation of a perpendicular collisionless shock.
△ Less
Submitted 8 February, 2019;
originally announced February 2019.
-
On the Boris solver in particle-in-cell simulation
Authors:
Seiji Zenitani,
Takayuki Umeda
Abstract:
A simple form of the Boris solver in particle-in-cell (PIC) simulation is proposed. It employs an exact solution of the Lorentz-force part, and it is equivalent to the Boris solver with a gyrophase correction. As a favorable property for stable schemes, this form preserves a volume in the phase space. Numerical tests of the Boris solvers are conducted by test-particle simulations and by PIC simula…
▽ More
A simple form of the Boris solver in particle-in-cell (PIC) simulation is proposed. It employs an exact solution of the Lorentz-force part, and it is equivalent to the Boris solver with a gyrophase correction. As a favorable property for stable schemes, this form preserves a volume in the phase space. Numerical tests of the Boris solvers are conducted by test-particle simulations and by PIC simulations. The proposed form provides better accuracy than the popular form, while it only requires few additional computation time.
△ Less
Submitted 31 October, 2018; v1 submitted 12 September, 2018;
originally announced September 2018.
-
Some remarks on the diffusion regions in magnetic reconnection
Authors:
Seiji Zenitani,
Takayuki Umeda
Abstract:
The structure of the diffusion regions in antiparallel magnetic reconnection is investigated by means of a theory and a Vlasov simulation. The magnetic diffusion is considered as relaxation to the frozen-in state, which depends on a reference velocity field. A field-aligned component of the frozen-in condition is proposed to evaluate a diffusion-like process. Diffusion signatures with respect to i…
▽ More
The structure of the diffusion regions in antiparallel magnetic reconnection is investigated by means of a theory and a Vlasov simulation. The magnetic diffusion is considered as relaxation to the frozen-in state, which depends on a reference velocity field. A field-aligned component of the frozen-in condition is proposed to evaluate a diffusion-like process. Diffusion signatures with respect to ion and electron bulk flows indicate the ion and electron diffusion regions near the reconnection site. The electron diffusion region resembles the energy dissipation region. These results are favorable to a previous expectation that an electron-scale dissipation region is surrounded by an ion-scale Hall-physics region.
△ Less
Submitted 28 March, 2014; v1 submitted 14 March, 2014;
originally announced March 2014.
-
Dynamics and microinstabilities at perpendicular collisionless shock: A comparison of large-scale two-dimensional full particle simulations with different ion to electron mass ratio
Authors:
Takayuki Umeda,
Yoshitaka Kidani,
Shuichi Matsukiyo,
Ryo Yamazaki
Abstract:
Large-scale two-dimensional (2D) full particle-in-cell simulations are carried out for studying the relationship between the dynamics of a perpendicular shock and microinstabilities generated at the shock foot. The structure and dynamics of collisionless shocks are generally determined by Alfven Mach number and plasma beta, while microinstabilities at the shock foot are controlled by the ratio of…
▽ More
Large-scale two-dimensional (2D) full particle-in-cell simulations are carried out for studying the relationship between the dynamics of a perpendicular shock and microinstabilities generated at the shock foot. The structure and dynamics of collisionless shocks are generally determined by Alfven Mach number and plasma beta, while microinstabilities at the shock foot are controlled by the ratio of the upstream bulk velocity to the electron thermal velocity and the ratio of the plasma-to-cyclotron frequency. With a fixed Alfven Mach number and plasma beta, the ratio of the upstream bulk velocity to the electron thermal velocity is given as a function of the ion-to-electron mass ratio. The present 2D full PIC simulations with a relatively low Alfven Mach number (M_A ~ 6) show that the modified two-stream instability is dominant with higher ion-to-electron mass ratios. It is also confirmed that waves propagating downstream are more enhanced at the shock foot near the shock ramp as the mass ratio becomes higher. The result suggests that these waves play a role in the modification of the dynamics of collisionless shocks through the interaction with shock front ripples.
△ Less
Submitted 30 January, 2014; v1 submitted 23 January, 2014;
originally announced January 2014.
-
Microinstabilities at perpendicular collisionless shocks: A comparison of full particle simulations with different ion to electron mass ratio
Authors:
Takayuki Umeda,
Yoshitaka Kidani,
Shuichi Matsukiyo,
Ryo Yamazaki
Abstract:
A full particle simulation study is carried out for studying microinstabilities generated at the shock front of perpendicular collisionless shocks. The structure and dynamics of shock waves are determined by Alfven Mach number and plasma beta, while microinstabilities are controlled by the ratio of the upstream bulk velocity to the electron thermal velocity and the plasma-to-cyclotron frequency. T…
▽ More
A full particle simulation study is carried out for studying microinstabilities generated at the shock front of perpendicular collisionless shocks. The structure and dynamics of shock waves are determined by Alfven Mach number and plasma beta, while microinstabilities are controlled by the ratio of the upstream bulk velocity to the electron thermal velocity and the plasma-to-cyclotron frequency. Thus, growth rates of microinstabilities are changed by the ion-to-electron mass ratio, even with the same Mach number and plasma beta. The present two-dimensional simulations show that the electron cyclotron drift instability is dominant for a lower mass ratio, and electrostatic electron cyclotron harmonic waves are excited. For a higher mass ratio, the modified two-stream instability is dominant and oblique electromagnetic whistler waves are excited, which can affect the structure and dynamics of collisionless shocks by modifying shock magnetic fields.
△ Less
Submitted 10 April, 2012;
originally announced April 2012.
-
Relativistic Electron Shock Drift Acceleration in Low Mach Number Galaxy Cluster Shocks
Authors:
Shuichi Matsukiyo,
Yutaka Ohira,
Ryo Yamazaki,
Takayuki Umeda
Abstract:
An extreme case of electron shock drift acceleration in low Mach number collisionless shocks is investigated as a plausible mechanism of initial acceleration of relativistic electrons in large-scale shocks in galaxy clusters where upstream plasma temperature is of the order of 10 keV and a degree of magnetization is not too small. One-dimensional electromagnetic full particle simulations reveal th…
▽ More
An extreme case of electron shock drift acceleration in low Mach number collisionless shocks is investigated as a plausible mechanism of initial acceleration of relativistic electrons in large-scale shocks in galaxy clusters where upstream plasma temperature is of the order of 10 keV and a degree of magnetization is not too small. One-dimensional electromagnetic full particle simulations reveal that, even though a shock is rather moderate, a part of thermal incoming electrons are accelerated and reflected through relativistic shock drift acceleration and form a local nonthermal population just upstream of the shock. The accelerated electrons can self-generate local coherent waves and further be back-scattered toward the shock by those waves. This may be a scenario for the first stage of the electron shock acceleration occurring at the large-scale shocks in galaxy clusters such as CIZA J2242.8+5301 which has well defined radio relics.
△ Less
Submitted 31 August, 2011;
originally announced September 2011.
-
Cross-scale coupling at a perpendicular collisionless shock
Authors:
Takayuki Umeda,
Masahiro Yamao,
Ryo Yamazaki
Abstract:
A full particle simulation study is carried out on a perpendicular collisionless shock with a relatively low Alfven Mach number (M_A=5). In the present study, we have performed a two-dimensional (2D) electromagnetic full particle simulation with a "shock-rest-frame model". The simulation domain is taken to be larger than the ion inertial length in order to include full kinetics of both electrons a…
▽ More
A full particle simulation study is carried out on a perpendicular collisionless shock with a relatively low Alfven Mach number (M_A=5). In the present study, we have performed a two-dimensional (2D) electromagnetic full particle simulation with a "shock-rest-frame model". The simulation domain is taken to be larger than the ion inertial length in order to include full kinetics of both electrons and ions. The present simulation result has confirmed the transition of shock structures from the cyclic self-reformation to the quasi-stationary shock front. During the transition, electrons and ions are thermalized in the direction parallel to the shock magnetic field. Ions are thermalized by low-frequency electromagnetic waves (or rippled structures) excited by strong ion temperature anisotropy at the shock foot, while electrons are thermalized by high-frequency electromagnetic waves (or whistler mode waves) excited by electron temperature anisotropy at the shock overshoot. Ion acoustic waves are also excited at the shock overshoot where the electron parallel temperature becomes higher than the ion parallel temperature. We expect that ion acoustic waves are responsible for parallel diffusion of both electrons and ions, and that a cross-scale coupling between an ion-scale mesoscopic instability and an electron-scale microscopic instability is important for structures and dynamics of a collisionless perpendicular shock.
△ Less
Submitted 1 April, 2011;
originally announced April 2011.
-
Nonlinear evolution of parallel propagating Alfven waves: Vlasov - MHD simulation
Authors:
Y. Nariyuki,
T. Umeda,
T. Kumashiro,
T. Hada
Abstract:
Nonlinear evolution of circularly polarized Alfvén waves are discussed by using the recently developed Vlasov-MHD code, which is a generalized Landau-fluid model. The numerical results indicate that as far as the nonlinearity in the system is not so large, the Vlasov-MHD model can validly solve time evolution of the Alfvénic turbulence both in the linear and nonlinear stages. The present Vlasov-…
▽ More
Nonlinear evolution of circularly polarized Alfvén waves are discussed by using the recently developed Vlasov-MHD code, which is a generalized Landau-fluid model. The numerical results indicate that as far as the nonlinearity in the system is not so large, the Vlasov-MHD model can validly solve time evolution of the Alfvénic turbulence both in the linear and nonlinear stages. The present Vlasov-MHD model is proper to discuss the solar coronal heating and solar wind acceleration by Alfveń waves propagating from the photosphere.
△ Less
Submitted 5 October, 2009;
originally announced October 2009.
-
Electron acceleration at a low-Mach-number perpendicular collisionless shock
Authors:
Takayuki Umeda,
Masahiro Yamao,
Ryo Yamazaki
Abstract:
A full particle simulation study is carried out on the electron acceleration at a collisionless, relatively low Alfven Mach number (M_A=5), perpendicular shock. Recent self-consistent hybrid shock simulations have demonstrated that the shock front of perpendicular shocks has a dynamic rippled character along the shock surface of low-Mach-number perpendicular shocks. In this paper, the effect of…
▽ More
A full particle simulation study is carried out on the electron acceleration at a collisionless, relatively low Alfven Mach number (M_A=5), perpendicular shock. Recent self-consistent hybrid shock simulations have demonstrated that the shock front of perpendicular shocks has a dynamic rippled character along the shock surface of low-Mach-number perpendicular shocks. In this paper, the effect of the rippling of perpendicular shocks on the electron acceleration is examined by means of large-scale (ion-scale) two-dimensional full particle simulations. It has been shown that a large-amplitude electric field is excited at the shock front in association with the ion-scale rippling, and that reflected ions are accelerated upstream at a localized region where the shock-normal electric field of the rippled structure is polarized upstream. The current-driven instability caused by the highly-accelerated reflected ions has a high growth rate to large-amplitude electrostatic waves. Energetic electrons are then generated by the large-amplitude electrostatic waves via electron surfing acceleration at the leading edge of the shock transition region. The present result suggests that the electron surfing acceleration is also a common feature at low-Mach-number perpendicular collisionless shocks.
△ Less
Submitted 10 December, 2008;
originally announced December 2008.
-
Two-dimensional full particle simulation of a perpendicular collisionless shock with a shock-rest-frame model
Authors:
Takayuki Umeda,
Masahiro Yamao,
Ryo Yamazaki
Abstract:
A two-dimensional (2D) shock-rest-frame model for particle simulations is developed. Then full kinetic dynamics of a perpendicular collisionless shock is examined by means of a 2D full particle simulation. We found that in the 2D simulation there are fewer nonthermal electrons due to surfing acceleration which was seen in the previous 1D simulations of a high Mach number perpendicular shock in a…
▽ More
A two-dimensional (2D) shock-rest-frame model for particle simulations is developed. Then full kinetic dynamics of a perpendicular collisionless shock is examined by means of a 2D full particle simulation. We found that in the 2D simulation there are fewer nonthermal electrons due to surfing acceleration which was seen in the previous 1D simulations of a high Mach number perpendicular shock in a low-beta and weakly magnetized plasma. This is because the particle motion along the ambient magnetic field disturbs the formation of coherent electrostatic solitary structures which is necessary for electron surfing acceleration.
△ Less
Submitted 16 June, 2008; v1 submitted 1 June, 2008;
originally announced June 2008.
-
Full particle simulation of a perpendicular collisionless shock: A shock-rest-frame model
Authors:
Takayuki Umeda,
Ryo Yamazaki
Abstract:
The full kinetic dynamics of a perpendicular collisionless shock is studied by means of a one-dimensional electromagnetic full particle simulation. The present simulation domain is taken in the shock rest frame in contrast to the previous full particle simulations of shocks. Preliminary results show that the downstream state falls into a unique cyclic reformation state for a given set of upstrea…
▽ More
The full kinetic dynamics of a perpendicular collisionless shock is studied by means of a one-dimensional electromagnetic full particle simulation. The present simulation domain is taken in the shock rest frame in contrast to the previous full particle simulations of shocks. Preliminary results show that the downstream state falls into a unique cyclic reformation state for a given set of upstream parameters through the self-consistent kinetic processes.
△ Less
Submitted 29 October, 2006; v1 submitted 25 July, 2006;
originally announced July 2006.