Методика моделирования

Для проведения расчетов необходимо задать соотношение масс бомбардирующей частицы m ₁ и атома мишени m ₂, потенциал взаимодействия между ними U (r), где r – расстояние между частицами, а также расположение атомов внутри твердого тела.

Работа выполняется в вычислительном центре кафедры физических проблем материаловедения.

Приближения, используемые при моделировании

Для упрощения алгоритма в настоящей работе выбран потенциал взаимодействия твердых шаров (потенциал «жесткой сердцевины»):

U (r) = ¥, при r £ R

U (r) = 0, при r > R,

где R – радиус твердого шара. При этом траектория движения частицы внутри твердого тела представляется отрезками прямых, симметричных относительно нормали к поверхности шара в точке столкновения.

Вторым упрощением является замена объемной (трехмерной) кристаллической решетки плоской (двумерной). Такая замена позволяет не только обойтись при описании каждого акта рассеяния вычислением единственного значения угла (что ускоряет процесс моделирования), но и дает более наглядное представление о траектории при ее изображении.

В работе использован простейший тип кристаллической решетки твердого тела – простая квадратная решетка (двумерный аналог простой кубической решетки). Параметр решетки а выбран равным 0,286 нм. Внешняя граница твердого тела (положение наружных атомных слоев мишени) определяется прямыми: х = 0, x = 25 a и y = -25 а, y = 25 a.

Начальные и граничные условия

Частица начинает движение из точки с координатами í– 25 a, 0ý под углом a к горизонтали, выбираемым случайным образом в диапазоне – 30° £ a £ 30°. Она считается выбывшей из процесса взаимодействия в двух случаях: когда частица покинула пределы кристаллита (в этом случае на мониторе появляется соответствующее сообщение) и когда энергия частицы упала ниже заданной величины. Начальная энергия частицы в работе принимается равной 30 МэВ, а конечная – 1 МэВ (значения взяты условно и не отвечают какой-либо конкретной экспериментальной ситуации).