Research Letters in Physics

Volume 2008, Article ID 746892, 4 pages

http://dx.doi.org/10.1155/2008/746892

## Computer Modeling of Displacement Cascades in Beryllium Irradiated with Intensive Neutron Flux

Institute for Nuclear Research and Nuclear Energy, 72 Tzarigradsko Chaussee Boulevard, Sofia 1784, Bulgaria

Received 21 October 2007; Accepted 3 January 2008

Academic Editor: William Dorland

Copyright © 2008 T. Troev et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

#### Abstract

Computer simulations of the radiation defects created in beryllium irradiated by fast neutrons ( MeV) using the Geant4 and SRIM packages were carried out. The atom cascade displacements in Be at a neutron fluence of n/c were determined to be 0.06 dpa and the helium concentration was calculated to be 168 appm. The concentration of has been estimated to be 5% in comparison to the He concentration. Nanoscale calculations were done in nm cube of fast neutron-irradiated Be. A correlation between the Be primary knock-on atom (PKA) energies and the damage cascades has been established. The final defect distributions of single vacancies, divacancies, and small vacancy clusters were examined. Our results indicate that the damages caused by He atoms are about 3 times less than damages caused by Be primary knock-on atoms (PKAs).

#### 1. Introduction

The purpose of our study is the computer modeling of the damage of a Be sample caused by intensive fusion neutron irradiation. In order to achieve our goal, we have adopted similar approach used by Becquart [1]. First, we have simulated the interaction of neutrons with a macroscopic Be target at high energy ( MeV). For this purpose, the fusion neutron spectrum data was taken from the experiment as an input for the simulation. We have recorded the initial energy and position of Be primary knock-on atoms (PKAs) and of all isotopes created in the nuclear reactions at the output of the calculation. The results obtained allowed us to calculate the concentration of He and Li isotopes and use it as an input for the next step in our simulation. Using the Monte-Carlo (MC) SRIM [2], code we have calculated the atom displacements caused by the Be and He PKAs. At the end of the simulation, the final positions of the atoms involved in damage cascades have been recorded. Finally, the atom displacements and the number of produced vacancies and interstitials as a function of the PKA energy for the isotopes created were calculated.

#### 2. The Methodology of the Simulations

We created a computer Monte-Carlo (MC) code based on the Geant4 CERN package [3] to simulate
ion and neutron interactions at high energies ( MeV) in a
Be homogeneous target. We took the initial neutron energy distribution for the
simulation from [4]. Neutron elastic and inelastic scattering, neutron capture
with and without emission of charged particles were taken into consideration.
Gamma ray and secondary electron (SE) emission as well as all types of
electromagnetic interaction of charged particles and photons with matter were
included in the present model. The high-precision database for low-energy
neutron cross-section G4NDL [3] based on the ENDF/B-VI [5] database has been
used with the purpose to achieve high accuracy of the neutron interaction in Be
at energies between 0.1 MeV and 14 MeV. In our numerical calculations, we have
considered the production of ^{3}He, ^{3}He, and ^{6}Li. The ^{7}Li production has a threshold
of 12 MeV, and therefore it takes place in a fusion reactor. Our computer
simulations were carried out for a macroscopic cylindrical Be target with
diameter and length of 5 cm. We should point out that the Be PKAs for the
damage cascades due to neutron scattering and capture were considered for a
neutron fluence of about n/cm^{2} ( MeV) and for small cubic volume with dimensions nm inside the target. We have adopted the procedure, described in
[4] for estimation of the damage in beryllium due to intensive neutron
irradiation. The results obtained for dpa (displacement per atom) cross-section
versus various neutron energy are shown in Figure 1 where
it is seen that the elastic neutron scattering ()
cross-section is more than an order of magnitude higher than the cross-section
of the neutron inelastic scattering (). The
magnitude of the neutron capture cross-section in Be with gamma ray
emission is relatively small but it is important for the tritium production.
Our calculations show that the threshold for the ()-induced dpa
is about 2 MeV and for the () is about 1
MeV. The dpa due to elastic scattering was present at over whole energy range
(0.1–14 MeV) in the
present simulation. The dpa cross-sections for () and () grow rapidly up to about 4 MeV, after that they reach saturation. The next step in our model
was to carry out calculations of damages in Be using the results from the
described high-energy neutron interactions. We used the binary collision
approximation (BCA) and the SRIM-MC code [2] which relies on the universal pair
potential of Ziegler et al. [6] and Robinson [7]. The SRIM code allows [1] the
calculation of the stopping of the PKA at energies higher than 10 keV/u with an
accuracy of 10%. In this approach, the crystal structure of Be has been
neglected. Every PKA created by the high-energy neutron interaction with Be
atoms was tracked in SRIM to its stopping in the Be target. According to the
model, the target atom which is hit by PKA and starts a recoil cascade is
identified by its recoiling energy. The number of displacement collisions
records how many target atoms were set in motion in the cascade with energies
above their displacement energy which is specified at the input. The Be target
vacancies are the next items which were recorded. The replacement collisions on
the other hand may reduce the total vacancies up to 30%. In the model, the
average displacement energy has been chosen to be 40 eV.
The lattice binding energy was taken to be 3.32 eV equal to the cohesive energy
of the Be crystal. At the end of every cascade, the ion type, the PKA energy,
the recoil energy, and the atom positions have been recorded for further
processing.

#### 3. Results and Discussion

In the present simulations, the He content was estimated to be 168 appm. The ^{6}Li atom concentration was
estimated to be 5% of the He content for the neutron fluence of n/cm^{2} ( MeV) used
in the simulations. If we scale the result for the neutron flux of n/cm^{2}/s expected for the fusion
reactor first wall (FW), our data is in agreement
with the available one in the literature [8] although our results are obtained
by a different procedure. It was established in our calculations that the total
damage of 0.06 dpa for a neutron fluence of n/cm^{2} scales well with the already
reported results for n/cm^{2} [4]. We have established that
the calculated number of the replacement collisions which lead to the
recombination of vacancies and self-interstitials was only 1.8% compared to the
total number of displacements. This result for the first time indicates that
high-energy recoils generate atomic collision cascades in Be in which a
fraction of the defects recombines. The calculated results of the Be subcascades as a function of the PKA energy are
shown in Figure 4. They indicate a correlation between PKA energies and the
number of subcascades. It is seen in Figure 4 that PKAs with higher energies
produce a large number of subcascades. It is seen that Be PKA energies in the range
20–120 keV have a higher probability to induce
subcascades. Obviously, the most probable number of subcascades is around 50.
Similar results but for -Fe are discussed in [9, 10]. The authors pointed out that for PKA with an energy
higher than some critical value in the range of
10–40 keV, the formation of subcascades is more
probable. On the other hand, PKA with higher energies do not initiate as large
a number of subcascades as one could expect. This is probably due to the fact
that the most efficient energy transfer in Be occurs at lower energies and thus
produces the greatest amount of damages. From our simulation results shown in
Figure 2, it is evident that low-energy cascades ( keV) are
responsible for about 70% of the defects production. The accumulation of He
atoms in the Be sample is due to the stopped He atoms losing their energy in
the cascade collisions. The distribution of vacancies created by beryllium atom
displacements and helium atoms created by neutron capture in Be is shown in
Figure 3. We considered the He atoms as
interstitials which contribute to the interstitial clusters concentration.

#### 4. Conclusions

Computer simulations of the creation of point defects in beryllium, irradiated by
intensive fast ( MeV) neutron fluence up to n/cm^{2}, have been performed. Our
numerical model includes two steps. In the first step, precise calculations of
the PKA energies and positions for Be, He, and ^{6}Li recoils have been made. The
damage in Be due to the intensive neutron irradiation was calculated to be 0.06
dpa. The He content was determined to be 168 appm which is in agreement with
the available results. The ^{6}Li content which is important
for the tritium production was calculated to be 5% of the He content. The
simulation results show that 70% of the damages occur at Be PKA energies below
100 keV. The number of subcascades increases with the PKA energy and reach
saturation at about 110 keV. The most probable subcascade number is 50 for PKA
energies close to 50 keV. The final Be and He atom distributions confirm that
the most likely defects are single vacancies, divacancies, and small-vacancy
nanovoids containing more than 5 atoms.

#### Acknowledgments

The authors gratefully acknowledge K. Berovski, S. Peneva, and P. Staykov for their participation in discussions and their contribution to the technical preparation of the manuscript. This paper was supported by the EURATOM/ INRNE CSA Contract no. 801499.

#### References

- C. S. Becquart, “RPV steel microstructure evolution under irradiation: a multiscale approach,”
*Nuclear Instruments and Methods in Physics Research B*, vol. 228, no. 1–4, pp. 111–121, 2005. View at Publisher · View at Google Scholar - SRIM-2006, http://www.srim.org.
- C. Inguimbert and R. Gigante, “NEMO: a code to compute NIEL of protons, neutrons, electrons, and heavy ions,”
*IEEE Transactions on Nuclear Science*, vol. 53, no. 4, part 1, pp. 1967–1972, 2006. View at Publisher · View at Google Scholar - S. Shimakawa, N. Nojiri, and N. Sekimura, “NPRIM computer code of neutron dose calculation for radiation damage estimation,” in
*REACTOR DOSIMETRY IN THE 21ST CENTURY Proceedings of the 11th International Symposium on Reactor Dosimetry*, pp. 84–92, World Scientific, River Edge, NJ, USA, 2002. View at Publisher · View at Google Scholar - ENDF/B, http://www-nds.iaea.org/endf.
- D. S. Gelles, G. A. Sernyaevb, M. D. Donnec, and H. Kawamura, “Radiation effects in beryllium used for plasma protection,”
*Journal of Nuclear Materials*, vol. 212–215, part 1, pp. 29–38, 1994. View at Publisher · View at Google Scholar - J. Ziegler, J. P. Biersack, and U. Littmark,
*The Stopping and Range of Ions in Solids*, vol. 1, Pergamon, New York, NY, USA, 1985. - M. T. Robinson, “Computer simulation studies of high-energy collision cascades,”
*Nuclear Instruments and Methods in Physics Research B*, vol. 67, no. 1–4, pp. 396–400, 1992. View at Publisher · View at Google Scholar - Materials Assessment Report G 74 MA 10 01-07-11 W 0.2, Ch 2.1, pp 29.
- R. E. Stoller, “The role of cascade energy and temperature in primary defect formation in iron,”
*Journal of Nuclear Materials*, vol. 276, no. 1–3, pp. 22–32, 2000. View at Publisher · View at Google Scholar - D. J. Bacon, F. Gao, and Y. N. Osetsky, “The primary damage state in fcc, bcc and hcp metals as seen in molecular dynamics simulations,”
*Journal of Nuclear Materials*, vol. 276, no. 1–3, pp. 1–12, 2000. View at Publisher · View at Google Scholar