Моделирование углеродных наноструктур и их свойств методом молекулярной динамики тема автореферата и диссертации по физике, 01.04.07 ВАК РФ

На правах рукописи

ымг

БЕЛОВА Елена Эдуардовна

МОДЕЛИРОВАНИЕ УГЛЕРОДНЫХ НАНОСТРУКТУР И ИХ СВОЙСТВ МЕТОДОМ МОЛЕКУЛЯРНОЙ ДИНАМИКИ

Специальность 01.04 07 -физика конденсированного состояния

Автореферат

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

ООЗОБ45ВЗ

Москва-2007

003064583

Работа выполнена в Институте биохимической физики им Н.М. Эмануэля Российской академии наук

Научный руководитель

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

Чернозатонский Леонид Александрович

Официальные оппоненты:

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

Елецкий Александр Валентинович

кандидат физико-математических наук Осадчий Александр Валентинович

Ведущая организация:

Институт спектроскопии РАН

Защита состоится « 21» сентября 2007 г в 10 часов на заседании Диссертационного Совета Д 002.231.01 при Институте радиотехники и электроники РАН по адресу. 125009, Москва, ГСП-9, ул. Моховая, д 11, корп. 7

С диссертацией можно ознакомиться в библиотеке Института радиотехники и электроники РАН

Автореферат разослан «45"» августа 2007 г Ученый секретарь диссертационного совета

доктор физико-математических наук, профессор С.Н Артеменко

ОБЩАЯ ХАРАКТЕРИСТИКА РАБОТЫ

Актуальность работы состоит в том, что в качестве объекта исследования были выбраны углеродные наноструктуры, открытые в течение последних 20 лет, а именно углеродные нанотрубы, их многотерминальные соединения и графены Интерес к этим наноструктурам был вызван тем, что вследствие молекулярного масштаба они обладают новыми необычными физическими и химическими характеристиками В настоящий момент благодаря своим уникальным свойствам они считаются одними из наиболее перспективных материалов наноэлектроники и наномеханики

Реальная структура углеродных нанотруб была открыта в 1991 году Ииджимой (Iijima, 1991) Однако еще в 1952 году Радушкевич и Лушкинович (Радушкевич, 1952), а позднее в 1976 году, независимо, Эндо с коллегами (Endo, 1976) наблюдали эти нанообъекты в электронный микроскоп К настоящему моменту уже известно, что углеродные нанотрубы обладают уникальной прочностью, жесткостью, теплопроводностью и электропроводностью При этом, в зависимости от особенностей молекулярной симметрии и диаметра, они могут быть как полупроводниками, так и обладать металлической проводимостью В настоящее время углеродные нанотрубы начинают применять в нано-композитах и полевых транзисторах, источниках света и элементах нелинейной оптики Между тем количество их новых потенциальных применений растет с каждым днем

В начале 90-х экспериментаторам удалось получить не обычные углеродные нанотрубы, а их ветвистые аналоги (Zhou, 1995) Примерно в это же время, теоретически было показано, что такие многотерминальные соединения из углеродных нанотруб с различной проводимостью представляют собой идеальные наноэлектронные устройства Например, соединение из металлической и полупроводниковой нанотрубы выполняет функции выпрямительного нано-диода (Chico, 1996, Yao, 1999, Lambm, 2000) Область применения трехтерми-нальных нанотрубных соединений может быть уже гораздо больше, поскольку третья нанотруба будет выполнять функции переключателя или усилителя, преобразуя нанотрубное соединение в транзистор (Andnotis, 2001) К настоящему моменту технология синтеза сложных нанотрубных соединений достаточно развита, и их синтез воспроизводим, контролируем, обладает высоким процентом выхода

Совсем недавно, в 2004 году смогли выделить еще один углеродный на-ноструктурный материал - графен (Geim, 2004), который представлял собой отдельный лист графита Несмотря на то, что впервые такая 2D структура теоретически была рассчитана еще 60 лет назад (Wallace, 1947), получить и идентифицировать его удалось лишь на современном этапе развития технологий Интерес к исследованию этого материала можно объяснить сразу несколькими причинами Во-первых, наноразмерный графен является одним из наиболее перспективных материалов для наноэлектроники Согласно проведенным исследованиям в нем проявляется эффект баллистического транспорта при ком-

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

Зачастую в качестве одного из методов исследования наноструктур используют численное моделирование Преимущество этого метода, в связи с высокой ценой исследуемых материалов и необходимостью применения для определения их параметров передовых инструментальных средств, очевидно На сегодняшний день существуют две основные категории методов моделирования наносистем метод молекулярной динамики (МД) и ab initio методы Обычно, точность результатов, получаемых методами ab initio выше, чем дает МД, но следует учитывать, что они также и более громоздки с точки зрения затрачиваемого компьютерного времени Метод МД позволяет проводить исследования структуры, механических и тепловых свойств молекул, состоящих из нескольких тысяч атомов В то же время использование разработанных не так давно эмпирических потенциалов с квантово-механической концепцией типа связи между атомами позволяет существенно повысить точность получаемых результатов

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

Основными задачами исследования являлись

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

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

3 Классификация ковалентных соединений углеродных нанотруб

4 Изучение реакции димеризации пересекающихся под прямым углом углеродных нанотруб с учетом межмолекулярного взаимодействия между ними, определение условий реакции и ее механизма

5 Изучение механических свойств нанотрубных Y-соединений несимметричного типа и поиск возможности их использования в качестве структурных элементов наномеханических устройств

Основные положения, выносимые на защиту

1 Разработана программа МД с многочастичным модифицированным потенциалом Бреннера для вычисления геометрических (длина связи, валентные углы) и энергетических (энергия химических связей) параметров углеродных и углеводородных наноструктур Эту программу можно использовать для решения задачи по геометрической оптимизации структур, состоящих из нескольких тысяч атомов и для изучения изотермических процессов

2 Показано, что графены с периодически расположенными вакансиями являются устойчивыми наноструктурами Определены их геометрические и энергетические параметры

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

4 Ковалентные соединения нанотруб, образованные по механизму (2+2)-циклоприсоединение, классифицированы по типам I, V, Т, У и X

5 Реакция димеризации пересекающихся под прямым углом нанотруб с образованием Х-соединения протекает по механизму (2+2)-циклоприсоедине-ние при высоких давлениях и температурах

6 Несимметричные нанотрубные У-соединения с углом между нанотрубой-стволом и нанотрубой-веткой 30° обладают широкой областью упругих из-гибных деформаций ветки Ширина этой области не зависит от диаметра и длины составляющих соединение нанотруб, а связана с величиной угла между стволом и веткой Такие соединения могут работать в режиме наноп-ружин

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

Диссертационная работа выполнена в Институте биохимической физики им НМ Эмануэля РАН в соответствии с планом научно-исследовательских работ Института в рамках Российской научной программы «Низкоразмерные квантовые структуры» и Российской научно-технической программы «Актуальные направления в физике конденсированных сред» (направление «Фулле-рены и атомные кластеры») Диссертационная работа также поддержана Российским фондом фундаментальных исследований (грант № 05-02-17443)

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

1 предложены новые наноструктуры на основе графена, представляющие ценность для наноэлектроники,

2 проведена классификация ковалентных соединений углеродных нанотруб, образованных по механизму (2+2)-циклоприсоединение,

3 изучена реакция димеризации углеродных нанотруб с учетом межмолекулярного взаимодействия между ними;

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

5 установлена возможность перехода таких соединений при механическом воздействии в новое устойчивое состояние,

6 предложено новое наномеханическое устройство типа молекулярного храповика, продемонстрирована принципиальная схема его работы

Достоверность полученных результатов подтверждена частично путем сравнения с известными в настоящее время экспериментальными данными Достоверность также проверялась при помощи сравнения некоторых результатов расчета по нашей программе МД с результатами расчетов других авторов, полученных методами аЪ initio

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

Личный вклад соискателя. Все результаты представленные в работе получены соискателем лично, либо в соавторстве при его непосредственном участии

Апробация работы. Материалы работы докладывались на молодежных конференциях ИБХФ РАН - ВУЗы «Биохимическая физика» (Москва, 2004, 2005), 2-ой и 5-ой Международной конференции «Углерод фундаментальные проблемы науки, материаловедение, технология» (Москва, 2003, 2006), 6th and

7 Biennal International Workshop "Fullerenes and Atomic Clusters" (С -Петербург, 2003, 2005), NATO Advanced Research Workshop "Hydrogen Materials Science and Chemistry of Carbon Nanomatenals" (Судак, Украина, 2003), NATO Advanced Study Institute "Carbon Nanotubes From Basic Research to Nanotechnology" (Созополь, Болгария, 2005), E-MRS 2007 Spring Meeting (Strasbourg, France, 2007)

Результаты работы также обсуждались на научных семинарах в следующих организациях Институт биохимической физики им H M Эммануэля РАН, Физический институт им П H Лебедева, Институт радиотехники и электроники РАН, Московский государственный университет (химический факультет)

Публикации. Основные результаты опубликованы в 11 печатных работах, в том числе в 5 статьях в рецензируемых журналах и 6 тезисах докладов на международных конференциях

Объем и структура диссертации. Диссертация состоит из введения, 5 глав, заключения и списка литературы Работа изложена на 132 страницах, иллюстрирована 53 рисунками и 6 таблицами Библиография включает список из 125 работ

СОДЕРЖАНИЕ ДИССЕРТАЦИИ

Во введении обоснованы актуальность и научная новизна работы, сформулированы цель исследования и основные положения, выносимые на защиту, представлена практическая ценность диссертации

Первая глава содержит обзор современной литературы по теме диссертации А именно, в обзоре описаны известные в настоящее время экспериментальные способы получения многотерминальных нанотрубных соединений и возможные виды этих соединений Большое внимание уделено описанию структуры таких соединений и ее отличию от структуры одиночных нанотруб Кроме того, в обзоре приведены все известные к настоящему моменту данные о механических свойствах многотерминальных нанотрубных соединений В отдельный пункт в обзоре выделена информация о ковалентных нанотрубных соединениях, так как они также являлись объектом исследования в диссертации Приведены современные данные об условиях, при которых между отдельными углеродными нанотрубами возможно образование ковалентных связей, и предполагаемые механизмы протекания этой реакции В главе также приводятся сведения по более изученной реакции полимеризации фуллеренов, протекающей по механизму (2+2)-циклоприсоединение Кроме того, в обзоре уделено внимание совсем новому наноматериалу - графену, а именно его структуре, свойствам, современным методам синтеза и идентификации

Вторая глава посвящена теоретическим основам разработанной программы МД Суть метода МД заключается в решении уравнений движения атомов, полученных из второго закона Ньютона

d2r dt

где mn и rn - масса и пространственные координаты атома п, V - потенциал системы, а Vn - знак пространственного градиента

В качестве потенциала использовался эмпирический многочастичный потенциал Бреннера (Brenner, 1990) в модификации с включенным в него ван-дер-ваальсовым взаимодействием (Sinnott, 1998) Наш выбор потенциала был обусловлен тем, что согласно утверждению авторов и уже имеющейся информации по его применению, этот потенциал является реактивным, т е позволяет описывать процессы, сопровождающиеся образованием новых ковалентных связей и разрывом уже существующих связей Этот потенциал имел вид

^ = к )- В„УА (г,)+ У^ (гщ)] (2)

^ ' 1

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

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

/»Г -1)ехр(-,/2^('-, -*<")) (3)

/Д^ч -^Л-^щ^-К'")) (4)

где , и Д, - параметры морзовской функции, а именно глубина ямы, равновесное расстояние и постоянная соответственно, ^ - поправка, введенная Бреннером Отметим, что при ^ =2, потенциалы принимают вид обычных функций Морзе Функция /с(г9) - функция обрезания потенциалов, которая позволяет учитывать взаимодействие лишь ближайших соседей

Согласно Бреннеру величина поправки Вук1 рассчитывается по формуле

вди = (5)

где Ъук и Ьр1 - функции от валентных углов, а в функции учитывается наличие связей у атомов ближайшего окружения

Потенциал ван-дер-ваальсового взаимодействия имеет форму степенного потенциала Леннарда-Джонса, и гладко сшивался в программе с потенциалами ковалентного взаимодействия при помощи кубического сплайна Этот потенциал принимает ненулевое значение только в том случае, когда между атомами отсутствует ковалентное взаимодействие Кроме того, для сокращения времени вычисления потенциал занулялся, когда расстояние между атомами превышало 1 нм

Из формулировки МД (1) следует необходимость вычисления пространственного градиента от потенциала V Полученное аналитическое выражение для силы, действующей на атом п со стороны его локального окружения согласно приведенному выше формализму, имеет вид

К (г„,)+ V ПВП]И УА {гп])+ Вщк1Чп УА (г„,)- V „¥гАГ (г„ ))

1

+ (6)

Здесь учтено, что в силу симметричности определений (3)-(5) УК(г:1)= Уя{гр), Вчк1 = В]М, а, следовательно, такие же равенства будут распространяться и на градиенты от них.

Из-за высокой сложности уравнения (1) для его решения были использованы численные методы, а именно многошаговый алгоритм Бимана (Веетап, 1976), имеющий глобальную погрешность третьего порядка для координаты и второго порядка для скорости

*,♦,=*,+V,* + (4*;-^У/би (7)

V,« = V, + + - ^ )й / 6т, (8)

Недостатком этого алгоритма является его несамостартуемость, т е для его запуска необходим набор положений частиц не только в начальный момент времени, но и на следующем шаге Для получения необходимых начальных данных на первом этапе использовалась схема Рунге-Кутты четвертого порядка В качестве критерия выбора временного шага Ъ мы использовали эмпирическое правило, согласно которому флуктуации полной энергии системы не должны превышать нескольких процентов от флуктуации ее потенциальной энергии (Хеерман, 1990) Шаг интегрирования в нашем случае составлял 10"2, что соответствовало реальному времени ~ 10~16с

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

±^Лжвт (9)

здесь N — это количество атомов в наносистеме, кв - постоянная Больцмана, Т-абсолютная температура наносистемы, слева в равенстве записана кинетическая энергия наносистемы

Для решения задач по геометрической оптимизации углеродных наноструктур в программе использовался метод отжига (КлгкраШск, 1983) Этот метод основан на аналогии с процессом образования кристалла из расплавленного вещества Сначала выбранную систему «расплавляют», а затем проводят постепенное понижение ее температуры При этом при каждой температуре моделирование, проводимое в нашем случае методом МД, длилось достаточно долго, для того, чтобы система достигла устойчивого состояния

Для вычисления различных характеристик наносистем при постоянной температуре в программе использовался метод изокинетической МД (Хеерман, 1990) Согласно формализму которого кинетическая энергия наносистемы не меняет своего значения на протяжении всего процесса моделирования Для реализации этого условия на каждом шаге МД моделирования требовалось проводить перенормировку скоростей на коэффициент/?, имеющий следующий вид

здесь Т, - это температура теплового резервуара

В следующих главах диссертации приведены результаты, полученные при помощи программы МД

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

Обычно синтезированный графен располагается на подложке В компьютерном эксперименте ее функции выполнял идеальный графитовый лист При этом подложка и исследуемый графен были сдвинуты друг относительно друга по схеме Бернала А-В Учет подложки необходим, так как показало проведенное МД моделирование и реальный эксперимент (Shioyama, 2001, Meyer, 2007), поверхность одиночного графена, т е без подложки, уже не является плоской. Поскольку в эксперименте имеют дело с массивными (по сравнению с исследуемым графитовым листом) подложками, положение последнего при оптимизации было зафиксировано После получения, в результате оптимизации, состояния модифицированного графена с наименьшим значением потенциальной энергии, вычислялась энергия С-С связей графена, приходящаяся на один углеродный атом

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

возникновении вакансий не меняется При этом величина валентных углов вблизи вакансии отклоняется от значения 120°, характерного для идеального графена, а именно а = 132°, р = 127°, у = 113° (см Рис 1) Исследование таких графенов при температуре 1000 К показало устойчивость имеющихся в них вакансий, которые даже при таких высоких температурах сохраняли свою форму и никуда не смещались Из-за наличия вакансий, в

Рис. 1. Графен с одной ва- структуре графена имеются углеродные атомы с кансией Серьм цветом вы- ,

у двумя химическими связями (по три атома на каж-

делены атомы с двумя свя- J . г

зями а, р, у - углы, нерав- ДУЮ вакансию), что приводит к снижению суммарные 120° н°й энергии связи в графене по сравнению с его

бездефектным аналогом Некоторые из полученных при помощи программы МД оптимальных графеновых структур, с периодически расположенными вакансиями представлены на Рис 2 Здесь пары (n,l) - (m,k) определяют периоды вакансий в двух перпендикулярных направлениях На одну периодическую ячейку этих графенов приходится одна вакансия Количество атомов в ячейке графенов равно 35

атомов для наноструктуры с Рис 2(a) и 79 атомов для наноструктуры с Рис 2(6) Энергия связи углеродного атома в модифицированных графенах (без

учета взаимодействия с подложкой) равна 7,159 эВ и 7,278 эВ соответственно Для сравнения, вычисленная в рамках используемой программы МД энергия связи углеродного атома в идеальном графене составляет 7,373 эВ В результате взаимодействия графена с подложкой (графитовым листом) энергия каждого его атома увеличивается еще на 0,034 эВ и становится равной 7,407 эВ Последний результат хорошо согласуется с данными, полученными методами ab initio, согласно которым энергия связи углеродного атома в графите равна 7,40 эВ (Yin, 1984)

Полученные результаты по геометрической оптимизации графенов с периодически расположенными вакансиями были в дальнейшем использованы для получения их электронной структуры методом функционала локальной электронной плотности Так как требование наличия подложки существенно увеличило размер наносистемы, то провести геометрическую оптимизацию графенов методом функционала локальной электронной плотности не представлялось возможным В совместной работе (Чернозатонский, 2006) было показано, что линии из вакансий на графене сильно изменяют его спектр - в зависимости от расположения этих линий, а именно кратности трем разности индексов |л — /| (|т — Ц ), графеновые полосы между ними можно считать полуметаллическими или полупроводящими, при распространении носителей в направлении, перпендикулярном линиям При этом сами линии будут являться металлическими нановолноводами Создать подобные вакансии на графитовой поверхности возможно бомбардировкой листа ионами гелия или водорода (Han, 2003, Lehtinen, 2004)

Также в данной главе представлены результаты геометрической оптимизации графенов с линиями адсорбированных Н-атомных пар Установлено, что в связи с появлением в графеновой структуре дополнительных С-Н связей и вызванной этим sp2-sp3 трансформацией некоторых атомов углерода, имеет ме сто искривление поверхности графена, и зр3-гибридизированные атомы располагаются примерно на 0,1 нм выше остальных Длина С-С связи, в которой участвует хотя бы один зр3-гибридизированный атом составляет 1,54 А Такую же длину имеет и химическая связь в алмазе Длина связи С-Н равна 1,08 А Для сравнения, длина С-Н связи в этане согласно литературным источникам равна

Рис. 2. Графены с периодически расположенными вакансиями (а) (3,0)-(-3,3), (б) (5,0)-(-4,4) Атомы вокруг вакансий выделены серым цветом В рамках -элементарные прямоугольные ячейки

1,10 Â Также изменяются и валентные углы вблизи С-Н связей На Рис 3 (а) показаны углы ССС, неравные 120°, а именно углы а = 109° и р = 116° Углы

НСС имеют следующие значения 0 = 106° или ср = 102° (см Рис 3) Кроме того, было показано, что при наличии в графене нескольких линий адсорбированных Н-атомных пар, находящихся друг от друга на расстоянии не менее 0,5 нм, их взаимное влияние можно исключить

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

Графены с линиями адсорбированных Н-атомных пар представлены на Рис 4 Водородные линии в них находятся на расстоянии друг от друга 1,7 нм (а), 1 нм и 1,5 нм (б) и 3 нм (в) Для обозначения этих графенов можно использовать индексы элементарных ячеек (n,m) при помощи которых

(а)

(б)

(в)

Рис. 3. Графен с одной линией адсорбированных пар атомов водорода Атомы водорода выделены серым цветом а, р - углы ССС, неравные 120° 9 и <р - углы НСС

принято характеризовать углеродные нанотрубы На каждую периодическую ячейку графенов (7,0) и (7,7) приходится по одной- паре атомов водорода, а в ячейке более сложного (6,0)+(5,0) графена их две Количество углеродных атомов в элементарных ячейках представленных на Рис 4 графенов равно 28,44 и 28, соответственно.

Результаты геометрической оптимизации модифицированных водородом графенов (Рис 4), полученные при помощи программы МД, были в дальнейшем использованы для определения их электрон-

Рис. 4. Графены с линиями адсорбированных Н-атомных пар (а) (7,0), (б) (6,0)+(5,0), (в) (7,7) Атомы водорода выделены серым цветом В рамках - элементарные прямоугольные ячейки

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

спектр Например, графен (7,0) является полупроводником с шириной запрещенной зоны 0,86 эВ, а графен (7,7), напротив, является металлом Более сложная структура - графен (6,0)+(5,0) согласно расчетам оказывается металлом с диэлектрическими прослойками Недавно была продемонстрирована возможность экспериментального получения таких графитовых поверхностей с химически адсорбированными Н-атомными парами (Ногпекэег, 2006)

Четвертая глава диссертации посвящена ковалентным соединениям углеродных нанотруб Подобные соединения одиночных нанотруб между собой при помощи ковалентных связей могут служить аналогами многотерминальных нанотрубных соединений и широко использоваться в наноэлектронике

Вообще говоря, существует много способов соединить две пары атомов, принадлежащих соседним нанотрубам Но согласно литературным данным, наиболее ожидаемым механизмом их соединения, является (2 + 2) - циклоприсоединение, когда между нанотрубами образуется четырех-звенное углеродное кольцо (Рис 5)

По аналогии с принятой классификацией многотерминальных нанотрубных топологических соединений в главе были определены и представлены основные типы ковалентных нанотрубных соединений — I, V, Т, У и X (Рис 6)

Рис. 5. Четырехзвенное углеродное кольцо (выделено серым цветом), возникающее при полимеризации нанотруб по механизму (2+2)-циклоприсоединения

(а)

Рис. 6. Типы ковалентных нанотрубных соединений, образованных по механизму (2 + 2) -циклоприсоединения (а) димер 1-типа, (б) димер У-типа, (в) тример Т-типа, (г) тример У-типа на примере нанотруб (8,0), (д) димер Х-типа из нанотруб (9,0) и (5,5) Серым цветом в соединениях (а)-(г) выделено четырехзвенное углеродное кольцо, в (д) - нижняя нанотруба

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

собой прямоугольник, в котором длина межтрубных связей составляет 1,551,57 Á, а вторая пара сторон в зависимости от типа соединения принимает значения из промежутка 1,57-1,63 Á

Для каждого из соединений была вычислена когезионная энергия как разницу между энергией соединения Е и суммарной энергией свободных на-нотруб, входящих в его состав АЕ = Е-Ес Результаты расчетов для различных типов ковалентных соединений нанотруб (Рис 6) приведены в Табл 1 из которой видно, что ни одно из соединений выигрыша в энергии не даёт А, следовательно, для их образования требуется дополнительная энергия Небольшая энергетическая разница (~ 2 эВ) была зафиксирована у нанотрубных соединений I-, V- и Т-типов, а значит, их синтез наиболее вероятен Предполагается, что получить такие многотерминальные нанотрубные соединения можно при помощи электронного луча (Terrones, 2002) или при высоких давлениях и температурах (Чернозатонский, 2001)

Таблица 1. Энергетические характеристики ковалентных нанотрубных соединений (Рис 6)

тип соединения (число атомов) энергия отдельных нанотруб Ес ,эВ энергия (2 + 2) -соединения Е,эВ когезионная энергия соединения АЕ,эВ

1-тип (316) 2238,4 2236,9 -1,5

V-тип (316) 2238,4 2237,7 -0,7

Т-тип (446) 3172,5 3170,1 -2,4

Х-тип (806) 5786,7 5783,4 -3,3

Y-тип (474) 3357,6 3352,8 -4,8

В качестве примера в главе детально исследуется реакция димеризации пересекающихся под прямым углом двух металлических нанотруб (9,0) и (5,5)

с учетом межмолекулярного взаимодействия между ними Расположение углеродных нанотруб в начале компьютерного эксперимента представлено на Рис 7 Расстояние между фрагментами нанотруб длиной 5 нм в начальный момент времени составляло 0,4 нм Так как в реальном эксперименте

(а) ~ (б)

Рис. 7. Расположение углеродных нанотруб в начале компьютерного эксперимента Серыми стрелками показано направление действия постоянной внешней силы Р на нанотрубу (9,0)

нанотрубы лежат на подложке, то она также была включена в модель. Ее роль

выполнял квадратный кусочек графитового листа со стороной 2 нм Такие размеры подложки позволяли предотвращать проскальзывания нанотрубы (5,5) в ходе компьютерного эксперимента На верхнюю нанотрубу (9,0) на протяжении всего компьютерного эксперимента действовала постоянная внешняя сила Г Направление ее действие изображено на Рис 7. Непосредственно в процессе моделирования это воздействие реализовывалось добавлением константы / к ^-координате силы, действующей на атомы нанотрубы (9,0) Величина константы / была связана со значением внешней силы Т7 как / = Л7 / п, где п -это количество атомов в нанотрубе (9,0) Площадь поверхности этой нанотрубы определялась по формуле Б = лО ,где £> - ее диаметр Так как ¿) = 7,05 А, то 5 = 22,15 А Тогда давление Р, оказываемое на нанотрубу, можно вычислить по формуле Р = .Р / 5

В ходе проведения компьютерного эксперимента оказалось, что для исследуемой реакции характерны давления порядка нескольких ГПа С другой стороны размеры (9,0) были относительно небольшие, поэтому силовому воздействию подвергались не ряд ее верхних атомов, а все В противном случае в нанотрубе наблюдались значительные деформации Компьютерный эксперимент проводился при постоянной температуре по методу изокинетической МД Температура теплового резервуара выбиралась из области 11001900 К, так как такой температурный режим уже использовался для исследования этой реакции без учета межмолекулярного взаимодействия (Чер-нозатонский, 2001) Максимальное время наблюдения за наносистемой в компьютерном эксперименте составляло 3000 МД-шагов, т е 0,3 пс

Из полученной в результате моделирования Т-Р-диаграммы реакции (Рис 8) следует, что реакция протекает в температурном режиме 1600 — 1900 К при давлении 4-7 ГПа При этом наиболее вероятным механизмом ее протекания является (2+2)-циклоприсоединение в результате которого образуется ковалентное нанотруб-ное соединение Х-типа Кроме того, было установлено, что при давлениях свыше 5 ГПа зависимость температуры реакции от давления носит линейный характер Аппроксимация этой кривой до комнатных температур (300 К) дает значение давления в 23,7 ГПа, которое очень хорошо согласуется с экспери-

3,5 4,0 4,5 5,0 5,5 6,0 6,5 7,0 7,5 Р, ГПа

Рис. 8. Фазовая Т-Р-диаграмма димеризации пересекающихся под прямым углом нанотруб (9,0) и (5,5) с учетом межмолекулярного взаимодействия между ними На вставке представлено полученное в результате реакции димерное нанотрубное соединение Х-типа

ментальными данными по сверхтвердой фазе однослойных углеродных нанот-руб, обнаруженной при давлениях в 24 ГПа (Popov, 2002) Эту фазу, таким образом, можно объяснить ковалентным соединением нанотруб

В пятой главе диссертации исследуются механические свойства нано-трубных несимметричных Y-соединений, которые в настоящее время успешно синтезируются (Terrones, 2002, Gothard, 2004), и их возможные наномеханиче-ские приложения

В нанотрубном соединении типа несимметричный Y от нанотрубы - основного ствола отходит нанотруба — ветка (см Рис 9) Согласно правилу Эйлера для получения топологического соединения из углеродных нанотруб необходимо внедрение дефектов (октагонов, гептагонов, пентагонов) в совершенную гексагональную решетку В работе рассматривались соединения, в которых ветка со стволом образует острый угол ~ 30° В главе представлены результаты исследования их поведения под внешней нагрузкой

Схема проведенного модельного эксперимента изображена на Рис 9. Поскольку обычно нанотрубные соединения закрепляются, то при моделировании концы нанотрубы-ствола были зафиксированы В качестве внешней нагрузки, создающей деформацию, использовался графеновый кусочек Сначала он был расположен параллельно оси ствола на расстоянии 0,5 нм, а затем передвигался по направлению к стволу Из-за молекулярного взаимодействия кусочка с соединением, ветка отталкивалась от графена, как только расстояние между ними было меньше 0,32 нм Размеры кусочка были выбраны таким образом, чтобы, во-первых, избежать проскальзывания нанотрубы-ветки и добиться ее движения только в направлении, перпендикулярном плоскости графена, а во-вторых, минимизировать компьютерное время эксперимента Шаг сдвига графена составлял 0,06 нм После каждого перемещения наносистему (нанотрубное соединение и графеновый кусочек) оптимизировали методом отжига, фиксируя при этом текущее положение графенового кусочка, под которое и подстраивалась ветка При таком режиме удавалось избегать деформации кончика ветки и изменения диаметра самой ветви на протяжении всего компьютерного эксперимента Диаметр нанотруб в изучаемых соединениях составлял ~ 1 нм, это значение близко к экспериментальным величинам (Terrones, 2002)

Вначале были изучены особенности поведения под внешней нагрузкой соединения (10,0),(4,4)ВН401 Согласно принятой терминологии эта запись обозначает, что ствол соединения образован нанотрубой (10,0), ветка - (4,4), а область его ответвления содержит один октагон и четыре гептагона Полученные результаты показали, что под действием внешней нагрузки происходит

Рис. 9. Несимметричное У-соедине-ние в начале компьютерного эксперимента Серым цветом выделены зафиксированные области, а именно концы нанотрубы-ствола и графеновый кусочек На вставке изображена область ответвления с выделенными серым цветом дефектами

деформация соединения На Рис 10 приведен график изменения энергии соединения (10,0),(4,4)ВН401 в зависимости от деформации Энергия недефор-мированного соединения была принята за нуль и от этого уровня были отложены значения энергии соединения в различных напряженных состояниях По

оси абсцисс на графике отложен относительный изгиб е Если расстояние между кончиком ветки и стволом в недеформированном состоянии принять за г0 (см Рис 9), то величину относительного изгиба е можно определить как е = Аг/г0, где Дг=г0-г Здесь г- расстояние между кончиком ветки и стволом в напряженном состоянии Было установлено, что при относительно малых деформациях (я < 0,5) зависимость энергии от деформации е носит квадратичный характер

Е = \Е"е2, (11)

где коэффициент пропорциональности £"=14 эВ В этом случае нанотрубное соединение ведет себя подобно пружине, а именно, при удалении графенового кусочка ветка возвращается в свое первоначальное состояние Никаких остаточных деформаций после снятия внешней нагрузки в соединении не наблюдалось

Исходя из полученных результатов, был оценен эффективный модуль Юнга нанопружины (10,0),(4,4)ВН401 По закону Гука напряжение сг при упругой деформации пропорционально относительной деформации

<7 = Уе#Дг/г0=У^, (12)

где Ус#- эффективный модуль Юнга нанопружины С другой стороны, по определению, напряжение о - это сила, приходящаяся на единицу площади а = Г/5", где сила Р= -с1Е1(1г, 5- площадь действия силы Тогда эффективный модуль Юнга нанопружины можно выразить следующим образом

У0 (13)

Для случая нанопружины (10,0),(4,4)ВН401 было оценено, что силовое воздействие на ветку соединения создается, в основном, графеновыми атомами, заключенными в квадрат со стороной 0,8 нм Таким образом, было получено, что эффективный модуль Юнга нанопружины (10,0),(4,4)ВН401 составляет около 1,1 ГПа

При дальнейшем увеличении напряжения в соединении (10,0),(4,4)ВН401 имеет место необычный для макромира эффект, а именно при деформациях е ~ 0,76 происходит самопроизвольное прилипание ветки к

Рис. 10. Зависимость энергии деформации соединения (10,0),(4,4)ВН401 от величины относительного изгиба е

стволу за счет молекулярного притяжения (Рис 11) Было установлено, что при

критической деформации £сг=0,76 величина ван-дер-ваальсова взаимодействия между веткой и стволом составляет 4,091 эВ/нм Новое состояние, в котором ветка расположена параллельно стволу, является стабильным Более того, это новое состояние оказывается энергетически более выгодным, чем начальное недеформированное состояние Разница в энергии составила 0,4 эВ Расстояние между стволом и веткой в этом состоянии равно 0,32 нм Это значение равно величине ван-дер-ваальсовой щели в нанотрубных жгутах (Thess, 1996)

Таким образом, было установлено, что в результате внешнего механического воздействия происходит переход соединения (10,0),(4,4)ВН401 в новое устойчивое состояние На основании литературных данных (Lee, 2004) в главе предполагается возможность обратного перехода под действием электрического поля В таком случае несимметричные нанотрубные Y-соединения, пребывающие в зависимости от внешних факторов в одном из двух своих возможных состояний, могли бы послужить основой для механизмов памяти

Далее, в главе было оценено влияния различных характеристик соединения, таких как диаметр ветви и ствола, их длина, наличие внутренних нанотрубных структур, на его поведение под нагрузкой Было установлено, что ширина области упругих деформаций не зависит от диаметра и длины, составляющих соединение нанотруб, а связана с величиной угла между стволом и веткой соединения Для рассмотренных в главе соединений с углом 30°, причем как однослойных, так и многослойных, величина относительного изгиба ветки достигала 0,5 Кроме того, по мере увеличения геометрических размеров соединений снижается их жесткость При этом многослойные нанотрубные соединения являются более жесткими, чем их однойслойные аналоги Эффективный модуль Юнга рассмотренных нанотрубных соединений находился в интервале 0,3-1,1 ГПа Переход соединения в состояние с веткой, расположенной параллельно стволу, может не происходить, если ветка соединения не достаточно длинна Например, в соединении (19,0),(9,9)ВН6 длина ветки должна превышать 6,2 нм В противном случае для соединения характерны упругие деформации вплоть до величины относительного изгиба 0,85 При достаточной длине ветки переход соединения в новое состояние имеет место при величине относительного изгиба более 0,7, причем это значение не зависит ни от геометрии соединения, ни от количества в нем слоев

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

Рис. 11. Нанотрубное соединение (10,0),(4,4)ВН401 в состоянии с веткой, расположенной параллельно стволу

В предлагаемом в диссертации молекулярном храповике в качестве зубчатого колеса использовалась нанотруба (12,12) с четырьмя несимметричными

зубьями (Рис 12) Образовать на поверхности нанотрубы такие зубья могут, например, осажденные частицы катализатора (Gothard, 2004) В таком случае количество зубьев на колесе будет напрямую зависеть от диаметра используемой нанотрубы Регулировать вращение такого молекулярного зубчатого колеса можно было бы при помощи несимметричного Y-соединения При вращении колеса в одном направлении ветка соединения приподнимается, не мешая, тем самым, свободному ходу колеса После прохождения зуба ветка соединения возвращается в свое исходное положение (см Рис 13) При этом движение колеса в обратном направлении вследствие несимметричной формы нанотрубного соединения оказывается невозможным

В качестве кандидата на роль собачки храповика в диссертации было рассмотрено нанотрубное соединение (10,0),(4,4)ВН401 Из-за большого размера молекулярного зубчатого колеса в компьютерном эксперименте использовалась только его элементарная часть, а именно кусочек графена с одним несимметричным зубом Также для удобства система отсчета была перенесена в центр зубчатого колеса, и проводилось моделирование движение нанотрубного соединения в двух противоположных направлениях путем его сдвига вдоль выбранной оси с шагом 0,15 нм Такой режим обеспечил движение нанотрубного соединения (10,0),(4,4)ВН401 относительно элемента колеса с линейной скоростью 0,05 нм/пс, что в другой системе отсчета соответствовало вращению зубчатого колеса из нанотрубы (12,12) с угловой скоростью 0,03 пс"1 Такая величина угловой скорости вполне типична для моделей наномеханических устройств на основе нанотруб (Han, 1997)

Рис. 12. Молекулярное зубчатое колесо, представляющее собой нанотрубу (12,12) с четырьмя несимметричными зубьями на поверхности

Рис. 13. Схема движения нанотрубного соединения (10,0),(4,4)ВН401 относительно элемента колеса в положительном направлении оси Схема построена на основании результатов компьютерного эксперимента

На Рис 13 схематически представлено движения нанотрубного соединения (10,0),(4,4)ВН401 относительно элемента колеса в положительном направлении оси, что в другой системе координат соответствовало вращению молекулярного колеса против часовой стрелки Видно, что ветка соединения после прохождения элемента колеса возвращается в свое начальное недеформиро-ванное состояние Движение соединения в обратном направлении, что соответствует вращению молекулярного зубчатого колеса по часовой стрелке, блокируется, когда сила Т7, действующая на его ветку со стороны элемента колеса, равна критическому значению 10 нН Если же значение силы превышает эту критическую величину, то происходит резкое отталкивание ветки соединения, сопровождающееся частичными разрывами в ней химических связей Таким образом, движение молекулярного зубчатого колеса возможно только в одну сторону, поскольку его движение в обратном направлении ограничивается на-нотрубным соединением, выполняющим функции собачки храпового механизма

В заключении перечислены основные результаты и выводы диссертации

ОСНОВНЫЕ РЕЗУЛЬТАТЫ РАБОТЫ

1 Разработана программа молекулярной динамики, которая может применяться для моделирования различных углеродных и углеводородных наноструктур с большим количеством атомов

2 Предложены новые графеновые структуры с периодически расположенными вакансиями Определена их устойчивость, а также геометрические (длины химических связей, валентные углы) и энергетические (энергия химических связей) параметры

3 Рассчитаны новые графеновые системы с адсорбированными на их поверхность Н-атомными парами Показана их устойчивость и получены геометрические (длины химических связей, валентные углы) и энергетические (энергия химических связей) характеристики

4 Проведена классификация ковалентных соединений нанотруб, образованных по механизму (2+2)-циклоприсоединение, определены их возможные типы I, V, Т, У и X Установлена устойчивость всех типов соединений и для каждого типа определены параметры межтрубных связей и когезионная энергия

5 Детально изучена реакция димеризации углеродных нанотруб (9,0) и (5,5), пересекающихся под прямым углом, при высоких давлениях (4-7 ГПа) и температурах (1600 — 1900 К) с учетом межмолекулярного взаимодействия между ними. Получена Т-Р-диаграмма реакции.

6 Изучены механические характеристики несимметричных нанотрубных У-соединений и определено влияние на них геометрических параметров соединения и его слойности Показано, что в соединениях при внешней нагрузке на ветвь ее деформация является упругой вплоть до сокращения рас-

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

7 Показано, что при подобном механическом воздействии возможен переход несимметричного нанотрубного Y-соединения в устойчивое состояние с веткой, расположенной параллельно стволу Определены условия этого перехода

8 Предложено новое наномеханическое устройство - молекулярный храповик, в котором в качестве зубчатого колеса используется нанотруба с несимметричными зубьями, а собачкой является несимметричное нанотрубное Y-соединение Показана принципиальная схема работы этого прибора

СПИСОК ПУБЛИКАЦИЙ ПО ТЕМЕ ДИССЕРТАЦИИ

1 Mikheeva Е, Chernozatonskn L, Andriotis А N, Menon М Cap-polymerized nanotube junctions structures and transport properties // Fullerenes and Atomic Clusters book of abstracts, 6th Biennal International Workshop, St Petersburg, June 30- July 4, 2003 - p 99

2 Михеева E Э, Чернозатонский JIA, Астахова T Ю Механизм образования ковалентных соединений из углеродных нанотруб // Углерод фундаментальные проблемы науки, материаловедение, технология тез докл 2 межд конф , МГУ, Москва, 15-17 октября 2003г - С 153

3 Mikheeva Е Е, Chernozatonskn L А, Astahova Т Yu Covalent-bmdmg carbon nanotube simulation of formation mechanisms and energy characteristics // Hydrogen Materials Science and Chemistry of Carbon Nanomaterials Edited by T N Veziroglu, S Yu Zagmaichenko, D V Schur et al - Kluwer Academic Publishers, Dordrecht/Boston/London - 2004 - V 172-p 279-282

4 Белова E Э, Чернозатонский Л А Упругие свойства трехтерминальных на-нотрубных соединений // Биохимическая физика тез докл 4 молод конф ИБХФ РАН - ВУЗЫ, Москва, 25 - 26 октября 2004г - С 21

5 Belova Е , Chernozatonskn L A "Spring" behavior of "bough" CNT junctions // Fullerenes and Atomic Clusters book of abstracts, 7th Biennal International Workshop, St Petersburg, June 27- July 1, 2005 - p 256

6 Белова E Э, Чернозатонский Л А Компьютерное моделирование синтеза многотерминальных наноструктур путем ковалентного соединения углеродных нанотруб // Биохимическая физика тез докл 5 межд молод конф ИБХФ РАН - ВУЗЫ, Москва, 14-16 декабря 2005г - С 36-37

7 Belova Е, Chernozatonskn L A Mechanical properties of three-termmal nanotube junction determined from computer simulations // Carbon Nanotubes From Basic Research to Nanotechnology Edited by V N Popov, Ph Lambin -Springer - 2006-V 222-p 215-217

8 Чернозатонский Л A, Сорокин П Б , Белова Е Э Электронные свойства графена и дефектных сверхрешеток на его основе // Углерод фундаментальные проблемы науки, материаловедение, технология тез докл 5 межд конф , МГУ, Москва, 18-20 октября 2006г - С 40

9 Чернозатонский JIА , Сорокин П Б , Белова Е Э , Брюнинг Й, Федоров А С Сверхрешетки металл - полупроводник (полуметалл) на графитовом листе с вакансиями //ПисьмаЖЭТФ-2006-Т 84-С 141-145

10 Чернозатонский Л А, Сорокин П Б , Белова Е Э , Брюнинг Й, Федоров А С Сверхрешетки, состоящие из «линий» адсорбированных пар атомов водорода на графене //Письма ЖЭТФ-2007-Т 85-С 84-89

11 Belova Е, Chernozatonsku L A Mechanical properties of carbon nanotube bough junctions A theoretical study //Phys Rev B-2007-V 75-p 073412-1-4

Отпечатано в копицентр.е «СТ ПРИНТ» Москва, Ленинские горы, МГУ, 1 Гуманитарный корпус www stprint ш e-mail' zakaz@stpnnt ш тел 939-33-38 Тираж 100 экз Подписано в печать 13 08 2007 г

Введение.

Глава 1. Обзор литературы.

1.1 Синтез многотерминальных нанотрубных соединений.

1.2 Структура многотерминальных нанотрубных соединений и их классификация.

1.3 Графен.

1.4 Ковалентно-связанные углеродные нанотрубы.

1.5 Механические свойства многотерминальных нанотрубных соединений

Глава 2. Метод молекулярной динамики.

2.1 Многочастичный модифицированный потенциал Бреннера.

2.2 Расчет силы, действующей на атом в поле многочастичного модифицированного потенциала Бреннера.

2.3 Численный алгоритм Бимана. Начальные условия.

2.4 Геометрическая оптимизация методом отжига.

2.5 Метод изокинетической МД.

2.6 Блок-схема алгоритма.

Глава 3. Новые графеновые структуры.

3.1 Схема компьютерного эксперимента.

3.2 Графены с вакансиями.

3.3 Графены с адсорбированными Н-атомными парами.

3.4 Выводы главы.

Глава 4. Ковалентные соединения углеродных нанотруб.

4.1 Основные типы ковалентных нанотрубных соединений.

4.2 Моделирование реакции димеризации пересекающихся под прямым углом нанотруб (9,0) и (5,5).

4.2.1 Схема компьютерного эксперимента.

4.2.2 Результаты моделирования.

4.3 Выводы главы.

Глава 5. Механические свойства нанотрубных несимметричных Yсоединений и их приложения.

5.1 Поведение нанотрубных несимметричных Y-соединений под внешней нагрузкой.

5.1.1 Схема компьютерного эксперимента.

5.1.2 Нанотрубное соединение (10,0),(4,4)ВН401.

5.1.3 Нанотрубное соединение (19,0),(9,9)ВН6.

5.1.4 Двойное нанотрубное соединение

10,0),(4,4)ВН401 @( 19,0),(9,9)ВН6.

5.2 Применение нанотрубного несимметричного Y-соединения в качестве собачки молекулярного храповика.

5.2.1 Молекулярный храповик.

5.2.2 Схема компьютерного эксперимента.

5.2.3 Результаты компьютерного эксперимента.

5.3 Выводы главы.

Несколько последних десятилетий общественной жизни были ознаменованы появлением и бурным развитием нанотехнологии, науки занимающейся изучением свойств объектов и разработкой устройств порядка нанометра (10"9 м). Впервые о возможности создания механизмов атомарных размеров еще в 1960 г. объявил Ричард Фейнман в своей знаменитой лекции «Внизу полным-полно места.» [1]. Его идеи, казавшиеся тогда современникам фантастическими, на нынешнем этапе развития технологии выглядят вполне реальными. Появление таких методов исследования как атомно-силовая микроскопия и туннельная электронная микроскопия позволили не только увидеть атомную структуру материалов, но и проводить различные манипуляции. В настоящее время принято считать, что нанотехнология является следующим логическим этапом развития электроники и других наукоёмких производств.

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

Реальная структура углеродных нанотруб была открыта в 1991 году Ииджимой [2]. Однако еще в 1952 году Радушкевич и Лушкинович [3], а позднее в 1976 году, независимо, Эндо с коллегами [4, 5] наблюдали эти нанообъекты в электронный микроскоп. К настоящему моменту уже известно, что углеродные нанотрубы обладают уникальной прочностью, жесткостью, теплопроводностью и электропроводностью [6-8]. При этом, в зависимости от особенностей молекулярной симметрии и диаметра, они могут быть как полупроводниками, так и обладать металлической проводимостью. В настоящее время углеродные нанотрубы начинают применять в нанокомпозитах и полевых транзисторах, источниках света и элементах нелинейной оптики [6-8]. Между тем количество их новых потенциальных применений растет с каждым днем.

В начале 90-х экспериментаторам удалось получить не обычные углеродные нанотрубы, а их ветвистые аналоги [9]. Примерно в это же время, теоретически было показано, что такие многотерминальные соединения из углеродных нанотруб с различной проводимостью представляют собой идеальные наноэлектронные устройства. Например, соединение из металлической и полупроводниковой нанотрубы выполняет функции выпрямительного нанодиода [10-12]. Область применения трехтерминальных нанотрубных соединений может быть уже гораздо больше, поскольку третья нанотруба будет выполнять функции переключателя или усилителя, преобразуя нанотрубное соединение в транзистор [13-17]. Отметим, что к настоящему моменту технология синтеза сложных нанотрубных соединений достаточно развита, и их синтез воспроизводим, контролируем, обладает высоким процентом выхода.

Совсем недавно, в 2004 году смогли выделить еще один углеродный наноструктурный материал - графен [18], который представлял собой отдельный лист графита. Несмотря на то, что впервые такая 2D структура теоретически была рассчитана еще 60 лет назад [19], получить и идентифицировать его удалось лишь на современном этапе развития технологий.

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

Зачастую в качестве одного из методов исследования наноструктур используют численное моделирование. Преимущество этого метода, в связи с высокой ценой исследуемых материалов и необходимостью применения для определения их параметров передовых инструментальных средств, очевидно. На сегодняшний день существуют две основные группы методов моделирования наноструктур: метод молекулярной динамики (МД) и ab initio методы. Вообще говоря, методы ab initio позволяют получать более точные результаты, но они имеют ограничения на размер системы. Для решения представленных в работе задач по структурам с большим количеством атомов (~1000) наиболее приемлемым оказался метод МД. Для повышения точности получаемых результатов нами использовался эмпирический потенциал, относящийся к группе потенциалов с квантово-механической рассмотрением типа связи между атомами.

Научная новизна диссертационной работы состоит в следующем:

1. предложены новые наноструктуры на основе графена, представляющие ценность для наноэлектроники;

2. впервые проведена классификация ковалентных соединеий углеродных нанотруб, образованных по механизму (2+2)-циклоприсоединение;

3. впервые изучена реакция димеризации углеродных нанотруб с учетом межмолекулярного взаимодействия между ними;

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

5. установлена возможность перехода таких соединений при механическом воздействии в новое устойчивое состояние;

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

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

Основные положения, выносимые на защиту:

1. Разработана программа МД с многочастичным модифицированным потенциалом Бреннера для вычисления геометрических (длина связи, валентные углы) и энергетических (энергия химических связей) параметров углеродных и углеводородных наноструктур. Эту программу можно использовать для решения задачи по геометрической оптимизации структур, состоящих из нескольких тысяч атомов и для изучения изотермических процессов.

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

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

4. Ковалентные соединения нанотруб, образованные по механизму (2+2)-циклоприсоединение, классифицированы по типам: I, V, Т, Y и X.

5. Реакция димеризации пересекающихся под прямым углом нанотруб с образованием Х-соединения протекает по механизму (2+2)-циклоприсоедине-ние при высоких давлениях и температурах.

6. Несимметричные нанотрубные Y-соединения с углом между нанотрубой-стволом и нанотрубой-веткой 30° обладают широкой областью упругих изгибных деформаций ветки. Ширина этой области не зависит от диаметра и длины составляющих соединение нанотруб, а связана с величиной угла между стволом и веткой. Такие соединения могут работать в режиме нанопружин.

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

8. Модель нового наномеханического устройства типа молекулярного храповика, основными элементами которого являются нанотруба с зубьями и несимметричное нанотрубное Y-соединение.

Представленные в диссертации результаты были опубликованы в работах [20-24].

5.3 Выводы главы

В данной главе при помощи программы МД были впервые исследованы механические свойства несимметричных нанотрубных Y-соединений. Нами было показано, что в таких соединениях при внешней нагрузке на ветвь её деформация является упругой вплоть до сокращения расстояния между концом ветки и стволом в два раза. Указанная область упругих деформаций характерна для соединений с углом 30°. Причем как однослойных, так и многослойных. По мере увеличения геометрических размеров соединений снижается их жесткость. При этом многослойные нанотрубные соединения являются более жесткими, чем их однойслойные аналоги. Эффективный модуль Юнга рассмотренных нанотрубных соединений находится в интервале 0,3-1,1 ГПа.

С увеличением деформации в соединении, как правило, происходит переход в новое состояние. В этом состоянии ветка соединения располагается параллельно стволу и удерживается в таком положении благодаря ван-дер-ваальсову взаимодействию между ней и стволом. Переход в состояние с веткой, расположенной параллельно стволу, может не происходить, если ветка соединения не достаточно для такого соединения длинна. Например, в соединении (19,0),(9,9)ВН6 длина ветки должна превышать 6,2 нм. В противном случае для соединения характерны упругие деформации вплоть до величины относительного изгиба 0,85. Отметим, что даже при таких высоких деформациях sp -гибридизация атомов соединения не меняется. При достаточной длине ветки переход соединения в новое состояние имеет место при величине относительного изгиба более 0,7, причем это значение не зависит ни от геометрии соединения, ни от количества в нём слоёв.

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

Достаточно широкая область упругих деформаций нанотрубных соединений типа несимметричный Y позволяет применять их в качестве структурных элементов разнообразных наномеханических устройств. Например, в данной главе была показана принципиальная возможность их использования в качестве собачки молекулярного храповика. В основе оригинального молекулярного храпового механизма лежит нанотруба (12,12) с несимметричными зубьями, вращение которой с угловой скоростью 0,03 пс"1 возможно только против часовой стрелки. Движение нанотрубы в обратном направлении ограничивается нанотрубным соединением (10,0),(4,4)ВН401. Заметим, что при вращении нанотрубы против часовой стрелки это несимметричное соединение свободному ходу колеса нисколько не мешало.

Заключение

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

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

2. Предложены новые графеновые структуры с периодически расположенными вакансиями. Определена их устойчивость, а также геометрические (длины химических связей, валентные углы) и энергетические (энергия химических связей) параметры.

3. Рассчитаны новые графеновые системы с адсорбированными на их поверхность Н-атомными парами. Показана их устойчивость и получены геометрические (длины химических связей, валентные углы) и энергетические (энергия химических связей) характеристики.

4. Проведена классификация ковалентных соединений нанотруб, образованных по механизму (2+2)-циклоприсоединение, определены их возможные типы: I, V, Т, Y и X. Установлена устойчивость всех типов соединений и для каждого типа определены параметры межтрубных связей и когезионная энергия.

5. Детально изучена реакция димеризации углеродных нанотруб (9,0) и (5,5), пересекающихся под прямым углом, при высоких давлениях (4-7 ГПа) и температурах (1600 - 1900 К) с учетом межмолекулярного взаимодействия между ними. Получена Т-Р-диаграмма реакции.

6. Изучены механические характеристики несимметричных нанотрубных Y-соединений и определено влияние на них геометрических параметров соединения и его слойности. Показано, что в соединениях при внешней нагрузке на ветвь её деформация является упругой вплоть до сокращения расстояния между концом ветки и стволом в два раза. Вычислен эффективный модуль Юнга нанопружины из такого соединения.

7. Показано, что при подобном механическом воздействии возможен переход несимметричного нанотрубного Y-соединения в устойчивое состояние с веткой, расположенной параллельно стволу. Определены условия этого перехода.

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

1. Фейнман Р. Внизу полным-полно места: приглашение в новый мир физики / Р. Фейнман // Рос. Хим. Ж.- 2002.- № 46.- С. 4-8.

2. Iijima S. Helical microtubules of graphitic carbon / S. Iijima // Nature-1991.- v. 354.-p. 56.

3. Радушкевич JI.В. О структуре углерода, образующегося при термическом разложении окиси углерода на железном контакте / JI.B. Радушкевич, В.М. Лушкинович // Ж. Физ. Химии 1952.- т. 26.- С. 8893.

4. Endo М. Structural improvement of carbon fibers prepared from benzene / M. Endo, T. Koyama, Y. Hishima // Jpn. J. Appl. Phys. 1976,-v. 15,- p. 2073-2076.

5. Oberlin A. High resolution electron microscope observations of grafitized carbon fibers / A. Oberlin, M. Endo, T. Koyama // Carbon 1976.- v. 14.- p. 133-135.

6. Carbon Nanotubes: synthesis, structure, properties and applications / eds. M.S. Dresselhaus, G. Dresselhaus.- Springer, Berlin, 2001.

7. Елецкий A.B. Углеродные нанотрубки и их эмиссионные свойства / А.В. Елецкий // УФН- 2002.- т. 172.- С. 401-438.

8. Елецкий А.В. Сорбционные свойства углеродных наноструктур / А.В. Елецкий//УФН-2004.-т. 174.-С. 1191-1231.

9. Zhou D. Complex branching phenomena in the growth of carbon nanotubes / D. Zhou, S. Seraphin//Chem. Phys. Lett. 1995,- v. 238-p. 286-289.

10. Structural and electronic properties of bent carbon nanotube / Ph. Lambin, A. Fonseca, J.P. Vigneron et al. // Chem. Phys. Lett. 2000.- v. 245 - p. 8589.

11. Pure carbon nanoscale devices: nanotube heterojunctions / L. Chico, V.H. Crespi, L.X. Benedict et.al. // Phys. Rev. Lett. 1996.- v. 76 - p. 971-974.

12. Carbon nanotube intramolecular junctions / Zh. Yao, H.W.Ch. Postma, L. Balents et. al. //Nature. 1999.- v. 402.-p. 273-276.

13. Menon M. Carbon nanotube "T junctions": nanoscale metal-semiconductor-metal contact / M. Menon, D. Srivastava // Phys. Rev. Lett. 1997.- v. 79-p. 4453-4456.

14. Ballistic switching and rectification in single wall carbon nanotube Y junctions / A.N. Andriotis, M. Menon, D. Srivastava et al. // Appl. Phys. Lett. 2001.- v. 79.- p. 266-268.

15. Rectification properties of carbon nanotube "Y junctions" / A.N. Andriotis, M. Menon, D. Srivastava et al. //Phys. Rev. Lett. 2001.- v. 87- p. 0668021-4.

16. Transport properties of single-wall carbon nanotube Y junctions / A.N. Andriotis, M. Menon, D. Srivastava et al. // Phys. Rev. В 2002.- v. 65 - p. 165416-1-13.

17. Novel electrical switching behaviour and logic in carbon nanotube Y-junctions / P.R. Bandaru, C. Daraio, S. Jin et al. // Nat. Mater. 2005.- v. 4,-p. 663-666.

18. Electric field effect in atomically thin carbon films / K.S. Novoselov, A.K. Geim, S.V. Morozov et al. // Science 2004.- v. 306 - p. 666-669.

19. Wallace P.R. The band theory of graphite / P.R. Wallace // Phys. Rev. -1947.- v. 71- p. 622-634.

20. Сверхрешетки металл полупроводник (полуметалл) на графитовом листе с вакансиями / JI.A. Чернозатонский, П.Б. Сорокин, Е.Э. Белова и др. // Письма ЖЭТФ - 2006.- т. 84.- С. 141-145.

21. Сверхрешетки, состоящие из «линий» адсорбированных пар атомов водорода на графене / J1.A. Чернозатонский, П.Б. Сорокин, Е.Э. Белова и др. // Письма ЖЭТФ 2007.- т. 85.- С. 84-89.

22. Belova E. Mechanical properties of carbon nanotube bough junctions: A theoretical study / E. Belova, L.A. Chernozatoskii // Phys. Rev. В 2007.- v. 75.-p. 073412-1-4.

23. Arc-grown Y-branched carbon nanotubes observed by scanning tunneling microscopy (STM) / Z. Osvath, A.A. Koos, Z.E. Horvath et al. // Chem. Phys. Lett. -2002.- v. 365-p. 338-342.

24. Li J. Growing Y-junction carbon nanotubes / J. Li, C. Papadopoulos, J. Xu // Nature. 1999.- v. 402,-p. 253-254.

25. Y-junction carbon nanotubes / B.C. Satishkumar, P.J. Thomas, A. Govindaraj et al. // Appl. Phys. Lett.- 2000.- v. 77.- p. 2530-2532.

26. Deepak F.L. Synthetic strategies for Y-junction carbon nanotubes / F.L. Deepak, A. Govindaraj, C.N.R. Rao // Chem. Phys. Lett. 2001.- v. 345 - p. 5-10.

27. Li W.Z. Straight carbon nanotube Y junctions / W.Z. Li, J.G. Wen, Z.F. Ren // Appl. Phys. Lett.- 2001.- v. 79 p. 1879-1881.

28. Y-junction carbon nanotubes grown by in situ evaporated copper catalyst / B. Gan, J. Ahn, Q. Zhang et al. // Chem. Phys. Lett. 2001.- v. 333 - p. 2328.

29. Ting J.-M. Multijunction carbon nanotube network / J.-M. Ting, C.-C. Chang // Appl. Phys. Lett.- 2002.- v. 80.- p. 324-325.

30. Controlled growth of Y-junction nanotubes using Ti-doped vapor catalyst / N. Gothard, C. Daraio, J. Gaillard et al. // NanoLett.- 2004.- v. 4 p. 213217.

31. Chernozatonskii L. Three-terminal junctions of carbon nanotubes: synthesis, structures, properties and applications / L. Chernozatonskii // J. Nanopartical Res.-2003.- v. 5.-p. 473-484.

32. Choi Y.C. Synthesis of Y-junction single-wall carbon nanotubes / Y.C. Choi, W. Choi // Carbon.- 2005.- v. 43.-p. 2737-2741.

33. Three-terminal carbon nanotube junctions: Current-voltage characteristics / L.W. Liu, J.H. Fang, L. Lu et al. // Phys. Rev. В 2005.- v. 71.- p. 1554241-4.

34. Y-branching of single walled carbon nanotubes / P. Nagy, R. Ehlich, L.P. Biro et al. // Appl. Phys. A- 2000.- v. 70.- p. 481-483.

35. Molecular junctions by joining single-walled carbon nanotubes / M. Terrones, F. Banhart, N. Grobert et al. // Phys. Rev. Lett. 2002.- v. 89- p. 075505-1-4.

36. Heyning O.T. A low cost method for the direct synthesis of highly Y-branched nanotubes / O.T. Heyning, P. Bernier, M. Glerup // Chem. Phys. Lett.-2005.- v. 409.-p. 43-47.

37. Crespi V.H. Relations between global and local topology in multiple nanotube junctions / V.H. Crespi // Phys. Rev. В 1998.- v. 58,-p. 12671.

38. Chernozatonskii L.A. Classification of three-terminal nanotube junctions / L.A. Chernozatonskii, S.V. Lisenkov // Fullerenes, nanotubes, and carbon nanostructures 2004.- v. 12 -p. 105-109.

39. Fabrication and electric-field-dependent transport measurements of mesoscopic graphite devices / Y. Zhang, J.P. Small, W.V. Pontius et al. // Appl. Phys. Lett.- 2005.- v. 86.-p. 073104-1-3.

40. Electronic confinement and coherence in patterned epitaxial grapheme / C. Berger, Z. Song, X. Li, et al. // Science.- 2006,- v. 312 p. 1191-1196.

41. Shioyama H. Cleavage of graphite to grapheme / H. Shioyama // J. Mat. Sci. Lett.- 2001.- v. 20.- p. 499-500.

42. Viculis L.M. A chemical route to carbon nanoscroll / L.M. Viculis, J.J. Mack, R.B. Kaner // Science.- 2003.- v. 299 p. 1361.

43. Raman fingerprint of grapheme / A.C. Ferrari, J.C. Meyer, V.Scardaci et al. // http://arxiv.org/abs/cond-mat/0606284.-20Q6.

44. Chernozatonskii L.A. Polymerized nanotube structures new zeolites? / L.A. Chernozatonskii // Chem. Phys. Lett. - 1998.- v. 297.- p. 257-260.

45. A theoretical study of buckminsterfiillerene reaction products: Сбо+Сбо / D.L. Strout, R.L. Muny, C. Xu et al. // Chem. Phys. Lett. 1993.- v. 214.-p. 576-582.

46. New phases of Сбо synthesized at high pressure / Y. Iwasa, T. Arima, R.M. Fleming et al. // Science. 1994.- v. 264- p. 1570-1572.

47. A.M. Rao, P. Zhou, K.A. Wang et al. // Science. 1993.- v. 259 - p. 955.

48. Single-crystalline (KC6o)n: A conducting linear alkali fulleride polymer / S. Pekker, A. Janossy, L. Mihaly et al. // Science. 1994.- v. 265 - p. 10771078.

49. Eklund P.C. Fullerene polymers and fullerene polymers composites / P.C. Eklund, A.M. Rao. Springer, New York, 2000.- 395 p.

50. Макарова Т.А. Электронная структура фуллеренов и фуллеритов / Т.А. Макарова, И.Б. Захарова СПб.: Наука, 2001. - 70 с.

51. A new phase transition in the T-P diagram of C6o fullerite / I.O. Bashkin, V.I. Rashchupkin, A.F. Gurov et al. // J. Phys.: Condens. Matter. 1994.- v. 6.-p. 7491-7498.

52. Углеродные системы из полимеризованных нанотруб: кристаллическая и электронная структуры / JT.A. Чернозатонский, М. Менон, Т.Ю. Астахова и др. // Письма в ЖЭТФ 2001.- Т. 74 - С. 523-527.

53. Superhard phase composed of single-wall carbon nanotubes / M. Popov, M. Kyotani, R.J. Nemanich et al. // Phys. Rev. В 2002.- v. 65 - p. 033408-1-4.

54. Brenner D.W. Empirical potential for hydrocarbons for use in simulating the chemical vapor deposition of diamond films / D.W. Brenner // Phys. Rev. В 1990.- v. 42.-p. 9458-9471.

55. Mechanics of carbon nanotubes / D. Qian, G.J. Wagner, W.K. Liu et al. // Appl. Mech. Rev. 2002.- v. 55.-p. 495-533.

56. Srivastava D. Nanomechanics of carbon nanotubes and composites / D.

57. Srivastava, C. Wei, K. Cho // Appl. Mech. Rev. 2003.- v. 56 - p. 215-230. 61.Елецкий A.B. Механические свойства углеродных наноструктур и материалов на их основе / А.В. Елецкий // УФН - 2007.- Т. 177.- С. 233274.

58. Чернозатонский J1.A. Эффект слипания ветвей Y-соединения углеродных нанотруб / J1.A. Чернозатонский, И.В. Пономарева // Письма в ЖЭТФ 2003.- Т. 78,- С. 777-781.

59. A second-generation reactive empirical bond order (REBO) potential energy expression for hydrocarbons / D.W. Brenner, O.A. Shenderova, J. A. Harrison et al. //J. Phys.: Condens. Matter. 2002.- v. 14.-p. 783-802.

60. Tsai P.-C. Coalescence, melting, and mechanical characteristics of carbon nanotube junctions / P.-C. Tsai, Y.-R. Jeng, T.-H. Fang // Phys. Rev. В -2006.- v. 74.-p. 045406-1-10.

61. Попл Д.А. Квантово-химические модели / Д.А. Попл // УФН 2002.- Т. 172.-С. 349-356.

62. Hohenberg P. Inhomogeneous electron gas / P. Hohenberg, W. Kohn // Phys. Rev.- 1964,-v. 136-p. B864-B871.

63. Kohn W. Self-consistent equations including exchange and correlation effects / W. Kohn, L.J. Sham // Phys. Rev.- 1965,- v. 140.- p. A1133-A1138.

64. K0H В. Электронная структура вещества волновые функции и функционалы плотности / В. Кон // УФН - 2002.- Т. 172.- С. 336-348.

65. Koster G.F. Wave functions for impurity levels / G.F. Koster, J.C. Slater // Phys. Rev.- 1954.- v. 95.- p. 1167-1176.

66. Menon M. Nonorthogonal tight-binding molecular-dynamics scheme for silicon with improved transferability / M. Menon, K.R. Subbaswamy // Phys. Rev. В 1997.- v. 55.-p. 9231-9234.

67. Menon M. Generalized tight-binding molecular dynamics scheme for heteroatomic systems: Application to SimCn clusters / M. Menon // J. Chem. Phys.-2001.-v. 114,-p. 7731-7735.

68. Хеерман Д.В. Методы компьютерного эксперимента в теоретической физике, Пер. с англ. / Под ред. С.А. Ахманова М: Наука. Гл. Ред. физмат. лит., 1990.- 176 с.

69. Гулд X. Компьютерное моделирование в физике: В 2-х частях. Часть 1: Пер. с англ./ X. Гулд, Я. Тобочник М.: Мир, 1990. - 349 с.

70. Abel G.C. Empirical chemical pseudopotential theory of molecular and metallic bonding /G.C. Abel //Phys. Rev. В 1985.- v. 31.-p. 6184-6196.

71. Tersoff J. New empirical-model for the structural properties of silicon / J. Tersoff // Phys. Rev. Lett. 1986.- v. 56.-p. 632-635.

72. Tersoff J. New empirical-approach for the structure and energy of covalent systems/J. Tersoff//Phys. Rev. В 1988.- v. 37.-p. 6991-7000.

73. Tersoff J. Empirical interatomic potential for carbon, with applications to amorphous carbon / J. Tersoff// Phys. Rev. Lett. 1988.- v. 61.- p. 28792882.

74. Tersoff J. Modeling solid-state chemistry interatomic potentials for multicomponent systems / J. Tersoff // Phys. Rev. В - 1989.- v. 39 - p. 5566-5568.

75. Mechanical properties of nanotubule fibers and composites determined from theoretical calculations and simulations / S.B. Sinnott, O.A. Shenderova, C.T. White et al. //Carbon 1998.- v. 36.-p. 1-9.

76. Robertson D.H. On the way to fullerenes: molecular dynamics study of the curling and closure of graphitic ribbons / D.H. Robertson, D.W. Brenner, C.T. White//J. Chem. Phys. 1992.- v. 96.-p. 6133-6135.

77. Robertson D.H. Temperature-dependent fusion of colliding Сбо fullerenes from molecular dynamics simulations / D.H. Robertson, D.W. Brenner, C.T. White //J. Chem. Phys. 1995.- v. 99.-p. 15721-15724.

78. Ion pickup of large, surface-adsorbed molecules: a demonstration of the Eley-Rideal mechanism / E.R. Williams, G.C. Jones Jr., L. Fang et al. // J. Am. Chem. Soc. 1992.-v. 114.-p. 3207-3210.

79. Alfonso D.R. Hydrocarbon adsorption on a diamond (100) stepped surface / D.R. Alfonso, S.E. Ulloa, D.W. Brenner // Phys. Rev. В 1994.- v. 49.- p. 4948-4952.

80. Robertson D.H. Energetics of nanoscale graphitic tubules / D.H. Robertson, D.W. Brenner, J.W. Mintmire // Phys. Rev. В 1992.- v. 45.- p. 1259212595.

81. Garg A. Effect of chemical functionalization on the mechanical properties of carbon nanotubes / A. Garg, S.B. Sinnott // Chem. Phys. Lett. 1998.- v. 295.-p. 273-278.

82. Mao Z. Molecular dynamics simulations of the filling and decorating of carbon nanotubules / Z. Mao, A. Garg, S.B. Sinnott // Nanotechnology -1999.-v. 10.-p. 273-277.

83. Ni B. Chemical functionalization of carbon nanotubes through energetic radical collisions / B. Ni, S.B. Sinnott // Phys. Rev. В 2000.- v. 61 - p. R16343- R16346.

84. A combined computational and experimental study of ion-beam modification of carbon nanotube bundles / B. Ni, R. Andrews, D. Jacques et al. // J. Phys. Chem. В -2001.-v. 105.-p. 12719-12725.

85. Garg A. Interactions of carbon-nanotubule proximal probe tips with diamond and grapheme / A. Garg, J. Han, S.B. Sinnott // Phys. Rev. Lett. -1998.- v. 81.- p. 2260-2263.

86. Kang J.W. The electroemission of endo-fullerenes from a nanotube / J.W.Kang, H.J. Hwang // Nanotechnology 2004.- v. 15-p. 1825-1830.

87. Stable configurations of C2o and C2g encapsulated in single wall carbon nanotubes / L. Zhou, Z.Y. Pan, Y.X. Wang et al. // Nanotechnology 2006.-v. 17-p. 1891-1894.

88. Los J.H. Intrinsic long-range bond-order potential for carbon: performance in Monte Carlo simulations of graphitization / J.H. Los, A. Fasolino // Phys. Rev. В 2003.- v. 68.-p. 024107-1-14.

89. Verlet L. Computer "experiments" on classical fluids. I. Thermodynamical properties of Lennard-Jones molecules / L. Verlet // Phys. Rev. 1967.- v. 159.-p. 98-103.

90. Beeman D. Some multistep methods for use in molecular dynamics calculations / D. Beeman // J. Сотр. Phys. 1976.- v. 20,- p. 130-139.

91. Корн Г. Справочник по математике (для научных работников и инженеров)/ Г. Корн, Т. Корн. М.: Наука, 1973. - 832 с.

92. Яворский Б.М. Справочник по физике для инженеров и студентов вузов / Б.М. Яворский, А.А. Детлаф. М.: Наука, 1965. - 848 с.

93. Kirkpatrick K.S. Optimization by simulated annealing / K.S. Kirkpatrick, C.D. Gelatt Jr., M.P. Vecchi//Science. 1983.- v. 220-p. 671-680.

94. Ma J. Simulated annealing using the classical density distribution / J. Ma, J.E. Straub//J. Chem. Phys. 1994.- v. 101.-p. 533-541.

95. Ballone P. Simulated annealing of carbon clusters / P. Ballone, P. Milani // Phys. Rev. В 1990.- v. 42,- p. 003201-1-4.

96. Tomanek D. Growth regimes of carbon clusters / D. Tomanek, M. A. Schluter // Phys. Rev. Lett. -1991.- v. 67.-p. 002331-1-4.

97. Hollow diamonds: stability and elastic properties / G. Benedek, E. Galvani, S. Sanguinetti et al. // Chem. Phys. Lett. 1995,- v. 244,- p. 339344.

98. Kohler Th. Molecular-dynamics study of nitrogen impurities in tetrahedral amorphous carbon / Th. Kohler, G. Jungnickel, Th. Frauenheim // Phys. Rev. В 1999.- v. 60.-p. 010864-1-8.

99. Lu Z-Y. Structures and dynamical properties of Cn, Sin, Gen, and Snn clusters with n up to 13 / Z-Y. Lu, C-Z. Wang, K.-M. Ho // Phys. Rev. В -2000.- v.61.-p. 002329-1-6.

100. Hodak M. Ordered phases of fullerene molecules formed inside carbonnanotubes/M. Hodak, L.A. Girifalco//Phys. Rev. В 2003.- v. 67-p. 075419-1-4.

101. Allinger N.L. Conformational analysis. 130. MM2. A hydrocarbon force field utilizing VI and V2 torsional terms / N.L. Allinger // J. Am. Chem. Soc. 1977.- v. 99.-p. 8127-8134.

102. Haile J.M. Extensions of the molecular dynamics simulation method. II. Isothermal systems / J.M. Haile, S. Gupta // J. Chem. Phys. 1983.- v. 79.-p. 3067-3076.

103. Irradiation-induced magnetism in graphite: A density functional study /P.O. Lehtinen, A.S. Foster, Y. Ma et al. //Phys. Rev. Lett. 2004.- v. 93-p. 187202-1-4.

104. Ferromagnetic spots in graphite produced by proton irradiation / K.H. Han, D. Spemann, P.Esquinazi et al. // Ads. Mater. 2003.- v. 15 - p. 17191722.

105. Metastable structures and recombination pathways for atomic1. V Vhydrogen on the graphite (0001) surface / L. Hornekser, Z. Sljivancanin, W. Xu et al. // Phys. Rev. Lett. 2006.- v. 96.-p. 156104-1-4.

106. The structure of suspended graphene sheets / J.C. Meyer, A.K. Geim, M.I. Katsnelson et al. // Nature 2007,-v. 446,- p. 60-63.

107. Yin M. T. Structural theory of graphite and graphitic silicon / M. T. Yin, M. L. Cohen // Phys. Rev. В 1984.- v. 29.- p. 6996 - 6998.

108. Sofo J.O. Graphane: A two-dimensional hydrocarbon / J.O. Sofo, A.S. Chaudhari, G.D. Barber// Phys. Rev. В 2007.- v. 75.-p. 153401-1-4.

109. Localization and quantization in covalently bonded carbon nanotube junctions / F. Cleri, P. Keblineski, I. Jang et al. // Phys. Rev. В 2004.- v. 69.-p. 121412(R)-121415(R).

110. Electrical transport through carbon nanotube junctions created by mechanical manipulation / H.W.Ch. Postma, M. Jonge, Z. Yao et al. // Phys. Rev. В 2000.-v. 62,-p. R10653-R10656.

111. Catalytic growth of single-wall carbon nanotubes from metal particles / J.H. Hafner, M.J. Bronikowski, B.R. Azamian et al. // Chem. Phys. Lett. -1998.-v. 296.-p. 195-202.

112. Gas-phase catalytic growth of single-walled carbon nanotubes from carbon monoxide / P. Nikolaev, M.J. Bronikowski, R.K. Bradley et al. // Chem. Phys. Lett. 1999.-v. 313.- p. 91-97.

113. Лазерно-индуцированные эффекты в спектрах комбинационного рассеяния света в одностенных углеродных нанотрубках / С.Н. Бокова, В.И. Конов, Е.Д. Образцова и др. // Квантовая электроника 2003.- т. 33.- С. 645-650.

114. Моделирование плотности одноэлектронных состояний для одностенных нанотрубок из углерода и нитрида бора / А.В. Осадчий, Е.Д. Образцова, С.В. Терехов и др. // Письма в ЖЭТФ 2003.- т. 77.- С. 479-484.

115. Menon M. Covalent bonding between fiillerenes p. 229-240 from Fullerene Polymers and Fullerene Polymer Composites / M. Menon, K.R. Subbaswamy // Springer, 2000.

116. Crystalline Ropes of Metallic Carbon Nanotubes / A. Thess, R. Lee, P. Nikolaev et al. // Science 1996.- v. 273- p. 483-487.

117. A three-terminal carbon nanorelay / S.W. Lee, D.S. Lee, R. E. Morjan et al. // NanoLett. 2004.- v. 4.- p. 2027-2030.

118. Храповой механизм — Википедия Электронный ресурс. -Режим доступа: http://en.wikipedia.org/wiki

119. Molecular dynamics simulations of carbon nanotube-based gears / J. Han, A. Globus, R. Jaffe et al. // Nanotechnology 1997.- v. 8.- p. 95-102.