“ЖУРНАЛ РАДИОЭЛЕКТРОНИКИ” N 11, 2011

оглавление

УДК 537.311.322

квантовохимическое Моделирование структуры кремниевых фуллеренов

 

В. В. Филиппов, А Н. Власов

Липецкий государственный педагогический университет, кафедра физики

 

Получена 10 ноября 2011 г.

 

Аннотация. Представлены результаты оптимизации атомной и электронной структуры кластеров кремния со структурой подобной углеродным фуллеренам. Определены основные квантово-энергетические параметры, характеризующие электронные свойства исследуемых наноформ. Для исследования устойчивости кремниевых фуллеренов определены энергии связей гидрированных и инкапсулированных атомами металлов наноструктур.

Ключевые слова: кремний, наночастица, фуллерен, NDDO PM3.

Abstract. The results of optimization of the atomic and electronic structure of silicon clusters with a structure similar to carbon fullerenes. Determined the basic quantum-energy parameters characterizing the electronic properties of the investigated nano-forms. To investigate the stability of silicone fullerenes determine the binding energy of the hydrogenated and the encapsulated metal atoms nanostructures.

Keywords: silicon, nanoparticles, fullerenes, NDDO PM3.

 

Введение

 

Известно, что одним из перспективных классов материалов нанотехнологий являются фуллерены [1]. После открытия фуллерена С60 в 1985 г. [2] начали появляться обоснованные предположения, что кремний и германий также могут образовывать энергетически устойчивые сфероидальные геометрические структуры. Относительно недавно японские ученые сообщили о том, что им удалось получить устойчивые кластеры из 12 и 18 атомов кремния с замкнутой структурой [3]. Стабилизации кремниевых наноструктур удавалось добиться путем помещения внутрь сфероидальной наночастицы атома переходного металла. Структура полученных кластеров представляла собой шестиугольную призму или икосаэдр в случае кластера МеSi12 ,  в случае кластера Me2Si18 – две шестиугольные призмы с общим основанием.

Необходимо отметить, что кремний является самым распространенным материалом современной электроники. Естественно полагать, что различные наноформы кремния можно использовать в качестве структурных элементов наноэлектронных транзисторов, выпрямляющих элементов, а также одноэлектронных приборов. Соответственно, у инженеров и ученых наблюдается особый интерес, направленный на экспериментальные исследования кремниевых кластеров [4-6], а также на теоретическое моделирование известных кластеров и прогнозирование их свойств [7-10]. В настоящее время уже получены устойчивые модификации некоторых кремниевых кластеров [4-6], а также фуллереноподобные структуры [11-13], кремниевые нанотрубки [1415] и однослойные пленки силицена [16]. Однако, физико-химические свойства указанных кремниевых наноформ исследованы не достаточно.

В данной работе представлены результаты оптимизации атомной и электронной структуры кластеров кремния со структурой подобной углеродным фуллеренам. Определены основные квантово-энергетические параметры, характеризующие электронные свойства исследуемых наноформ.

 

Методика расчета 

 

В качестве расчетных методов выбраны метод силового поля ММ+ (для оптимизации структуры остова) [17] и квантовохимический метод NDDO РМ3-UHF (для расчета электронных свойств) [18, 19]. Выбор именно данных методов обусловлен хорошей степенью их надежности и воспроизводимости результатов [20-24]. Для исследования устойчивости кремниевых фуллеренов определены также энергии связей гидрированных и инкапсулированных структур атомами переходных металлов. Во всех рассмотренных ниже случаях мультиплетность структуры определялась как минимальная, поскольку в приведенных выше экспериментальных работах по кремниевым наноформам не отмечалось о проявлении ими каких-либо магнитных свойств.

В методе силового поля ММ+ энергии заряд диполь и диполь-диполь оцениваются с точностью до членов . Вычисление потенциальной энергии растяжения валентных связей в ММ+ использует разложение до кубического члена:

          .                                       (1)

Энергия деформации валентных углов вычисляется в ММ+ с включением члена шестого порядка:

          .             (2)

Поправка, учитывающая изменение связи при изменении валентного угла вычисляется виде суммы по всем связям и валентным углам по формуле:

          .                                      (3)

Поправка, учитывающая внеплоскостные деформации вычисляется в виде суммы по всем группам атомов, образующих тригональную плоскую структуру по формуле:

          .                     (4)

Видно, что для каждого угла функция энергии аналогична функции изгиба валентных углов, но в качестве независимой переменной в данном случае используется угол, образуемый связью с плоскостью, в которой расположены атомы в недеформированном состоянии.

Энергия ван-дер-ваальсовых взаимодействий, описывающая притяжение атомов на больших расстояниях и отталкивание на близких в ММ+ определяется выражением:

          ,                        (5)

где , а  и  – значения ван-дер-ваальсовых радиусов, а суммирование производится по всем валентно не связанным атомам.

В указанных выше выражениях константы , , , , , , ,  ‑ экспериментальные параметры, определяемые с точностью до постоянного множителя, зависящего от выбора системы единиц.

В методе РМ3 полная энергия кластера вычисляется следующим образом:

,                       (6)

          ,                                               (7)

          ,                             (8)

          .                                                (9)

Здесь – матрица электронной плотности,  – оператор кинетической энергии электрона и энергии его притяжения к ядрам.

Заряд на атомах определялся по методу Малликена:

                                          (10)

где  – заряд ядра атома ,  – электронная заселенность орбитали  на атоме ,  – интеграл перекрывания орбиталей  и  атомов  и ,  – половина от электронной заселенности перекрывания орбиталей  и  атомов  и .

 

Результаты расчетов

 

В приведенной таблице 1 представлены геометрические структуры оптимизированных методом ММ+ некоторых кремниевых фуллеренов, указаны классы симметрии структур, а также рассчитанные средние длины связей. В таблице 2 отражены результаты квановохимических расчетов энергетических параметров исследуемых структур.

 

Таблица 1. Атомная структура оптимизированных кремниевых фуллеренов.

Si20 , cимм. – С1, <а>=2.35

Si28, cимм. – С1, <а>=2.34

Si36, cимм. – С2, <а>=2.33

Si40, cимм. – С1, <а>=2.34

Si50, cимм. – D5H, <а>=2.33

Si60, cимм. – СI, <а>=2.34

 

Таблица 2. Рассчитанные квантово-энергетические параметры кремниевых фуллеренов.

Структура

2S+1

Есвязи, eV

Есвязи/ат., eV

ЕHOMO, eV

DЕ, eV

qMe, e

Si20

1

-77.23

-3.86

-7.46

2.66

-

H20Si20

1

-143.29

-3.58

-9.18

6.71

-

CuSi20

2

-90.94

-4.33

-7.90

3.26

+0.2

CaSi20

1

-110.40

-5.26

-7.71

3.11

-0.46

Si28

1

-114.65

-4.09

-7.96

3.72

-

H28Si28

1

-199.03

-3.55

-8.80

6.21

-

CuSi28

2

-120.91

-4.17

-7.55

3.04

+0.56

CaSi28

1

-141.79

-4.89

-7.51

3.03

-0.05

Si36

1

-149.56

-4.15

-7.97

3.72

-

H36Si36

1

-254.06

-3.53

-8.75

6.07

-

CuSi36

2

-154.97

-4.18

-8.00

3.8

-0.03

CaSi36

1

-170.72

-4.62

-7.62

3.00

+0.2

Si40

1

-166.18

-4.15

-8.03

3.82

-

H40Si40

1

-281.70

-3.52

-8.59

5.89

-

Si50

1

-209.94

-4.20

-7.67

3.06

-

H50Si50

1

-350.08

-3.50

-8.57

5.79

-

Si60

1

-255.25

-4.25

7.98

3.88

-

H60Si60

1

-427.00

-3.56

8.57

5.73

-

 

 

Выводы и результаты

 

Выполненные расчеты по оптимизации атомной структуры кремниевых фуллеренов и определению квантово-энергетических параметров позволяют сделать следующие нижеуказанные выводы.

Атомный остов составлен на основе пятиугольников и шестиугольников с типом связи sp2, не характерным для кристаллических кремниевых структур. Среднее межатомное расстояние кремниевых фуллеренов примерно соответствует длине связи в кристаллическом кремнии.

Все представленные наноструктуры являются энергетически устойчивыми, поскольку имеют отрицательное значение энергии связи. С увеличением числа атомов кремния наблюдается незначительный рост удельной энергии связи. Увеличения энергии связи можно добиться, помещая в центр структуры металлические атомы или насыщая поверхностные атомы кремния водородом.

Атом металла, помещенный в центр фуллерена, может проявлять как донорные, так и акцепторные свойства. В зависимости от структуры заполнения электронами атомных орбиталей. Так для меди (Cu) 3-й энергетический уровень заполнен практически полностью, один валентный электрон слабо связан с атомным остовом, вследствие чего медь проявляет донорные свойства. Напротив, для кальция (Ca) возможно проявление акцепторных свойств, в связи с полным заполнением 4s подуровня, очевидно, при малом расстоянии относительно кремния начинает заполняться 4р подуровень; с увеличением расстояния от кремния кальций начинает проявлять донорные свойства. Характерно, что энергия связи для CuSin меньше чем для структур CaSin, это можно объяснить большей разницей электроотрицательностей между кремнием (1.95) и кальцием (1.0) по отношению к кремнию (1.95) и меди (1.7) [25].

Инкапсулирование исследуемых структур металлическими атомами приводит к повышению , однако, для структуры с 36 атомами кремния при добавлении меди данная зависимость нарушается поскольку в этом случае металл практически не связан с кремниевым остовом. При гидрировании кремниевого фуллерена вакантные связи насыщаются, что способствует значительному понижению верхнего заполненного уровня и увеличению разницы  и возрастанию работы выхода.

Таким образом, меняя структуру фуллерена, можно получать материалы с заведомо необходимыми свойствами, в частности величину работы выхода, электронного сродства и ширину  щели т.е. аналога запрещенной зоны объемного материала.

 

 

Литература

 

1.       Пул Ч., Оуэнс Ф. Нанотехнологии.  – М.: Техносфера, 2006. – 336 с.

2.       Керл Р.Ф. Истоки открытия фуллеренов: эксперимент и гипотеза // Успехи физич. наук. – 1998. – Т. 168. – № 3 – С. 331-342.

3.       Huira H., Miyazaki T., Kanayama T. Formation of Metal-Encapsulating Si Cage Clusters // Phys. Rev. Lett. – 2001. – Т. 86. – №9. – P. 1733-1736.

4.       Hoffmann M.A., Wrigge G., Issendorff B.V., Müller J., Ganteför G., Haberland H. Ultraviolet photoelectron spectroscopy of  to  // Eur. Phys. J. D. – 2001. – V. 16. ‑ № 1-3. – P. 9-11.

5.       Beck T.L., Berry R.S. The interplay of structure and dynamics in the melting of small clasters // J. Chem. Phys. – 1998. – V. 88. ‑ № 6. – P. 3910–3922.

6.       Sugano S. Microcluster Physics. // Berlin: Springer, 1991. –375 pp.

7.       Курганский С.И., Борщ Н.А. Геометрическая структура и спектральные характеристики электронных состояний кремниевых наночастиц // Физика и техн. полупр. – 2004. – Т. 38. – №5. – С. 580–584.

8.       Борщ Н.А., Переславцева Н.С., Курганский С.И. Атомная и электронная структура кремниевых и кремний-металлических наночастиц Si20, Si20, NaSi20 и KSi20 // Физика и техн. полупр. – 2006. – Т. 40. – № 12. – С. 1457–1462.

9.       Мелешко В.П., Мороков Ю.Н., Швейгерт В.А. Структура водородосодержащих кремниевых кластеров. Малые кластеры // Журн. структ. хим. – 1999. – Т. 40. – № 1. – С. 13-20.

10.   Yoo S., Zeng X.C. Structures and relative stability of medium-sized silicon clusters. IV. Motif-based low-lying clusters Si21 – Si30 // J. Chem. Phys. – 2006. – V. 124.№ 5. – P. 054304.

11.   Kumar V. Novel caged clusters of silicon: Fullerenes, Frank–Kasper polyhedron and cubic // Bull. Mater. Sci. ‑ 2003. ‑V. 26. ‑ № 1. ‑ P. 109–114.

12.   Pei Y., Gao Yi, Zeng X.C. Exohedral silicon fullerenes: SiNPtN/2 20≤ N ≤ 60 // J. Chem. Phys. – 2007. – V. 127. №4. – P. 044704.

13.   Gao Yi., M4Si28 (M = Al,Ga): Metal-encapsulated tetrahedral silicon fullerene // J. Chem. Phys. – 2005. – V. 123. ‑ № 20. – P. 204325.

14.   Park M.-H., Kim M.G., Joo J., Kim K., Kim J., Ahn S., Cui Yi, Cho J. Silicon nanotube battery anodes // Nano Lett. – 2009. – V. 9. – №. 11. – Р. 3844-3847.

15.   Мазуренко Е.A., Дорошенко М.Н., Герасимчук А.И. Синтез, свойства и моделирование кремниевых и германиевых нанотрубок // Украинский химический журнал. – 2008. – Т. 74. – № 11. – С. 3-15.

16.   Powell D., Silicene: It could be the new graphen // Science News. – 2011. ‑ V. 179. ‑ № 9. – P. 14.

17.   Соловьев М.Е., Соловьев М.М. Компьютерная химия. – М.: Солон-Пресс, 2005. – 536 с.

18.   Stewart J.J.P. Optimization of parameters for semiempirical methods // J. Comput. Chem. – 1989. – V. 10. – №2. – P. 209-264.

19.   Stewart J.J.P. Optimization of Parameters for Semi-Empirical Methods III-Extension of PM3 to Be, Mg, Zn, Ga, Ge, As, Se, Cd, In, Sn, Sb, Te, Hg, Tl, Pb, and Bi // J. Comput. Chem. – 1991. – V. 12. – №3. – P. 320-341.

20.   Чернозатонский Л.А. Новый класс диоксидных нанотруб МO2 = Si, Ge, Sn, Pb) из "квадратных" решеток атомов – их структура и энергетические характеристики // Письма в ЖЭТФ. – Т. 80. ‑ № 10. – С. 732-736.

21.   Курганский С.И., Борщ Н.А. Геометрическая и электронная структура кремниевых и кремниево-металлических наночастиц. // Изв. РАН. Сер. физич. – 2004. – Т. 68. – № 7. – С. 1023-1025.

22.   Филиппов В.В., Переславцева Н.С., Курганский С.И. Квантовохимическое моделирование структуры напряженных нанокристаллов кремния на германиевой подложке // Изв. РАН. Сер. физич. – 2008. – Т. 72. ‑ №9. – С. 1314-1316.

23.   Филиппов В.В., Власов А.Н. Моделирование электронных свойств кремниевых наночастиц с плотной атомной упаковкой // Известия вузов. Физика. – 2010. – № 1 – С. 70-75.

24.   Власов А.Н., Филиппов В.В. Квантово-энергетические и кинетические свойства материалов кремниевой наноэлектроники на основе кластеров Si2-Si10 [Электронный ресурс] // Журнал радиоэлектроники (электронный журнал). – 2011. – № 8. – Режим доступа: http://jre.cplire.ru/jre/aug11/6/text.html.

25.   Бацанов С.С. Структурная химия. Факты и зависимости. – М.: Диалог–МГУ, 2000. – 292 с.