Наноэлектроника лит-ра / dragunov
.pdf
контр-электрода. Затем между наконечником и контр-электродом подавалось переменное напряжение, что вызывало колебания жгута. При этом регистрировалась резонансная частота и амплитуда колебаний. Для оценки механических характеристик авторы [72] пользовались упрощенным соотношением, полагая, что резонансная частота колебаний
|
2 |
|
|
|
|
|
|
|
|
Ed |
|
|
|
fi |
i |
D |
. |
|||
8 L2 |
|
|||||
|
|
|
|
|
|
|
Результаты измерений резонансной частоты для жгутов УНТ с различными параметрами приведены в табл. 5.4. Размеры жгутов определялись с использованием просвечивающего электронного микроскопа.
|
|
|
|
|
|
|
Т а б л и ц а |
5.4 |
|
|
|
|
|
|
|
|
|
||
Жгуты незаполненных нанотрубок |
Жгуты заполненных нанотрубок |
||||||||
|
|
|
|
|
|
|
|
|
|
D, нм |
L, мкм |
f, МГц |
(Ed/ )1/2, м/c |
D, нм |
|
L, мкм |
f, МГц |
(Ed/ )1/2, м/c |
|
16 |
1.635 |
12.46 |
15069 |
14 |
|
1.693 |
13.2 |
19465 |
|
|
|
|
|
|
|
|
|
|
|
16 |
2.354 |
6.82 |
16572 |
16 |
|
2.405 |
7.96 |
20962 |
|
17 |
1.705 |
14.94 |
17842 |
18 |
|
2.616 |
8.04 |
22603 |
|
|
|
|
|
|
|
|
|
|
|
17 |
1.819 |
9.84 |
13452 |
19 |
|
1.556 |
21.56 |
19847 |
|
19 |
2.429 |
6.6 |
14649 |
20 |
|
2.682 |
6.32 |
16607 |
|
|
|
|
|
|
|
|
|
|
|
23 |
1.653 |
17.44 |
15071 |
20 |
|
2.733 |
7.66 |
20766 |
|
|
|
|
|
|
|
|
|
|
|
24 |
3.004 |
5.16 |
13754 |
21 |
|
3.442 |
4.56 |
18215 |
|
27 |
2.006 |
11.52 |
12318 |
23 |
|
5.916 |
1.52 |
16796 |
|
|
|
|
|
|
|
|
|
|
|
53 |
7.027 |
1.42 |
9391 |
35 |
|
6.831 |
1.64 |
15757 |
|
67 |
10.273 |
0.9 |
10046 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
68 |
11.615 |
0.52 |
7372 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Среднее значение |
13230 |
Среднее значение |
19002 |
|
|||||
|
|
|
|
|
|||||
Стандартное отклоне- |
3187 |
Стандартное отклоне- |
2307 |
|
|||||
|
ние |
|
|
|
ние |
|
|
||
|
|
|
|
|
|
|
|
||
|
|
|
|
|
|
|
|||
Из таблицы следует, что среднее значение параметра (Ed/ )1/2 |
для |
||||||||
незаполненных УНТ составляет 13230 |
3187 м/c, в то время как для |
||||||||
«стручков» (Ed/ )1/2 = 19002 2307 м/c. Эти данные показывают, что в результате заполнения УНТ молекулами фуллеренов с учетом изменения плотности среднее значение модуля упругости по отношению к изгибу Ed повышается от 240 105 до 650 156 ГПа. Обратим также
211
внимание на тенденцию к снижению модуля упругости заполненных жгутов с ростом их диаметра, которая, по-видимому, обусловлена эффектом скольжения нанотрубок, составляющих жгут, относительно друг друга при изгибе. Такая тенденция отмечалась также и в [73], где изучали жгуты однослойных нанотрубок.
Упругие механические свойства УНТ являются объектом и многочисленных детальных расчетов. Можно выделить два общих подхода к расчету механических характеристик нанотрубок. Первый – основан на современных методах квантовой химии и молекулярной динамики, учитывающих реальное расположение атомов углерода на цилиндрической поверхности нанотрубки и потенциал взаимодействия этих атомов. Объектом исследования при этом, как правило, является сегмент трубки, состоящий из двух атомов. В этом случае из первых принципов вычисляется изменение потенциальной энергии системы атомов углерода, регулярным образом размещенных на поверхности цилиндра, при различных деформациях. Такая процедура представляется наиболее последовательной, но в силу большой неопределенности формы потенциала взаимодействия между атомами, а также значительного числа атомов, включаемых в рассмотрение, получаемые результаты носят скорее качественный характер. При этом результаты расчетов, выполненных разными авторами на основании различных модельных представлений, могут заметно различаться. Однако в рамках данного подхода имеется возможность установления зависимости механических характеристик УНТ от таких их структурных особенностей, как хиральность и наличие дефектов. В этом случае с успехом применяются полуэмпирические [74] и эмпирические модели [75, 77], а также метод молекулярной динамики [74–82].
Другой подход к анализу механических свойств УНТ основан на использовании современных представлений механики сплошных сред и теории упругости [83, 84]. В этом случае поверхность нанотрубки, составленная из атомов углерода, заменяется сплошной оболочкой, механические характеристики которой определяются исходя из соответствия результатам измерений либо результатам последовательных расчетов.
Такой подход, как всякий эмпирический метод, не является последовательным, однако в настоящее время он является, пожалуй, единственным средством рассмотрения сложных задач, относящихся, например, к установлению механических характеристик многослойных УНТ [11].
212
Одним из первых примеров детального расчета упругих характеристик однослойных УНТ может служить работа [75], в основе которой лежит эмпирическая модель силовых постоянных. При этом сами силовые постоянные, определяющие потенциал взаимодействия атомов, определялись эмпирически из условия согласования расчета и экспериментальных значений упругих констант и фононных частот. Результаты расчета параметров однослойных УНТ различной структуры [75] приведены в табл. 5.5.
Т а б л и ц а 5.5
(n, m) |
|
|
R, нм |
С11 |
С33 |
B, ТПа |
E, ТПа |
M, ТПа |
|
(5.5) |
|
|
0.34 |
0.397 |
1.054 |
0.191 |
0.971 |
0.436 |
0.280 |
|
|
|
|
|
|
|
|
|
|
(6.4) |
|
|
0.34 |
0.397 |
1.054 |
0.191 |
0.972 |
0.437 |
0.280 |
(7.3) |
|
|
0.35 |
0.397 |
1.055 |
0.190 |
0.973 |
0.454 |
0.280 |
(8.2) |
|
|
0.36 |
0.397 |
1.057 |
0.190 |
0.974 |
0.452 |
0.280 |
(9.1) |
|
|
0.37 |
0.396 |
1.058 |
0.191 |
0.974 |
0.465 |
0.280 |
(10.0) |
|
|
0.39 |
0.396 |
1.058 |
0.190 |
0.975 |
0.451 |
0.280 |
|
|
|
|
|
|
|
|
|
|
(10.10) |
|
0.68 |
0.398 |
1.054 |
0.191 |
0.972 |
0.457 |
0.278 |
|
|
|
|
|
|
|
|
|
|
|
(50.50) |
|
3.39 |
0.399 |
1.054 |
0.192 |
0.972 |
0.458 |
0.277 |
|
|
|
|
|
|
|
|
|
|
|
(100.100) |
|
6.78 |
0.399 |
1.054 |
0.192 |
0.972 |
0.462 |
0.277 |
|
|
|
|
|
|
|
|
|
|
|
(200.200) |
|
13.5 |
0.399 |
1.054 |
0.192 |
0.972 |
0.478 |
0.277 |
|
|
|
|
|
|
|
|
|
||
Графит, вдоль гек- |
|
|
|
|
|
|
|||
сагональной |
плос- |
1.06 |
|
0.0083 |
1.02 |
0.44 |
0.16 |
||
кости |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Графит, |
|
перпенди- |
|
|
|
|
|
|
|
кулярно |
|
|
гексаго- |
|
0.036 |
0.0083 |
0.0365 |
0.004 |
0.012 |
нальной плоскости |
|
|
|
|
|
|
|||
Алмаз, |
вдоль куби- |
1.07 |
1.07 |
0.442 |
1.063 |
0.5758 |
0.10415 |
||
ческой оси |
|
||||||||
|
|
|
|
|
|
|
|||
|
|
|
|
|
|
|
|
|
|
Здесь (n, m) – индексы хиральности нанотрубки; R – ее радиус; E – модуль Юнга; M – модуль сдвиговой деформации (кручения относительно оси нанотрубки); B – объемный модуль упругости, определяемый соотношением
p B |
|
V |
|
|
|
|
V |
|
213 |
|
|
где p – однородное давление, которому подвергается объект, а V / V – относительное изменение объема, вызванное этим воздействием; – коэффициент Пуассона, величина которого определяется как отношение относительного поперечного сжатия
к относительному про-
дольному растяжению ( = |
/ ). Для изотропных материалов связь |
||
между объемным модулем упругости и модулем Юнга имеет вид |
|||
|
|
E |
|
B |
|
|
. |
|
3(1 2ν) |
||
Результаты расчета, представленные в табл. 5.5, показывают, что значения модулей упругости ОУНТ практически не зависят от ее диаметра и индексов хиральности, что, однако, может быть естественным следствием расчетов, выполненных на основании эмпирической модели силовых постоянных. Так как в этом случае механические свойства реальной оболочки УНТ моделируются непрерывной оболочкой, свойства которой слабо зависят от угла между направлением внешнего воздействия и ориентацией графитовой плоскости. Данные расчеты также показывают, что показатели жесткости нанотрубок в продольном и радиальном направлениях имеют близкие значения.
Как уже отмечалось, результаты определения модулей упругости нанотрубок весьма чувствительны к выбору значения толщины стенок цилиндрической оболочки, которой моделируется рассматриваемый объект. Такая чувствительность является причиной значительного расхождения результатов расчетов, полученных при использовании различных предположений о значении этого параметра.
Указанных трудностей, связанных с неопределенностью толщины стенок нанотрубки, возникающих в расчетах, использующих модели непрерывных оболочек, удается избежать, применяя альтернативное (5.3.1) определение модуля Юнга [85]:
E |
|
1 |
|
2W |
|
, |
(5.3.5) |
V0 |
2 |
|
|||||
|
|
|
|
||||
|
|
|
|
0 |
|
||
|
|
|
|
|
|
|
|
где V0 – равновесный объем, |
а W – энергия растяжения. Результаты |
||||||
расчета параметров однослойных УНТ различной хиральности и диаметра [85], выполненного методом молекулярных орбиталей в приближении сильной связи [85], приведены в табл. 5.6.
214
|
|
|
Т а б л и ц а 5.6 |
|
|
|
|
Индексы хиральности |
Диаметр, нм |
E, ТПа |
|
(n, m) |
|
||
|
|
|
|
|
|
|
|
(10.0) |
0.791 |
1.22 |
0.275 |
(6.6) |
0.820 |
1.22 |
0.247 |
|
|
|
|
(10.5) |
1.034 |
1.25 |
0.265 |
|
|
|
|
(10.10) |
1.360 |
1.24 |
0.256 |
|
|
|
|
(20.0) |
1.571 |
1.26 |
0.270 |
|
|
|
|
(15.15) |
2.034 |
1.25 |
0.256 |
Сопоставление данных, приведенных в табл. 5.5 и 5.6, показывает, что отличие составляет около 30 %. Эта величина и характеризует возможности современных подходов к моделированию механических свойств УНТ. Отметим также, что, несмотря на значительные расхождения в определении абсолютных значений модуля Юнга, качественный вывод о слабой чувствительности этого параметра к величине диаметра нанотрубки и ее индексам хиральности прослеживается в обоих расчетах. Заметим однако, что условием применимости данных моделей является выполнение неравенства d / L 0.1 (здесь d – диаметр, а L – длина нанотрубки) [86], так как согласно оценкам только в этом случае можно не учитывать влияние конечного продольного размера нанотрубки на ее упругие свойства.
Как уже отмечалось, рассмотренные подходы не позволяют достаточно достоверно установить зависимости модулей упругости от хиральности, диаметра и длины нанотрубки. Однако уникальные механические характеристики (большая механическая прочность в сочетании с ажурным строением и необыкновенной легкостью), предопределяющие широкое использование УНТ в композиционных материалах, продолжают стимулировать все более детальные исследования их упругих свойств. В работе [87] проводилось теоретическое изучение зависимостей модулей Юнга и кручения тонких ОУНТ типа zigzag и armchair с использованием их квантовой модели от диаметра d
идлины L нанотрубок, а также влияния открытых концов. В результате было установлено, что модуль кручения увеличивается с ростом диаметра и уменьшается с ростом длины, в то время как модуль Юнга с увеличением длины нанотрубки увеличивается.
Результаты расчетов, выполненных в [87], приведены на рис. 5.28
и5.29.
215
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
, ТПа |
|
ТПа |
|
|
|
|
|
|
|
|
E, |
|
|
|
|
|
|
E |
|
|
|
|
|
|
|
|
|
L, A |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
N |
||||
|
|
а |
|
|
|
|
б |
|
Рис. 5.28. Изменение модуля Юнга:
а – с увеличением числа атомов УНТ (6, 0), б – с увеличением длины УНТ [87]. Кривые 1 и 2 соответствуют двум разным моделям
L, A |
d, A |
а |
б |
Рис. 5.29. Изменение модуля кручения УНТ с увеличением длины (а); изменение модуля кручения с увеличением диаметра УНТ при неизменной длине (б); для zigzag УНТ L = 15.6 A, для armchair УНТ L = 23.31 A [87]; f E – расчет из энергии кручения; fY – расчет через модуль Юнга
Изменение параметров упругости трубок типа zigzag и armchair с увеличением их диаметра отражено в табл. 5.7. Из расчетов следует, что модули Юнга и кручения нанотрубок стремятся к насыщению при отношении диаметра к длине ~ 0.3. Можно также считать, что коэффициент Пуассона
при изменении диаметра трубки остается постоянным. При этом согласно расчетам для трубок zigzag ν 0.43 , а для трубок armchair ν 0.45 .
216
Т а б л и ц а 5.7
|
УНТ типа zigzag |
|
|
УНТ типа armchair |
|
||||
d, нм |
N |
(n, m) |
E, ТПа |
|
d, нм |
N |
(n, m) |
E, ТПа |
|
0.483 |
96 |
(6, 0) |
0.671 |
0.445 |
0.418 |
120 |
(3, 3) |
0.683 |
0.439 |
0.560 |
112 |
(7, 0) |
0.692 |
0.433 |
0.551 |
160 |
(4, 4) |
0.721 |
0.444 |
0.637 |
128 |
(8, 0) |
0.709 |
0.435 |
0.685 |
200 |
(5, 5) |
0.740 |
0.447 |
0.714 |
144 |
(9, 0) |
0.722 |
0.433 |
0.819 |
240 |
(6, 6) |
0.751 |
0.448 |
0.792 |
160 |
(10, 0) |
0.730 |
0.431 |
0.954 |
280 |
(7, 7) |
0.758 |
0.449 |
0.869 |
176 |
(11, 0) |
0.733 |
0.429 |
|
|
|
|
|
Наряду с модулями упругости важной характеристикой нанотрубок является предельная прочность на разрыв, определяемая как растягивающее усилие (в расчете на единицу площади), при превышении которого происходят необратимые изменения в структуре нанотрубок.
Из определения предельной прочности на разрыв следует, что величина этого параметра не должна зависеть от числа нанотрубок, используемых при измерениях. Таким образом, при экспериментальном измерении этого параметра можно использовать жгуты нанотрубок и тем самым избежать ряда проблем, возникающих при работе с весьма миниатюрными образцами индивидуальных УНТ. Данный подход был успешно реализован в работе [88]. Для приложения к жгуту механической нагрузки авторы использовали наконечник атомного силового микроскопа. Максимальное значение деформации, полученное в результате обработки измерений, выполненных на четырех образцах, составило 5.8 0.9 %. Это значение соответствует величине прочности на разрыв жгута нанотрубок 45 7 ГПа, что примерно в 20 раз превышает соответствующее значение для высокопрочных сталей. Следует, однако, отметить, что многослойные УНТ, выращенные методом CVD, обладают значительно более низкой прочность на разрыв. Согласно [69] прочность на разрыв таких трубок составляет всего лишь 1.72
0.64 ГПа. Таким образом, прочность УНТ, полученных методом химического осаждения в парах, примерно на порядок ниже, чем при синтезе традиционным электродуговым методом, что связано с появлением большого числа дефектов при использовании метода CVD.
Подробные исследования механических характеристик УНТ в широком интервале деформаций были проведены в работах [76, 89]. В рамках метода молекулярной динамики авторы рассчитали зависимости напряжение–деформация в упругой, упруго-пластической и пластической областях для многослойных нанотрубок с различным числом
217
Механическое напряжение, ГПа
расчет [88] расчет [76]
Одноосная деформация растяжения
Рис. 5.30. Зависимость напряжения от одноосной деформации для ОУНТ (12, 12)
Механическое напряжение, ГПа Механическое напряжение, ГПа
150
125
100
75
50
25
00.00 0.05 0.10 0.15 0.20 0.25
Одноосная деформация растяжения
а
160
120
80
60
00 |
|
|
|
|
|
0.00 |
0.05 |
0.10 |
0.15 |
0.20 |
0.25 |
Одноосная деформация растяжения
в
Механическое напряжение, ГПа Механическое напряжение, ГПа
180
140
100
60
20
0.00 |
0.05 |
0.10 |
0.15 |
0.20 |
0.25 |
Одноосная деформация растяжения
б
160
120
80
60
00 |
|
|
|
|
|
0.00 |
0.05 |
0.10 |
0.15 |
0.20 |
0.25 |
Одноосная деформация растяжения
г
Рис. 5.31. Зависимость механического напряжения от одноосной деформации для УНТ с L/d = 4.5:
а – ОУНТ (10, 10), б – двухслойная УНТ (5 ,5)&(10,10), в – трехслойная УНТ
(5 ,5)&(10,10)&(15, 15), г – трехслойная УНТ (5 ,5)&(10,10)&(15, 15)&(20,20) [76]
218
слоев. Результаты расчетов приведены на рис. 5.30–5.32, а также в табл. 5.8.
Дополнительные сведения о строении, свойствах и возможных применениях углеродных нанотрубок можно найти в работах [90–93].
Механическое напряжение, ГПа Механическое напряжение, ГПа
160
120
80
60
00 |
|
|
|
|
|
0.00 |
0.05 |
0.10 |
0.15 |
0.20 |
0.25 |
Одноосная деформация растяжения
а
160
120
80
60
00 |
|
|
|
|
|
0.00 |
0.05 |
0.10 |
0.15 |
0.20 |
0.25 |
Одноосная деформация растяжения
в
Механическое напряжение, ГПа Механическое напряжение, ГПа
180
140
100
60
20
0.00 |
0.05 |
0.10 |
0.15 |
0.20 |
0.25 |
Одноосная деформация растяжения
б
160
120
80
60
00 |
|
|
|
|
|
0.00 |
0.05 |
0.10 |
0.15 |
0.20 |
0.25 |
Одноосная деформация растяжения
г
Рис. 5.32. Зависимость механического напряжения от одноосной деформации для УНТ с L/d = 9.1:
а – ОУНТ (10, 10), б – двухслойная УНТ (5 |
,5)&(10,10), в – трехслойная УНТ |
(5 ,5)&(10,10)&(15, 15), г – трехслойная УНТ |
(5 ,5)&(10,10)&(15, 15)&(20,20) [76] |
219
|
|
|
|
|
|
|
Т а б л и ц а 5.8 |
||
|
|
|
|
|
|
|
|
|
|
|
|
|
Упругая |
Упруго-пласти- |
Предельные па- |
||||
Тип на- |
|
Модуль |
область |
ческая область |
раметры |
||||
L/D |
|
|
|
|
|
|
|||
нотрубки |
Юнга E |
напря- |
дефор- |
напря- |
дефор- |
напря- |
дефор- |
||
|
|||||||||
|
|
|
жение, |
жение, |
жение, |
||||
|
|
|
мация |
мация |
мация |
||||
|
|
|
ГПа |
ГПа |
ГПа |
||||
|
|
|
|
|
|
||||
Одно- |
4.5 |
1.043 |
61.03 |
0.0585 |
136.9 |
0.231 |
140.4 |
0.280 |
|
слойная |
|||||||||
|
|
|
|
|
|
|
|
||
трубка |
|
|
|
|
|
|
|
|
|
9.1 |
1.031 |
62.71 |
0.0594 |
142.1 |
0.236 |
148.5 |
0.279 |
||
(10, 10) |
|||||||||
|
|
|
|
|
|
|
|
||
Двух- |
4.5 |
1.161 |
72.31 |
0.0627 |
161.6 |
0.247 |
162.4 |
0.279 |
|
слойная |
|||||||||
|
|
|
|
|
|
|
|
||
трубка |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
(5, 5), |
9.1 |
1.175 |
72.87 |
0.0621 |
163.3 |
0.242 |
168.4 |
0.281 |
|
(10, 10) |
|
|
|
|
|
|
|
|
|
Трехслой- |
4.5 |
1.000 |
60.68 |
0.0605 |
143.0 |
0.238 |
143.4 |
0.281 |
|
ная трубка |
|||||||||
(5, 5), |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
(10, 10) |
9.1 |
0.972 |
56.45 |
0.0611 |
138.1 |
0.246 |
141.4 |
0.282 |
|
(15, 15) |
|||||||||
|
|
|
|
|
|
|
|
||
Четырех- |
4.5 |
0.932 |
60.75 |
0.0654 |
134.3 |
0.235 |
138.2 |
0.281 |
|
слойная |
|||||||||
|
|
|
|
|
|
|
|
||
трубка |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
(5, 5), |
|
|
|
|
|
|
|
|
|
(10, 10), |
9.1 |
0.872 |
57.84 |
0.0633 |
127.8 |
0.241 |
132.7 |
0.280 |
|
(15, 15) |
|||||||||
|
|
|
|
|
|
|
|
||
(20, 20) |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
5.4.Фононный спектр углеродных нанотрубок
Исследование фононных мод УНТ проводилось в ряде работ [49, 90, 94, 95]. Согласно [90, 96] частоты колебательных мод ОУНТ могут быть вычислены из спектра двумерного графенового слоя с помощью уравнения
1D (k) 2D kK2 K1 , |
0, 1, 2, ..., N 1 , |
где 1D – частота колебательной моды для одномерной нанотрубки, 2D – частота для двумерного графенового слоя, k – волновой вектор
220
