No Arabic abstract
This paper continues our previous work on a nanostructural evolution model for Fe-C alloys under irradiation, using Object Kinetic Monte Carlo modeling techniques. We here present a number of sensitivity studies of parameters of the model, such as the carbon content in the material, represented by generic traps for point defects, the importance of traps, the size dependence of traps and the effect of the dose rate.
This work extends our Object Kinetic Monte Carlo model for neutron irradiation-induced nanostructure evolution in Fe-C alloys to consider higher irradiation temperatures. The previous study concentrated on irradiation temperatures < 370 K. Here we study the evolution of vacancy and self-interstitial atom (SIA) cluster populations at the operational temperature of light water reactors, by simulating specific reference irradiation experiments. Following our previous study, the effect of carbon on radiation defect evolution can be described in terms of formation of immobile complexes with vacancies, that in turn act as traps for SIA clusters. This dynamics is simulated using generic traps for SIA and vacancy clusters. The traps have a binding energy that depends on the size and type of the clusters and is also chosen on the basis of previously performed atomistic studies. The model had to be adapted to account for the existence of two kinds of SIA clusters, <111> and <100>, as observed in electron microscopy examinations of Fe alloys neutron irradiated at the temperatures of technological interest. The model, which is fully based on physical considerations and only uses a few parameters for calibration, is found to be capable of reproducing the experimental trends, thereby providing insight into the physical mechanisms of importance to determine the type of nanostructural evolution undergone by the material during irradiation.
Neutron irradiation induces in steels nanostructural changes, which are at the origin of the mechanical degradation that these materials experience during operation in nuclear power plants. Some of these effects can be studied by using as model alloy the iron-carbon system. The Object Kinetic Monte Carlo technique has proven capable of simulating in a realistic and quantitatively reliable way a whole irradiation process. We have developed a model for simulating Fe-C systems using a physical description of the properties of vacancy and self-interstitial atom (SIA) clusters, based on a selection of the latest data from atomistic studies and other available experimental and theoretical work from the literature. Based on these data, the effect of carbon on radiation defect evolution has been largely understood in terms of formation of immobile complexes with vacancies that in turn act as traps for SIA clusters. It is found that this effect can be introduced using generic traps for SIA and vacancy clusters, with a binding energy that depends on the size of the clusters, also chosen on the basis on previously performed atomistic studies. The model proved suitable to reproduce the results of low (<350 K) temperature neutron irradiation experiments, as well as the corresponding post-irradiation annealing up to 700 K, in terms of defect cluster densities and size distribution, when compared to available experimental data from the literature. The use of traps proved instrumental for our model.
The solid inner core of the Earth is predominantly composed of iron alloyed with several percent Ni and some lighter elements, Si, S, O, H, and C being the prime candidates. There have been a growing number of papers investigating C and H as possible light elements in the core, but the results are contradictory. Here, using ab initio simulations, we study the Fe-C and Fe-H systems at inner core pressures (330-364 GPa). Using the evolutionary structure prediction algorithm USPEX, we have determined the lowest-enthalpy structures of possible carbides (FeC, Fe2C, Fe3C, Fe4C, FeC2, FeC3, FeC4 and Fe7C3) and hydrides (Fe4H, Fe3H, Fe2H, FeH, FeH2, FeH3, FeH4) and have found that Fe2C (Pnma) is the most stable iron carbide at pressures of the inner core, while FeH, FeH3 and FeH4 are stable iron hydrides at these conditions. For Fe3C, the cementite structure (Pnma) and the Cmcm structure recently found by random sampling are less stable than the I-4 and C2/m structures found here. We found that FeH3 and FeH4 adopt chemically interesting thermodynamically stable structures, in both compounds containing trivalent iron. The density of the inner core can be matched with a reasonable concentration of carbon, 11-15 mol.percent (2.6-3.7 wt.percent) at relevant pressures and temperatures. This concentration matches that in CI carbonaceous chondrites and corresponds to the average atomic mass in the range 49.3-51.0, in close agreement with inferences from the Birchs law for the inner core. Similarly made estimates for the maximum hydrogen content are unrealistically high, 17-22 mol.percent (0.4-0.5 wt.percent), which corresponds to the average atomic mass in the range 43.8-46.5. We conclude that carbon is a better candidate light alloying element than hydrogen.
The manipulation of the antiferromagnetic interlayer coupling in the epitaxial Fe/Cr/Fe(001) trilayer system by moderate 5 keV He ion beam irradiation has been investigated experimentally. It is shown that even for irradiation with very low fluences (10^14 ions/cm^2) a drastic change in strength of the coupling appears. For thin Cr-spacers (below 0.6 - 0.7 nm) the coupling strength decreases with fluence, becoming ferromagnetic for fluences above (2x10^14 ions/cm^2). The effect is connected with the creation of magnetic bridges in the layered system due to atomic exchange events caused by the bombardment. For thicker Cr spacers (0.8 - 1.2 nm) an enhancement of the antiferromagnetic coupling strength is found. A possible explanation of the enhancement effect is given.
We calculate the sink strength of dislocations and toroidal absorbers using Object Kinetic Monte Carlo and compare with the theoretical expressions. We get good agreement for dislocations and loop-shaped absorbers of 3D migrating defects, provided that the volume fraction is low, and fair agreements for dislocations with 1D migrating defects. The master curve for the 3D to 1D transition is well reproduced with loop-shaped absorbers and fairly well with dislocations. We conclude that, on the one hand, the master curve is correct for a wide range of sinks and that, on the other, OKMC techniques inherently take correctly into account the strengths of sinks of any shape, provided that an effective way of appropriately inserting the sinks to be studied can be found.