E-ISSN:2583-9152

Author Introduction

Photonic materials

Journal of Condensed Matter

2023 Volume 1 Number 2 Jul-Dec
Publisherwww.thecmrs.in

A Unique Approach to Exactly Solve Optical Pulses in Nonlinear Meta-materials

L Nanda1*
DOI:10.61343/jcm.v1i02.44

1* L Nanda, Panchayat College, Bargarh, Odisha 768028, India.

Ultra-short pulse propagation in nonlinear NRM has been investigated where a wide class of solutions for bright and dark solitons has been analyzed for distinct parameter ranges. Here the pulse propagation in nonlinear meta-material has been analytically investigated by solving the nonlinear Schrödinger's equation (NLSE) in composite media expressing frequency dispersion in the dielectric permittivity (ε) and the magnetic permeability (µ). The solutions are exactly shown to be of trigonometric & localized types. The analytical and simulation-based result has been utilized to study the intensity variation in case of a nonlinear meta-material which typically behaves as a negative refractive medium (NRM), for which both ε and µ exhibit frequency dispersion and are negative in nature. The peak of the pulse-intensity curve slowly decreases as the frequency increases towards the magnetic plasma frequency. The stability of the solitonic solutions has also been established.

Keywords: Meta-materials, NRM, NLSE, SRRs, Solitons.

Corresponding Author How to Cite this Article To Browse
L Nanda, Panchayat College, Bargarh, Odisha 768028, India.
Email:
L Nanda, A Unique Approach to Exactly Solve Optical Pulses in Nonlinear Meta-materials. J.Con.Ma. 2023;1(2):151-155.
Available From
https://jcm.thecmrs.in/index.php/j/article/view/44

Manuscript Received Review Round 1 Review Round 2 Review Round 3 Accepted
2023-11-10 2023-11-15 2023-11-20 2023-11-25 2023-12-01
Conflict of Interest Funding Ethical Approval Plagiarism X-checker Note
None Nil Yes 21.98%

© 2023by L Nandaand Published by Condensed Matter Research Society. This is an Open Access article licensed under a Creative Commons Attribution 4.0 International License https://creativecommons.org/licenses/by/4.0/ unported [CC BY 4.0].

Introduction

Nonlinear wave propagation in optics has led to innumerable innovations in several fields. The study of optical solitons in meta-materials has been a new and exciting field of research which opens the door towards numerous scientific theoretical and experimental investigations, the most notable part of which is the study of surface dynamics [1] and fiber optics [2]. In recent years meta-materials [3] have been a subject of immense theoretical and practical interest in infrared and optical frequencies due to their several applications ranging from super-resolution to cloaking. These are artificially structured materials where both the electric and the magnetic responses can be obtained at any desired frequency regime. Most meta-materials exhibit a linear response. However, nonlinearity in meta-materials can be introduced with the addition of a stack of thin wires and split ring resonators inside a nonlinear dielectric. Such type of structured materials has various applications like switching material properties from negative refractive to positive refractive, negative refraction photonic crystals and so on.

For several years significant efforts have been given by several groups to solve the higher order nonlinear Schrodinger's equation (NLSE). There have been a wide class of solutions including bright and dark solitons phase locked with the sources [4] identified for distinct parameter ranges where the solutions can be periodic as well as localized. Further investigations were carried out to study nonlinear propagation of ultra-short pulses in dispersive and negative refractive index media where a wide class of solitary waves has been obtained [5]. In any homogeneous nonlinear meta-materials, the ultra-short electromagnetic pulse propagation is often explained with the help of a system of coupled NLSE which are extracted from the Maxwell's equations. In this article a unique method has been used to exactly solve the NLSE [5] and study the solutions in a medium exhibiting frequency dependent dielectric permittivity and magnetic permeability. The solutions have also been analyzed inside a dispersive negative refractive medium (NRM).

Solution of the Required Nonlinear Schrödinger's Equation 

For ultra-short-wave propagation,

the group velocity dispersion (GVD) is of utmost importance. Hence ignoring the second order temporal derivatives while keeping the linear derivatives up to second order, the Generalized nonlinear Schrödinger equation (NLSE) for relatively transparent materials has been written [5] as,    jcm_44_1.jpg
The Drude-Lorentz model is used to express the dispersion of the dielectric permittivity(ε) and the magnetic permeability(µ) asjcm_44_2.jpg    respectively, where jcm_44_3.jpg   and   Here ωp and ωm  are respectively the plasma frequencies(electric and magnetic respectively) and Λp signifies the electric plasma wavelength. The refractive index is given byjcm_44_4.jpg  The medium dispersion is shown in Fig. 1(a) forjcm_44_4_a.jpg  It allows propagating waves for a positive refractive material (PRM) whenjcm_44_4_b.jpg  for a negative refractive material (NRM) when and jcm_44_4_c.jpgabsorptive elsewhere. The group velocity is represented as jcm_44_4_d.jpgin units of  Eq. (1) can be re-written as


jcm_44_4_e.jpg

Here E(z, t)  denotes the complex electric field with the subscripts  and  corresponding to the partial derivatives for space and time respectively.

The different coefficients in Eq. (3) arejcm_44_5.jpg

where -a1 signifies the GVD. The self-phase modulation (SPM) is given by

jcm_44_6.jpg

Figure 1: (a) The medium dispersion with the refractive index (n(ω))), and real parts of ε(ω) and μ(ω) shown as a function of (b) the different coefficients in the NLSE are shown as a function of ω.
jcm_44_7.jpg
whereas the self-steepening (SS) term is

represented by -a where a =jcm_44_8.jpg
The Fresnel diffraction (FD) coefficient is given by jcm_44_9.jpgwhich mostly defines the pseudo-quintic nonlinearity. Furtherjcm_44_10.jpg  where represents the mixed spatio-temporal second order term. The various coefficients of the NLSE are shown in Fig. 1(b). Here the complex form of the electric field is used which is function of the travelling coordinate jcm_44_11.jpgwhere η signifies the width. The electric field is represented byjcm_44_12.jpg where  stands for the amplitude and  signifies the phase where both are taken as real. Substituting the electric field in Eq. (3), and writing down all the terms together, we get,
jcm_44_15.jpg
where prime stands for the differentiation with respect to ξ. The R.H.S in Eq. (4) is zero means that both components of the complex expression have to separately vanish. They are represented by a pair of coupled equations

jcm_44_16.jpg

where k1 is the constant of integration. Substituting the value of X' from Eq. (7) in Eq. (5) and simplifying


further we get
,jcm_44_17.jpg
After integrating Eq. (8) further, the conserved value of  is obtained which is a constant of integration and is given as
  jcm_44_18.jpg              
Here,jcm_44_19.jpg where is taken as a real function with f(ξ) as a constant parameter. Taking, Eq. (9) is expressed in terms of, which after successive differentiation gives:jcm_44_20.jpg
By making the coefficient of vanish, the above equation converges into a cubic equation which after simplifying givesjcm_44_22.jpgjcm_44_23.jpg

The solution of Eq. (11) is mostly of rational type which is given as f=jcm_44_24.jpg, where there are various possible solutions of, for example f1Cn(ξ,m), with m taken as the modulus parameter. Substituting the expression of  in the cubic equation given by Eq. (11), a polynomial of Cn(ξ,m) is obtained. The consistency conditions can be extracted by putting the coefficients for all terms individually

equal to zero. In this article the solution for m=0 is shown where, Cn(ξ,m=0)=Cos(ξ) Substituting A=0, B≠0 (which are the allowed values) in the consistency conditions, the coefficients obtained are, B=-C3/2η2, D=-2/3 and 2η2=-C1/4 with the consistency condition c12=(-27C2C32/2)1/3 which implies that c1is negative. The obtained solutions of Eq. (11) are periodic and are represented by,jcm_44_25.jpg
Solution in a Meta-material

In the above analysis, the analytical and simulation-based solution of an optical pulse traversing through a nonlinear dispersive medium has been obtained. The results have been analyzed for the case of artificially structured meta-materials which exhibit material dispersion. Further the optical wave propagation has been investigated in the nonlinear NRM which can be achieved within a particular frequency range (Fig.1). The pulse-intensity has been plotted with respect to the scaled frequency () and (ξ) (Fig.2). It is apparent that the peaks corresponding to the pulse-intensity curve decreases with the increase in frequency towards.

jcm_44_26.jpgFigure 2: Pulse-intensity plotted as a function of (ξ) and the scaled frequency () in the negative refractive medium (NRM).

Conclusions

A unique method has been used to obtain an exact analytical and simulation-based solution of an optical pulse traversing through a nonlinear medium having frequency dispersion. It has been shown that the intensity is dependent on frequency when the medium is taken to be dispersive. The solutions have been investigated in an NRM which exhibits dispersive dispersion. The solutions have been


obtained within a particular frequency range. It has been observed that the pulse-intensity decreases with the increase in frequency towards. The stability of the solution is also well established.

Acknowledgments

I heartily thank Prof. S. A. Ramakrishna for his valuable suggestions. Also, my sincere thanks to Prof. P. Panigrahi, and Dr. S. Shaw for helpful discussions. I thank Panchayat College, Bargarh for support

References

1. I. Shadrivov, A. A. Sukhorukov, and Y. S. Kivshar, Phys. Rev. E., 69, 016617 (2004).

2. G. P. Agrawal, Nonlinear Fiber Optics(Academic Press. San Diego, 1995).

3. J. B. Pendry A. J. Jolden, D. J. Robbins, and W. J. Stewart, IEEE Trans. MTT.47, 2075 (1999).

4. V. M. Vyas, P. Patel, P. K. Panigrahi, C. N. Kumar, and W. Greiner, J. Phys. A: Math Gen, 39, 9151 (2006).

5. M. Scalora et. al., Phys. Rev. Lett., 95, 013902 (2005).