THE CALCULATION OF THE ELECTRICAL CONDUCTIVITY OF GRAPHENE NANOFIBERS BY DISCRETE FOURIER TRANSFORMATION METHOD

УДК 538.9

DOI: 18.101/2306-2363-2019-2-3-35-38

THE CALCULATION OF THE ELECTRICAL CONDUCTIVITY OF GRAPHENE NANOFIBERS BY DISCRETE FOURIER TRANSFORMATION METHOD

Andrey Nomoev

Doctor of Physical and Mathematical Sciences

Head of the Laboratory of Physics of Composite Materials

Institute of Physical Materials Science SB RAS

Buryat State University

E-mail: nomoevav@mail.ru

Natalya Yumozhapova

Buryat State University

E-mail: ynat@mail.ru

Graphenenanoribbons, as a typical two-dimensional nanometer material, has shown its unique application potential in electrical characteristics, thermal properties, and thermoelectric properties by virtue of its novel electronic structure. Through studies of graphene&s electrical properties in recent years, progress has been made not only in experimental testing, but also in theoretical calculation. This aspect of graphene reviewed in this paper.

Для цитирования

Andrey Nomoev, Natalya Yumozhapova The calculation of the electrical conductivity graphene nanoribbons by discrete Fourier transformation method // Вестник Бурятского государственного университета. Химия. Физика. 2019. Вып. 2-3. С. 39-43.

Graphenenanoribbons (GNRs) are narrow strips of a one-atomic-layer thick two-dimensional material made of carbon atoms arranged in a honeycomb structure. Its fascinating electrical, optical, and mechanical properties ignited enormous interdisciplinary interest from the physics, chemistry, and materials science fields [1, 2]. These nanomaterials boast potential solutions to downsizing of future electronics in both size and energy consumption, as well as opening the door to faster processing.

We have calculated the electronic properties of graphene within the density functional theory framework. We used the SIESTA code (Spanish Initiative for Electronic Simulation with Thousand of Atoms) [3], which employs linear combinations of numerical atomic orbitals to solve the Kohn-Sham equations with three dimensional periodic boundary conditions within a supercell [4]. This is convenient for treating infinite systems such as crystals, but can also handle finite systems by making the supercell sufficiently large separating the objects. The framework in the supercell approach is Bloch&s theorem, which states that for a

ВЕСТНИК БУРЯТСКОГО ГОСУДАРСТВЕННОГО УНИВЕРСИТЕТА. ХИМИЯ. ФИЗИКА

periodic system the electronic wave function can be written as a product of a wavelike part and a cell-periodic part, i.e.

where n is a discrete band index and k is a reciprocal lattice vector belonging to the first Brillouin zone corresponding to the supercell. The theorem allows for mapping the Kohn-Sham eigenvalue problem into the reciprocal space, where one can separately obtain for each k-point a discrete set of eigenvalues and eigen states of the Hamiltonian. The expectation value of some one-body operator d is then calculated as

where for practical purposes the integral over the first Brillouin zone with volume Qbz is approximated by a sum over k-points with weight factors ©k (adding up to one). The numerical solution of the Kohn-Sham equation implies a finite basis set for the representation of the wave functions. The SIESTA method requires the use of atomic-like localized orbitals which guarantee the Hamiltonian and overlap matrix to be sparse. For each atom i positioned at Ri a set of atom-centered orbitals is defined as a product of a spherical harmonic and a radial function

where the distances are conveniently written in terms of ri = r — Ri and the angular momentum is labeled by l, m. The SIESTA basis implies that the calculation of the overlap matrix and most matrix elements of the Hamiltonian are two center integrals, which are effectively calculated within Fourier space where the convolution becomes a simple product. However, the density and some Hamiltonian matrix elements are calculated on a real-space grid, whose fineness is conveniently described by an energy grid cutoff. In SIESTA different parametrizations for the local density approximation and the generalized gradient approximation to the exchange-correlation functional are implemented. For our calculations we selected the generalized gradient approximation functional as described by the Perdew-Burke-Ernzerhof [5] which is proved to be one of the functionals that provide the smallest differences between calculated and experimental values in our field of interest. The current through the contact region has been calculated using Landauer-Buttiker formula [6]

Where Go = 2(e2/h) is the unit of quantum conductance and T(E, Vb) is the transmission probability of electrons incident at an energy E through the device under the potential bias Vb. The electrochemical potential difference between the left and right electrodes is eVb = — ^r.

The calculation of the electrical conductivity using the TranSiesta package consists of the following steps:

systems. This calculation is used as input to the self-consistency procedure in the Green function method.

Fig. 1. Represented the response current for an electrode temperature of 300 K.

fc -500 d

■1000 -1500 -2000

-2 -1.5 -1 -0.5 0 0.5 1 1.5

Bias, Volt

Fig. 2. Current-voltage characteristics of the armchair GNR.

The figure below represented the response current for an electrode temperature of 300 K. The low-bias current-voltage (I-V) characteristic of the device are linear (fig. 2). Because the main contribution to the conductivity of GNRs is made by electrons. In a defectless armchair GNRs electrons pass through the nanoribbon without scattering, which corresponds to the quantum ballistic mechanism of conductivity.

Basic science has much to gain from studies of graphene. The electrical, optical, mechanical, and thermal properties of graphene achieve extreme values that are hardly met by any other materials. In this paper, we studied the electrical properties of graphenenanoribbon. We computed the current-voltage curve by method that come from the first principles methods.

This work was carried out according to the State assignment of the Institute of Physical Materials Science SB RAS and according to the RFBR project No. 18-42030004.

ВЕСТНИК БУРЯТСКОГО ГОСУДАРСТВЕННОГО УНИВЕРСИТЕТА. ХИМИЯ. ФИЗИКА

References

— P. 183-191.

— 2012. V. 490. — P. 192-200.

РАСЧЕТ ЭЛЕКТРОПРОВОДНОСТИ ГРАФЕНОВЫХ НАНОВОЛОКОН МЕТОДОМ ДИСКРЕТНОГО ПРЕОБРАЗОВАНИЯ ФУРЬЕ

Андрей Номоев

доктор физико-математических наук

заведующий лабораторией физики композитных материалов Институт физического материаловедения СО РАН 670047, Улан-Удэ, ул. Сахьяновой, 6 Бурятский государственный университет 670000, Улан-Удэ, ул. Смолина, 24а Е-mail: nomoevav@mail.ru

Наталья Юможапова Бурятский государственный университет 670000, Улан-Удэ, ул. Смолина, 24а E-mail: ynat@mail.ru

Графенановолокна, как типичный двумерный нанометровый материал, показал свой уникальный потенциал применения в электрических характеристиках, тепловых свойствах и термоэлектрических свойствах благодаря своей новой электронной структуре. Благодаря исследованиям электрических свойств графена в последние годы был достигнут прогресс не только в экспериментальных испытаниях, но и в теоретических расчетах. Этот аспект графена рассмотрен в данной статье. Ключевые слова: графеновые нановолокна; метод теории функциональной плотности; вольт-амперные характеристики.