1 Introduction
With the development of laser technology, especially the invention of the chirped-pulse amplification (CPA) technique[ Reference Strickland and Mourou 1 ], laser peak power and focused intensity have been increasing dramatically[ Reference Danson, Haefner, Bromage, Butcher, Chanteloup, Chowdhury, Galvanauskas, Gizzi, Hein, Hillier, Hopps, Kato, Khazanov, Kodama, Korn, Li, Li, Limpert, Ma, Nam, Neely, Papadopoulos, Penman, Qian, Rocca, Shaykin, Siders, Spindloe, Szatmári, Trines, Zhu, Zhu and Zuegel 2 , Reference Radier, Chalus, Charbonneau, Thambirajah, Deschamps, David, Barbe, Etter, Matras, Ricaud, Leroux, Richard, Lureau, Baleanu, Banici, Gradinariu, Caldararu, Capiteanu, Naziru, Diaconescu, Iancu, Dabu, Ursescu, Dancus, Ur, Tanaka and Zamfir 3 ]. Simultaneously, intense laser–matter interactions have attracted great interest since they not only promise wide potential applications but also give birth to rich physical phenomena[ Reference Mourou, Tajima and Bulanov 4 , Reference Tanaka, Spohr, Balabanski, Balascuta, Capponi, Cernaianu, Cuciuc, Cucoanes, Dancus, Dhal, Diaconescu, Doria, Ghenuche, Ghita, Kisyov, Nastasa, Ong, Rotaru, Sangwan, Sóderstróm, Stutman, Suliman, Tesileanu, Tudor, Tsoneva, Ur, Ursescu and Zamfir 5 ]. With the increase of laser power, however, the generation and manipulation of intense laser pulses have become more and more challenging due to the limited damage threshold of traditional optical devices[ Reference Stuart, Feit, Herman, Rubenchik, Shore and Perry 6 ]. In contrast, plasma-based optical elements can sustain ultra-high light intensity and hence provide an alternative approach to the generation and manipulation of intense lasers.
In the last decades, plasma optics has developed into an attractive methodology for the manipulation and further amplification of intense lasers. For example, plasma mirrors have been widely used to enhance the contrast of intense laser pulses[ Reference Thaury, Quéré, Geindre, Levy, Ceccotti, Monot, Bougeard, Réau, d’Oliveira, Audebert, Marjoribanks and Martin 7 ] or tightly focus laser pulses toward extreme intensities[ Reference Nakatsutsumi, Kon, Buffechoux, Audebert, Fuchs and Kodama 8 ]. More significantly, both stimulated Raman and Brillouin scatterings in plasmas could be employed to amplify intense laser pulses[ Reference Malkin, Shvets and Fisch 9 – Reference Edwards, Mikhailova and Fisch 12 ]. Moreover, magnetized plasmas may have some specific advantages in the control and amplification of intense laser pulses[ Reference Weng, Zhao, Sheng, Yu, Luan, Chen, Yu, Murakami, Mori and Zhang 13 – Reference Li, Wu, Zuo, Zeng, Wang, Wang, Mu, Hu and Su 16 ]. Further, holographic plasma lenses have recently been proposed as a novel plasma optical device for broad applications[ Reference Dodin and Fisch 17 – Reference Edwards, Munirov, Singh, Fasano, Kur, Lemos, Mikhailova, Wurtele and Michel 19 ]. In particular, laser-induced plasma density gratings (PDGs) have been extensively studied and proposed for wide potential applications, including the temporal compression, polarization control and manipulation of intense laser pulses[ Reference Lehmann and Spatschek 20 – Reference Lehmann and Spatschek 25 ]. Analogously, laser-induced plasma density modulations can tune laser power distribution to achieve a symmetric inertial confinement fusion (ICF) implosion in the so-called crossed-beam energy transfer[ Reference Glenzer, MacGowan, Michel, Meezan, Suter, Dixit, Kline, Kyrala, Bradley, Callahan, Dewald, Divol, Dzenitis, Edwards, Hamza, Haynam, Hinkel, Kalantar, Kilkenny, Landen, Lindl, LePape, Moody, Nikroo, Parham, Schneider, Town, Wegner, Widmann, Whitman, Young, Van Wonterghem, Atherton and Moses 26 – Reference Kruer, Wilks, Afeyan and Kirkwood 28 ].
In this paper, we investigate the spatial and temporal distributions of the polarization state of a probe laser pulse after it passes through a PDG. As shown in Figure 1, the PDG is driven by two intersecting laser pulses that have Gaussian intensity profiles in both the longitudinal and transverse directions. Due to the birefringence caused by the PDG, it could be employed as a waveplate for an ultrashort laser pulse. In this study, however, the formed PDG is not only nonuniform in space but also varies over time. Therefore, the refractive index modulation induced by such a PDG will be nonuniform and time-dependent. As a result, the probe laser pulse may be depolarized by such a dynamic PDG, that is, its polarization state becomes variable in both space and time. The laser beams with varying polarization states may play an important role in mitigating parametric instabilities in laser–plasma interactions relevant to ICF[ Reference Barth and Fisch 29 , Reference Ma, Li, Weng, Yew, Kawata, Gibbon, Sheng and Zhang 30 ]. The paper is organized as follows. The formation and evolution of the PDG are simulated and analyzed in Section 2. The laser depolarization by the dynamic PDG is illustrated in Section 3. Finally, the conclusion and some discussions are presented in Section 4.
Schematic of laser depolarization by a PDG. The PDG driven by intersecting laser pulses #1 and #2 will be nonuniform in the
$y$
direction and also time-dependent. After the probe laser pulse passes through such a PDG, its polarization state will become nonuniform and time-dependent.
2 Dynamics of the plasma density grating
To understand the depolarization of laser beams, we first investigate the dynamics of the PDG, which can be well described by the fluid model in the linear stage in 1D cases[
Reference Sheng, Zhang and Umstadter
31
,
Reference Ma, Weng, Li, Li, Wang, Yew, Chen, McKenna and Sheng
32
]. For simplicity, two driver laser beams in the theoretical analysis are assumed to have constant intensities, and they propagate oppositely along the positive or negative
$y$
-axis into a homogeneous plasma. Using the normalized laser vector potentials
${\boldsymbol{a}}_{1,2}=e{\boldsymbol{A}}_{1,2}/{mc}^2$
, the total vector potential of the driver laser beams becomes
$$\begin{align}\boldsymbol{a}={a}_1\cos \left({k}_1y-{\omega}_1t\right)\;{\boldsymbol{e}}_x+{a}_2\cos \left({k}_2y+{\omega}_2t\right)\;{\boldsymbol{e}}_x.\end{align}$$
Here, the two driver laser beams have the same frequency
${\omega}_{1,2}={\omega}_0$
, and the same wave number in the plasma
${{k}_{1,2}=\sqrt{1-{n}_0/{n}_{\mathrm{c}}}{\omega}_0/c}$
, where
${n}_0$
is the initial plasma density,
${n}_{\mathrm{c}}={\omega}_0^2{\varepsilon}_0{m}_{\mathrm{e}}/{e}^2$
is the critical plasma density corresponding to the laser frequency
${\omega}_0$
,
${\varepsilon}_0$
is the permittivity of free space and
${m}_{\mathrm{e}}$
and
$e$
are the electron mass and charge, respectively. For a linearly polarized laser, its normalized vector potential
$a$
is related to its intensity
$I$
as
$I$
[W/cm2] =
$1.37\times {10}^{18}\cdot {a}^2/{\lambda}^2$
[μm2], where
$\lambda$
is the laser wavelength in a vacuum. The intersection of two driver laser pulses will result in a periodic modulation in the total light amplitude and thus induce the following periodic ponderomotive force:
$$\begin{align} \begin{array}{l}{\boldsymbol{F}}_{\mathrm{p}}={m}_{\mathrm{e}}{c}^2{a}_1{a}_2{k}_1\sin \left(2{k}_1y\right)\;{\boldsymbol{e}}_y.\end{array}\end{align}$$
Driven by such a pondermotive force, the electron density will be modified firstly. Then, the ions will follow the electrons due to the charge-separation field. Finally, a periodic plasma density structure that is termed as the PDG will be induced. In the linear growth stage, the density perturbation of the PDG can be estimated as follows:
$$\begin{align}\delta {n}_{\mathrm{e}}=-\left(2{k}^2{c}^2/{\omega}_{\mathrm{p}}^2\right){a}_1{a}_2\cos \left(2{k}_1y\right)\left[1-\cos \left({\omega}_{\mathrm{p}}t\right)\right],\end{align}$$
where
${\omega}_{\mathrm{p}}$
is the plasma frequency. After the linear growth stage, however, the saturation, collapse and regrowth of the PDG will take place[
Reference Ma, Weng, Li, Li, Wang, Yew, Chen, McKenna and Sheng
32
,
Reference Peng, Riconda, Grech, Su and Weber
33
]. Further, the formed PDG will be nonuniform, since the driver laser beams always have finite durations and spot sizes in reality.
To illuminate the complicated dynamics of the PDG, a 2D particle-in-cell (PIC) simulation is conducted using the Osiris code[
Reference Fonseca, Silva, Tsung, Decyk, Lu, Ren, Mori, Deng, Lee, Katsouleas and Adam
34
]. The simulation box is
$\left[-200\lambda \le x\le 200\lambda \right]\times \left[-160\lambda \le y\le 160\lambda \right]$
, which is divided into
$6400\times 5120$
cells. The plasma is located in
$\left[-200\lambda \le x\le 200\lambda \right]\times \left[-150\lambda \le y\le 150\lambda \right]$
and it has a uniform electron number density
${n}_{\mathrm{e}}=0.3{n}_{\mathrm{c}}$
in the region
$-195\lambda \le x\le 195\lambda$
and a
$5\lambda$
rising or falling ramp at each side in the
$x$
direction. For the plasma region, 16 macro-particles per cell are allocated. The plasma is assumed initially to be cold and the ions are the protons with
${m}_{\mathrm{p}}=1836{m}_{\mathrm{e}}$
. Two Gaussian driver laser pulses are incident along the positive or negative
$y$
-axis, respectively. They are p-polarized (the electric fields along the
$x$
direction) with a wavelength
$\lambda$
= 1 μm. These two driver pulses are assumed to have the same normalized amplitude
${a}_1={a}_2=0.02$
, the same spot size (full width at half maximum (FWHM) in intensity) of 100 μm and the same FWHM duration of
$150{T}_0$
. It is assumed that the peaks of the driver laser pulses will arrive at the simulation box center at
$t=0$
if there is no plasma, and the simulation begins at
$t=-460{T}_0$
.
The PDG structure at
$t=590{T}_0$
is displayed in Figure 2(a), where an obvious modulation in the plasma density is evidenced for the laser intersecting region. Zooming into the center region
$\left[-100\lambda \le x\le 100\lambda \right]\times \left[-5\lambda \le y\le 5\lambda \right]$
, Figure 2(b) illustrates that the density modulation has a clear periodicity along the
$y$
direction with a wavelength close to the laser wavelength in the plasma. However, Figure 2(c) shows that the density modulation in the
$x$
direction will be nonuniform and confined in a finite region, since the driver laser beams have finite spot sizes. Here, the peak and trough of the PDG are roughly achieved at
$y=\mp \lambda /4$
, respectively. Furthermore, Figure 2(d) shows that the density modulation in the
$y$
direction is also nonuniform.
The electron density distribution of (a) the overall plasma region and (b) the center region
$\left[-100\lambda \le x\le 100\lambda \right]\times \left[-5\lambda \le y\le 5\lambda \right]$
at
$t=590{T}_0$
, respectively. (c) The corresponding electron density profiles along the x direction at
$y=-\lambda /4,0,\lambda /4$
, respectively. Here,
$y=\pm \lambda /4$
are roughly along the plasma density trough and peak, respectively. (d) The corresponding electron density profiles along the
$y$
direction at
$x=0$
, in which the inset displays the enlarged density profile in the region
$-3\lambda <y<3\lambda$
. The upper and lower envelopes of this density profile are also outlined by the red and blue curves, respectively.
More importantly, the plasma density modulation will be time-dependent. The time evolution of the plasma density profile along the
$y$
-axis is displayed in Figure 3. It is illustrated that the peak density of the PDG increases gradually to a maximum value in the first development stage. Then the PDG will saturate and collapse due to the ion wave breaking[
Reference Ma, Weng, Li, Li, Wang, Yew, Chen, McKenna and Sheng
32
,
Reference Peng, Riconda, Grech, Su and Weber
33
]. Interestingly, the formation of the PDG will restart again after its collapse. As shown in Figure 3, the formation, saturation and collapse of the PDG could take place for several rounds. Due to its spatial nonuniformity and temporal evolution, the PDG will induce a complicated modulation upon the laser pulse passing through it. The laser depolarization by such a complicated and changeable PDG will be illuminated in the following paragraphs.
The time evolution of the electron density profile along the
$y$
-axis. Note that the PDG experiences a time periodic process of formation, saturation and collapse. The simulation parameters are the same as those in Figure 2.
3 Depolarization of the laser pulse
The propagation of an electromagnetic wave through the PDG is similar to the electron movement in crystalline solids. Therefore, the theory of energy bands in solid physics can be employed for the theoretical analysis of the laser propagation through the PDG. For simplification, the PDG is assumed to be composed of alternating high- and low-density layers. Each high-density layer has a uniform density
${n}_{\mathrm{h}}$
and a thickness
${l}_{\mathrm{h}}$
, and each low-density layer has a uniform density
${n}_{\mathrm{l}}$
and a thickness
${l}_{\mathrm{l}}$
. For keeping the plasma quasi-neutral, the relationship
${n}_0l={n}_{\mathrm{h}}{l}_{\mathrm{h}}+{n}_{\mathrm{l}}{l}_{\mathrm{l}}$
is satisfied, where
$l={l}_{\mathrm{h}}+{l}_{\mathrm{l}}$
is the periodic length of the PDG. For a spatially infinite PDG, the dispersion relation for the s-polarized (transverse electric (TE) mode) wave is given by the following:
$$\begin{align} \begin{array}{l}\cos (kl)=\cos \left({k}_{\mathrm{h}}{l}_{\mathrm{h}}\right)\cos \left({k}_{\mathrm{l}}{l}_{\mathrm{l}}\right)\\[4pt] {}\kern4em -\dfrac{1}{2}\left(\dfrac{k_{\mathrm{l}}}{k_{\mathrm{h}}}+\dfrac{k_{\mathrm{h}}}{k_{\mathrm{l}}}\right)\sin \left({k}_{\mathrm{h}}{l}_{\mathrm{h}}\right)\sin \left({k}_{\mathrm{l}}{l}_{\mathrm{l}}\right),\end{array}\end{align}$$
while for the p-polarized (transverse magnetic (TM) mode) wave it is given by the following:
$$\begin{align} \begin{array}{l}\cos (kl)=\cos \left({k}_{\mathrm{h}}{l}_{\mathrm{h}}\right)\cos \left({k}_{\mathrm{l}}{l}_{\mathrm{l}}\right)\\[4pt] {}\kern4em -\dfrac{1}{2}\left(\dfrac{n_{\mathrm{l}}^2{k}_{\mathrm{l}}}{n_{\mathrm{h}}^2{k}_{\mathrm{h}}}+\dfrac{n_{\mathrm{h}}^2{k}_{\mathrm{h}}}{n_{\mathrm{l}}^2{k}_{\mathrm{l}}}\right)\sin \left({k}_{\mathrm{h}}{l}_{\mathrm{h}}\right)\sin \left({k}_{\mathrm{l}}{l}_{\mathrm{l}}\right),\end{array}\end{align}$$
where
$k$
is the laser wave number in the PDG when it propagates along the layers,
${k}_{\mathrm{h}}=\left(\omega /c\right)\sqrt{1-{n}_{\mathrm{h}}/{n}_{\mathrm{c}}}$
and
${k}_{\mathrm{l}}=\left(\omega /c\right)\sqrt{1-{n}_{\mathrm{l}}/{n}_{\mathrm{c}}}$
.
Due to their different dispersion relations, the s- and p-polarized electromagnetic waves will have different phase velocities. More importantly, the phase velocity difference will vary in both time and space as long as the PDG is spatially nonuniform and temporally variable. Using the electron density profile presented in Figure 2(d), the phase velocity of the s- and p-polarized light waves along the
$y$
-axis can be obtained according to Equations (4) and (5), respectively. The results are displayed in Figure 4, in which the phase velocity difference is clearly illustrated and it is nonuniform along the
$y$
-axis. Consequently, the final phase difference between the s- and p-polarized light waves will be proportional to the integral of the phase velocity difference over time, as follows:
$$\begin{align} \begin{array}{l}\Delta \phi ={\int}_0^lk\frac{v_{{\mathrm{TM}}}(x)-{v}_{{\mathrm{TE}}}(x)}{v_{{\mathrm{TE}}}(x)} \mathrm{d}l.\end{array}\end{align}$$
The phase velocities of the s-polarized (
${v}_{{\mathrm{TE}}}$
) and p-polarized (
${v}_{{\mathrm{TE}}}$
) light waves obtained from Equations (4) and (5), respectively, in which the electron density profile presented in Figure 2(d) is employed.
Since the phase velocity difference varies in time and space, the resultant phase difference between the s- and p-polarized light components of a probe laser pulse will also change with time and space. Therefore, the final polarization state of the probe laser pulse that initially has both non-zero s- and p-polarized light components will vary in both space and time. That is to say, the probe laser pulse will be depolarized by the dynamic PDG.
To illuminate the laser depolarization by the dynamic PDG, some additional PIC simulations are conducted. In the first simulation, the parameters of two driver laser pulses and the initial plasma are the same as those used in Figure 3. In addition to two driver laser pulses, a probe laser pulse whose polarization plane is oriented at 45° with respect to the
$x{-}y$
plane is incident along the
$x$
-axis into the PDG. The probe laser pulse has Gaussian intensity profiles in both the longitudinal and transverse directions, it has the peak intensity
${a}_0=0.05$
, the wavelength
$\lambda$
= 1 μm, the focal FWHM spot size of 100 μm and the FWHM duration
$\tau =200{T}_0$
. The peak of the probe laser pulse is assumed to arrive at
$\left(x,y\right)=\left(0,0\right)$
at
$t=515{T}_0$
if there is no plasma, that is, it is delayed by
$515{T}_0$
in time with respect to the driver laser pulses.
To describe the laser polarization state, the Stokes parameters
$I$
,
$Q$
,
$U$
and
$V$
are introduced[
Reference Tinbergen
35
]. Here
${I={\left|{E}_y\right|}^2+{\left|{E}_z\right|}^2}$
denotes the intensity regardless of polarization,
${Q={\left|{E}_y\right|}^2-{\left|{E}_z\right|}^2}$
denotes the linear polarization along the
$y\left(+\right)$
- or
$z\left(-\right)$
-axis,
$U=2\operatorname{Re}\left\{{E}_y^{\ast }{E}_z\right\}$
denotes the linear polarization at
$+45{}^{\circ}$
(+) or −45°
$\left(-\right)$
from the
$y$
-axis and
$V=2\operatorname{Im}\left\{{E}_y^{\ast }{E}_z\right\}$
denotes the right-handed
$\left(+\right)$
or left-handed
$\left(-\right)$
circular polarization, where
${E}_y$
(
${E}_z$
) is the complex amplitude of the electric field along the
$y(z)$
-axis and
${E}_y^{\ast }$
is the complex conjugate of
${E}_y$
. It is worth pointing out that the probe laser pulse initially has
$U=I$
and
$Q=V=0$
, since it is linearly polarized at +45° with respect to the
$y$
-axis. This implies that the s- and p-polarized light components of the probe laser are in phase and of the same amplitude at the beginning.
During the probe laser propagation in the PDG, its s- and p-polarized light components gradually become out of phase due to their different phase velocities. More importantly, the phase difference between them varies in time and space. As a result, the Stokes parameter distributions of the probe laser pulse become very complicated after it passes through the PDG, as shown in Figure 5. Above all, the parameters
$Q$
and
$V$
are no longer zero. Moreover, there are periodic variations in the parameters
$U$
and
$V$
. Further, the longitudinal (at
$y=0$
) and transverse (at
$x=425\lambda$
) profiles of the Stokes parameters are displayed in Figures 6(a) and 6(b), respectively. All of these illustrate that the polarization state of the probe laser becomes variable in both time and space after it passes through the PDG.
The spatial distributions of the Stoke parameters (a)
$I$
, (b)
$Q$
, (c)
$U$
and (d)
$V$
of the probe laser pulse at
$t=940{T}_0$
after it passes through the PDG. Here, all Stokes parameters are normalized to the instantaneous maximum laser intensity
${I}_{{\mathrm{max}}}$
. The simulation parameters are given in the text.
(a) Longitudinal profiles of the Stokes parameters at
$y=0$
and (b) transverse profiles of the Stokes parameters at
$x=425\lambda$
. (c) Longitudinally averaged polarization degree
${P}_{\mathrm{l}}$
and (d) transversely averaged polarization degree
${P}_{\mathrm{t}}$
. The simulation parameters are the same as those in Figure 5.
For the Stokes parameters, they always satisfy
${Q}^2+{U}^2+ {V}^2\leqslant {I}^2$
, where the equality is achieved for fully polarized lights. For partially polarized lights,
$0<\left({Q}^2+{U}^2+{V}^2\right)/{I}^2<1$
. It is worth pointing out that
${Q}^2+{U}^2+{V}^2={I}^2$
is satisfied at every infinitely small point in Figure 5. If the parameters
$I$
,
$Q$
,
$U$
and
$V$
are averaged over the longitudinal or transverse directions, however, the following averaged polarization degree
$$\begin{align*}{P}_{{\mathrm{t}},{\mathrm{l}}}=\sqrt{{\left\langle Q\right\rangle}_{{\mathrm{l}},{\mathrm{t}}}^2+{\left\langle U\right\rangle}_{{\mathrm{l}},{\mathrm{t}}}^2+{\left\langle V\right\rangle}_{{\mathrm{l}},{\mathrm{t}}}^2}/{\left\langle I\right\rangle}_{{\mathrm{l}},{\mathrm{t}}}\end{align*}$$
will be obviously smaller than one, where
${\left\langle \dots \right\rangle}_{{\mathrm{l}},{\mathrm{t}}}$
represents the average over the longitudinal (l) or transverse (t) directions. Figure 6(c) shows that the longitudinally averaged polarization degree
${P}_{\mathrm{l}}$
is lower than 0.5 in the near-axis region
$\mid y\mid \leqslant 50\lambda$
, while Figure 6(d) shows that the transversely averaged polarization degree
${P}_{\mathrm{t}}$
is less than 0.4 in the whole pulse duration. These indicate that the temporal and spatial variations of the polarization state induced by the PDG are equivalent to a decrease of the averaged polarization degree. Moreover, Figure 6(c) indicates that the longitudinally averaged polarization degree
${P}_{\mathrm{l}}$
is nonuniform in the transverse direction since the depth of the PDG is nonuniform. Further, the transverse distribution of the PDG depth varies from time to time, which strongly depends on the durations of the driver laser pulses as well as the time delay between the driver laser pulses and the probe laser pulse.
In the above analysis, the PDG is induced by two counter-propagating driver laser pulses, and the probe laser is incident along the direction orthogonal to the driver lasers. Actually, the PDG can also be induced by two intersecting driver laser pulses with an intersection angle
$\varphi$
, as shown in Figure 7. In this case, the interference term of two intersecting driver laser pulses can be written as follows:
$$\begin{align} {\mathit{\boldsymbol a}}_1\cdot {\mathit{\boldsymbol a}}_2&=\dfrac{a_1{a}_2\cos \varphi }{2}\Big[\cos \left(2\omega t-2 kx\cos \dfrac{\varphi }{2}\right)\nonumber\\[8pt] &\quad +\cos \left(2 ky\sin \dfrac{\varphi }{2}\right)\Big], \end{align}$$
Laser depolarization by the PDG that is induced by two intersecting laser pulses with an intersection angle
$\varphi$
.
where the second term on the right-hand side induces the formation of the PDG. According to Equation (7), the periodic length of the PDG is related to the angle as
${l=\lambda /\left[2\sin \left(\varphi /2\right)\right]}$
, except that the PDG generated in this case is similar to that generated by two counter-propagating lasers, as shown in Figures 2 and 3. The formed PDG in this case can also be employed to depolarize the probe laser pulse. In addition, this scenario may spontaneously take place in the beam-crossing region in either indirect- or direct-drive ICF schemes if the polarizations and incident angles of the laser beams are arranged properly[
Reference Glenzer, MacGowan, Michel, Meezan, Suter, Dixit, Kline, Kyrala, Bradley, Callahan, Dewald, Divol, Dzenitis, Edwards, Hamza, Haynam, Hinkel, Kalantar, Kilkenny, Landen, Lindl, LePape, Moody, Nikroo, Parham, Schneider, Town, Wegner, Widmann, Whitman, Young, Van Wonterghem, Atherton and Moses
26
,
Reference Michel, Divol, Williams, Weber, Thomas, Callahan, Haan, Salmonson, Dixit, Hinkel, Edwards, Macgowan, Lindl, Glenzer and Suter
27
].
To investigate the laser depolarization in the beam-crossing region, an additional PIC simulation is performed. In the simulation, a plasma with a uniform density
${n}_0=0.27{n}_{\mathrm{c}}$
is located in
$\left[-75\lambda \le x\le 75\lambda \right]\times \left[-400\lambda \le y\le 400\lambda \right]$
. The simulation box located within
$\left[-125\lambda \le x\le 425\lambda \right]\times \left[-400\lambda \le y\le 400\lambda \right]$
is resolved with 16 cells per wavelength, and 16 macro-particles are initially allocated per cell. Two driver laser pulses #1 and #2 are incident into the plasma at angles of
$\varphi /2=\pm \pi /4$
with respect to the
$x$
-axis, respectively. Two driver laser pulses are p-polarized (the electric fields are oriented in the
$x{-}y$
plane). A probe laser pulse whose polarization plane is oriented at
${45}^{\circ }$
with respect to the
$x{-}y$
plane is incident along the
$x$
-axis into the plasma. All laser pulses have spatio-temporal Gaussian intensity profiles, the same wavelength of
$\lambda$
= 1 μm and the same normalized amplitude
${a}_0={a}_1={a}_2=0.0468$
. Two driver laser pulses are assumed to have the same FWHM spot size of
$200\lambda$
and the same FWHM duration of
$200{T}_0$
. The probe laser pulse is assumed to have the FWHM spot size of
$100\lambda$
and the FWHM duration of
$150{T}_0$
. It is assumed that if there is no plasma the peaks of the driver laser pulses and the probe laser pulse will arrive at
$\left(x,y\right)=\left(0,0\right)$
at
$t=0$
and
$t=370{T}_0$
, respectively. The simulation begins at
$t=-305{T}_0$
.
Analogous to Figure 5, Figure 8 shows that the Stokes parameter distributions of the probe laser pulse become very complicated after it passes through the PDG in this case. The periodic variations in parameters U and V are clearly illustrated in Figures 8(c) and 8(d), respectively. Such periodic variations in the Stokes parameters are also shown by the longitudinal (at
$y=0$
) and transverse (at
$x=295\lambda$
) profiles of the Stokes parameters in Figures 9(a) and 9(b), respectively. Moreover, Figure 9(c) shows that the longitudinally averaged polarization degree
${P}_{\mathrm{l}}$
is lower than 0.2 in the region
$\mid y\mid \le 100\lambda$
, whereas Figure 9(d) shows that the transversely averaged polarization degree
${P}_{\mathrm{t}}$
is about 0.6 in the whole pulse duration. This demonstrates that the PDG can be induced not only by two counter-propagating laser pulses but also generally by two intersecting laser pulses with an intersection angle
$\varphi$
. The latter may be spontaneously encountered for the crossed laser beams in ICF plasmas. The formed PDG will in turn modulate the polarization state of the laser beams, and consequently reduce their polarization degree.
The spatial distributions of the Stoke parameters (a) I, (b) Q, (c) U and (d) V of the probe laser pulse at
$t=650{T}_0$
after it passes through the PDG that is induced by two intersecting laser pulses with an intersection angle
$\varphi$
. Here, all Stokes parameters are normalized to the instantaneous maximum laser intensity
${I}_{{\mathrm{max}}}$
. The simulation parameters are given in the text.
(a) Longitudinal profiles of the Stokes parameters at y = 0 and (b) transverse profiles of the Stokes parameters at
$x=295\lambda$
. (c) Longitudinally averaged polarization degree
${P}_{\mathrm{l}}$
and (d) transversely averaged polarization degree
${P}_{\mathrm{t}}$
. The simulation parameters are the same as those in Figure 8.
4 Discussion and conclusions
Since the normalized parameters are used in the simulations, such a dynamic PDG can in principle be used to depolarize the intense laser pulses at any wavelength. However, the density of the employed background plasma should increase for a shorter laser wavelength, which is often used in ICF. Correspondingly, collisional absorption may become the dominant process in a denser plasma. Fortunately, laser energy deposition via the collisional absorption usually is of benefit to ICF[ Reference Atzeni and Meyer-ter-Vehn 36 ]. In addition, we find that the plasma temperature will increase gradually. Correspondingly, the peak density of the PDG will decrease due to the enhanced thermal pressure[ Reference Ma, Weng, Li, Li, Wang, Yew, Chen, McKenna and Sheng 32 , Reference Peng, Riconda, Grech, Su and Weber 33 ]. This may finally set a limitation on the use of dynamic PDGs in hot plasmas. On the other hand, it is worth noting that the theoretical growth rates of parametric instabilities will be reduced due to the enhanced Landau damping in plasmas with high temperatures[ Reference Liu 37 ].
For the application of dynamic PDGs, the time-scale in which the polarization of the probe laser pulse becomes altered is a crucial factor. As shown in Figure 6(a), the dynamic PDG can make the polarization of the probe laser pulse vary obviously within a few tens of laser wave periods (
${\sim}0.2$
ps), while the stimulated Raman scattering (SRS) growth time in a homogenous plasma is also about 0.2 ps for the laser intensities of the order of
${10}^{14}$
W/cm2 that are usually employed in ICF[
Reference Montgomery
38
]. In general, the growth time of the stimulated Brillouin scattering (SBS) is much longer than that of the SRS. Therefore, one can expect that the laser depolarization via the dynamic PDG will play an important role in mitigating parametric instabilities such as SRS and SBS in ICF.
The time-scale in which the probe laser polarization is varied strongly depends on the typical variation time of the PDG. Further, it would decrease with the increase of the depth of the PDG (or the maximal achievable peak ion density). Moreover, it would also decrease with the increase of the width of the PDG, that is, the span of the PDG in the x-direction in the simulations. The latter is normally determined by the spot size of the driver laser pulses. However, it is difficult to give a formula for the time-scale in which the probe laser polarization is varied.
Concerning the variation time of the PDG, we have studied its dependence on the laser and plasma parameters by a series of 1D PIC simulations. As shown in Figure 10(a), we find that the saturation time
${T}_{\mathrm{s}}$
of the PDG decreases obviously with the increasing laser intensity
${a}_0$
for a given initial plasma density
${n}_0=0.3{n}_{\mathrm{c}}$
. Here, the saturation time
${T}_{\mathrm{s}}$
is defined as the time required to achieve the maximal peak ion density. However, the maximal achievable peak ion density
${n}_{{\mathrm{max}}}$
is nearly unchanged with the increasing
${a}_0$
. On the other hand, Figure 10(b) shows that the saturation time
${T}_{\mathrm{s}}$
of the PDG increases slightly with the increasing plasma density
${n}_0$
for a given laser intensity
${a}_0=0.02$
, while the maximal achievable peak ion density
${n}_{{\mathrm{max}}}$
increases obviously with the increasing plasma density
${n}_0$
.
The saturation time
${T}_{\mathrm{s}}$
(black solid lines) and the maximal achievable ion density
${n}_{{\mathrm{max}}}$
(red solid lines) as functions of (a) the laser intensity
${a}_0$
for a given initial plasma density
${n}_0=0.3{n}_{\mathrm{c}}$
and (b) the initial plasma density
${n}_0$
for a given laser intensity
${a}_0=0.02$
. Except for the laser intensities and initial plasma densities, other laser–plasma parameters are the same as those used in Figure 5.
Based on the above analysis, one can see that the intensities and spot sizes of the driver laser pulses and the initial plasma density will combine to set a minimum time-scale in which the probe laser polarization can be varied. In general, the time-scale in which the probe laser polarization is varied is expected to decrease with the increase of the initial background plasma density as well as the increase of the intensity and spot size of the driver laser pulses.
In summary, we have shown that the PDG induced by two intersecting driver laser pulses is not only nonuniform in space but also varies over time. Such a dynamic PDG can induce a modulation upon the polarization state of the probe laser pulse passing through it. The time-scale in which the probe laser polarization is varied depends on the initial background plasma density, as well as the intensities and spot sizes of the driver laser pulses. The polarization degree of the probe laser that is averaged over either the longitudinal or transverse direction will be greatly reduced, that is, the probe laser pulse is depolarized by the dynamic PDG. The depolarized laser beams may play an important role in mitigating parametric instabilities relevant to laser-driven ICF. More importantly, the scenario of laser depolarization may spontaneously take place for crossed laser beams that are employed as the drivers in either indirect- or direct-drive ICF schemes.
Acknowledgement
This work was supported by the National Natural Science Foundation of China (Grant Nos. 11975154, 11991074, 12005287 and 12135009), the Strategic Priority Research Program of Chinese Academy of Sciences (Grant Nos. XDA25050100 and XDA25010100), the Presidential Foundation of the China Academy of Engineering Physics (Grant No. YZJJLX2016008) and the Science Challenge Project (Grant No. TZ2018005). Simulations have been carried out on the Pi supercomputer at Shanghai Jiao Tong University.
Open access
