Если бы решетка была стационарной, как это считалось при рассмотрении, например, зонной теории, то потенциальная энергия кристалла была бы равна сумме вкладов от всех различных пар:
U = 1/2 | (8.1) |
Если учесть тепловые колебания, то
U = 1/2 | (8.2) |
где r(Ri), r'(R') - мгновенные положения в узле (Ri), S(Ri) = r(Ri) - Ri - смещение i-того атома из узла. Разлагая в ряд Тейлора:
f(r + a)
= f(r) + a
f(r)
+ 1/2 (a
)2f(r)
+ 1/3! (a
)3f(r)+…
полагая r = R-R' и a = S(Ri) - S(Rj) = Sij
| U = N/2 | (8.3) |
Первый член:
N/2
iV(Ri)
= ueq - есть потенциальная энергия кристалла в равновесии (каждый атом
жестко закреплен в равновесном положении).
Второй член: 1/2
ijSij
V(Rij),
здесь
V(Rij)
=
V(Ri - Rj)
- есть сила, действующая на атом i со стороны всех атомов j в равновесном положении.
Поскольку эта сила в равновесии = 0, то 1-й неисчезающей поправкой в потенциальной
энергии является 3-й член.
Третий член: 1/4
ij(Sij
)2V(Rij)
0.
Приближение, в котором ограничиваются членом ~ (Sij)2
в разложении в степенной ряд, называется гармоническим: U
Ueq
+ Uharm, где
Uharm = 1/4 | (8.4) |
4-й и следующие члены дают ангармонические поправки. С учетом ангармонизма
U = Ueq + Uharm +
Uanh | (8.5) |
На практике очень часто останавливаются на учете кубического члена, пренебрегая остальными, хотя, строго говоря, этого недостаточно.
Проявления ангармонизма:
| На главную сайта Dvoika.net |