Hostname: page-component-cd9895bd7-p9bg8 Total loading time: 0 Render date: 2024-12-27T02:35:33.550Z Has data issue: false hasContentIssue false

Molecular-Dynamics Simulation of Displacement Cascades in Zircon (ZrSiO4)

Published online by Cambridge University Press:  10 February 2011

J-P. Crocombette
Affiliation:
DTA/DECM/SRMP C.E. Saclay, 91191 Gif/Yvette Cedex, France, [email protected]
D. Ghaleb
Affiliation:
DCC/DRDD/SCD C.E.A. Valrho-Marcoule, Bp 171-30207 Bagnols/cèze Cedex, France, [email protected]
Get access

Abstract

Zircon (ZrSiO4) is of great interest for the nuclear industry as it is one of the new crystalline waste form considered for the disposal of actinides, for example weapons' plutonium in USA. In this study the effects of displacement cascade due to α-decay has been modelled from an atomistic point of view by molecular dynamics simulation using Born-Mayer-Huggins empirical potential. The numerical values of the parameters of this potential have been fitted on structural equilibrium properties of the crystal and on atomic arrangements.

Displacement cascades are reproduced by accelerating one of the atoms of the cell, thus modelling the effect of the α-decay recoil nucleus. Kinetic energies up to 2 keV have been introduced. The unfolding of the cascades and the final structures have been studied in detail. The centre of the displacement cascade exhibits an amorphous zone where the zircon structure is completely lost. It contains an assembly of distorted SiO4 tetrahedra and disordered zirconium polyhedra. The zirconium ions (originally surrounded by 8 oxygen atoms) exhibit a decrease in their coordination number to 7 or 6 in agreement with what is observed for zirconium ions in amorphous zircon, zirconia or glasses. The size of the amorphous zone and the number of atoms displaced have been estimated for different recoil energies.

The energy stored during the cascade has been calculated. It exhibits an overall good agreement with the available experimental data at complete amorphization.

The existence of an amorphous track in our calculated cascades shows that the correct model for the amorphisation process should take into account the existence of a direct impact amorphous zone.

Type
Research Article
Copyright
Copyright © Materials Research Society 1998

Access options

Get access to the full version of this content by using one of the access options below. (Log in options will check for institutional or personal access. Content may require purchase if you do not have access.)

References

REFERENCES

1 Ewing, R. C., Lutze, W. and Weber, W. J., J. Mater. Res. 10, (1995), 243 Google Scholar
2 Weber, W. J., Ewing, R.C. and Wang, L., J. Mater. Res. 9, 688 (1994)Google Scholar
3 Doan, N. V. and Tietze, H., Nucl. Instr. and Meth. in Phys. Res B 102, 58 (1995)Google Scholar
4 Delaye, J. M. and Ghaleb, D., J. Non-Cryst Solids 195, 239 (1996).Google Scholar
5 Delaye, J. M. and Ghaleb, D., J. Nuclear Mat. 244, 22 (1997)Google Scholar
6 Robinson, K., Gibbs, G. V. and Ribbe, P. H. S., Amer Mineral. 56, 782 (1971)Google Scholar
7 Ozkan, H., Cartz, L., Jamiero, J. C., J. Applied Phys. 45, 556 (1974)Google Scholar
8 Zhu, H., Averback, R.S., Nastasi, M. Phil. Mag. A 71, 735 (1995)Google Scholar
9 Farges, F. Phys. Chem. Minerals 20, 504 (1994)Google Scholar
10 Wyckoff, R., Crystal Structures, Interscience, New York (1951)Google Scholar
11 Petit-Maire, D., PhD thesis, University of Paris VI (1988)Google Scholar
12 Webb, R. and Carter, G., Rad. Effects 42, 159 (1979)Google Scholar
13 Weber, W. J. J. Mater. Res. 5, 2687 (1995)Google Scholar