Abstract

Paraxial theory of relativistic self-focusing of Gaussian laser beams in plasmas for arbitrary magnitude of intensity of the beam has been presented in this paper. The nonlinearity in the dielectric constant arises on account of relativistic variation of mass. An appropriate expression for the nonlinear dielectric constant has been used to study laser beam propagation for linearly/circularly polarized wave. The variation of beamwidth parameter with distance of propagation, self-trapping condition, and critical power has been evaluated. The saturating nature of nonlinearity yields two values of critical power of the beam ( and ) for self-focusing. When the beam diverges. When the beam first converges then diverges and so on. When the beam first diverges and then converges and so on. Numerical estimates are made for linearly/circularly polarized wave applicable for typical values of relativistic laser-plasma interaction process in underdense and overdense plasmas. Since the relativistic mechanism is instantaneous, this theory is applicable to understanding of self-focusing of laser pulses.

1. Introduction

The interaction of ultrahigh-power laser beams with plasmas is not only of technological importance but also rich in a variety of nonlinear phenomena. These phenomena become particularly interesting and involved when the laser power is sufficiently intense to cause the electron oscillation (quiver) velocity to become relativistic [1]. An important process that can affect the size of the focused spot of the radiation is self-focusing. In recent years the field of laser plasma interaction in the relativistic regime has been identified as an emerging area and is often referred to as high-field science. Several complex phenomena are included in the area of high-field science. These phenomena (on account of their nonlinear nature) are all significantly affected by the field distribution in the beam, and hence, self-focusing occupies a unique position in the field as it affects all other phenomena.

Relativistic electron motion in a plasma due to an intense laser pulse modifies the refractive index and leads to two effects: relativistic induced transparency and relativistic self-focusing. In dense plasma with , light cannot propagate and is reflected from the surface. However, for relativistic intensities generating large factors, the plasma becomes transparent. The dependence of the electron mass on intensity causes significant change of characteristic properties of the plasma and its nonlinear processes; the increase in effective electron mass decreases the electron plasma frequency and makes the penetration of electromagnetic waves into the overdense area of plasma, which is known as self-induced transparency or electromagnetically induced transparency, and is the cause of relativistic self-focusing. A combination of the relativistic self-focusing and relativistic induced transparency enables transmission of laser energy deep into plasmas which is useful for fast ignition of inertial fusion [2]. This so-called penetration sensitively depends on the focal position of the laser intensity due to the inhomogeneous density profile of the plasma and convergence of the laser pulse by final focusing optics.

For ultrafast laser pulses lasting the order of a picosecond or less, the drift velocity of electrons in a plasma can be comparable to the velocity of light, causing a significant increase in the mass of the electron and consequently in the effective dielectric constant of the plasma. This leads to self-focusing of the laser beam, as pointed out [3]. The fact that the relativistic mechanism is the only mechanism of self-focusing that can manifest itself for subpicosecond pulses makes the understanding of this mechanism very important [4, 5]. The earlier analysis concerned with relativistic self-focusing of laser beams in plasmas ignored the effect of the saturating nature of the nonlinear dielectric constant and was essentially a perturbation treatment, based on the quadratic dependence of the dielectric constant on the electric field of the beam [6]. Hence, such an analysis has limited applicability to understand the self-focusing of beams with arbitrarily large electric field and consequently arbitrarily large nonlinear part of the dielectric constant.

Relativistic focusing of radiation beams in plasma is applicable to pulse length in the range , where is the ion-acoustic speed [7], and can arise owing to the creation of a density depression in the plasma as well as to an increase in the electron mass due to relativistic effects. The time scale for self-focusing due to a density depression is longer than that for self-focusing due to relativistic mass increase. The latter time scale is , where is the radiation frequency and is assumed to be much greater than the electron plasma frequency. The analysis presented here will be concerned only with relativistic self-focusing on a time scale sufficiently short that the plasma density profile does not evolve significantly under the influence of the radiation beam. This implies that the pulse length of the radiation beam must be short compared with (the time scale for the density depression to occur) and, of course, long compared with a radiation period . We make use of the steady-state paraxial-ray theory, as developed by Akhamanov et al. [8] and Sodha et al. [9].

The organization of the paper is as follows. In Section 2 the nonlinear dielectric constant due to relativistic variation of mass for a time-harmonic plane wave (i.e., a linearly/circularly polarized wave) is found. Section 3 is concerned with the relativistic self-focusing equation relating the variation of beamwidth parameter with distance of propagation. The self-trapping condition and the critical power are evaluated in Section 4. Results and discussion are made in Section 5.

2. Nonlinear Dielectric Constant of the Plasma

The relativistic equation of motion for a point charge is where is the rest mass of a point charge, is the velocity of light in vacuum, E is the electric field, and B is the wave magnetic field. For a homogeneous plane wave traveling in the direction of a unit vector n, the electric and magnetic vectors are perpendicular to n and are related by Substituting for B on the right-hand side of the equation of motion, (1) gives Considering the scalar product of (3) with n yields Equation (4) together with the energy equation establishing that is constant. We now eliminate in favour of the phase variable Here is the angular frequency of the wave, and r is the positive vector of the particle. Then where the prime denotes differentiation with respect to , and using (6) Substitution of (8) and (9) into (3) enables the factor to be cleared and leaves which is a second-order linear differential equation for r as a function of . The dielectric constant of the plasma is given by where   is the relativistic mass, is the relativistic factor, and is the plasma frequency in the absence of the beam. For a circularly polarized wave, the components of (10) are where   is the phase, and and approximately satisfy the linear dispersion relation Using (12), the average relativistic factor is where for self-focusing considering the order of electric field up to second order and that the value of is unity for circularly polarized wave [10] and half for linearly polarized wave. The dielectric constant can be written as [9] where is the linear part of dielectric constant, is the nonlinear relativistic term, and The effective dielectric constant for arbitrary large nonlinearity can be written as where is the dimensionless beamwidth parameter and and represent, respectively, the linear and nonlinear dielectric constant for arbitrary nonlinearity.

3. Relativistic Self-Focusing Equation

Consider the propagation of Gaussian laser beam of frequency along -direction; at the intensity distribution of the beam is given by and here, “” is the radial coordinate of the cylindrical coordinate system, “” is the initial beamwidth, and is the axial amplitude. The wave equation governing the electric vector of beam in plasmas with the effective dielectric constant (20) is To solve the prevoius equation, the Wentzel-Kramers-Brillouin (WKB) approximation has been used. In the WKB approximation, the second term of (22) has been neglected, which is justified when . The electric vector of the main beam can be expressed as In the presence of the modified background electron concentration due to relativistic nonlinearity, the intensity distribution of the beam inside the plasma can be obtained by using the WKB and paraxial ray approximations. Following Akhmanov et al. [8] and Sodha et al. [9] the solution for E can be written as where Substituting for E and from (24) and (20) in (23), one obtains Putting in (26) and separating real and imaginary parts we get where is called the eikonal and related to the curvature of the wavefront. For a slightly converging/diverging beam the solution for a Gaussian laser beam (21) can be written as where It can be seen previously that is the inverse of radius of curvature of the wavefront and is the width of the beam. In the geometrical optics approximation, represents a ray in a plane containing the -axis. On substituting for into (27), using the paraxial ray approximation, such that , equating the coefficients of on both sides of the resulting equations and substituting for into (29), the dimensionless beamwidth parameter is given by Transforming the coordinate and the initial beamwidth to dimensionless forms we get the characteristic beam propagation equation as Physically, (34) governs the variation of beamwidth parameter with distance of propagation. Substituting for from (24) and , , and from (29)–(31) into (20) and using the paraxial ray approximation, can be written as correct to terms in , where Equation (35) can put into the convenient form of (20), with [11] Consider second order of nonlinearity, for circularly polarized wave; we can rewrite the previous set of equations as On substituting values of and from (38) in (34), the equation governing the beamwidth parameter is In (39) the first term on right-hand side is responsible for diffraction divergence of the beam, and the second term corresponds to convergence, arising due to nonlinearity through the relativistic factor of linearly or circularly polarized beam, which contributes to focusing.

4. Self-Trapping Condition

When the two terms on the right-hand side of (39) balance each other the beam propagates without convergence or divergence, which is referred to as uniform waveguide propagation or self-trapping of the beam. For an initial plane wavefront of the beam the initial conditions on are and . As both the terms on right-hand side cancel each other hence at , ; since is also zero and at , for all values of ; in other words beam can propagate without convergence or divergence. Therefore, the condition for self-trapping is for linearly polarized beam and

for circularly polarized beam. The critical power of the beam is

5. Results and Discussion

The differential equation (39) is to study the propagation of linearly and circularly polarized beams through plasma. We have chosen parameters accessible for interaction of such beams with moderately underdense plasma to slightly overdense plasma. The numerical values are chosen for laser intensity ranging between I = 1016 to 1018 W/cm2, electron density between 0.1 to 1.5 of critical density, laser frequency  sec−1 and initial radius of the beam laying between 1 to 3 μm. The dielectric function given by (20) and beam propagation equation (39) together with boundary conditions, are numerically solved. Further we consider the propagation of the beam in axially inhomogeneous plasma by including linear density variation in plasma.

Figure 1 shows the variation of relativistic factor versus intensity for linearly and circularly polarized laser beam. From the graph it is evident that the growth in relativistic factor for circularly polarized laser beam is more as compared to linearly polarized beam. Figure 2 shows plot of dielectric function versus axial intensity for linearly and circularly polarized beams in both underdense and overdense plasmas corresponding to different concentration of electron density. It is seen that initially there is an increase in dielectric function and after a certain value it attains a constant value. Because of very strong fields, the nonlinearity saturates. As a result the electromagnetic wave drives all of the plasma out of the regions of large field intensity and establishes equilibrium with radiation inside a vacuum channel surrounded by plasma. As we move from underdense region to overdense region the saturation effect slows down.

The critical curve (i.e., a plot of normalized radius of the beam “ ” against the axial intensity proportional to critical beam power) (40a) and (40b) is represented by Figure 3 for linearly and circularly polarized beams at . Equations (40a) and (40b) have two critical powers and , corresponding to two different values of electric field strength (or intensity of electric field), say and such that . The beam can be self-focused when its power lies between the two critical values , and the range increases rapidly with increasing normalized radius. The region above the critical curve is the self-focusing region, and region below the critical curve is known as the defocusing region. The effect of different concentration of electron density can be clearly seen from these critical curves. As concentration of electron density increases, that is, if we move from underdense region to overdense region the critical curve shifts downward.

Figures 4(a) and 4(b) show the variations of dimensionless beamwidth parameter with dimensionless distance of propagation at relativistic intensity in plasma for different values of electron density in underdense and overdense plasmas. From the critical curves for = 0.7 () and = 4, the value for linearly/plane polarized beam, and for circularly polarized beam. At () and , for linearly/plane polarized beam, and for circularly polarized beam. The beam shows focusing/defocusing for different values of . This is because as we change the strength of the electric field, the maxima and minima get shifted. For a given value of when , the second term on the RHS of (39) dominates the first term at , and is negative; hence, decreases with the distance of propagation.

For certain value of , the two terms on RHS of (39) get cancelled resulting in = 0, and for larger values of , RHS becomes positive. However, the beam still continues to converge on account of the curvature it has already gained, though becomes less and less negative. At certain , ( and . Beyond this point, the diffraction term dominates over the nonlinear term, and hence exhibits the oscillatory behavior. From the figures it is evident that due to relativistic mass effect, increase in effective electron mass decreases the electron plasma frequency and makes the penetration of the beam into the overdense region of plasma, which is known as self-induced transparency. At higher concentration of electron density beam gets more focused as compared to low concentration of electron density. The variation of dimensionless axial intensity of the beam with distance of propagation corresponding to Figures 4(a) and 4(b) for self-focusing region is depicted in Figures 5(a) and 5(b). The intensity peaks for higher concentration of electron densities are higher with a reduced distance between two successive peaks.

Figure 6(a) illustrates the variations of dimensionless beamwidth parameter with distance of propagation for different concentrations of electron densities in self-focusing region with linear increase in electron density. Careful observation of the results from figure indicates that for linearly increasing density plasma, minimum value of the normalized beamwidth parameter for the second and higher orders decreases continuously in underdense and overdense plasmas. In underdense plasma it propagates long distance as compared to overdense plasma, but its focusing is less. The corresponding graph for variation of dimensionless axial intensity (/) of the beam versus distance of propagation is shown in Figure 6(b). From the numerical analysis and graphs it is evident that both linearly polarized and circularly polarized beams can propagate through underdense plasma as well as they penetrate overdense plasma or say they exhibit self-induced transparency effect.

In underdense plasma there is not much difference in propagation of linearly polarized and circularly polarized beams, while the propagation is affected in case of overdense plasma and focal point changes.

A threshold exists for onset of self-focusing, as this effect must overcome the spreading of the beam due to diffraction. A light beam is focused to be self-trapped at any arbitrary diameter and will thus not spread. Further, self-trapping occurs at a critical power level, which is independent of the beam diameter. The theory presented here has a much wide range of application than the previous one, because it is not a perturbation treatment. Furthermore, it is applicable at all intensities—high and low.