Главная страница
Навигация по странице:

  • GEANT 4

  • 3.2. УСТАНОВКА И ЭКСПЛУАТАЦИЯ GEANT 4

  • 3.3. ОПИСАНИЕ ОСНОВНЫХ ПАРАМЕТРОВ ЗАДАЧИ

  • 3.4. РЕЗУЛЬТАТЫ ЧИСЛЕННОГО МОДЕЛИРОВАНИЯ

  • Энергетические спектры.

  • Вычисление методом Монте-Карло. ВЫЧИСЛЕНИЕ МЕТОДОМ МОНТЕ-КАРЛО КОЛИЧЕСТВА ПЕРВИЧНЫХ ДЕФЕКТОВ В К. вычисление методом монтекарло количества первичных дефектов в кристаллическом кремнии, создаваемых протонами и нейтронами различных энергий


    Скачать 0.62 Mb.
    Названиевычисление методом монтекарло количества первичных дефектов в кристаллическом кремнии, создаваемых протонами и нейтронами различных энергий
    АнкорВычисление методом Монте-Карло
    Дата13.02.2022
    Размер0.62 Mb.
    Формат файлаdoc
    Имя файлаВЫЧИСЛЕНИЕ МЕТОДОМ МОНТЕ-КАРЛО КОЛИЧЕСТВА ПЕРВИЧНЫХ ДЕФЕКТОВ В К.doc
    ТипПояснительная записка
    #360814
    страница2 из 6
    1   2   3   4   5   6


    3. ЧИСЛЕННОЕ МОДЕЛИРОВАНИЕ ВЗАИМОДЕЙСТВИЯ НЕЙТРОНОВ И ПРОТОНОВ РАЗЛИЧНЫХ ЭНЕРГИЙ С КРЕМНИЕМ С ПОМОЩЬЮ GEANT4.

     

    3.1. ФИЗИЧЕСКИЕ ПРОЦЕССЫ, ВКЛЮЧЕННЫЕ В GEANT4

     

    В данной работе проводилось моделирование взаимодействия протонов и нейтронов с энергиями 1 МэВ, 14 МэВ, 63 МэВ. В этом энергетическом диапазоне учитывались следующие физические процессы:

     compt – комптоновское рассеяние

    conv – образование e+e- пар

    eBrem – тормозное излучение

    eIoni – ионизационные потери электронов

    msc – многократное рассеяние

    annihil – аннигиляция легких частиц

    hIoni – ионизационные потери адронов

    LElastic – упругое взаимодействие адронов

    ProtonInelastic – неупругое ядерное взаимодействие протонов

    NeutronInelastic – неупругое ядерное взаимодействие нейтронов

    LCapture – захват нейтрона

    Fission – распадный процесс, вызываемый нейтронами

     К ядерным реакциям относятся ядерные взаимодействия упругого и неупругого типа, захват нейтрона, а так же распадный процесс, вызываемый нейтронами.

    Данные о сечениях взаимодействия нейтронов в Geant4 содержатся в библиотеке G4NDL, которая позволяет моделировать прохождение нейтронов с энергиями от тепловых до 20 МэВ. В диапазоне энергий до 20 МэВ в Geant4 применяются экспериментальные данные из библиотеки G4NDL, а в диапазоне выше 20 МэВ применяется следующая экстраполяция:

     

    =F(En)p12ln(N)[1+A1/3-p2(1-1/A1/3)] ,     (8)

     

    где F – корректирующий коэффициент,  N –количество нейтронов, En – энергия нейтронов в ГэВ, A – массовое число; p1 и p2 – параметры выбранной физической модели [4]. В данной работе использована модель промежуточного ядра, т.к. она наиболее точно описывает взаимодействие нейтронов с кремнием в используемом энергетическом диапазоне.

    3.2. УСТАНОВКА И ЭКСПЛУАТАЦИЯ GEANT4

    Подробная инструкция по установке Geant4  находится по адресу http://geant4.web.cern.ch/geant4/G4UsersDocuments/UsersGuides/InstallationGuide/html/index.html. Тем не менее, существует ряд нетривиальных особенностей, касающихся установки и настройки Geant4, на которые следует обратить внимание.

    1. Существует два основных варианта установки Geant4: под Windows и под Linux (или любую другую Unix-подобную систему). Для установки Geant4 под Windows потребуется среда программирования Microsoft Visual C++ 6 с установленным Service Pack 5, а так же оболочка Cygwin, эмулирующая работу Unix-терминала в Windows. Cygwin можно скачать с сайтаhttp://www.cygwin.com. После установки Cygwin процедура установки Geant4 под Windows ничем не отличается от аналогичной процедуры под Linux за исключением настройки некоторых переменных окружения, описание которых содержится в инструкции по установке Geant4.

    2. Для нормального функционирования Geant4 потребуется библиотека CLHEP(A Class Library for High Energy Physics), включающая описания генераторов случайных чисел, необходимых для работы Geant4. Найти её исходные модули и готовые для использования библиотеки можно по адресу http://cern.ch/clhep/.

    3. Далее необходимо скачать и распаковать исходный код Geant4 и файлы физических данных для ядерных и электромагнитных взаимодействий (G4NDL, G4EMLOW, PhotonEvaporation,RadiativeDecay). Последние версии всегда можно найти по адресу http://cern.ch/geant4/.

    4. После того, как все необходимые исходные коды скачаны и распакованы, можно приступать к установке Geant4 – т.е. к компиляции библиотек, входящих в состав исходного кодаGeant4. В комплекте с исходными кодами Geant4 поставляется утилита Configure, облегчающая установку Geant4. Для того чтобы начать установку, необходимо ввести команду

    ./Configure –install

    после чего утилита конфигурации в интерактивном режиме начнет сбор сведений о конфигурации системы. После того, как процесс подготовки завершен, начинается компиляция библиотекGeant4, которая занимает достаточно продолжительный промежуток времени (1-5 часов в зависимости от быстродействия ЭВМ). В процессе компиляции на экране могут появляться различные предупреждения (warning), но как правило они не оказывают отрицательного влияния на ход установки. Если процесс компиляции успешно завершён, т.е. не возникло ошибок (error), то можно считать установку пакета Geant4 завершенной.

    Порядок работы с Geant4 состоит в следующем: сначала составляется план эксперимента, изучаются основные свойства и характер взаимодействия экспериментальных частиц. Затем составляется программа на объектно-ориентированном языке C++, использующая библиотеки, входящие в состав Geant4. После этого исходные тексты программ компилируются, программа выполняется, и происходит анализ полученных результатов. По мере  анализа полученных данных в программу вносятся изменения, затем процесс повторяется до тех пор, пока не будут получены удовлетворительные результаты.

    Моделирование физического эксперимента с помощью Geant4 основано на методе Монте-Карло. Это предполагает ход физических процессов в эксперименте с наперед заданными вероятностями, чтобы обеспечить максимальную реалистичность эксперимента.

    Моделирование прохождения нейтронов в веществе включает в себя: описание геометрии детектора, описание химического состава входящих в эксперимент материалов, описание типов частиц, участвующих в эксперименте, описание физических процессов и первичного события (вылет первичной частицы из источника).

     

    ● Описание химического состава материалов, входящих в эксперимент:

     В качестве материала для детектора был использован природный  кремний с плотностью 2,33 г/см3 и молярной массой 28,09 г/моль. Состав природного кремния: Si28 – 92,23%, Si29 – 4,68%, Si30 – 3,09%.

     

    ● Описание геометрии детектора:

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

     

    ● Описание типов частиц, участвующих в эксперименте:

    Для того, чтобы в эксперименте участвовали какие-либо частицы, их необходимо вначале описать. В нашем эксперименте будем рассматривать протоны, нейтроны, гамма-кванты, электроны, позитроны, μ+, μ- , а также альфа-частицы, дейтроны и другие частицы.

     

    Рассмотрим пример выходных данных эксперимента прохождения нейтрона с энергией 15 МэВ в кремниевом кубике 20х20х20 см3. Данные представлены в текстовом виде. Каждая частица имеет свой номер трека (Track ID) и номер родительской частицы (Parent ID). На каждом шаге (Step) указываются координаты частицы (X, Y, Z), её энергия (KinE), энергетические потери (dEStep),  длина пробега на текущем шаге (StepLeng), полная длина пробега (TrackLeng), место нахождения (Volume) и тип взаимодействия (Process).

    ***********************************************************************************    G4Track  Information:      Particle = neutron,          Track ID = 1,           Parent ID = 0

    **********************************************************************************

    Step#      X             Y            Z        KineE      dEStep    StepLeng   TrakLeng  Volume     Process

        0     -10 cm        0 fm      0 fm     15 MeV     0 eV       0 fm             0 fm        World      initStep

        1    -2.49 cm      0 fm      0 fm     15 MeV     0 eV     7.51 cm        7.51 cm    World   NeutronInelastic

     

    Визуализация осуществлялась с помощью программного средства FukuiRenderer DAWN, позволяющего отображать траектории частиц и составляющие части детектора в виде трехмерного изображения. Трехмерное изображение представлено в виде проекции на плоскость (Рис. 4) с нанесенными на неё траекториями и энергиями частиц, участвующих во взаимодействиях, произведенных первичной частицей.

     



    Рис. 4

     

    Примечание: полная история движения и взаимодействий вторичных частиц, соответствующая этому рисунку, представлена в приложении 1.

    3.3. ОПИСАНИЕ ОСНОВНЫХ ПАРАМЕТРОВ ЗАДАЧИ

     

    Численное моделирование проводится методом Монте-Карло с помощью пакета Geant4. Процедура моделирование прохождения нейтронов и протонов в кремнии представляет повторяющийся процесс запуска начальной частицы (протона или нейтрона) с заданной энергией. После того, как частица вылетела из своего начального положения (источника первичных частиц) и провзаимодействовала в объеме детектора, регистрируются координаты и энергии вторичных частиц, заполняются угловые и энергетические спектры. Затем процесс продолжается до тех пор, пока не будет накоплено нужное количество ядер отдачи.

    Блок схема для программы моделирования и общая картина эксперимента представлены ниже.





    3.4. РЕЗУЛЬТАТЫ ЧИСЛЕННОГО МОДЕЛИРОВАНИЯ

     

    Моделирование проводилось на операционной системе Linux Mandrake 9.0. Для вычислений использовалась ЭВМ с процессором 1,3 МГц, оперативной памятью 256 Мб. Версия пакетаGeant4 4.4.1 с патчем N1. Последнюю версию Geant4 можно скачать с сайта http://cern.ch/geant4/.

     

    Энергетические спектры.

     

    В результате моделирования прохождения протонов и нейтронов в природном кремнии были получены энергетические спектры ядер отдачи и угловые распределения. В таблице 1 представлены условия и результаты моделирования.

     

    Таблица 1.

     

    Толстая пластина, 50 см

    Тонкая пластина, 1мм

    Час-тица

    n

    p

    n

    p

    E

    1 МэВ

    14 МэВ

    63 МэВ

    14 МэВ

    63 МэВ

    1 МэВ

    14 МэВ

    63 МэВ

    14 МэВ

    63 МэВ

    n

    3,5•104

    1,3•105

    105

    3•108

    1,1•106

    1,7•106

    106

    3•105

    106

    1,5•105

    N

    106

    106

    106

    4•104

    2•105

    104

    104

    104

    104

    104

    t

    4

    10

    25

    10

    18

    5

    3

    5

    3

    9

    ΔE

    32029

    131117

    283965

    6075

    103456

    388

    4271

    9212

    1684

    9224

    ζ

    91,5

    7,2

    4,5

    0,0001

    0,15

    0,023

    0,03

    0,05

    0,01

    0,1

    D

    3660

    4034

    11358

    0,08

    376

    0,9

    17,1

    123

    7

    246

     

     

    E – энергия налетающих частиц, n – количество налетевших частиц,

    N-количество образованных ядер отдачи, t-время моделирования

    (минуты), ΔE-полная энергия образованных ядер отдачи (МэВ),

    ζ – доля неионизационных потерь энергии по отношению к полным потерям энергии (%),

    D – количество дефектов от одной налетающей частицы.

     

     

    Доля неионизационных потерь энергии вычислялась по формуле:

    .          (9)

    Количество дефектов от одной налетающей частицы определяется выражением:

    .      (10)

     



    Рис.7



    Рис.8



    Рис.9



    Рис.10



    Рис.11



    Рис.12



    Рис.13

     



    Рис.14

     



    Рис.15



    Рис.16

     

    Примечание: на рис. 10 приведены только неупругие взаимодействия. Количество упругих  взаимодействий равно 26272. Все упругие взаимодействия образуют ядра отдачи  c энергией до 10 кэВ. На рис. 11 не показаны ядра отдачи  с энергией меньше 10 кэВ. В первый бин спектра входят 100258 ядер от упругих взаимодействий. На рис. 14, рис. 15, рис. 16  приведены результаты только неупругих взаимодействий. Упругие взаимодействия составляют  

    60% от общего числа взаимодействий и входят в первый бин спектра.
    1   2   3   4   5   6


    написать администратору сайта