

Тепловое трение и тепловое движение межузельных атомов в облучённом веществе.
Выполнил:
студент группы Т9-32а Говрас Е.А.
Москва, 2009 |
1 |
|

Модель Френкеля – Конторовой.
Лагранжиан цепочки атомов:
|
2 |
|
|
|
L |
mxn |
|
(xn 1 xn a)2 V (x1,..., xn ,...) [1] |
|
2 |
|
|||
n |
|
2 n |
α – константа взаимодействия атомов в цепочке a – равновесное положение между атомами
Потенциал окружающей решётки (д.б. периодическим с периодом a):
|
|
2 |
|
2 |
|
|
|
xn |
|
V (x1 |
,..., xn ,...) |
m 0 a |
|
sin |
2 |
||||
|
|
|
|
|
|
|
|||
2 |
2 |
|
|
a |
|||||
|
|
|
|
m |
|
|
|
ω0 - вариационный параметр амплитуды потенциала
* номер атома в цепочке n не обязательно совпадает с номером минимумов потенциала m
[1] – S.L. Dudarev, Journal of Nuclear Materials, 307-311 (2002) 881-885 |
2 |
|

Решение в виде солитона.
Пусть в цепочке есть один дополнительный атом:
|
|
u |
|
1 |
|||||||||
и поле смещений атомов цепочки слабо меняется на периоде |
a |
|
|
||||||||||
x |
|||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
||
Тогда |
|
|
2a |
|
|
|
0 (x X ) |
|
|
||||
|
|
u(x) |
|
|
|
||||||||
|
|
|
arctan exp |
c |
|
|
|
||||||
|
|
|
|
|
|
|
|
|
|||||
c |
a |
2 |
|
|
|
|
|
|
|
|
|
|
|
m |
- скорость звука в веществе |
|
|
|
|
|
|
|
|
||||
|
|
|
|
|
|
|
|
|
|
|
|
X – координата центра дефекта
Форму солитона характеризует единственный безразмерный параметр:
N |
c |
|
- эффективное число атомов, участвующих в формировании солитона |
|
0a |
||||
|
|
3

Сравнение результатов.
Сравниваем решения для одно-, двух-, трёхатомных дефектов с результатами моделирования (105 атомов)
4

Броуновское движение.
Если кластер сохраняет форму в процессе теплового движения, то движение центра кластера описывается уравнением Ланжевена:
m* X m* X (t)
m* - эффективная масса кластера
- коэффициент теплового трения
(t) |
- случайная сила: (t) T 0 |
(t) (t') |
* |
T 2mT (t t') [T]=эрг |
D mT* - коэффициент тепловой диффузии
Lb |
c |
- средняя длина баллистической траектории кластеров |
|
|
|||
|
|
5

Коэффициент трения.
Уравнение Фоккера-Планка для функции распределения движущихся дефектов:
F |
|
p F |
|
|
|
|
* |
F |
|
|
|
|
|
|
|
|
|
0 |
|
|
|
|
|
|
|
|
pF m T |
||
t |
* |
|
|
||||||
|
m x |
|
p |
|
|
p |
Выводя это уравнение из уравнений движения дефекта, можно получить выражение для коэффициента трения.
В случае наличия теплового возбуждения фононов положения минимумов решётки возле дефекта будут испытывать случайные смещения. В этом случае уравнения для смещений атомов цепочки:
m2ut2
(x,t)
|
2 2u |
|
m 02a |
|
2 [u(x,t) (x,t)] |
|||
a |
|
2 |
|
|
|
sin |
|
|
x |
2 |
2 |
a |
|||||
|
|
|
|
|
|
- поле акустических фононов решётки
(x,t) |
T |
0 - усреднение ведётся по состоянию термодинамического равновесия |
|
|
6

Коэффициент трения.
Решение ищем в виде: u(x,t) u0[x X (t)] (x,t)
u0[x X (t)] решение в виде солитона с произвольно движущимся центром
После вычислений:
1.06 4 |
a |
3 |
|
|
T |
|
ln |
|
2 c |
||||
|
|
|
|
|
|
|
2 |
|
|
||||
|
0 |
c |
|
|
|
|
0a |
||||||
|
|
|
|
mc |
|
|
|
|
|||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
0 |
|
Ta |
4 |
Ta |
|
|
|
|||||
|
|
|
|
|
|||||||||
5.53 |
|
|
|
|
c |
|
, |
c |
1 |
||||
|
|
|
|||||||||||
|
amc |
|
|
|
|
|
|
, Ta |
1 |
|
||
|
c |
|
Если с помощью этого коэффициента посчитать среднюю длину произвольного теплового смещения дефекта, получим:
o
LT 100K 15 - так как эта величина значительно превышает постоянную решётки, становится возможным объяснить «прыжки» дефектов, наблюдаемые в ММД.
7

Этапы вычисления коэффициента трения.
I. Уравнение движения цепочки атомов в присутствии теплового возбуждения фононов надо вторично проквантовать.
II. Ввести фононные переменные и получить эффективный гамильтониан взаимодействия дефекта с фононами.
III. С помощью этого гамильтониана получить кинетическое уравнение для матрицы плотности дефекта.
IV. Перейти к пределу Фоккера-Планка этого уравнения. Из вида этого уравнения для дефекта найти коэффициент трения.
F |
|
p F |
|
|
|
|
* |
F |
|
||
|
|
|
|
||||||||
|
|
|
|
|
|
|
|
|
0 |
|
|
|
|
|
|
|
|
|
|
pF m T |
|
||
|
* |
|
|
|
|||||||
t |
|
m x |
|
p |
|
|
|
p |
8 |
||
|
|
|
|
9