A Study of Beam Parameters Using NLSE in Chalcogenide Glass Through Variational Method with a Gaussian Trial Function
Chironjit Hazarika1, *, Abhijeet Das2, Subrata Hazarika2
1Department of Physics, Mariani College, Jorhat, Assam, India
2Department of Physics, Assam University, Diphu Campus, Karbi Anglong, Assam, India
To cite this article:
Chironjit Hazarika, Abhijeet Das, Subrata Hazarika. A Study of Beam Parameters Using NLSE in Chalcogenide Glass Through Variational Method with a Gaussian Trial Function. American Journal of Optics and Photonics.Vol.3, No. 4, 2015, pp. 43-47. doi: 10.11648/j.ajop.20150304.11
Abstract: Using variational method for an elliptical Gaussian optical beam trial function, self –action in bulk chalcogenide glass (Kerr media) is investigated. Emphasis is laid on the study of variation in beam width, curvature, phase and intensity of the beam with propagation distance. Solutions predict stationary self-focusing of the elliptical beam and an effective beam collapse at 10Pcr input power. These study is significant in the choice of parameters in optical communications.
Keywords: Self- Focusing, Kerr Nonlinearity, Gaussian Laser Beam
Self- focusing, an optical self- action effect, occurs when the propagating electromagnetic field induces a change in the refractive index of the non-linear medium through which the field propagates. Self-focusing and self-trapping effects of laser pulses and beam can lower the threshold of other non-linear processes like self-phase modulation, optical damage, regulate design of high energy lasers sources required for long – distance optical communication system , controlled laser fusion [2,3], coherent beam combiningand molecular dynamic diagnostics [8,9] etc. These consequences of self-action effects on the propagation of intense laser pulses and beams have attracted attention of researchers for more than fifty years now. Intense researches on laser propagation in nonlinear media and associated self-action effects lead to formulation of different theoretical concepts and mathematical models of self-focusing, beam self-trapping, filamentation, supercontinuum generation and filament plasma defocusing [5,6,7,10]. Starting with Chiao et. al.’s paper on self–trapped model of an intense laser beam in di-electric waveguide , which was an explanation to Hercher’s experimental observation of first ever ‘optical damage’ in 1964 caused by a Q-switch laser in glass , more than 1500 research papers involving self-focusing and self-trapping in one way or other  have been published since the concept of laser beam self-focusing and self-trapping was proposed by Askaryan in 1962 . Theoretically, the temporal and spatial electromagnetic-field behaviour of a propagating beam in nonlinear bulk media were investigated by several methods that includes beam propagation method (BPM) [14,15], nonlinear Schrödinger equation (NLSE) [16,17], finite-difference time-domain (FDTD) method [18,19] etc. Currently, frontier theoretical research on laser beam propagation in nonlinear medium deals with formulations and creative solutions of nonlinear Schrodinger equation (NLSE) that describes laser beam collapse, self-trapping, dispersion, filamentation, modulation instability, pulse splitting and other extraordinary particularities of nonlinear beam propagation [20,21,22]. In this communication propagation of intense elliptical Gaussian beam in chalcogenide glass with refractive index (=2.8) is investigated as the beam approches the self focusing distances (zsf) by variational method. Using Gaussian trial function, equations are derived for beam width, curvature, phase shift as a function of propagation distance to study their variations and interpretions of self action effects.
The present study characterizes its significance as currently non-linear effects are exploited at low optical powers in chalcogenide glasses. Due to their large intensity dependent refractive index and their low linear absorptions at optical communication wavelength makes these glasses attractive materials for use in a variety of optical devices [23,24,27].Hence Study of the behaviour of the intense beam in chalcogenide glass at self-focusing distance can facilitated the choice of power, wavelength etc. and fabricate symmetric waveguides for efficient optical transmission.
This paper is organized as follows: in section 2, the theoretical formalism of the problem is outlined, section 3 is devoted to results and discussion and paper is concluded in section 4.
2. Basic Formulation
The propagation of a light beam in a Kerr medium under the paraxial approximation can be described by the nonlinear Schrodinger equation,
Where k is the linear wave vector, x and y are the transverse co-ordinates, z is the longitudinal co-ordinates,is the third order nonlinear coefficient and is the linear refractive index of the medium. This equation describes the slowly varying field inside the nonlinear medium. The first term and second term of equation (1) represents the diffraction of the beam in transverse direction and the last term works as a source to generate different modes of the electric field due to the nonlinear response of the medium. Kerr nonlinearity is one of the important nonlinear response occurring in glasses, gases etc. When the diffraction term comes into balance with the focusing due to Kerr nonlinearity then self-focusing of optical beam occur.
The solution of the variational problem
also solves the nonlinear Schrodinger equation (1) with an appropriate Lagrangian L. Equation (1), is reformulated as a variational problem using the Lagrangian
where ‘ * ’ denotes complex conjugation.
As a trial solution, the Gaussian field distribution function is considered for the procedure of the variational approach
E(x, y, z) = (4)
where (z) and are the amplitude and phase of the complex field (x, y, z), and are the width parameters of the Gaussian beam in x and y direction respectively, (z) and (z) are the curvature parameters in x and y direction.
Inserting the trial function into equation (3), and integration with respect to x and y yield the effective Lagrangian.
The effective Lagrangian can be written as,
Now from the reduced variational principle,dz =0, i.e considering the beam propagation along z direction, it is possible to derive variational equations with respect to ,,(z),(z) , Using the procedure of [25,26] and after some rearrangement, the following differential equations for α, β, (z),(z) and are obtained, which describe the evolution of the Gaussian beam in Kerr medium.
(z) = (6)
(z) = (7)
only depends on z, So we can write ≈ and ≈ .
Differentiating equation (6) and (7) and putting the values of in equation (9) and (10) yield,
Equations (12) and (13) will be used for numerical analysis.
3. Result and Discussion
The finite difference scheme based on fixed point iteration method is used to obtain numerical solutions of equations (12) and (13) that describe the spatial evolution of beam width. Here the critical power is calculated using the relation .
The solutions are plotted in Fig. 1 and are for a beam propagating in chalcogenide glass of refractive index (n0)=2.8 and non-linear refractive index)of 13x /W at three input powers viz. 1.5, 3 and 10 times critical power Pcr. which is 9kw for the wavelength λ= 1500nm.Using the boundary condition α= for values of predict an exponential decrease in α and β which approaches zero over a finite propagation distance. i.e. when > the width of an elliptical beam collapses over a finite propagation distance which for 1.5, 3 and 10 times Pcr is 8.3cm, 6.3cm and 4 cm, respectively. However, from Eqns. (12) and (13) a necessary condition for self-trapping can be written as
at z = (14)
Where is self –focusing length. It implies that at , α = β which is a condition at self-focusing. The values of thus determined from the derived condition are 7.1 cm, 5.6 cm and 3.0 cm for 1.5, 3 and 10 times Pcr, respectively. In these solutions value of k was taken as
Figure 1. Variation of beam widths and β with propagation distance for input power 1.5, 3, 10 times the critical power for =.02, = .015.
1.4[28,29] and the wavelength of propagation as 1500 nm - the‘eye-safe region’ wavelenth in optical communication.It is observed from Fig. 1 that the beam width collapses more rapidly with input power and the collapse distance decreases with increase in input power of the Gaussian beam. The collapse distance scales as for P>> but at higher powers collapse scales as [31,33].Collapse occurs when the whole beam shrinks to a point. Further, a collapsing beam exhibits sharp rise in intensity .The intensity of the propagating beam at three input powers in the present case increases on reaching zsf . But, the rise in intensity of the beam with input power 10only bears signature of a collapsing beam during self- focusing. The variation of intensity with propagation distance is presented in Figure 2 and is plotted using Eqn.(11). The values of H0 for the three input powers are 0.1974x, 0.3948x and 0.1316x(in c.g.s).
Figure 2. The variation of with propagation distance at input power 1.5and 10.
The threshold power for collapse of a beam for an elliptical beam is approximated as times the threshold power for circular beam,.Here e = β/α .The threshold power for collapse increases with input beam ellipticity and the threshold power for collapse is above the critical power (35).
Using the boundary condition at z=0, α= and β= , the solutions for Equations (6) and (7) are given by,
The variation of x and y direction curvature and with propagation distance at input power 1.5,3 is shown in figre(3) and (4).
Figure 3. Variation of curvature with distance of propagation for input power 1.5, 3, 10 times Pcr..
Figure 4. Variation of curvature with distance of propagation for input power 1.5, 3, 10 times the critical power.
Evolution of both the curvatures c1 and c2 shows a decreasing trend with propagation distance and falls of rapidly around half the collapse distance. The variations are accentuated as input power increases. The curvatures diminish three folds more for 8.5 times increase in input power on propagating 4.3 cm less. The negative sign in the magnitude of c1 and c2 is indicative of the decrease in curvatures with propagation distance i.e convergent point, the radius of the wave front of the transmitted beam is negative.
Again, using the boundary condition ,at z=0, α= and β= , the solution of equation (8) is given by,
ϕ(z) =[ + ] (17)
Figure 5 displays the variation of longitudinal phase ϕ as a function of propagation distance for input power 1.5, 3, 10 times Pcr. It is observed that phase ϕ(z) is negative and for the beam with Pin =10 phase changes abruptly at zsf . It should be noted that the longitudinal phase may be positive or negative depending on the value of intensity parameters.
The evolution of the beam width,curvature and phase was analysed by variational approach. The solutions obtained using finite difference scheme based on fixed point iteration predict stationary self-focussing for elliptical Gaussian beam propagating in chalcogenide glass with distinct beam collapse evident at 10 Pcr input power. The variations in curvature, phase and intensity were more accentuated for beam with input power of 10Pcr compared to those in input power of 1.5Pcrand 3Pcr and occurs as the beam approches self-focusing length or at self-focus length.It is indicative that self-action effects are more effective at input power of 10 Pcr .The change of phase is negative with distance of propagation.The present study is significant for optical communication leading to transfer of energy in narrow beams without losses at low optical power(for silica fibre =2.1MW).Further,the study of beam parameters at self focusing distances also helps in material selection and material modification in the fibre optic communications,because properly selecting beam parameters, like power,wavelength and Kerr nonlinear material can prevent the damage of fibre and waveguide etc.