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

UptoLike

ВУЗ: 

ТПУ | Томск

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

Рубрика: 

Print [ "заряд атома z2=",z2," заряд частицы z1=",z1];
Print [ "масса атома,а.е.м. m2=",m2, " масса частицы,а.е.м. m1="
,m1];
Print [ "пост. решетки, анг cons=",cons, "энергия частиц ,мэв el=",el
];
Print ["толщина, ч.пос.реш. nlayers=",nlayers,"число частиц nI=",
nI];
Print [ "нач. угол влета, рад tetin=", tetin];
Print [ "шаг между атомами, манг dimz=",dimz];
Print [ "радиус захвата, манг rbm=", rbm];
Print [ "радиус томаса-ферми, манг iatf=", N[iatf] ," epsi=",
N[epsi]];
Print [ "крит. угол psikr=", psikr, " рад psigr= ", N[psigr] ,"
град"];
Print [ "крит.приц.параметр rbbkr=", rbbkr];
Ion=Table[ Random[Real, maxc ],{nI},{2}];(* начальные координаты ионов
*)
tetaxy=Table[0, {nI}, {2}]; (* начальное угловое распределение *)
ListPlot[Ion,Frame->True,PlotRange->{{0,maxc},{0,maxc}},
AspectRatio -> 1,PlotLabel -> init]; (*график распределения нач.потока
*)
(* массив расстояний X и Y до атома *)
Action[at_]:=Module [ {r,dr},rxry=Table[Ion[[i]]-chanal[[at]],{i,nI}];
r=Round[ Table[Hyp[rxry[[i,1]],rxry[[i,2]]],{i,nI}] ]; (*гипотенуза *)
(*
приращение попер. расст. dx dy *)
Table[If[r[[i]]<rbm, tetaxy[[i]]=tetaxy[[i]]+an[[ r[[i]] ]]*rxry[[i]],
tetaxy[[i]]=tetaxy[[i]] ], {i,nI}];
Ion=Ion+dimz*tetaxy;
(* новые координаты без учета выхода из канала
*)
(* Анализ выхода из канала и модификация координат *)
Do[ Which [ Ion[[i,1]]<0, {Ion[[i,1]]=-Ion[[i,1]],
tetaxy[[i,1]]=-tetaxy[[i,1]] } ],{i,nI} ];
Do[ Which [ Ion[[i,2]]<0, {Ion[[i,2]]=-Ion[[i,2]],
tetaxy[[i,2]]=-tetaxy[[i,2]] } ] ,{i,nI}];
Do[ Which [ Ion[[i,1]]>maxc, {Ion[[i,1]]=2*maxc-
Ion[[i,1]], tetaxy[[i,1]]=-tetaxy[[i,1]] } ] ,{i,nI}];
Do[ Which [ Ion[[i,2]]>maxc, {Ion[[i,2]]=2*maxc-
Ion[[i,2]], tetaxy[[i,2]]=-tetaxy[[i,2]] } ] ,{i,nI}];
]
Ячейка 7: Основной цикл
Do [
{ Do[ Action[i],{i,4}];
ListPlot[Ion,Frame->True,PlotRange->{{0,maxc},{0,maxc}},
AspectRatio -> 1, PlotLabel -> n];
} ,{n,nlayers}]
▪ Возможные результаты
Выборочные результаты работы программы приведены на рис. 3.6. Но не для германия,
а для арсенида галлия GaAs. В таблице Менделеева Ga
31
, Ge
32
и As
33
расположены рядом и
их заряды, а, следовательно, и воздействия на пролетающие ионы, не сильно отличаются
друг от друга. Исследуемый поток ионов прослеживался до глубины L=668
нм.
Первоначальный поток распределен равномерно (0 нм). С ростом глубины частицы
начинают удаляться от атомных рядов к центру канала.
52

Страницы