
2
features of CTS in two plasmas. We theoretically cal-
culate the CTS spectra in the two-plasma state in Sec.
II. When two electron distribution functions overlap with
each other, spectral asymmetry arises due to the differ-
ent rates of electron Landau damping for the two peaks.
We perform numerical simulations to calculate the CTS
spectra in Sec. III. In the presence of the electron two-
stream-type instability, we observe a new spectrum corre-
sponding to the excited wave and the spectral asymmetry
of the ion-acoustic feature due to the change in the elec-
tron distribution function. We provide a discussion and
summary of our research in Sec. IV.
II. DYNAMIC STRUCTURE FACTOR IN A
TWO-PLASMA STATE
We theoretically calculate the CTS spectrum in a two-
plasma state. We consider the electron and ion distribu-
tion functions to be the sum of two Maxwell distributions:
fe,i(v) = X
j
nj
e,i
ne,i
fj
e,i(v),(1)
where fj
e,i(v) and nj
e,i are the distribution function and
number density of the j-th plasma species, respectively.
The total number density is given by ne,i =Pjnj
e,i. The
individual distribution functions are written as fj
e,i(v) =
exp[−(v−vdj )2/v2
te,i]/(π1/2vte,i), where vdj is the j-th
drift velocity. We assume there are no net currents, so
J=Pj[(qenj
e+qinj
i)vdj ] = 0, where qe,i is the charge
of an electron or ion. The thermal velocity is given by
vte,i = (2kBTe,i/me,i)1/2, where kB,Te,i, and me,i are
the Boltzmann constant, temperature, and mass, respec-
tively. In this paper, we assume the temperatures of the
two components (electrons and ions) are the same for
simplicity. We set the ion-to-electron mass ratio to 100.
The reduced mass is employed to compare the dynamic
structure factor with the simulated spectrum in Sec. III.
The dynamic structure factor with a realistic mass ratio is
presented in Appendix A. The dynamic structure factor
that describes the spectral shape of Thomson scattering
near the equilibrium plasmas is expressed as7
S(k, ω) = 2π
k1−χe
ε
2feω
k+Z
χe
ε
2fiω
k.
(2)
The first and second terms on the right-hand side in
Eq. (2) are associated with the non-collective scattering
or Langmuir (electron plasma) waves and ion-acoustic
waves, respectively. The susceptibility χe,i is given by
χe,i =4πq2
e,ine,i
me,ik2Z∞
−∞
∂fe,i (v)
∂v
ω
k−vdv
=−4πq2
e,i
me,ik2X
j"nj
e,i
v2
te,i
Z0ω
k−vdj
vte,i #,
(3)
where Z0(ζ) is the derivative of the plasma dispersion
function. The wavenumber and frequency of the ob-
served density fluctuation are given by k=kS−kIand
ω=ωS−ωI, where kS,kI,ωS, and ωIare the scattered
wavenumber, incident wavenumber, scattered frequency,
and incident frequency, respectively. The dielectric func-
tion satisfies ε= 1+χe+χi. The peaks of the CTS spec-
trum are mainly determined by ε∼0, which represents
the resonant condition of parametric resonance. The dis-
persion relations of the Doppler-shifted ion-acoustic and
electromagnetic waves are given by (ω−k·vd)2=c2
Sk2
and ω2=ω2
p+c2k2, where vd,cS, and ωpare the flow
velocity, sound velocity, and plasma frequency, respec-
tively. Considering the dispersion relations and the res-
onant condition, we find the peak wavenumbers of the
CTS spectra kS±:
k2
S±'c2+ (vd±cS)2
c2−(vd±cS)2kIsin θ
2+2(vd±cS)
c2−(vd±cS)2ωI2
+k2
Icos2θ
2,
(4)
where θis the angle between the incident and scattered
wavenumbers. The peaks are determined not only by ε∼
0 but also by the electron susceptibility in the numerator
in Eq. (2). Substituting the dielectric function into the
absolute value term in Eq. (2), we get
χe
ε
2=
1
1 + x
2
=1
(1 + Re[x])2+ (Im[x])2,(5)
where x= (1 + χi)/χe. The peak is determined to be
at the wavenumber with the smallest denominator in Eq.
(5) where Re[x]∼ −1 and Im[x]∼0 are satisfied.
We plot Eqs. (1), (2), and (5) in Fig. 1with the param-
eters Te/(mec2)=2×10−3,Te1=Te2= 10Ti1= 10Ti2,
n1
e,i/n2
e,i = 1, Z= 1, vd1= 0, ckI/ωpe = 5, and θ= 90°.
We compare the spectra obtained from f=f1+f2
(S1+2) and the sum of the spectra obtained indepen-
dently from f1(S1) and f2(S2) by changing the drift
velocity vd2/vte = 0, 1, and 3. Figures 1(a)–1(c) show
the scattered spectrum in terms of the wavenumber shift,
∆k=kS−kI, normalized to ωpe/c, where ωpe is the
electron plasma frequency. The black, green dashed, red
dotted, and blue dotted curves in Figs. 1(a)–1(c) show
the CTS spectra of S1+2,S1+S2,S1, and S2, respec-
tively. The red and blue curves in Figs. 1(d)–1(f) are
the real and imaginary parts of (1 + χi)/χein terms of
the wavenumber shift, respectively. The resonant condi-
tion of Re[(1 + χi)/χe]∼ −1 is shown by the horizontal
red dashed line, and that of Im[(1 + χi)/χe]∼0 is shown
as the horizontal blue dashed line. We plot Eq. (4) as
the vertical dotted lines in Figs. 1(a)–1(f) as a criterion
of the peaks determined by Eq. (5). The insets in Figs.
1(b), 1(c), 1(e), and 1(f) are enlarged versions of the left
spectrum for ease of view. Figures 1(g)–1(i) represent the
electron (blue) and ion (red) distribution functions. The