Hostname: page-component-7b9c58cd5d-dlb68 Total loading time: 0 Render date: 2025-03-15T12:25:34.700Z Has data issue: false hasContentIssue false

Nonlinear electromagnetic Eigen modes of a self created magnetized plasma channel and its stimulated Raman scattering

Published online by Cambridge University Press:  15 December 2011

Updesh Verma*
Affiliation:
Govt. Degree College Ballarpur, Rampur, U.P., India
A.K. Sharma
Affiliation:
Center for Energy Studies, Indian Institute of Technology Delhi, New Delhi, India
*
Address correspondence and reprint requests to: Updesh Verma, Center for Energy Studies, Indian Institute of Technology Delhi, New Delhi-110016, India. E-mail: updeshv@gmail.com
Rights & Permissions [Opens in a new window]

Abstract

A theoretical formalism is developed to obtain the mode structure of right circularly polarized nonlinear laser Eigen mode in a self created plasma channel in the presence of an axial magnetic field. The nonlinearity in electron response arises due to relativistic mass effect and ponderomotive force induced density redistribution. The Eigen mode is seen to be unstable to stimulated Raman backscattering involving an electrostatic quasi-mode and a scattered electromagnetic wave. The growth rate increases with ambient magnetic field.

Type
Research Article
Copyright
Copyright © Cambridge University Press 2011

1. INTRODUCTION

The propagation of high power laser of finite spot size in plasmas is influenced by diffraction and self focusing effects. The former causes divergence of the beam on the scale of a Rayleigh length R d = kr 02, where k is the wave number of the laser and r 0 is its spot size. The self focusing arises due to nonlinear permittivity that may be caused by a variety of sources, e.g., ohmic nonlinearity, ponderomotive nonlinearity and relativistic mass nonlinearity (Sodha et al., Reference Sodha, Khanna and Tripathi1974a, Reference Sodha, Mittal, Kumar and Tripathi1974b, Reference Sodha, Ghatak and Tripathi1976; Gill et al., Reference Gill, Mahajan and Kaur2010; Saini & Gill, Reference Saini and Gill2004; Singh & Tripathi, Reference Singh and Tripathi2009; Verma & Sharma, Reference Verma and Sharma2011). The theories of self focusing are mostly limited to paraxial ray approximation where eikonal of the wave is expanded up to the second power of r, the radial (transverse) coordinate. Sodha et al. (Reference Sodha, Khanna and Tripathi1974a, Reference Sodha, Mittal, Kumar and Tripathi1974b) have surveyed the early studies on self focusing and filamentation in unmagnetized plasmas with primary focus on long time scale Ohmic and ponderomotive nonlinearities. Esaray et al. (1997) have studied the short-pulse laser propagation in gases undergoing ionization. They also observed optical guiding of a laser beam in plasmas by taking relativistic effects, ponderomotive effects, and preformed density channels. Kurki-Suonio et al. (Reference Kurki-suonio, Morrison and Tajima1989) have studied the propagation of a circularly polarized laser in underdense plasma, including relativistic and ponderomotive nonlinearities, and showed that it remains intact as it propagates inside the plasma. Liu and Tripathi (Reference Liu and Tripathi2000) have developed a nonparaxial theory of laser self defocusing in a tunnel ionizing gas by expanding the eikonal up to fourth power in r.

Pathak and Tripathi (Reference Pathak and Tripathi2006) have developed a formalism to study the nonlinear Eigen modes of plasma without resorting to any expansion of the eikonal. The wave equation is solved numerically. The laser intensity profile in an unmagnetized plasma under relativistic and ponderomotive nonlinearity, turns out to be close to a Lorentzian. The Eigen modes in high density plasmas are stable to stimulated Raman scattering (SRS). However, in low density plasmas, laser Eigen modes are unstable to SRS. This study is limited to unmagnetized plasma.

In this paper, we study the nonlinear Eigen modes of a laser in a self created magnetized plasma channel (Verma & Sharma, Reference Verma and Sharma2009a; Panwar & Sharma, Reference Panwar and Sharma2009), under the combined effects of ponderomotive and relativistic mass nonlinearities, and examine their stability to stimulated Raman scattering. This study is relevant to intense short pulse laser produced plasmas where multi mega-gauss magnetic fields have been observed (Cai et al., Reference Cai, Yu, Zhu and Zhou2007; Borghesi et al, Reference Borghesi, MacKinnon, Barringer, Gaillard, Gizzi, Meyer, Willi, Pukhov and Meyer-ter-Vehn1997; Verma & Sharma, Reference Verma and Sharma2009b). Horovitz et al. (Reference Horovitz, Eliezer, Henis, Paiss, Moshe, Ludmirsky, Werdiger, Arad and Zigler1998) experimentally measured the generation of axial magnetic fields. They performed experiment with 1.06 µm laser of 7 ns pulse duration and 109 to1014 W/cm2 intensity range. They observed axial fields from 500 G to 2.17 MG. Naseri et al. (Reference Naseri, Bychenkov and Rozmus2010) reported that inverse Faraday Effect is the main mechanism of axial magnetic field generation in inhomogeneous plasmas. Pukhov and Meyer-ter-Vehn (Reference Pukhov and Meyer-ter-Vehn1998) in their three-dimensional particle-in-cell simulations of intense short pulse laser gas jet interaction have observed 100 MG magnetic fields in one-tenth critical density plasmas at ~1019 W/cm2 intensity at 0.8 µm wavelength. Frolov (Reference Frolov2009) observed quasi-static magnetic fields generated in low-density plasma in the interaction of counter propagating laser pulses. The magnetic field introduces anisotropy and cyclotron resonance effects in the propagation of electromagnetic waves and opens up a rich variety of electrostatic modes. Parametric coupling (Liu & Tripathi, Reference Liu and Tripathi1986; Porkolab & Chang, Reference Porkolab and Chang1972; Short & Simon Reference Short and Simon2004) between these modes has been a fascinating subject of extensive study for quiet some time. Clark and Fisch (Reference Clark and Fisch2005) proposed a scheme that utilizes the stimulated Raman backscattering (SRBS) in plasma to amplify a short and frequency downshifted seed pulse. Sharma and Kourakis (Reference Sharma and Kourakis2009) used SRBS to investigate the pulse compression and pulse amplification mechanisms. Trines et al. (Reference Trines, Kamp, Schep, Leemans, Esarey and Sluijter2004) did PIC simulation and show that suppression of SRBS increases the high energy electron yield. Sajal et al. (Reference Sajal and Tripathi2004, Reference Sajal, Panwar and Tripathi2006) have studied relativistic stimulated Raman forward scattering (SRFS) and shown that the Langmuir wave in the presence of a strong azimuthal magnetic field goes over to a localized lower hybrid wave, and the growth rate of the stimulated Raman backward scattering process is decreased. Liu and Tripathi (Reference Liu and Tripathi1996) have studied the stimulated forward Raman scattering of a laser in the channel for none relativistic case. The localization effect reduces the growth rate of Raman process. Liu and Tripathi (Reference Liu and Tripathi2001) also studied the relativistic laser guiding in a strong azimuthal magnetic field in plasma in the paraxial ray approximation. A linearly polarized laser with circular spot size is seen to acquire an elliptical spot size as it is self focused (Liu & Tripathi, Reference Liu and Tripathi1994).

In Section 2, we deduce the nonlinear wave equation governing the propagation of a right circularly polarized mode in plasma in the presence of axial magnetic field, allowing for a transverse amplitude variation. The nonlinearity is taken to arise due to relativistic mass effect and ponderomotive force induced electron density redistribution. The ions are taken to be stationary. The wave equation is solved numerically to obtain the nonlinear Eigen modes. In Section 3, we study the stimulated Raman backscattering of a large amplitude Eigen mode. In Section 4, we discuss the results.

2. NONLINEAR ELECTROMAGNETIC PLASMA EIGEN MODE IN PRESENCE OF AXIAL MAGNETIC FIELD

Consider pre-formed plasma of electron density n 0 embedded in a static magnetic field $B_s \hat z$. A circularly polarized electromagnetic mode propagates through it with t, z variation of the electric field as

(1)
\vec E_0=A_0 \lpar x\rpar \lpar \hat x+i\hat y\rpar e^{ - i\lpar {\rm \omega} _0 t - k_0 z\rpar }, \eqno\lpar 1\rpar

where A 0 has a slow dependence on x. The magnetic field of laser is

(2)
{\rm \omega} =i{\rm \omega} _0 g^{1/2} .\eqno\lpar 2\rpar

The laser imparts an oscillatory velocity to electrons,

(3)
{\rm \vec v}_0=\displaystyle{{e\vec E_0 } \over {mi{\rm \gamma} _0 \lpar {\rm \omega} _0 - {\rm \omega} _c /{\rm \gamma} _0 \rpar }}\comma \; \eqno\lpar 3\rpar

where ωc = eB s/m, −e and m are the electron charge and mass and ${\rm \gamma} _0=\displaystyle{1 \over {\lpar 1 - {\rm v}_0^2 /c^2 \rpar ^{1/2} }}$ is the gamma factor. For |v0|2/c 2 ≪ 1, simplifies to

(4)
{\rm \gamma} _0=1+\displaystyle{{a_0^2 } \over {2\lpar 1 - {\rm \omega} _c /{\rm \omega} _0 \rpar ^2 }} - \displaystyle{{\lpar 1+3{\rm \omega} _c /{\rm \omega} _0 \rpar } \over {8\lpar 1 - {\rm \omega} _c /{\rm \omega} _0 \rpar ^5 }}a_0^4\comma \; \eqno\lpar 4\rpar

where a 02 = e 2A 02/m 2ω02c 2.

The laser exerts a direct current ponderomotive force on electrons,

(5)
\vec F_p=- \displaystyle{{m{\rm \gamma} _0 } \over 2}\lpar {\rm \vec v}_{\rm 0}\cdot \nabla \rpar {\rm \vec v}_{\rm 0}^{\;\ast } - \displaystyle{e \over 2}\lpar {\rm \vec v}_{\rm 0} \times \vec B_0^{\ast}\rpar \comma \; \eqno\lpar 5\rpar
=\displaystyle{{ - e^2 \lsqb \vec E_0 \times \lpar \nabla \times \vec E_0 \rpar \rsqb } \over {2mi{\rm \gamma} _0 \lpar {\rm \omega} _0 - {\rm \omega} _c /{\rm \gamma} _0 \rpar i{\rm \omega} _0 }}+\displaystyle{{e^2 \lsqb \lpar \vec E_0 \cdot\nabla \rpar \vec E_0 \rpar \rsqb } \over {2m{\rm \gamma} _0 \lpar {\rm \omega} _0 - {\rm \omega} _c /{\rm \gamma} _0 \rpar ^2 }}.

Using identities $\vec A \times \vec B \times \vec C=\vec B\lpar \vec A.\vec C\rpar - \vec C\lpar \vec A.\vec B\rpar $ and $\nabla \lpar \vec A.\vec B\rpar =\vec A \times \lpar \nabla \times \vec B\rpar +\vec B \times \lpar \nabla \times \vec A\rpar +\lpar \vec A.\nabla \rpar \vec B+\lpar \vec B.\nabla \rpar \vec A$, we get $\vec F_p={{ {&- e^2 {\rm \omega} _0 \lsqb \lpar 1 - {\rm \omega} _c /2{\rm \omega} _0 {\rm \gamma} _0 \rpar \, \nabla \lpar \vec E_0 .\vec E_0 \rpar \cr & - \lpar {\rm \omega} _c /2{\rm \omega} _0 {\rm \gamma} _0 \rpar \lpar \vec E_0 \times \nabla \times \vec E_0 \rpar \rsqb} \over \displaystyle{{2m{\rm \gamma} _0 {\rm \omega} _0 \lpar {\rm \omega} _0 - {\rm \omega} _c /{\rm \gamma} _0 \rpar ^2 }}$. But we assumed ωc0γ0 < 1, hence one can neglect the second term from above equation and hence we get

\vec F_p=\displaystyle{{e^2 \lpar 1 - {\rm \omega} _c /2{\rm \omega} _0 {\rm \gamma} _0 \rpar \, \nabla \lpar \vec E_0 .\vec E_0 \rpar } \over {2m{\rm \gamma} _0 \lpar {\rm \omega} _0 - {\rm \omega} _c /{\rm \gamma} _0 \rpar ^2 }}=e\nabla {\rm \phi} _p\comma \;

hence

(6)
{\rm \phi} _p=\displaystyle{{ - e} \over {4m{\rm \gamma} _0 \lpar {\rm \omega} _0 - {\rm \omega} _c /{\rm \gamma} _0 \rpar ^2 }}\lpar 2 - {\rm \omega} _c /{\rm \omega} {\rm \gamma} _0 \rpar \lpar A_0 \cdot A_0^{\ast}\rpar \comma \; \eqno\lpar 6\rpar

The ponderomotive force pushes the electrons away from the laser axis, creating a low electron density channel near the laser axis. For laser pulse duration shorter than the ion response time, the ions may be taken to be stationary. The ion space charge field $\vec E_s=- \vec \nabla {\rm \phi} _s\comma \; $ pulls the electrons back. In a time on the order of ωp−1, where ωp2 = n 0e 2/mɛ0, a quasi-steady-state is realized, when the ponderomotive potential balances the space charge potential and further lowering of electron density at the axis stops. The time scale of plasma channel formation is of the order of electron plasma period (~1/ωp). Hence in the quasi-steady-state

(7)
{\rm \phi} _s=- {\rm \phi} _p=\displaystyle{e \over {4m{\rm \gamma} _0 \lpar {\rm \omega} _0 - {\rm \omega} _c /{\rm \gamma} _0 \rpar }}\lpar 2 - {\rm \omega} _c /{\rm \omega} _0 {\rm \gamma} _0 \rpar .\eqno\lpar 7\rpar

Using ϕs in Poisson's equation, ∇2ϕs = e(n e − n 0)/ɛ0, we obtained the modified electron density

(8)
n_e=n_0 - \displaystyle{{{\rm \varepsilon} _0 } \over e}\nabla ^2 {\rm \phi} _p .\eqno\lpar 8\rpar

With the relativistic modification of electron mass and ponderomotive modification of density, the non-zero component of the conductivity tensor of plasma can be written as

{\rm \sigma} _{xx}={\rm \sigma} _{yy}=\displaystyle{{in_e e^2 } \over {m{\rm \omega} _0 {\rm \gamma} _0 \lpar 1 - {\rm \omega} _c^2 /{\rm \omega} _0^2 {\rm \gamma} _0^2 \rpar}} \comma\; {\rm \sigma} _{xy}=- {\rm \sigma} _{yx}=\displaystyle{{n_e e^2 {\rm \omega} _c } \over {m{\rm \omega} _0^2 {\rm \gamma} _0^2 \lpar 1 - {\rm \omega} _c^2 /{\rm \omega} _0^2 {\rm \gamma} _0^2 \rpar}}\comma \;{\rm \sigma}_{zz}=\displaystyle{{in_e e^2 \over {n_0{\rm \gamma}_0 {\rm \varepsilon}_0m{\rm \omega} _0^2}}.
(9)
\hskip3.5pc{\rm \sigma} _{xz}={\rm \sigma} _{zx}={\rm \sigma} _{yz}={\rm \sigma} _{zy}=0. \fleqno\lpar 9\rpar

Using third and fourth Maxwell's equations, we deduce the wave equation as

(10)
\nabla ^2 \vec E_0 - \nabla \lpar \nabla .\vec E_0 \rpar +\displaystyle{{{\rm \omega} _0^2 } \over {c^2 }}\lpar \mathop {\rm \varepsilon} \limits_=.\vec E_0 \rpar =0.\eqno\lpar 10\rpar

Effective dielectric constant for right circularly polarized extraordinary mode in magnetized plasma can be written as

(11)
{\rm \varepsilon} _ \pm={\rm \varepsilon} _{xx} \pm {\rm \varepsilon} _{xy}=1 - \displaystyle{{n_e {\rm \omega} _p^2 } \over {n_0 {\rm \gamma} _0 {\rm \omega} _0^2 \lpar 1 \mp {\rm \omega} _c /{\rm \omega} {}_0{\rm \gamma} _0 \rpar }}\comma \; \eqno\lpar 11\rpar

where

{\rm \varepsilon} _{xx}={\rm \varepsilon} _{yy}=1 - \displaystyle{{n_e {\rm \omega} _p^2 } \over {n_0 {\rm \gamma} _0 {\rm \omega} _0^2 \lpar 1 - {\rm \omega} _c^2 /{\rm \omega} _0^2 {\rm \gamma} _0^2 \rpar }}\comma \;
{\rm \varepsilon} _{xy}=- {\rm \varepsilon} _{yx}=\displaystyle{{in_e {\rm \omega} _c {\rm \omega} _p^2 } \over {n_0 {\rm \gamma} _0^2 {\rm \omega} _0^3 \lpar 1 - {\rm \omega} _c^2 /{\rm \omega} _0^2 {\rm \gamma} _0^2 \rpar }}\comma \;\;\;\;
\hskip4.6pc{\rm \varepsilon} _{zz}=1 - \displaystyle{{n_e {\rm \omega} _p^2 } \over {n_0 {\rm \gamma} _0 {\rm \omega} _0^2 }}\fleqno
(12)
\hskip3.2pc\hbox{and} \,{\rm \varepsilon} _{xx}={\rm \varepsilon} _{yy}={\rm \varepsilon} _{yz}={\rm \varepsilon} _{zy}=0.\fleqno\lpar 12\rpar

Using $\partial E_z /\partial z=- \lpar 1/{\rm \varepsilon} _{zz} \rpar \lsqb \lpar \partial /\partial x\rpar \lpar {\rm \varepsilon} _{xx} E_x+{\rm \varepsilon} _{xy} E_y \rpar +\lpar \partial /\partial y\rpar \lpar \!-\! {\rm \varepsilon} _{xy} E_x+{\rm \varepsilon} _{xx} E_y \rpar \rsqb $, and following Sodha et al. (Reference Sodha, Khanna and Tripathi1974a, Reference Sodha, Mittal, Kumar and Tripathi1974b, Reference Sodha, Ghatak and Tripathi1976) the wave equation governing the propagation of right circularly polarized wave can be written as

(13)
\displaystyle{{\partial ^2 A_+} \over {\partial z^2 }}+\displaystyle{1 \over 2}\left({1+\displaystyle{{{\rm \varepsilon} _ - } \over {{\rm \varepsilon} _{zz} }}} \right)\nabla _ \bot ^2 A_++ \displaystyle{{{\rm \omega} _0^2 } \over {c^2 }}{\rm \varepsilon} _ - A_ +=0\comma \; \eqno\lpar 13\rpar

when ωc ≪ ω0, then $\left({1+\displaystyle{{{\rm \varepsilon} _ - } \over {{\rm \varepsilon} _{zz} }}} \right)\approx 2$ and hence the above equation reduces to well-known result as derived by Liu and Tripathi (Reference Liu and Tripathi1994). Here A + = E x + iE y is the extraordinary mode of the laser. Let A + = A +0 (x)e ik 0z, where A +0 = 2A 0. In an under dense plasma ωp2 ≪ ω02 and $\left({1+\displaystyle{{{\rm \varepsilon} _ - } \over {{\rm \varepsilon} _{zz} }}} \right)\approx 2$. Using these values, Eq (13) becomes

(14)
\displaystyle{{d^2 a_0 } \over {dx^2 }}=\left({k_0^2 - \displaystyle{{{\rm \omega} _0^2 } \over {c^2 }}{\rm \varepsilon} _ - } \right)a_0 .\eqno\lpar 14\rpar

Put the value of modified electron density from Eq. (8), ɛ_ can be written as

(15)
{\rm \varepsilon} _ -=1 - {{{{\rm \omega} _p^2 } \over {{\rm \omega} _0^2 \lpar {\rm \gamma} _0+{\rm \omega} _c /{\rm \omega} _0 \rpar }}} - \displaystyle{{\lpar 2 - {\rm \omega} _c /{\rm \omega} _0 {\rm \gamma} _0 \rpar c^2 \lpar \partial ^2 a_0^2 /\partial x^2 \rpar } \over {4{\rm \gamma} _0 {\rm \omega} _0^2 \lpar {\rm \gamma} _0+{\rm \omega} _c /{\rm \omega} _0 \rpar \lpar 1 - {\rm \omega} _c /{\rm \omega} _0 {\rm \gamma} _0 \rpar ^2 }}.\eqno\lpar 15\rpar

Put the value of ɛ_ in the above equation, it becomes

(16)
\eqalign{ \displaystyle{{d^2 a_0^2 } \over {dX^2 }} &=\displaystyle{1 \over {2a_0^2 \lpar 1 - 2{\rm \alpha} a_0^2 \rpar }}\left({\displaystyle{{da_0^2 } \over {dX}}} \right)^2\cr &\quad +\displaystyle{{2a_0^2 } \over {\lpar 1 - 2{\rm \alpha} a_0^2 \rpar }}\left[{\left({\displaystyle{{k_0 c} \over {{\rm \omega} _0 }}} \right)^2 - 1} \right]\cr &\quad +\displaystyle{{2\lpar {\rm \omega} _p^2 /{\rm \omega} _0^2 \rpar a_0^2 } \over {\lpar 1 - 2{\rm \alpha} a_0^2 \rpar \lpar {\rm \gamma} _0+{\rm \omega} _c /{\rm \omega} _0 \rpar }}} \eqno\lpar 16\rpar

where a 0 = eA 0/mω0c, X = xω0/c and

(17)
{\rm \alpha}=\displaystyle{{\lpar 2 - {\rm \omega} _c /{\rm \omega} _0 {\rm \gamma} _0 \rpar } \over {4{\rm \gamma} _0 \lpar {\rm \gamma} _0+{\rm \omega} _c /{\rm \omega} _0 \rpar \lpar 1 - {\rm \omega} _c /{\rm \omega} _0 {\rm \gamma} _0 \rpar ^2 }}.\eqno\lpar 17\rpar

We have solved Eq. (16) numerically to obtain the mode structure of the Eigen mode for parameters: ωp0 = 0.4, k 0c0 = 0.91652 and ωc0 = 0.4.

Figure 1 shows the mode structure of a laser in the presence of magnetic field. It resembles Gaussian profile with half width x 0 = 0.436c0, a 0 = a 00ψ(x), ${\rm \psi} \lpar x\rpar =\left({\displaystyle{1 \over {{\rm \pi} ^{1/4} }}} \right)e^{ - x^2 /2x_0^2 } $. This is different from the Pathak and Tripathi, Lorentzian profile that have obtained in unmagnetized plasma. Figure 2 shows the variation of ω0p with k 0cp for different values of magnetic fields. The value of ω0p is minimum at minimum k0cp and increases monotonically with increasing k 0cp. The minimum value of ω0p increases with increasing the value of ωcp.

Fig. 1. (Color online) Electromagnetic mode amplitude square as a function of normalized transverse distance for ωp202 = 0.25, a 00 = 0.4, and ωc0 = 0.4.

Fig. 2. (Color online) Variation of normalized frequency ω0p with normalized wave vector k 0cp for a high power laser in plasmas. The other parameters are: a 00 = 0.4, ωcp = 0.2 and 0.4.

3. STIMULATED RAMAN BACK SCATTERING

Now we examine the parametric decay of a large amplitude Eigen mode into a reactive quasi-mode of potential

(18)
{\rm \phi} _{\rm \omega}=A\lpar x\rpar \exp \lsqb \!- i\lpar {\rm \omega} t - kz\rpar \rsqb \comma \; \eqno\lpar 18\rpar

and a backscattered electromagnetic wave with electric and magnetic fields

(19)
\vec E_1=\vec A_1 \lpar \hat x - i\hat y\rpar e^{ - i\lpar {\rm \omega} _1 t - k_1 z\rpar } \comma\; \vec B_1=\displaystyle{{\vec k_1 \times \vec E_1 } \over {{\rm \omega} _1 }}\comma \; \eqno\lpar 19\rpar

where ω1 = ω−ω0, k 1 = k − k 0. In an underdense plasma, ω ≪ ω0 hence |ω1| ≈ ω0, |k 1| ≈ k 0 and k ≈ 2k 0. The scattered wave imparts oscillatory velocity to electrons

(20)
{\rm \vec v}_1=\displaystyle{{e\vec E_1 } \over {mi{\rm \gamma} _0 \lpar {\rm \omega} _1+{\rm \omega} _c /{\rm \gamma} _0 \rpar }}\comma\eqno\lpar 20\rpar

The pump and scattered waves exert a ponderomotive force on electrons at $\lpar {\rm \omega}\comma \; {\vec k}\rpar $, $\vec F_{\,pz\comma {\rm \omega} }=e\nabla {\rm \phi} _{\,p{\rm \omega} }=ik{\rm \phi} _{\,p{\rm \omega} } \hat Z$, where

(21)
{\rm \phi} _{\,p{\rm \omega} }=\displaystyle{e \over {m{\rm \gamma} _0\, k}}\left[{\displaystyle{{k_1 {\rm \omega} _0 \lpar {\rm \omega} _1+{\rm \omega} _c /{\rm \gamma} _0 \rpar +k_0 {\rm \omega} _1 \lpar {\rm \omega} _0 - {\rm \omega} _c /{\rm \gamma} _0 \rpar } \over {{\rm \omega} _0 {\rm \omega} _1 \lpar {\rm \omega} _0 - {\rm \omega} _c /{\rm \gamma} _0 \rpar \lpar {\rm \omega} _1+{\rm \omega} _c /{\rm \gamma} _0 \rpar }}} \right]A_0 A_1 .\eqno\lpar 21\rpar

The electron response to ϕω and ϕpω can be obtained by solving the equation of motion,

m{\rm \gamma} _0 \displaystyle{{\partial {\rm \vec v}} \over {\partial t}}=- e\vec E_{\rm \omega}+\vec F_{\rm \omega} - e\lpar {\rm \vec v} \times \vec B_s \rpar
(22)
\vec{\rm v}=\displaystyle{{ - e\nabla \lpar {\rm \phi} _{\rm \omega}+{\rm \phi} _{\,pz} \rpar } \over {mi{\rm \omega} {\rm \gamma} _0 }}\comma \; \eqno\lpar 22\rpar

using this value in the continuity equation we get the modified density of electrons

n_{\rm \omega}=\displaystyle{1 \over {i{\rm \omega} }}n_o \lpar \nabla .{\rm \vec v}\rpar .

Using Poisson's equation ∇2ϕω = en ω0, and defining eϕω/mc 2 = φω, we get

(23)
\eqalign{& \displaystyle{{\partial ^2 {\rm \varphi} _{\rm \omega} } \over {\partial x^2 }}\left[{1 - \displaystyle{{{\rm \omega} _p^2 } \over {{\rm \gamma} _0 \lpar {\rm \omega} ^2 - {\rm \omega} _c^2 /{\rm \gamma} _0^2 \rpar }}} \right]- k_z^2 \lpar 1 - {\rm \omega} _p^2 /{\rm \omega} ^2 {\rm \gamma} _0 \rpar {\rm \varphi} _{\rm \omega}\cr &\quad= - \, \, \displaystyle{{{\rm \omega} _p^2\, k_z {\rm \omega} _0^2 a_0 a_1 } \over {{\rm \omega} ^2 {\rm \gamma} _0^2 }}\cr &\qquad\times\left[{\displaystyle{{k_1 {\rm \omega} _0 \lpar {\rm \omega} _1+{\rm \omega} _c /{\rm \gamma} _0 \rpar +k_0 {\rm \omega} _1 \lpar {\rm \omega} _0 - {\rm \omega} _c /{\rm \gamma} _0 \rpar } \over {{\rm \omega} _0 {\rm \omega} _1 \lpar {\rm \omega} _1+{\rm \omega} _c /{\rm \gamma} _0 \rpar \lpar {\rm \omega} _0 - {\rm \omega} _c /{\rm \gamma} _0 \rpar }}} \right]\comma} \eqno\lpar 23\rpar

Where a 1 = eA 1/mω0c. For ω > ωc, ωp, Eq. (23) can be written as

(24)
\eqalign{\displaystyle{{\partial ^2 {\rm \varphi} _{\rm \omega} } \over {\partial x^2 }} - k_z^2 {\rm \varphi} _{\rm \omega}&=- \displaystyle{{{\rm \omega} _p^2\, k_z {\rm \omega} _0^2 a_0 a_1 } \over {{\rm \omega} ^2 {\rm \gamma} _0^2 \lpar 1 - {\rm \omega} _p^2 /{\rm \omega} ^2 {\rm \gamma} _0 \rpar }}\cr&\quad \times \left[{\displaystyle{{k_1 {\rm \omega} _0 \lpar {\rm \omega} _1+{\rm \omega} _c /{\rm \gamma} _0 \rpar +k_0 {\rm \omega} _1 \lpar {\rm \omega} _0 - {\rm \omega} _c /{\rm \gamma} _0 \rpar } \over {{\rm \omega} _0 {\rm \omega} _1 \lpar {\rm \omega} _1+{\rm \omega} _c /{\rm \gamma} _0 \rpar \lpar {\rm \omega} _0 - {\rm \omega} _c /{\rm \gamma} _0 \rpar }}} \right]\comma \; \eqno\lpar 24\rpar

the current density at the side band $\lpar {\rm \omega} _1\comma \; {\vec k}_1 \rpar $ can be written as

(25)
\vec J_1=- n_e e{\rm \vec v}{}_1 - \displaystyle{1 \over 2}n_{\rm \omega} e{\rm \vec v}_0^{\,\ast}.\eqno\lpar 25\rpar

The normalized wave equation for the side band, on assuming $\nabla .\vec E_1 \approx 0$, can be written as

(26)
\eqalign{& \displaystyle{{d^2 a_1 } \over {dx^2 }}+\Bigg[ \[\displaystyle{{{\rm \omega} _1^2 } \over {c^2 }} - k_1^2 - \displaystyle{{{\rm \omega} _p^2 } \over {c^2 {\rm \gamma} _0 \lpar 1+{\rm \omega} _c /{\rm \omega} _1 {\rm \gamma} _0 \rpar }} \cr &\quad - \displaystyle{{\lpar 2 - {\rm \omega} _c /{\rm \omega} _0 {\rm \gamma} _0 \rpar \lpar \partial ^2 a_0^2 /\partial x^2 \rpar } \over {4{\rm \gamma} _0^2 \lpar 1+{\rm \omega} _c /{\rm \omega} _0 {\rm \gamma} _0 \rpar \lpar 1 - {\rm \omega} _c /{\rm \omega} _0 {\rm \gamma} _0 \rpar ^2 }} \Bigg ] a_1 \cr &\quad =- \displaystyle{{n_{\rm \omega} e^2 {\rm \omega} _1 a_0^{\ast}} \over {2{\rm \varepsilon} _0 mc^2 {\rm \gamma} _0 \lpar {\rm \omega} _0 - {\rm \omega} _c /{\rm \gamma} _0 \rpar }}.} \eqno\lpar 26\rpar

The left-hand side of this equation is the same as that of Eq. (14) with ω0 replaced by $\left\vert {{\rm \omega} _1 } \right\vert $. Hence, in the absence of the parametric coupling the mode structure of a 1 is the same as that of the pump Eigen mode,

(27)
a_1=a_{10} {\rm \psi} \lpar x\rpar .\eqno\lpar 27\rpar

In low density plasma when the pump amplitude is large, the growth rate of plasma wave may become larger than the plasma frequency, ω2≫ωp2 turning into a reactive quasi-mode. Substituting for a 1 in Eq. (26), multiplying the resulting equation by ψdx′, x′ = x/x 0 , and integrating over x′ from −∞ to +∞ we obtain

(28)
\eqalign{& {\rm \omega} _1^2 - k_1^2 c^2=\displaystyle{{c^2 } \over {x_0^2 }}\left[1 - \mathop{\scale190%\vint}_\limits{ - \infty }_^\infty {x^{\prime 2} } {\rm \psi} ^2 dx^{\prime} \right]\cr &\quad - \displaystyle{{{\rm \omega} _p^2 {\rm \omega} _1\, k^2 c^2 \left\vert {a_{00} } \right\vert ^2 } \over {2{\rm \omega} ^2 {\rm \omega} _0 }}\mathop{\scale190%\vint}\limits_{ - \infty }^\infty {\displaystyle{{{\rm \psi} ^4 dx^{\prime}} \over {{\rm \gamma} _0^3 \lpar 1 - {\rm \omega} _c /{\rm \omega} _0 {\rm \gamma} _0 \rpar ^2 }}} \cr &\quad +\mathop{\scale190%\vint}_{ - \infty }^\infty \left[\displaystyle{{{\rm \omega} _p^2 } \over {{\rm \gamma} _0 \lpar 1+{\rm \omega} _c /{\rm \omega} _1 {\rm \gamma} _0 \rpar }}\right. \cr &\quad \left. +\displaystyle{{\lpar 2 - {\rm \omega} _c /{\rm \omega} _0 {\rm \gamma} _0 \rpar c^2 } \over {4{\rm \gamma} _0^2 \lpar 1+{\rm \omega} _c /{\rm \omega} _1 {\rm \gamma} _0 \rpar \lpar 1 - {\rm \omega} _c /{\rm \omega} _1 {\rm \gamma} _0 \rpar ^2 }}\displaystyle{{\partial ^2 } \over {\partial x^2 }}\left\vert {a_0 } \right\vert ^2 \right]{\rm \psi} ^2 dx^{\prime}.}\eqno\lpar 28\rpar

One may approximate the left-hand side as 2ωω0. Since ω3 < ω2, we can neglect ω3 terms. Hence Eq. (28) can be written as

(29)
{\rm \omega} = i {\rm \omega} g^{1/2}\comma \eqno \lpar 29\rpar

where

(30)
g=\displaystyle{{\lpar \Omega _p^2 {\rm \eta} _1^2 /4\rpar \mathop{\scale150%\vint}\limits_{ - \infty }^\infty {\displaystyle{{a_{00}^2 {\rm \psi} ^4 dx^{\prime}} \over {{\rm \gamma} _0^3 \lpar 1 - \Omega _c /{\rm \gamma} _0 \rpar ^2 }}} } \over {\left({\displaystyle{{0.278} \over {{\rm \eta} _0^2 }}} \right)+\mathop{\scale150%\vint}\limits_{ - \infty }^\infty {\displaystyle{{\Omega _p^2 {\rm \psi} ^2 dx^{\prime}} \over {2{\rm \gamma} _0 \lpar 1 - \Omega _c /{\rm \gamma} _0 \rpar }}}+\mathop{\scale150%\vint}\limits_{ - \infty }^\infty {\displaystyle{{\lpar 2 - \Omega _c /{\rm \gamma} _0 \rpar \displaystyle{{\partial ^2 a_0^2 } \over {\partial X^2 }}{\rm \psi} ^2 dx^{\prime}} \over {8{\rm \gamma} _0^2 \lpar 1 - \Omega _c /{\rm \gamma} _0 \rpar ^3 }}} }}\comma \; \eqno\lpar 30\rpar

other normalized quantities are as follows, Ωp = ωp0, Ωc = ωc0, Ω = ω/ω0, ${\rm \eta} _0=\displaystyle{{x_0 {\rm \omega} _0 } \over c}$, ${\rm \eta} _1=\displaystyle{{kc} \over {{\rm \omega} _0 }}$, x′ = x/x 0 and $X=\displaystyle{{x{\rm \omega} _0 } \over c}$. Solving Eq. (29) numerically, one obtains the growth rate Γ = Imω.

Figure 3 shows the schematic diagram of the laser plasma interaction. Figure 4 shows the variation of normalized growth rate as a function of normalized axial laser amplitude a 00 for ωc0 = 0.4 and ωp0 = 0.1, 0.2 and 0.4. As ω0 is a constant, different value of ωp0 imply different value of background electron density. Therefore, Figure 4 shows the importance of background electron density and showed that as the plasma frequency (density of electrons) increases, Γ/ω0 increases slowly. For larger values of laser intensity, the SRBS instability grows more rapidly. This shows that SRBS is significantly influenced by the plasma frequency. Figure 5 shows the variation of growth rate with axial laser amplitude. As a 00 increases from 0 to 0.4, Γ/ω0 varies from 0.001 to 0.024 for ωc0 = 0.1. But as ωc0 increases from 0.1 to 0.4, the value of Γ/ω0 increases from 0.024 to 0.035. This shows that magnetic field plays an important role in the growth of SRBS, and instability become stronger with increasing the magnetic field.

Fig. 3. (Color online) Schematic diagram of laser plasma interaction.

Fig. 4. (Color online) Normalized growth rate Γ/ω0 as a function of normalized axial laser amplitude a 00 for ωc0 = 0.4 and ωp0 = 0.1, 0.2 and 0.4.

Fig. 5. (Color online) Normalized growth rate as a function of normalized axial laser amplitude for a 00 for ωp0 = 0.4 and ωc0 = 0.1,0.2 and 0.4.

4. DISCUSSION

The presence of magnetic field modifies the Eigen frequency of the fundamental laser Eigen modes. The mode structure resembles a Gaussian distribution. This is different from unmagnetized plasma where Pathak and Tripathi obtain the mode structure resembling a Lorentizian profile. The growth rate of stimulated Raman backscattering in the reactive quasi-mode regime, relevant at high intensity, increases with the magnetic field. This is caused due to the cyclotron resonance effect though the theory is not applicable close to cyclotron resonance. It also increases with plasma density. As the laser intensity increases the width of nonlinear Eigen modes decreases. This is due to the fact that at higher laser intensity nonlinear self focusing effect becomes stronger, hence to compensate for it stronger diffraction effects, caused by reduced spot size, are required. Our theory is limited to a 0 ≤ 0.5 i.e., for weakly relativistic cases.

ACKNOWLEDGMENTS

The authors are grateful to Prof. V. K. Tripathi, IIT Delhi for fruitful discussions. The author is very thankful to the Govt. Degree College Bilaspur, Rampur, U.P. for encouragement and financial support.

References

REFERENCES

Borghesi, M., MacKinnon, A.J., Barringer, L., Gaillard, R., Gizzi, L.A., Meyer, C., Willi, O., Pukhov, A. & Meyer-ter-Vehn, J. (1997). Relativistic channeling of a picosecond laser pulse in a near-critical preformed plasma. Phys. Rev. Lett. 78, 879882.CrossRefGoogle Scholar
Cai, H-B., Yu, W., Zhu, S. & Zhou, C. (2007). Generation of strong quasistatic magnetic fields in interactions of ultraintense and short laser pulses with overdense plasma targets. Phys. Rev. E 76, 036403/7.Google ScholarPubMed
Clark, D.S. & Fisch, N.J. (2005). Raman laser amplification in preformed and ionizing plasmas. Laser Part. Beams 23, 101106.Google Scholar
Esarey, E., Sprangle, P., Krall, J. & Ting, A. (1997). Self-focusing and guiding of short laser pulses in ionizing gases and plasmas. IEEE J. Quan. Electron. 33, 18791914.Google Scholar
Frolov, A.A. (2009). Generation of quasistatic magnetic fields in the interaction of counter propagating laser pulses in a low-density plasma. Plasma Phys. Rept. 35, 668676.Google Scholar
Gill, T.S., Mahajan, R. & Kaur, R. (2010). Relativistic and ponderomotive effects on evolution of laser beam in a non-uniform plasma channel. Laser Part. Beams 28, 1120.Google Scholar
Horovitz, Y., Eliezer, S., Henis, Z., Paiss, Y., Moshe, E., Ludmirsky, A., Werdiger, M., Arad, B. & Zigler, A. (1998). The inverse Faraday effect in plasmas produced by circularly polarized laser light in the range of intensities 109–1014 W/cm2. Phys. Lett. A 246, 329334.CrossRefGoogle Scholar
Kurki-suonio, T., Morrison, P.J. & Tajima, T. (1989). Self-focusing of an electromagnetic beam in a plasma. Phys. Rev. A 40, 32303239.Google Scholar
Liu, C.S. & Tripathi, V.K. (1994). Interaction of Electromagnetic Waves with Electron Beams and Plasmas. Singapore: World Scientific.Google Scholar
Liu, C.S. & Tripathi, V.K. (1996). Stimulated Raman scattering in a plasma channel. Phys. Plasmas 3, 34103413.CrossRefGoogle Scholar
Liu, C.S. & Tripathi, V.K. (2001). Relativistic laser guiding in an azimuthal magnetic field in a plasma. Phys. Plasmas 8, 285288.Google Scholar
Liu, C.S. & Tripathi, V.K. (1986 ). Parametric instabilities in a magnetized plasma. Phys. Rept. 130, 143216.Google Scholar
Liu, C.S. & Tripathi, V.K. (2000). Laser frequency upshift, self-defocusing, and ring formation in tunnel ionizing gases and plasmas. Phys. Plasmas 7, 43604363.CrossRefGoogle Scholar
Naseri, N., Bychenkov, V.Yu. & Rozmus, W. (2010). Axial magnetic field generation by intense circularly polarized laser pulses in underdense plasmas. Phys. Plasmas 17, 083109/10.CrossRefGoogle Scholar
Panwar, A. & Sharma, A.K. (2009). Self-phase modulation of a laser in self created plasma channel . Laser Part. Beams 27, 249253.CrossRefGoogle Scholar
Pathak, V. B. & Tripathi, V.K. (2006). Nonlinear electromagnetic plasma Eigen modes and their stability to stimulated Raman scattering. Phys. Plasmas 13, 082105/4.Google Scholar
Porkolab, M. & Chang, R.P H. (1972). Instabilities and induced scattering due to nonlinear landau damping of longitudinal plasma waves in a magnetic field. Phys. Fluids 15, 283296.Google Scholar
Pukhov, A. & Meyer-ter-Vehn, J. (1998). Relativistic laser-plasma interaction by multi-dimensional particle-in-cell simulations. Phys. Plasmas 5, 18801886.Google Scholar
Saini, N.S. & Gill, T.S. (2004). Enhanced Raman scattering of a rippled laser beam in a magnetized collisional plasma. Laser Part. Beams 22, 3540.CrossRefGoogle Scholar
Sajal, V. & Tripathi, V.K. (2004). Stimulated Raman scattering of a laser beam in a plasma with azimuthal magnetic field. Phys. Plasmas 11, 42064212.Google Scholar
Sajal, V., Panwar, A. & Tripathi, V.K. (2006). Relativistic forward stimulated Raman scattering of a laser in a plasma channel. Phys. Scr. 74, 484488.CrossRefGoogle Scholar
Sharma, A. & Kourakis, I. (2009). Laser pulse compression and amplification via Raman backscattering in plasma. Laser Part. Beams 27, 279285.Google Scholar
Short, R.W. & Simon, A. (2004). Theory of three-wave parametric instabilities in inhomogeneous plasmas revisited. Phys. Plasmas 11, 53355340.Google Scholar
Singh, R. & Tripathi, V.K. (2009). Filamentation of laser in a magnetized plasma under relativistic and pondromotive nonlinearities. Phys. Plasmas 16, 052108/5.CrossRefGoogle Scholar
Sodha, M., Ghatak, A. & Tripathi, V.K. (1976). Self focusing of laser beams in plasmas and semiconductors. Prog. Opt. 13, 169265.CrossRefGoogle Scholar
Sodha, M., Khanna, R. & Tripathi, V.K. (1974 a). The self-focusing of electromagnetic beams in a strongly ionized magnetoplasma. J. Phys. D: Appl. Phys. 7, 21882197.CrossRefGoogle Scholar
Sodha, M., Mittal, R., Kumar, S. & Tripathi, V.K. (1974 b). Self-focusing of electromagnetic waves in a magneto plasma. Opto-Electron. 6, 167180.CrossRefGoogle Scholar
Trines, R.M.G.M., Kamp, L.P.J., Schep, T.J., Leemans, W.P., Esarey, E.H. & Sluijter, F.W. (2004). Enhancement of high-energy electron generation through suppression of Raman backscattering. Europhys. Lett. 66, 492498.CrossRefGoogle Scholar
Verma, U. & Sharma, A.K. (2009 a). Effect of self focusing on the prolongation of laser produced plasma channel. Laser Part. Beams 27, 3339.CrossRefGoogle Scholar
Verma, U. & Sharma, A.K. (2009 b). Laser second harmonic generation in a rippled density plasma in the presence of azimuthal magnetic field. Laser Part. Beams 27, 719724.CrossRefGoogle Scholar
Verma, U. & Sharma, A.K. (2011). Laser focusing and multiple ionization of Ar in a hydrogen plasma channel created by a pre-pulse. Laser Part. Beams 29, 219225.Google Scholar
Figure 0

Fig. 1. (Color online) Electromagnetic mode amplitude square as a function of normalized transverse distance for ωp202 = 0.25, a00 = 0.4, and ωc0 = 0.4.

Figure 1

Fig. 2. (Color online) Variation of normalized frequency ω0p with normalized wave vector k0cp for a high power laser in plasmas. The other parameters are: a00 = 0.4, ωcp = 0.2 and 0.4.

Figure 2

Fig. 3. (Color online) Schematic diagram of laser plasma interaction.

Figure 3

Fig. 4. (Color online) Normalized growth rate Γ/ω0 as a function of normalized axial laser amplitude a00 for ωc0 = 0.4 and ωp0 = 0.1, 0.2 and 0.4.

Figure 4

Fig. 5. (Color online) Normalized growth rate as a function of normalized axial laser amplitude for a00 for ωp0 = 0.4 and ωc0 = 0.1,0.2 and 0.4.