Добавил:
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:
Скачиваний:
1
Добавлен:
11.06.2024
Размер:
200.39 Кб
Скачать

Тепловое трение и тепловое движение межузельных атомов в облучённом веществе.

Выполнил:

студент группы Т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