Theoretical research of crooked graphene Professor Olga E. Glukhova - - PowerPoint PPT Presentation

theoretical research of crooked graphene
SMART_READER_LITE
LIVE PREVIEW

Theoretical research of crooked graphene Professor Olga E. Glukhova - - PowerPoint PPT Presentation

Theoretical research of crooked graphene Professor Olga E. Glukhova Saratov State University e-mail: glukhovaoe@info.sgu.ru The energy of a system of ion cores and valence electrons is written as E E E (1) tot bond rep Here the


slide-1
SLIDE 1

Theoretical research of crooked graphene

Professor Olga E. Glukhova Saratov State University

e-mail: glukhovaoe@info.sgu.ru

slide-2
SLIDE 2

Tight-binding method (метод сильной связи)

О.Е. Глухова, СГУ , glukhovaoe@info.sgu.ru

2

The energy of a system of ion cores and valence electrons is written as

rep bond tot

E E E  

(1) Here the term

bond

E

is the bond structure energy that is calculated as the sum

  • f energies of the single-particle occupied states. Those single-particle

energies become known by solving the following equation

  

n n |

|

n

   H

, (2) where H is the one-electron Hamiltonian,

n

 is the energy of the nth single-

particle state. The wave functions

n

|

can be approximated by linear combination

  

   

l l l |

|

n n

C

, (3) where 

n is an orthogonal basis set, l is the quantum number index and

 labels the ions.

slide-3
SLIDE 3

О.Е. Глухова, СГУ , glukhovaoe@info.sgu.ru

3

The matrix elements in equation (2) are calculated after fitting a suitable database obtained from the. The overlap matrix elements, which takes into account four types of interaction ssσ, spσ, ppσ и ppπ and the pair repulsive potential are calculated by the formulas:

                                          

4 4 1

2 3 2 1 3

exp ) (

p p p ij

p p p r p r p V r V

 

, (4) where i and j are orbital moments of a wave function, α presents the bond type ( or  ). The bond structure energy is determined by the formula

n n bond

E  2 (5) This expression is the sum of the energies of molecular orbitals obtained by diagonalization of the Hamiltonian. The parameter n is the number

  • f occupied orbitals, and εn is the energy of single-particle orbitals.
slide-4
SLIDE 4

The phenomenon energy (феноменологическая энергия)

О.Е. Глухова, СГУ , glukhovaoe@info.sgu.ru

4

Term

rep

E

in equation (9) is phenomenon energy that is a repulsive potential. It can be expressed as a sum of two-body potentials as

 

    ,

r V E

rep rep

, (6) where

rep

V is a pair potential between atoms at  and  . This two-body potential

describes an interaction between bonded and non-bonded atoms.

                                          

4 4 6

2 3 2 6 3 5

exp ) (

p p p rep

p p p r p r p p r V

, (7) where i and j are orbital moments of wave function,  presents the bond type ( or

 ). The values of the parameters



V , the atomic terms and pn for carbon compounds

are given in table 1.

slide-5
SLIDE 5

О.Е. Глухова, СГУ , glukhovaoe@info.sgu.ru

5

Table 1. Values of the parameters

s

 , эВ

p

 , эВ

ss

V  , эВ

sp

V

 , эВ pp

V

, эВ pp

V

, эВ

  • 10,932
  • 5,991
  • 4,344

3,969 5,457

  • 1,938

p1 p2, Å p3, Å p4 p5, эВ p6 2,796 2,32 1,54 22 10,92 4,455

Optimization of atomic structure is implemented by entire system energy minimization on atomic coordinates. The study of the compression process was implemented with the algorithm presented earlier. The parameters were fitted from the experimental data for fullerenes and carbon nanotubes. Transferability to other carbon compounds was tested by comparison with ab initio calculations and experiments.

slide-6
SLIDE 6

The TB parameters for carbon nanoclusters

О.Е. Глухова, СГУ , glukhovaoe@info.sgu.ru

6

Our transferable tight-binding potential can reproduce changes correctly in the electronic configuration as a function of the local bonding geometry around each carbon atom.

slide-7
SLIDE 7

The Hamiltonian

О.Е. Глухова, СГУ , glukhovaoe@info.sgu.ru

7

slide-8
SLIDE 8

The interaction of 𝞀-orbitals

О.Е. Глухова, СГУ , glukhovaoe@info.sgu.ru

8

All S- and P-orbitals are given in the real Cartesian coordinates system. To reproduce changes in the electronic configuration of the local bonding geometry around each atom correctly we defined P-

  • rbital as an axial vector. Each axial vector makes the

angle with direction Rij (α, β, θ ) and it may be written as the geometrical sum of the two vectors:

 

x xD x

P P P   

,

 

y yD y

P P P   

,

 

z zD z

P P P   

Here

xD

P 

,

yD

P 

,

zD

P 

are projections to an interatomic direction,

 x

P  , etc are projections to an orthogonal direction. So, to describe the interaction between Pz and Px we must write: bonding

  • bonding

  

   

z x zD xD z x

P P P P P P      

slide-9
SLIDE 9

О.Е. Глухова, СГУ , glukhovaoe@info.sgu.ru

9

The angle between projections

xD

P  and

zD

P  is

equal to zero, but the other angle between projections to an orthogonal direction is not zero and it is equal to γ. As a result

  • f

some mathematical transformations we can write the expressions for cosγ and the energy of the interaction between Pz and Px in the following way:      sin sin cos cos cos     ,

 

) ( ) ( cos cos ) (

ij ij ij PxPz

r V r V r V

PxPz PxPz

 

    

slide-10
SLIDE 10

О.Е. Глухова, СГУ , glukhovaoe@info.sgu.ru

10

As well known, the expression for the energy of the interaction between S and P-orbitals can be defined simply:

cos ) ( ) (  

ij ij SPz

r V r V

SPz

slide-11
SLIDE 11

The rehybridization

О.Е. Глухова, СГУ , glukhovaoe@info.sgu.ru

11

The presented scheme provides the consideration and calculation

  • f

the rehybridization between σ- and π-orbitals to reproduce the electronic configuration and the local bonding geometry around each atom. In figure 6 we can see that the atom in sp2 hybridization becomes one in sp2+Δ hybridization because of a curvature of the topological

  • network. The degree of rehybridization is

defined by the pyramidalization angle. This angle is calculated by the following formula:

2   

 

p

Our own program was used to research the nanoribbons with the help of the tight-binding

  • method. It provides the calculation of the total

energy of nanostructures, which consist of 50- 5000 atoms. We adapted our TB method to be able to run the algorithm on a parallel computing machine (a computer cluster).

slide-12
SLIDE 12

The electron spectra

О.Е. Глухова, СГУ , glukhovaoe@info.sgu.ru

12

slide-13
SLIDE 13
  • 1. Optimization of the initial atomic structure

(optimized Hooke-Jeeves method)

Block scheme of the modified Hooke-Jeeves method Block scheme of the subprogram of investigation by sample

О.Е. Глухова, СГУ , glukhovaoe@info.sgu.ru

13

slide-14
SLIDE 14

Efficiency of the parallel calculations

At optimization of the structure with transition from four block cluster on eight block cluster the efficiency of calculations increases in 1.4 times. At carrying out of calculations by means of the parallel manner in four flow the efficiency of calculations increases in average of 1.8 times in comparison with sequential method.

О.Е. Глухова, СГУ , glukhovaoe@info.sgu.ru

14

slide-15
SLIDE 15

Research of the local stress field of the atomic grid of graphene nanoribbons and prediction of the appearance of defects in compression process

О.Е. Глухова, СГУ , glukhovaoe@info.sgu.ru

15

The method of the calculation of the local stress of atomic network include following steps:

  • 1. the optimization of the atomic structure of the stable unstrained graphene sheet by the

minimization of the total energy by the coordinates of atoms;

  • 2. the calculation of distribution of the bulk energy density on atoms of the stable

unstrained graphene;

  • 3. the optimization of the atomic network of the deformed structure;
  • 4. the calculation of distribution of the bulk energy density on atoms;
  • 5. the calculation of the local stress field of atomic.
slide-16
SLIDE 16

О.Е. Глухова, СГУ , glukhovaoe@info.sgu.ru

16

The bulk energy density of the bamboo-like nanotube was calculated by formula:

 

i i j ij VdW i j j i k k j i l ijkl tors i j ij A ij ij R i

V r V V r V B r V w                            

    

     ) ( , , , ) (

) ( ) ( ) ( ) ( 

(8) where ) ( ij

r VR

and ) ( ij

r VA

  • the pair potentials of repulsion and attraction between chemically bonded atoms which are

determined by the atoms type and the distance between them.

ij

r is the distance between atoms i and j , i and j - are the

numbers of interaction atoms;

ij

B is the multiparticle term correcting interaction energy of the atoms pair i – j considering

specificity of interaction of σ-and π-electron clouds;

) (

ijkl tors

V 

is the torsional potential which is the function of the linear dihedral angle

ijkl

constructed on the basis of atoms with an edge on bond i – j (

l k, are the atoms forming chemical bonds

with i , j ). ) ( ij

r V

VdW

is the van der Waals interaction potential between the chemically unbounded atoms;

3

3 4 r

Vi 

is the

  • ccupancy volume of the atom i ,

r is the Van der Waals radius of the carbon atoms which is equal of 1.7 Å

The strain of the atomic framework near the atom with number i is calculated as:

i i i

w w   

where

i

w is the bulk energy density of atom i of the graphene sheet which is in equilibrium;

i

w is the bulk energy density of

the CBNT. The value of the bulk energy density

i

w in the centre of the graphene sheet is equal of -58.609 GPa. At the edges of the graphene sheet the bulk energy density is more because the atoms of the edges have only two links with other carbon atoms. It is equal -41.546 GPa. It is suggested that stress is equal to zero on the atoms in center and on edges.

slide-17
SLIDE 17

Approbation of method (prediction of the graphene destruction with defect of hydrogenation )

Q.X. Pei Y.W. Zhang, V.B. Shenoy, CARBON 48 ( 2010) 898–904

О.Е. Глухова, СГУ , glukhovaoe@info.sgu.ru

17

slide-18
SLIDE 18

Occurrence of the defect

О.Е. Глухова, СГУ , glukhovaoe@info.sgu.ru

18

slide-19
SLIDE 19

Prediction of the defects occurrence at strain

  • f the atomic grid

О.Е. Глухова, СГУ , glukhovaoe@info.sgu.ru

19

slide-20
SLIDE 20

The absorption of H-atom on the atomic network

О.Е. Глухова, СГУ , glukhovaoe@info.sgu.ru

20

slide-21
SLIDE 21

The total energy of the structure depends on the distance between the hydrogen atom and the carbon atom

О.Е. Глухова, СГУ , glukhovaoe@info.sgu.ru

21

90

i i

С Кривизна  

slide-22
SLIDE 22

Investigation of the atomic grid curvature influence on the adsorption properties

О.Е. Глухова, СГУ , glukhovaoe@info.sgu.ru

22

Кривизна 6,9% Кривизна 8,6% Graphene Enthalpy of reaction The planar

  • 23.31 kcal/mol

Curvature of 6,9 %

  • 28.13 kcal/mol

Curvature of 8,6 %

  • 31.59 kcal/mol
slide-23
SLIDE 23

Influence of the curvature on the potential barrier depth

О.Е. Глухова, СГУ , glukhovaoe@info.sgu.ru

23

Approximation formula: Dependence

  • f

the difference between potential energy minima for the planar and strained graphene nanoribbon from the surface curvature

) 16 . exp( 1 . C h  

slide-24
SLIDE 24

Compression of graphene nanoribbons

Graphene nanoribbons

Modeling of the compression process

О.Е. Глухова, СГУ , glukhovaoe@info.sgu.ru

24

slide-25
SLIDE 25

Comression of graphene nanoribbons

Nanoribbons and nanoparticles Compressed nanoribbon: 96% of initial length

О.Е. Глухова, СГУ , glukhovaoe@info.sgu.ru

25

Velocity of the compression 10 m/sec

slide-26
SLIDE 26

Change of electron structure at an axial compression

Compression of armchair ribbon (646 atoms)

DOS of π-electron system

О.Е. Глухова, СГУ , glukhovaoe@info.sgu.ru

26

98,1 % 90,1%

slide-27
SLIDE 27

Change of electron structure at an axial compression

Compression of zigzag ribbon (550 atoms)

DOS of π-electron system

О.Е. Глухова, СГУ , glukhovaoe@info.sgu.ru

27

98,0% 90,0%

slide-28
SLIDE 28

Change of the pyramidalization angle along the strained atomic framework Compression of ribbon

О.Е. Глухова, СГУ , glukhovaoe@info.sgu.ru

28

slide-29
SLIDE 29

Change of the half-waves quantity at compression (molecular dynamics, tight-binding method)

Compression of armchair ribbon: 98%

О.Е. Глухова, СГУ , glukhovaoe@info.sgu.ru

29

L=71Ǻ, λ/2=35.5Ǻ, A=2.2Ǻ 9 hexagons

L –length of nanoribbon , λ/2 – length of half-wave, A – amplitude of half-wave

L=181.7Ǻ, λ/2=60.5Ǻ, A=5.3Ǻ L=258.4Ǻ, λ/2=64.6Ǻ, A= 5.65Ǻ L=335.12Ǻ, λ/2=66.2Ǻ, A= 5.4Ǻ

14-15 hexagons

slide-30
SLIDE 30

Change of the half-waves quantity at compression (molecular dynamics, tight-binding method)

Compression of zigzag ribbon: 98%

О.Е. Глухова, СГУ , glukhovaoe@info.sgu.ru

30

L=65Ǻ, λ/2=32.5Ǻ, A=2.8Ǻ

12 hexagons

L=165.18Ǻ, λ/2=55.06Ǻ, A=5.4Ǻ L=198.7Ǻ, λ/2=49.6Ǻ, A=5.6Ǻ

L –length of nanoribbon , λ/2 – length of half-wave, A – amplitude of half-wave

20 hexagons

slide-31
SLIDE 31

Maps of stress

О.Е. Глухова, СГУ , glukhovaoe@info.sgu.ru

31

slide-32
SLIDE 32

Map of stress

О.Е. Глухова, СГУ , glukhovaoe@info.sgu.ru

32

slide-33
SLIDE 33

Nanoribbons armchair (552 atoms), zigzag (646 atoms)

Ionization potential Energy gap

О.Е. Глухова, СГУ , glukhovaoe@info.sgu.ru

33

Armchair, width 2,21 nm Zigzag, width 1,99 nm Zigzag, width 1,99 nm Armchair, width 2,21 nm

slide-34
SLIDE 34

Investigation of the strain and prediction of the defects formation and destruction of graphen nanoribbon

Tension of graphene nanoribbon

О.Е. Глухова, СГУ , glukhovaoe@info.sgu.ru

34

slide-35
SLIDE 35

Спасибо за внимание!

Саратовский государственный университет имени Н.Г. Чернышевского

О.Е. Глухова, СГУ , glukhovaoe@info.sgu.ru

35