Open access peer-reviewed chapter

Chiral Waves in Graphene Medium and Optical Simulation with Metamaterial

By H. Torres-Silva

Submitted: February 24th 2012Reviewed: July 4th 2012Published: November 14th 2012

DOI: 10.5772/51328

Downloaded: 1252

1. Introduction

In modern age Engineers have paved the way for a new generation of faster, more powerful cell phones, computers and other electronics by developing a practical technique to replace silicon with carbon on large surface. The capability of silicon, the material at the heart of computer chips has been harnessed beyond its limits by engineers and carbon has come up as an integrating replacement for the same. The material called“Graphene” which is a single layer of atoms arranged in honeycomb lattice could let electronics to process information and produce radio transmission 10 times better than silicon based devices.

For theorists, such a system is also of great interest because it provides a physical realization of two-dimensional field theories with quantum anomalies. Indeed, the continuum limit of the effective theory describing the electronic transport in graphene is that of two-dimensional massless Dirac fermions. The reported and predicted phenomena include the Klein paradox (the perfect transmission of relativistic particles through high and wide potential barriers), the anomalous quantum Hall effect induced by Berry phases and its corresponding modified Landau levels and the experimental observation of a minimal conductivity.

From the point of view of its electronic properties, graphene is a two-dimensional zero-gap semiconductor with the cone energy spectrum, and its low-energy quasiparticles are formally described by the Dirac-like Hamiltonian [1, 2].


Where vF106ms1is the Fermi velocity and σ=(σx,σy)are the Pauli matrices. The fact that charge carriers are described by the Dirac-like equation, rather than the usual Schrödinger equation, can be seen as a consequence of graphene’s crystal structure, which consists of two equivalent carbon sublattices [1, 2]. Quantum mechanical hopping between the sublattices leads to the formation of two cosine-like energy bands, and their intersection near the edges of the Brillouin zone yields the conical energy spectrum. As a result, quasiparticles in graphene exhibit the linear dispersion relationEG=E=kvF, as if they were massless relativistic particles with momentum k (for example, photons) but the role of the speed of light is played here by the Fermi velocityvFc/300. Owing to the linear spectrum, it is expected that graphene’s quasiparticles will behave differently from those in conventional metals and semiconductors where the energy spectrum can be approximated by a parabolic (free-electron-like) dispersion relation.

From a crystallographic point of view, the graphene is a triangular Bravais lattice with a diamond-shaped unit tile consisting of two sites so one gets the honeycomb structure. The very unique feature of the graphene band structure is that the two lowest-energy bands, known as the valence and the conduction bands, touch at two isolated points located at the corners of the Brillouin zone. In the immediate vicinity of these degeneracy points, known as the Dirac points, the band structure is a cone. In natural graphene samples, there is exactly one electron per site, and thus, at zero temperature, all levels in the valence band are filled (a situation known as half-filling). As a result, the energy of the last occupied level precisely slices the band structure at the Dirac points. The low-energy excitations of this system are then described by the massless two-dimensional Weyl-Dirac equation and their energy dispersion relationω=vFkis that of relativistic massless fermions with particle-hole symmetry. In graphene these massless fermions propagate with a velocityvF. The maths is simple but the principles are deep. We will review the formulation of graphene’s massless Dirac Hamiltonian, under the chiral electromagnetism approach, like a metamaterial media, hopefully demystifying the material’s unusual chiral, relativistic, effective theory. The novel result here is that in our theory we do not make cvF[1, 2, 27-35], but we obtain vFas vF=c(1k0T)ifk0T>0orvF=c(1+k0T)ifk0T<0. These results are derived of the Chiral Electrodynamics with T as the chiral parameter andk0=ω/cwhen the electric wave Eis quasi parallel to the magnetic waveH.

With this approach we can review different phenomena such: The Dirac point with a double-cone structure for optical fields, an optical analogy with Dirac fermions in graphene, can be realized in optically homogenous metamaterials. The condition for the realization of Dirac point in optical systems is the varying of refractive index from negative to zero and then to positive.

Also we give a support to the similitude of the band structure of a macroscopic photonic crystal with the electronic band structure of grapheme, which is experimentally much more difficult to access, allows the experimental study of various relativistic phenomena. With our analytical and numerical analysis we hope to verify that, similar behaviors exist to electrons in graphene treated as mass-less particle, ie, electron wave propagation.

This chapter presents a short review on the chiral propagation of electron waves in monolayer graphene and optical simulation with optical field in the negative-zero-positive index metamaterial NZPIM and its close connection. Section II presents an enhanced vector diagram of Maxwell’s equations for chiral media with quasi parallel electromagnetic fields,EH. Chiral waves in graphene acting as metamaterial media are discussed in section III. In section IV, two component equations and tunneling rate of Dirac electron in graphene are derived. Section V is about Zitterbewegung of optical pulses near the Dirac point inside a negative-zero-positive index metamaterial, and shows that the chiral field near the Dirac point becomes a diffusive wave The last sections describe the theoretical description of minimal conductivity in graphene under chiral approach and absorption of light by quasi 2D Dirac fermions.

2. A vector diagram of Maxwell’s equations for chiral media with EH

The idea of representing Maxwell’s time-harmonic equations in homogeneous isotropic media by vector diagram as put forward by Wilton [3] and by S. Uckun [4] deserves consideration. All the common relations between field and potential quantities implied by Maxwell’s equations can be represented by a diagram. It is started that the diagram not only illustrates Maxwell’s equations, but also many of the methods for constructing diagram are based on the formal similarity between many theorems of vector calculus and those of vector algebra.

An isotropic chiral medium is a macroscopically continuous medium composed of equivalent chiral objects that are uniformly distributed and randomly oriented. A chiral object is a three-dimensional body that cannot be brought into agreement with its mirror image by translation and rotation. An object of this sort has the property of handedness and must be either left-handed or right-handed. An object that is not chiral is said to be achiral, and thus all objects are either chiral o achiral. Due to their novel properties and wide applications in microwave and radar engineering, chiral media has been undergoing extensive research during the last years. That is why this study aims to cover chiral medium for the representation of Maxwell’s equations in vector diagram form. In a chiral media a cross coupling between electric and magnetic filed exists. Thus, the vector diagram has vectors along all three coordinate axes whereas the vector diagram presented by Wilton [3] for achiral media has vectors only in one plane with H vector normal to it.

2.1. Vector diagram construction

Assumingexp(iωt)time dependence, Maxwell’s time-harmonic equations [5, 6] for isotropic, homogeneous, linear media are:


Chirality is introduced into the theory by defining the following constitutive relations to describe the isotropic chiral medium [5]


Where the chirality admittance ωεTindicates the degree of chirality of the medium, and the εandμare permittivity and permeability of the chiral medium, respectively. Since Dand Eare polar vectors and Band Hare axial vectors, it follows that εand μare true scalars and ωεTis a pseudoscalar. This means that when the axes of a right-handed Cartesian coordinate system are reversed to form a left-handed Cartesian coordinate system, ωεTchanges in sign whereas εand μremain unchanged.

For a graphical representation of the above relationships, following Wilton’s procedure [3], let us assume vector-differential operator,is an ordinary vector and treat the divergence and curl operations in equations (1) to (4) as ordinary scalar (dot) and vector (cross) products, respectively. Equation (3) implies that is perpendicular to Band the vector ×Bmust be perpendicular to both andB.

As shown in Figure 1, three transverse coordinate axes are chosen as, from Equation (1) B=×E/(jω)and, ×B/(jωεμ')=××E/k2wherek2=ω2εμ'.

Since B=0always, this conditions will hold identically if Bis expressed as the curl of a vector potential Asince the divergence of the curl of a vector is identically zero, thus

andAmust be perpendicular to both and Band lie in and ×Bplane. However, Ais not unique since only its components perpendicular to contribute to the cross product. Therefore,A, the component of Aparallel to, must be specified. The curl equation forE, as in Equation (1), and Equation (7) give ×(E+jωA)=0where the quantity in parentheses should be parallel to and the curl of the gradient of a scalar function ϕis identically zero; so the above equation is E+jωA=ϕor

That is shown in Figure 1.

Figure 1.

Diagram of the full Maxwell System for a chiral media with Lorentz gauge μ ε 1 − k o 2 T 2 = μ ' ε ' .

Following the Uckun’s approach [5], we substitute Equation (7) into Equation (6) having


Substituting Equation (9) and Equation (5) into Equation (2) gives ××A+jωμεT1ko2T2×E=jωμε1ko2T2Eω2μεT1ko2T2B+μ1ko2T2J

××A+jωμεT1ko2T2×E=jωμε1ko2T2Eω2μεT1ko2T2B+μ1ko2T2J, playing the value of×Efrom Equation (1) into the above equation ××A+2μω2εT1ko2T2B=jωμε1ko2T2E+μ1ko2T2Jand using the vector identity ××A=(A)2Aenables us to write the above equation as 2A+2ω2μεT1ko2T2×A=(·A)jωμε1ko2T2Eμ1ko2T2J2A+2ω2μεT1ko2T2×A=(A)jωμε1ko2T2Eμ1ko2T2J, so using Equation (8), we have

Here Ais arbitrary, so in order to specifyA, for uniqueA, we may choose


And eliminate the term in parentheses in Equation (10). The choice in Equation (11) can be known as a chiral Lorentz gauge. Then Equation (10) will be simplified to


Divide both sides of Equation (12) by (jωμε1ko2T2)and reorganize it to get


( Figure 1, shows this vectorial equation).

The difference between our approach and the Uckun’s procedure [4], is that we take the chiral media characterized by D=ε(E+T×E)andB=μ(H+T×H). In this form we can obtain the condition of quasi spatial parallel condition between Band Ewhere the main equation is like a Beltrami equation which is important for the numerical simulation of graphene systems.

Placing the value of Bfrom Equation (1) into Equation (5) D=ε(E+T×E)will be obtained. Placing the value ofH, from Equation (6), and Dinto Equation (2) will give


by rearranging this equation we obtain


will be obtained as shown in Figure 1. In this figure we putμμ/(1k02T2).

Taking divergence of Equation (5) and using Equations (3) and (4) in it


will be derived. To find the projection of Eonto, from Equation(15) E=ρ/ε, take the gradient of both sides and divide by scalar value 2to normalize to a unit vector. So


Similarly, getting gradient of both sides of Equation(11), using the vector identify ××A+2Afor (A)and normalizing by 2we have:


So it will be obtained as parallel component of A to coordinate.

By using Equations (8), (13), (14), (16), and (17) the vector diagram of Lorenz gauge can be completed as shown in Figure 1, where all Maxwell’s relations and potential quantities appear.

Now let us examine derivation of some relations from the diagram. For example, it is seen that the component of Eand J/(jωε)parallel to must be equal and opposite. By taking the divergence of Equation (14) and using E=ρ/εcan be shown that


Taking the gradient in both sides of Equation (18) and dividing it by scalar value 2will give the same value as Equation (16) with opposite sign. From the right side of Equation (18) it is seen that


This is the known continuity equation. Since the divergence of the curl of any vector is identically zero, the divergence of Equation (2) yields.0=jωD+J. Using Equation (4) convert this immediately into continuity equation as, expected. Again, as seen in Figure 1, 2T×Eand jω2T×Αare equal and opposite vectors. From Equation(8), taking curl of both side and using the vector identity ×ϕ=0will show that


as expected. By using the vector calculus a few possible equations from the vector diagram can be written as follows


For example, adding Equations (13) and (14) side by side and using Equation (8) will give Equation (24) which shows the correctness of the equation derived from the diagram 1. Instead of Lorenz gauge we can choose Coulomb’s gauge.


In Equation (10) so that it will take the form


where the subscript "c"is used it indicate Coulomb’s gauge. Using Equation (8) and (20) in the above equation


will be obtained. Placing the values of Equations (5) and (6) into Equation (2) will give


and value of ×Efrom Equation (1) will give


Combining these equations with Equation (26) and using Equation (7) we have


By using the same coordinates axes, Band (1ko2T2)×B/(jωμε)and Equation (1), (8), (16), (26) and (27) for the Coulomb gauge. It is clear from Equation (25) that the component of the vector A parallel tois equal to zero.

As seen in Figure 1, Lorenz gauge are the best choice because these make Aeither parallel or perpendicular to any of the other vectors and simplify its relationship to those vectors. In Figure 1a, if the chirality factorT, goes to zero, pointK, Land Rapproach pointM, P, and Nrespectively, in which case the diagram will be the same as in Reference [3] for linear, homogeneous, isotropic achiral medium. If (1ko2T2)0then E is quasi parallel to B, and parallel toAso all vectors remain in an only plane. (see Figure 2).

Figure 2.

Two dimensional condition when ( 1 − k o 2 T 2 ) → 0 for a graphene system, the thickness is about one carbon atom

This Beltrami condition is useful to numerical calculations in graphene. We apply this approach to a two dimensional chiral graphene slab. This result cannot be obtained with the Uckun’s approach [4]. In terms of chiral magnetic potentialAcwithJ=0, andρ=0we have the wave equation

Ask=kx2+ky2+kz2, ifkx=ksinθ, ky=0, andkz=kcosθ, we have the matrix:

The dispersion relation of the transversal wave is (k2(1k02T2)+k02)24k04k2T2(sin2θ+cos2θ)=0k=k±=±k0/(1±k0T).

That is


So we have k±=±ωc(1±k0T)1=ω/vFwith vF=c(1±k0T)if we consider a bi dimensional graphene system.

The novel result here is that in our chiral theory we do not make cvFbut we obtain vFas vF=c(1k0T)ifk0T>0orvF=c(1+k0T)ifk0T<0. These results are derived of the Chiral Electrodynamics with T as the chiral parameter andk0=ω/c[5, 6, 7]. This situation corresponds to a model of chiral electron where the Efield is almost parallel to theHfield, the movement of electrons is helical in the valence band of the membrane graphene. Semi classically, for a hydrogen atom, this condition was shown by Huang [36] for Ewave perpendicular to Hwave with radiation and by Torres-Silva for Eparallel toH, without radiation. See refs ([5, 37, 38]). Here, we extend this result to the electron bound to carbon atom.

In the next section we study the situation when the refractive index is negative.

3. Chiral waves in graphene acting as metamaterial media

Metamaterials are composite materials in which both permittivity and permeability possess negative values at some frequencies has recently gained considerable attention [see e.g., [8-12]. This idea was originally initiated by Veselago in 1967, who theoretically studied plane wave propagation in a material whose permittivity and permeability were assumed to be simultaneously negative [11]. Recently Shelby, Smith, and Schultz constructed such a composite medium for the microwave regime, and experimentally showed the presence of anomalous refraction in this medium [10]. Previous theoretical study of electromagnetic wave interaction with omega media using the circuit-model approach had also revealed the possibility of having negative permittivity and permeability in omega media for certain range of frequencies [9].

The anomalous refraction at the boundary of such a medium with a conventional medium, and the fact that for a time-harmonic monochromatic plane wave the direction of the Poynting vector is antiparallel with the direction of phase velocity, can lead to exciting features that can be advantageous in design of novel devices and components. For instance, as a potential application of this material, compact cavity resonators in which a combination of a slab of conventional material and a slab of metamaterial with negative permittivity and permeability are possible. The problems of radiation, scattering, and guidance of electromagnetic waves in metamaterials with negative permittivity and permeability, and in media in which the combined paired layers of such media together with the conventional media are present, can possess very interesting features leading to various ideas for future potential applications such as phase conjugators, unconventional guided-wave structures, compact thin cavities, thin absorbing layers, high-impedance surfaces, to name a few. In this section, we will first present a brief overview of electromagnetic properties of the media with negative permittivity and permeability, and we will then discuss some ideas for potential applications of these materials.

Such a medium is therefore termed left-handed medium [12]. In addition to this ‘‘left-handed’’ characteristic, there are a number of other dramatically different propagation characteristics stemming from a simultaneous change of the signs of εandμ, including reversal of both the Doppler shift and the Cerenkov radiation, anomalous refraction, and even reversal of radiation pressure to radiation tension. This configuration exhibit also chirality and a rotation of the polarization so the analysis of metamaterial presented by several authors provides a good but not exact characterization of the metamaterial [9]. The evidence of chirality behavior suggests that if it is included in the conditions to obtain a metamaterial behavior of a medium futher progress will be obtained. In this short paper, we propose to investigate the conditions to obtain a metamaterial having simultaneously negative εand negative μand very low eddy current loss. As a initial point, we consider a media where the electric polarization depends not only on the electric fieldE, and the magnetization depends not only on the magnetic fieldH, and we may have, for example, constitutive relations given by the Born-Federov formalism [12].


The pseudoscalar Trepresents the chirality of the material and it has length units. In the limitT0, the constitutive relations (31) and (32) for a standard linear isotropic lossless dielectric with permittivity εand permeability μare recovered.

According to Maxwell’s equations, electromagnetic waves propagating in a homogeneous dielectric magnetic material are either positive or negative transverse circularly polarized waves, and can be expressed as


whereE0±=E0(x^±iy^), and×E±(r,t)=k±E±, k±0is the chiral wave number.

If the phase velocity and energy flow are in the same directions, and from Maxwell’s equation, one can see that the electricEand magnetic field and Hand the wave vector kwill form a right-handed triplet of vectors. This is the usual case for right-handed materials. In contrast, if the phase velocity and energy flow are in opposite directions, andE, H, and kwill form a left-handed triplet of vectors. This is just the peculiar case for left handed materials where the effective permittivity ε(1+T×)and the effective permeability μ(1+T×)are simultaneously negative. So, for incident waves of a given frequencyω, we can determine whether wave propagation in the composite is right handed or left handed through the relative sign changes ofk.

To advance in our propose we considerer other more popular representation to describe a chiral medium, [12] as


in which electromagnetic coupling terms are added to the basic terms. Bi-isotropy or bianisotropy is used for calling such constitutive equations, according to the parameters to be scalars or tensors. If κ=0andχ0, it is the Tellegen medium; if κ0andχ=0, as the requirement of reciprocity, it is the Pasteur medium:


There is a long dispute on strong chiral medium since it was introduced theoretically. Traditional electromagnetic conclusions have limited us to understand strong chirality, i.e. k+T1orμPεPκ2, until we see the fact that artificial Veselago’s medium [9] was successfully realized in certain frequency bands [3]. For the case, EH, we have to ask the following question: can strong chiral medium exist?

In Ref. [14], the reason for traditional restriction of chirality parameters was concluded

as: 1) The wave vector of one eigenwave will be negative; 2) The requirement of a positive definite matrix to keep positive energy:


With the exploration of backward-wave medium, we know that negative wave vector, or opposite phase and group velocities, are actually realizable. And there is an unfortunate mathematical error in the second reason: in linear algebra, only if it is real and symmetric, positive definite matrix is equivalent to that all eigenvalues should be positive. The matrix (14) is a complex one, making the analysis on restriction of positive energy meaningless.

Actually, in a strong bi-isotropic medium with constitutive relations as Eqs. (1) and (2), the energy can be drawn as


Mathematically, the amount of energy density propagated is proportional to the magnitude of the Poynting vectorS, whereS=E×H. With the conditionEH, S(r,t)=0so we find right circularly polarized wave or left circularly polarized wave.

The concept of parallel fields is important in the theoretical formulation of: Space electromagnetism and vacuum, the classical and quantum gravitational fields, the study of elementary particles, operator and Dirac matrices, fields and chiral electrodynamics [15, 16]. If we putE=αBand solving the Maxwell equations with the Born Fedorov relations we findα=i/με/4E=iηH=iμ/2ε/2Hsoε(ω), μ(ω)and T are transformed as:


In this case the total density energy is may be


In this special case where the energy propagated in one direction is equal to that propagated in the opposite direction, there is no net energy flow in the medium and the sum of the two TEM waves form what is generally known as a standing wave. The condition for a standing wave is that the time average of Svanishes. This can be achieved if Sis zero all the time everywhere in the region of space under consideration, i.e., S(r,t)=0(see Fig. 4). Examination of Eq. (39) shows thatS=0ifEH. In this last case, a particular solution of Eq. (31,3 2) is whenk02T2=1, where we have the conditionEH, andE=iηH, so we find the Beltrami force free equation E+2T×E=0and the vector Poynting S(r,t)vanishes [15-16].

In terms of Eqs. (41-42), ifμPεP<κ2for whole frequency range, the energy will still keep positive as long as the permittivity and permeability are positive, under the weak spatial dispersion condition. Therefore, the real reason for traditional strong-chirality limitation is neither negative wave vector nor energy conversation. Next we will point out two other important reasons.

First, with the assumption thatεP>0, μP>0, κ>0andκ>μPεP, we easily show thatε,μand T turn to negative from the transformation between Pasteur constitutive relations and BF relations shown in Eqs. (41)-(43). This is absolutely unacceptable before people realizing Veselago’s medium. Actually, strong chiral medium can be equivalent to Veselago’s medium for the right circularly polarized wave [9-16]. The negativeεand μhave shown such a point. Hence the negative sign in the BF model is not strange at all, since we realize effective double-negative with strong chirality parameter instead of simultaneously frequency resonances. For a limiting case, the chiral nihility [10], in which μP0and εP0whileκ0, the parameters in DBF representation becomeε,μandT=1/ωκ, remaining a finite value after a simple mathematical analysis. There is no evidence that strong chirality cannot exist in this aspect.

Second, it is the effectiveness of linear models. Similar to the case that linear optical and electromagnetic models can no longer deal with very strong optical intensity and electromagnetic field, we introduce nonlinear optics to take into account the higher order terms of polarization.

If the spatial dispersion is strong enough, the higher order coupling terms cannot be neglected as before. People used to mistake strong chirality with strong spatial dispersion, hence adding a limitation to chirality parameter,κμPεP. We believe that this is the most probable reason.

However, the strong spatial dispersion is embodied in the BF model, e.g. the value of T, while the strong chirality is represented by the Pasteur model, e.g. the ratio of κtoμPεP. That is to say, strong chirality does not necessarily lead to strong spatial dispersion.

Based on Eqs. (41)-(43), we have computed T and ε/εPor μ/μPversusκ/μPεP, as shown in Figs. 2 and 3. When κis very close toμPεP, the value of T is quite large, indicating a strong spatial dispersion. Hence the singular point is the very point of traditional limitation. However, withκcontinuously increasing, the spatial dispersion strength falls down very quickly. Therefore, if κis not aroundμPεP, e.g. κ<0.8μPεPorκ>1.2μPεP, we need not take nonlinear terms into consideration at all. Hence the strong spatial dispersion and nonlinearity cannot put the upper limitation to chirality parameters either.

Figure 3.

The strength relationship of chirality and spatial dispersion. ω T / c versus κ / μ P ε P The point of κ / μ P ε P = 1 is singularity, corresponding infinite spatial dispersion coefficient T. When κ / μ P ε P 1 , T becomes negative for keeping the positive rotation term coefficients with negative μ ∥ and ε ∥ . The 2-D grapheme system may be modeled when ω T / c ∼ 0.9 − 1.1 and κ / μ P ε P ∼ 1

Figure 4.

ε ∥ / ε p , μ ∥ / μ p versus κ / μ P ε P . With chirality strength increases, ε ∥ and μ ∥ reduces quickly from ε p and μ p to − ∞ . The case 2-D graphene corresponds to κ / μ P ε P ∼ 1 .

If ω|T|/c=0.9927orω|T|/c=1.0073, we have that the Fermi velocity is given byvF=c(1±k0T), so it appears a connection between the chiral electrodynamics and αthe fine structure constant.

From Figure 4. we see that near of ε/εp,μ/μp0the graphene have a linear behavior so the Beltrami equation is


Solving Beltrami's equation with boundary conditions, we can demonstrate the optical Zitterbewegung effect by means of electromagnetic pulses propagating through a negative-zero-positive index metamaterial (NZPIM) if we make


Thus in optics, the Beltrami’s equation for electromagnetic waves can be reduced to the Helmholtz equation.

In later sections we discuss some effects such as Zitterbewegung of optical pulses, diffusion phenomenon and tunneling rate of dirac electron in graphene.

4. Two component equations and tunneling rate of dirac electron in graphene

The usual choice of an orthogonal set of four plane-wave solutions of the free-particle Dirac equation does not lend itself readily to direct and complete physical interpretation except in low energy approximation. A different choice of solutions can be made which yields a direct physical interpretation at all energies. Besides the separation of positive and negative energy states there is a further separation of states for which the spin is respectively parallel or antiparallel to the direction of the momentum vector. This can be obtained from the Maxwell’s equation without charges and current in the EHconfiguration. Dirac's four-component equation for the relativistic electron is [5-7]. (see our chapter “Chiral Transverse Electromagnetic standing waves with EHin the Dirac Equation and the spectra of the Hydrogen Atom”, section 2,Intech 2011, Behaviour of electromagnetic waves in different media and structures, edited by Ali Akdagli). Here we considerer a bidimensional graphene system so the Dirac's four-component equation for the relativistic electron is:




and I is the two-by-two identity matrix and the Fermi velocity vFis deduced from the chiral electrodynamics withvF=c(1±k0T), where T is the chiral parameter in a metamaterial condition. This result is capital to our approach because we fin a contact point between the graphene system and optical metamaterial makingvF=c(1±k0T), no making cvFas other authors do it. In Figure 5 we have1.0073>k0T>0.9927.

The Hamiltonian commutes with the momentum vectorp^. In order to resolve this degeneracy we seek a dynamical variable which commutes with both H andp^. Such a variable isσ^p^, whereσ^is the matrix Pauli. The eigenfunctions of the commuting variables p andσ^p^are simultaneous:


Thus for a simultaneous eigenstate of p^andσ^p^, the value of σ^p^will be +p or –p, corresponding to states for which the spin is parallel or antiparallel, respectively, to the momentum vector like a graphene system.

A simultaneous eigenfunction of H and p will have the form of a plane wave


where the ψjare the four components of the state function and ujfour numbers to be determined. Then E can have either of the two values.


We now demand that ψjbe also an eigenfunction of σ^p^belonging to one of the eigenvaluespE, say, wherepE=±p.

The eigenvalue equation is


Since Wcan be given either of the two values ±εandpE, the two values±p, we have found for given p four linearly independent plane wave solutions. It is easily verified that they are mutually orthogonal.

Figure 5.

Diagram where the graphene is located in the plane (x, y), the chiral photon is in either direction. The pontos black figure on the left indicates the position of carbon atoms.

The physical interpretation of the solutions is now clear. Each solution represents a homogeneous beam of particles of definite momentum p, of definite energy, either±ε, and with the spin polarized either parallel or antiparallel to the direction of propagation. From here we can obtain the well known equation for grapheneH0=ivFσ.

In this section, we study the tunneling rate of Dirac electrons in graphene through a barrier with an intense electromagnetic field. A one transport phenomenon in graphene is the chiral tunneling [1, 2, 17, 18]. In mono layer graphene a perfect transmission through a potential barrier in the normal direction is expected. This tunneling effect is due to the chirality of the Dirac electrons, which prevents backscattering in general. This kind of reflectionless transmission is independent of the strength of the potential, which limits the development of graphene-based field-effect transistors (FET). The perfect transmission can be suppressed effectively when the chiral symmetry of the Dirac electrons is broken by a laser field, when the n-p junctions in graphene are irradiated by an electromagnetic field in the resonant condition [19,20].We consider a rectangular potential barrier with heightH0, width D in the X direction, and infinite length in the Y direction. The Fermi level lies in the valence band in the barrier region and in the conduction band outside the barrier. The gray filled areas indicate the occupied states. The optical field under chiral condition is propagated perpendicular to the layer surface and it is circularly polarized along the Z,Y direction with a detuningΔ0=2Ebω. We choose Δ0>0to ensure that there is no inter-band absorption inside the barrier. Meanwhile, 2Ekωis used to guarantee that the influence of the optical field outside the barrier can be neglected. Thus, neglecting the scattering between different valleys, the scattering process of Dirac electrons in K point is described by the time-dependent Dirac equation, Eq. (45) withm=0, so ψ=ψg=(ψE,ψH)twhere t means transposed wave function.


where ψg(r,t)=[ϕA(r,t),ϕB(r,t)]tis the wave function, He=vFσ^p=ivFσ^is the unperturbed Dirac Hamiltonian obtained from the chiral electrodynamics [17], σ^=(σxσy)are the Pauli matrices, vF106ms1is the Fermi velocity, H0(r)is the height of the potential barrier, I is the unit matrix, and Hintis the interaction Hamiltonian with a chiral electromagnetic potential.


where e is the electron charge and the chiral potential vector is [A]=[Axeiωt,iAyeiωt]with i=1and |Ax|=|Ay|so we have circular polarized wave of the electromagnetic field. When the Dirac electrons perpendicularly collide with the barrier perpendicularly, we can rewrite Eq. (53) as a set of partial differential equations


Since the tunneling time is of order of sub-picosecond and the potential H12(t)and H21(t)vary as fast as the frequency of incident light beams, this scattering process is strongly time-dependent.

In order to study such a strongly time-dependent scattering process, we employ the finite- difference time-domain (FDTD) method to solve Eq. (47) and Eq. (48) numerically in the time-domain [21]. In the traditional FDTD method, the Maxwell’s equations are discretized by using central-difference approximations of the space and time partial derivatives. As a time-domain technique, the FDTD method can demonstrate the propagation of electromagnetic fields through a model in real time. Similar to the discretization of Maxwell’s equations in FDTD, we denote a grid point of the space and time as (i,k)=(iΔx,kΔt)[22]. For the nodal variables we can apply the usual Lee discretization method in 2-D systems withk02T21(See Figure 6).

Figure 6.

Lee discretization. Here H x = ( ∇ × A → ) x E x = − ( ∂ ( 1 + T ∇ × ) A → / ∂ t ) x and with J → = 0 (See [21, 22, 24]).

For any function of space and timeG(iΔx,kΔt)=Gk(i), the first order in time or space partial differential can be expressed as:


These eqs. can be replaced by a finite set of finite differential equations like:


For computational stability, the space increment and the time increment need to satisfy the relation Δx>vFΔt[21]. Furthermore, the space increment Δxmust far smaller than the wavelength of electronsΔx<λe/10, and the time increment must be far smaller than the period of the electromagnetic fieldTl.

At the boundary, one-dimensional Mur absorbing boundary conditions are used [21-22]. To compare our results with [25] which use linear polarization for the vector potentialAwhere the chiral parameter T is zero, we considerer at the input boundary, a normalized Gaussian electronic wave packet, wheretgand τgdenote the peak position and the pulse width, respectively. Thus, by solving Eq. (57) directly in the time domain we can demonstrate the propagation of a wave packet through a barrier in real time.

Numerical simulations are shown in Fig. 7. The following parameters are used in our calculation: the peak position t0= 1.5 ps, the pulse width τg= 1.0 ps, the space increment Δx= 0.1 nm, the time incrementsΔt=5×105ps, and the height of the potential barrier H0= 500 meV.

Figure 7.

a) Numerical simulations of a wave packet given by N = | ϕ A | 2 + | ϕ B | 2 , tunneling through a barrier without pump beams. Figures (b)-(d) show the time sequence of a wave packet tunneling through a barrier with pump intensity I ω = 3 M W / c m 2 , Δ 0 = 5meV, and D = 350 nm. The light grey shows the barrier area.

When there is no pump beams, a perfect chiral tunneling can be found [see Fig. 7 (a)]. This result is consistent with that of Ling et al. [25]. But when the sample is irradiated by an intense non resonant laser beam, a reflected wave packet appears [see Fig. 7 (d)]. The perfect transmission is suppressed. By analyzing the transmitted wave packet and the reflected wave packet, we can obtain the tunneling rate.

To explain the suppression of chiral tunneling, we first investigate the chiral potential wave in the barrier within a rotating wave approximation [23, 24]. Figure 7 (a) shows the renormalized band as a function of momentum k with intensityIω=3MW/cm2.

Here, the important point is that we make±k0T1+α1.0073, where α=e2/(4πε0c)is the fine structure constant, so we verify that our chiral theory is correct. Without external beams, (fig. 7 (a)), we verify the known results on chiral tunneling [11]. Figures 7 (b)-(d) show the time sequence of a wave packet tunneling through a barrier with pump intensityIω=3MW/cm2, Δ0= 5meV, and D = 350 nm.

Figure 8.

Left: The reflectance R(circles) and the transmittance T (squares) of the barrier as a function of the detuning for I ω = 3 M W / c m 2 and D = 350 nm. Right: The transmittance versus I ω ( M W / c m 2 ) , having D as parameter. Here, k 0 2 T 2 ≈ 1 , that is the wave electric field is almost parallel to the wave magnetic field within the graphene device.

Under intense light beams, the dressed states are strongly mixed with valence states and conduction states. Therefore, the chiral symmetry of Dirac electrons in graphene can be broken and perfect chiral tunneling is strongly suppressed. Numerical results are shown in Fig. 8 (left) with pump intensity Iω=3MW/cm2and D = 350 nm. From Fig. 8(left) we can find that the transmission is strongly suppressed, even with laser detuning (e.g., Δ0= 10 meV, the transmittance is about 0.03).

Figure 8 (left), show that the reflectance decreases, and the transmittance increases as Δ0increases. The strong laser field can enhance band mixing and reduce the transmittance. If D increases we can see that the wide barrier can prolong the interaction time between electrons and photons, reduce the tunneling rate, and lower the threshold of the pump laser power (Figure 8 right). Our results can be compared with [25] which use linear polarization ofHint, however we think that inside of a plate of graphene, the eigenvectors of A have to be characterized as chiral waves, because the electrons have a chiral nature in a graphene device. Is for this reason that our numerical results correspond to half times the results shown in [25].

5. Zitterbewegung of optical pulses near the Dirac point inside a negative-zero-positive index metamaterial

This optical analog of the Zitterbewegung effect is similar to two-dimensional (2DPC) photonic crystal.

The Zitterbewegung effect can appear in solids, because the electron at the touching point between two interacting bands in a solid also obeys the massless Dirac equation [26]. In recent years, a growing attention has been devoted to the simulations or demonstrations of the Zitterbewegung effects in controllable physical systems, such as, ultra cold atoms, superconductors, semiconductor nanostructures with spin-orbit coupling, a single trapped ion, and graphene.

For a homogenous medium, when a light field is polarized in the z direction, the Beltrami’s equation×E±(r,t)=k±E±can be transformed as a Helmholtz equation written as


with a wave number k(ω)This equation can be written in the form of the Dirac equation


In the case of grapheme near ofε/εp,μ/μp0, there are two bands touch each other forming a double-cone structure with a linear dispersion, which can also be called Dirac Point. Near the Dirac point the light transport obeys the massless Dirac equation if we makevF=c(1±k0T)=vD. Also naturally, it would be of great interest to find out what is the condition to have the Dirac dispersion for the light field in a homogenous medium, the wave vector of a medium, k(ω), can be expanded as k(ω)=k(ωD)+(ωωD)/ ifk(ωD)=0, then we have a linear dispersionk(ω)=(ωωD)/vD, and the Dirac equation is


whereEz1(x,y,ω), Ez2(x,y,ω)are two eigenfunctions of the electrical fields (polarized in the z-direction) corresponding to the same frequencyω, and vDis the group velocity at the Dirac point (ωD). For media satisfying Eq. (58) withωD>0, the wave number, k(ω), varies from negative to zero and then to positive with the frequency, and so does the refractive index of the medium with zero refractive index atωD.

Note that the positive and negative branches of the band structure coexist. Equation (58) is the massless Dirac equation of the light fields in homogenous materials, which is the same as that of electrons in graphene [27]. Therefore, for the NZPI metamaterial [satisfying Eq. (58)], we will have the Dirac point with a double-cone structure for the light field at frequencyωD. It is expected that the propagation of the light field is analogous to that of the electron in graphene. In fact, in the 2DPC [28], the effective refractive index varies from negative to positive near this Dirac point. AtωD, Eq. (45) takes the form of diffusion equation,


which are the same as that of the massless Dirac equation at zero energy. [29]. Therefore, it can be predicted that the behavior of light fields at ωDhas the diffusive properties inside the medium of Eq. (2), like electrons at Dirac point of graphene [27-29]. At or nearωD, because ofk2=kx2+ky20, thekycomponent becomes a pure imaginary number for any realkx, thus the fields along theydirection between a given interval L have the total energy transmittance


which tells us that the propagation of light field at (or near) ωDexhibits the L / 1 scaling, a main characteristic of the diffusion phenomenon. This effect can be obtained if ε/εp,μ/μp0with low damping rates.

Here, the Dirac point with the double-cone structure for the light field can also be realized in a homogenous negative-zero-positive index (NZPI) medium [26], in which the two-dimensional Helmholtz equation could be written as the two-dimensional massless Dirac equation. The condition for the realization of the Dirac point in the homogenous optical medium is the index varying from negative to zero and then to positive with frequency, and the light field also obeys eq. (45) near the Dirac point.. Thus, in this work, we show that the Zitterbewegung effect with optical pulses appears near the Dirac point in NZPIM slabs. With the realization of the NZPIMs in experiments, we believe that it already has a great possibility to observe Zitterbewegung effect with optical pulses in the GHz region.

6. Theoretical description of minimal conductivity in graphene under chiral approach

The first results on the theoretical combination of graphene electron transport properties and a characteristic property of Dirac chiral fermions were obtained by Katsnelson et al, [2]. The intrinsic nature of Dirac fermions gives rise to minimal conductivity even for an ideal crystal, that is, without any scattering processes. The simplest way for the theoretical consideration is the Landauer approach [29]. Assuming that the sample is a ring of length Lyin y direction; author used the Landauer formula to calculate the conductance in the x direction (Fig. 9).

Figure 9.

Geometry of the sample. The direction of the current is parallel to the x axis [27].

The convenient boundary conditions are not physical, but to get finite transparency one should chooseLxLy.

In the coordinate representation the Dirac equation (59) at zero energy takes the form (Ψi=vDEzi)


The solutions of these equations areΨ1(xiy),Ψ2(x+iy), which are complex conjugated functions which are general solutions.

Due to periodicity in y direction both wave functions should be proportional toexp(ikyy)whereky=2πn/Ly,, n = 012,... This means that the dependence on the x is also fixed:

the wave functions are proportional toexp(2πnx/Ly). The introduced boundary conditions at the sample edges are x=0x = 0 andx=Lx. The assumption is that the sample is doped graphene with the Fermi energyEF=vkF=V0. The wave functions in the sample are supposed to have the same y-dependence, that is,Ψ1,2(x,y)=Ψ1,2(x)exp(ikyy).

Requiring continuity of the each wave function at the edges of sample, one can find the transmission coefficient:




Further, one should assume thatkFLx1and put ϕ=0in equation. Thus, the trace of the transparency which is just the conductance per valley per spin is


The conductance then equalsσLy/Lxand the conductivity ise2/πh. Experimentally, it is close to e2/h [31], that is in πtimes larger than present estimation.

Zitterbewegung − circular motion of elementary particles caused by an interference between positive and negative energy states − leads to the fluctuation of the position of an electron. This relativistic “jittering” of an electron in graphene could be interpreted in terms of classical physics as an interaction of electron with some potential caused by the presence of disorder in crystal. Therefore, the Zitterbewegung plays a role of “intrinsic” disorder in the system which appears in the presence of minimal conductivity of the ideal crystal (without scattering) even at zero temperatures.

Under this chiral approach, we show that the optical transparency of suspended graphene can be defined by the fine structure constant, α=e2/c, the parameter that describes coupling between light and relativistic electrons and is traditionally associated with quantum electrodynamics rather than condensed matter physics. Despite being only one atom thick, graphene is found to absorb a significant (α=0.78%toπα=2.3%) fraction of incident white light, which is a consequence of graphene’s unique electronic structure. This value translates directly into universal dynamic conductivity G=e2/4πwithin a few % accuracy.

There is a small group of phenomena in condensed matter physics, which are defined only by the fundamental constants and do not depend on material parameters. Examples are the resistivity quantum h/e2that appears in a variety of transport experiments, including the quantum Hall effect and universal conductance fluctuations, and the magnetic flux quantum h/eplaying an important role in the physics of superconductivity (here h is the Planck constant and e the elementary charge). By and large, it requires sophisticated facilities and special measurement conditions to observe any of these phenomena. Here, we show that such a simple observable as the visible transparency of graphene [1] is defined by the fine structure constant,α. Our results are in agreement with the theory of ideal two-dimensional (2D) Dirac fermions [2,3] and its recent extension into visible optics [4], which takes into account the triangular warping and nonlinearity of graphene’s electronic spectrum, but our result is π1times the result of [2,3].

Optical properties of thin films are commonly described in terms of dynamic or optical conductivity G. For a 2D Dirac spectrum with a conical dispersion relation ε=vF|k|(kthe wavevector). G was theoretically predicted [32, 33] to exhibit a universal valueG=e2/4, our result is G0=e2/4πif the photon energy E is much larger than both temperature and Fermi energyεF. Both conditions are stringently satisfied in our visible-optics experiments. The universal value of G also implies that all optical properties of graphene (its transmittance T, absorption Aabsand reflection R) can be expressed through fundamental constants only (T, Aabsand R are unequivocally related to G in the 2D case). In particular, it was noted by Kuzmenko et al [34] that Ttrasp=(1+2πG0/c)2(1+0.5πα)1παfor the normal light incidence. Our result is given byTtrasp1α. We emphasize that – unlike G – both T and R are observable quantities that can be measured directly by using graphene membranes.

7. Absorption of light by quasi 2D Dirac fermions

Here, following [35], we show how the universal value of graphene’s opacity can be understood qualitatively, without calculating its dynamic conductivity. Let a light wave with electric field Eand frequency ω fall perpendicular to a graphene sheet of a unit area. The incident energy flux is given by Wi=c/4π(|E|2)

Taking into account the momentum conservationkfor the initial |iand final |fstates, only the excitation processes pictured in Fig. 10 contribute to the light absorption. The absorbed energy Wa=nωis given by the number n of such absorption events per unit time and can be calculated by using Fermi’s golden rule as n=(2π/ƛ)(|L|2)Dwhere L is the matrix element for the interaction between light and Dirac fermions, and D is the density of states at ε=E/2=ω/2(see Fig. A1). For quasi 2D Dirac fermions, D(ω/2)=ω/π2vF2and is a linear function of ε.

The interaction between light and Dirac fermions is generally described by the Hamiltonian


where the first term is the standard Hamiltonian for 2D Dirac quasiparticles in graphene [1] and

HA=vFσ^(e/c)A=vFσ^(e/iω)Edescribes their interaction with electromagnetic field.

Here A=(ic/ω)Eis the vector potential and σthe standard Pauli matrices. Averaging over all initial and final states and taking into account the valley degeneracy, our calculations yield


This results in Wa=(e2/4)|E|2and, consequently, absorption Aabs=Wa/Wi=πe2/c)=παboth of which are independent of the material parameter vFthat cancels out in the calculations ofWa. Also note that the dynamic conductivity G=Wa/|E|2is equal toe2/4. Our result isAabsWa/Wi=πe2/πc)=αin correspondence with vF=c(1k0T)Because graphene practically does not reflect light (R<<1 as discussed above), its opacity (1 – T) is dominated by the derived expression forAabs.

In the case of a zero-gap semiconductor with a parabolic spectrum (e.g., bilayer graphene at low ε), the same analysis based on Fermi’s golden rule yieldsAabs=2πα. This shows that the fact that the optical properties of graphene are defined by the fundamental constants is related to its 2D nature and zero energy gap and does not directly involve the chiral properties of Dirac fermions.

On a more general note, graphene’s Hamiltonian H has the same structure as for relativistic electrons (except for coefficient vF instead of the speed of light c). The interaction of light with relativistic particles is described by a coupling constant, a.k.a. the fine structure constant. The Fermi velocity is only a prefactor for both HamiltoniansH0andHAand, accordingly, one can expect that the coefficient may not change the strength of the interaction, as indeed our calculations show.

Figure 10.

Excitation processes responsible for absorption of light in graphene.

Thus, we have found that the visible opacity of suspended graphene is given by π α within a few percent accuracy and increases proportionally to the number of layers N for few-layer graphene. Its dynamic conductivity at visible frequencies is remarkably close to the universal value of(e2/4)N. The agreement between the experiment and theory is particularly striking because it was believed the universality could hold only for low energies (< 1eV) beyond which the electronic spectrum of graphene becomes strongly warped and nonlinear and the approximation of Dirac fermions breaks down.

Electrons from the valance band (bottom) are excited into empty states in the conduction band with conserving their momentum and gaining the energy E= h ω.

The approximation of 2D Dirac fermions is valid for graphene only close to the Dirac point and, for higher energies ε, one has to take into account such effects in graphene’s band structure as triangular warping and nonlinearity.

Most theories suggest4e2/πh, which is about π times smaller than the typical values observed experimentally. The reason of disagreement is still an open question − “mystery of a missing pie”. A value of conductivity e2/hat the Dirac point can emerge in case of preservation of chiral symmetry by disorder or when the dominant disorder does not scatter electrons between two valleys. Furthermore, macroscopic inhomogeneity (on the scale larger than the mean free path) is also important for measurements ofσmin.

8. Conclusions

In this chapter we presented a short review on the chiral propagation of electron waves in monolayer graphene and optical simulation with optical field in the negative-zero-positive index metamaterial NZPIM and its close connection. Section II presented an enhanced vector diagram of Maxwell’s equations for chiral media with quasi parallel electromagnetic fields,EH. Chiral waves in graphene acting as metamaterial media were discussed in section III. In section IV, two equations and tunneling rate of Dirac electron in graphene were derived. In Section V, Zitterbewegung of optical pulses near the Dirac point inside a negative-zero-positive index metamaterial, showed that the chiral field near the Dirac point becomes a diffusive wave. The last sections described the theoretical description of minimal conductivity in graphene under chiral approach and absorption of light by quasi 2D Dirac fermions.

We reviewed the formulation of graphene’s massless Dirac Hamiltonian, under the chiral electromagnetism approach, like a metamaterial medium, hopefully demystifying the material’s unusual chiral, relativistic, effective theory. The novel result here was that in our theory we did not make cvF[1, 2, 27-35], but we obtained vFas vF=c(1k0T)ifk0T>0orvF=c(1+k0T)ifk0T<0. These results were derived of the Chiral Electrodynamics with T as the chiral parameter andk0=ω/c[5, 6, 7]. This situation corresponds to chiral electrons where the Efield is almost parallel to theHfield, and the movement of electron is helical in the valence band of the membrane graphene. Semi classically, for an electron in hydrogen atom, this condition was shown by Huang (2006) for Ewave perpendicular to Hwave with radiation (E=H) and by Torres-Silva for Eparallel toH, without radiation (E=iH), see refs ([5, 36, 37]). Here, we have extended this result to the electron bound to a carbon atom.

There is a small group of phenomena in condensed matter physics, which are defined only by the fundamental constants and do not depend on material parameters. Examples are the resistivity quantum h/e2that appears in a variety of transport experiments, including the quantum Hall effect and universal conductance fluctuations, and the magnetic flux quantum h/eplaying an important role in the physics of superconductivity. In this chapter we founded that under this chiral approach, the optical transparency of suspended graphene can be defined by the fine structure constant,α=e2/c, that isTtrasp1α, andvF/c=α=(1k0T).

Also we give a support to the similitude of the band structure of a macroscopic photonic crystal with the electronic band structure of graphene, which is experimentally much more difficult to access, allows the experimental study of various relativistic phenomena. With our analytical and numerical analysis we hope to verify that, similar behaviors exist to electrons in graphene treated as mass-less particle, ie, electron wave discovered only 8 years ago graphene is already one of the most studied carbon allotropes. But this material still poses a lot of theoretical and experimental questions, which have to be answered.


I wish to thank to colleague Jorge Benavides Silva of the EIEE for many useful discussions on particles and photons

How to cite and reference

Link to this chapter Copy to clipboard

Cite this chapter Copy to clipboard

H. Torres-Silva (November 14th 2012). Chiral Waves in Graphene Medium and Optical Simulation with Metamaterial, Solutions and Applications of Scattering, Propagation, Radiation and Emission of Electromagnetic Waves, Ahmed Kishk, IntechOpen, DOI: 10.5772/51328. Available from:

chapter statistics

1252total chapter downloads

2Crossref citations

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

Solutions and Applications of Scattering, Propagation, Radiation and Emission of Electromagnetic Waves

Edited by Ahmed Kishk

Next chapter

Numerical Modeling of Electromagnetic Wave Propagation Through Bi-Isotropic Materials

By I. Barba, A. Grande, A.C.L. Cabeceira, A. Gómez, J.A. Pereda and J. Represa

Related Book

Advancement in Microstrip Antennas with Recent Applications

Edited by Ahmed Kishk

First chapter

Design Techniques for Conformal Microstrip Antennas and Their Arrays

By Daniel B. Ferreira, Cristiano B. de Paula and Daniel C. Nascimento

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