A&A 468, 1075-1081 (2007)
DOI: 10.1051/0004-6361:20066977
J. I. Sakai - Y. Tanaka
Laboratory for Plasma Astrophysics, Faculty of Engineering, University of Toyama, 3190 Gofuku, Toyama, 930-8555, Japan
Received 20 December 2006 / Accepted 1 March 2007
Abstract
Aims. We investigate what happens when a fast magnetosonic shock wave associated with a coronal mass ejection (CME) collides obliquely with a coronal streamer with a stable current sheet.
Methods. A two-dimensional relativistic and fully electromagnetic Particle-In-Cell (PIC) code is used.
Results. It is shown that the fast magnetosonic shock with initial Alfvén Mach
compresses the current sheet, resulting in strong deformation of the current sheet. In the later stage a secondary fast magnetosonic shock wave can be generated almost perpendicular to the current sheet, and it propagates away to the opposite side of the original shock. This newly generated shock wave may emit a type II radio burst. The simulation results may be applied to a split of electromagnetic wave emissions when a shock wave associated with CMEs collides obliquely with a coronal streamer. For weak initial Alfvén Mach
,
a secondary shock wave does not appear, while the current sheet can be deformed and become unstable for tearing-like modes associated with magnetic reconnection.
Key words: plasmas - shock waves - radiation mechanisms: non-thermal - methods: numerical - Sun: radio radiation
Sheeley et al. (2000) show examples of the fast CMEs associated with the deflection of coronal streamers, and van der Holst et al. (2002) observe that when the fast CME collides with a helmet streamer, a split occurs in the type II emission. They modeled this CME-streamer interaction by using a 2D MHD simulation. Mancuso & Abbo (2004) investigate the observed bifurcation of the radio emission lanes by modeling the interaction of a piston-driven spherical shock with a vertical current sheet located above the active region.
Previous theoretical and simulation models of the interaction between a current sheet and a fast magnetosonic shock wave are based on MHD models. To understand the physical processes of the emissions of type II radio bursts during the shock-streamer collision, we must take plasma kinetic processes into account by using particle-in-cell (PIC) simulation. Recently, Haruki et al. (2006) investigated how the emission of electromagnetic waves can be enhanced when a fast magnetosonic shock wave associated with CMEs collides perpendicular to a coronal streamer with a stable current sheet, by using a two-dimensional relativistic and fully electromagnetic PIC code. It was shown that the ions in front of the shock can be accelerated by the shock-surfing acceleration mechanism (Sagdeev & Shapiro 1973; Lee et al. 1996; Lipatov & Zank 1999; Hoshino & Shimada 2002; Lee et al. 2005; Dieckmann et al. 2006). This shock compresses the current sheet, resulting in the local electron temperature anisotropy. Electron Bernstein waves are generated by the local electron temperature anisotropy and are converted into electromagnetic waves (X-mode) through the linear mode conversion due to density inhomogeneity (Murtaza & Shukla 1984; Yin et al. 1999). The nonlinear excitation process of the electron Bernstein waves in a pinching current sheet was investigated by Haruki & Sakai (2001). There, the electromagnetic waves are observed in both forward and backward regions of the shock. These simulation results were applied to the enhancement of electromagnetic wave emission when a shock wave associated with a CME collides with a coronal streamer.
In the present paper we extend the above work by Haruki et al. (2006) to investigate what happens when a fast magnetosonic
shock wave associated with a CME collides obliquely with a coronal streamer with a
stable current sheet. It is shown that the fast magnetosonic shock for strong initial Alfvén Mach
compresses the current sheet, resulting in
strong deformation of the current sheet.
In the later stage a secondary fast magnetosonic shock wave can be generated almost perpendicular to the current sheet, and it propagates to the
opposite side of the original shock. This newly generated shock may emit the type II radio burst, as shown by Sakai et al. (2005). For weak initial Alfvén Mach
a secondary shock wave does not appear, while the current sheet can be deformed and become unstable for tearing-like modes associated with magnetic reconnection.
This paper is organized as follows. In Sect. 2 we present our simulation model and in Sect. 3 we show our simulation results. In the final section we summarize our results.
![]() |
Figure 1: Schematic picture of the coronal mass ejection (CME) associated with a shock wave that collides obliquely with Harris type current sheet and the coordinate system used in the simulation. |
Open with DEXTER |
![]() |
Figure 2:
The initial magnetic field structure of a Harris type current sheet and proton density. The square high-density region moves
with Alfvén Mach speed
![]() |
Open with DEXTER |
Figure 1 shows a schematic picture of our simulation model. It is assumed that a CME originating in the
lower corona produces a shock and that this shock undergoes a oblique collision with a coronal streamer that
contains a current sheet as shown. The initial magnetic field and plasma number
density for a Harris type current sheet are the same as in Haruki et al. (2006). Figure 2 shows the initial magnetic field
and the proton density distribution. As shown in Fig. 2 we impose a dense plasma cloud (square region) far from the center of the
current sheet to excite a fast magnetosonic shock wave by imposing the plasma velocity with Alfvén Mach speed
on the cloud. The initial velocity of the cloud is
.
The code used here is a 2D relativistic and fully electromagnetic PIC code, modified from the 3D TRISTAN
code (Buneman 1993). The system size in 2D is
,
where
is the simulation grid size. It is assumed that the physical
quantities are constant in z, i.e.,
.
The average number density of background
electrons and ions is 100 per cell. The given particles fill the whole domain. Charge neutrality is initially kept because of having the
same position for the electron and proton. Periodic boundary conditions in the x-direction and free boundary conditions
in the y-direction are imposed on particles and fields.
First, we impose the following Harris type current sheet (Harris 1962; Pritchett 2001) that is characterized by the magnetic field, current,
and plasma particle number density, namely,
The other parameters used for this simulation are as follows. The simulation time step is
,
where
is used as a normalization of time. The Debye length is
.
The electron skin depth is
,
which is also used as a normalization of space. The mass
ratio of electron and proton is
,
whose value is artificial due to computer memory. The temperatures of the
electrons and protons are the same,
.
The Larmor radii of electrons and protons are
and
,
respectively.
We performed another simulation case without the dense plasma cloud. We confirmed that the
current sheet is very stable until
.
For this case, we did not observe any wave generation.
![]() |
Figure 3:
Time development of proton number density ![]() ![]() |
Open with DEXTER |
![]() |
Figure 4:
Time development of proton phase space,
![]() ![]() ![]() |
Open with DEXTER |
![]() |
Figure 5:
Time development of proton density distribution at a)
![]() ![]() ![]() ![]() |
Open with DEXTER |
![]() |
Figure 6:
a) Proton density distribution at
![]() ![]() ![]() ![]() ![]() |
Open with DEXTER |
In Fig. 6 we present the spacial structure associated with the secondary generated shock wave at
.
Figures 6a, c, and d show proton density distribution, magnetic field Bx, and electric field Ey near the the shock front, respectively. Figure 6b shows the proton phase space plot (
)
at
.
Some protons near Y=70 can be accelerated to the positive z-direction.
These protons are accelerated by the shock surfing acceleration mechanism (Sagdeev & Shapiro 1973) near the secondary generated shock wave. We can estimate the maximum proton velocity obtained near the shock front as
.
We obtain
by using the field data of Ey=0.14 and Bx=0.76 at X=36.6 and Y=72.9. This value is almost consistent with the simulation result seen in Fig. 6b. Therefore we may expect that the secondary generated fast magnetosonic shock wave can emit the type II radio bursts shown in Sakai et al. (2005). To see the wave emission in this simulation, we need larger simulation box.
![]() |
Figure 7:
Time development of proton density distribution for
![]() ![]() ![]() ![]() |
Open with DEXTER |
![]() |
Figure 8:
Current Jz for
![]() ![]() |
Open with DEXTER |
![]() |
Figure 9:
Time development of vector plots of magnetic fields for
![]() ![]() ![]() ![]() |
Open with DEXTER |
Finally we present simulation results with a different Mach number that is smaller than the above case.
We show the case where the initial driving flow velocity is
.
Figure 7 shows the time development of proton number density at a)
,
b)
,
and c)
.
As seen in Fig. 7b, the plasma cloud associated with the shock moves along the current sheet, resulting in the bending of the current sheet. In this case there appears no secondary shock wave as observed in the previous case. During the later stage (see Fig. 7c), the current sheet breaks into two current filaments. Figure 8 shows the spacial distribution of current Jz at
,
corresponding to Fig. 7c. In Fig. 9 we show the time development of vector plots of magnetic fields corresponding to Fig. 7. From Fig. 9c we see the formation of magnetic islands corresponding to the current filaments seen in Fig. 7c. From all this, we conclude that for a weak Mach number like
the current sheet can be deformed and become unstable for tearing-like modes associated with magnetic reconnection.
By using a 2D relativistic and fully electromagnetic PIC code, we investigated the oblique collision between a fast magnetosonic shock and a Harris type current sheet.
It was shown that the fast magnetosonic shock with Mach number
compresses the current sheet, resulting in
strong deformation of it.
In the later stage a secondary fast magnetosonic shock wave can be generated almost perpendicular to the current
sheet and it propagates away to the
opposite side of the original shock. This newly generated shock wave may emit the type II radio burst. The simulation results may
be applied to a split of electromagnetic wave emissions when a shock wave associated with CMEs collides
obliquely with a coronal streamer.
While for a weak shock with
,
the secondary shock wave does not appear. But the current sheet can be deformed and become unstable for tearing-like modes associated with magnetic reconnection.
Acknowledgements
The authors thank the anonymous referee for useful comments that improved our original paper.