Open access peer-reviewed chapter

The Nonlinear Analysis of Chiral Medium

By Andrey Nikolaevich Volobuev

Submitted: June 28th 2018Reviewed: July 9th 2018Published: November 5th 2018

DOI: 10.5772/intechopen.80067

Downloaded: 147

Abstract

The principle of calculation of a plate from a metamaterial with inductive type chiral inclusions is submitted. It is shown that distribution of an electromagnetic wave to such substance can be investigated with the help of introduction of a chiral parameter and on the basis of a detailed method of calculation. By comparison of two methods the dependence of chiral parameter from frequency of electromagnetic radiation falling on a plate is found. With the help of a detailed method the nonlinear equation for potential on the chiral plate is found. It is shown that this equation has solutions as solitary and standing waves but not running waves. The analysis of the received solutions of the nonlinear equation is carried out.

Keywords

  • metamaterial
  • chiral medium
  • chiral parameter
  • nonlinear equation
  • detailed method
  • solitary waves
  • standing waves

1. Introduction

Now the metamaterials (Greek “meta” outside), i.e. composite materials with the various inclusions allocated both chaotically, and periodically are widely applied in particular in a radio engineering, at designing space devices, in medicine, etc. [1, 2, 3]. Due to these inclusions the received materials have many useful physical, electric, optical and other properties which are not present at natural substances. Among metamaterials the substances with chiral properties [4] which capable to rotate a polarization plane of electromagnetic waves are distinguished. In optics as analogue of similar substances are optical active substances, for example, quartz, a solution of glucose etc.

However the methods of metamaterials calculation are enough limited [5]. Basically all calculations are based on the decision of the Maxwell’s equations and the material equations selected according to a problem.

The existing method has restrictions since are usually used only averaged characteristics of metamaterials for example chiral parameter.

In the present work attempt of more detailed approach to properties of the chiral inclusions into metamaterials is made also the analysis of these properties on interaction of chiral elements with the electromagnetic wave falling on a plate from a metamaterial is carried out.

2. Standard method of calculation of metamaterial with electromagnetic wave interaction

At research of metamaterials with chiral inclusions on the basis of Maxwell’s equations usually use the material equations including so-called chiral parameter χ. In [6] the material equations in the following kind are offered:

D=εaEiχVH,E1
B=μaH±iχVE,E2

where D and B there are induction of electric and magnetic fields in the electromagnetic wave propagating in chiral medium, E and H—the strength of electric and magnetic components of wave, εaand μa—absolute electric penetrance and magnetic permeability of chiral medium, V—velocity of an electromagnetic wave in chiral medium, χ—chiral parameter, in this case dimensionless size.

In [6], it is shown that the material Eqs. (1) and (2) can be written down in more simple kind:

D=1χεaE,E3
B=1±χμaH.E4

In formulas (1)(4) top signs concern to right-handed rotation chiral element bottom to left-handed rotation.

Using (3) and (4) it is possible to show [6] that if a chiral medium has only reactance the electromagnetic wave in it submits to the wave equations:

ΔD=1±χV22Dt2,E5
ΔB=1±χV22Bt2,E6

where t there is a time.

Further us the Eq. (5) will interest only. Substituting (3) in (5) and passing to scalar potential φ[7] we shall find:

Δφ=1±χV22φt2.E7

Let us search for the decision of the Eq. (7) as:

φφ0=φrexpiωt,E8

where φ0there is initial a potential level, r—set of spatial coordinates, ω—a cyclic frequency of electromagnetic wave.

Substituting (8) in (7), we have:

Δφr+1±χ2k2φr=0,E9

where k=ωVthere is a module of the wave vector of a falling electromagnetic wave.

Solving Eq. (9) with use of initial and boundary conditions it is possible to investigate processes of reflection, refraction, diffraction of an electromagnetic wave in the metamaterial.

3. Detailed method of calculation of metamaterial with electromagnetic wave interaction

Let us consider a plate of the metamaterial with chiral inclusions of the inductive type. The plate will consist of the dielectric in which are included the current-carrying chiral elements as spirals which axis is directed across a plate.

On Figure 1 the irradiation of a plate by an electromagnetic wave is shown. We assumed as before that chiral inclusions have no active resistance.

Figure 1.

The plate of metamaterials irradiated by an electromagnetic wave.

Feature of a plate is the capacity distributed on its surfaces at dot inductive inclusions. Therefore to examine the interaction of separate chiral element having inductance and capacity with an electromagnetic wave is incorrectly.

At the irradiation on the plate there is a potential difference submitting to the Eq. (7). The density of a current through plate will look like:

jm=Cmφt+φφ0gm,E10

where Сmthere is capacity of the plate area unit, φ—potential on a plate concerning an initial level φ0, gm—electrical conductivity of units of the plate area due to an inductive component.

The first term (10) reflects a capacitor bias current, the second term—an inductive current through the chiral elements.

For spiral chiral element it is possible to write down the equation of a voltage balance:

LiSijit=φφ0,E11

where jithere is density of a current through the i-th chiral element, Li—inductance i-th chiral element, Si—the area of plate, falling one chiral element having inductive electrical conductivity gi.

The density of a current jithrough the chiral element depends on a potential difference on a plate and electrical conductivity it chiral element giunder the formula of the Ohm’s law:

jiSi=giφφ0.E12

Substituting (12) in (11) we shall find:

gi=φφ0Liφt.E13

Electrical conductivity falling unit area of a plate it is equal:

gm=φφ0SiLiφt,E14

where it is taken into account gi=gmSi.

Having substituted (14) in (10) we shall find:

jm=Cmφtφφ02SiLiφt.E15

Using Ci=CmSi—capacity of the plate falling one chiral element, and designating ω02=1CiLiwe shall find:

jmCmφt=φt2φφ02ω02.E16

Let us consider a plate consisting of chiral elements one lines, Figure 2.

Figure 2.

The single-row chiral plate.

Along this plate the inductive current flows.

The law of an electromagnetic induction for this current looks like:

LIXt=φφ0,E17

where IX=γXSlφφ0there is a longitudinal inductive current, γX—specific inductive electrical conductivity a single-row plate, L—its inductance, S—the area of cross-section of a single-row plate, l—its length.

Hence:

γXSL1φt=φφ0,E18

where L1=Llthere is inductance of a single-row plate unit of length.

Under the Ohm’s law for density of a longitudinal current we have:

jX=γXφX.E19

Hence:

djX=γX2φX2dX.E20

Having divided (20) on (18), and having reduced on γXwe shall find:

djX=φφ0SL1φt2φX2dX.E21

On the other hand taking into account that a longitudinal current is defined only presence of a cross-section current (or on the contrary) we have:

SdjX=jmbdX,E22

where b there is width of a single-row plate.

Substituting (21) in (22) we have:

jm=φφ0L1bφt2φX2.E23

Further substituting (23) in (16) we shall find:

φφ0CmL1b2φX2=φt2φφ02ω02.E24

Taking into account С1=Сmb—capacity of a single-row plate unit of length and V2=1C1L1—a square of an electromagnetic field along a plate velocity we have:

V2φφ02φX2=φt2φφ02ω02.E25

The nonlinear Eq. (25) can be transformed to a kind correct for spatial geometry:

V2Δφ+ω02φφ0=1φφ0φt2.E26

Linearization of Eqs. (26) can be carried out by a ratio (8):

Δφr+kS2φr=0,E27

where is kS2=ω02+ω2V2=k02+k2.

4. Various kinds of the equation solution of a metamaterial and an electromagnetic wave interaction

Eqs. (27) and (9) reflect the same physical process – propagation of electromagnetic fluctuations on the chiral plate. Distinction consists that at a deduction (27) as against (9) was not necessity to use the material Eqs. (1)(4) i.e. the chiral parameter was not used.

On the basis of Eqs. (27) and (9) identity it is possible to put down:

kS2=k02+k2=1±χ2k2.E28

Hence the chiral parameter can be written down as:

±χ=1+k02k21.E29

If k0<<kor ω0<<ωformula (29) becomes simpler:

χ=±k022k2=±ω022ω2.E30

Let us notice that quantum calculation of an optical active substance [8, 9] results in the formula for chiral parameter:

χ=23ћω0jω0j2ω2,E31

where ћthere is Planck’s reduced constant, δ- size proportional to product of the real parts electric and magnetic dipole moments of an optical active molecule power transition excited by a light of the wave given length, ω0j—in this case the frequency corresponding to power transition 0j[10].

The increase in a degree of frequency dependence ω0up to square-law in formula (30) in comparison with (31) is characteristic at transition from quantum area in classical.

On Figure 2 the illustrative graph of the potential fluctuations on the chiral plate is shown according to the oscillatory solutions satisfying Eqs. (9) and (27). Character of fluctuations will be investigated below.

The nonlinear Eq. (25) has at least one more solution as a solitary wave:

φφ0=φmaxexpk0XX0±ω0tt022.E32

whereas before k0=ω0V, φmaxthere is a peak value of potential, t0—an initial of a time reading, X0—coordinate of the chiral element center. The sign a minus concerns to a wave spreading from left to right, a sign plus from right to left.

The illustrative graph of the solution (32) as function of coordinate X also is shown on Figure 2. Growth of potential above chiral inclusions is caused by proportionality of the chiral inclusions reactance their inductivities φφ0XLi=ωLi.

From the analysis of both curves it is possible to conclude that the top curve, Figure 2, concern to often enough inclusions of the chiral elements in a plate, and bottom to more rare. Therefore into the solution (32) to enter a chiral parameter it is irrational.

Obviously for the nonlinear Eqs. (25) or (26) there should be a multiwave solution. However to find such solution it is extremely difficult. Multiwave solutions are found for very much limited circle of the nonlinear wave equations [11, 12]. The multiwave solution should depend on concentration of the chiral elements in a plate. Only with its help it is possible to understand under what conditions it is possible is proved to enter the chiral parameter, i.e. to understand borders of the material Eqs. (1)(4) applicability.

Let us consider in more detail a kind of the wave arising on single-row chiral plate at falling on it of an electromagnetic wave.

The nonlinear Eqs. (25) and (26) can be solved a method of variables division [13]. We search for the solution of Eq. (25) as:

φφ0=φXTt.E33

where φXthere is function only coordinates X, Tt—a function only time t.

Having substituted (33) in (25) we shall find:

V2φXT2t2φXX2=φXTtt2φ2XT2tω02.E34

Let us divide both parts of the equation on φ2XT2t. In result we shall receive:

V21φX2φXX2+ω02=1TtTtt2=α2.E35

where αthere is a constant.

Eq. (35) breaks up to two independent equations. The equation dependent on X looks like:

2φXX2+k02+α2V2φX=0.E36

Comparing (36) and (27) we notice that kS2=k02+α2V2. Hence k2=α2V2, and hence α=ω.

The solution of Eq. (36) we shall write down as:

φX=φ0expikSX.E37

where φ0there is value of function φXin the beginning of coordinates.

The second equation of equality (35) looks like:

Ttt=iωTt.E38

Solving this equation we shall find:

Tt=T0expiωt,E39

where T0there is initial value of function Tt.

Using (33) and the real parts of solutions (37), (39) we shall find the solution of Eq. (25):

φφ0=φАexpiωtexpikSX=φАcosωtcoskSX=φАcosωtcos2πXλ.E40

where it is designated φA=T0φ0—a peak value of potential on a plate, λ—a wave length.

Formula (40) represents the equation of a standing wave.

Condition of the nodes occurrence in a standing wave Xns=±2n+1λ4, where n=0,1,2,.

On the ends of the single-row chira plate, Figure 2, should be nodes of a standing wave. If excitation of a wave occurs in the center of a plate the number of the maximal node can be found under the formula ±l2=±2nmax+1λ4or nmax=lλ12.

It is necessary to note that running waves φφ0=φА2coskSX±ωtare not the solution of Eq. (25) therefore formula (40) from the physical point of view cannot be presented as the sum of the direct, and reflected from borders plate waves though mathematical this procedure is simple for making. It is consequence of Eq. (25) nonlinearity.

5. Conclusion

Distribution of potential to a plate from a metamaterial with inductive chiral inclusions is investigated as with use of the material equations together with the Maxwell’s equations, and on the basis of a detailed method of calculation of the chiral elements and an electromagnetic wave interaction. Comparison of two approaches has allowed to find out that introduction of the chiral parameter is correct only at enough high concentration of the chiral inclusions. On the basis of comparison of two methods results the frequency dependence of chiral parameter is found. At use of a detailed method of calculation the nonlinear equation for the potential having solutions as standing waves and solitary waves is received. Running waves are not the solution of this equation. Necessity of the multiwave solution existence of the nonlinear equation which should depend on concentration of the chiral elements in a metamaterial is marked.

© 2018 The Author(s). Licensee IntechOpen. This chapter is distributed under the terms of the Creative Commons Attribution 3.0 License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

How to cite and reference

Link to this chapter Copy to clipboard

Cite this chapter Copy to clipboard

Andrey Nikolaevich Volobuev (November 5th 2018). The Nonlinear Analysis of Chiral Medium, Chirality from Molecular Electronic States, Takashiro Akitsu, IntechOpen, DOI: 10.5772/intechopen.80067. Available from:

chapter statistics

147total chapter downloads

More statistics for editors and authors

Login to your personal dashboard for more detailed statistics on your publications.

Access personal reporting

Related Content

This Book

Next chapter

Chirality Properties of Modeling Water in Different Aqueous Systems

By Khakhalin Andrey Vladimirovich and Gradoboeva Olga Nikolaevna

Related Book

First chapter

Introductory Chapter: Crystallography

By Takashiro Akitsu

We are IntechOpen, the world's leading publisher of Open Access books. Built by scientists, for scientists. Our readership spans scientists, professors, researchers, librarians, and students, as well as business professionals. We share our knowledge and peer-reveiwed research papers with libraries, scientific and engineering societies, and also work with corporate R&D departments and government entities.

More About Us