1. INTRODUCTION
Self-phase-modulation (SPM) of laser is an important nonlinear process in a variety of media with intensity dependent index of refraction. Max et al. (Reference Max, Arons and Bruce Langdon1974) observed the nonlinear frequency shift of a strong electromagnetic wave in plasma, due to weak relativistic effects. Yablonovitch (Reference Yablonovitch1974a, Reference Yablonovitch1974b) investigated SPM and short pulse generation from laser breakdown plasmas. Tsintsadze et al. (Reference Tsintsadze, Tskhakaya and Stenflo1979) obtained the generalized dispersion relation for modulational instabilities due to relativistic electron mass variations. The SPM is accompanied by self focusing/defocusing of the laser beams and is significantly modified by those effects. Gill and Saini (Reference Gill and Saini2007) observed the enhancement of Raman scattering in collisional plasma by the interaction of rippled laser beam with upper hybrid mode and also observed the focusing of the upper hybrid wave. Willi et al. (Reference Willi, Toncian, Borghesi, Fuchs, D'Humieres, Antici, Audebert, Brambrink, Cecchetti, Pipahl and Romagnani2007) presented a novel technique for focusing and energy selection of MeV protons by employing a hollow micro-cylinder, which irradiated at the outer wall by a high intensity ultra-short laser pulse. Purohit et al. (Reference Purohit, Chauhan and Sharma2008) studied the excitation of an upper hybrid wave by a relativistic laser beam in the presence of perpendicular static magnetic field and also introduced the relativistic electron mass nonlinearity as well as relativistic self-focusing effect. Liu and Tripathi (Reference Liu and Tripathi2000) have developed an unified formalism of SPM and self defocusing in a tunnel ionizing plasma. The self defocusing limits the frequency upshift of the laser. Liu and Tripathi (Reference Liu and Tripathi2001) studied the frequency downshift of a Gaussian beam due to relativistic self focusing of the laser.
Watts et al. (Reference Watts, Zepf, Clark, Tatarakis, Krushelnik, Dangor, Alott, Clarke, Neely and Norreys2002) observed the relativistic SPM in the interaction of a high-intensity laser pulse (1 ps, 1.053 µm, 80 J) with the plasma. Singh et al. (Reference Singh, Konar and Sharma2001) observed the resonant cross modulation of two laser beams in a semiconductor slab. Saini and Gill (Reference Saini and Gill2006) studied the self focusing and SPM of an elliptic Gaussian laser beam in collisionless magnetized plasma, by using the variational approach considering the effects of nonlinearity and diffraction. Gupta and Suk (Reference Gupta and Suk2007) obtained the electron acceleration by two crossing chirped lasers of same amplitude and frequency, at an arbitrary angle, and causing modulation of laser intensity. Recently, Nickles et al. (Reference Nickles, Ter-Avetisyan, Schnuerer, Sokollik, Sandner, Schreiber, Hilscher, Jahnke, Andreev and Tikhonchuk2007) experimentally observed the interactions of a laser pulse at intensities above 1019 W/cm2, and covered the rear and front side acceleration mechanisms, particle dynamics inside the dense target, proton source characteristics, strong modulations in proton and deuteron emission spectra, and generation of quasi-monoenergetic deuteron bursts.
Liu and Jetendra (Reference Liu and Jetendra2006) studied the self-defocusing/focusing of a right circularly polarized laser, in a preexisting density channel, and studied the effect of ponderomotive force driven plasma wave on the growth of modulational instability. Borghesi et al. (Reference Borghesi, Kar, Romagnani, Toncian, Antici, Audebert, Brambrink, Ceccherini, Cecchetti, Fuchs, Galimberti, Gizzi, Grismayer, Lyseikina, Jung, Macchi, Mora, Osterholtz, Schiavi and Willi2007) studied the probing technique to provide the maps of impulsive electrostatic fields with high spatial and temporal resolution by the high-intensity laser matter interactions. The dynamics of ponderomotive channeling in underdense plasmas is also observed by the processes of Debye sheath formation and MeV ion front expansion at the rear of laser-irradiated thin metallic foils. Hafizi et al. (Reference Hafizi, Tang, Sprangle and Hubbard2000) studied the relativistic focusing and ponderomotive channeling of intense laser beams—ponderomotive channeling modified the effective potential of the laser spot size—and obtained the envelope equation by using the source-dependent expansion method with Laguerre-Gaussian Eigen functions. Experimental observations of relativistic focusing and ponderomotive channeling have been reported in the literature (Faenov et al., Reference Faenov, Magunov, Pikuz, Skobelev, Gasilov, Stagira, Calegari, Nisoli, De Silvestri, Poletto, Villoresi and Andreev2007; Torrisi et al., Reference Torrisi, Margarone, Laska, Krasa, Velyhan, Pfeifer, Ullschmied and Ryc2008; Chessa et al., Reference Chessa, Mora and Antonsen1998; Sun et al., Reference Sun, Ott, Lee and Guzdar1987; Kurki-Suonio et al., Reference Kurki-Suonio, Morrison and Tajima1989; Borisov et al., Reference Borisov, Borovskiy, Korobkin, Prokhovrov, Rdodes and Shiryaev1990, Reference Borisov, Borovskiy, Shiryaev, Korobkin, Prokhorov, Solem, Luk, Boyer and Rhodes1992a, Reference Borisov, Borovskiy, Korobkin, Prokhorov, Shiryaev, Shi, Luk, Mcpherson, Solem, Boyer and Rhodes1992b; Konar & Manoj, Reference Konar and Manoj2005; Abramyan et al., Reference Abramyan, Litvak, Mironov and Sergeev1992; Chen & Sudan, Reference Chen and Sudan1993; Annou et al., Reference Annou, Tripathi and Srivastava1996; Tzeng & Mori, Reference Tzeng and Mori1998; Monot et al., Reference Monot, Auguste, Gibbon, Jakober, Mainfray, Dulieu, Louis-Jacquet, Malka and Miquel1995).
In this paper, we study the effect of self focusing on SPM of a laser pulse in self-channeled plasma. The self focusing, caused by relativistic mass shift as well as radial ponderomotive force, modifies the electron density. In Section 2, we obtain the nonlinear current density due to a high power laser. In Section 3, we drive coupled equations for amplitude and phase of the laser pulse. In Section 4, we solve these equations in near-axis approximation and by expanding the Eikonal up to second order in radial coordinates r.
2. NONLINEAR CURRENT DENSITY
Consider the propagation of a laser pulse in a plasma of electron density n 00 along the z-axis. At z = 0, the electric field of the laser is
where g(t) is the temporal shape of the pulse For z > 0, we may write,
where A(t, z, r) is the complex amplitude, φ (t, z) is the fast phase of the wave, and one may take ω = ∂φ/∂t and k = −∂φ/∂z.
The laser pulse imparts an oscillatory velocity to the electrons,
where −e and m are the electronic charge and mass, respectively, and c is the velocity of light in vacuum. It also exerts a relativistic ponderomotive force on the electrons ,
The ponderomotive force expels the electrons away from the region of higher electric field, while the ions remain stationary due to their heavy mass. On the time scale longer than the electron plasma period and shorter than the ion plasma period (ωp−1 < t < ωpi−1), one may take the static space charge potential, caused by the displacement of electrons, φs ≈ −φp. Using this in the Poisson equation, ∇2φs = 4π e(n e − n i), we get
The nonlinear electron current density at the frequency of the laser can now be written as,
3. COUPLED EQUATIONS FOR AMPLITUDE AND PHASE
The wave equation governing the propagation of the laser in underdense plasma is
We get
where ωp2 = (4πn 00e 2/m) is the plasma frequency. The second term on left-hand-side arises due to relativistic mass variation and the third is due to the ponderomotive force.
On substituting from Eq. (2) in the wave equation and assuming the wave amplitude A to be a slowly varying function of t and z (Wentzel-Kramer-Brillouin approximation) we obtain,
where γ0 = γ at r = 0.
Differentiating Eq. (9) with respect to t by use of ∂k/∂t = −∂ω/∂z,
where is the group velocity.
We define a function F = (ω/ω0)1/2A, t / = t − z/c and z / = z, Eqs (10) and (11) become
and
We may write F = F 0 exp(iS), where F 0 (t /, z /, r) and S(t /, z /, r) are real and separate the real and the imaginary parts of Eq. (11),
and
Equations (13), (14), and (15) are the coupled phase and amplitude equations.
4. FREQUENCY DOWN SHIFT
We solve Eqs (14) and (15) by expanding S as S = S 0 + S 2r 2/r 02 and introducing a function f(z /) such that S 2 = (ωr 02/2c)(1/f)(∂f/∂z /). Then Eq. (15), for an initially Gaussian profile of laser intensity, gives
where (r 0f) is the modified radius of the laser beam and the beam width parameter f has implicit dependence on t /.
Using Eq. (16) in Eq. (14) and collecting the coefficient of r 2, we obtain
where v0/c = eE 00/mω0c is the normalized oscillatory velocity and ξ = z //R d, Ω = ω/ω0, R d = (ω0/c)r 02.
The boundary conditions at ξ = 0 for an initially plane wavefront are f = 1, ∂f/∂ξ = 0, Ω = 1, and S 0 = 0. The first term on the right-hand-side of Eq. (17) represents the diffraction divergence, the second term gives the beam convergence due to relativistic mass non-linearity, whereas the third and forth terms describe the convergence due to ponderomotive nonlinearity. We choose the temporal profile g(t /) = tanh (t //τ) for t / > 0 and zero otherwise, where τ is the pulse rise time, and solve the coupled Eqs (17) and (18) numerically for the parameters: r 0ωp/c = 5, ω0τ = 40, v0/c = 1, γ0 = √2, T ≅ (t − z/c)/τ = 0 − 3, z/R d = 0 − 3.
Figure 1 shows the variation of the beam width parameter f with the normalized distance of propagation ξ. At T = 0, f increases monotonically with ξ due to diffraction. At higher T, the nonlinear convergence exceeds the diffraction divergence, and f decreases with ξ, i.e., focusing occurs. After certain distance, the electrons are significantly depleted from the axial region, and the nonlinear convergence weakens. At a certain value of ξ, the beam starts to diverge due to the predominance of diffraction divergence. As the beam acquires a large spot size, the nonlinear convergence again predominates beyond ξ = 0.7. Figure 2 shows the frequency of the laser pulse as a function of retarded time T at z/R d = 1. The front of pulse undergoes frequency downshift at T = 0. With increasing T, ω/ω0 increases monotonically and attains a saturation value. In the simultaneous presence of relativistic mass nonlinearity and ponderomotive nonlinearity, self focusing is stronger; hence red shift remains fairly constant with retarded time throughout the channel. Figure 3 shows the variation of the normalized axial intensity |E|r=02/E 002 as a function of T at z = R d and z = 0. At z = 0, the normalized axial intensity remains tanh (T) with retarded time, but at z = R d due to periodic self focusing, it deviates from tanh (T) and increases monotonically with the retarded time. For the normalized retarded time T = 1.5, the variation of the radial intensity |E|2/|E|r=02 with (r/r 0) at z = R d and z = 0 is shown in Figure 4. Due to periodic self focusing, the radial intensity is also more narrowly peaked at z ~ R d than at z = 0, but its shape remains Gaussian at z = R d and z = 0.
5. DISCUSSION
A fast rising laser pulse, with tanhyperbolic temporal profile and Gaussian radial profile, propagating through preformed plasma undergoes downshift in frequency. Maximum downshift occurs at t / = (t − z/vg) ~ 0, i.e., at the foot of the pulse. However, self focusing significantly modifies the temporal profile of the pulse, hence influences the frequency shift. When relativistic mass nonlinearity alone is considered, the frequency shift diminishes as one approaches the higher intensity of the pulse. With the inclusion of ponderomotive nonlinearity the frequency downshift continues to be significant up to longer times due to the temporal variation of electron density in the axial regions. Normalized axial intensity deviates from tanh (T) form. The deviation increases with retarded time and is maximum at z ~ R d (Fig. 3).