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

UptoLike

ВУЗ: 

ТПУ | Томск

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

Рубрика: 

В качестве изменяемых параметров в нашем случае берутся число ионов в потоке и
число слоев. Слой в данном случае для канала <100> равен постоянной решетки, в нем
последовательно рассматривается взаимодействие потока ионов с 4-мя атомами, которые
образуют <100>-канал в элементарной ячейке кристалла.
Далее следует ячейка с операторами печати параметров, операторами, задающими
равномерно распределенный начальный поток ионов с помощью метода Монте-Карло.
Также в эту ячейку входят операторы, формирующие ключевую функцию Action для
расчета траекторий ионов. После прохода каждого атома осуществляется анализ выхода
частицы из канала и модификация пространственных и угловых координат (возврат в
«эффективный» канал) (рис 3.3). После выполнения этой ячейки происходит печать
параметров и графический вывод координатного распределения начального потока.
И, наконец, задаем основной цикл для функции Action с графическим выводом
распределения потока после каждого слоя.
▪ Программа
z2=32; (* заряд атома *) z1=2; (* заряд частицы *)
Ячейка 1: Задание параметров
m2=72.59 масса атома а.е.м. *) m1=4.003; (* масса частицы,а.е.м.*) ;(*
cons=5.657; (* пост. решетки, А *) el=1.8; (* энергиЯ частиц ,Мэв *)
elkev=el*1000; (* энергиЯ частиц, кэв*)
tetin=0; ч. угол влета, рад *) (* на
dimz=1414.245; (* шаг между соседними атомами в канале , мA*)
maxc=rbx=rby=2000; (* граница канала , мA*)
chanal={{0,0},{0,2000},{2000,2000},{2000,0}}; (* координаты атомов в
канале *)
(* Трансф. коэффициенты *)
radgr=1/Degree; (* рад в град *) angma=1000; (* ангстрем в мА *)
maang=0.001; (* мА в А *)
atf=0.468345/Sqrt[ z1^(2/3) + z2^(2/3) ]; (* радиус Томаса-Ферми, A *)
Ячейка 2: Рабочие константы
psikr=psich=Sqrt[ z1*z2*2.8804*10^(-5)/(cons*el)]; (* критический угол,
рад *)
psigr=psichg=psich*radgr; (* критический угол, град *)
teting=tetin*radgr; iatf=atf*angma; epsi=z1*z2*14.4*10^(-6)/(el*atf);
Hyp[n_,m_]:=N[Sqrt[n^2+m^2]]; (* гипотенуза *)
Ячейка 3: Функции
(* угол рассеяния *)
ra=rb/atf;
Ansin[rb_,atf_,epsi_]= epsi* (0.6*BesselK[1, 6*ra]+
0.66*BesselK[1, 1.2*ra]+0.105*BesselK[1, 0.3*ra]);
(* критический прицельный параметр *)
rbbkr=rb/.FindRoot[ rb-Ansin[ rb,atf,epsi]==0, {rb,10^(-10),cons}];
rbm=3000; (* радиус охвата рассеянием , мА *) b=istep*maang; (* A
*)
Ячейка 4: Массив углов рассеяния
an=N[ Table[Ansin[ b,atf,epsi]/istep,{istep,rbm}] ] ; (* рад/мА *)
nI=1000 ; (* число ионов *)
Ячейка 5: Изменяемые параметры
nlayers=20; (* толщина , число. постоянных решетки *)
dte=Date[];
Ячейка 6: Расчет распределения потока
Print [ " ",dte[[4]],":",dte[[5]]," ",
dte[[3]],"-",dte[[2]],"-",dte[[1]] ];
51

Страницы