®
www.volsu.ru
DOI: https://doi.org/10.15688/mpcm.jvolsu.2017.4.9
УДК 541.138 ББК 24.5
ВЛИЯНИЕ АДСОРБЦИИ АТОМОВ И МОЛЕКУЛ КИСЛОРОДА НА ЭЛЕКТРОННОЕ СТРОЕНИЕ ГРАФЕНОВОЙ НАНОЛЕНТЫ
Елена Николаевна Шамина
Кандидат физико-математических наук, доцент кафедры математики и информатики, Волгоградский государственный медицинский университет shamina.alena@gmail.com
площадь Павших Борцов, 1, 400131 г. Волгоград, Российская Федерация
Николай Геннадьевич Лебедев
Доктор физико-математических наук,
профессор кафедры теоретической физики и волновых процессов, Волгоградский государственный университет шкоку. lebedev@volsu. ги, tf@volsu. ги
просп. Университетский, 100, 400062 г. Волгоград, Российская Федерация
Аннотация. В работе представлены результаты квантово-химических полуэмпирических исследований процессов адсорбции атомов и молекул кислорода на поверхности графеновой наноленты. Расчеты выполнены на основе модели молекулярного кластера с граничными псевдоатомами с использованием полуэмпирических расчетных схем MNDO и РМ3. Определены энергетические характеристики процессов адсорбции.
Ключевые слова: модель молекулярного кластера, квантовая химия, MNDO, РМ3, энергия адсорбции.
Введение
Графеновые наноленты (ГНЛ) рассматриваются в литературе в качестве перспективного материала для создания таких устройств, как полевые транзисторы, коммутационные устройства, сенсоры и др. Одним из преимуществ их использования является зависимость электронных свойств от характерных размеров наноленты и геометрии краев [7]. Электронные свойства ГНЛ ^ могут быть модифицированы посредством внесения в структуру нанолент дефектов, таких как о вакансий [6] или посредством адсорбции атомов или молекул [4; 5; 8].
г Ранее в ряде работ был изучен вопрос о зависимости электронных свойств оксидированных Д углеродных наноструктур от их геометрии и положения адатома. В частности, в работах [2-3] на Ц основе полуэмпирических исследований авторами была показана осциллирующая зависимость энер-$ гий адсорбции процессов оксидирования от диаметра и хиральности углеродной нанотрубки. Чередо-„ вание минимальных и максимальных значений энергий с ростом диаметра трубок позволило определить оптимальные диаметры хиральных УНТ, на которых адсорбция атомарного и молекулярного д кислорода реализуется наиболее эффективно в том смысле, что формируется энергетически более | устойчивая система. В работе [9] авторами, в рамках теории функционала плотности, был изучен д вопрос о зависимости электронных свойств молекулярно окидированного графена от положения ад© сорбированной молекулы кислорода относительно поверхности. Было показано, что конфигурация, в
которой связь O-O параллельна поверхности, характеризуется наибольшим значением энергии адсорбции (особенно над центром шестиугольника), чем та, где связь O-O перпендикулярна поверхности.
В данной работе представлены результаты расчета энергетических характеристик процессов адсорбции атома кислорода и диссоциативной адсорбции молекулы кислорода на поверхности гра-феновых нанолент типа «arm-chair». Проведено исследование электронных свойств таких структур в зависимости от концентрации атомов на поверхности и их геометрического расположения.
Расчеты электронного строения данных структур выполнены в рамках модели молекулярного кластера с граничными псевдоатомами с использованием квантово-химических полуэмпирических методов MNDO и PM3 [1]. Модель молекулярного кластера позволяет сравнительно легко описывать структурные дефекты, функциональные группы на поверхности кристаллов и т. д. Это способствует широкому ее применению в большинстве хемосорбционных и каталитических задач.
Выбор модели
В качестве геометрических моделей изучаемых нанолент выбраны кластеры, содержащие 3 шестиатомных цикла (гексагона) по ширине и от 3 до 6 элементарных ячеек (ЭЯ) по длине. Граничные разорванные химические связи замыкались атомами водорода, которые выбирались в качестве граничных псевдоатомов.
В работе рассмотрено три варианта положения атомарного и молекулярного кислорода (рис. 1). Положение адатомов кислорода в одном случае выбрано над центром С-С-связи, так как он является двухвалентным при формировании ионно-ковалентных связей (рис. 1, а), а в другом случае непосредственно над атомами углерода поверхности наноленты (рис. 1, в). Положение адсорбирующейся молекулы кислорода выбрано параллельно С-С-связи (рис. 1, б).
а б в
Рис. 1. Фрагменты структуры графеновой наноленты с указанием положения адсорбирующихся частиц:
а - атом кислорода над центром С-С-связи; б - молекула кислорода параллельно С-С-связи; в - атомы кислорода над атомами углерода
Кроме того, адсорбирующиеся частицы помещались в центральной части наноленты, чтобы уменьшить влияние граничных псевдоатомов.
Обсуждение результатов
Рассчитаны следующие основные электронно-энергетические характеристики: энергии граничных орбиталей ЕВЗМО и ЕшМО; энергия адсорбции Еад и ширина запрещенной щели Еg и изменение ширины запрещенный зоны ЛЕ^, обусловленное адсорбцией.
В результате адсорбции атом кислорода образует две химические связи с соседними атомами углерода, молекула кислорода диссоциирует на атомы, которые образуют по две химической связи с поверхностными активными центрами С, при этом происходит деформация поверхности наноленты так, что атомы углерода смещаются с положений равновесия наружу наноленты. В результате такой деформации связи С-С углеродных гексагонов, над которыми происходило присоединение О2, удлинились по сравнению с невозмущенными значениями в среднем на 9 %.
Анализ результатов показывает, что увеличение длины наноленты приводит к изменению энергетических характеристик процессов оксидирования углеродных нанолент. Величины энер-
гий верхней занятой и нижней вакантной МО зависят от длины наноленты: ЕВЗМО монотонно уменьшается, а ЕНВМО - увеличивается (рис. 2). С увеличением длины наноленты разница между энергиями граничных орбиталей убывает, что свидетельствует об увеличении перекрывания между валентной зоной и зоной проводимости и увеличении реакционной способности системы. То есть адсорбированные на поверхности частицы увеличивают сродство ГНЛ к другим атомам или молекулам.
о -1 -2 -3
ы
- J Fhbmo
L, ЭЯ
а
б
' MNDO ■РМЗ
ь,эя
в
Рис. 2. Энергия граничных орбиталей ГНЛ с адсорбирующимися частицами как функция длины:
а - один адатом кислорода; б - молекула кислорода; в - четыре адатома кислорода
Адсорбированные частицы можно рассматривать как дефект структуры - квантовую точку, поляризующую поверхность ленты. Квантовая точка на поверхности создает дополнительные локализованные квантовые состояния электронов, которые могут заполняться электронами системы. Согласно расчетам величина энергии запрещенной щели Еу ГНЛ с одним атомом кислорода возрастает в пределах ~1,1 ^ 1,4 эВ (MNDO) и ~0,89 ^ 1,32 эВ (РМ3). При адсорбции молекулы кислорода и ее диссоциации значения величины энергии запрещенной зоны колеблются в пределах значений ~1,05 ^ 1,68 эВ (MNDO) и ~0,79 ^ 1,55 эВ (РМ3). Величина энергии запрещенной щели Е^ у ГНЛ с четырьмя атомами кислорода составляет ~1,21 ^ 1,67 эВ (MNDO) и ~1,27 ^ 1,61 эВ (РМ3). Полученные данные свидетельствуют о полупроводниковых свойствах изучаемых нанолент.
5
L, ЭЯ
а
1,8
1,6
1,4
1,2
ее 1
eS 0,8 Ы
0,6 0,4 0,2 О
MNDO ■РМЗ
L, ЭЯ
Рис. 3. Энергия запрещенной зоны оксида ГНЛ как функция длины:
а - один адатом кислорода; б - молекула кислорода; в - четыре адатома кислорода
б
в
Следует заметить, что с увеличением длины наноленты энергия Eg уменьшается (см. рис. 3, а, б). Это говорит об увеличении «металличности» изучаемой структуры, что соответствует экспериментальным данным: бесконечная графеновая поверхность обладает проводящими свойствами. Однако следует отметить, что при адсорбции четырех атомов кислорода наблюдается отсутствие корреляции численных результатов, полученных разными полуэмпирическими методами (рис. 3, в). Данное расхождение, на наш взгляд, связано с параметризацией использованных полуэмпирических методов. При этом метод PM3 использует модифицированную формулу для межатомного взаимодействия [1], что приводит к увеличению энергии взаимодействия между молекулами кислорода, адсорбированными на поверхности ГНЛ. Сильная корреляция между молекулами O2 влияет на энергию первого электронного перехода, которая сопоставляется с величиной Eg
Зависимости энергии химической связи (энергии адсорбции) от длины наноленты представлены на рисунке 4. Для структуры, содержащей n адсорбированных атомов кислорода, энергия химической связи E (n), приходящаяся на один адсорбированный атом, рассчитывалась следующим образом:
Еад(п) = (E
total
egnr neo) / n,
где Etotal - энергия графеновой наноленты или ее фрагмента (с адсорбированными атомами кислорода на поверхности); EGNR - полная энергия чистой ГНЛ; ЕО - энергия изолированного атома кислорода.
Кривые на рисунке 4 наглядно демонстрируют, что энергия связывания атома кислорода колеблется около некоторого постоянного значения. При адсорбции одного атома кислорода энергия химической связи составляет в среднем -4,5 эВ в методе MNDO и -3,7 эВ в методе PM3. Отрицательные значения энергии адсорбции указывают на то, что система находится в устойчивом состоянии. При адсорбции четырех атомов кислорода энергия химической связи составляет в среднем -4,2 эВ в методе MNDO и -3,7 эВ в методе PM3.
3 4 5 6 7
ь, эя
б
Рис. 4. Энергия адсорбции кислорода как функция длины ГНЛ:
а - один адатом кислорода; б - молекула кислорода; в - четыре адатома кислорода (см. также с. 100)
а
2 3 4 5 6 7
-3,0 -,-1-,-,—
S "3,5 -
-4,5 -
-5,0 -
L/ЭЯ
в
Рис. 4. Окончание
Но при этом наблюдается слабое монотонное убывание величины Еад (рис. 4, в) с ростом длины ленты, что свидетельствует об увеличении устойчивости образованного оксида ГНЛ и уменьшении влияния границ ленты на процесс адсорбции. При адсорбции молекулы кислорода энергия химической связи составляет в среднем -0,03 эВ в методе MNDO и 0,25 эВ в методе PM3. То есть такое положение молекулы над поверхностью графеновой наноленты является энергетически менее выгодным. Данные результаты хорошо согласуются с аналогичными расчетами с помощью неэмпирических методов [9].
Заключение
На основе модели молекулярного кластера с использованием квантово-химических полуэмпирических методов MNDO и PM3 изучены электронные свойства графеновых нанолент с адсорбированными атомами и молекулами кислорода. Показана зависимость энергетических характеристик процесса адсорбции от конфигурации адатома и концентрации адатомов. В частности, наибольшая эффективность адсорбции наблюдается в тех случаях, когда адатом кислорода располагается над центром С-С-связи. Поэтому адсорбированные на поверхности ГНЛ молекулы кислорода диссициируют на отдельные адатомы, которые адсорбируются над центрами С-С-связей.
Из анализа результатов расчета следует, что энергия адсорбции адатома кислорода и молекулы кислорода слабо зависит от длины наноленты, а двух молекул кислорода монотонно возрастает с увеличением длины ГНЛ.
СПИСОК ЛИТЕРА ТУРЫ
1. Степанов, Н. Ф. Квантовая механика и квантовая химия / Н. Ф. Степанов. - М. : Мир, 2001. - 519 с.
2. Шамина, Е. Н. Влияние хиральности углеродных нанотрубок на диссоциативную адсорбцию молекулярного кислорода на поверхности / Е. Н. Шамина, Н. Г. Лебедев // Журнал физической химии. - 2015. -№5(89). - С. 806-811.
3. Шамина, Е. Н. Хиральный эффект адсорбции атомарного кислорода на поверхности углеродных нанотрубок / Е. Н. Шамина, Н. Г. Лебедев // Вестник Волгоградского государственного университета. Серия 1, Математика. Физика. - 2013. - № 1 (18). - С. 90-97.
4. Electronic, structural, and transport properties of Ni-doped graphene nanoribbons / V A. Rigo, T. B. Martins, J. R. Antonio da Silva [et al.] // Phys. Rev. B. - 2009. - Vol. 79. - Р. 075435.
5. Enhanced half-metallicity in edge-oxidized zigzag graphene nanoribbons / O. Hod, V. Barone, J. Peralta [et al.] // Nano Lett. - 2007. - Vol. 7 (8). - Р. 2295-2299.
6. Miller, M. Defect induced distortion of armchair and zigzag graphene and boron nitride nanoribbons / M. Miller, F. J. Owens // Chem. Phys. Lett. - 2013. - Vol. 570. - P. 42-45.
7. Son, Y. W. Energy Gaps in Graphene Nanoribbons / Y. W. Son, M. L. Cohen, S. G. Louie // Phys. Rev. Lett. -2006. - Vol. 9 (21). - P. 216803-216806.
8. Wang, Z. Adsorption of transition metal atoms (Co and Ni) on zigzag graphene nanoribbon / Z. Wang, J. Xiao, M. Li // Appl. Phys. A. - 2013. - Vol. 110 (1). - P. 235-239.
9. Yu, S. S. Oxidation of Graphene Nanoribbon by Molecular Oxygen / S. S. Yu, W. T. Zheng, J. Qing // IEEE Trans. Nanotechnol. - 2008. - Vol. 7 (5). - P. 628-635.
REFERENCES
1. Stepanov N.F. Kvantovaya mekhanika i kvantovaya khimiya [Quantum mechanics and quantum chemistry]. Moscow, Mir Publ., 2001. 519 p.
2. Shamina E.N., Lebedev N.G. Vliyanie khiralnosti uglerodnykh nanotrubok na dissotsiativnuyu adsorbtsiyu molekulyarnogo kisloroda na poverkhnosti [Chiral effect of the dissociative adsorption of molecular oxygen on the carbon nanotube surface]. Zhurnalfizicheskoy khimii [Russian Journal of Physical Chemistry A], 2015, no. 5 (89), pp. 806-811.
3. Shamina E.N., Lebedev N.G. Khiralnyy effekt adsorbtsii atomarnogo kisloroda na poverkhnosti ugleronykh nanotrubok [The chiral adsorption effect of atomic oxygen on the carbon nanotube surface]. Vestnik Volgogradskogo gosudarstvennogo universiteta. Seriya 1, Matematika. Fizika [Science Journal of Volgograd State University. Mathematics. Physics], 2013, no. 1 (18), pp. 90-97.
4. Rigo V.A., Martins T.B., Antonio da Silva J.R. et al. Electronic, structural, and transport properties of Ni-doped graphene nanoribbons. Phys. Rev. B, 2009, vol. 79, pp. 075435.
5. Hod O., Barone V., Peralta J. et al. Enhanced half-metallicity in edge-oxidized zigzag graphene nanoribbons. Nano Lett, 2007, vol. 7 (8), pp. 2295-2299.
6. Miller M., Owens F.J. Defect induced distortion of armchair and zigzag graphene and boron nitride nanoribbons. Chem. Phys. Lett., 2013, vol. 570, pp. 42-45.
7. Son Y. W., Cohen M.L., Louie S.G. Energy Gaps in Graphene Nanoribbons. Phys. Rev. Lett, 2006, vol. 9 (21), pp. 216803-216806.
8. Wang Z., Xiao J., Li M. Adsorption of transition metal atoms (Co and Ni) on zigzag graphene nanoribbon. Appl. Phys. A., 2013, vol. 110 (1), pp. 235-239.
9. Yu S.S., Zheng W.T., Qing J. Oxidation of Graphene Nanoribbon by Molecular Oxygen. IEEE Trans. Nanotechnol., 2008, vol. 7 (5), pp. 628-635.
INFLUENCE OF ADSORPTION OF ATOMS AND MOLECULES OF OXYGEN ON THE ELECTRONIC STRUCTURE OF GRAPHENE NANORIBBONS
Elena Nikolaevna Shamina
Candidate of Physical and Mathematical Sciences, Associate Professor, Department of Mathematics and Computer Science, Volgograd State Medical University shamina.alena@gmail.com
Pavshikh Bortsov Sq., 1, 400131 Volgograd, Russian Federation
Nikolay Gennadyevich Lebedev
Doctor of Physical and Mathematical Sciences, Professor, Department of Theoretical Physics and Wave Phenomena, Volgograd State University, nikolay. lebedev@volsu. ru, tf@volsu. ru
Prosp. Universitetsky, 100, 400062 Volgograd, Russian Federation
Abstract. The article presents the results of calculations of the energy characteristics of adsorption of oxygen atoms and molecules on graphene nanoribbons surface. The calculations
of electronic structures were carried out by semi-empirical MNDO and PM3 quantum-chemical methods using the molecular cluster model with boundary pseudoatoms. The main geometrical and electronic-energy characteristics of graphene nanoribbons with adsorbed oxygen atoms and molecules were calculated.
It is revealed that the changes of the frontier energies ^HOMO and £LUMO indicate that the properties of nanoribbons depend on the result of adsorption, in particular, on the increase of these systems reactivity.
It is also revealed that an oxygen molecule is weakly bound to the GNR surface with low adsorption energy. The configuration where the oxygen atom is located above the center of the C-C bond, or directly above the carbon atoms of the surface of the GNR has a stronger binding between an oxygen and surface.
Key words: molecular cluster model, quantum chemistry, MNDO, PM3, adsorption energy.