Generation of Dynamic Terahertz (THz) Pulse Profiles for Fourth Order Harmonics in A Fiber – Like Nonlinear Uniaxial Crystal.

DOI : 10.17577/IJERTV3IS060583

Download Full-Text PDF Cite this Publication

Text Only Version

Generation of Dynamic Terahertz (THz) Pulse Profiles for Fourth Order Harmonics in A Fiber – Like Nonlinear Uniaxial Crystal.

Timtere Pascal

Department of Physics, Modibbo Adama University of Technology, Yola

Adamawa State, Nigeria.

Adam Usman

Department of Physics, Modibbo Adama University of Technology, Yola

Adamawa State, Nigeria.

Yakubu, Hammanjoda. Jagami Department of Pure Science, Taraba State Po;ytechnic, Jalingo Taraba State, Nigeria.

Abstract In This work, the Maxwells equations were used to obtain a transverse wave equation. A trial solution of the wave equation that has both the amplitude and the field distribution function of the terahertz (THz) radiation were used to generate a dimensionless nonlinear Schrödinger equation (NLSE) for the fourth order harmonics. The numerical simulation of the NLSE for the THz radiation was tested on a fiber like nonlinear uniaxial crystal based on three parameters. The nonlinear

spatial solitary waves based on the nonlinear Schrödinger equations associated with them. For propagation along the z axis, the transverse dimension for the temporal solitary wave

is in terms of time t, and the dimensionless pulse envelope U

(z,t) evolves according to the NLSE as [5]

coefficient, the dimensionless pulsewidth, and the order of the

U 2U 2

solitary wave. The dynamic pulse profiles show that THz could

i 2

U U 0


have an application in optical communications, just like Laser in

z 2

t 2 NL

optical fiber communication.

keywords terahertz, uniaxial crystal, Maxwells equations,



d 2 k

d 2

is the second derivative of the linear

nonlinear Schrödinger equation, pulse profiles


    The THz radiation (~ 1 x 1012 Hz) is broadly applied to sub millimeter wave energy that fills the wavelength between 0.001 m to 0.0001 m, corresponding to frequency range between 0.3 THz 3 THz [1]. Spectroscopists had coined the term for emission frequencies that fall below the far infra red (FIR) and the microwave regions of the electromagnetic spectrum. The ability of THz radiation to penetrate packaging, similar to X rays but without ionizing the contents makes THz potentially useful for law enforcement applications such as detecting explosives and illegal drugs [2]. Nonlinear optical phenomena occur typically at high optical intensities [3], and most of the wave equations involved are governed by the Nonlinear Schrödinger Equations (NLSE). Optical solitons have been regarded as the next generation technology for high capacity optical communications, chiefly because of their promise to transmit signals over long distances. Agrawal [4] describes the soliton as the particle like properties of pulse envelopes in dispersive nonlinear medium that under certain conditions, the envelope not only propagates undistorted but survives collisions just as particles do. There are basically two types of solitary waves. Xie [5] distinguished the temporal and the

    wave number k with respect to angular frequency , and vNL is the nonlinear coefficient. 2 is called the group velocity dispersion (GVD) that changes sign as one moves from normal to anomalous dispersion. It also requires a material medium and/or a waveguide in order to exist. For 2 = 0, the pulse will experience the nonlinear change as a self phase modulation (SPM). So, when nonlinearity counteracts or balances dispersion in nonlinear media, one obtains temporal solitary wave. Maxwells laws govern the interaction of bodies which are electrically charged. The bodies and their charges may either be stationary or in motion [6]. Materials in which two of the components of the dielectric constant are equal (such as x y z ) are termed uniaxial crystals.

    This class includes trigonal, tetragonal and hexagonal crystals. Syms and Cozens [7] explained that for uniaxial crystals, two of the refractive indices are identical (nx = ny = no), and nz = ne, and, that this particular ellipsoid can be characterized by two values of refractive indices; no (the ordinary index), which corresponds to nx or ny, and ne (extraordinary index) corresponding to nz. In this work, the Nonlinear Schrödinger Equation (NLSE) for Fourth Order is numerically simulated and stable dynamic pulse profiles for THz radiation are obtained.


    E(x, y, z,t) F(x, y)A(z,t) ei(t z)


    The Nonlinear Schrödinger Equation (NLSE) is an example of a universal nonlinear model that describes many physical nonlinear scenarios [8]. The mathematical description of solitons requires the solution of nonlinear Schrödinger equation (NLSE) in a dispersive medium. The Maxwells

    where A(z,t) is the amplitude of the pulse envelope, is a normalization constant, and F(x,y) is the field distribution in the (x,y) plane and corresponds to the mode structure. Timtere [14] obtained the dimensionless NLSE for the fourth order harmonics as

    equations provide the most fundamental description of electric and magnetic fields. According to Wang [9], THz

    i U



    U U U 3U 0

    pulses generated from an optoelectronic crystal can be calculated using the Maxwells equations. For a neutral dielectric medium (one with no free charges, i.e 0 ), the


    2 2 o o 1

    four Maxwells equations are;

    . D = 0 (2)

    . B = 0 (3)

    E = B (4)



    where U is the dimensionless pulse amplitude (height), is

    the dimensionless pulse distance, v is the nonlinear coefficient, and is the dimensionless time. = ± 1 defines

    normal or anomalous dispersion propagation. For o = ± 1 implies n1 > 0 or n1 < 0, and 1 = ± 1 implies n3 > 0 or n3 < 0. For most nonlinear applications, it is the anomalous dispersion (2 < 0) that is applicable [12]. This is because in

    H =



    the study of nonlinear effects, it is in the anomalous dispersion regime that optical fibers support solitons through a balance between the dispersive and nonlinear effects. This

    where E is the electric field , B is the magnetic field

    ,D E is the electric flux and H is the magnetic intensity. For a non magnetic media,

    means that in this work = – 1. However, in equation (10), o and 1 take positive values of +1. This is because the fourth and the fifth terms that correspond to the second ( 2 ) and

    the fourth ( 4 ) order nonlinearities are positive for the

    B o H


    three uniaxial crystals studied in a work by [15].

    where o is the permeability of free space (~ 4 x10-7 N/A2). Duffin [10] defines the electric flux as



    D o E P


    One requires a numerical simulation of the dimensionless 4th


    is the permittivity of free space (~ 8.85×10-12

    order NLSE to test the authenticity of the NLSE developed.

    2 2 o The other major reason for the simulation is to establish the

    C /N.m ), P is the polarization, and is the dielectric dipole moment per unit volume of the medium. The Maxwells equations are used to obtain the transverse fields wave equation of the form

    balance between the nonlinear coefficient, v in the fifth term of equation (10) and the group velocity dispersion (GVP) in the second term of equation (10). The explicit finite difference scheme is a nuerical approach that is often

    2 E 1

    c 2

    2 E 1



    t 2 2

    2 P

    t 2


    necessary for an understanding of the nonlinear effects in fibre like waveguides [11]. This scheme was used to solve the NLSE generated from equation (10). The values of the

    Equation (8) is the fundamental electromagnetic wave equation of which the THz radiation is an example, and it shows how the electric field depends upon the dipole moment, P of the medium. According to Agrawal [11], the NLSE can be solved exactly using inverse scattering method (ISM) that is similar to the Fourier Transform Method. However, the NLSE belongs to a class of non-integrable partial differential equations that are not amenable to ISM, it is the variational method (VM) that proves to be a more suitable substitute for the ISM [12]. The stationary solitary wave solution for NLSE was used to generate a stationary pulse profile for the THz radiation. The propagation of THz wave in a nonlinear medium is governed by the wave equation derived from Maxwells equation given as in equation (8). The solution of equation (8.) is of the form [13]

    nonlinear coefficient v were selected from a work by [13], where the allowed values for the nonlinear coefficient, v saturates at 1. The values of o were chosen based on the fact that good dynamic profiles, and very small dispersion and nonlinear effects are mostly observed during short pulse duration [4]. The selection of the order of the solitary wave is also based on the fact that small numerical instabilities are pronounced when N 2 . In this work, the input pulse for all

    the dynamic pulse profiles is the Gaussian shape, which is a component of the input THz pulse generated by [16].


    One should note that pulses emitted from many lasers (where THz radiation is generated) are approximated by Gaussian shapes, but other pulse shapes such as the hyperbolic secant pulse shape that occur naturally in relation to optical solitons could be considered [4]. Figs. 1 4 represent the dynamic pulse evolution of THz radiation at different values of nonlinear coefficient v, dimensionless pulsewidth o, and the order of the solitary wave N

    Fig. 3: The dynamic THz pulse profile for nonlinear coefficient = 2, the dimensionless pulsewidth is of magnitude o = 0.8, and the solitary wave order N = 1.






    Fig. 1: The dynamic THz pulse profile for nonlinear coefficient = 0.004, the dimensionless pulsewidth is of magnitude o = 1 unit, and the solitary wave

    order N = 1.

    Fig. 4: The dynamic THz pulse profile for nonlinear coefficient = 1, the dimensionless pulsewidth is of magnitude o = 0.5, and the solitary wave order N = 1.

    In Figure 1, the dynamic pulse profile is numerically stable throughout the dimensionless distance. This could be as a result of the value of v used, which is within the allowed

    range of values as obtained by [13]. In Figure 2, the dynamic pulse profile is also numerically stable throughout the

    dimensionless distances. The change in the value of the pulsewidth from 1 to 2, and the order of the solitary wave


    -4 -8 from 1 to 2 has not caused any significant change in the pulse


    4 profile. This could be as a result of the value of v used, which is within the allowed range. In Figs. 3 and 4, the dynamic

    Fig. 2: The dynamic THz pulse profile for nonlinear coefficient = 0.004, the dimensionless pulsewidth is of magnitude o = 2 units, and the solitary wave order N = 2.

    pulse profiles produce regular periods. These periods could be due to the increase in the value of the nonlinear coefficient. This periodicity is one of the characteristics of a solitary wave. The periodicity could be desirable if it is not much pronounced. In general, dispersion and nonlinearity act together along the length of the fibre like THz waveguide. Within this length, there is a great deal of smooth temporal THz profile, which is very useful for optical communication systems [11].


In this work, we have generated temporal THz pulses that could be propagated in a nonlinear uniaxial crystal. The results show that stable dynamic THz pulse profiles could be obtained within the allowed range of values of the nonlinear coefficient. The work has also shown that THz radiation could be used in optical telecommunications in a uniaxial crystal fibre like waveguides as a carrier wave, just as lasers in silica optical fibers.


  1. P.H. Siegel, Terahertz technology. Invited paper. IEEE transactions on microwave theory and techniques. .vol. 50. No.3. pp 910 928. 2002

  2. A. M., Leigh; W. Shi,,, J. Zong, Z., Yao, S, Jiang,. and N. Peyghambarian, Narrowband pulsed THz sources using eyesafe region fiber lasers and a nonlinear crystal. IEEE Photonics Technology letters. Vol. 21. No. 1. Pp

    27 29. 2009

  3. S. Gupta,, L. K. Moore,, W. K. Murch, and M. D. Stamper Kurn, Cavity Nonlinear optics at low photon numbers from collective atomic

    motion. arXiv: 0706. 1052V1. pp 1 4. 2007

  4. P.G. Agrawal, Fiber optic communication systems. John Wiley and Sons, Inc. New York. U.S.A. 1992

  5. K.. Xie, Spatial soliton interactions. Unpublished doctoral Thesis

    . University of Salford. U.K. 1993

  6. E.G. Thomas, and A.J. Meadows, Maxwells equation and their applications. Adam Hilger Ltd. Bristol England. 1985

  7. R. Syms, and J. Cozens, Optical guided waves and device. McGraw

    Hill Book Company London. 1993.

  8. J. Belmonto Beitia,, M.V. Perez Garcia, and V. Brazhuyi, Solitary waves in coupled Schrodinger equations with spatially inhomogeneous nonlinearities. arXiv:1001. 2425V1 [nlin – ps]. pp 1 17. 2010.

  9. Z. Wang, Generation of Terahertz Radiation via nonlinear Optical methods. IEEE Transactions on geosciences and remote sensing.

    Vol.1 No. 1. pp 1 5. 2002

  10. W.J. Duffin, Electricity and magnetism. 3rd ed. McGraw- Hill Book Company (UK) Limited. 1981.

  11. P.G. Agrawal, Nonlinear Fiber Optics. Academic Press Inc. Sandiego.

    U.S.A. 1989.

  12. A. Usman., Studies of nonlinear Schrodinger equations for electromagnetic pulse propagation in optical fibers. Unpublished doctoral Thesis, University of Science of Malaysia. 2000.

  13. P. Timtere, A, Usman and. J.C. Ododo. Variational solution for fourth order nonlinear Schrödinger equation for Terahertz (THz) pulses. Nigerian Journal of Physics vol. 21 (2) Pp 23 29. 2010.

  14. P. Timtere, A theoretical study for generation of temporal terahertz (THz) pulses in nonlinear uniaxial crystals. Unpublished doctoral Thesis, Modibbo Adama University of Technology, Yola. Nigeria. 2014

  15. P. Timtere U. Adam, and J. C. Ododo, Evaluation of some nonlinear optical parameters for terahertz (THz) radiating uniaxial crystals. J. Nig. Ass. Math. Physics. Vol. 22. pp 129 142. 2012.

  16. P. Timtere; U. Adam, and J. C. Ododo, Generation of stationary pulse profile for terahertz (THz) radiation in a nonlinear uniaxial medium. IJERT. Vol. 2. Issue 5. Pp 841 848. 2013.

Leave a Reply