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 is the Fermi velocity and 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 relation, 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 velocity. 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 relationis that of relativistic massless fermions with particle-hole symmetry. In graphene these massless fermions propagate with a velocity. 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 [1, 2, 27-35], but we obtain as iforif. These results are derived of the Chiral Electrodynamics with T as the chiral parameter andwhen the electric wave is quasi parallel to the magnetic wave.
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,. 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
The idea of representing Maxwell’s time-harmonic equations in homogeneous isotropic media by vector diagram as put forward by Wilton  and by S. Uckun  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  for achiral media has vectors only in one plane with H vector normal to it.
2.1. Vector diagram construction
Chirality is introduced into the theory by defining the following constitutive relations to describe the isotropic chiral medium 
Where the chirality admittance indicates the degree of chirality of the medium, and the andare permittivity and permeability of the chiral medium, respectively. Since and are polar vectors and and are axial vectors, it follows that and are true scalars and is 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, changes in sign whereas and remain unchanged.
For a graphical representation of the above relationships, following Wilton’s procedure , 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 and the vector must be perpendicular to both and.
Since always, this conditions will hold identically if is expressed as the curl of a vector potential since the divergence of the curl of a vector is identically zero, thus
That is shown in Figure 1.
Here is arbitrary, so in order to specify, for unique, we may choose
Divide both sides of Equation (12) by and reorganize it to get
( Figure 1, shows this vectorial equation).
The difference between our approach and the Uckun’s procedure , is that we take the chiral media characterized by and. In this form we can obtain the condition of quasi spatial parallel condition between and where the main equation is like a Beltrami equation which is important for the numerical simulation of graphene systems.
by rearranging this equation we obtain
will be obtained as shown in Figure 1. In this figure we put.
will be derived. To find the projection of onto, from Equation(15) , take the gradient of both sides and divide by scalar value to normalize to a unit vector. So
Similarly, getting gradient of both sides of Equation(11), using the vector identify for and normalizing by we have:
So it will be obtained as parallel component of A to coordinate.
Now let us examine derivation of some relations from the diagram. For example, it is seen that the component of and parallel to must be equal and opposite. By taking the divergence of Equation (14) and using can be shown 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.. Using Equation (4) convert this immediately into continuity equation as, expected. Again, as seen in Figure 1, and are equal and opposite vectors. From Equation(8), taking curl of both side and using the vector identity will 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
and value of from Equation (1) will give
By using the same coordinates axes, and 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 either parallel or perpendicular to any of the other vectors and simplify its relationship to those vectors. In Figure 1a, if the chirality factor, goes to zero, point, and approach point, , and respectively, in which case the diagram will be the same as in Reference  for linear, homogeneous, isotropic achiral medium. If then E is quasi parallel to B, and parallel toso all vectors remain in an only plane. (see Figure 2).
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 . In terms of chiral magnetic potentialwith, andwe have the wave equation
The dispersion relation of the transversal wave is .
So we have with if we consider a bi dimensional graphene system.
The novel result here is that in our chiral theory we do not make but we obtain as iforif. These results are derived of the Chiral Electrodynamics with T as the chiral parameter and[5, 6, 7]. This situation corresponds to a model of chiral electron where the field is almost parallel to thefield, 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  for wave perpendicular to wave with radiation and by Torres-Silva for parallel to, 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 . 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 . 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 .
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 . 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 . 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 field, and the magnetization depends not only on the magnetic field, and we may have, for example, constitutive relations given by the Born-Federov formalism .
The pseudoscalar represents the chirality of the material and it has length units. In the limit, 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
where, and, is 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 electricand magnetic field and and the wave vector will 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, and, , and will form a left-handed triplet of vectors. This is just the peculiar case for left handed materials where the effective permittivity and the effective permeability 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 of.
To advance in our propose we considerer other more popular representation to describe a chiral medium,  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 and, it is the Tellegen medium; if and, 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. or, until we see the fact that artificial Veselago’s medium  was successfully realized in certain frequency bands . For the case, , we have to ask the following question: can strong chiral medium exist?
In Ref. , 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.
Mathematically, the amount of energy density propagated is proportional to the magnitude of the Poynting vector, where. With the condition, so 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 putand solving the Maxwell equations with the Born Fedorov relations we findso, 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 vanishes. This can be achieved if is zero all the time everywhere in the region of space under consideration, i.e., (see Fig. 4). Examination of Eq. (39) shows thatif. In this last case, a particular solution of Eq. (31,3 2) is when, where we have the condition, and, so we find the Beltrami force free equation and the vector Poynting vanishes [15-16].
In terms of Eqs. (41-42), iffor 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, , and, we easily show thatand 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 negativeand 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 , in which and while, the parameters in DBF representation become,and, 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,. 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. That is to say, strong chirality does not necessarily lead to strong spatial dispersion.
Based on Eqs. (41)-(43), we have computed T and or versus, as shown in Figs. 2 and 3. When is very close to, the value of T is quite large, indicating a strong spatial dispersion. Hence the singular point is the very point of traditional limitation. However, withcontinuously increasing, the spatial dispersion strength falls down very quickly. Therefore, if is not around, e.g. or, 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.If or, we have that the Fermi velocity is given by, so it appears a connection between the chiral electrodynamics and the fine structure constant.
From Figure 4. we see that near of the 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 makewith
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 configuration. Dirac's four-component equation for the relativistic electron is [5-7]. (see our chapter “Chiral Transverse Electromagnetic standing waves with in 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 is deduced from the chiral electrodynamics with, 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 making, no making as other authors do it. In Figure 5 we have.
The Hamiltonian commutes with the momentum vector. In order to resolve this degeneracy we seek a dynamical variable which commutes with both H and. Such a variable is, whereis the matrix Pauli. The eigenfunctions of the commuting variables p andare simultaneous:
Thus for a simultaneous eigenstate of and, the value of 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 are the four components of the state function and four numbers to be determined. Then E can have either of the two values.
We now demand that be also an eigenfunction of belonging to one of the eigenvalues, say, where.
The eigenvalue equation is
Since can be given either of the two values and, the two values, we have found for given p four linearly independent plane wave solutions. It is easily verified that they are mutually orthogonal.
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 graphene.
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 height, 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. We choose to ensure that there is no inter-band absorption inside the barrier. Meanwhile, 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) with, so where t means transposed wave function.
where is the wave function, is the unperturbed Dirac Hamiltonian obtained from the chiral electrodynamics , are the Pauli matrices, is the Fermi velocity, is the height of the potential barrier, I is the unit matrix, and is the interaction Hamiltonian with a chiral electromagnetic potential.
where e is the electron charge and the chiral potential vector is with and 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 and 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 . 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 . For the nodal variables we can apply the usual Lee discretization method in 2-D systems with(See Figure 6).
For any function of space and time, 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 . Furthermore, the space increment must far smaller than the wavelength of electrons, and the time increment must be far smaller than the period of the electromagnetic field.
At the boundary, one-dimensional Mur absorbing boundary conditions are used [21-22]. To compare our results with  which use linear polarization for the vector potentialwhere the chiral parameter T is zero, we considerer at the input boundary, a normalized Gaussian electronic wave packet, whereand denote 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 = 1.5 ps, the pulse width = 1.0 ps, the space increment = 0.1 nm, the time increments, and the height of the potential barrier = 500 meV.
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. . 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 intensity.
Here, the important point is that we make, where 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 . Figures 7 (b)-(d) show the time sequence of a wave packet tunneling through a barrier with pump intensity, = 5meV, and D = 350 nm.
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 and D = 350 nm. From Fig. 8(left) we can find that the transmission is strongly suppressed, even with laser detuning (e.g., = 10 meV, the transmittance is about 0.03).
Figure 8 (left), show that the reflectance decreases, and the transmittance increases as increases. 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  which use linear polarization of, 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 .
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 . 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 equationcan be transformed as a Helmholtz equation written as
with a wave number This equation can be written in the form of the Dirac equation
In the case of grapheme near of, 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 make. 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, , can be expanded as so if, then we have a linear dispersion, and the Dirac equation is
where, are two eigenfunctions of the electrical fields (polarized in the z-direction) corresponding to the same frequency, and is the group velocity at the Dirac point (). For media satisfying Eq. (58) with, the wave number, , 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.
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 . 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. It is expected that the propagation of the light field is analogous to that of the electron in graphene. In fact, in the 2DPC , the effective refractive index varies from negative to positive near this Dirac point. At, Eq. (45) takes the form of diffusion equation,
which are the same as that of the massless Dirac equation at zero energy. . Therefore, it can be predicted that the behavior of light fields at has the diffusive properties inside the medium of Eq. (2), like electrons at Dirac point of graphene [27-29]. At or near, because of, thecomponent becomes a pure imaginary number for any real, thus the fields along thedirection between a given interval L have the total energy transmittance
which tells us that the propagation of light field at (or near) exhibits the L / 1 scaling, a main characteristic of the diffusion phenomenon. This effect can be obtained if with 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 , 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, . 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 . Assuming that the sample is a ring of length in y direction; author used the Landauer formula to calculate the conductance in the x direction (Fig. 9).
The convenient boundary conditions are not physical, but to get finite transparency one should choose.
In the coordinate representation the Dirac equation (59) at zero energy takes the form ()
The solutions of these equations are,, which are complex conjugated functions which are general solutions.
Due to periodicity in y direction both wave functions should be proportional towhere,, n = 0,±1,±2,... This means that the dependence on the x is also fixed:
the wave functions are proportional to. The introduced boundary conditions at the sample edges are x = 0 and. The assumption is that the sample is doped graphene with the Fermi energy. The wave functions in the sample are supposed to have the same y-dependence, that is,.
Requiring continuity of the each wave function at the edges of sample, one can find the transmission coefficient:
Further, one should assume thatand put in equation. Thus, the trace of the transparency which is just the conductance per valley per spin is
The conductance then equalsand the conductivity is. Experimentally, it is close to /h , 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, , 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 (to) fraction of incident white light, which is a consequence of graphene’s unique electronic structure. This value translates directly into universal dynamic conductivity 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 that appears in a variety of transport experiments, including the quantum Hall effect and universal conductance fluctuations, and the magnetic flux quantum playing 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  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 , which takes into account the triangular warping and nonlinearity of graphene’s electronic spectrum, but our result is times 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 (the wavevector). G was theoretically predicted [32, 33] to exhibit a universal value, our result is if the photon energy E is much larger than both temperature and Fermi energy. 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 and reflection R) can be expressed through fundamental constants only (T, and R are unequivocally related to G in the 2D case). In particular, it was noted by Kuzmenko et al  that for the normal light incidence. Our result is given by. 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 , 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 and frequency ω fall perpendicular to a graphene sheet of a unit area. The incident energy flux is given by
Taking into account the momentum conservationfor the initial and final states, only the excitation processes pictured in Fig. 10 contribute to the light absorption. The absorbed energy is given by the number n of such absorption events per unit time and can be calculated by using Fermi’s golden rule as where L is the matrix element for the interaction between light and Dirac fermions, and D is the density of states at (see Fig. A1). For quasi 2D Dirac fermions, and 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  anddescribes their interaction with electromagnetic field.
Here is 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 and, consequently, absorption both of which are independent of the material parameter that cancels out in the calculations of. Also note that the dynamic conductivity is equal to. Our result isin correspondence with Because graphene practically does not reflect light (R<<1 as discussed above), its opacity (1 – T) is dominated by the derived expression for.
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 yields. 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 Hamiltoniansandand, accordingly, one can expect that the coefficient may not change the strength of the interaction, as indeed our calculations show.
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. 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 suggest, 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 at 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.
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,. 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 [1, 2, 27-35], but we obtained as iforif. These results were derived of the Chiral Electrodynamics with T as the chiral parameter and[5, 6, 7]. This situation corresponds to chiral electrons where the field is almost parallel to thefield, 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 wave perpendicular to wave with radiation () and by Torres-Silva for parallel to, without radiation (), 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 that appears in a variety of transport experiments, including the quantum Hall effect and universal conductance fluctuations, and the magnetic flux quantum playing 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,, that is, and.
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