A computer simulation study of the accommodation and diffusion of He in uranium- and plutonium-doped zircon (ZrSiO4)

We report the results of an atomistic computational study of He accommodation and diffusion in the Pu⁴⁺- and U⁴⁺-doped zircon (ZrSiO₄). The He-cation potentials derived for this work give results of comparable accuracy to DFT calculations. We have calculated the structural features of doped lattices as well as He solution energies in interstitial sites in the perfect and doped zircon and its diffusion in these lattices. The mode of He accommodation in the perfect zircon is influenced mainly by the topological features of the lattice, promoting site preference of He towards accommodation in the interstitial sites present in the middle of c cylinder channels, whereas the presence of Pu⁴⁺ and U⁴⁺ dopants in the zircon lattice significantly affects the energetics of He accommodation and diffusion in the lattice. Doping causes strong local structural distortions, extending to next nearest-neighbour atoms of the dopants to a radius of up to ∼4 Å, in agreement with experimental results. The presence of dopants in the vicinity of He enhances the solubility of He in the lattice compared to the perfect lattice. The mechanism of diffusion is also affected, where the dopants can create a He trap along the most energetically favourable pathway in the (0 0 1) direction, which may slow down the movement of He along the c direction. The dopants also lower the energy barriers by ∼50% in the octahedral sites.     

Li同位素在斜顽辉石、镁铝榴石及镁铝尖晶石晶格内扩散与分馏的理论计算研究

斜顽辉石、镁铝榴石和镁铝尖晶石作为辉石族、石榴石族以及尖晶石族中的重要端元,是地球上地幔主要组成矿物。Li同位素是重要的地幔地球化学示踪剂,其在橄榄石、辉石和石榴石等地幔矿物中的扩散分馏的性质对理解Li同位素作为地幔地球化学示踪剂非常重要。我们通过经典力场经验势方法,从原子尺度上计算研究了不同温压条件下Li同位素在斜顽辉石、镁铝榴石和镁铝尖晶石晶格中分别通过不同的填隙机制和取代空位机制迁移的活化能和其在不同晶格位上的分馏效应。我们发现Li同位素是通过取代空位机制在斜顽辉石、镁铝榴石和镁铝尖晶石中进行迁移扩散。Li同位素在不同晶格位上的分馏作用计算表明,在三种矿物中重同位素7Li会更多地进入晶格填隙位中,而6Li则相对更多进入Mg位。温度是影响这种分馏作用的一个关键因素,相应的结果可用来解释地幔Li同位素组成特征及冷却条件下的同位素分馏等科学问题。     

Geothermochronology based on noble gases: II. Stability of the (U-Th)/He isotope system in zircon

The kinetics of He migration from zircon of variable degree of metamictization was investigated. The migration parameters of He were experimentally determined, the influence of radiation damage and the degree of metamictization on the stability of the (U-Th)/He isotope system was evaluated, the mechanisms of noble gas escape from zircon were investigated, new data on the kinetics of He migration were obtained and compared with previous results for the kinetics of Xe migration from zircon of the same geologic objects. It was shown that He occurs in two energy positions in the zircon lattice: the main position (more than 80% He) with an activation energy of ∼39 kcal/mol and k ₀ = 10¹¹ yr⁻¹ and the second position with an activation energy for migration of 5–10 kcal/mol and k ₀ ∼ 10⁶ yr⁻¹. It was concluded that He migration from the main energy position is better described by a single-jump mechanism. The migration of He from the second energy position is consistent with the diffusion mechanism. It was shown that deviations from the linear dependence in the lnln(He₀/He_t)-1/T coordinates are probably related to the destruction of volume defects containing He atoms at high temperatures (more than 1000°C on the experimental time scale) resulting in a change from the single-jump to diffusion mechanism and the presence of atoms migrating via the diffusion mechanism. It was shown that the peak width in the spectrum of radiogenic He release and the appearance of a second peak also depend on the fraction of atoms migrating in accordance with the diffusion mechanism. It was found that the low activation energy for He release from the second energy position indicates the existence of continuous He loss from the zircon lattice.     

Low-temperature anisotropic diffusion of helium in zircon: Implications for zircon (U-Th)/He thermochronometry

In the last decade the zircon (U-Th)/He (ZHe) thermochronometer has been applied to a variety of geologic problems. Although bulk diffusion coefficients for He in zircon are available from laboratory step-heating experiments, little is known about the diffusion mechanism(s) and their dependence on the crystallographic structure of zircon. Here, we investigate the diffusion of He in perfectly crystalline zircon using atomistic simulation methods that provide insights into the structural pathways of He migration in zircon. Empirical force fields and quantum-mechanical calculations reveal that the energy barriers for He diffusion are strongly dependent on structure. The most favorable pathway for He diffusion is the [0 0 1] direction through the open channels parallel to the c-axis (, activation energy for tracer diffusion of a He atom along [0 0 1]). In contrast, energy barriers are higher in other directions where narrower channels for He diffusion are identified, such as [1 0 0], [1 0 1], and [1 1 0] (ΔE^∗ of 44.8, 101.7, and 421.3 kJ mol⁻¹, respectively). Molecular dynamics simulations are in agreement with these results and provide additional insight in the diffusion mechanisms along different crystallographic directions, as well as the temperature dependence. Below the closure temperature of He in zircon [T_c ∼ 180 °C, Reiners P. W., Spell T. L., Nicolescu S., and Zanetti K. A. (2004) Zircon (U-Th)/He thermochronometry: He diffusion and comparisons with Ar-40/Ar-39 dating. Geochim. Cosmochim. Acta68, 1857-1887], diffusion is anisotropic as He moves preferentially along the [0 0 1] direction, and calculated tracer diffusivities along the two most favorable directions differ by approximately five orders of magnitude (D_[001]/D_[100] ∼ 10⁵, at T = 25 °C). Above this temperature, He atoms start to hop between adjacent [0 0 1] channels, along [1 0 0] and [0 1 0] directions (perpendicular to the c-axis). The diffusion along [1 0 0] and [0 1 0] is thermally activated, such that at higher temperatures, He diffusion in zircon becomes nearly isotropic (D_[001]/D_[100] ∼ 10, at T = 580 °C). These results suggest that the anisotropic nature of He diffusion at temperatures near the closure temperature should be considered in future diffusivity experiments. Furthermore, care should be taken when making geologic interpretations (e.g., exhumation rates, timing of cooling, etc.) from this thermochronometer until the effects of anisotropic diffusion on bulk ages and closure temperature estimates are better quantified.     

Radiation induced lattice defects in natural zircon (ZrSiO4) observed at atomic resolution

Fission tracks and point defects in natural zircon are directly observed by a 1 MV electron microscope at atomic resolution for three types of samples adjusted to the 100 orientation. Lattice planes intersecting the fission tracks at high angles are distorted in a rather irregular manner over a wide region up to more than 100 Å wide. Diameter of the tracks, ranging from 25 Å to 40Å, is much narrower than those so far reported for the U-doped synthetic zircon (100–200 Å), UO₂ thin film (100 Å), mica (66 Å, 240 Å) or fluorophlogopite (150 Å). The fact that fairly long tracks thousands of angstroms in length are observed in thin 100-oriented sample hundreds of angstroms in thickness and that some of them are nearly parallel to a low index lattice plane suggest a possible occurrence of channelling in the process of track formation. Parallel tracks often observed in chemically etched specimens support the idea of channelling. Slightly bent tracks are sometimes observed. It is concluded from computer simulation that many contrast anomalies of bright and dark spots in the lattice image are due to point defects of vacancies and interstitial atoms, mainly produced by the direct atomic collision with α-particles or by passage of ionizing nuclear particles. Optimum conditions of the observation of point defects with highest contrast are studied. One interstitial Zr atom or one Zr ion vacancy will give very low contrast and will be not detectable unless the crystal is less than two unit cells thick. A pair of Zr ion vacancies, however, yields extended detectable limit of thickness. Some of the observed defects are in good accordance with those simulated.     

Zircon (U-Th)/He thermochronometry: He diffusion and comparisons with Ar/Ar dating

(U-Th)/He chronometry of zircon has a wide range of potential applications including thermochronometry, provided the temperature sensitivity (e.g., closure temperature) of the system be accurately constrained. We have examined the characteristics of He loss from zircon in a series of step-heating diffusion experiments, and compared zircon (U-Th)/He ages with other thermochronometric constraints from plutonic rocks. Diffusion experiments on zircons with varying ages and U-Th contents yield Arrhenius relationships which, after about 5% He release, indicate E_a = 163-173 kJ/mol (39-41 kcal/mol), and D₀ = 0.09-1.5 cm²/s, with an average E_a of 169 ± 3.8 kJ/mol (40.4 ± 0.9 kcal/mol) and average D₀ of 0.46^+0.87_−0.30 cm²/s. The experiments also suggest a correspondence between diffusion domain size and grain size. For effective grain radius of 60 μm and cooling rate of 10°C/myr, the diffusion data yield closure temperatures, T_c, of 171-196°C, with an average of 183°C. The early stages of step heating experiments show complications in the form of decreasing apparent diffusivity with successive heating steps, but these are essentially absent in later stages, after about 5-10% He release. These effects are independent of radiation dosage and are also unlikely to be due to intracrystalline He zonation. Regardless of the physical origin, this non-Arrhenius behavior is similar to predictions based on degassing of multiple diffusion domains, with only a small proportion (<2-4%) of gas residing in domains with a lower diffusivity than the bulk zircon crystal. Thus the features of zircon responsible for these non-Arrhenius trends in the early stages of diffusion experiments would have a negligible effect on the bulk thermal sensitivity and closure temperature of a zircon crystal.We have also measured single-grain zircon (U-Th)/He ages and obtained ⁴⁰Ar/³⁹Ar ages for several minerals, including K-feldspar, for a suite of slowly cooled samples with other thermochronologic constraints. Zircon He ages from most samples have 1 σ reproducibilities of about 1-5%, and agree well with K-feldspar ⁴⁰Ar/³⁹Ar multidomain cooling models for sample-specific closure temperatures (170-189°C). One sample has a relatively poor reproducibility of ∼24%, however, and a mean that falls to older ages than predicted by the K-feldspar model. Microimaging shows that trace element zonation of a variety of styles is most pronounced in this sample, which probably leads to poor reproducibility via inaccurate α-ejection corrections. We present preliminary results of a new method for characterizing U-Th zonation in dated grains by laser-ablation, which significantly improves zircon He age accuracy.In summary, the zircon (U-Th)/He thermochronometer has a closure temperature of 170-190°C for typical plutonic cooling rates and crystal sizes, it is not significantly affected by radiation damage except in relatively rare cases of high radiation dosage with long-term low-temperature histories, and most ages agree well with constraints provided by K-spar ⁴⁰Ar/³⁹Ar cooling models. In some cases, intracrystalline U-Th zonation can result in inaccurate ages, but depth-profiling characterization of zonation in dated grains can significantly improve accuracy and precision of single-grain ages.     

He diffusion systematics in minerals: Evidence from synthetic monazite and zircon structure phosphates

Helium diffusivity was measured in synthetic rare-earth-element orthophosphates with systematically varying properties to evaluate potential controls on He transport in minerals. In the zircon structure phosphates (in this study, the phosphates of Tb, Dy, Ho, Er, Tm, Yb, and Lu as well as synthetic xenotime, YPO₄), He diffusion is strongly anisotropic. Transport apparently proceeds preferentially through channels aligned with the c-axis. The activation energy for diffusion is almost the same (122 ± 6 kJ/mol) in all members of this family, but there is a monotonic decrease in D_o with atomic number from TbPO₄ (∼10⁵ cm²/s) to LuPO₄ (∼10 cm²/s). The c-parallel channels become increasingly constricted in the same sequence, likely accounting for the systematically decreasing diffusivity. The He closure temperature (r = 1 cm, dT/dt = 10 °C/Myr) increases with atomic number from 44 °C for TbPO₄ to 88 °C for LuPO₄. Diffusion of radiogenic helium from natural zircon and xenotime is much slower than these synthetic analogs predict, suggesting that coupled substitution of REE and P for Zr and Si and/or radiation damage profoundly modify the energetics of interstitial He diffusion. In particular, α-recoil may play a key role by damaging the continuity and integrity of the channels.Monazite structure phosphates (here La, Ce, Pr, Nd, Sm, and Gd phosphate) are far more He retentive than those of the zircon structure. Activation energies increase smoothly with atomic number from LaPO₄ (183 kJ/mol) to NdPO₄ (224 kJ/mol) then decrease again to GdPO₄ (198 kJ/mol). D_o values mimic this pattern, spanning a range from ∼10⁻¹ cm²/s (GdPO₄) to 10⁴ cm²/s (NdPO₄). Nevertheless, He closure temperatures increase monotonically with atomic number, from 300 °C in LaPO₄ to 410 °C in GdPO₄. No evidence was obtained bearing on diffusion anisotropy, but the monazite structure lacks through-going channels so it is not expected. Diffusion parameters for radiogenic helium in natural monazite are similar to those obtained on the synthetic analogs.Ionic porosity is not the primary control on He diffusion in the orthophosphates. Within a given structure and with limited elemental substitution, ionic porosity and He closure temperature are negatively correlated, as predicted. However, differences between crystal structures are far more important than ion packing density: at comparable ionic porosity the monazite structure phosphates have He closure temperatures ∼300 °C higher than the xenotime structure phosphates. Modifications to the structures by radiation damage likely play a similarly significant role in controlling He diffusion.     

Carbon in solid solution in forsterite—a key to the untractable nature of reduced carbon in terrestrial and cosmogenic rocks

Carbon is an incompatible element in oxide and silicate lattices. Until now it has been believed to exist only in the form of CO^2?₃ ions, molecularily dissolved CO₂ or graphitic inclusions. Recently it has been shown that carbon can dissolve in refractory oxides, like MgO and CaO, in the form of carbon atoms.The experimental results obtained with carbonaceous MgO are reviewed and new results are presented which demonstrate that synthetic forsterite and natural olivines can also take up atomic carbon in solid solution. The incorporation of the carbon atoms is treated thermodynamically. Near the melting point they probably occupy cation vacancies, but with decreasing temperature they are progressively transferred on interstitial sites. On these sites they are very mobile and tend to segregate into the elastically relaxed subsurface zone, but exsolution to graphite is prevented by the strain fields surrounding each carbon atom.Upon heating, however, the atomic carbon may react with lattice oxygen to give CO₂ and with co-dissolved hydrogen to give a wide variety of hydrocarbons. The underlying reaction mechanisms, involving the formation and decay of O^? ions, are discussed in view of the so-called ‘carbonatic carbon’ and ‘reduced carbon’ in magmatic minerals and meteorites, in view of the diamond genesis and also in view of the reversible CO₂ solubility in silicate magmas at high pressures and temperatures.     

Neutron irradiation and post-irradiation annealing of rutile (TiO2−x ): effect on hydrogen incorporation and optical absorption

Neutron irradiation and post-irradiation annealing under oxidising and reducing conditions have been used to investigate H incorporation in, and the optical properties of, reduced (TiO_2−x ) rutile. Optical absorption in rutile is mainly due to a Ti³⁺ Ti⁴⁺ intervalence charge transfer effect. The main mechanism for H incorporation in rutile involves interstitial H not coupled to other defects, which has important implications for the rate of H diffusion, and possibly also on the electrical properties of rutile. Additional minor OH absorption bands in IR spectra indicate that a small amount of interstitial H is coupled to defects such as Ti³⁺ on the main octahedral site, and indicates that more than one H incorporation mechanism may operate. Concentration of oxygen vacancies has a controlling influence on the H affinity of rutile.     

NMR applied to minerals. V. The localization of vacancies in the octahedral sheet of aluminous biotites

The localization of vacancies in the octahedral sheet of aluminous biotites has been investigated by the use of ¹H nuclear magnetic resonance (NMR) and infrared (IR) spectroscopy. The polytype to which the samples belonged was determined by x-ray diffraction methods. The joint use of all three techniques was essential in determining unambiguously the exact position of the vacancies and orientation of the associated O-H and Fe-H vectors. It was found that vacancies are located in the M1 (pseudo-centrosymmetric) sites and that contrary to what is usually assumed the OH dipoles are not, in this case, oriented directly towards the vacancy site.     

Atomic ordering around the oxygen vacancies in sillimanite

Computer simulation is used to investigate the short range ordering around an isolated oxygen vacancy in sillimanite. The static lattice energy with the use of empirical potentials is calculated, for different Al/Si distributions around a vacancy in a supercell of sillimanite. A parametrisation of the total energy is built up and used to deduce the best Al/Si ordering around the oxygen vacancy. It is found that a strong ordering about the vacancy occurs. In the ab-plane two sets of aluminium cluster are found besides the vacancy, surrounded above and below by silicon atoms, a configuration that promotes local charge balance. By placing two vacancies on sites directly adjacent to the same oxygen, the central oxygen site is bonded to four cations: this situation is found to be energetically unfavourable.     

Kinetics of zircon dissolution and zirconium diffusion in granitic melts of variable water content

The experimental dissolution of zircon into a zircon-undersaturated felsic melt of variable water content at high pressure in the temperature range 1,020° to 1,500° C provides information related to 1) the solubility of zircon, 2) the diffusion kinetics of Zr in an obsidian melt, and 3) the rate of zircon dissolution. Zirconium concentration profiles observed by electron microprobe in the obsidian glass adjacent to a large, polished zircon face provide sufficient information to calculate model diffusion coefficients. Results of dissolution experiments conducted in the virtual absence of water (<0.2% H₂O) yield an activation energy (E) for Zr transport in a melt ofM=1.3 [whereM is the cation ratio (Na+K+2Ca)/(Al·Si)] of 97.7±2.8 kcal-mol^?1, and a frequency factor (D ₀) of 980 _?580 ^+1,390 cm²-sec^?1. Hydrothermal experiments provide an E=47.3±1.9 kcal-mol^?1 andD ₀=0.030 _?0.015 ^+0.030 cm²-sec^?1. Both of these results plot close to a previously defined diffusion compensation line for cations in obsidian. The diffusivity of Zr at 1,200° C increases by a factor of 100 over the first 2% of water introduced into the melt, but subsequently rises by only a factor of five to an apparent plateau value of ～2×10^?9 cm²-sec^?1 by ～6% total water content. The remarkable contrast between the wet and dry diffusivities, which limits the rate of zircon dissolution into granitic melt, indicates that a 50 μm diameter zircon crystal would dissolve in a 3 to 6% water-bearing melt at 750° C in about 100 years, but would require in excess of 200 Ma to dissolve in an equivalent dry system. From this calculation we conclude that zircon dissolution proceeds geologically instantaneously in an undersaturated, water-bearing granite. Estimates of zircon solubility in the obsidian melt in the temperature range of 1,020° C to 1,500° C confirm and extend an existing model of zircon solubility to these higher temperatures in hydrous melts. However, this model does not well describe zircon saturation behavior in systems with less than about 2% water.     

Water incorporation in synthetic and natural MgAl2O4 spinel

The solubility and incorporation mechanisms of water in synthetic and natural MgAl₂O₄ spinel have been investigated in a series of high-pressure/temperature annealing experiments. In contrast to most other nominally anhydrous minerals, natural spinel appears to be completely anhydrous. On the other hand, non-stoichiometric Al-rich synthetic (defect) spinel can accommodate several hundred ppm water in the form of structurally-incorporated hydrogen. Infrared (IR) spectra of hydrated defect spinel contain one main O-H stretching band at 3343-3352 cm⁻¹ and a doublet consisting of two distinct O-H bands at 3505-3517 cm⁻¹ and 3557-3566 cm⁻¹. IR spectra and structural refinements based on single-crystal X-ray data are consistent with hydrogen incorporation in defect spinel onto both octahedral and tetrahedral O-O edges. Fine structure of O-H bands in IR spectra can be explained by partial coupling of interstitial hydrogen with cation vacancies, or by the effects of Mg-Al disorder on the tetrahedral site. The concentration of cation vacancies in defect spinel is a major control on hydrogen affinity. The commercial availability of large single crystals of defect spinel coupled with high water solubility and similarities in water incorporation mechanisms between hydrous defect spinel and hydrous ringwoodite (Mg₂SiO₄) suggests that synthetic defect spinel may be a useful low-pressure analogue material for investigating the causes and consequences of water incorporation in the lower part of Earth’s mantle transition zone.     

On the mechanisms for H and Al incorporation in stishovite

The solubility and incorporation mechanisms of hydrogen in synthetic stishovite as a function of Al₂O₃ content have been investigated. Mechanisms for H incorporation in stishovite are more complex than previously thought. Most H in stishovite is incorporated via the Smyth et al. (Am Mineral 80:454–456, 1995) model, where H docks close to one of the shared O–O edges, giving rise to an OH stretching band in infrared (IR) spectra at 3,111–3,117 cm⁻¹. However, careful examination of IR spectra from Al-stishovite reveals the presence of an additional OH band at 3,157–3,170 cm⁻¹. All H is present on one site, with interstitial H both coupled to Al³⁺ substitutional defects on adjacent octahedral (Si⁴⁺) sites, and decoupled from other defects, giving rise to two distinct absorption bands. Trends in IR data as a function of composition are consistent with a change in Al incorporation mechanism in stishovite, with Al³⁺ substitution for Si⁴⁺ charge-balanced by oxygen vacancies at low bulk Al₂O₃ contents, and coupled substitution of Al³⁺ onto octahedral (Si⁴⁺) and interstitial sites at high bulk Al₂O₃ contents. Trends in OH stretching frequencies as a function of Al₂O₃ content suggest that any such change in Al incorporation mechanism could alter the effect that Al incorporation has on the compressibility of stishovite, as noted by Ono et al. (Am Mineral 87:1486–1489, 2002).     

First-principles study of OH defects in zircon

The infrared spectroscopic properties of selected OH defects in zircon are investigated by first-principles calculations. The explicit treatment of the coupled nature of OH motions in the stretching modes, together with the calculation of the intensity and polarization of absorption bands, makes it possible to directly compare theoretical and experimental data. The bands observed at 3,420 cm^?1 (polarization parallel to c axis) and 3,385 cm^?1 (polarization perpendicular to c axis) in natural and synthetic samples correspond to the IR-active vibrational modes of the hydrozircon defect, that is, fully protonated Si vacancy. The broad band observed at 3,515 cm^?1 in the spectrum of zircon crystals grown in F-rich environments is consistent with the occurrence of composite (OH,F) tetrahedral defects. Calculations also show that the band observed at 3,200 cm^?1 in the spectrum of synthetic undoped samples can be ascribed to fully protonated Zr vacancies. The theoretical values of integrated absorption coefficients indicate that general correlations can be reasonably used to determine the concentration of OH groups in zircon.     

Crystal chemistry of Sc-bearing synthetic diopsides

Scandium substitution in the diopside structure was studied by single-crystal X-ray diffraction on a series of synthetic diopside samples. These diopsides were doped with increasing amounts of Sc³⁺ through a coupled substitution involving ^M1(Sc³⁺)₁ ^M2[] ₁ ^T (B)₁ ^M1(Mg²⁺)_?1 ^M2(Ca²⁺)_?1 ^T(Si⁴⁺)_?1 exchange, whereby charge compensation is achieved by vacancies at the M2 sites and B at the tetrahedral sites. The substitution of scandium for magnesium at the M1 site results in an increase in volume and distortion of the M1O₆ polyhedron. The accompanying creation of vacancies at the M2 sites causes an increase in the M2 polyhedral volume. The modifications of the M1 and M2 polyhedra result in an increase of the polyhedral strip along the b lattice direction and a straightening of the tetrahedral chain. The geometrical modification of the M1 polyhedron due to scandium incorporation is comparable to those observed when similar amounts of Ti³⁺/Ti⁴⁺ substitute for Mg in the diopside structure, suggesting a structural control on the solubility of Sc and Ti in diopside that may influence the extent of the solid solutions between the Sc and Ti end-members.     

Application of lattice imagery to radiation damage investigation in natural zircon

High-, intermediate-, and low-type zircon crystals of natural origin were investigated using a 1,000 kV high-resolution electron microscope. The lattice images obtained successfully for high zircon were in good accordance with computer simulated ones, and 1.5 Å separations, the nearest distance between zirconium atoms projected along the a axis, were clearly resolved under a certain instrumental condition. The images of fission tracks and surrounding areas show nearly perfect lattice images and that within the fission tracks, with a width of 20 ～ 30 Å and length of ca. 1,000 Å, the structure is heavily disordered, almost amorphous; that both sides of the tracks the lattices are displaced or dislocated, and that in the area adjacent to the tracks, bright and dark spot images occur, corresponding to vacancies and their interstitial atoms. In low zircon, the structure is completely destroyed to show an entirely amorphous state, whereas an intermediate type consists of domains of the order of 50 ～ 100 Å across with nearly regular lattices, along whose boundaries strongly disordered areas with widths of few tens of angstroms appear, but the relative orientations of the neighbouring domains are almost continuous. Thus a whole process of metamictization is visualized on a lattice scale. Metamictization proceeds principally by the formation of fission tracks, the direct result of fast movement of nuclear particles; recoil nuclei therefrom seem to play a less important role in the destruction of the structure.     

Atomistic simulation of the mechanisms of noble gas incorporation in minerals

Atomistic simulations have been carried out to investigate the mechanisms of noble gas incorporation in minerals using both the traditional two-region approach and the “supercell” method. The traditional two-region approach has been used to calculate defect energies for Ne, Ar, Kr and Xe incorporation in MgO, CaO, diopside and forsterite in the static limit and at one atmosphere pressure. The possibilities of noble gas incorporation via both substitution and interstitial mechanisms are studied. The favored mechanism varies from mineral to mineral and from noble gas to noble gas. In all minerals studied, the variation of the solution energies of noble gas substitution with atomic radius appears approximately parabolic, analogous to those for 1⁺, 2⁺, 3⁺ and 4⁺ trace element incorporation on crystal lattice sites. Noble gas solution energies thus also fall on a curve, similar to those previously observed for cations with different charges, but with much lower curvature.The “supercell” method has been used to investigate the pressure dependence of noble gas incorporation in the same systems. Results indicate a large variation of the solubility of the larger noble gases, Kr and Xe with pressure. In addition, explicit simulation of incorporation at the (0 0 1) surface of MgO shows that the solubility of the heavier noble gases may be considerably enhanced by the presence of interfaces.     

Changes in the structural and magnetic properties of Ni–substituted hematite prepared from metal oxinates

Samples prepared by the novel method based on the thermal decomposition at 700 °C in air of mixed Fe and Ni oxinates were characterized by thermogravimetric analysis, X-ray diffraction, magnetization measurements, and variable-temperature Mössbauer spectroscopy. It is found that the combustion treatment produces Ni-hematite and trevorite, the fraction of the latter increasing with the increment of Ni in the metal oxinates. Results indicate that the substitution of Ni²⁺ for Fe³⁺ in the hematite structure, which was found to be less than 5.3 mol%, is accompanied by the presence of oxygen and structural vacancy sites. Both the metal replacement and the presence of defects cause the a and c cell hematite parameters to decrease. The Néel temperature and the difference between the saturation fields for the antiferromagnetic AF and the weakly ferromagnetic WF phases are also found to decrease with nickel content. These effects are due to the fact that the magnetic behavior of hematite depends on the presence of Ni, vacancy sites, and lattice distortion. The large decrease in the Néel temperature with Ni doping suggests that structural vacancies are also present On the contrary, Ni incorporation does not appreciably affect the Morin temperature and the temperature range in which both AF and WF phases coexist. It is suggested that these effects are probably due to the differing effects that Ni and defects can produce on the magnetic dipolar and the single ion anisotropies.     

Contribution of interstitial OH groups to the incorporation of water in forsterite

Water incorporation in forsterite samples synthesized under low to medium silica-activity conditions mostly occurs via a substitutional mechanism in which a Si vacancy is compensated by four protons. Corresponding IR absorption spectra display a cluster of narrow and weakly anharmonic OH-stretching bands at wavenumbers above 3,500 cm^?1. However, this diagnostic spectrum is often superimposed to one broader absorption band, rarely two, displaying pronounced temperature-dependent properties and tentatively assigned to H atoms in interstitial position (Ingrin et al. in Phys Chem Miner 40:499–510, 2013). Here, we investigate the structural and vibrational properties of selected interstitial H-bearing defects in forsterite using a first-principles modeling approach. We show that the broad bands discussed by Ingrin et al. (Phys Chem Miner 40:499–510, 2013) are most likely related to interstitial OH groups in the vacant octahedral sites alternating with the M2 sites along the c axis of the forsterite structure. The corresponding OH defects lead to the formation of fivefold coordinated Si species. Their peculiar thermal properties stem from the vibrational phase relaxation due to the anharmonic coupling of the high-energy local OH-stretching mode with a low-energy vibrational mode. This “exchange mode” corresponds to the hindered longitudinal translation of the OH group. These results suggest that at high pressure, hydrogen incorporation in forsterite is dominated by coexisting interstitial OH groups and (4H)_Si defects.