первое слагаемое – энергия ненапряженного изогнутого бруса, два других суть вклад энергии деформации, обусловленный действием механического напряжения Т0 + Т, где Т0 – остаточное напряжение, которое могло возникнуть, например, в процессе роста трубки, Т – напряжение, связанное с удлинением трубки под действием кулоновских сил со стороны управляющего электрода:
T = |
ES |
∫0L z′2 dx. |
(8.5) |
|
2L |
||||
|
|
|
Приближенно считая управляющий электрод бесконечной плоскостью, расположенной на расстоянии R от нанотрубки, для емкости на единицу ее длины с(z) (см. рис. 74) имеем
c(z) = |
|
|
1 |
|
≈ |
1 |
+ |
z(x) |
(8.6) |
|
|
|
|
|
|
|
, |
||||
2 ln |
2 |
(R − z ) |
r |
2ln [2R r] |
2 ln2 [2R r] |
|||||
|
|
|
|
|
|
|
|
|
|
|
где разложение в ряд Тейлора верно при условии z << R. В этом случае силами Ван-дер-Ваальса между нанотрубкой и подложкой можно пренебречь. Выражение для электростатической энергии системы довольно громоздко. Мы ограничимся здесь рассмотрением случая, когда емкости потенциальных барьеров CL и CR, соединяющих нанотрубку с левым и правым электродами, пренебрежимо малы (более общая ситуация рассмотрена в [64]).). Полагая также Т0 = 0, получим для WQ:
|
Q |
( |
|
) |
|
(ne)2 |
G |
|
|
|
|
2CG [z] |
|
||||
|
W z |
|
x |
|
= |
|
− neV |
≈ |
|
|
|
|
|
|
|
|
(8.7) |
≈ |
(ne)2 ln (2R r ) |
− |
(ne2 )2 ∫0L z (x)dx − neVG . |
|||||
|
L |
|
|
|
|
L R |
|
|
где VG – потенциал управляющего электрода (потенциал V источника (электрода L) в рассматриваемом приближении не входит в выражение для электростатической энергии, а потенциал стока (электрода R) выбран за начало отсчета), пе – избыточный (квантованный) заряд нанотрубки. В предположении, что заряд распределен по поверхности нанотрубки равномерно, для ее емкости относительно электрода G можно написать:
C |
G |
= |
L c z |
( |
x |
dx. |
(8.8) |
|
|
∫0 |
|
) |
|
126
Варьируя функционал полной энергии системы
W z |
( |
x |
= W z |
( |
x |
+W z |
( |
x |
|
||
n |
|
) |
el |
|
) |
Q |
|
) |
|||
относительно функции z(x) при постоянном числе избыточных электронов на трубке, получим уравнение для формы нанотрубки, [60]:
IEz′′′′−Tz′′ = K0 ≡ |
(ne)2 |
|
z |
|
(8.9) |
|
2 |
+ O |
|
|
, |
||
|
||||||
|
L R |
|
R |
|
|
|
где K0 – кулоновская сила, действующая на единицу длины нанотрубки (при z/R << 1 ее можно приближенно считать постоянной). Предполагая, что напряжение T есть константа (мы определим ее позже из условия самосогласования (8.5)), нетрудно найти решение уравнения (8.9) с понятными граничными условиями (z(0) = z(L) = 0, z′(0) = z′(L) = 0):
|
|
|
sh (ξL) |
|
zn (x) = |
K0 L |
|||
|
|
|
||
2T ξ |
ch (ξL) −1 |
|||
|
|
|||
|
|
|
|
|
ξ =
|
|
x |
2 |
|
|
ch(ξx) −1 − sh (ξx) +ξx − ξ |
|
|
, |
||
|
|
|
|||
|
|
L |
|
||
|
|
|
(8.10) |
||
|
|
|
|
|
|
EIT .
Подставляя (8.10) в (8.5), найдем связь между напряжением Т в нанотрубке и линейной плотностью внешней силы K0. Эту связь легко получить в двух наиболее интересных предельных случаях:
K0 |
L S |
(60480EI |
|
), T << EI L , |
|
|
2 |
6 |
|
2 |
2 |
T = |
(ES 24)1 3 (K0 L)2 3 |
(8.11) |
|||
|
, T >> EI L2 . |
||||
|
|
|
|
|
|
Верхняя строка в (8.11) соответствует малому изгибу трубки, |
|||||
когда z ≤ r; нижняя – |
случаю, когда трубка изогнута сильно, т.е. |
||||
r < z << R, L. Для каждого их двух этих режимов можно определить также смещение нанотрубки в ее средней точке – величину
znmax = zn (L
2):
|
|
|
(ne) |
2 |
2 |
|
|
|
|
|
|
||
|
0,003 |
|
L |
, T << |
EI |
|
|
||||||
|
E r4 |
R |
L2 |
|
|
||||||||
max |
|
|
|
|
|
|
|
||||||
zn |
= |
|
|
(ne)2 3 L2 3 |
|
EI |
|||||||
|
|
|
|
|
|||||||||
|
|
0,24 |
|
|
|
|
|
|
|
, T >> |
|
||
|
|
E |
1 3 |
r |
2 3 |
R |
1 3 |
2 |
|||||
|
|
|
|
|
|
|
|
|
|
L |
|||
|
Er5 R |
|
||
n << |
|
|
|
; |
e |
2 2 |
|||
|
L |
|
(8.12) |
|
|
|
|
|
|
n >> Er5 R .
e2 L2
127
Для одностенной нанотрубки с r = 0,65 нм, Е = 1,25 ТПа, L = = 500 нм, R = 100 нм (назовем ее Е-нанотрубкой) кроссовер –
переход |
от режима малого изгиба к режиму большого изгиба |
(T EI |
L2 ) происходит уже при n ~ 5 ÷ 10. В режиме большого |
изгиба максимальное отклонение формы Е-нанотрубки от прямолинейной составляет znmax = 0,24n2
3 (нм). Заметим, что этот режим
недостижим ни в одном из существующих субмикронных устройств на основе кремния: все они функционируют только в режиме малого изгиба.
8.2.2. Равновесный заряд нанотрубки
Результатов, полученных выше, недостаточно для сравнения с экспериментом – необходимо установить соотношение между количеством избыточных электронов на трубке и напряжением на управляющем электроде VG. Для этого следует минимизировать
равновесную полную энергию системы |
W eq = W z |
n ( |
x |
|
по отно- |
|
|
n |
n |
|
) |
|
|
шению к параметру n. В обоих рассмотренных выше предельных случаях мы можем выписать функцию Wneq в явном виде:
W eq |
≡W |
st |
− δW = (ne)2 |
ln |
2R − neV − |
|
||||||||||
|
n |
|
|
|
|
|
|
|
L |
|
|
r |
G |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
(8.13) |
||
|
|
|
|
|
|
|
|
4 |
|
4 |
2 |
|
|
|
2 |
|
0,0009(ne) |
|
|
(Er R |
|
), |
T << EI |
L ; |
|
||||||||
|
|
|
|
|
||||||||||||
− |
|
|
( |
|
|
) |
|
|
( |
|
|
1 3 |
|
|
|
|
|
|
|
ne |
8 3 |
|
Er2 R4 L |
) |
, T >> EI |
L2 . |
|
||||||
0,08 |
|
|
|
|
|
|
||||||||||
Два первых слагаемых в (8.13) отвечают электростатической энергии неизогнутой трубки, а третье слагаемое есть вклад механических степеней свободы – напряжения, изгиба и изменения емкости CG при изгибе трубки. Как правило, это нелинейное наномеханическое слагаемое мало по сравнению с остальными: в Е- нанотрубке оно достигает того же порядка величины, что и вклад от Wst, лишь при п ≈ 3000, когда приближение (8.6) уже неверно. Отрицательный знак наномеханического слагаемого легко объясним: трубка подвижна, и минимум энергии системы при изменении напряжения VG достигается не только за счет изменения заряда трубки, но и посредством изменения ее формы.
128
Из (8.13) для величины п, минимизирующей энергию, получим
|
|
VG L |
|
1 |
|
|
(8.14) |
n = |
Int |
|
+ |
|
+ δn |
, |
|
|
|
||||||
|
|
2e ln (2R r ) |
|
2 |
|
|
|
где символ Int (...) |
|
|
|
|
|
||
предполагает вычисление целой части выраже- |
|||||||
ния в круглых скобках, а δ п – малая поправка, которая в режиме большого изгиба ведет себя как VG5
3.
На рис. 75 показана зависимость смещения средней точки Е- нанотрубки от напряжения на управляющем электроде VG при r =
= 0,65 нм, Е = 1,25 ТПа, |
L = 500 нм, R = 100 нм. |
При VG ≈ 0,5 В |
происходит переход от |
режима малого изгиба |
трубки, когда |
zmax VG2 , к режиму большого изгиба с zmax VG2 3 |
(по данным рас- |
|
чета работы [64]). Видно, что при изменении VG поперечное смещение нанотрубки zmax изменяется ступенчатым образом.
Рис. 75. Зависимость смещения средней точки Е-нанотрубки от напряжения на управляющем электроде. Пунктирная кривая получена без учета кулоновских эффектов и влияния механического напряжения.
8.2.3. Собственные частоты механических колебаний нанотрубки
Спектр собственных механических колебаний – важная характеристика системы, которая поддается непосредственному измерению в эксперименте. Как увидим из дальнейшего, кулоновское взаимодействие в рассматриваемой NEMS сильно влияет на упругие свойства входящей в ее состав нанотрубки – настолько сильно, что качественно изменяет ее поведение.
Для того чтобы найти спектр собственных колебаний нанотрубки, можно добавить к постоянному напряжению управляющего
129
электрода малое возмущение, периодически изменяющееся во времени. Теперь поперечное смещение произвольной точки нанотрубки зависит не только от ее продольной координаты, но и от време-
ни, z = z (x,t), поэтому в левой части уравнения (8.9) появляется внешняя сила (−ρS z ), где ρ – объемная плотность трубки, ρ = 1,35
г/см3. Понятно, что по-прежнему форма нанотрубки в основном определяется большой и постоянной во времени компонентой управляющего напряжения и в среднем по времени описывается решением (8.10) статического уравнения (8.9), в котором напряжение Т – по-прежнему константа (8.11), которая определяется из условия самосогласования (8.5). Другими словами, именно большая постоянная составляющая VG определяет величину механического напряжения в нанотрубке.
Уравнение для малого возмущения поперечного отклонения, δz (x,t ), которое развивается на фоне большого статического сме-
щения (8.10), с учетом сказанного имеет следующий вид
IE δz′′′′ −T δz′′ − ρSω2 δz = 0, |
(8.15) |
а частоты поперечных собственных мод нанотрубки в системе определяются из требования, что это уравнение с граничными условиями δz (0) = δz (L) = 0, δz′(0) = δz′(L) = 0 имеет нетривиальное
решение. В итоге получаем уравнение для собственной частоты:
ch ( y1 )cos( y2 ) − |
y2 |
− y |
2 |
sh ( y1 )sin ( y2 ) =1, |
||||||
1 |
|
2 |
||||||||
2 y1 y2 |
|
|||||||||
|
|
|
|
|
|
|
(8.16) |
|||
|
L |
( ξ4 + 4λ |
2 ± ξ2 )1 2 |
|
ρS |
|||||
y1,2 = |
, λ = |
ω, |
||||||||
2 |
EI |
|||||||||
|
|
|
|
|
|
|
|
|||
а параметр ξ определен в (8.10). По причинам, которые мы обсудили ранее, в задачах наноэлектромеханики наиболее интересна так
называемая фундаментальная (наинисшая) собственная мода ω0. В двух предельных случаях, рассмотренных выше, решения уравнения (8.16) для этой моды таковы:
|
EI |
|
−2 |
+ 0.28ξ |
2 |
, ξL <<1; |
|
ω0 = |
22,38L |
|
(8.17) |
||||
|
|
|
|
||||
|
ρS |
πξL−1 + 2πL−2 , |
|
ξL >>1. |
|
||
130