Читать далее - Нанотехнологическое общество России

advertisement
О ПОВЕРХНОСТНОЙ ЭНЕР ГИИ И ПОВЕРХНОСТНОМ ДАВЛЕНИИ
В НАНОКРИСТАЛЛАХ АЛМ АЗА И КРЕМНИЯ
М.Н. М АГО МЕД ОВ
Институт проблем геотермии Дагестанского научного центра РАН, Махачкала, Россия
Поверхностные свойства играют основную роль как в процессе возникновения нанокристалла, так и
при росте кристалла вообще. В связи с этим в данной работе без каких-либо подгоночных параметров
(исходя только из потенциала межатомного взаимодействия) произведен аналитический расчет
поверхностных свойств (удельной поверхностной энергии, поверхностного натяжения и поверхностного
давления) нанокристаллов алмаза и кремния. Изучена эволюция этих поверхностных свойств при
изменении, как размера, так и формы поверхности нанокристалла при различных температурах.
Результаты показали, что поверхностное давление для нанокристалла: Psf(N)  N –1/3, оказалось намного
меньшим, чем давление Лапласа: Pls(N)  N – 1/3, для такого же нанокристалла при данных значениях
температуры T, плотности N/V и числе атомов N. На примере нанокристаллов алмаза и кремния показано,
что с ростом температуры от 20 K до плавления изотерма Psf(N) понижается намного заметнее, чем
изотерма Pls(N), а при высоких температурах изотерма Psf(N), изменив форму зависимости от N, заходит в
область растяжения малых нанокристаллов. Показано, что параметр решетки может, как сжиматься, так
и растягиваться при уменьшении размера нанокристалла, либо при деформации его формы.
Как и в наших предыдущих работах [1-3] положим, что нанокристалл со свободной поверхностью
имеет вид прямоугольного параллелепипеда с квадратным основанием, ограненный гранями (100).
Величина f = Nps / Npo – это параметр формы, который определяется отношением числа атомов на
боковом ребре Nps к числу атомов на ребре основания Npo . Число атомов в нанокристалле: N = f Npo3 /  ,
изменяется в пределах: 23 /   N  , где  =  / (6 kp), kp – коэффициент упаковки структуры.
Как было показано в работах [1-3], зависимость нормированного значения среднего
координационного числа kn* от числа атомов N и параметра формы f определяется выражением:
kn*(N, f) = kn(N, f) / kn(N=) = 1 – Zs(f) (2 / N ) 1/3 ,
(1)
Zs(f) = (1 + 2 f) / (3 f 2/3) .
(2)
Функция формы Zs(f) достигает минимума равного единице при f = 1, т.е. для формы куба: Zs(f = 1) = 1.
Для пластинчатых (f < 1) или стержневидных (f > 1) форм значение Zs(f) больше единицы. Поэтому
функция kn*(f) при любом N имеет максимум при f = 1, т.е. для наиболее термодинамически устойчивой
(среди всевозможных форм прямоугольного параллелепипеда) – кубической формы нанокристалла.
Объем V и площадь поверхности  для прямоугольного параллелепипеда равны:
V = Npo3 f [c(N, f)]3 = N  [c(N, f)]3 ,
 = 6 [c(N, f)]2 s ( N )2/3 Zs(f) ,
(3)
где c(N, f) – среднее (по всему объему нанокристалла) расстояние между центрами ближайших атомов, s
– коэффициент, учитывающий плотность упаковки атомов на грани нанокристалла: s  2/3. Легко
видеть, что объем нанокристалла V в выражении (3) зависит от формы системы (т.е. от величины f)
только через зависимость c(N, f), т.е. через расстояние между центрами ближайших атомов.
Пусть взаимодействие атомов описывается парным потенциалом Ми-Леннарда-Джонса:
(r) = [ D / ( b – a ) ][ a ( ro / r) b – b ( ro / r) a ] ,
(4)
где D и ro – глубина и координата минимума потенциальной ямы, b и a – параметры, определяющие
жесткость и дальнодействие парного потенциала: b > a.
Как показано в [1, 2], удельную поверхностную энергию , ее производную по температуре, и
поверхностное давление Psf можно найти из выражений:
(T, c, N, f) = – [kn(N=) D / (12 s ro2)] R2 Le(N, f) ,
(5)
(  /  T)c, N, f = – {  kb / [2 s (b + 2) ro2 kn*(N, f)]} R2 CE(e/T) ,
(6)
Psf(T, c, N, f) = Pls ( 1 + p ) ,
p = (1/2)( ln  /  ln c)T, N, , f .
(7)
Здесь R = ro/c, kb – постоянная Больцмана, e – температура Эйнштейна,  и q первый и второй
параметры Грюнайзена (метод расчета данных функций описан в работе [4]),
p = – 1 – (1/2){ U(R) – 54 [ / (b + 2)](q –  ty)[kb e /D kn(N, f)] Ew(y) } / Le(N, f) ,
Le(N, f) = U(R) + 18 [ / (b + 2)][kb e(N, f) / D kn(N, f)] Ew(y) ,
y = e/T ,
(8)
ty(y) = 1 – {2 y exp(y) / [exp(2 y) – 1]},
Ew(y) = 0.5 + [exp(y) – 1] – 1,
CE(y) = 3 y2 exp(y) / [exp(y) – 1] 2,
U(R) = R [ U(R) /  R] = a b (R b – R a) / (b – a) .
В формуле (7) первый сомножитель – это так называемое «давление Лапласа»:
Pls =  (  /  V)N, , f = (2/3)( / V)  = 4 s Zs(f)  / [c (N )1/3] = 4 s [1 – kn*(N, f)]  / (c ).
(9)
При Psf > 0 поверхностное давление сжимает, а при Psf < 0 – растягивает нанокристалл.
2
Значения параметров межатомного потенциала (4) для алмаза (m = 12.01 a.m.u., kn = 4, kp = 0.3401, 3
= 2.25) из [5], которые были использованы для расчетов, и получаемые из них при R = ro / c = 1 значения
характеристических температур и параметров Грюнайзена равны:
ro = 0.1545 nm, Ds / kb = 42702.72 K = 3.68 eV, b = 3.79, a = 2.21, Aw(R=1)3 = 13.427 K,
e(R=1) = 1060.83 K,
(R=1) = 1414.44 K,
(R=1) = 0.956,
q(R=1)103 = 9.16.
В таблице 1 представлены результаты расчета удельной поверхностной энергии (5) и ее изохорной
производной по температуре (6) для макро-алмаза (N = ) при различных значениях температуры от 20
до Tm = 4300 K и R = 1. Температура 4300 K была определена в [8] как теоретическая температура
плавления алмаза. Сравнение рассчитанных значений  с теоретическими оценками (100)est для алмаза
из работ [6, 7] показало хорошее согласие.
T
(100), (5)
T / e
(100)est – (  /  T)V , (6) – p , (8)
ten(100)
– ten / T
K
10 – 3 J / m2 10 – 3 J/ m2
10 – 6 J / m2 K
10 – 3 J / m2
10 – 6 J / m2 K
20
0.019
6118.2
0
0.986772
80.928
6.49
150
0.141
6118.1
4.57
0.986910
80.084
6513 [6]
81.05
300
0.283
6114.8
41.51
0.988891
67.927
84.62
600
0.566
6094.7
83.33
0.993020
42.541
5710 [7]
40.75
1000
0.943
6057.8
97.90
0.995668
26.240
14.95
2000
1.885
5955.3
104.95
0.998104
11.292
6.18
3000
2.828
5849.5
106.32
0.999125
5.117
3.69
4300
4.054
5710.8
106.89
0.999943
0.326
Формулу (7) можно переписать в виде: Psf = ten (  /  V)N, , f , где ten – поверхностное натяжение,
которое определяется уравнением Шаттлвоса [9–11]: ten =  (1 + p). Рассчитанные величины ten и ее
изменение с температурой для грани (100) приведены в таблице 1.
1,00
14
20 K
Diamond
0,95
12
20 K
0,90
(d/dT)R, f*
*
10
8
4300 K
0,85
6
0,80
4
4300 K
0,75
2
0
Рис. 1.
1
2
3
lg(N)
4
5
6
На рис. 1 показаны рассчитанные изотермо-изоморфные (f = 1: форма куба) размерные зависимости
нормированных (на значение для макро-кристалла) величин: удельной поверхностной энергии (5)
(сплошные линии, левая шкала) и изохорно-изоморфной производной удельной поверхностной энергии
по температуре (6) (пунктир, правая шкала) для нано-алмаза при R = 1 и температурах равных 20 K
(линии с кубиками) и 4300 K (линии с кружками). Символы на изоморфах указывают положение
разрешенных (при данном f) значений N в нанокристалле: N = INT[f Npo3 / ] при Npo = 2, 3, 4, … .
Функция INT[X] округляет X до целого. Число атомов в минимальном (Npo = 2) нано-алмазе равно: Nmin =
INT[23 / 1.54] = 5. Любое отклонение формы нанокристалла от кубической (стержневидная: f > 1, или
пластинчатая: f < 1) приводит к изменению величины Nmin и положения разрешенных значений N в
изоморфных (при данном f) нанокристаллах. Размерные изотермо-изоморфные зависимости для
«некубичных» нанокристаллов лежат на рис. 1 для функции * ниже, а для функции (  /  T)c, N, f* выше,
чем соответствующие зависимости для «кубических» нанокристаллов.
3
Diamond
Psf [kbar]
K
800
43
600
30
20
Pls [kbar]
1000
00
25
20
K
15
400
10
20 K
200
5
4300 K
0
0
Рис. 2.
1
2
lg(N)
3
4
5
На рис. 2 представлены рассчитанные изотермо-изоморфные (f = 1) размерные зависимости
давления Лапласа (пунктир с кубиками, левая шкала) и поверхностного давления (сплошные линии с
треугольниками, правая шкала) для нано-алмаза при R = 1 и T = 20  4300 K. Из рис. 2 видно, что
величина Psf из (7) намного меньше, чем значение Pls рассчитанное из (9). Более того, если с ростом
температуры от 20 K до 4300 K изотерма Pls(f=1, N) сравнительно мало понижается, то изотерма
Psf(f=1,N) понижается довольно заметно и даже, изменив форму зависимости от N, заходит в область
растяжения малых нанокристаллов, где Psf < 0. Изотерма Psf(3000 K, f=1, N) при Npo = 3 (т.е. при N = 17)
имеет максимум: Psf(3000 K, f=1, 17)max = 384.6 bar. Изотерма Psf(4300 K, f=1, N) между Npo = 15 и 16 (т.е.
при N = 2192  2660) переходит из области сжатия в область растяжения, а при Npo = 30 (т.е. при N =
17537) она имеет максимум: Psf(4300 K, f=1, 17537)max = 1.9 bar.
Значения параметров межатомного потенциала (4) для кремния (m = 28.09 a.m.u., kn = 4, kp = 0.3401,
3 = 2.25) из [5], которые были использованы для расчетов, и получаемые из них при R = ro / c = 1
значения характеристических температур и параметров Грюнайзена равны:
ro = 0.2351 nm,
Ds / kb = 26921.28 K = 2.32 eV,
b = 4,
a = 2.48,
Aw(R=1)3 = 3.186 K,
e(R=1) = 411.6 K,
(R=1) = 549.06 K,
(R=1) = 0.994,
q(R=1)103 = 5.802.
В таблице 2 представлены результаты расчета удельной поверхностной энергии (5) и ее изохорной
производной температуре (6) для макро-кремния (N = ) при различных значениях температуры от 20 до
Tm = 1685 K и R = 1. Сравнение рассчитанных значений  с теоретическими оценками (100)est из [6, 7,
12] показало хорошее согласие с наиболее свежими данными из [7].
T
(100), (5)
– (  /  T)V
– p , (8)
T / e
(100)est
ten(100)
– ten / T
K
10 – 3 J / m2
10 – 3 J / m2 10 – 6 J / m2K
10 – 3 J / m2 10 – 6 J / m2K
20
0.049
1671.69
0
0. 991496
14.216
2165 [6]
7.84
70
0.170
1671.64
4.517
0. 99173
13.824
42.69
150
0.364
1670.37
25.743
0. 993769
10.409
1340 
29.27
300
0.729
1665.18
39.897
0. 996386
6.018
1490 [7]
10.17
600
1.457
1652.25
44.783
0. 998204
2.968
3.57
1000
2.428
1634.05
45.916
0. 999057
1.541
1140 [12]
1.89
1300
3.157
1620.23
46.181
0. 999399
0.974
1.29
1685
4.092
1602.42
46.338
0. 999703
0.476
В работах [10, 11] был подвергнут сомнению сам факт существования поверхностного давления для
нанокристаллов. Причем, в качестве экспериментального подтверждения «отсутствия давления Лапласа»
приводились экспериментальные результаты, обнаружившие, как увеличение параметра решетки
нанокристаллов при уменьшении их размера, так и независимость межатомного расстояния от числа
4
атомов N в нанокристалле. Как следует из полученных результатов, поверхностное давление для
ограненного нанокристалла отлично от нуля, но оно много меньше того поверхностного давления,
которое получается по формуле Лапласа (9). Из полученных результатов видно, что данный эффект
обусловлен поведением функции p, которая для жидкой фазы равна нулю, так как: (  /  )T, N = 0.
Исходя из наших результатов, размерное изменение расстояния между центрами ближайших атомов
в нанокристалле (или среднего по нанокристаллу параметра решетки) можно понять из следующих
очевидных соотношений:
V / V = 3 c / c = – P / BT = – Psf / BT = – Pls ( 1 + p ) / BT ,
(10)
где BT = – V ( P /  V)T – изотермический модуль упругости, который вне области фазовых переходов
первого рода (плавление, или полиморфные переходы) всегда больше нуля.
В работе [13] отмечалось, что для металлов выполняется следующее корреляционное соотношение:
 / BT  0.05 nm. Для алмаза при T = 0 K получено: BT=0K = 4430 kbar [14]. Используя значение из табл. 1:
(0 K) = 6118.210 – 3 J / m2, получим:
(0 K) / BT=0K  0.014 nm, или (Cdia) = (0 K) / (ro BT=0K)  0.0906.
(11)
Для кремния: BT=0K = 977 kbar [14], а из табл. 2: (0 K) = 1617.710 – 3 J / m2. Это дает:
(0 K) / BT=0K  0.017 nm, или (Si) = (0 K) / (ro BT=0K)  0.0728.
(12)
Допустим, что полученные для макрокристалла корреляционные соотношения (11) и (12)
выполняются и для нанокристалов. Тогда из (9)-(12) получим:
c/c = – Pls (1 + p) / (3 BT )  – (4/3) [(0 K) / BT=0K] (s / ) [1 – kn*(N, f)](1 + p) / c 
 – (4/3)  (s / ) [1 – kn*(N, f)](1 + p) / (c / ro) .
(13)
Из (13) видно, что для макрокристалла (т.е. при N  , и когда в соответствии с (1): kn* 1) имеем:
c/c  0. Что касается нанокристалла, то здесь величина c/c определяется не только «размерным
аргументом» kn*, но и величиной p < 0, а знак разности c определяется знаком функции: 1 + p.
На рис. 3 и 4 показаны изотермо-изоморфные зависимости параметра решетки c/c (в %) от размера,
оцененные по формуле (13) для нанокристаллов алмаза и кремния при R = 1 и двух значений параметра
формы: f = 1 (сплошные линии и сплошные фигуры) и f = 5 (пунктирные линии и открытые фигуры).
Расчеты показали, что с ростом температуры функция c меняет знак для малых нанокристаллов. При
низких температурах величина c/c отрицательна и ее абсолютная величина тем больше, чем меньше
значение размерного аргумента kn*. С ростом температуры функция c сначала зануляется, а потом
переходит в положительную область, где величина c/c тем больше, чем меньше значение kn*: т.е. чем
меньше число атомов N вдоль изоморфы (т.е. при f = const), либо чем больше величина параметра формы
f отклонена от единицы вдоль изомерной зависимости (т.е. при N = const).
c/c x 100
Diamond
4300
K
0,00
-0,02
-0,04
f=1
f=5
K
00
0
1
f=1
f=5
20
K
-0,06
-0,08
lg(N)
Рис. 3.
-0,10
1
2
3
4
На рис. 3 две верхние спадающие линии с кружками – изотермы 4300 K, две средние
возрастающие кривые с треугольниками – изотермы 1000 K, две нижние возрастающие кривые с
квадратами – изотермы 20 K. Верхние линии с открытыми и сплошными кружками (изотермы 4300 K) –
изоморфы f = 1 и f = 5 в данном масштабе сливаются.
5
c/c x 100
0,00
168
5K
K
1000
f=1
Silicon
f=5
-0,01
Рис. 4.
20
K
f=5
15
0K
30
0
K
f=1
lg(N)
-0,02
1
2
3
4
На рис. 4 возрастающие кривые: с ромбиками – изотермы 1000 K, с треугольниками – изотермы 300
K, со звездочками – изотермы 150 K, с кружками – изотермы 20 K. с треугольниками – изотермы 1000 K,
две нижние возрастающие кривые с квадратами – изотермы 20 K. Верхняя линия с квадратами –
изотермы 1685 K: изоморфы f = 1 и f = 5 в данном масштабе сливаются.
Из рис. 3 и 4 видно, что наиболее заметное размерное сжатие решетки нанокристаллов алмаза и
кремния получаются при относительно низких температурах: при T < 1000 K или при T / e(R=1) < 1.
Работа выполнена при финансовой поддержке Программы Президиума РАН (проект № 12.1.19) и
РФФИ (грант № 09 – 08 – 96508-р-юг-а).
1.
2.
3.
4.
5.
6.
7.
8.
9.
10.
11.
12.
13.
14.
СП И С О К ЛИ Т ЕР А ТУ РЫ
Магомедов М.Н. О зависимости поверхностной энергии от размера и формы нанокристалла // Физика
Твердого Тела. 2004. Т. 46, № 5. С. 924 – 937.
Магомедов М.Н. Об изменении «поверхностного» давления в нанокристалле с температурой //
Теплофизика Высоких Температур. 2005. Т. 43, № 6. С. 870 – 879.
Магомедов М.Н. О зависимости температуры фазового перехода кристалл-жидкость от размера и
формы простого нанокристалла // Теплофизика Высоких Температур. 2009. Т.47, № 1. С.49–60.
Магомедов М.Н. Об определении температуры Дебая из экспериментальных данных // Физика
Твердого Тела. 2003. Т. 45, № 1. С. 33 – 36.
Магомедов М.Н. О природе ковалентной связи в кристаллах подгруппы углерода // Журнал
Неорганической Химии. 2004. Т. 49, № 12. C. 2057 – 2067.
Oshcherin B.N. On Surface Energies of ANB8–N Semiconducting Compounds // Physica Status Solidi (a).
1976. V. 34, № 2. P. K181 – K186.
Stekolnikov A.A., Bechstedt F. Shape of Free and Constrained Group-IV Crystallites: Influence of Surface
Energies // Phys. Rev. B. 2005. V. 72, № 12. P. 125326.
Шульман Л.А. Температуры плавления алмаза и BNсф и критерий Линдеманна // Сверхтвердые
материалы. 1993, № 4. С. 58 – 59.
Shuttleworth R. The Surface Tension of Solids // Proceedings of the Physical Society. Section A. (London).
1950. V. 63. Part 4, № 364A. P. 444 – 457.
Зубов В.И. Некоторые размерные эффекты и свойства ультрадисперсных систем // Журнал
Всесоюзного Химического Общества им. Д.И. Менделеева. 1991. Т. 36, № 2. С. 133 – 137.
Нагаев Э.Л. Малые металлические частицы // Успехи Физ. Наук. 1992. Т. 162, № 9. С. 49 – 124.
Jiang Q., Luand H.M., Zhao M. Modeling of Surface Energies of Elemental Crystals // Journal of Physics:
Condensed Matter. 2004. V. 16, № 4. P. 521 – 530.
Комник Ю.Ф. Физика металлических пленок. М.: Атомиздат, 1979. 264 с.
Стишов С.М. Энергия, сжимаемость и ковалентность в подгруппе углерода // Письма в Журнал
Экспериментальной и Теоретической Физики. 2000. Т. 71, № 1. С. 25 – 27.
Download