Компьютерное моделирование физических явлений. Малютин В.М - 46 стр.

UptoLike

ВУЗ: 

ТПУ | Томск

Составители: 

Рубрика: 

22
2
2
3
1
36
D
D
T
T
M
kT T
σ
⎡
⎤
=+
⎢
⎥
⎣
⎦
h
, (3.5)
где T – температура кристалла, M – атомная масса, k – постоянная
Больцмана, T
D
– температура Дебая.
В модели можно разыгрывать только поперечные компоненты
смещений (по x и y), так как из-за малости углов между направлением
движения частиц и осью кристалла продольными смещениями можно
пренебречь.
Окончательные координаты атома:
0
0
,
,
T
T
xx x
yy y
=
+
=+
(3.6)
где x
0
и y
0
– координаты равновесного положения атома;
2
2
2 ln( ) cos(2 ),
2ln()sin(2)
T
T
x
y ,
σ
ξπξ
σ
ξπ
=−
=−
ξ
(3.7)
где ξ — число, равномерно распределенное случайным образом на
интервале [0;1].
На рисунке 2.1 показан «эффективный канал», в который
помещается весь поток ионов при условии, что ионы между собой не
взаимодействуют. Здесь мы используем основное свойство
монокристаллов — ближний и дальний порядок.
Если в процессе многократного рассеяния ион покинул пределы
эффективного канала, но не деканалировал (угол вектора импульса иона
с кристаллографической осью меньше критического угла
каналирования), можно считать возможным искусственно поместить его
в эффективный канал при условии сохранения характеристик движения
иона в соседнем канале (его фазовые координаты относительно
ближайшего окружения атомов). Такой подход позволяет значительно
сократить время расчета на компьютере. Изменения координат и углов
при таком подходе указаны на рис 3.3.
46

Страницы