UDC 539.19
DOI: 10.18698/1812-3368-2019-1-27-36
DEVELOP A MODEL TO STUDY THE ENERGY DISTRIBUTION OF CASCADES OF ATOMIC COLLISIONS
A.A. Aleksandrov1 V.A. Akatev1 E.V. Metelkin2 E.Yu. Baryscheva1
[email protected] [email protected] [email protected] [email protected]
1 Bauman Moscow State Technical University, Moscow, Russian Federation
2 Russian State Social University, Moscow, Russian Federation
Abstract
In this paper, based on the solution of the Boltzmann kinetic equation, we determine the energy distribution function describing the steady-state deceleration of the cascade of moving atoms taking into account their multiplication at the power interaction potential U ~1/ rn. A new approach to the solution of the kinetic equation based on the extended concept of primary knocked-on atoms (PKA) is used for its calculation. One of the advantages of using the power interaction potential is that in this case it is possible to obtain simple analytical formulas for the distribution function of the cascade of slowing-down atoms taking into account their multiplication and demonstrate the simplicity and convenience of the proposed new approach to the solution of the kinetic equation. On the other hand, based on the obtained results it is possible to estimate the accuracy of various approximate solutions. It is shown that this approach will be applicable to other interatomic interaction potentials, if the average PKA energy loss in individual collisions decreases mono-tonically with decreasing energy, and the relative PKA energy loss in individual collisions will be small
Keywords
Kinetic equation, collisions, cascade of atoms, slowing-down of atoms, interaction potential
Received 24.04.2018 © Author(s), 2019
The development of nuclear reactors and thermonuclear installations involves the selection of radiation-resistant materials as their bodies and individual elements must withstand long-term exposure to radiation. Irradiation of solids by fast particles (neutrons, fast ions) leads to the fact that the atoms of the crystal lattice received from the incoming particle energy greater than a certain threshold
value fly out of their equilibrium positions. Further, the exchange of energy between moving atoms and atoms located in the nodes of the crystal lattice leads to the emergence of new generations of knocked-on atoms. Thus, there is a so-called cascade of atomic collisions. As a result of the cascade development the primary point defects (vacancies, interstitial atoms) are formed in the solid body, which determine the accumulation of radiation damage in irradiated materials and subsequent changes in their physical and mechanical properties [1, 2].
For the correct description of these processes, it is necessary the information on the nature of the energy distribution of the decelerating cascade of moving atoms. A significant number of works are devoted to the study of the development of atomic collision cascades in a solid, for example [3-9]. To be noticed is that the authors of works [7, 8] found a distribution function that describes the energy distribution of a cascade of moving atoms taking into account their multiplication for the simplest special case of elastic spherically symmetric scattering in the system of the center of mass without taking into account the binding energy of atoms in the lattice sites. There are no previous results for more complex interaction potentials of simple analytical solutions describing the cascade of decelerating atoms.
In this paper, based on the solution of the Boltzmann kinetic equation we determine the energy distribution function describing the steady-state deceleration of the cascade of moving atoms taking into account their multiplication at the power interaction potential U ~1/ rn [10]. A new approach to the solution of the kinetic equation based on the extended concept of primary knocked-on atoms (PKA) is used for its calculation. Henceforward, the PKA is that atom, which after the collision has more energy [6-9]. An atom with lower energy will be called the target atom (or knocked-on atom).
One of the advantages of using the power potential of interaction is the possibility of obtaining simple analytical formulas for the distribution function of the cascade of decelerating atoms taking into account their multiplication, which makes it possible to demonstrate the simplicity and convenience of the proposed below new approach to the solution of the kinetic equation. Using the results that we obtained it is also possible to estimate the accuracy of various approximate solutions.
Note that in the work [11] the authors determined the spatial and time characteristics of the primary knocked-on relativistic electrons (PKRE) decelerating in the substance due to ionization losses: the dependence of the PKRE energy on time, the dependence of the energy of the decelerating PKRE on the distance covered and the dependence of the distance covered by PKRE on the deceleration time.
Problem statement. Consider a solid body consisting of atoms of the same type, whose binding energy at the nodes of the lattice is ed. Let us study in more detail the interaction between a moving atom and an atom located in a lattice node. Let the moving atom with energy E' as a result of collision with an atom of the crystal lattice becomes a state with energy E, by transmitting energy to the atom located in the node of the crystal lattice, E' - E -ed. In this case, there is only one of two possible events can be realized: 1) the energy of a moving atom E after elastic scattering, there may be more energy of the knocked-on atom E' - E - ; 2) the energy knocked-on of an atom E' - E -ed there may be more energy from the moving atom after the impact.
The kinetic equation describing the unsteady energy distribution of moving atoms (both PKA and knocked-on) from an instantaneous uniformly distributed over the space of a monoenergy source without taking into account their collisions with each other is written as [12]
Eo
+ л(Ес-E-sd) j dE'Z(E' ^E'-Esd)Ф(Е',t) + N08(E ^Eo)S(t), (1)
where f (E, t) dE is number of atoms with energy E in the range dE at time t per unit volume; 2 (E) is full macroscopic cross section of atomic scattering; 0( E, t) = vf (E, t) is the flow of moving atoms, v is the atoms' speed; 2 (E' ^ E) is macroscopic cross section of atomic scattering; 5(x) is the Dirac Delta function; N 0 is total number of atoms; E0 is initial energy of moving atoms; n (x) is Heaviside unit function,
The first integral on the right side of the kinetic equation (1) describes the transition of a moving atom with energy E' to the state with energy E. In this case, the knocked-on atom receives kinetic energy E' - E -ed. The second integral describes the formation of a knocked-on atom with energy E, when a moving atom moved into a state with energy E' - E -ed.
According to the terminology adopted above, we consider PKA to be the one that has more energy after each act of interaction. Obviously, contributions to the kinetic Eq. (1) from both the first and the second integral should be taken into account when the condition is met E < E' - E -ed or
f (E, t) dt
Eo
+ E(E) Ф (E, t) = j dE' Z (E ' ^ E) Ф (Et) +
E
E+&d
E' < 2E + sd. In this case, the kinetic equation describing the energy distribution of PKA will be written as [7, 8]:
. t-w E° -sd
- at E >-
2
MM + E(E) O (E, t) = j° dE' E (E ' ^ E) O (E ', t) +
dt e
E°
+ ^(E°-ESd) j dE'E(E'^E-E-zd)O(E,t) + tf0S(t)S(E-E°); (2a)
E+sd
. -o -sd
- at E <-
2
f(E t) 2E+sd
ff(Eli) + E(E) O(E, t) E= j dE' E(E ' ^ E) O(E', t) +
dt E
2E+sd
+ j dE ' E (E ' ^ E ' - E -sd) O (E', t). (2b)
E+sd
Based on the law of conservation of energy from Eq. (2) in the works [7, 8] the authors found that the average energy loss of PKA in a separate collision at E ' > sd described by the expression for
E' (E' -sd )/2
A°(E') = j dE (E' - E)P(E' ^ E) + j dE (sd + E)P (E' ^ E), (3)
(E' -sd V2 °
where P (E ' ^ E) = E (E' ^ E^E (E' ) is the probability of transition of a moving atom in interaction with a stationary atom from a state with energy E in a state with energy E in a single interval of energies.
The result (3) has the following interpretation. Broken atoms can be formed when E' >sd. If atom with energy E', colliding with the lattice atom acquires energy E lying in the interval from (E '-sd)/2 till E', then it knocked-on the atom has less energy, and then it is not considered. In this case, the energy loss of the moving atom will be (E ' - E), which reflects the first term in (3). If the incoming atom with energy E' as a result of the interaction, it moved into a state with energy E, lying in the interval 0...(E '-sd)/2, the great energy acquired knocked-on atom, which are followed in the future. The atom itself with energy E excluded from consideration, so that the loss of energy is (E + sd), and what reflects the second term in the formula (3).
Exact analytical solution of Eq. (2) for a simple special case of elastic, spherically symmetric in the system of the center of mass of scattering of parti-
cles with the same mass at sd = 0 obtained in the works [7, 8]. The solution of Eq. (2) for the power potential of interaction between atoms will be obtained below and the distribution function describing the steady-state deceleration of the cascade of moving atoms taking into account their multiplication will be determined on its basis.
Characteristics of the deceleration process of moving atoms. For the power potential of interaction between atoms U ~ 1/ rn the differential cross section of the scattering of a moving atom on the lattice atom is written as [10]
E (E' ^ E) = E0
1
'm/ p' p\1+m '
E'm (E' - E)
(4)
where E 0 = const; m = 1/ n.
Using (4), we find the complete macroscopic cross section
E'-Ed
E (E') = { E (E' ^ E) dE =
mE 'mEm
1 -I Ed E'
(5)
In obtaining the formula (4) it was assumed that the energy of a moving atom exceeds the binding energy of the atom at the lattice site E ' > sd.
Using (4), (5), the formula (3) can determine the average energy loss of PKA energy in a separate collision. In a dimensionless form with A(s) = A0jsd is
- at m = 1
A(E) =
1 -1
ln||l + 2ln f 1 +1 ] +1 - 1
(6)
- at m < 1
A(E) =
m e
1-m
1 -Ii
1
m (1 - m)
1-m
2m I 1 + -I -1
(7)
Here s = E'j sd — the dimensionless energy.
The relative energy losses of the moving PKA in a separate collision are determined by the ratio
A(s)
5 (E) = -
(8)
Energy dependence of the average energy losses of PKA in individual collisions, constructed by formulas (6), (7) for various exponents n (respectively
ь
E
m = 1; 2/3; 1/2; 2/5; 1/3; 1/4), shown in Figure, a. The average energy losses when decelerating PKA decreases and have a maximum value when the initial energy. This value can be found by formulas (6), (7), taking e »1: - at m = 1
Am = ln ^ +1; (9)
- at m < 1
e0-m (2m -1) , x
Am --. (10)
1 - m
It should be noted that the calculations carried out by the formulas (9), (10) are in good agreement with the results shown in the Figure, a.
According to the accepted terminology, PKA transmit energy to knockedon atoms, the latter appear only at energies e< A m. Consequently, in the field of energy e> A m there are only PKA and in this area their distribution function describes the entire cascade of decelerating atoms. Thus, in the field of energy
Am <e<80 (11)
the energy distribution of the cascade of decelerating atoms taking into account their multiplication is described by the PKA distribution function, which is the solution of Eq. (2).
It should be noted that in the works [7, 8] the authors found exact solutions of Eq. (1), (2) for a simple special case of elastic, spherically symmetric scattering in the system of the center of mass without taking into account the binding energy of atoms in the lattice nodes. In this case, the PKA distribution function (solution of Eq. (2)) coincided in the energy region (11) with the energy distribution of the cascade of decelerating atoms taking into account their multiplication (solution of Eq. (1)).
The energy dependence of the relative energy losses of PKA in individual collisions, built on the formula (8) using the results (6), (7), is shown in the Figure, b. The relative energy losses of the decelerating PKA in individual collisions are small, which makes it possible to use the approximation of continuous energy losses to determine the distribution function of the decelerating PKA [12].
The stationary distribution function of decelerating the primary knocked-on atoms. Since the relative energy losses of the decelerating PKA in
Energy dependences of average energy losses of PKA in separate collisions, constructed by formulas (6), (7) (a) and (8) using the results (6), (7) (b), class of PKA distribution functions (c) at value n = 1 (1), 1,5 (2), 2 (3), 2,5 (4), 3 (5), 4 (6)
Ф(Е)
103 102 10
1
10
,-1
10" 1(T
2 104 4-104 6-Ю4 8-104 e
a separate collision are small, the approximation of continuous energy losses can be used to determine the stationary solution of Eq. (2) [12]. In this approximation, the stationary distribution function of decelerating PKA is written as [12]
= Sc8(£-E0)+_J^, (12)
V ' Z (E) A0(E)I(E) V '
where S0 is the number of PKA emitted by a source per unit time per unit volume.
In expression for (12), the first term describes the scattered radiation emitted by the source, i.e., atoms with energy E0, and the second-the scattered radiation, i.e., atoms slowed down further to lower energies.
Using (5)-(7), the distribution function (12) can be represented in dimen-sionless form:
- at m = 1
Ф'(е) =
- at m < 1
ms
i-il
S (s-so) +
lni -| + 2ln f 1 +1 1 +1 -1
(13)
Ф'(е) =
ms
1 - 1
S (s-so) +
,2m-1
1
m (1 - m)
2m i 1 + -
1-m
-1
(14)
J_ ms
where 0' = OEo/(S0sdm_1).
Family of PKA distribution functions (second term in (13), (14)) for several exponents n shown in Figure, c. In the energy range (11), this result is the energy distribution of the cascade of decelerating atoms taking into account their multiplication.
Summary. In this paper, based on the solution of the Boltzmann kinetic equation we determine the energy distribution function describing the steady-state deceleration of the cascade of moving atoms taking into account their multiplication at the power interaction potential U ~1/rn [10]. A new approach to the solution of the kinetic equation based on the extended concept of primary knocked-on atoms (PKA) is used for its calculation. This approach can also be applied to other interatomic interaction potentials if the average PKA energy losses in individual collisions decrease monotonically with decreasing energy, and the relative PKA energy losses in individual collisions are small. In this case, the distribution function describing the steady-state deceleration of the cascades of moving atoms taking into account their multiplication in the energy region (11) will be determined by the expression for (12), in which the values of Z(E) and A0(E) are found by expressions for (3), (5).
Translated by K. Ivanova
REFERENCES
[1] Lehmann Chr. Interaction of radiation with solids and elementary defect production. Kernforschungsanlage Juelich, 1977.
[2] Was Gary S. Fundamentals of radiation materials science. Metals and alloys. Springer, 2010.
[3] Ryazanov A.I., Metelkin E.V. Concerning the theory of radiation cascades of atomic collisions in a solid with an arbitrary interatomic interaction potential. Radiat. Eff., 1980, vol. 52, iss. 1-2, pp. 15-23. DOI: 10.1080/00337578008210012
[4] Satoh Y., Kojima S., Yoshiie T., et al. Criterion of subcascade formation in metals from atomic collision calculation. J. Nucl. Mater., 1991, vol. 179-181, part 2, pp. 901-904. DOI: 10.1016/0022-3115(91)90234-X
[5] Satoh Y., Yoshiie T., Kiritani M. Binary collision calculation of subcascade structure and its correspondence to observe subcascade defects in 14 MeV neutron irradiated copper. J. Nucl. Mater, 1992, vol. 191-194, part B, pp. 1101-1105.
DOI: 10.1016/0022-3115(92)90645-2
[6] Metelkin E.V., Ryazanov A.I. Threshold energy of formation of subcascades. At. Energy, 1997, vol. 83, iss. 3, pp. 653-657. DOI: 10.1007/BF02415246
[7] Metelkin E.V., Ryazanov A.I., Semenov E.V. Developing new theoretical models of the formation of atomic collision cascades and subcascades in irradiated solids. J. Exp. Theor. Phys., 2008, vol. 107, no. 3, pp. 394-404. DOI: 10.1134/S1063776108090070
[8] Metelkin E.V., Ryazanov A.I., Semenov E.V. Ispol'zovanie kineticheskogo uravneniya Bol'tsmana dlya opisaniya ansamblya dvizhushchikhsya atomov, vybitykh iz uzlov kristal-licheskoy reshetki, dlya issledovaniya obrazovaniya kaskadov i subkaskadov atomnykh stolknoveniy v obluchaemykh tverdykh telakh [Using Boltzmann kinetic equation for description of moving atoms assembly, knocked-out from lattice sites for study on formation of atomic collision cascades and subcascades in radiated solid bodies]. Preprint Kurchstov Institute. Moscow, NRC Kurchatov Institute Publ., 2008 (in Russ.).
[9] Ryazanov A.I., Metelkin E.V., Semenov E.V. Modeling of cascade and subcascade formation at high PKA energies in irradiated fusion structural materials. J. Nucl. Mater., 2009, vol. 386-388, pp. 132-134. DOI: 10.1016/j.jnucmat.2008.12.071
[10] Lindhard J., Vibeke Nielsen, Scharff M. Approximation method in classical scattering by screened Coulomb fields. Kopenhagen, Munksgaard, 1968.
[11] Akat'yev V.A., Metelkin E.V., Savinov A.M. Space-time characteristics of relativistic electrons moderated in matter. At. Energy, 2017, vol. 122, iss. 5, pp. 364-368.
DOI: 10.1007/s10512-017-0280-6
[12] Isakov A.I., Kazarnovskiy M.V., Medvedev Yu.A., et al. Nestatsionarnoe zamedlenie neytronov. Osnovnye zakonomernosti i nekotorye prilozheniya [Non-stationary moderation of neutrons. Main laws and some applications]. Moscow, Nauka Publ., 1984.
Aleksandrov A.A. — Dr. Sc. (Eng.), Professor, Rector of Bauman Moscow State Technical University, Head of Department of Ecology and Industrial Safety, Bauman Moscow State Technical University (2-ya Baumanskaya ul. 5, str. 1, Moscow, 105005 Russian Federation).
Akatev V.A. — Dr. Sc. (Eng.), Professor, Department of Ecology and Industrial Safety, Bauman Moscow State Technical University (2-ya Baumanskaya ul. 5, str. 1, Moscow, 105005 Russian Federation).
Metelkin E.V. — Dr. Sc. (Phys.-Math.), Professor, Department of Technosphere Safety and Ecology, Russian State Social University (Vil'gel'ma Pika ul. 1, str. 4, Moscow, 129226 Russian Federation).
Baryscheva E.Yu. — Graduate Student, Department of Ecology and Industrial Safety, Bauman Moscow State Technical University (2-ya Baumanskaya ul. 5, str. 1, Moscow, 105005 Russian Federation).
Please cite this article as:
Aleksandrov A.A., Akatev V.A., Metelkin E.V., et al. Develop a Model to Study the Energy Distribution of Cascades of Atomic Collisions. Herald of the Bauman Moscow State Technical University, Series Natural Sciences, 2019, no. 1, pp. 27-36. DOI: 10.18698/1812-3368-2019-1-27-36
В Издательстве МГТУ им. Н.Э. Баумана вышел в свет учебник авторов
В.И. Ванько, О.В. Ермошиной, Г.Н. Кувыркина
«Вариационное исчисление и оптимальное управление»
Наряду с изложением основ классического вариационного исчисления и элементов теории оптимального управления рассмотрены прямые методы вариационного исчисления и методы преобразования вариационных задач, приводящие, в частности, к двойственным вариационным принципам. На примерах из физики, механики и техники показана эффективность методов вариационного исчисления и оптимального управления для решения прикладных задач. Для студентов и аспирантов технических университетов, а также для инженеров и научных работников, специализирующихся в области прикладной математики и математического моделирования.
По вопросам приобретения обращайтесь:
105005, Москва, 2-я Бауманская ул., д. 5, стр. 1 +7 (499) 263-60-45 [email protected] http://baumanpress.ru
В.И. Ванько, О.В. Ермошина, Г.Н. Кувыркин
Вариационное исчисление и оптимальное управление