**Electronic and Vibrational Properties of Adsorbed and Embedded Graphene and Bigraphene with Defects**

Alexander Feher, Eugen Syrkin, Sergey Feodosyev, Igor Gospodarev, Elena Manzhelii, Alexander Kotlar and Kirill Kravchenko

Additional information is available at the end of the chapter

http://dx.doi.org/10.5772/50562

## **1. Introduction**

34

[30] Schwinger, J., Tsai, W. Y. and Erber, T. (1976). Classical and quantum theory of synergic

[31] Schwinger, J. (1970).Particles, Sources and Fields Vol. I (Addison-Wesley, Reading,

[32] Schwinger, J. (1988). Particles, sources and fields, (Addison-Wesley Publishing

[33] Sokolov, A. A.; Loskutov, Yu. M. and Ternov, I.M. (1962). Quantum mechanics,

[34] Tamm, I. E. and Frank, I. M. (1937). The coherent radiation of a fast electron in a

[35] Yildiz, A. (1977). Meissner effect in gauge fields *Phys. Rev. D* 16, No. 12 , 3450.

ˇ *Annals of Physics* (NY) 96, 303.

synchrotron-Cerenkov radiation,

Mass.).

Company, Inc.).

134 New Progress on Graphene Research

(Moscow). (in Russian).

medium, *Dokl. Akad. Nauk SSSR* 14, 109.

The perennial interest in studying the physical properties of nanofilms has increased sub‐ stantially over the last few years due to the development of nanotechnologies and the syn‐ thesis of new compounds – especially those based on carbon, which are extremely interesting for both fundamental research and potential applications.

An important feature of carbon nanofilms (including those with defects) is a close relation between the electronic and phonon properties, which is exhibited, for example, in the gra‐ phene-based systems with superconducting properties [1,2].

It is well known that graphene monolayers cannot exist as planar objects in the free state, because in flat 2D-crystals the mean-square amplitudes of the atoms in the direction normal to the layer plane diverge even at *T* =0 (see, e.g., [3]). So we can study and practically apply only such graphene, which is deposited on a certain substrate providing the stability of the plane carbon nanofilms (see, e.g., [4-6]). Only small flakes can be detached from the sub‐ strate and these flakes immediately acquire a corrugated shape [7]. When studying the elec‐ tronic properties of graphene a dielectric substrate is often used. The presence of the substrate greatly increases the occurrence of various defects in graphene and carbon nano‐ films. Our investigations make it possible to predict the general properties of phonon and electron spectra for graphene and bigraphene containing different defects.

This chapter consists of three sections: first section is devoted to the calculation of local dis‐ crete levels in the electron spectra of graphene with different defects. In the second section

© 2013 Feher et al.; licensee InTech. This is an open access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/3.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited. © 2013 Feher et al.; licensee InTech. This is a paper distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/3.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

we describe the electronic properties of bilayered graphene and, finally, the third section deals with the influence of defects on electron spectra of bigraphene.

the symmetric and antisymmetric combinations of wave functions in the two sublattices in‐ tersect at the edge of the Brillouin zone, which leads to a cone-shaped energy spectrum near the *K* and *K*´ points of the first Brillouin zone. The electrons obey the linear dispersion law

Electronic and Vibrational Properties of Adsorbed and Embedded Graphene and Bigraphene with Defects

The electronic spectrum of graphene can be described by a strong coupling approximation, and it is sufficient to consider the interaction between nearest neighbors only (see, e.g.,

where *i* and *j* are the labels of the nodes of the two-dimensional lattice, *ε<sup>i</sup>* is the energy of

Curve 1 in Figure 1 shows the density of electronic states of graphene as calculated using the method of Jacobi matrices [16, 17]. In a perfect graphene the local Green's function

lattices. Peculiarity of the density of states at *ε* =*ε*(*K*) (the value *ε*(*K*) corresponds to the Fer‐ mi energy *εF* in graphene) determines the behavior of the real part of the Green's function near *εF* . For a wide class of perturbations caused by defects we can find, using the Lifshitz equation [21], quasilocalized states in the interval −*ε*(*M* ), *ε*(*M* ) (in this model *ε*(*M* )= *J* ). This equation, which determines the energy of these states, can be written as (see, e.g.,

^ )−<sup>1</sup> <sup>|</sup>*i* coincides with the total Green's function

^ )−<sup>1</sup> <sup>|</sup>*i* because of the physical equivalence of the atoms of both sub‐

Re*G*(*ε*)=*S*(*ε*, *Λik* ) (2)

Im*G*(*ε* + *iγ*, *i*) of impurity atoms are calculated in

^

( *Λik* are matrix


http://dx.doi.org/10.5772/50562

137


is the so-called overlap integral.

where the *S*(*ε*, *Λik* ) function is determined by the perturbation operator *Λ*

*<sup>π</sup>* lim *γ*↓0

and with the overlap integral *Ji*<sup>0</sup> =(1 + *η*)*J* , the function *S*(*ε*, *ε*˜, *η*) has the form

[6]. For an isolated substitutional impurity with the energy *ε*<sup>0</sup> =*ε*˜ of the impurity node *i* =0

elements of this operator on defined basis).

The local spectral densities *ρ*(*ε*, *<sup>i</sup>*)<sup>≡</sup> <sup>1</sup>

(in ordinary metals and semiconductors the dispersion law is parabolic).

*H* ^ <sup>=</sup>∑ *i εi*

[5,6,18-20]). The corresponding Hamiltonian is

electron at node *i*, and *Jij*

^ −*H*

1 *<sup>N</sup>* ∑ *i*=1 *N i* |(*εI* ^ −*H*

*G*(*ε*, *i*)=*i* |(*εI*

*G*(*ε*)= lim *N* →*∞*

[3,17])

#### **2. Impurity levels in the electron spectra of graphene**

The exceptionality of graphene is manifested in the phonon and electron properties. Gra‐ phene is a semimetal whose valence and conduction bands touch at the points K and K' of the Brillouin zone [8,9]. In the pure graphene unique electronic properties are manifested by the charge carriers behaving as massless relativistic particles - the dependence of energy on the momentum is linear rather than - as in ordinary solids - quadratic. Thus, the lower-di‐ mensionality affects the formation of phonon localized states [10] and also the formation of localized states in the electronic spectrum. Absent gap between the valence and conduction bands is a consequence of the symmetry between two equivalent sublattices in graphene [11]. Presence of impurities lowers the symmetry. The influence of vacancies placed into one of the graphene sublattices was investigated in [12], where it was shown that the equiva‐ lence of the sublattices is broken.

In this section we describe the characteristics of localized and local states present in gra‐ phene due to the impurities of nitrogen and boron, respectively. The presence of the sub‐ strate greatly increases the possibility to introduce various defects into graphene. For example, in the graphene deposited on silicon, vacancies can occur [13, 14], whereas in graphite (a set of weakly interacting graphene monolayers) vacancies heal and form a stack‐ ing fault with local fivefold symmetry axis [15]. Impurity atoms embedded in graphene may lead to the appearance of impurity states outside the band of quasi-continuous spectrum. At low impurity concentrations (when impurity is considered as an isolated defect) these states appear in the form of local discrete levels (LDL).

Although such levels in various quasiparticle spectra have been known and studied over 60 years, an adequate description is still absent, even in the harmonic approximation for suffi‐ ciently realistic models of the crystal lattice. The dependence of the appearance conditions and characteristics of LDL on the parameters of a perfect lattice and defect was identified only in the most general terms. However, LDL may be used as an important source of infor‐ mation about the defect structure and force interactions in real crystals. To extract such in‐ formation it is useful to have analytical expressions that relate main characteristics of LDL to the parameters of both the defect and the host lattice.

Here we present the results of our calculations and analyses of the characteristics of the elec‐ tronic local discrete levels for substitutional impurities in graphene, especially for a boron substitutional impurity, using an analytical approximation based on the Jacobi matrices method [16, 17].

The fact that the charge carriers in graphene are formally described by the Dirac equation and not by the Schrödinger equation is due to the symmetry of the crystal lattice of gra‐ phene, which consists of two equivalent carbon sublattices. Electronic subbands formed by the symmetric and antisymmetric combinations of wave functions in the two sublattices in‐ tersect at the edge of the Brillouin zone, which leads to a cone-shaped energy spectrum near the *K* and *K*´ points of the first Brillouin zone. The electrons obey the linear dispersion law (in ordinary metals and semiconductors the dispersion law is parabolic).

we describe the electronic properties of bilayered graphene and, finally, the third section

The exceptionality of graphene is manifested in the phonon and electron properties. Gra‐ phene is a semimetal whose valence and conduction bands touch at the points K and K' of the Brillouin zone [8,9]. In the pure graphene unique electronic properties are manifested by the charge carriers behaving as massless relativistic particles - the dependence of energy on the momentum is linear rather than - as in ordinary solids - quadratic. Thus, the lower-di‐ mensionality affects the formation of phonon localized states [10] and also the formation of localized states in the electronic spectrum. Absent gap between the valence and conduction bands is a consequence of the symmetry between two equivalent sublattices in graphene [11]. Presence of impurities lowers the symmetry. The influence of vacancies placed into one of the graphene sublattices was investigated in [12], where it was shown that the equiva‐

In this section we describe the characteristics of localized and local states present in gra‐ phene due to the impurities of nitrogen and boron, respectively. The presence of the sub‐ strate greatly increases the possibility to introduce various defects into graphene. For example, in the graphene deposited on silicon, vacancies can occur [13, 14], whereas in graphite (a set of weakly interacting graphene monolayers) vacancies heal and form a stack‐ ing fault with local fivefold symmetry axis [15]. Impurity atoms embedded in graphene may lead to the appearance of impurity states outside the band of quasi-continuous spectrum. At low impurity concentrations (when impurity is considered as an isolated defect) these states

Although such levels in various quasiparticle spectra have been known and studied over 60 years, an adequate description is still absent, even in the harmonic approximation for suffi‐ ciently realistic models of the crystal lattice. The dependence of the appearance conditions and characteristics of LDL on the parameters of a perfect lattice and defect was identified only in the most general terms. However, LDL may be used as an important source of infor‐ mation about the defect structure and force interactions in real crystals. To extract such in‐ formation it is useful to have analytical expressions that relate main characteristics of LDL to

Here we present the results of our calculations and analyses of the characteristics of the elec‐ tronic local discrete levels for substitutional impurities in graphene, especially for a boron substitutional impurity, using an analytical approximation based on the Jacobi matrices

The fact that the charge carriers in graphene are formally described by the Dirac equation and not by the Schrödinger equation is due to the symmetry of the crystal lattice of gra‐ phene, which consists of two equivalent carbon sublattices. Electronic subbands formed by

deals with the influence of defects on electron spectra of bigraphene.

**2. Impurity levels in the electron spectra of graphene**

lence of the sublattices is broken.

136 New Progress on Graphene Research

appear in the form of local discrete levels (LDL).

the parameters of both the defect and the host lattice.

method [16, 17].

The electronic spectrum of graphene can be described by a strong coupling approximation, and it is sufficient to consider the interaction between nearest neighbors only (see, e.g., [5,6,18-20]). The corresponding Hamiltonian is

$$\stackrel{\frown}{H} = \sum\_{i} \varepsilon\_{i} \mid \text{ii} \parallel - \sum\_{i,j} J\_{ij} \mid i \not\mid \tag{1}$$

where *i* and *j* are the labels of the nodes of the two-dimensional lattice, *ε<sup>i</sup>* is the energy of electron at node *i*, and *Jij* is the so-called overlap integral.

Curve 1 in Figure 1 shows the density of electronic states of graphene as calculated using the method of Jacobi matrices [16, 17]. In a perfect graphene the local Green's function *G*(*ε*, *i*)=*i* |(*εI* ^ −*H* ^ )−<sup>1</sup> <sup>|</sup>*i* coincides with the total Green's function *G*(*ε*)= lim *N* →*∞* 1 *<sup>N</sup>* ∑ *i*=1 *N i* |(*εI* ^ −*H* ^ )−<sup>1</sup> <sup>|</sup>*i* because of the physical equivalence of the atoms of both sub‐ lattices. Peculiarity of the density of states at *ε* =*ε*(*K*) (the value *ε*(*K*) corresponds to the Fer‐ mi energy *εF* in graphene) determines the behavior of the real part of the Green's function near *εF* . For a wide class of perturbations caused by defects we can find, using the Lifshitz equation [21], quasilocalized states in the interval −*ε*(*M* ), *ε*(*M* ) (in this model *ε*(*M* )= *J* ). This equation, which determines the energy of these states, can be written as (see, e.g., [3,17])

$$\text{Re}G(\varepsilon) = \text{S}\begin{pmatrix} \varepsilon \ \Lambda \end{pmatrix} \tag{2}$$

where the *S*(*ε*, *Λik* ) function is determined by the perturbation operator *Λ* ^ ( *Λik* are matrix elements of this operator on defined basis).

The local spectral densities *ρ*(*ε*, *<sup>i</sup>*)<sup>≡</sup> <sup>1</sup> *<sup>π</sup>* lim *γ*↓0 Im*G*(*ε* + *iγ*, *i*) of impurity atoms are calculated in [6]. For an isolated substitutional impurity with the energy *ε*<sup>0</sup> =*ε*˜ of the impurity node *i* =0 and with the overlap integral *Ji*<sup>0</sup> =(1 + *η*)*J* , the function *S*(*ε*, *ε*˜, *η*) has the form

**Figure 1.** Electronic density of states of perfect graphene (curve 1) and the local density of states for an isolated boron substitutional impurity (curve 2).

$$S(\varepsilon, \ \tilde{\varepsilon}, \ \eta) = \frac{(1 + \eta)^2}{\tilde{\varepsilon} + \varepsilon \eta (2 + \eta)} \tag{3}$$

**Figure 2.** Graphical solution of equation (2) for boron substitutional impurity in graphene. Curve 1 is the electronic

Electronic and Vibrational Properties of Adsorbed and Embedded Graphene and Bigraphene with Defects

Such expressions were obtained in [23] for localized vibrations in the phonon spectrum of a three-dimensional crystal. Authors proposed an analytical approximation of the basic char‐ acteristics of local vibrations based on the rapid convergence of the real part of the Green's function outside the band of quasi-continuous spectrum using the method of Jacobi matri‐

Let us, briefly, to the extent necessary to understand the use of the classification of the eigen‐ functions of Hamiltonian (1), to present the basics of the method of Jacobi matrices. This method allows, without finding the dispersion laws, to calculate directly the local partial Green's functions of the system, corresponding to the perturbation of one or more atoms.

space of electronic excitations of atoms. Its dimension is *qN* , where *N* is the number of atoms in the system, and *q* is the dimension of the displacement of a single atom (*q = 1, 2, 3)*. Vectors of the space *H* are denoted by an arrow above the symbol, and "ordinary" *q*-dimen‐

If, using the generating vector (*p* is the number of excited atoms) and the

^ *nh* →

vectors of this sequence forms, in the *H* space, a cyclic subspace invariant to the operator *H*

This subspace contains, within itself, all the atomic displacements generated by the vector

<sup>0</sup> . The corresponding partial Green's function is determined as a matrix element

→

<sup>0</sup>}*n*=0 *<sup>∞</sup>* , then the linear envelope covering the

is the eigenvalue. Quantity *ρ*(*λ*)≡*π* <sup>−</sup><sup>1</sup>

<sup>0</sup>∈*H* , where *H* is the

^ .

Im*G*00(*λ*) is

 and 2′ are

http://dx.doi.org/10.5772/50562

139

density of states of ideal graphene, curve 2 is the corresponding real part of the Green's function. Curves 1′

the "approximations of two moments" of these functions. Curve 3 represents the function *S*(ε) .

This perturbation is described by the so-called generating vector *h*

0) , where *<sup>λ</sup>* <sup>≡</sup>*<sup>ω</sup>* <sup>2</sup>

ces [16,17].

*h* →

*G*00(*λ*)≡(*h*

→ <sup>0</sup>, *εI* ^ −*H* ^ <sup>−</sup><sup>1</sup> *h* →

sional vectors are in bold italics.

Hamiltonian (1), we construct the sequence {*H*

For a nitrogen impurity *ε*˜ −*ε*(*K*)≈ −0.525*J* and *η* ≈ −0.5 (according to [4]). As shown in [22], equation (2) has a solution for both interval −*ε*(*M* ), *ε*(*K*) and interval *ε*(*K*), *ε*(*M* ) .

The local density of states of the nitrogen substitutional impurity calculated in [11] has qua‐ si-local maxima in both intervals. For an boron substitutional impurity ( *η* ≈0.5 ) [6], quasilocalized states are absent in the −*ε*(*M* ), *ε*(*M* ) interval [12]. Figure 2 shows the graphical solution of the Lifshitz equation (2) for a given impurity atom. In this case the Lifshitz equa‐ tion has no solutions in interval −*ε*(*M* ), *ε*(*M* ) (corresponding dependences *S*(*ε*) are shown as curves 3 in Figure 2). The local Green's function of the boron impurity (curve 2 in Figure 1) has two poles outside the band of quasi-continuous spectrum, which are called local dis‐ crete levels and which are also solutions of equation (2). As is clearly seen in Figure 1 the area under the curve 2 is smaller (by the sum of the residues at these poles) than the area under the curve 1.

Local discrete levels can be an important source of information about defective structure and force interactions in real crystals. To extract this useful information we should have ana‐ lytical expressions that relate the main characteristics of LDL (primarily their energy) to the parameters of the defect and the host lattice.

Electronic and Vibrational Properties of Adsorbed and Embedded Graphene and Bigraphene with Defects http://dx.doi.org/10.5772/50562 139

**Figure 2.** Graphical solution of equation (2) for boron substitutional impurity in graphene. Curve 1 is the electronic density of states of ideal graphene, curve 2 is the corresponding real part of the Green's function. Curves 1′ and 2′ are the "approximations of two moments" of these functions. Curve 3 represents the function *S*(ε) .

Such expressions were obtained in [23] for localized vibrations in the phonon spectrum of a three-dimensional crystal. Authors proposed an analytical approximation of the basic char‐ acteristics of local vibrations based on the rapid convergence of the real part of the Green's function outside the band of quasi-continuous spectrum using the method of Jacobi matri‐ ces [16,17].

**Figure 1.** Electronic density of states of perfect graphene (curve 1) and the local density of states for an isolated boron

For a nitrogen impurity *ε*˜ −*ε*(*K*)≈ −0.525*J* and *η* ≈ −0.5 (according to [4]). As shown in [22],

The local density of states of the nitrogen substitutional impurity calculated in [11] has qua‐ si-local maxima in both intervals. For an boron substitutional impurity ( *η* ≈0.5 ) [6], quasilocalized states are absent in the −*ε*(*M* ), *ε*(*M* ) interval [12]. Figure 2 shows the graphical solution of the Lifshitz equation (2) for a given impurity atom. In this case the Lifshitz equa‐ tion has no solutions in interval −*ε*(*M* ), *ε*(*M* ) (corresponding dependences *S*(*ε*) are shown as curves 3 in Figure 2). The local Green's function of the boron impurity (curve 2 in Figure 1) has two poles outside the band of quasi-continuous spectrum, which are called local dis‐ crete levels and which are also solutions of equation (2). As is clearly seen in Figure 1 the area under the curve 2 is smaller (by the sum of the residues at these poles) than the area

Local discrete levels can be an important source of information about defective structure and force interactions in real crystals. To extract this useful information we should have ana‐ lytical expressions that relate the main characteristics of LDL (primarily their energy) to the

*<sup>ε</sup>*˜ <sup>+</sup> *εη*(2 <sup>+</sup> *<sup>η</sup>*) (3)

*<sup>S</sup>*(*ε*, *<sup>ε</sup>*˜, *<sup>η</sup>*)= (1 <sup>+</sup> *<sup>η</sup>*)<sup>2</sup>

equation (2) has a solution for both interval −*ε*(*M* ), *ε*(*K*) and interval *ε*(*K*), *ε*(*M* ) .

substitutional impurity (curve 2).

138 New Progress on Graphene Research

under the curve 1.

parameters of the defect and the host lattice.

Let us, briefly, to the extent necessary to understand the use of the classification of the eigen‐ functions of Hamiltonian (1), to present the basics of the method of Jacobi matrices. This method allows, without finding the dispersion laws, to calculate directly the local partial Green's functions of the system, corresponding to the perturbation of one or more atoms. This perturbation is described by the so-called generating vector *h* → <sup>0</sup>∈*H* , where *H* is the space of electronic excitations of atoms. Its dimension is *qN* , where *N* is the number of atoms in the system, and *q* is the dimension of the displacement of a single atom (*q = 1, 2, 3)*. Vectors of the space *H* are denoted by an arrow above the symbol, and "ordinary" *q*-dimen‐ sional vectors are in bold italics.

If, using the generating vector (*p* is the number of excited atoms) and the Hamiltonian (1), we construct the sequence {*H* ^ *nh* → 0}*<sup>n</sup>*=0 *<sup>∞</sup>* , then the linear envelope covering the vectors of this sequence forms, in the *H* space, a cyclic subspace invariant to the operator *H* ^ . This subspace contains, within itself, all the atomic displacements generated by the vector *h* → <sup>0</sup> . The corresponding partial Green's function is determined as a matrix element *G*00(*λ*)≡(*h* → <sup>0</sup>, *εI* ^ −*H* ^ <sup>−</sup><sup>1</sup> *h* → <sup>0</sup>) , where *<sup>λ</sup>* <sup>≡</sup>*<sup>ω</sup>* <sup>2</sup> is the eigenvalue. Quantity *ρ*(*λ*)≡*π* <sup>−</sup><sup>1</sup> Im*G*00(*λ*) is called the spectral density generated by the initial displacement *h* → 0 . In the basis {*h* → *<sup>n</sup>*}*n*=0 *∞* which is obtained by the orthonormalization of the sequence {*H* ^ *nh* → 0}*<sup>n</sup>*=0 *<sup>∞</sup>* , the operator (1) is represented in the form of a tridiagonal Jacobi matrix (or J-matrix). This matrix has a simple spectrum, what greatly simplifies finding the partial Green's functions and spectral densi‐ ties. As can be seen, this method does not use explicitly the translational symmetry of the crystal, making it extremely effective for treating systems in which such symmetry is bro‐ ken. The method of Jacobi matrices is particularly effective for treating systems with a sim‐ ply connected quasi-continuous band of spectrum *D*. In this case, with increasing rank of the J-matrix ( *n* →*∞* ), its diagonal elements *an* converge to *a* corresponding to the middle of the bandwidth *D*, and nondiagonal elements *b <sup>n</sup>* converge to *b* corresponding to the one-quarter of the bandwidth *D*.

In practice, it is usually possible to calculate the Jacobi matrix of the Hamiltonian of a finite

Electronic and Vibrational Properties of Adsorbed and Embedded Graphene and Bigraphene with Defects

is called analytical approximation of LGF. All dependences in Figure 1 and curves 1 and 2 in Figure 2 were calculated by the formula (8), using the Jacobi matrix of the Hamiltonian (1) with rank *n* =600 . If we count the energy from the Fermi energy level, then all diagonal ele‐ ments of Jacobi matrices are zero ( *an* =*a* =0; *b* =*ε*<sup>0</sup> / 2 , where *ε*<sup>0</sup> =3*J* is the half-width of the quasi-continuous spectrum). A good accuracy of the approximations shown in figures is confirmed also by the fact that they show the nonanalyticity effects corresponding to the densities of states of systems with the dimension larger than unity (so-called van Hove sin‐ gularities). In the vicinity of these singularities the expression (8) slowly converges to the

of the LGF calculated by formula (8) with *n*=1. As can be seen in the band of the quasi-con‐ tinuous spectrum these relationships have very little in common with curves 1 and 2. Thus,

singularities at the edges of the band of quasi-continuous spectrum, characteristic for the 2*D* systems are absent. However, outside the band of quasi-continuous spectrum (also in the

cide, and if we put in the Lifshitz equation (2) instead of LGF its approximation (8) for *n*=1, the obtained solutions give the energies of LDL with quite high accuracy. Moreover, these solutions can be easily found analytically. The LGF approximation by formula (8) for *n*=1 was named the approximation of two moments in [23]. Indeed, it follows from the orthonor‐

2

→ 0) =(*h* → 0, *εI* ^ −*H* ^ <sup>−</sup>*<sup>Λ</sup>* ^ <sup>−</sup><sup>1</sup> *h* → 0) . In

Finding the characteristics of LDL is more convenient without using equation (2), looking

the approximation of two moments for the subspace generated by the excitation of an im‐

area of intersection of the real part of LGF with curve (3)) curves 2 and 2′

(*ε*)*K∞*(*ε*)

in Figure 2 show local density of states and their corresponding real parts

does not even hint at the *V*-shaped "Dirac" singularity at *ε* =*ε*(*K*)=*εF* , and on

in the interval −*ε*(*M* ), *ε*(*M* ) both non-monotonous parts and the logarithmic

(*ε*)*K∞*(*ε*) (8)

http://dx.doi.org/10.5772/50562

141

practically coin‐

*Qn*(*ε*)−*bn*−1*Qn*−<sup>1</sup>

*Pn*(*ε*)−*bn*−1*Pn*−<sup>1</sup>

*G*(*ε*, *h* → 0) ≈

true values of the really and imaginary parts of LGF.

mality of the polynomials defined in [16,17] that

−*ε*<sup>0</sup>

(*<sup>ε</sup>* <sup>2</sup> <sup>−</sup>*a*<sup>0</sup> 2 )

*b*0

*ερ*(*ε*)*dε* =*M*1,

for them as the poles of the LGF perturbed Hamiltonian *<sup>G</sup>*˜(*<sup>ε</sup>*, *<sup>h</sup>*

<sup>2</sup> *ρ*(*ε*)*dε* =1⇒*b*<sup>0</sup> = *M*<sup>2</sup> −*M*<sup>1</sup>

−*ε*<sup>0</sup>

rank. The expression

Curves 1′

the curve 1′

the curve 2′

*∫* −*ε*<sup>0</sup>

*∫* −*ε*<sup>0</sup>

−*ε*<sup>0</sup> *P*1

−*ε*<sup>0</sup> *P*1 and 2′

(*ε*)*ρ*(*ε*)*dε* =0⇒*a*<sup>0</sup> = *∫*

−*ε*<sup>0</sup>

−*ε*<sup>0</sup>

2(*ε*)*ρ*(*ε*)*dε* = *∫*

purity atom, we get

For the local Green's function (LGF), corresponding to the excitations of one or more atoms, which are determined by the generating vector *h* → <sup>0</sup> , we get following expression using the *J*matrix method

$$G(\varepsilon, \vec{h}\_0) = \begin{pmatrix} \vec{h}\_0 \ \mathbf{f}\_\varepsilon \hat{\mathbf{I}} - \hat{H} \mathbf{J}^{-1} \vec{h}\_0 \end{pmatrix} = \lim\_{n \to \omega} \frac{Q\_n(\varepsilon) - b\_{n-1} Q\_{n-1}(\varepsilon) K\_n(\varepsilon)}{P\_n(\varepsilon) - b\_{n-1} P\_{n-1}(\varepsilon) K\_n(\varepsilon)} \tag{4}$$

where *I* ^ is the unit operator and polynomials *Pn*(*ε*) and *Qn*(*ε*) are determined by the follow‐ ing recurrence relations

$$b\_n\{P,Q\}\_{n+1}(\varepsilon) = (\varepsilon - a\_n)\{P,Q\}\_n(\varepsilon) - b\_{n-1}\{P,Q\}\_{n-1}(\varepsilon) \tag{5}$$

The initial conditions are *P*−<sup>1</sup> (*ε*)=*Q*<sup>0</sup> (*ε*)≡0 *, P*<sup>0</sup> (*ε*)≡1 *, Q*<sup>0</sup> (*ε*)≡0 , and *Q*<sup>1</sup> (*ε*)≡*b*<sup>0</sup> −1 , function *K∞*(*ε*) corresponds to the LGF operator, with all elements of its J-matrix being equal to their limit values *a* and *b*:

$$K\_{\circ}(\varepsilon) = 2b^{-2} \left[ \varepsilon - a - Z(\varepsilon) \cdot \sqrt{(\varepsilon - a + 2b)(\varepsilon - a - 2b)} \right] \tag{6}$$

$$\mathcal{Z}\left(\varepsilon\right) \equiv \Theta\left(a - \mathfrak{D}b - \varepsilon\right) + i\Theta\left(\varepsilon - a + \mathfrak{D}b\right) \cdot \Theta\left(a + \mathfrak{D}b - \varepsilon\right) - \Theta\left(\varepsilon - a - \mathfrak{D}b\right) \tag{7}$$

The method of Jacobi matrices can treat as a regular singular perturbation a much larger num‐ ber of perturbations of the phonon spectrum due to the presence of various crystal defects than the traditional methods [19, 20]. In addition, perturbations do not change the band‐ width of the quasi-continuous spectrum, and consequently, the asymptotic values of the ele‐ ments of the J-matrix can be regarded as an asymptotically degenerated regular perturbation [23]. This type of perturbations covers virtually all perturbations of the phonon spectrum caused by local defects. The calculation of vibration characteristics of such systems is per‐ formed, using the method of J-matrices, with the same accuracy as for the initial ideal system. In practice, it is usually possible to calculate the Jacobi matrix of the Hamiltonian of a finite rank. The expression

called the spectral density generated by the initial displacement *h*

represented in the form of a tridiagonal Jacobi matrix (or J-matrix). This matrix has a simple spectrum, what greatly simplifies finding the partial Green's functions and spectral densi‐ ties. As can be seen, this method does not use explicitly the translational symmetry of the crystal, making it extremely effective for treating systems in which such symmetry is bro‐ ken. The method of Jacobi matrices is particularly effective for treating systems with a sim‐ ply connected quasi-continuous band of spectrum *D*. In this case, with increasing rank of the J-matrix ( *n* →*∞* ), its diagonal elements *an* converge to *a* corresponding to the middle of the bandwidth *D*, and nondiagonal elements *b <sup>n</sup>* converge to *b* corresponding to the one-quarter

For the local Green's function (LGF), corresponding to the excitations of one or more atoms,

→

^ is the unit operator and polynomials *Pn*(*ε*) and *Qn*(*ε*) are determined by the follow‐

*<sup>n</sup>*(*ε*)−*bn*−<sup>1</sup>

(*ε*)≡1 *, Q*<sup>0</sup>

*Z*(*ε*)≡*Θ*(*a* −2*b* −*ε*) + *iΘ*(*ε* −*a* + 2*b*) ⋅*Θ*(*a* + 2*b* −*ε*)−*Θ*(*ε* −*a* −2*b*) (7)

*<sup>K</sup>∞*(*ε*)=2*<sup>b</sup>* <sup>−</sup><sup>2</sup> *<sup>ε</sup>* <sup>−</sup>*<sup>a</sup>* <sup>−</sup>*Z*(*ε*) <sup>⋅</sup> (*<sup>ε</sup>* <sup>−</sup>*<sup>a</sup>* <sup>+</sup> <sup>2</sup>*b*)(*<sup>ε</sup>* <sup>−</sup>*<sup>a</sup>* <sup>−</sup>2*b*) (6)

*Qn*(*ε*)−*bn*−1*Qn*−<sup>1</sup>

*Pn*(*ε*)−*bn*−1*Pn*−<sup>1</sup>

{*P*, *Q*} *n*−1

which is obtained by the orthonormalization of the sequence {*H*

which are determined by the generating vector *h*

*bn*{*P*, *Q*}

*n*+1

(*ε*)=(*ε* −*an*){*P*, *Q*}

(*ε*)≡0 *, P*<sup>0</sup>

*K∞*(*ε*) corresponds to the LGF operator, with all elements of its J-matrix being equal to their

The method of Jacobi matrices can treat as a regular singular perturbation a much larger num‐ ber of perturbations of the phonon spectrum due to the presence of various crystal defects than the traditional methods [19, 20]. In addition, perturbations do not change the band‐ width of the quasi-continuous spectrum, and consequently, the asymptotic values of the ele‐ ments of the J-matrix can be regarded as an asymptotically degenerated regular perturbation [23]. This type of perturbations covers virtually all perturbations of the phonon spectrum caused by local defects. The calculation of vibration characteristics of such systems is per‐ formed, using the method of J-matrices, with the same accuracy as for the initial ideal system.

(*ε*)=*Q*<sup>0</sup>

*G*(*ε*, *h* → 0) =(*h* → 0, *εI* ^ −*H* ^ <sup>−</sup><sup>1</sup> *h* → <sup>0</sup>) =lim *n*→*∞*

ing recurrence relations

limit values *a* and *b*:

The initial conditions are *P*−<sup>1</sup>

of the bandwidth *D*.

140 New Progress on Graphene Research

matrix method

where *I*

→

<sup>0</sup> , we get following expression using the *J*-

(*ε*)*K∞*(*ε*) (4)

(*ε*) (5)

(*ε*)≡*b*<sup>0</sup> −1

, function

(*ε*)*K∞*(*ε*)

(*ε*)≡0 , and *Q*<sup>1</sup>

^ *nh* → 0}*<sup>n</sup>*=0

0 . In the basis {*h*

*<sup>∞</sup>* , the operator (1) is

→ *<sup>n</sup>*}*n*=0 *∞*

$$G(\varepsilon, \vec{h}\_0) \approx \frac{Q\_n(\varepsilon) - b\_{n-1} Q\_{n-1}(\varepsilon) K\_\infty(\varepsilon)}{P\_n(\varepsilon) - b\_{n-1} P\_{n-1}(\varepsilon) K\_\infty(\varepsilon)}\tag{8}$$

is called analytical approximation of LGF. All dependences in Figure 1 and curves 1 and 2 in Figure 2 were calculated by the formula (8), using the Jacobi matrix of the Hamiltonian (1) with rank *n* =600 . If we count the energy from the Fermi energy level, then all diagonal ele‐ ments of Jacobi matrices are zero ( *an* =*a* =0; *b* =*ε*<sup>0</sup> / 2 , where *ε*<sup>0</sup> =3*J* is the half-width of the quasi-continuous spectrum). A good accuracy of the approximations shown in figures is confirmed also by the fact that they show the nonanalyticity effects corresponding to the densities of states of systems with the dimension larger than unity (so-called van Hove sin‐ gularities). In the vicinity of these singularities the expression (8) slowly converges to the true values of the really and imaginary parts of LGF.

Curves 1′ and 2′ in Figure 2 show local density of states and their corresponding real parts of the LGF calculated by formula (8) with *n*=1. As can be seen in the band of the quasi-con‐ tinuous spectrum these relationships have very little in common with curves 1 and 2. Thus, the curve 1′ does not even hint at the *V*-shaped "Dirac" singularity at *ε* =*ε*(*K*)=*εF* , and on the curve 2′ in the interval −*ε*(*M* ), *ε*(*M* ) both non-monotonous parts and the logarithmic singularities at the edges of the band of quasi-continuous spectrum, characteristic for the 2*D* systems are absent. However, outside the band of quasi-continuous spectrum (also in the area of intersection of the real part of LGF with curve (3)) curves 2 and 2′ practically coin‐ cide, and if we put in the Lifshitz equation (2) instead of LGF its approximation (8) for *n*=1, the obtained solutions give the energies of LDL with quite high accuracy. Moreover, these solutions can be easily found analytically. The LGF approximation by formula (8) for *n*=1 was named the approximation of two moments in [23]. Indeed, it follows from the orthonor‐ mality of the polynomials defined in [16,17] that

$$\begin{aligned} \int\_{-\epsilon\_0}^{-\epsilon\_0} P\_1(\epsilon) \rho(\varepsilon) d\varepsilon &= 0 \Rightarrow a\_0 = \int\_{-\epsilon\_0}^{-\epsilon\_0} \varepsilon \rho(\varepsilon) d\varepsilon = M\_1, \\ \int\_{-\epsilon\_0}^{-\epsilon\_0} P\_1^2(\varepsilon) \rho(\varepsilon) d\varepsilon &= \int\_{-\epsilon\_0}^{-\epsilon\_0} \frac{\left(\varepsilon^{-2} - a\_0^2\right)}{b\_0^2} \rho(\varepsilon) d\varepsilon = 1 \Rightarrow b\_0 = \sqrt{M\_2 - M\_1^2} \end{aligned}$$

Finding the characteristics of LDL is more convenient without using equation (2), looking for them as the poles of the LGF perturbed Hamiltonian *<sup>G</sup>*˜(*<sup>ε</sup>*, *<sup>h</sup>* → 0) =(*h* → 0, *εI* ^ −*H* ^ <sup>−</sup>*<sup>Λ</sup>* ^ <sup>−</sup><sup>1</sup> *h* → 0) . In the approximation of two moments for the subspace generated by the excitation of an im‐ purity atom, we get

$$\widetilde{G}(\varepsilon) = \frac{1}{\varepsilon - a\_0 - b\_0^2 K\_{\text{ov}}(\varepsilon)}\tag{9}$$

*Gmn*(*ε*, *h* → 0) =(*h* → *<sup>m</sup>*, *εI* ^ −*H* ^ <sup>−</sup><sup>1</sup> *h* →

*μn* (±) =*Pn* 2 (*εd* (±) ) ⋅*μ*<sup>0</sup> (±)

The intensities *μn*

spectrum.

(in this case *ε*˜*ε*<sup>0</sup>

ure 2 merge at |*ε* | >*ε*0 .

ed by arrows.

(±)

*<sup>n</sup>*) = −*Pm*(*ε*)*Qn*(*ε*) + *Pm*(*ε*)*Pm*(*ε*)*G*(*ε*, *h*

This implies that the damping of LDL, i.e. the decay of its intensity with the increasing dis‐ tance from the impurity atom (i.e. with the increase of *n*) follows the equation

*a*0

decay with increasing *n* according to *μn*<sup>0</sup>

from *n*=1 they form an infinitely decreasing geometric progression whose denominator

<sup>2</sup> <sup>+</sup> (2*<sup>γ</sup>* <sup>−</sup>1)*ε*<sup>0</sup>

*n*=0 *∞ μn* (+) =∑ *n*=0 *∞ μn* (−)

LDL is the formation of one quasi-particle outside the band of the quasi-continuous electron

Formulas (10, 12-14, 16) give simple analytical expressions of the local conditions of the exis‐

(in this case *ε*˜*ε*<sup>0</sup>

It is seen that such levels exist in a very wide range of variables *ε*˜ and *η* . The absence of LDL is possible only in a narrow range of values *ε*˜ at *η* 3 / 2−1 . If *η* 3 / 2−1 at least one local discrete level exists. In fact, the lines delineating the area of existence of LDL in graphene in the plane {*ε*˜, *η*} must pass through the origin of coordinates, since Re*G*(*ε*)→*∞* for *ε* → ± *ε*<sup>0</sup> . However, since this divergence is logarithmic, for any appreciable splitting of LDL from the

boundary of quasi-continuous spectrum there is a certain threshold. Curves 2 and 2′

(*<sup>γ</sup>* <sup>−</sup>1)=*ε*<sup>0</sup> 2(1 <sup>+</sup> *<sup>η</sup>*)<sup>2</sup> / <sup>3</sup>−1 ) lie above and under curves in Figure 3, as indicat‐

tence of discrete levels due to the presence of a substitutional impurity in graphene.

(+)

*ε*0 (2*γ* −1)

<sup>2</sup> <sup>∓</sup> *<sup>a</sup>*<sup>0</sup>

2

*Pn*(*ε<sup>d</sup>* (±)

) = 2*γ* ⋅ ±

*q* (±)

Summing these progressions we see that ∑

Regions of the LDL existence for *ε<sup>d</sup>*

<sup>=</sup> *<sup>a</sup>*<sup>0</sup>

. Using the method of mathematical induction we can prove that

Electronic and Vibrational Properties of Adsorbed and Embedded Graphene and Bigraphene with Defects

<sup>2</sup> <sup>+</sup> (2*<sup>γ</sup>* <sup>−</sup>1)*ε*<sup>0</sup>

*ε*0 (2*γ* −1)

<sup>2</sup> <sup>∓</sup> *<sup>a</sup>*<sup>0</sup>

*n*

(±)

=2*<sup>γ</sup>* <sup>⋅</sup> *<sup>q</sup> <sup>n</sup>* <sup>⋅</sup>*μ*<sup>0</sup>

(±)

=1 , that is the formation of each

(1−*γ*)=*ε*<sup>0</sup> <sup>1</sup>−2(1 <sup>+</sup> *<sup>η</sup>*)<sup>2</sup> / <sup>3</sup> ) and for *ε<sup>d</sup>*

→

0) (*m*<*n*).

http://dx.doi.org/10.5772/50562

(15)

143

(16)

(−)

in Fig‐

, that is, starting

For the case of an isolated substitutional impurity

$$a\_0 = \tilde{\varepsilon}; b\_0 = \sqrt{3}(1 + \eta)f = \frac{1 + \eta}{\sqrt{3}}\varepsilon\_0 \tag{10}$$

from where we get

$$\widetilde{G}(\varepsilon) = \frac{(1-\gamma)\varepsilon - a\_0 + Z(\varepsilon)\gamma\sqrt{\left|\varepsilon\right|^2 - \varepsilon\_0^2\left|}}{R(\varepsilon)}\tag{11}$$

where *R*(*ε*)=(1−2*γ*)*<sup>ε</sup>* <sup>2</sup> <sup>−</sup>2(1−*γ*)*a*0*<sup>ε</sup>* <sup>+</sup> *<sup>a</sup>*<sup>0</sup> <sup>2</sup> + *γ* <sup>2</sup> *ε*0 2 and

$$\gamma \equiv b\_0^2 \mid 2b^2 = 2(1+\eta)^2/3\tag{12}$$

Local discrete levels are poles (11), i.e. the roots of *R*(*ε*) are

$$
\varepsilon\_d^{(\pm)} = \frac{(\gamma - 1)a\_0 \pm \gamma \sqrt{a\_0^2 + (2\gamma - 1)\varepsilon\_0^2}}{2\gamma - 1} \tag{13}
$$

Residues at these poles *μ*<sup>0</sup> (±) =*r e* ‵ *s ε*=*ε<sup>d</sup>* (±) *<sup>G</sup>*˜(*ε*) are called intensities of LDL and they determine the relative LDL "amplitude" on the impurity atom: *μ*<sup>0</sup> (±) =1−*π* <sup>−</sup><sup>1</sup> *∫* −*ε*<sup>0</sup> *ε*0 Im*<sup>G</sup>*˜(*ε*)*d<sup>ε</sup>* . The condition that the intensity differs from zero defines the existence region of LDL. In this case

$$\begin{split} \mu\_{0}^{(\star)} &= \frac{\gamma \, a\_{0} + (\gamma - 1) \sqrt{a\_{0}^{2} + (2\gamma - 1) \varepsilon\_{0}^{2}}}{(2\gamma - 1) \sqrt{a\_{0}^{2} + (2\gamma - 1) \varepsilon\_{0}^{2}}} \cdot \Theta \Big( \frac{a\_{0}}{\varepsilon\_{0}} - 1 + \gamma \Big); \\ \mu\_{0}^{(\star)} &= \frac{-\gamma \, a\_{0} + (\gamma - 1) \sqrt{a\_{0}^{2} + (2\gamma - 1) \varepsilon\_{0}^{2}}}{(2\gamma - 1) \sqrt{a\_{0}^{2} + (2\gamma - 1) \varepsilon\_{0}^{2}}} \cdot \Theta \Big( \gamma - 1 - \frac{a\_{0}}{\varepsilon\_{0}} \Big). \end{split} \tag{14}$$

It was shown in [23] that

Electronic and Vibrational Properties of Adsorbed and Embedded Graphene and Bigraphene with Defects http://dx.doi.org/10.5772/50562 143

$$G\_{mn}(\varepsilon,\ \vec{h}\_{\ \!\!\/)} = \left(\vec{h}\_{\ \!\!\/ \!\/ \!\/ \!\/}\,\mathsf{T}\,\hat{\mathsf{I}} - \hat{\mathsf{H}}\,\Box^{-1}\vec{h}\_{\ \!\!\/ \!\/}\right) = -P\_m(\varepsilon)Q\_n(\varepsilon) + P\_m(\varepsilon)P\_m(\varepsilon)G\{\varepsilon,\ \vec{h}\_{\ \!\!\/)}\tag{m$$

*<sup>G</sup>*˜(*ε*)= <sup>1</sup> *ε* −*a*<sup>0</sup> −*b*<sup>0</sup> 2

*<sup>a</sup>*<sup>0</sup> <sup>=</sup>*ε*˜; *<sup>b</sup>*<sup>0</sup> <sup>=</sup> 3(1 <sup>+</sup> *<sup>η</sup>*)*<sup>J</sup>* <sup>=</sup> <sup>1</sup> <sup>+</sup> *<sup>η</sup>*

(1−*γ*)*<sup>ε</sup>* <sup>−</sup>*a*<sup>0</sup> <sup>+</sup> *<sup>Z</sup>*(*ε*)*<sup>γ</sup>* <sup>|</sup>*<sup>ε</sup>* <sup>2</sup> <sup>−</sup>*ε*<sup>0</sup>

<sup>2</sup> + *γ* <sup>2</sup> *ε*0 2 and

(*γ* −1)*a*<sup>0</sup> ± *γ a*<sup>0</sup>

the intensity differs from zero defines the existence region of LDL. In this case

*γa*<sup>0</sup> + (*γ* −1) *a*<sup>0</sup>

(2*γ* −1) *a*<sup>0</sup>

<sup>=</sup> <sup>−</sup>*γa*<sup>0</sup> <sup>+</sup> (*<sup>γ</sup>* <sup>−</sup>1) *<sup>a</sup>*<sup>0</sup>

(2*γ* −1) *a*<sup>0</sup>

2*γ* −1

<sup>2</sup> <sup>+</sup> (2*<sup>γ</sup>* <sup>−</sup>1)*ε*<sup>0</sup>

<sup>2</sup> <sup>+</sup> (2*<sup>γ</sup>* <sup>−</sup>1)*ε*<sup>0</sup>

<sup>2</sup> <sup>+</sup> (2*<sup>γ</sup>* <sup>−</sup>1)*ε*<sup>0</sup>

<sup>2</sup> <sup>+</sup> (2*<sup>γ</sup>* <sup>−</sup>1)*ε*<sup>0</sup>

*γ* ≡*b*<sup>0</sup>

Local discrete levels are poles (11), i.e. the roots of *R*(*ε*) are

*εd* (±) =

(±) =*r e* ‵ *s ε*=*ε<sup>d</sup>* (±)

relative LDL "amplitude" on the impurity atom: *μ*<sup>0</sup>

*μ*0 (+) =

*μ*0 (−)

It was shown in [23] that

3

<sup>2</sup> <sup>+</sup> (2*<sup>γ</sup>* <sup>−</sup>1)*ε*<sup>0</sup>

(±)

2

<sup>2</sup> <sup>⋅</sup>*Θ*( *<sup>a</sup>*<sup>0</sup> *ε*0

2

<sup>2</sup> <sup>⋅</sup>*Θ*(*<sup>γ</sup>* <sup>−</sup>1<sup>−</sup>

2

=1−*π* <sup>−</sup><sup>1</sup> *∫*

2 |

For the case of an isolated substitutional impurity

*<sup>G</sup>*˜(*ε*)=

where *R*(*ε*)=(1−2*γ*)*<sup>ε</sup>* <sup>2</sup> <sup>−</sup>2(1−*γ*)*a*0*<sup>ε</sup>* <sup>+</sup> *<sup>a</sup>*<sup>0</sup>

Residues at these poles *μ*<sup>0</sup>

from where we get

142 New Progress on Graphene Research

*<sup>K</sup>∞*(*ε*) (9)

*<sup>R</sup>*(*ε*) (11)

<sup>2</sup> / 2*b* <sup>2</sup> =2(1 + *η*)<sup>2</sup> / 3 (12)

*<sup>G</sup>*˜(*ε*) are called intensities of LDL and they determine the

Im*<sup>G</sup>*˜(*ε*)*d<sup>ε</sup>* . The condition that

−*ε*<sup>0</sup>

−1 + *γ*);

*a*0 *ε*0 ).

*ε*0

*ε*<sup>0</sup> (10)

(13)

(14)

This implies that the damping of LDL, i.e. the decay of its intensity with the increasing dis‐ tance from the impurity atom (i.e. with the increase of *n*) follows the equation *μn* (±) =*Pn* 2 (*εd* (±) ) ⋅*μ*<sup>0</sup> (±) . Using the method of mathematical induction we can prove that

$$P\_n(\varepsilon\_d^{(\pm)}) = \sqrt{2\gamma} \cdot \left[ \pm \frac{\sqrt{a\_0^2 + (2\gamma - 1)\varepsilon\_0^2} \mp a\_0}{\varepsilon\_0(2\gamma - 1)} \right]^n \tag{15}$$

The intensities *μn* (±) decay with increasing *n* according to *μn*<sup>0</sup> (±) =2*<sup>γ</sup>* <sup>⋅</sup> *<sup>q</sup> <sup>n</sup>* <sup>⋅</sup>*μ*<sup>0</sup> (±) , that is, starting from *n*=1 they form an infinitely decreasing geometric progression whose denominator

$$q^{(\omega)} = \left[\frac{\sqrt{a\_0^2 + (2\gamma - 1)\varepsilon\_0^2} \mp a\_0}{\varepsilon\_0(2\gamma - 1)}\right]^2 \tag{16}$$

Summing these progressions we see that ∑ *n*=0 *∞ μn* (+) =∑ *n*=0 *∞ μn* (−) =1 , that is the formation of each LDL is the formation of one quasi-particle outside the band of the quasi-continuous electron spectrum.

Formulas (10, 12-14, 16) give simple analytical expressions of the local conditions of the exis‐ tence of discrete levels due to the presence of a substitutional impurity in graphene.

Regions of the LDL existence for *ε<sup>d</sup>* (+) (in this case *ε*˜*ε*<sup>0</sup> (1−*γ*)=*ε*<sup>0</sup> <sup>1</sup>−2(1 <sup>+</sup> *<sup>η</sup>*)<sup>2</sup> / <sup>3</sup> ) and for *ε<sup>d</sup>* (−) (in this case *ε*˜*ε*<sup>0</sup> (*<sup>γ</sup>* <sup>−</sup>1)=*ε*<sup>0</sup> 2(1 <sup>+</sup> *<sup>η</sup>*)<sup>2</sup> / <sup>3</sup>−1 ) lie above and under curves in Figure 3, as indicat‐ ed by arrows.

It is seen that such levels exist in a very wide range of variables *ε*˜ and *η* . The absence of LDL is possible only in a narrow range of values *ε*˜ at *η* 3 / 2−1 . If *η* 3 / 2−1 at least one local discrete level exists. In fact, the lines delineating the area of existence of LDL in graphene in the plane {*ε*˜, *η*} must pass through the origin of coordinates, since Re*G*(*ε*)→*∞* for *ε* → ± *ε*<sup>0</sup> . However, since this divergence is logarithmic, for any appreciable splitting of LDL from the boundary of quasi-continuous spectrum there is a certain threshold. Curves 2 and 2′ in Fig‐ ure 2 merge at |*ε* | >*ε*0 .

So, a good agreement between the results of numerical calculation of the LDL characteris‐ tics using Jacobi matrices of high rank, and their analytical description by the Green's func‐ tion (9), which relates these characteristics to the parameters of the defect (13) (14), makes it relatively easy to extract the defect parameters from the known characteristics of LDL. Exper‐

Electronic and Vibrational Properties of Adsorbed and Embedded Graphene and Bigraphene with Defects

ing (10), (12) and (13), to the determination of the parameters *ε*˜ and *η* and this might represent a significant advance in creating nanomaterials with predetermined spectral characteristics.

is, the impurity levels can not be completely localized on the impurities, but they also ap‐ pear in the spectra of surrounding carbon atoms. This greatly increases the probability of ex‐

Bilayer graphene is a carbon film consisting of two graphene monolayers, linked together by (as in bulk graphite) van der Waals forces. Since the distance between the layers (film thick‐ ness) *h* ~3.5 Å the bilayer graphene can be considered not a nanofilm but a subnanofilm. The constants of the interatomic interaction of bilayer graphene were determined and its phonon density of states and partial contributions to this quantity from the atomic displacements along different crystallographic directions were calculated [7,24]. On the basis of the analy‐ ses of the mean-square amplitudes of atomic displacements calculated using data from the spectral densities, we have shown that the flat shape of a free bilayer graphene remains sta‐ ble up to the temperatures much higher than the room temperature, which makes this com‐ pound promising for nanoelectronics. In this section we calculate and analyze the electronic spectrum of a defect-free bilayer graphene. Naturally of greatest interest is its behavior in the energy range close to *εF* where there are characteristic Dirac points on the spectrum of

The unit cell of graphene contains two physically equivalent atoms and therefore local Green's function and the local density of states (LDOS) of the atoms of different sublattices are identical. On the other hand, bilayer graphene unit cell contains four atoms, and atoms of different sublattices of a single graphene layer interact differently with the atoms of the

The electronic spectrum of bilayer graphene, as well as the electronic spectrum of graphene can be described in the strong coupling approximation. Corresponding Hamiltonian has form (1). For graphene and bigraphene we assume that the electron hopping within the lay‐ er is possible only between nearest neighbors *Jij* = *J* ≈2.8 eV (see for example [25]). Electron hopping between layers is also assumed to be possible only between nearest neighbors from different layers, that is, between those which lie at a distance *h* from each other. Denote the

such neighbors (sublattice *AI* and *AII*, see Figure 1). In sublattices *BI* and *BII* no such neigh‐

. Note that only half of the bilayer graphene atoms have

(±)

http://dx.doi.org/10.5772/50562

(+) →*μ*<sup>0</sup> (−)

should lead, us‐

145

→1 / 2 . That

imental measurement (e.g. by scanning tunneling microscopy) of values *ε<sup>d</sup>*

As can be seen from Figure 4, with increasing *η* the intensity of LDL *μ*<sup>0</sup>

**3. The electronic spectrum of bilayer graphene**

graphene monolayer (whose plane shape is unstable).

corresponding hopping integral *J* ′

other layer and their physical equivalence is disrupted (Figure 5a).

perimental detection of such levels, even at low concentrations of impurities.

**Figure 3.** Regions of the existence of discrete levels for substitutional impurities in graphene.

Figure 4 shows, for *ε*˜ ≈0.525*J* (for boron, from [6]), the dependences of energies, the LDL in‐ tensities at the boron impurity and the damping parameters of the value *η* that characterizes the change in the overlap integral of the boron impurity atom (10). Solid lines show the char‐ acteristics of LDL, calculated using the approximation of the Green's function (9), i.e. accord‐ ing to the analytical formulas (13), (14) and (16). Open circles show the results of numerical calculations of dependences *ε<sup>d</sup>* (±) (*η*) and *μ*<sup>0</sup> (±) (*η*) , using the Green's function in the form of (8), calculated by the Jacobi matrices of the *n* ≥100 rank.

**Figure 4.** The basic characteristics of LDL in the presence of a boron substitutional impurity in graphene.

It is seen that at the threshold values *ε<sup>d</sup>* (±) = ± *ε*0 of the LDL formation, the intensities of LDL equal zero and the parameters of damping are equal to unit. Further increase of |*ε<sup>d</sup>* (±) | is accompanied by the increase of *μ*<sup>0</sup> (±) and by the strengthening of the containment level.

So, a good agreement between the results of numerical calculation of the LDL characteris‐ tics using Jacobi matrices of high rank, and their analytical description by the Green's func‐ tion (9), which relates these characteristics to the parameters of the defect (13) (14), makes it relatively easy to extract the defect parameters from the known characteristics of LDL. Exper‐ imental measurement (e.g. by scanning tunneling microscopy) of values *ε<sup>d</sup>* (±) should lead, us‐ ing (10), (12) and (13), to the determination of the parameters *ε*˜ and *η* and this might represent a significant advance in creating nanomaterials with predetermined spectral characteristics. As can be seen from Figure 4, with increasing *η* the intensity of LDL *μ*<sup>0</sup> (+) →*μ*<sup>0</sup> (−) →1 / 2 . That is, the impurity levels can not be completely localized on the impurities, but they also ap‐ pear in the spectra of surrounding carbon atoms. This greatly increases the probability of ex‐ perimental detection of such levels, even at low concentrations of impurities.

#### **3. The electronic spectrum of bilayer graphene**

**Figure 3.** Regions of the existence of discrete levels for substitutional impurities in graphene.

(*η*) and *μ*<sup>0</sup>

(±)

**Figure 4.** The basic characteristics of LDL in the presence of a boron substitutional impurity in graphene.

(±)

(±)

equal zero and the parameters of damping are equal to unit. Further increase of |*ε<sup>d</sup>*

(±)

calculated by the Jacobi matrices of the *n* ≥100 rank.

calculations of dependences *ε<sup>d</sup>*

144 New Progress on Graphene Research

It is seen that at the threshold values *ε<sup>d</sup>*

accompanied by the increase of *μ*<sup>0</sup>

Figure 4 shows, for *ε*˜ ≈0.525*J* (for boron, from [6]), the dependences of energies, the LDL in‐ tensities at the boron impurity and the damping parameters of the value *η* that characterizes the change in the overlap integral of the boron impurity atom (10). Solid lines show the char‐ acteristics of LDL, calculated using the approximation of the Green's function (9), i.e. accord‐ ing to the analytical formulas (13), (14) and (16). Open circles show the results of numerical

(*η*) , using the Green's function in the form of (8),

= ± *ε*0 of the LDL formation, the intensities of LDL

and by the strengthening of the containment level.

(±) | is Bilayer graphene is a carbon film consisting of two graphene monolayers, linked together by (as in bulk graphite) van der Waals forces. Since the distance between the layers (film thick‐ ness) *h* ~3.5 Å the bilayer graphene can be considered not a nanofilm but a subnanofilm. The constants of the interatomic interaction of bilayer graphene were determined and its phonon density of states and partial contributions to this quantity from the atomic displacements along different crystallographic directions were calculated [7,24]. On the basis of the analy‐ ses of the mean-square amplitudes of atomic displacements calculated using data from the spectral densities, we have shown that the flat shape of a free bilayer graphene remains sta‐ ble up to the temperatures much higher than the room temperature, which makes this com‐ pound promising for nanoelectronics. In this section we calculate and analyze the electronic spectrum of a defect-free bilayer graphene. Naturally of greatest interest is its behavior in the energy range close to *εF* where there are characteristic Dirac points on the spectrum of graphene monolayer (whose plane shape is unstable).

The unit cell of graphene contains two physically equivalent atoms and therefore local Green's function and the local density of states (LDOS) of the atoms of different sublattices are identical. On the other hand, bilayer graphene unit cell contains four atoms, and atoms of different sublattices of a single graphene layer interact differently with the atoms of the other layer and their physical equivalence is disrupted (Figure 5a).

The electronic spectrum of bilayer graphene, as well as the electronic spectrum of graphene can be described in the strong coupling approximation. Corresponding Hamiltonian has form (1). For graphene and bigraphene we assume that the electron hopping within the lay‐ er is possible only between nearest neighbors *Jij* = *J* ≈2.8 eV (see for example [25]). Electron hopping between layers is also assumed to be possible only between nearest neighbors from different layers, that is, between those which lie at a distance *h* from each other. Denote the corresponding hopping integral *J* ′ . Note that only half of the bilayer graphene atoms have such neighbors (sublattice *AI* and *AII*, see Figure 1). In sublattices *BI* and *BII* no such neigh‐ bors exist, since nearest neighbors from different layers are at a distance *h* <sup>2</sup> + *a* <sup>2</sup> ( *a* ≈1.415 Å is the distance between nearest neighbors in the layer plane). Since this distance is only by less than 10% greater than *h* , we can neglect the interaction with the atoms of the sublattices *BI* and *BII* and this does not lead to qualitative changes in the behavior of the spectra near *ε<sup>F</sup>* (see, e.g. [6]). Then the dispersion relation of each of the four branches of the electronic spec‐ trum of bilayer graphene can be written as

$$\begin{aligned} \varepsilon\_{1,2}(k) &= \pm \sqrt{\varepsilon\_0^2(k) + \frac{\overline{I}^{\prime\prime}}{2} - \overline{I}^{\prime\prime} \overline{\left(\varepsilon\_0^2(k) + \frac{\overline{I}^{\prime\prime}}{4}\right)}};\\ \varepsilon\_{3,4}(k) &= \pm \sqrt{\varepsilon\_0^2(k) + \frac{\overline{I}^{\prime\prime}}{2} + \overline{I}^{\prime\prime} \overline{\left(\varepsilon\_0^2(k) + \frac{\overline{I}^{\prime\prime}}{4}\right)}} \end{aligned} \tag{17}$$

Dispersion curves along highly symmetric directions ΓK, ΓM and KM for *J'* = 0.1*J* are shown in Figure 6. Region near the K-point *Κ* =(0,±4*π* / 3*a* 3,0) and *Κ* =(±2*π* / 3*a*, ±2*π* / 3*a* 3,0) is shown in the inset. The same inset shows also the dispersion curves for graphene (18). We clearly see the quasi-relativistic nature of the electronic spectrum of graphene as well as or‐ dinary quadratic dispersion curves *ε*1,2(*k*) for bigraphene. Spectral branches

Electronic and Vibrational Properties of Adsorbed and Embedded Graphene and Bigraphene with Defects

, *J* ′

ar (relativistic) dispersion relation. For electronic modes of bilayer graphene *ε*1,2(*k*) near the

2(*Κ* + *κ*))

2(*<sup>Κ</sup>* <sup>+</sup> *<sup>κ</sup>*) <sup>2</sup> } <sup>=</sup> *<sup>ε</sup>*<sup>0</sup>

9*J* <sup>2</sup> *a* 2 *κ* 2 4*J*<sup>1</sup>

2 −*ε*<sup>0</sup>

> *J*1 <sup>2</sup> <sup>+</sup> <sup>2</sup>*ε*<sup>0</sup>

4(*Κ* + *κ*)≈

4(*Κ* + *κ*)

*J*1 2 <sup>2</sup> <sup>+</sup> *<sup>ε</sup>*<sup>0</sup>

2(*<sup>Κ</sup>* <sup>+</sup> *<sup>κ</sup>*) <sup>≈</sup> <sup>81</sup>*<sup>J</sup>* <sup>4</sup>

<sup>2</sup> ) and putting *<sup>k</sup>* <sup>=</sup>*<sup>Κ</sup>* <sup>+</sup> *<sup>κ</sup>* ( *κ*1 ) we find *ε*<sup>0</sup>

) interval. Indeed, if *k* takes value along ΓK,

(*Κ* + *κ*)≈ ∓

2(*Κ* + *κ*)−

*a* 4 *κ* 4

16*J*<sup>1</sup> <sup>2</sup> ,

3*aκ*

http://dx.doi.org/10.5772/50562

<sup>2</sup> , i.e. a line‐

(19)

147

(20)

*<sup>ε</sup>*3,4(*k*)∉( <sup>−</sup> *<sup>J</sup>* ′

then *ε*<sup>0</sup>

, *J* ′

(*k*)= ± (1 + 2cos

*Κ* point, we can write

<sup>2</sup> (*Κ* + *κ*)≈

*J*1 2 <sup>2</sup> <sup>+</sup> *<sup>ε</sup>*<sup>0</sup>

*ε*1,2

<sup>−</sup> *<sup>J</sup>*<sup>1</sup> 2 <sup>2</sup> <sup>+</sup> *<sup>ε</sup>*<sup>0</sup> ) are determined in the ( − *J* ′

2(*<sup>Κ</sup>* <sup>+</sup> *<sup>κ</sup>*)<sup>−</sup> ( *<sup>J</sup>*<sup>1</sup>

2 *J*<sup>1</sup>

**Figure 6.** Dispersion curves along high-symmetry directions of bilayer graphene.

2 <sup>2</sup> <sup>+</sup> *<sup>ε</sup>*<sup>0</sup>

<sup>2</sup> / <sup>2</sup> <sup>+</sup> *<sup>ε</sup>*<sup>0</sup>

*ε*1,2(*Κ* + *κ*)≈ ±

4(*Κ* + *κ*)

*ak* 3

2(*<sup>Κ</sup>* <sup>+</sup> *<sup>κ</sup>*) <sup>⋅</sup> {1<sup>−</sup> *<sup>ε</sup>*<sup>0</sup>

that is, we get an ordinary quadratic dispersion law

where *ε*<sup>0</sup> (*k*) is the electronic spectrum of graphene, calculated in the strong coupling approx‐ imation:

$$\varepsilon\_0(k) = \pm J \sqrt{1 + 4\cos\left|k \cdot \frac{a\_1 \cdot a\_2}{2}\right|} \left[\cos\left|k \cdot \frac{a\_1 + a\_2}{2}\right) + \cos\left|k \cdot \frac{a\_1 \cdot a\_2}{2}\right|\right] \tag{18}$$

where *a*<sup>1</sup> =( *<sup>a</sup>* <sup>2</sup> , *<sup>a</sup>* <sup>3</sup> <sup>2</sup> ,0) and *a*<sup>2</sup> =( *<sup>a</sup>* <sup>2</sup> , <sup>−</sup> *<sup>a</sup>* <sup>3</sup> <sup>2</sup> ,0) are two-dimensional Bravais lattice vectors (see Figure 5b).

**Figure 5.** The structure of bilayer graphene (a); Bravais lattice and the first Brillouin zone (b).

Dispersion curves along highly symmetric directions ΓK, ΓM and KM for *J'* = 0.1*J* are shown in Figure 6. Region near the K-point *Κ* =(0,±4*π* / 3*a* 3,0) and *Κ* =(±2*π* / 3*a*, ±2*π* / 3*a* 3,0) is shown in the inset. The same inset shows also the dispersion curves for graphene (18). We clearly see the quasi-relativistic nature of the electronic spectrum of graphene as well as or‐ dinary quadratic dispersion curves *ε*1,2(*k*) for bigraphene. Spectral branches *<sup>ε</sup>*3,4(*k*)∉( <sup>−</sup> *<sup>J</sup>* ′ , *J* ′ ) are determined in the ( − *J* ′ , *J* ′ ) interval. Indeed, if *k* takes value along ΓK, then *ε*<sup>0</sup> (*k*)= ± (1 + 2cos *ak* 3 <sup>2</sup> ) and putting *<sup>k</sup>* <sup>=</sup>*<sup>Κ</sup>* <sup>+</sup> *<sup>κ</sup>* ( *κ*1 ) we find *ε*<sup>0</sup> (*Κ* + *κ*)≈ ∓ 3*aκ* <sup>2</sup> , i.e. a line‐ ar (relativistic) dispersion relation. For electronic modes of bilayer graphene *ε*1,2(*k*) near the *Κ* point, we can write

$$\begin{split} \varepsilon\_{1,2}^{2}(\mathbf{K}+\boldsymbol{\kappa}) &\approx \frac{\int\_{1}^{2} + \varepsilon\_{0}^{2}(\mathbf{K}+\boldsymbol{\kappa})}{\mathfrak{D}} - 4\left\{ \left[\frac{\boldsymbol{f}\_{1}^{2}}{2} + \varepsilon\_{0}^{2}(\mathbf{K}+\boldsymbol{\kappa})\right]^{2} - \varepsilon\_{0}^{4}(\mathbf{K}+\boldsymbol{\kappa}) \approx \frac{\int\_{1}^{2}}{2} + \varepsilon\_{0}^{2}(\mathbf{K}+\boldsymbol{\kappa}) - \\ &- \left[\frac{\int\_{1}^{2}}{2} + \varepsilon\_{0}^{2}(\mathbf{K}+\boldsymbol{\kappa})\right] \cdot \left[1 - \frac{\varepsilon\_{0}^{4}(\mathbf{K}+\boldsymbol{\kappa})}{2\left[\int\_{1}^{2}/2 + \varepsilon\_{0}^{2}(\mathbf{K}+\boldsymbol{\kappa})\right]^{2}}\right] = \frac{\varepsilon\_{0}^{4}(\mathbf{K}+\boldsymbol{\kappa})}{\int\_{1}^{2} + 2\varepsilon\_{0}^{2}(\mathbf{K}+\boldsymbol{\kappa})} \approx \frac{81f^{4}\operatorname{a}^{4}\operatorname{\kappa}^{4}}{16f\_{1}^{2}}, \end{split} \tag{19}$$

that is, we get an ordinary quadratic dispersion law

bors exist, since nearest neighbors from different layers are at a distance *h* <sup>2</sup> + *a* <sup>2</sup>

2(*k*) <sup>+</sup> *<sup>J</sup>* ′<sup>2</sup>

2(*k*) <sup>+</sup> *<sup>J</sup>* ′<sup>2</sup>

*a*<sup>1</sup> - *a*<sup>2</sup>

<sup>2</sup> , <sup>−</sup> *<sup>a</sup>* <sup>3</sup>

**Figure 5.** The structure of bilayer graphene (a); Bravais lattice and the first Brillouin zone (b).

<sup>2</sup> ) cos(*<sup>k</sup>* <sup>⋅</sup>

trum of bilayer graphene can be written as

146 New Progress on Graphene Research

where *ε*<sup>0</sup>

imation:

where *a*<sup>1</sup> =( *<sup>a</sup>*

Figure 5b).

*ε*0

<sup>2</sup> , *<sup>a</sup>* <sup>3</sup>

*ε*1,2(*k*)= ± *ε*<sup>0</sup>

*ε*3,4(*k*)= ± *ε*<sup>0</sup>

(*k*)= ± *J* 1 + 4cos(*k* ⋅

<sup>2</sup> ,0) and *a*<sup>2</sup> =( *<sup>a</sup>*

Å is the distance between nearest neighbors in the layer plane). Since this distance is only by less than 10% greater than *h* , we can neglect the interaction with the atoms of the sublattices *BI* and *BII* and this does not lead to qualitative changes in the behavior of the spectra near *ε<sup>F</sup>* (see, e.g. [6]). Then the dispersion relation of each of the four branches of the electronic spec‐

<sup>2</sup> <sup>−</sup> *<sup>J</sup>* ′ *<sup>ε</sup>*<sup>0</sup>

<sup>2</sup> <sup>+</sup> *<sup>J</sup>* ′ *<sup>ε</sup>*<sup>0</sup>

(*k*) is the electronic spectrum of graphene, calculated in the strong coupling approx‐

*a*<sup>1</sup> + *a*<sup>2</sup>

<sup>2</sup> ) <sup>+</sup> cos(*<sup>k</sup>* <sup>⋅</sup>

*a*<sup>1</sup> - *a*<sup>2</sup>

<sup>2</sup> ,0) are two-dimensional Bravais lattice vectors (see

<sup>2</sup> ) (18)

2(*k*) <sup>+</sup> *<sup>J</sup>* ′<sup>2</sup> 4 ;

2(*k*) <sup>+</sup> *<sup>J</sup>* ′<sup>2</sup> 4 , ( *a* ≈1.415

(17)

$$
\varepsilon\_{1,2}(\text{K} + \kappa) \approx \pm \frac{9f^2 a^2 \kappa^2}{4f\_1} \tag{20}
$$

**Figure 6.** Dispersion curves along high-symmetry directions of bilayer graphene.

The electron effective mass in considered branches, determined from the relation *ε* =(ℏ*κ*)<sup>2</sup> / 2*m*\* , is equal to

$$
\mu m^\* = \frac{2\hbar^2 f}{9f^2 a^2}^{'}\tag{21}
$$

Total electron DOS can be represented as a mean-arithmetic function of the two LDOS corre‐

*g*(*ε*)= *ρA*(*ε*) + *ρB*(*ε*) / 2. For each sublattice with perfect structure the LDOS may be written as

where index *s* is the designation of sublattice, index α is the designation of branch, and

In this figure we clearly see two-dimensional van-Hove peculiarities at energy values corre‐ sponding to points *Γ* and *Μ* of the first Brillouin zone ( see Figure 5b). These eigenvalues are given in the top inset in Figure 7. The bottom inset shows LDOS, in enlarged scale, near *εF* . It is seen that near the Fermi level the local density of states as well the total density of states are analytical and their dependences on energy are essentially nonlinear (for comparison the DOS of graphene is also shown in the same inset). Besides, *ρA*(*ε*) differs from both *ρB*(*ε*) and

(*α*, *κ*) are the eigenwave functions corresponding to atoms from sublattices. LDOS are cal‐

*δdφ* |*ψ<sup>s</sup>*

(*ε*)=*ρAII*

Electronic and Vibrational Properties of Adsorbed and Embedded Graphene and Bigraphene with Defects

(*α*, *κ*)| <sup>2</sup>

(*ε*)≡*ρA*(*ε*); *ρBI*

(*ε*)=*ρBII*

http://dx.doi.org/10.5772/50562

<sup>|</sup>∂*ε*1,2 / <sup>∂</sup>*<sup>κ</sup>* <sup>|</sup> (24)

(*ε*)≡*ρB*(*ε*) and

149

sponding to atoms of sublattices *A* and *B*: *ρAI*

*ψs*

*<sup>ρ</sup>s*(*ε<sup>F</sup>* ) <sup>≈</sup> *<sup>S</sup>*<sup>0</sup>

the total density of states and it approaches to zero for *ε* →0 .

**Figure 7.** LDOS for atoms of different sublattices of bigraphene.

(2*π*)<sup>2</sup> <sup>∑</sup> *α*=1 4 *∮ ε*(*δ*)=*ε*

culated by the method of Jacobi matrix [16, 17] and are shown in Figure 7.

In the case considered above *J* ′ =0.1*J* the effective mass is *m*\* ≈2.75⋅10−<sup>32</sup> *kg*, and if *J* ′ → *J* it tends to a value close to the free electron mass. Since *m*\* ~ *J* ′ , then by changing the interlayer hopping integral we can change the effective mass of charge carriers.

It should be noted that for *J* ′ →0 , that is for the transition from bilayer graphene to two noninteracting graphene monolayers the effective mass, *m*\* →0 , and formulas (19 - 21) cannot cor‐ rectly describe this transition, since they were obtained under the assumption *ε*<sup>0</sup> (*k*)< < *J* ′ .

Electron density of states for values of energy near ε<sup>F</sup> is determined by branches of *ε*1 and *ε*<sup>2</sup> only, and it follows from (3) that *g*<sup>1</sup> (*ε*)= *g*<sup>2</sup> (−*ε*) . Then

$$\chi\_{\mathcal{S}}(\varepsilon) = \frac{\Sigma\_0}{(2\pi)^2} \oint\_{\varepsilon(\kappa) = \varepsilon} \frac{d l\_{1,2}}{\lceil \partial \varepsilon\_{1,2} / \partial \kappa \rceil} \tag{22}$$

where *Σ*<sup>0</sup> =3*<sup>a</sup>* <sup>2</sup> <sup>3</sup> / 2 is the area of the two-dimensional Bravais cell, and integration is done along a closed isoenergetic line *ε*1,2(*k*)=*ε* . At *ε* =0 (Fermi level) the line contracts into a point and near εF the contour of integration is a circle. Taking into account (20) we may write

$$g\_{1,2}(\varepsilon\_F) \approx \frac{S\_0}{(2\pi)^2} \oint\_{\varepsilon(\delta) = \epsilon} \frac{\delta d\varrho}{|\operatorname{\bf }\varepsilon\_{1,2}/\partial\kappa \,|} = \frac{3a^2 \sqrt{3}}{2(2\pi)^2} \cdot 2\pi \frac{\delta}{\underbrace{18f^2 a^2 \kappa}\_{4f}} = \frac{f}{6\pi f^2} = const.\tag{23}$$
 
$$g\_{3,4}(\varepsilon\_F) = 0; \Rightarrow g(\varepsilon\_F) \approx \frac{f}{f^2 \sqrt{3}}.\tag{24}$$

This means that the electron density of states (DOS) is constant and different from zero. As follows from (20), DOS is an analytical function and has minimum at *ε* =0 and near *ε<sup>F</sup>* the function *g*(*ε*)~*ε* <sup>2</sup> .

Total electron DOS can be represented as a mean-arithmetic function of the two LDOS corre‐ sponding to atoms of sublattices *A* and *B*: *ρAI* (*ε*)=*ρAII* (*ε*)≡*ρA*(*ε*); *ρBI* (*ε*)=*ρBII* (*ε*)≡*ρB*(*ε*) and *g*(*ε*)= *ρA*(*ε*) + *ρB*(*ε*) / 2. For each sublattice with perfect structure the LDOS may be written as

The electron effective mass in considered branches, determined from the relation

=0.1*J* the effective mass is *m*\*

Electron density of states for values of energy near ε<sup>F</sup> is determined by branches of *ε*1 and *ε*<sup>2</sup>

where *Σ*<sup>0</sup> =3*<sup>a</sup>* <sup>2</sup> <sup>3</sup> / 2 is the area of the two-dimensional Bravais cell, and integration is done along a closed isoenergetic line *ε*1,2(*k*)=*ε* . At *ε* =0 (Fermi level) the line contracts into a point and near εF the contour of integration is a circle. Taking into account (20) we may write

2(2*π*)<sup>2</sup> <sup>⋅</sup>2*<sup>π</sup>*

This means that the electron density of states (DOS) is constant and different from zero. As follows from (20), DOS is an analytical function and has minimum at *ε* =0 and near *ε<sup>F</sup>* the

*J* <sup>2</sup> 3 .

*<sup>g</sup>*3,4(*ε<sup>F</sup>* ) =0; <sup>⇒</sup>*g*(*ε<sup>F</sup>* ) <sup>≈</sup> *<sup>J</sup>* ′

*δ* 18*J* <sup>2</sup> *a* 2 *κ* 4*J* ′

(−*ε*) . Then

*dl* 1,2 ~ *J* ′

→0 , that is for the transition from bilayer graphene to two non-

*<sup>a</sup>* <sup>2</sup> (21)

≈2.75⋅10−<sup>32</sup> *kg*, and if *J* ′

→0 , and formulas (19 - 21) cannot cor‐

<sup>|</sup>∂*ε*1,2 / <sup>∂</sup>*<sup>κ</sup>* <sup>|</sup> (22)

<sup>=</sup> *<sup>J</sup>* ′ <sup>3</sup>

<sup>6</sup>*π<sup>J</sup>* <sup>2</sup> <sup>=</sup>*const*

(*k*)< < *J* ′ .

, then by changing the interlayer

→ *J* it

(23)

*m*\* <sup>=</sup> <sup>2</sup>ℏ<sup>2</sup> *J* ′ 9*J* <sup>2</sup>

rectly describe this transition, since they were obtained under the assumption *ε*<sup>0</sup>

(*ε*)= *g*<sup>2</sup>

(2*π*)<sup>2</sup> *<sup>∮</sup> ε*(*κ*)=*ε*

*<sup>g</sup>*(*ε*)= *<sup>Σ</sup>*<sup>0</sup>

*δdφ* <sup>|</sup>∂*ε*1,2 / <sup>∂</sup>*<sup>κ</sup>* <sup>|</sup> <sup>=</sup> <sup>3</sup>*<sup>a</sup>* <sup>2</sup> <sup>3</sup>

*ε* =(ℏ*κ*)<sup>2</sup> / 2*m*\*

148 New Progress on Graphene Research

, is equal to

In the case considered above *J* ′

It should be noted that for *J* ′

only, and it follows from (3) that *g*<sup>1</sup>

*<sup>g</sup>*1,2(*ε<sup>F</sup>* ) <sup>≈</sup> *<sup>S</sup>*<sup>0</sup>

.

function *g*(*ε*)~*ε* <sup>2</sup>

(2*π*)<sup>2</sup> *<sup>∮</sup> ε*(*δ*)=*ε*

tends to a value close to the free electron mass. Since *m*\*

interacting graphene monolayers the effective mass, *m*\*

hopping integral we can change the effective mass of charge carriers.

$$\rho\_s(\varepsilon\_F) \approx \frac{S\_0}{(2\pi)^2} \sum\_{s=1}^4 \oint\_{\varepsilon(\delta) = \varepsilon} \frac{\delta d\rho \parallel \psi\_s(\alpha, \kappa) \parallel \,^2}{\vert \, \partial \varepsilon\_{1,2} / \partial \kappa \vert} \,\tag{24}$$

where index *s* is the designation of sublattice, index α is the designation of branch, and *ψs* (*α*, *κ*) are the eigenwave functions corresponding to atoms from sublattices. LDOS are cal‐ culated by the method of Jacobi matrix [16, 17] and are shown in Figure 7.

In this figure we clearly see two-dimensional van-Hove peculiarities at energy values corre‐ sponding to points *Γ* and *Μ* of the first Brillouin zone ( see Figure 5b). These eigenvalues are given in the top inset in Figure 7. The bottom inset shows LDOS, in enlarged scale, near *εF* . It is seen that near the Fermi level the local density of states as well the total density of states are analytical and their dependences on energy are essentially nonlinear (for comparison the DOS of graphene is also shown in the same inset). Besides, *ρA*(*ε*) differs from both *ρB*(*ε*) and the total density of states and it approaches to zero for *ε* →0 .

**Figure 7.** LDOS for atoms of different sublattices of bigraphene.

Indeed, putting zero eigenvalue in the equation of eigenfunctions of Hamiltonian (1) we get the values *ψAI* <sup>~</sup>*ψAII* <sup>=</sup>*O*(*<sup>κ</sup>* 2) , *ψBI* <sup>~</sup>*ψBII* ~1 . Therefore, near the Fermi level *ρA*(*ε*)~*<sup>ε</sup>* <sup>2</sup> *ρB*(*ε*) .

This figure also shows the dependences S(ε) appearing in the Lifshitz equation (2) and corre‐ sponding correspond to the presence of substitutional impurities of boron and nitrogen in

Electronic and Vibrational Properties of Adsorbed and Embedded Graphene and Bigraphene with Defects

For an isolated substitutional impurity that differs from the atom of the basic lattice by the values of energies at the impurity site and also by the overlap integral, function *S*(*ε*) has for each sublattice form (3). Also, as in the case of graphene, for the nitrogen impurity in the

states in this interval. For the boron impurity, equation (2) has two solutions outside the

LDOS of isolated impurity atom of nitrogen in sublattice *A* or *B* are shown in Figure 9 (we

As the inset shows, the nitrogen impurity does not forms a quasilocal maximum on LDOS and substantially changes it near the Fermi level. Figure 10 shows a LDOS of the boron im‐ purity in sublattices *A* and *B*. As in the case of boron impurity in graphene (see Figure 1), the area under this curve is less than unity. Outside the band of the quasi-continuous spec‐

**Figure 9.** LDOS of nitrogen impurity in the sublattices of bigraphene (red curve is for *A* sublattice, purple is for *B* sub‐ lattice). For comparison the figure shows the DOS of graphene (black dashed line) and LDOS (thin solid gray and dash‐

(*Μ*), *ε*<sup>3</sup>

remind that atoms *A*I and *A*II as well as *B*I and *B*II are physically equivalent).

(*Μ*) equation (2) has a solution and this impurity forms quasilocalized

(*Μ*) , corresponding to local discrete levels.

http://dx.doi.org/10.5772/50562

151

bigraphene.

interval *ε*<sup>2</sup>

ed gray, respectively).

(*Μ*), *ε*<sup>1</sup>

band of quasi-continuous spectrum *ε*<sup>4</sup>

trum local discrete energy levels are formed.

Peculiarities at both ends, i.e. at *ε* ≈ ± 0.1*J* = ± *J* ′ , originate from the contribution of modes ε3,4 which are not represented on interval *ε* ∈ − *J* ′ , *J* ′ . Peculiarities at *ε* ≈ ± 0.05*J* are due to the fact that beginning from these energies the anisotropy of isofrequency lines becomes essential.

So we can conclude that, in contrary to graphene, the bigraphene has an ordinary non-rela‐ tivistic form of the electronic dispersion law. The effective mass of electron in the bigra‐ phene strongly depends on the value of integral describing the hopping between two layers, and this value may be changed by external conditions (for example by pressure). Near the Fermi level the LDOS of atoms of different sublattices qualitively differ from each other. If the LDOS for the atoms of sublattice *A* at the Fermi level equals to zero and it slowly in‐ creases near this level ("a quasi-gap" appears), then the LDOS of sublattice *B* for the same energy values differs from zero and increases very quickly.

## **4. The influence of defects on the electron spectrum of bigraphene**

Some peculiarities in the behavior of the bigraphene electron spectrum near *ε* =*ε*(*Κ*) , that is near the Fermi level, indicate the possibility of a strong influence of various defects [26-28].

Let us first consider point defects [26]. Figure 8 shows the real parts of the local Green's functions Re*GA*(*ε*) and Re*GB*(*ε*) of atoms of the sublattices *A* and *B*.

**Figure 8.** Real parts of the local Green's functions for the atoms of the two sublattices of bigraphene. The dashed lines are *S*(ε) functions correspond to substitutional impurities of Boron (B-curves) and Nitrogen (N-curves).

This figure also shows the dependences S(ε) appearing in the Lifshitz equation (2) and corre‐ sponding correspond to the presence of substitutional impurities of boron and nitrogen in bigraphene.

Indeed, putting zero eigenvalue in the equation of eigenfunctions of Hamiltonian (1) we get

fact that beginning from these energies the anisotropy of isofrequency lines becomes essential.

So we can conclude that, in contrary to graphene, the bigraphene has an ordinary non-rela‐ tivistic form of the electronic dispersion law. The effective mass of electron in the bigra‐ phene strongly depends on the value of integral describing the hopping between two layers, and this value may be changed by external conditions (for example by pressure). Near the Fermi level the LDOS of atoms of different sublattices qualitively differ from each other. If the LDOS for the atoms of sublattice *A* at the Fermi level equals to zero and it slowly in‐ creases near this level ("a quasi-gap" appears), then the LDOS of sublattice *B* for the same

*ρB*(*ε*) .

, originate from the contribution of modes ε3,4

, *J* ′ . Peculiarities at *ε* ≈ ± 0.05*J* are due to the

the values *ψAI* <sup>~</sup>*ψAII* <sup>=</sup>*O*(*<sup>κ</sup>* 2) , *ψBI* <sup>~</sup>*ψBII* ~1 . Therefore, near the Fermi level *ρA*(*ε*)~*<sup>ε</sup>* <sup>2</sup>

**4. The influence of defects on the electron spectrum of bigraphene**

Some peculiarities in the behavior of the bigraphene electron spectrum near *ε* =*ε*(*Κ*) , that is near the Fermi level, indicate the possibility of a strong influence of various defects [26-28].

Let us first consider point defects [26]. Figure 8 shows the real parts of the local Green's

**Figure 8.** Real parts of the local Green's functions for the atoms of the two sublattices of bigraphene. The dashed lines

are *S*(ε) functions correspond to substitutional impurities of Boron (B-curves) and Nitrogen (N-curves).

Peculiarities at both ends, i.e. at *ε* ≈ ± 0.1*J* = ± *J* ′

150 New Progress on Graphene Research

which are not represented on interval *ε* ∈ − *J* ′

energy values differs from zero and increases very quickly.

functions Re*GA*(*ε*) and Re*GB*(*ε*) of atoms of the sublattices *A* and *B*.

For an isolated substitutional impurity that differs from the atom of the basic lattice by the values of energies at the impurity site and also by the overlap integral, function *S*(*ε*) has for each sublattice form (3). Also, as in the case of graphene, for the nitrogen impurity in the interval *ε*<sup>2</sup> (*Μ*), *ε*<sup>1</sup> (*Μ*) equation (2) has a solution and this impurity forms quasilocalized states in this interval. For the boron impurity, equation (2) has two solutions outside the band of quasi-continuous spectrum *ε*<sup>4</sup> (*Μ*), *ε*<sup>3</sup> (*Μ*) , corresponding to local discrete levels.

LDOS of isolated impurity atom of nitrogen in sublattice *A* or *B* are shown in Figure 9 (we remind that atoms *A*I and *A*II as well as *B*I and *B*II are physically equivalent).

As the inset shows, the nitrogen impurity does not forms a quasilocal maximum on LDOS and substantially changes it near the Fermi level. Figure 10 shows a LDOS of the boron im‐ purity in sublattices *A* and *B*. As in the case of boron impurity in graphene (see Figure 1), the area under this curve is less than unity. Outside the band of the quasi-continuous spec‐ trum local discrete energy levels are formed.

**Figure 9.** LDOS of nitrogen impurity in the sublattices of bigraphene (red curve is for *A* sublattice, purple is for *B* sub‐ lattice). For comparison the figure shows the DOS of graphene (black dashed line) and LDOS (thin solid gray and dash‐ ed gray, respectively).

**Figure 10.** LDOS of the boron impurity in different sublattices of bigraphene (designations are the same as in Figure 9).

The energies of these levels for the considered ratio between the overlap integrals J and J' are slightly different. As can be seen from Figure 8, the energies of the local discrete levels for the substituted boron atom in sublattice *A* or sublattice *B* can be calculated using the two-moment approximation (13). Near the Fermi level the LDOS of the boron atom, substi‐ tuting an atom of sublattice *B*, is considerably lower than the LDOS of the carbon atom of this sublattice. That is the boron impurity lowers the conductivity of bigraphene.

Because in bigraphene atoms of sublattices *A* and *B* are physically inequivalent, the influ‐ ence on their electron spectra by various defects is different. In the first part of this section we have described the influence of substitutional impurities on the electron density of states. However, the influence of vacancies in sublattices on the electron DOS is even more pro‐ found. Figures 11 and 12 show the LDOS of neighbors of vacancies in sublattices *A* (Figure 11) and *B* (Figure 12).

Neighboring atoms are in the same layer as a vacancy, either in sublattice *B* (top) or in sub‐ lattice *A* (bottom). Insets show the arrangement of atoms. Atoms are shown in the same col‐ or as the corresponding LDOS.

**Figure 11.** LDOS of bigraphene atoms which are neighbors of a vacancy in the sublattice *A*. Neighboring atoms are in the same layer as a vacancy, either in sublattice *B* (top) or in sublattice *A* (bottom). Insets show the arrangement of

Electronic and Vibrational Properties of Adsorbed and Embedded Graphene and Bigraphene with Defects

http://dx.doi.org/10.5772/50562

153

atoms. Atoms are shown in the same color as the corresponding LDOS.

Electronic and Vibrational Properties of Adsorbed and Embedded Graphene and Bigraphene with Defects http://dx.doi.org/10.5772/50562 153

**Figure 10.** LDOS of the boron impurity in different sublattices of bigraphene (designations are the same as in Figure 9).

The energies of these levels for the considered ratio between the overlap integrals J and J' are slightly different. As can be seen from Figure 8, the energies of the local discrete levels for the substituted boron atom in sublattice *A* or sublattice *B* can be calculated using the two-moment approximation (13). Near the Fermi level the LDOS of the boron atom, substi‐ tuting an atom of sublattice *B*, is considerably lower than the LDOS of the carbon atom of

Because in bigraphene atoms of sublattices *A* and *B* are physically inequivalent, the influ‐ ence on their electron spectra by various defects is different. In the first part of this section we have described the influence of substitutional impurities on the electron density of states. However, the influence of vacancies in sublattices on the electron DOS is even more pro‐ found. Figures 11 and 12 show the LDOS of neighbors of vacancies in sublattices *A* (Figure

Neighboring atoms are in the same layer as a vacancy, either in sublattice *B* (top) or in sub‐ lattice *A* (bottom). Insets show the arrangement of atoms. Atoms are shown in the same col‐

this sublattice. That is the boron impurity lowers the conductivity of bigraphene.

11) and *B* (Figure 12).

152 New Progress on Graphene Research

or as the corresponding LDOS.

**Figure 11.** LDOS of bigraphene atoms which are neighbors of a vacancy in the sublattice *A*. Neighboring atoms are in the same layer as a vacancy, either in sublattice *B* (top) or in sublattice *A* (bottom). Insets show the arrangement of atoms. Atoms are shown in the same color as the corresponding LDOS.

The reason for the specific evolution of the LDOS with increasing distance from the vacancy has been explained in [22]. In bigraphene there is an analogous situation. Electronic spectra of vacancy neighbors belonging to other sublattice also have sharp resonance peaks at ε = ε<sup>F</sup> (see upper parts Figure 11 and 12). So, if the vacancy is in the sublattice *A*, then the maxi‐ mum of the LDOS of the sublattice *B* atoms is sharp, and the maximum on the spectrum the sublattice *A* atoms is blurred. Therefore we can conclude that the vacancy in bigraphene

Electronic and Vibrational Properties of Adsorbed and Embedded Graphene and Bigraphene with Defects

http://dx.doi.org/10.5772/50562

155

Unique properties of both graphene and bigraphene are caused, above all, by an unusual symmetry leading to the absence of a gap between the valence and conduction bands. The quantum states of quasiparticle in these bands are described by the same wave function, i.e. quasiparticle in graphene and bigraphene have the so-called chiral symmetry. In graphene this leads to the fact that dispersion relation of electron spectrum is linear and is described by the Dirac equation, characteristic for massless ultrarelativistic quasiparticle in quantum electrodynamics. On other hand, in bigraphene the presence of chiral symmetry leads at low energies to an ordinary parabolic dispersion relation, i.e. quiasiparticle of a new type appear – massive chiral fermions having no analogy in quantum electrodynamics. At symmetry breaking a gap between the valence and conduction band appears, allowing to tune the con‐

In this chapter, using the method of Jacobi matrices, we analyzed how different impurities affect the energy gap width and the local density of states. The method of Jacobi matrices enables to investigate heterogeneous systems and to calculate the densities of states for each atom in different sublattices. Such analysis is necessary for correct comparison with experi‐ mental results. As defects we considered the vacancies and the impurity of nitrogen or bor‐ on. The presence of nitrogen leads to the formation of sharp resonance peaks (quasilocal states) inside the continuous spectrum; on the other hand, the boron impurity leads to states outside the continuous spectrum (local states). Both quasilocal and local states can be inves‐ tigated experimentally by, for example, a scanning tunneling microscope. In the presence of vacancies we have analyzed how the density of states in each of the sublattices A or B changes. Different situations were analyzed. For example, if the vacancy in this sublattice vanishes at Fermi level, whereas in sublattice B the density has maximum at this point. We also investigated the conditions for opening the gap and changing its width. Main attention was paid to the analysis of electronic properties of considered systems. Moreover, computa‐ tional method we used has also been successfully applied to the analysis of vibrational states. For example, increasing the overlap integral between the boron impurity and carbon atom leads to the strengthening of force interaction contacts between them. In addition, bor‐ on atom is 16% lighter than carbon atom, i.e. all necessary conditions are sent for the appear‐ ance of local vibrations in the phonon spectrum of both graphene and bigraphene with boron impurity. We hope that predicted peculiarities in electronic and vibrational spectra of perfect graphene and bigraphene as well as of graphene and bigraphene with defects will be

should have a more pronounced effect.

ducting properties of these materials.

detected experimentally.

**5. Conclusion**

**Figure 12.** LDOS of bigrafene atoms which are neighbors of a vacancy in the sublattice *B*.

The reason for the specific evolution of the LDOS with increasing distance from the vacancy has been explained in [22]. In bigraphene there is an analogous situation. Electronic spectra of vacancy neighbors belonging to other sublattice also have sharp resonance peaks at ε = ε<sup>F</sup> (see upper parts Figure 11 and 12). So, if the vacancy is in the sublattice *A*, then the maxi‐ mum of the LDOS of the sublattice *B* atoms is sharp, and the maximum on the spectrum the sublattice *A* atoms is blurred. Therefore we can conclude that the vacancy in bigraphene should have a more pronounced effect.

## **5. Conclusion**

**Figure 12.** LDOS of bigrafene atoms which are neighbors of a vacancy in the sublattice *B*.

154 New Progress on Graphene Research

Unique properties of both graphene and bigraphene are caused, above all, by an unusual symmetry leading to the absence of a gap between the valence and conduction bands. The quantum states of quasiparticle in these bands are described by the same wave function, i.e. quasiparticle in graphene and bigraphene have the so-called chiral symmetry. In graphene this leads to the fact that dispersion relation of electron spectrum is linear and is described by the Dirac equation, characteristic for massless ultrarelativistic quasiparticle in quantum electrodynamics. On other hand, in bigraphene the presence of chiral symmetry leads at low energies to an ordinary parabolic dispersion relation, i.e. quiasiparticle of a new type appear – massive chiral fermions having no analogy in quantum electrodynamics. At symmetry breaking a gap between the valence and conduction band appears, allowing to tune the con‐ ducting properties of these materials.

In this chapter, using the method of Jacobi matrices, we analyzed how different impurities affect the energy gap width and the local density of states. The method of Jacobi matrices enables to investigate heterogeneous systems and to calculate the densities of states for each atom in different sublattices. Such analysis is necessary for correct comparison with experi‐ mental results. As defects we considered the vacancies and the impurity of nitrogen or bor‐ on. The presence of nitrogen leads to the formation of sharp resonance peaks (quasilocal states) inside the continuous spectrum; on the other hand, the boron impurity leads to states outside the continuous spectrum (local states). Both quasilocal and local states can be inves‐ tigated experimentally by, for example, a scanning tunneling microscope. In the presence of vacancies we have analyzed how the density of states in each of the sublattices A or B changes. Different situations were analyzed. For example, if the vacancy in this sublattice vanishes at Fermi level, whereas in sublattice B the density has maximum at this point. We also investigated the conditions for opening the gap and changing its width. Main attention was paid to the analysis of electronic properties of considered systems. Moreover, computa‐ tional method we used has also been successfully applied to the analysis of vibrational states. For example, increasing the overlap integral between the boron impurity and carbon atom leads to the strengthening of force interaction contacts between them. In addition, bor‐ on atom is 16% lighter than carbon atom, i.e. all necessary conditions are sent for the appear‐ ance of local vibrations in the phonon spectrum of both graphene and bigraphene with boron impurity. We hope that predicted peculiarities in electronic and vibrational spectra of perfect graphene and bigraphene as well as of graphene and bigraphene with defects will be detected experimentally.

## **Acknowledgements**

This work is the result of the project implementation: Research and Education at UPJŠ – Heading towards Excellent European Universities, ITMS project code: 26110230056, sup‐ ported by the Operational Program Education funded by the European Social Fund (ESF). This work was also supported by the grant of the Ukrainian Academy of Sciences under the contract No 4/10-H and by the grant of the Scientific Grant Agency of the Ministry of Educa‐ tion, Science, Research and Sport of the Slovak Republic and Slovak Academy of Sciences under No. 1/0159/09.

[7] Meyer, J. C., Geim, A. K., Katsnelson, M. I., Novoselov, K. S., Booth, T. J., & Roth, S.

Electronic and Vibrational Properties of Adsorbed and Embedded Graphene and Bigraphene with Defects

http://dx.doi.org/10.5772/50562

157

[11] Slonchechevski, J. C., & Weiss, P. R. (1958). Band Structure of Graphite. *Physical Re‐*

[12] Feher, A., Syrkin, E. S., Feodos'ev, S. B., Gospodarev, I. A., & Kravchenko, K. V. (2011). Quasi-Particle Spectra on Substrate and Embedded Graphene Monolayers. In: Mikhailov S. (ed), *Physics and Applications of Graphene Theory*, Rijeka, InTech, 93-112.

[13] Wu, S., Jing, L., Li, Q., Shi, Q. W., Chen, J., Su, H., Wang, X., & Yang, J. (2008). Aver‐ age Density of States in Disordered Graphene Systems. *Physical Review B*, 77(19),

[14] Feher, A., Gospodarev, I. A., Grishaev, V. I., Kravchenko, K. V., Manzhelii, E. V., Syr‐ kin, E. S., & Feodos'ev, S. B. (2009). Effect of Defects on Quasi-particle Spectra of

[15] Chen, L., Zhang, Y., & Shen, Y. (2007). Self-healing in Defective Carbon Nanotubes: a

[16] Peresada, V. I. (1968). New Calculation Method in the Theory of Crystal Lattice. *Con‐ densed Matter Physics*, 4, ed B.I. Verkin (Kharkov: FTINT AN Ukr. SSR), 172, (in Russi‐

[17] Peresada, V. I., Afanas'ev, V. N., & Borovikov, V. S. (1975). On Calculation of Density of States of Single-Magnon Perturbations in Ferromagnetics. *Soviet Low Temperature*

[18] Skrypnyk, Yu. V., & Loktev, V. M. (2006). Impurity Effects in a Two-Dimensional

[19] Skrypnyk, Yu. V., & Loktev, V. M. (2008). Spectral Function of Graphene with Short-

[20] Bena, C., & Kivelson, S. A. (2005). Qusiparticle Scattering and Local Density of States

[21] Lifshits, I. M. (1945). On the Theory of regular Perturbations Report of AS SSSR (in

[22] Feher, A., Syrkin, E. S., Feodosyev, S. B., Gospodarev, I. A., Manzhelii, E. V., Kotlyar, A. V., & Kravchenko, K. V. (2011). The Features of Low Frequency Atomic Vibrations and Properties of Acoustic Waves in Heterogeneous Systems. In: Vila R. (ed.), *Waves*

(2007). The Structure of Suspended Graphene Sheets. *Nature*, 446, 60-63.

in Laminar Crystals. *Fizika Nizkih Temperatur*, 5(9), 1069-1073.

Graphite and Graphene. *Low Temperature Physics*, 35(7), 679-686.

Molecular Dynamic Study. *Journal of Physics: Cond. Matter*, 9, 38612-6.

Systems with the Dirac Spectrum. *Physical Review B*, 73, 241402-6.

Range Imourity Centers. *Low Temperature Physics*, 34(9), 818-825.

in Graphite. *Physical Review B*, 72(12), 125432-7.

*in Fluids and Solids*, Rijeka, InTech, 103-126.

*view*, 109(2), 272-279.

195411.

an).

*Physics*, 1(4), 227-232.

Russian). 48, 83-86.

[8] Wallace, P. (1947). The Band Theory of Graphite. *Physical Review*, 71(9), 622-634. [9] McClure, J. W. (1956). Diamagnetizm of Graphite. *Physical Review*, 104(3), 666-671. [10] Syrkin, E. S., & Feodosyev, S. B. (1979). The Phonon Spectrum and Local Vibrations

## **Author details**

Alexander Feher1\*, Eugen Syrkin2 , Sergey Feodosyev2 , Igor Gospodarev2 , Elena Manzhelii2 , Alexander Kotlar2 and Kirill Kravchenko2

\*Address all correspondence to: alexander.feher@upjs.sk

1 Institute of Physics, Faculty of Science, P. J. Šafárik University in Kosice, Slovakia

2 B.I.Verkin Institute for Low Temperature Physics and Engineering NASU, Ukraine

## **References**


**Acknowledgements**

156 New Progress on Graphene Research

under No. 1/0159/09.

**Author details**

Alexander Kotlar2

**References**

Alexander Feher1\*, Eugen Syrkin2

and Kirill Kravchenko2

Graphene. *Physical Review Letters*, 98, 14681-4.

ley-VCH, 10.1002/3527603085.

*EPL*, 80, 67007-67017.

Doped Graphene. *Physical Review Letters*, 104, 13683-4.

Dirac Fermions in Graphene. *Nature*, 438, 197-201.

\*Address all correspondence to: alexander.feher@upjs.sk

This work is the result of the project implementation: Research and Education at UPJŠ – Heading towards Excellent European Universities, ITMS project code: 26110230056, sup‐ ported by the Operational Program Education funded by the European Social Fund (ESF). This work was also supported by the grant of the Ukrainian Academy of Sciences under the contract No 4/10-H and by the grant of the Scientific Grant Agency of the Ministry of Educa‐ tion, Science, Research and Sport of the Slovak Republic and Slovak Academy of Sciences

, Sergey Feodosyev2

1 Institute of Physics, Faculty of Science, P. J. Šafárik University in Kosice, Slovakia

2 B.I.Verkin Institute for Low Temperature Physics and Engineering NASU, Ukraine

[1] Uchoa, B., & Castro Neto, A. H. (2007). Superconducting States of Pure and Doped

[2] McChesney, J. L., Bostwick, A., Ohta, T., Seyller, T., Horn, K., Gonza´lez, J., & Roten‐ berg, E. (2010). Extended van Hove Singularity and Superconducting Instability in

[3] Kosevich, A. M. (1999). *The Crystal Lattice (Phonons, Solitons, Dislocations)*, Berlin, Wi‐

[4] Novoselov, K. S., Gein, A. K., Morozov, S. V., Diang, D., Katsnelson, M. I., Grigorie‐ va, I. V., Dubonos, S. V., & Firsov, A. A. (2005). Two-Dimensional Gas of Massless

[5] Castro Neto, A. H., Guinea, F., Peres, N. M. R., Novoselov, K. S., & Geim, A. K. (2009). The Electronic Properties of Graphene. *Review of Modern Physics*, 81, 109-155.

[6] Peres, N. M. R., Klironomos, F. D., Tsai, S-W., Santos, J. R., Lopes dos Santos, J. M. B., & Castro Neto, A. H. (2007). Electron Waves in Chemistry Substituted Graphene.

, Igor Gospodarev2

, Elena Manzhelii2

,


[23] Kotlyar, A. V., & Feodosyev, S. B. (2006). Local Vibrational Modes in Crystal Lattices with a Simply Connected Region of the Quasi-continuous Phonon Spectrum. *Low Temperature Physics*, 32(3), 256-269.

**Section 2**

**Experimental Aspect**

