Next Article in Journal
X-ray and Neutron Study on the Structure of Hydrous SiO2 Glass up to 10 GPa
Next Article in Special Issue
Ultrafast Elemental Mapping of Platinum Group Elements and Mineral Identification in Platinum-Palladium Ore Using Laser Induced Breakdown Spectroscopy
Previous Article in Journal
Surface Subsidence Prognosis above an Underground Longwall Excavation and Based on 3D Point Cloud Analysis
Previous Article in Special Issue
Luminescence of Agrellite Specimen from the Kipawa River Locality
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

The In-Situ Quantification of Structural Radiation Damage in Zircon Using Laser-Induced Confocal Photoluminescence Spectroscopy

by
Christoph Lenz
1,2,3,*,
Elena Belousova
3 and
Gregory R. Lumpkin
1
1
Australian Nuclear Science and Technology Organisation, Lucas Heights, Sydney, NSW 2234, Australia
2
Institut für Mineralogie und Kristallographie, Universität Wien, 1090 Vienna, Austria
3
Australian Research Council (ARC) Centre of Excellence for Core to Crust Fluid Systems (CCFS) and National Key Centre for Geochemical Evolution and Metallogeny of Continents (GEMOC), Department of Earth and Environmental Sciences, Macquarie University, Sydney, NSW 2109, Australia
*
Author to whom correspondence should be addressed.
Minerals 2020, 10(1), 83; https://doi.org/10.3390/min10010083
Submission received: 11 December 2019 / Revised: 6 January 2020 / Accepted: 8 January 2020 / Published: 20 January 2020
(This article belongs to the Special Issue Modern Luminescence Spectroscopy of Minerals)

Abstract

:
We present a new methodology for laser-induced steady-state photoluminescence (PL) spectroscopy of Dy3+ that aims at a direct quantification of the amorphous fraction fa present in zircon (ZrSiO4), which undergoes a transition from a crystalline to a metamict state due to cumulative self-irradiation damage caused by the radioactive decay of substituted U and Th. Using state-of-the-art confocal spectrometers attached to optical microscopes, measurements may be performed non-destructively on the micrometre length-scale with the option to visualize radiation-damage patterns as revealed by hyperspectral PL maps. Zircon from the Ratnapura district (Sri Lanka, ~520 Ma), was used as reference material to substantiate the applicability of the proposed method. The accumulation of radiation damage in this material was investigated in detail and obtained fa values correlate with calculated α-doses in accordance to the direct impact model reported variously in the literature. The impact of chemically-induced, heterogeneous broadening of Raman and Dy3+ emission spectral bands is discussed on two examples from Mt. Malosa district, Malawi. A mean weighted U-Pb isotope age of 111 ± 1 Ma (pegmatitic-type) and a discordia age of 112 ± 1.6 Ma (hydrothermal-type) as obtained by LA-ICP-MS confirm their close genetic and temporal relationship. Studied zircon examples demonstrate that the amount of radiation damage present may have a substantial effect on the precision of LA-ICP-MS ages, but cannot be considered an exclusive cause for bias of obtained isotope ages.

1. Introduction

The transformation of an initially crystalline material into an amorphous state due to corpuscular irradiation is commonly referred to as “amorphisation” or “metamictisation” [1,2,3]. In nature, the accumulation of self-irradiation induced structural damage is created mainly during α-decay events in the 232Th, 235U, and 238U decay chains, by the nuclear interaction (atomic “knock-ons”) of recoiled heavy daughter [4,5]. Many accessory minerals, including zircon (ZrSiO4), incorporate variable amounts of actinides, whose radioactive decay creates structural disorder in their crystal structure over geological periods of time. In zircon this causes changes of its physical and solid-state properties such as refraction and birefringence, hardness, density, and elastic moduli [6,7,8,9], and a general decrease of the chemical resistance [10,11,12]. The generally increased susceptibility of radiation-damaged zircon to chemical alteration is of enormous importance, as these processes may negatively affect the ability of the mineral to retain radiogenic nuclei and therefore bias results of chemical and isotopic age determination applied in U-Pb geochronology [13,14,15]. Moreover, the investigation of gradually radiation damaged and metamict minerals, and their synthetic analogues, has increased appreciably over the past decades, stimulated by the potential use of mineral-like ceramics as waste forms for the immobilisation of reprocessed spent nuclear fuel and other radioactive waste [16,17]. Information obtained from studies on radiation-damage-assisted alteration in accessory minerals has important implications for the validation of the long-term performance of nuclear waste forms for the disposal in geological repositories, as radiation damage negatively affects the ability of analogue ceramics to immobilise incorporated radionuclides, e.g., Pu, Am, U and Th [18,19].
One of the major problems in the evaluation of effects caused by radiation-damage is the accurate quantification of the latter; i.e., the estimation of the fraction of the crystal structure which has been amorphised by self-irradiation. Since the early study by Holland and Gottfried [6], changes in zircon birefringence has been used to visualize and estimate the radiation damage in complexly zoned zircon [8,20,21,22,23]. Electron microscopic techniques, such as scanning electron microscopy and transmission electron microscopy are applied frequently to visualize radiation damage in various materials [24,25,26,27,28,29,30,31,32]. Those techniques, however, are mostly of qualitative nature and not appropriate to quantify the amorphous component. For the latter, different X-ray diffraction methods [4,26,33,34] and several spectroscopic techniques have been developed. For instance, Zhang and Salje [35] used infrared-absorption spectroscopy, and Farnan and Salje [36] applied 29Si nuclear magnetic-resonance spectroscopy to directly determine the amorphous fraction fa in radiation-damaged zircon. The major disadvantage of the latter techniques is that all of them probe bulk material, i.e., a few milligram of powdered or entire monolithic zircon crystals have been used, respectively. The samples to be analysed, hence, need to meet certain pre-conditions which include: (i) a sufficient amount of material available and (ii) an adequate chemical homogeneity to ensure that radiation damage has been accumulated homogeneously within the sample over geological periods of time. This is the major reason why those studies took advantage of relatively large, gem-quality zircon samples from gem-gravels in Ratnapura, Sri Lanka that have high chemical homogeneity. For the same reason, zircon samples from this location are thoroughly studied also with respect to their use as reference material in U-Pb geochronology (e.g., SL7: Kinny et al. [37]; CZ3: Pidgeon et al. [38]; M257: Nasdala et al. [39]; M127: Nasdala et al. [40]; BR266: Stern [41]; BB: Santos et al. [42]).
To estimate the degree of radiation-induced structural disorder in complexly zoned zircon on the µm-scale, Raman spectroscopy has been introduced by Nasdala et al. [22,43] using the band width of the symmetric stretching vibration B1g3) of zircon’s SiO4 tetrahedrally-coordinated anionic complex. Confocal photoluminescence (PL) spectroscopy of rare-earth elements (REE) has been introduced recently for the same purpose [44,45]. Trivalent REEs are common trace elements in zircon and many other accessory minerals, and are found in close geochemical relationship to the actinides U and Th due to their similar ionic radii and high valences [46]. The luminescence of REEs is excited most efficiently using appropriate light or laser sources [47,48,49,50,51] and is further known to have unique spectroscopic characteristics due to their shielded electronic intra ff transitions (see detailed discussion in Lenz et al. [52]). Emissions arising from those transitions consist of several narrow sublevel-bands (“Stark levels”) that strongly depend on the local structural environment the specific REE ion is incorporated in, and hence, may be used as site-specific structural probe. Very similar to the Raman spectroscopic approach, the band widths of narrow PL bands of REEs have been used as a measure of irradiation-induced structural disorder (i.e., Dy3+ and Nd3+ in zircon; cf., [44]). A major challenge for using linewidths of Raman or PL signals for a quantitative measure of radiation effects is its calibration using reference samples of known amorphous fraction fa or equivalent effective α-dose (for zircon see Palenik et al. [8]). Calibrations based on the α-dose, which is the theoretical number of α-events calculated from the concentration of U and Th, their half-life and the geological age of the sample [6], are often biased because of the insufficient knowledge of the samples thermal and, hence, annealing history (see the detailed discussion in Nasdala et al. [22]).
In addition, several Raman spectroscopic studies on zircon have revealed that Raman band widths occasionally may be affected strongly by the substitution and presence of extensive amounts of non-formula elements. Such examples comprise high-U zircon of the “Chernobyl lava” [53] and other U-rich samples [54,55], as well as hydrothermal zircon from the Mt. Malosa Massif (Malawi), that is characterised by extremely high concentrations of REEs [56]. This chemically-induced effect may be critical only for a few examples of zircon that grew in late-stage magmatic to hydrothermal environments, but have been documented to be important for other mineral systems that show higher chemical flexibility than zircon (e.g., monazite-type group minerals as reported by Ruschel et al. [57]).
In this study, we tested a new approach presented by Lenz et al. [58] that aims at a direct determination of the amorphous fraction fa from in-situ confocal PL measurements without the need for an indirect correlation via spectral band-widths. This methodology has been developed for artificially heavy-ion irradiated orthophosphate ceramics that showed a rather complex, non-linear band-width broadening upon damage accumulation [58]. Here, we demonstrate the applicability of this approach to naturally self-irradiation damaged zircon using well-studied reference material from the Ratnapura district (Sri Lanka). The basic principle and underlying interpretation of spectral PL data for a visualization of the amorphous fraction fa in hyperspectral maps are demonstrated in detail.

2. Materials and Methods

2.1. Samples

Previous studies on zircon samples from Sri Lanka have investigated the correlation between the accumulated amorphous fraction fa versus the calculated α-dose in very detail [33,35,36]. Consequently, we tested our PL methodology with samples from this location. Various gem-quality, homogeneous single crystals (see Table A1) and two brownish, non-transparent, large, cm-sized, but heterogeneously zoned, single crystals were investigated in detail (SL1–minID: P3H-LYV; SL2–minID: 4TL-XKA; see Table A2). The Sri Lankan zircon samples stem from well-known gemstone placers located in the central Sri Lankan granulite belt in the Ratnapura district. Zircon from the Ratnapura district are discussed to originate from granitoid pegmatites that have intruded into Al-rich sediments, but the original location of its source rocks is unknown [25].
In this study, we also included large, heterogeneous zircon crystals from a granitic, alkaline pegmatite near Mt. Malosa, located within the Zomba–Malosa Complex of the Chilwa Alkaline Province in Malawi, Africa. Pegmatites of this region have undergone multiple late-stage hydrothermal alteration and are characterised by unique mineral composition including aegirine, arfvedsonite, fluorite, Ca-REE-fluorocarbonates, Nb-Ta-Y oxides in large miarolitic cavities with Na–Be-silicates, xenotime-(Y), as well as zircon [56,59]. Three morphological types of zircon from this location have been distinguished [56]: (1) translucent reddish-brownish crystals, (2) opaque short prismatic brownish zircon and (3) translucent prismatic yellow crystals. We have studied zircon samples of type 2 (sample MalawiX–minID: GJ0-32C) and type 3 (MalawiA–minID: 8DE-AE0).
Following the approach of Nasdala et al. [60], the studied large, heterogeneous single-crystals were cut in half using a diamond-tipped stainless-steel wire saw (WELL Diamond Saw 3032; Ø 300 µm effective cutting). One half was subjected to heat treatment in furnace at ~1300°C for 3 days, whereas the other half remained untreated. Both halves were then embedded in one epoxy mount in close proximity, and gently polished. This preparation procedure was chosen to obtain sample mounts for direct comparison of two halves of one initial crystal which are widely similar in terms of their internal textures and chemical compositions, but significantly different in terms of their structural state. Note that one half is still in its natural, potentially radiation-damaged state, whereas its counterpart has undergone complete structural reconstitution.

2.2. Photoluminescence and Raman Spectroscopy

A confocal, dispersive LabRam Evolution HR800 spectrometer (Horiba, Kyoto, Japan) was used for Raman and laser-induced steady-state photoluminescence (PL) spectroscopic measurements. The spectrometer is equipped with a Si-based, Peltier-cooled charge-coupled device (CCD) detector and coupled to an Olympus BX41 optical microscope with a confocal hole in the analysing beam path achieving a maximum spatial resolution of ~1 µm lateral and ~2–3 µm in depth using a 100× objective (Olympus MplanN, NA = 0.9). We used a 473 nm diode-pumped solid-state laser with the beam path aligned through the microscope objective (quasi-backscattering configuration) to directly excite the most prominent luminescence emission in zircon (e.g., Dy3+, Figure 1a; see details in [44]). Raman measurements were performed using the 633 nm emission of a He-Ne gas laser. Laser power was adjusted to not exceed ca. 8–15 mW at the sample surface. PL measurements were performed using a 600 lines per millimetre (l/mm) grating that achieves a spectral resolution of ca. 2.5 cm−1 (in the blue to yellow spectral range of the main Dy emissions). High-resolution Raman spectra were recorded using a 1800 l/mm grating with a resolution of 0.8 cm−1 in the red spectral range. Spectral calibration of the spectrometer system was done using the 633 nm He-Ne laser Rayleigh peak, and its accuracy cross-checked regularly using neon-krypton lamp emissions in the respective wavelength ranges. Peak positions are reproducible within ±1.5 cm−1.
Raman bands were background-corrected and fitted using Gauss-Lorentz (pseudo-Voigt) functions. Note that for some samples (e.g., sample MalawiA) a two-fold Gauss-Lorentz function was used to account for strongly asymmetric band shapes. Obtained full-widths at half-maxima (FWHM) were corrected for the system’s spectral resolution according to Váczi [61].
Photoluminescence hyperspectral maps were built from raster of point measurements obtained with the help of a coupled, software-controlled Märzhauser x-y-z mechanical stage. Due to cm-sized dimensions of samples studied here, a somewhat lower spatial resolution of a long-distance 50× objective (Olympus LMPlanFLN, NA = 0.5) was accepted for comparably large hyperspectral maps. Data reduction and treatment for each measurement spot was done using a software-based, automated procedure that is described in more detail below (see Results and Discussion).

2.3. Trace-Element Analyses and U-Pb Geochronology

Trace-element (TE) concentrations and isotopic ratios for U-Pb age determination of various measurement spots on the large zircon crystals SL1, SL2, MalawiX and MalawiA were obtained using an Agilent 7700 quadrupole, inductively-coupled plasma mass spectrometer (ICP-MS) instrument (Agilent Technologies Inc., Santa Clara, CA, USA) equipped with a 193 nm Photon Machines excimer Laser Ablation (LA) system (Photon Machines Inc./Teledyne Cetac, Omaha, NE, USA) with 5 Hz repetition rate and an energy of around 0.06 µJ to 8J/cm2 (summary of data in Table A2, Table A3, Table A4, Table A5, Table A6, Table A7 and Table A8). A very fast scanning data acquisition protocol was employed to minimise signal noise. Data acquisition for each analysis took three minutes (one minute on background, two minutes on signal). Ablation was carried out in He gas to improve sample transport efficiency and to provide stable signals and reproducible Pb/U fractionation. Provided that constant ablation conditions are maintained, accurate correction for U/Pb fractionation can then be achieved using an isotopically homogeneous zircon standard. Sample analyses performed in this study were run along with analyses of the Gemoc GJ-1 zircon standard [62]. This standard is slightly discordant, and has a TIMS 207Pb/206Pb age of 608.5 Ma [63]. The other well-characterised zircon standard 91, 500 and Mud Tank were analysed within the measurement runs as an independent control on reproducibility and instrument stability. U-Pb ages were calculated from raw signal data using the software package Glitter [64]. Individual time-resolved data analysis provides isotopically homogeneous segments of the signal to be selected for integration. We correct the integrated ratios for ablation related fractionation and instrumental mass bias by calibration of each selected time segment against the identical time segments for the standard zircon analyses.
Quantitative results for trace elements reported here were obtained through calibration of relative element sensitivities using the NIST-610 standard glass as external calibration standard. The BCR-2g basaltic glass standard and zircon GJ-1 reference material were analysed along the measurement runs as independent control on reproducibility and instrument stability. Theoretical Zr content was used for internal calibration for unknown zircon samples. The precision and accuracy of the NIST-610 analyses are 1–2% for REE, Y, Sr, Nb, Hf, Ta, Th and U at the ppm concentration level, and from 5% to 10% for Ca, P and Ti (further details see in [65]).
Note, that the analysed concentration of trace dysprosium (Dy) that is considered relevant for its PL emission efficiency is low in the Sri Lankan samples (SL2: 6 to 41 ppm Dy; see Table A9) and comparably high in the Malawian samples (655 to 2123 ppm Dy; see Table A7 and Table A8).

3. Results and Discussion

3.1. Zircon’s Amorphous Fraction Decoded from Dy photoluminescence

The impact of self-irradiation-induced structural disorder on the PL emissions of REEs in zircon is exemplified in Figure 1. Excitation of zircon with a blue laser (λexc = 473 nm) typically results in emission of multiple trace REEs substituted in the crystal structure with that of Dy3+ and Sm3+ in the green-yellow spectral range being most intensive (Figure 1a; [48,50,51]). Figure 1 shows PL spectra from a Cambodian zircon sample from the Ratanakiri district that is characterised by very low degree of structural radiation-damage accumulation as detected with Raman spectroscopy (“Rata”, see details in [66]). Further, spectra of mildly (sample M257, and M146) and severely radiation damaged zircon (sample G3, and N17) of the Ratnapura district, Sri Lanka are shown (see description of samples above, and [44]). Crystal-field split sub-level bands of the Dy3+ (4F9/26H13/2) emission may be easily identified. These narrow bands increasingly broaden and decrease in intensity with progressive accumulation of radiation damage (Figure 1a; cf., [44]). The Dy3+ emission of a completely amorphised zircon sample (N17) is degenerated to a broad hump with a lack of the typical fine-splitting. Indeed, this degenerated broad-band emission of Dy3+ is typical for materials with absent long-range order such as reported for glasses [67,68]. A detailed investigation of spectra of partly radiation damaged samples reveal that fine-split bands not only broaden appreciably with an increasing degree of structural damage, but also superimpose a broad-band contribution that resembles the shape of the Dy emission of the completely amorphised zircon N17 (see red-hatched curve in Figure 1b).
We consequently interpret this contribution in Dy3+ spectra of partly radiation-damaged samples M146 and G3 as emissions arising from Dy ions in a structurally disordered “amorphised” cationic environment. This aspect is illustrated in more detail by Figure 2. A schematic energy level diagram of the 4f9 electronic configuration of Dy3+ is shown in Figure 2a. The most intensive PL transition 4F9/26H13/2 that is of interest here is indicated. Note, that each electronic level is split into fine sub-levels in dependence of the local anionic crystallographic coordination. Each Dy ion, hence, acts as an individual structural probe, in which subtle changes of the local crystal-field may have a strong effect on its luminescence emission. Figure 2b shows a graphical representation of a molecular-dynamics simulation of a single 30 keV recoil in zircon which typically decays from a single α-decay event (modified after Trachenko et al. [69]). For illustration, we graphically add random substitutions of Dy ions on Zr-sites that are used as luminescence probe in this study (blue balls in Figure 2b). Note, that the zircon structure is heavily altered on the local atomic scale upon a recoiling daughter nucleus and cause the Dy-probes to be located finally in various, structurally different sites. These may be fully ordered crystallographic environments which have not been affected by the impact of α-decay events. Other sites in close proximity to the impact area may be (strongly) stressed but still intact, whereas the structure is completely altered along recoil trajectories, i.e., with breakdown of atomic bonds and loss of regular coordination around Dy ions.
The detected PL of Dy3+ from the analysed volume of partly radiation-damaged zircon may, hence, be interpreted to represent a superposition of emissions from a multitude of individual Dy ions in sites with different degree of degradation. A high accumulation of radiation damage consequently result in PL emissions from high number of Dy ions in structurally altered sites. Given the assumption that Dy ions are distributed statistically within the analysed volume of our in-situ PL measurements and all of them contribute equally to the detected luminescence signal (i.e., the quantum conversion is the same and independent among different Dy sites), the relative integrated intensity of the deconvoluted amorphous component in relation to the overall integrated intensity, may give a reasonable estimate of the fraction of the amorphised material (amorphous fraction fa). In the examples of partly radiation-damaged samples M146 and G3 shown in Figure 1b, the red-hatched indicated component corresponds to ca. 45% and 73% of the signals surface integral, respectively.
The relation of the “amorphous” spectral component with respect to the overall signal integral is in principle independent of the absolute integrated intensity. For example, the presence of an increased number of statistically distributed Dy ions, due to the substitution of higher REE contents, may cause a higher emission intensity, but would not substantially alter the relation of Dy ions to be found in the radiation-damaged volume with respect to the total number of Dy ions present (i.e., if the assumption of a statistical distribution of Dy ions holds). Similarly, the estimation of the amorphous fraction by the latter approach may be accurate if other secondary PL quenching effects (e.g., effective absorption by broad-band irradiation-induced defect centres) may cause a substantial decrease of the overall luminescence intensity as typically reported for zircon e.g., in case of photo- and cathodoluminescence [44] or ionoluminescence [71].
A major challenge for the quantitative estimation of the amorphous fraction according to the aforementioned approach is the accurate spectra deconvolution of the broad-band contribution of the Dy3+ emission that is considered to represent the “amorphous” component. Figure 3 presents two possible data-reduction strategies that may be applied (after background correction as indicated in Figure 1a). A first option is the fitting of a multitude of individual peaks that are modelled by Gauss-Lorentz (i.e., pseudo-Voigt) functions with their three spectral parameters (peak position, band width and amplitude) freely adjustable using an iterative least-squares optimization (Figure 3a). The amorphous component is treated here as a single Gaussian function with a fixed FWHM of 540 cm−1 and a peak position at 17,360 cm−1 (~576 nm) according to fitting results of the completely amorphised zircon sample N17 (cf., Figure 1a).
The amplitude of the broad band is then refined together with the free parameters of the other Gauss-Lorentz functions that fit to the narrow, crystal-field split sublevel bands. The major drawback using latter fitting procedure is the large number of variables to be optimised per single spectra. Here, 11 Gauss-Lorentz functions are used, that comprise 3 free parameters each, in addition to the broad-bands amplitude (34 free variables per spectrum in total, cf., Figure 3a). This is generally unproblematic for fitting single spectra, but may cause extensive computing time if large data-sets of hyperspectral maps with 50,000 to 80,000 and more data points (with a single spectrum per measurement point) need to be treated. Furthermore, this fitting procedure is incapable of accurately resolve low or absent fractions of the amorphous component. This is obvious from fitting results of sample “Rata” (Figure 3a), that is considered to represent a highly crystalline zircon sample with no radiation damage accumulated [66]. This fitting effectively introduced a significant amorphous component as background to optimise the determination coefficient (R2) although this sample is considered free of any measurable structural radiation damage.
To eliminate latter problems, we tested an alternative least-squares fitting procedure that is based on multiple reference spectra (so-called “loadings”) with their individual amplitude (scaling) as free parameter optimised against experimental data (classic least-squares fitting, CLS). Loadings comprise spectra that are considered to represent (i) non-radiation damaged, crystalline samples (e.g., sample Rata or other annealed zircon samples); (ii) a reference that represent the amorphous component (e.g., the completely amorphised sample N17, cf., Figure 1a); and (iii) multiple references that represents slightly affected, stressed, but still crystalline volumes (cf., Figure 2b).
The latter non-amorphous, but affected sample volume causes narrow crystal-field split bands to effectively broaden with an increasing accumulation of radiation damage. This band broadening is accounted for in this fitting procedure by using residue spectra of partly-radiation damaged samples with broadened sub-level bands, but with the amorphous component subtracted according to the first fitting procedure described above (cf., Figure 3a). Note, however, that relative band heights especially in crystalline samples may differ due to polarization- and orientation-dependent effects [52]. To minimize errors in the fitting of samples with unknown orientation, it is therefore advisable to use loadings that represents samples of variable crystal orientations. With respect to the first mentioned fitting procedure, the number of free variables to be optimised in this alternative procedure is considerably reduced (i.e., effectively the number of loadings: e.g., 1 × amorphous component + 3 × residual references + 3 × crystalline references; cf., Figure 3b). This low number of free variables substantially minimize computing time of large data-sets, but may concomitantly results in a somewhat larger error because of less-accurate approximation of spectral components. Dependent on the number and type of loadings used, the amorphous fraction fa estimated by this method readily vary ± 5%. The major advantage, however, relies in the more accurate estimation of low and/or absent amorphous components by using suited “crystalline” reference loadings (cf., fitting results of sample Rata in Figure 3b).

3.2. Sri Lankan Zircon as Reference for the Verification of the PL-Based Estimation of the Amorphous Fraction

Zircon samples from the Ratnapura district (Sri Lanka), are ideal candidates to test and verify the applicability of the proposed PL methodology. Due to the prolonged interest in this sample material in the past, the amorphous fraction fa in Sri Lankan zircon may be easily inferred from the knowledge of the α-dose as calculated from U and Th concentrations and the age of the samples, respectively [26,35,36]. In addition to PL data on fa from a series of well-studied homogeneous Sri Lankan samples (see Table A1), we included two large, but strongly heterogeneously zoned Sri Lankan samples to demonstrate the potential of this technique to operate in-situ with a spatial resolution on the micrometre length-scale. Figure 4 shows PL hyperspectral maps of these two zircon samples SL1 (Figure 4a) and SL2 (Figure 4b). Both samples were cut in halves with one half annealed during a heat treatment at 1300 °C for three days, and thereafter prepared as pairs in close proximity (see inset sketch in Figure 4 for details on preparation). This was done to compare counterparts of the same sample that are chemically analogous to the greatest possible extent, but different in the structural state with the heat-treated half, cleared of the accumulated radiation damage. The samples were scanned using an instrument-coupled, software-controlled, mechanical x-y table with each pixel of the map being an individual laser-induced PL point-measurement. Each spectrum was then deconvoluted using an automated CLS fitting procedure according to the principle presented above and the result given colour-coded. As expected, the two halves of samples SL1 and SL2 differ strongly with respect to their structural state. No amorphous component was detected in spectra obtained from the heat-treated halves. In contrast, PL maps of the non-treated halves reveal internal growth zonation as well as secondary overgrowth features (see arrows as indicated in Figure 4b). Strong accumulation of radiation damage is found in zones with elevated concentrations of U and Th (see e.g., concentrations of U and Th indicated in Figure 4a). Element concentrations of U and Th, as well as isotope ratios for U-Pb age determination of various measurement spots on the zircon crystals SL1 and SL2 were obtained using LA-ICP-MS (see Table A2). Note, however, that we placed measurements in different zones of the heterogeneous samples, but just in those providing sufficient space for a laser-ablation spot (spot diameter 50 µm) and multiple PL measurements in close proximity.
Figure 5a shows a plot with the calculated α-dose of each measurement spot together with data from homogeneous Sri Lankan references against their amorphous fraction fa as obtained by PL (see detailed data and equation for calculation of the α-dose in Table A1). Data points resemble very well the correlation of the amorphous fraction fa vs. α-dose (red curve in Figure 5a) given by Zhang and Salje [35]. Latter authors have summarised existent data on the amorphous fraction fa of Sri Lankan zircon as obtained by three alternative bulk methods. Based on their data, a regression according to the direct impact model have been proposed as best-fit that describes the accumulation of the amorphous fraction fa via:
f a = 1 e B a × D α
where Dα is the α-dose and Ba the amount of amorphous material produced per α-recoil. A least-squares regression analysis of our data from Sri Lankan zircon using this model results in Ba = 2.87 × 10−19 g (R2 = 0.99), which is fairly close to a Ba = 2.73 × 10−19 g given by Zhang and Salje [35]. Note, however, that Sri Lankan zircon samples may have experienced a considerable thermal event in the geological history, and, hence the radiation damage now present may be systematically lower than expected if no annealing would have affected their structural state. Based on Raman data, Nasdala et al. [32] have discussed a correction factor of ×0.55 to the calculated α-dose, to account for geothermal annealing of Sri Lankan zircon (“effective” α-dose, see [32]). Within the context of the direct impact model (Equation (1)) this correction concomitantly results in a “true” Ba of ~4.96 × 10−19 g of amorphised material per α-recoil (i.e., ~2.73 × 10−19 g divided by 0.55; see discussion below along with Figure 9b).
Due to the obvious strong variance of fa among different growth zones in sample SL1 and SL2, we studied the impact of structural radiation damage on LA-ICP-MS isotopic age determination in more detail. We compared concordant 206Pb/238U ages obtained from laser-ablation measurement spots in chemically equivalent growth zones on both, the untreated radiation-damaged and the annealed crystal half with the crystallinity recovered (see indicated pairs of measurement spots in Figure 4). 206Pb/238U ages obtained from the untreated, radiation-damaged half of sample SL2 were found to be systematically older with increasing Th + U concentration and, hence, increasing fa present. Equivalent measurement spots on the annealed half of sample SL2 in contrast show no systematic correlation with the Th + U content, but a substantially lower mean 206Pb/238U age. Whereas the radiation-damaged half of SL2 gives a mean 206Pb/238U age of 531 ± 12 Ma (n = 8, MSWD = 17), the annealed half gives 514 ± 3 Ma (n = 8, MSWD = 1.1) with a 4-times improved analytical precision (see Table A4). Bias of the isotopic ages as caused by secondary loss/gain of U or Pb isotopes is rather unlikely in this context and can be excluded here due to the fact that very similar regions of the very same crystals were probed and compared. Here, the detected bias of LA-ICP-MS ages may be interpreted to result from the impact of radiation damage itself. Indeed, various authors reported a radiation-damage matrix-effect that is discussed to be caused by the altered physico-mechanical properties and, hence, laser-ablation behaviour of zircon in dependence of its damage accumulation [72,73,74]. Thus, an appropriate heat treatment of zircon samples prior to U-Pb isotope dating by LA-ICP-MS has been suggested to substantially improve its precision and accuracy (see discussion in [75,76]). Significant changes in physical and chemical properties of radiation-damaged zircon have been reported to typically result in more effective ablation rates that is evident in deeper ablation pits, and generally cause isotope ages to appear systematically older [73,74].
In contrast to the latter, we observed the exact opposite behavior for Sri Lankan sample SL1. The 206Pb/238U ages obtained from various measurement spots on unannealed zircon SL1 systematically appear younger, although an even higher amorphous fraction up to complete amorphisation was detected (cf., spots 1−3 on sample SL1 with fa up to 94%; indicated in Figure 5b, data see Table A2). Note that the trend from older to younger ages with increasing U + Th concentrations does not follow the expected age progression as evident from the growth zonation. Here, the core of zircon SL1 with the highest fa and U + Th concentrations yield the youngest LA-ICP-MS isotopic age (cf., spot SL1u_2 in Figure 5b). The annealing procedure applied to sample SL1 had a minor effect on the mean 206Pb/238U ages only and its precision was not improved substantially. The radiation-damaged half of SL1 gives a mean 206Pb/238U age of 528 ± 9 Ma (n = 11, MSWD = 22; cf., Table A3), whereas the annealed half gives 520 ± 7 Ma (n = 11, MSWD = 16; cf., Table A3).
Moreover, a third zircon example (sample MalawiX) from a completely different geological setting and with a much younger geological age shows no systematic age bias at all; neither in correlation with the U + Th concentration, nor with the amount of radiation damage present (see samples characteristics in Figure 7a below). Note, that this sample contains huge amounts of REEs, U and Th with growth zones that reaches an fa value up to 72% (cf., spot MalawiXu_13 in Table A2). These strongly radiation-damaged zones in Malawi zircon yield very consistent ages in comparison to zones that accumulated almost no damage (cf., Table A5). Accordingly, the difference of mean 206Pb/238U ages between unannealed (112 ± 1.3 Ma; n = 9, MSWD = 3.4) and annealed half (110 ± 1.6 Ma, n = 9, MSWD = 5.4) is small. The above contradicting results, however, indicate that the degree of structural radiation damage alone cannot fully explain the strong variance in isotope age bias of LA-ICP-MS measurements in general.

3.3. PL Photoluminescence Compared to Raman Spectroscopy

A detailed analysis of laser-induced PL spectra of Dy3+ in zircon revealed that crystal field fine-split sublevels obviously change in their spectral parameters with an increasing degree of structural radiation-damage accumulation (see bands labelled Dy1 and Dy2 in Figure 3a). Band widths appreciably broaden and peak positions shift to higher wavenumbers. Surprisingly, these changes correlate closely to changes in the Raman-band parameters of the ν3(SiO4) symmetric stretching vibration as obtained from the Sri Lankan samples studied here (Figure 6). This supports similar results by Lenz et al. [44]. Note, however, that in this study, PL sublevel bands of Dy3+ were fitted with Gauss-Lorentz functions after subtraction of the amorphous component according to the principle discussed above (cf., Figure 3b). With respect to latter PL model (cf., Figure 2b), changes in Dy3+ sublevel band widths and positions strictly speaking refer only to non-affected crystalline and/or marginally affected domains that are stressed and distorted, but not being completely structurally disordered (i.e., ‘amorphised’). Thus, the close correlation of band parameters of PL and Raman spectroscopic features may indicate that very similar information on the short-range order are obtained from two perspectives; i.e., structural distortion at the anionic SiO4-tetraheder as probed by Raman ν3 and that of cationic polyhedra substituted by Dy3+ probed by laser-induced PL.
Spectral band widths and their interpretation with respect to the degree of radiation damage, however, have been discussed to be biased by other factors that cause a random structural distortion and finally lead to inhomogeneous broadening. A major disturbing factor reported in this context is the substitution of especially heterovalent non-formula elements into the crystal structure. This has been documented e.g., for narrow PL bands of Nd3+ [52] and Raman bands [57] of monazite-(Ce). This chemically-induced band-broadening effect is generally considered less important in zircon due to the more restricted crystal-chemical flexibility. Studies by Geisler et al. [53] and Kempe et al. [56], however, revealed that the Raman band-widths of zircon indeed may systematically depend on elevated contents of U and/or REEs. In this study, we included zircon sample material from the Mt. Malosa massif, Malawi that is known for its inhomogeneous broadening of Raman ν3 bands due to the very high amount of REE substitutions [56,59]. The two samples studied represent two different generations of zircon found in pegmatites around the Mt. Malosa, Malawi (see sample description above). Large, prismatic, brownish, intransparent crystals (sample MalawiX) have been interpreted to be of pegmatitic origin that grew in miarolitic cavities [54,56,59], whereas yellowish, semi-transparent crystals (MalawiA) grew as a result of a hydrothermal overprint during and/or shortly after the formation of the pegmatites itself (~120 Ma, [56,59]). PL hyperspectral maps of cut-in-half samples reveal their complex growth and post-growth history (Figure 7a). Note, that the FWHMs of Dy emission bands are not used to estimate the radiation damage here, but analogously to PL maps presented in Figure 4, the amorphous fraction fa was determined using the PL deconvolution procedure described above (cf., Figure 3b).
Again, the annealed crystal halves do not show substantial structural damage left after heat treatment. In contrast, the untreated crystal-halves are characterised by multiple growth zones that accumulated considerable radiation damage in dependence of the amount of substituted U and Th. Crystal growth of the large pegmatite-type zircon crystal MalawiX started from rounded inherited zircon and quartz inclusions (see arrows in Figure 7a). The zircon inclusion is characterised by low U and Th trace-element concentrations of 142 and 80 ppm, respectively, and hence, accumulated almost no significant radiation damage over geological periods of time (spot 11, Figure 7a, Table A2). Simultaneous crystallisation of zircon and quartz indicates Si-saturated growth conditions. Crystallization from a highly fractionated, late-stage growth environment is further indicated by geochemical provenance discrimination based on Hf and Y concentrations after Shnukov [77] as presented in Figure 8b. Concentrations of Hf and Y of the Malawian zircon samples fall well into the discrimination field characteristic of “felsic rocks with high SiO2 content” (group IV, see details in [65]). Inclusions of sample MalawiX were progressively enclosed by a zircon ‘core’ region with dominant pyramidal faces (1 0 1), followed by multiple growth zones that lead to a morphology with dominant prism faces (1 0 0). Note, that zones in pyramidal growth direction are more enriched in U and Th and hence, accumulated more structural radiation damage than zones in growth direction of the prism faces with less damage present (see “sector zoning” indicated by arrows in Figure 7a). Further, outer growth zones of sample MalawiX were heavily affected by secondary alteration (see arrows in Figure 7a again). These altered regions have a fine-grained micro-texture that according to Soman et al. [59] has been discussed to result from aqueous-solution driven dissolution-reprecipitation processes during hydrothermal alteration (cf., their Figure 1). In accordance to observations by latter authors, we detected elevated concentrations of Ca and Fe in this area which is considered a typical indicator for hydrothermally-driven alteration (cf., measurement points 4, 5, 11, and 12 of TE analyses presented in Table A7). In contrast, sample MalawiA show very subtle accumulation of radiation damage only, but especially in overgrowth zones that adds to slightly curvilinear prism faces of the crystal (note adapted scale colour-coding, Figure 7a). No obvious features of secondary alteration comparable to those found in sample MalawiX are detected in sample MalawiA.
Zircon sample MalawiX gives a 206Pb/238U age of 111 ± 1 Ma (n = 18, MSWD = 5) as obtained from weighted mean of measurement spots with concordant isotopic ratios (cf., Table A5). Measurement spots located within heavily altered zones, however, exclusively yield discordant U-Pb isotope ratios with a discordia line having a lower concordia intercept at a similar age of 111 ± 1 Ma (Figure 7b). This discordia age is indistinguishable to the concordant mean 206Pb/238U age and confirms previous findings of a co-genetic and/or temporal closely related alteration during the formation of the pegmatite-type zircon [59]. Note that measurement spots within altered areas are likely to represent mixture analyses of various alteration phases, but no obvious spikes or disturbance of isotope ratios during time-resolved laser-ablation measurements were recorded. The latter, in principal, does not exclude the presence of very fine, homogeneously dispersed submicron inclusions within the ablated volume. The disturbance of U-Pb isotope ratios that leads to the systematic discordia here, however, is interpreted to most likely result from entry of common Pb during the hydrothermal overprint.
Preferential alteration in sample MalawiX took place along growth zones that are characterised by strong damage accumulation as seen in relict zones with high fa in the upper half of the crystal (Figure 7a). The upper half of the crystal was more likely exposed to circulating fluids as this part readily faced towards the open miarolitic cavity in comparison to the lower half close to the contact-wall (see left-side margin of the crystal). The lower half is obviously less altered but shares similar growth-zonation characteristics. Note, however, that no substantial radiation damage could have been accumulated at the time of alteration shortly after formation, and therefore, may not be considered as primary cause for the favored alteration here. We consider a preferential alteration due to the high amount of U, Th and REE substitution in this growth zones (that i.e., also lead to enhanced radiation damage over progressive geological time). The latter may be explained by an enhanced dissolution driving force as induced by the higher solubility of the trace element-rich zircon material in contrast to TE-poor zircon in hydrothermal solution (see discussion in [59]).
Interestingly, zircon sample MalawiA does not show characteristic alteration features, but LA-ICP-MS measurement spots placed across the crystal yielded relatively precise, but discordant U-Pb isotopic ratios (see Figure 7b). Only a few spots gave ratios that lie on the concordia within an error of 2σ (cf., Table A6). The constructed discordia gives a lower intercept age of 112.3 ± 1.6 Ma, that is, again, close to the age of formation and alteration of the pegmatitic zircon MalawiX. The consistent discordia trends of both zircon-types, altered MalawiX as well as MalawiA, underpins their close hydrothermal genetic relation with a very similar isotope disturbance as caused by incorporation of common Pb. Common Pb in yellowish “hydrothermal” type of zircon from Mt. Malosa have been confirmed also by U-Pb sensitive high-resolution ion microprobe analyses reported by Kempe et al. [56]. The hydrothermal nature of zircon sample MalawiA is, in addition, supported by geochemical indications that are illustrated in Figure 8. According to Hoskin [78], zircon which crystallize at hydrothermal conditions typically has elevated concentrations of especially LREEs and therefore gives substantially flattened chondrite-normalised REE concentration patterns in comparison to their magmatic analogue. This is clearly evident from REE patterns shown here (Figure 8a). Moreover, the less pronounced Ce/Ce* anomaly of sample MalawiA in comparison to that of sample MalawiX supports its hydrothermal nature [78]. Note, however, that overall REE concentrations and their patterns are approximately one order of magnitude higher than those reported by Hoskin [78] for their hydrothermal zircon samples. Here, total REE concentrations reach very high values of around 3 wt. %, whereby the REE totals surprisingly do not differ substantially between both Malawian sample types (cf., Table A7 and Table A8).
Their major difference is the strong enrichment of LREEs on the expense of HREE and the obvious lack of P for charge compensation in sample MalawiA (cf., <100 ppm P in Table A8). In accordance to results reported by Kempe et al. [56], this unusual crystal chemistry has a substantial impact on the Raman spectroscopic band widths. In Figure 9a, the correlation of Raman (ν3) FWHMs in dependence on the accumulated α-dose of reference samples are shown. Those references are conclusively discussed not to be affected by geothermal annealing in the geological past and, hence, have retained all of their structural damage [22]. Whereas Raman band-widths obtained from MalawiX scatter around latter calibration line, obtained FWHMs from measurement spots on MalawiA are systematically broader than expected from their calculated, relatively low α-dose (see blue circles in Figure 8a). Obtained FWHMs of 12 to 15 cm−1 typically indicate a rather substantial accumulation of radiation damage (~ 0.9–1.1 × 1018 events/g, [22]), but very low maximum doses of 0.16 × 1018 events/g were calculated from U and Th concentration data for sample MalawiA (Table A2). Raman band-widths of zircon that is (partially) annealed in the geothermal history typically plot at somewhat lowered FWHMs with respect to the calibration given by Nasdala [22]. Data points that plot above this line, hence, require cautious interpretation considering contradicting effects of band broadening. Kempe et al. [56] explained this band broadening by internal strain and distortion of the structure caused by the incorporation of large amounts of trivalent REEs that are widely charge-compensated by oxygen vacancies due to the insufficient quantity of other charge-balancing species such as pentavalent P on the Si-site or U, Nb and Ta, and others on the Zr-site (see discussion in [56]). An alternative charge-balancing mechanism that likewise may effectively compensate excess REEs in zircon involves the incorporation of OH- replacing O2- in the anionic tetrahedral SiO4-unit [79]. Charge balance is maintained here via H+ + REE3+ = Zr4+. Kempe et al. [56] explicitly mentioned that no water or hydroxyl groups were found by Raman spectroscopic analyses. The vast amount of REE substitutions and their luminescence emissions, however, strongly interfere in spectral regions typical of OH and H2O vibrational modes (2000–3500 cm−1). Using multiple excitation wavelengths, we found strong PL artefacts in Raman spectra of the Malawian zircon samples within the respective spectral region (473 nm: Dy3+; 488 nm: Ho3+ and Dy3+; 532 and 514 nm: Sm3+; and 633 nm: Nd3+). If hydroxyl or water species are present in these zircon samples, their Raman signal simply is strongly obscured by latter REE3+ emissions. The potential of Raman spectroscopy to detect OH/H2O species in zircon therefore needs to be questioned [51,80]. The exact nature of the substitution mechanisms that charge-balance the extraordinary amount of REEs in the Malawian zircon samples, hence, remain speculative.
As illustrated in Figure 9b, results of the PL-based determination of fa on the sample MalawiX are fairly in accordance to those of Raman band-width analysis (cf., Figure 9a). The fa vs. α-dose trend obtained for this sample resembles the typical accumulation behavior of zircon that is considered not to be affected by thermal annealing in the geological history. This calibration trend is given here as green curve in Figure 8b and was inferred from the direct impact model of Zhang and Salje [35] for Sri Lankan zircon (cf., Figure 5a), but with α-doses corrected for its geothermal annealing using a factor of ×0.55 as proposed by Nasdala [22]. This correction results in Ba = 4.96 × 10−19 g (see Equation (1) and discussion above). Note, that the latter correction of the direct-impact model by Zhang and Salje [35] yielded a marginally lower Ba than reported by Palenik et al. (Ba = 5.48 × 10−19 g; [8]). They have deduced Ba indirectly from an approximation of Raman band-widths of Sri Lankan zircon to the Raman calibration line by Nasdala et al. [22] using an “effective” damage accumulation time period of 375 Ma (see discussion in [8]).
A regression analysis of Equation (1) to PL data obtained from Malawi samples, however, yields a Ba = 4.86 × 10−19 g that is fairly close to the above mentioned back-calculated “true” Ba of 4.96 × 10−19 g. Values for fa obtained from PL measurement points on sample MalawiA also fall onto latter trend. Based on these results, we may conclude that the pegmatitic-type zircon (sample MalawiX) and also, the hydrothermal-type zircon generation was likely not affected by further subsequent geothermal events after their formation ~111 Ma ago. The expected low accumulation of radiation damage is, in contrast to results from Raman FWHMs, not overestimated here (see again Figure 9). This is ascribed to a major advantage of the CLS-fitting procedure of Dy3+ emission spectra presented above (cf., Figure 3b). A deconvolution that is based on input loadings allows to consider suited reference spectra from the annealed sample half of MalawiA, which restored complete crystallinity upon the applied heat treatment. Subtle changes in the spectral profile as caused by the amorphous component, hence, are detected most efficiently in direct comparison to this reference.

4. Summary and Conclusions

We documented a promising laser-based spectroscopic tool to non-destructively characterise radiation-damage in zircon in-situ on a micrometre length-scale. The proposed structural probe is based on the steady-state laser-induced photoluminescence of Dy3+ and opens up the opportunity to quantitatively visualize the amorphous fraction fa present in complexly zoned, heterogeneous samples using PL hyperspectral mapping. We found conclusive evidence that the fa obtained from Sri Lankan zircon using a classic least-square (CLS) component fitting of Dy3+ emissions is consistent with reported literature values from alternative bulk techniques [33,35,36]. The latter implies further evidence for a direct impact model of radiation-damage accumulation (i.e., the accumulation of fa) in zircon with a progressive α-dose. The proposed methodology, that comprise a deconvolution of the “amorphous” spectral component from Dy3+ luminescence emissions, has the advantage to provide a direct estimation of the amorphous fraction present without the need for calibration with related spectral changes, e.g., spectral band-width broadening of narrow crystal-field split REE3+ luminescence [44] or Raman bands [8], and hence may be used as a complementary in-situ spectroscopic technique. This is of specific surplus for occasionally observed samples that show chemically-induced band broadening of Raman and PL band-widths which biases individual calibrations of FWHM vs. damage accumulation, as exemplified in this study for hydrothermally-grown zircon from Mt. Malosa, Malawi.
Effective, laser-based excitation of Dy3+ luminescence in zircon is preferentially achieved using laser wavelengths in the UV to blue spectral region (e.g., with 488 or 473 nm as used in this study, [51]). Accurate spectral deconvolution of Dy3+ emission spectra requires a minimum spectral resolution of the applied spectrometer system to differentiate the broad-band component ascribed to result from Dy ions in structurally disordered regions from narrow crystal-field split bands that are interpreted to stem from slightly or non-distorted crystalline remnants. We recommend a minimum spectral resolution of better than 3 to 4 cm−1 (~0.1 nm in the yellow spectral range) because narrow Dy3+ bands of the 4F9/26H13/2 transition may have individual FWHMs down to ~10 cm−1 (0.3 nm) at room temperature (cf., Figure 6; [44]). Further, suited reference spectra from well-characterised sample material are needed to serve as loadings for a sound CLS fitting. Most importantly, those comprise spectra that (i) represents samples which underwent complete amorphisation, and (ii) samples that are considered to be fully crystalline samples and free of structural damage. The latter may be synthetic analogues or natural zircon samples that experienced very low α-doses or have been annealed artificially using a suitable heat treatment (e.g., 72 h at 1300 °C).
Thermal annealing, either induced artificially by heat treatments and/or naturally as caused by geothermal events, may have an important impact on the structural state of zircon, and causes the detected amorphous fraction fa to be substantially lower than expected from the calculated α-dose based on U + Th concentrations and the age of the sample. Misfits in obtained fa from a reference correlation, therefore, may provide important indications for potential geothermal events in the sample’s geological history. In this study, we found that investigated Malawian zircon samples did not experience a substantial impact from geothermal events after their formation ~111 Ma ago.
Careful comparison of LA-ICP-MS isotope ages of naturally, radiation damaged zircon with its chemically analogous, but annealed damage-free counterpart further reveal that the ablation behavior of altered zircon matrix may lead to a systematic bias of the isotope fractionation in dependence on the amount of accumulated radiation-damage. The mere quantity of damage alone, however, was found not to be an exclusive factor for obtained biases. Based on these results, we strongly suggest to initiate further studies that shed light on the role of structural radiation-damage affective to the analytical precision of LA-ICP-MS isotope dating; especially with focus on the impact of geothermal damage-annealing on the sub-microstructure of zircon.
We propose here that the concept of the presented PL methodology may be extended easily (a) to other accessory minerals and ceramic waste-form materials of interest, and (b) to other REE3+ species to be used as potential structural probes, especially those that are excitable/emissive in complementary spectral ranges (cf. for instance, Nd3+ in orthophosphates [58]). The estimation of fa in zircon and other mineral phases using laser-induced PL spectroscopy has numerous analytical advantages. One of the major benefits is the possibility to perform spot analyses non-destructively using highly confocal spectrometers coupled to optical probe-heads or microscopes with the latter operating on the µm length-scale. Further, this technique has great potential to be applied as a remote technique with minimal sample preparation needed, e.g., for minute in-situ monitoring of self-irradiation experiments with 238Pu-, 241Am, or 244Cm-doped analogue ceramics through transparent containers.

Author Contributions

All authors listed have made a substantial, direct, and intellectual contribution to the presented study. Conceptualization & project administration, C.L., G.R.L.; methodology, C.L.; E.B.; draft preparation, C.L.; review and editing, G.R.L., E.B. All authors have read and agreed to the published version of the manuscript.

Funding

C.L. and E.B. gratefully acknowledge the use of instrumentation funded through an honorary-associate agreement with the Arc Ccfs/Gemoc at Macquarie University, Sydney. The analytical data were obtained using instrumentation funded by DEST Systemic Infrastructure Grants, ARC LIEF, NCRIS/AuScope, industry partners and Macquarie University. This is contribution 1418 from the ARC Centre of Excellence for Core to Crust Fluid Systems (http://www.ccfs.mq.edu.au) and 1361 in the GEMOC Key Centre (http://www.gemoc.mq.edu.au). Financial support of C.L. and funding for publishing open access was provided by the Austrian Science Fund (FWF), through project J3662-N19 granted to C.L.

Acknowledgments

Thanks are due to Lu, K., and Palmer, T. (ANSTO) and Wagner, A. (University Vienna) for sample preparation. We thank Nasdala, L. (University Vienna) providing Sri Lankan zircon reference material, with sample M146 provided by Hofmeister, W. (Institute of Gemstone Research, Idar Oberstein and Mainz), and zircon G3 by Kennedy, A.K. (Curtin University, Perth). Sample MalawiA was made available by Kolitsch, U. (NHM Vienna; Sample reference-No.: N3848). We are indebted to Karatchevtseva, I. (ANSTO), Chanmuang Nasdala, C. and Talla, D. (University Vienna) for technical assistance and furnace heat-treatments of samples. We also gratefully acknowledge technical assistance with LA-ICP-MS instrumentation by Murphy, R. and Lai, Y.-J. (Macquarie University, Sydney).

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A

Table A1. Uranium and thorium concentrations, ages, calculated α-doses, and the amorphous fraction fa of various homogeneous, Sri Lankan zircon crystals and a damage-free reference from the Ratanakiri district, Cambodia (Rata).
Table A1. Uranium and thorium concentrations, ages, calculated α-doses, and the amorphous fraction fa of various homogeneous, Sri Lankan zircon crystals and a damage-free reference from the Ratanakiri district, Cambodia (Rata).
SampleU *
[ppm]
Th *
[ppm]
206Pb/238U Age *
[Ma]
Calc. α-dose †
[×1018 events/g]
Amorphous Fraction ‡ fa [%]
CZ3550 ± 1030 ± 2563.9 ± 1.31.05–1.1021.8
B188556 ± 2459 ± 4559 ± 81.03–1.1525.8
BR1796 ± 1339 ± 1558 ± 131.47–1.6036.2
BR231772 ± 10109 ± 2571 ± 41.53–1.5939.7
M257840 ± 27235 ± 20561.3 ± 0.31.66–1.7837.1
BR266909 ± 17201 ± 7559.0 ± 0.31.77–1.8837.5
M146923 ± 17411 ± 9567 ± 41.92–2.0344.6
M127923 ± 23439 ± 11524.3 ± 0.41.78–1.8843.4
OR11490 ± 70279 ± 18522 ± 32.62–2.9254.6
G1681499 ± 33257 ± 9547 ± 32.83–3.0058.4
G42355 ± 84330 ± 12564 ± 54.48–4.9176.0
G32572 ± 96585 ± 34542 ± 54.77–5.2576.3
N175568 ± 114344 ± 45551 ± 710.27–11.0297.4
Rata140 ± 10090 ± 801.2 ± 0.3±00.0
Notes: Errors quoted are estimates that include errors of individual measurements and sample heterogeneity (scatter among multiple analyses).
* Sources for U and Th concentrations, ages, calculated α-doses: Nasdala et al. [39] for M257, Nasdala et al. [40] for M127, Zeug et al. [66] for zircon from Ratanakiri, Cambodia (Rata). See Nasdala et al. [32] (and references therein) for all other samples.
† α-doses were calculated from U–Pb age and U and Th concentrations after: D = 8 · c u · N A · 0.9928 M 238 · 10 6 · ( e λ 238 · t 1 ) + 7 · c u · N A · 0.0072 M 235 · 10 6 · ( e λ 235 · t 1 ) + 6 · c T h · N A M 232 · 10 6 · ( e λ 232 · t 1 ) , where cu and cTh are the present actinide concentrations; NA is Avogadro’s constant; M238, M235 and M232 are the molecular weights of the parent isotopes; λ238, λ235 and λ232 are the respective decay constants and t is the age [6,22,25].
‡ The amorphous fraction as obtained with the PL method (cf., Figure 1b). Minimum errors of ± 5% need to be considered for scatter among multiple analyses and deviations due to different model loadings (cf., Figure 3b).
Table A2. Uranium and thorium concentrations, calculated α-doses, and the amorphous fraction fa obtained from multiple measurement spots on the unannealed (u), heterogeneous crystal halves of sample SL1 and SL2 from the Ratnapura district, Sri Lanka, and the samples MalawiA and MalawiX from the Mt. Malosa massif, Malawi.
Table A2. Uranium and thorium concentrations, calculated α-doses, and the amorphous fraction fa obtained from multiple measurement spots on the unannealed (u), heterogeneous crystal halves of sample SL1 and SL2 from the Ratnapura district, Sri Lanka, and the samples MalawiA and MalawiX from the Mt. Malosa massif, Malawi.
SamplesU [ppm]Th [ppm]Calc. α-dose †
[×1018 events/g]
fa ‡ [%]
Sample SL1
SL1u_01314625206.6684.5
SL1u_02418535198.9394.4
SL1u_03426936229.1390.2
SL1u_066791954.0633.3
SL1u_0711611073.4539.1
SL1u_0818071181.3058.6
SL1u_09180611822.1262.3
SL1u_10198812113.2869.1
SL1u_1115387533.7256.1
Sample SL2
SL2u_112717104.5369.8
SL2u_210103413.5460.7
SL2u_39295121.6738.2
SL2u_49404301.7742.8
SL2u_55382160.9727.6
SL2u_66603391.1431.2
SL2u_74351681.2134.5
SL2u_84441340.8524.4
MalawiX
MalawiXu_0135679340.8137.5
MalawiXu_0215110,0410.9238.8
MalawiXu_0336419,4651.8054.3
MalawiXu_0492160180.8535.4
MalawiXu_05454357042.1466.5
MalawiXu_0658613,9591.4145.8
MalawiXu_079292890.8332.7
MalawiXu_0814610,7880.9844.0
MalawiXu_0924410,1380.9638.9
MalawiXu_10474458412.2365.8
MalawiXu_11142800.060.0
MalawiXu_12515341742.2360.7
MalawiXu_13568271052.6871.8
MalawiA
MalawiAu_11123660.072.2
MalawiAu_21303200.087.4
MalawiAu_31133280.070.0
MalawiAu_41193690.080.1
MalawiAu_51063080.070.0
MalawiAu_615710210.155.7
MalawiAu_714612720.167.6
Notes: † α-doses were calculated from U and Th concentrations and the U-Pb age using equation given in Table A1 (see also [6,22,25]). Mean weighted 206Pb/238U ages as obtained from multiple LA-ICP-MS measurement spots on both, annealed and un-annealed halves were used for calculation (cf., Table A3, Table A4 and Table A5).
‡ The amorphous fraction as obtained with the PL method (cf., Figure 3b). Errors of ± 5% need to be considered for scatter among multiple analyses and background correction.
Table A3. U-Pb isotope ratios and 206Pb/238U ages of the large, heterogeneous Sri Lankan zircon crystal SL1 with multiple analyses spots on both, the radiation damage-affected unannealed half (u) and the heat-treated, annealed half (a) with its crystallinity recovered.
Table A3. U-Pb isotope ratios and 206Pb/238U ages of the large, heterogeneous Sri Lankan zircon crystal SL1 with multiple analyses spots on both, the radiation damage-affected unannealed half (u) and the heat-treated, annealed half (a) with its crystallinity recovered.
Analysis SpotIsotope RatiosAges [Ma]
207Pb/235U±2σ206Pb/238U±2σ206Pb/238U±2σ
SL1u_010.665380.00950.083580.0009517.55
SL1u_020.664480.00950.082650.0009511.95
SL1u_030.665480.00960.082150.00095095
SL1u_040.687750.01050.086350.0009533.96
SL1u_050.677410.01960.083710.0012518.37
SL1u_060.702310.01280.087940.0010543.46
SL1u_070.723130.01270.087960.0010543.56
SL1u_080.704120.01240.087560.0010541.16
SL1u_090.701720.01260.086160.0010532.86
SL1u_100.696880.01280.084930.0010525.56
SL1u_110.698290.01330.087360.0010539.96
Unannealed half (u) weighted mean age and 95% confidence
528 ± 9 Ma (n = 11, MSWD = 22)
SL1a_010.651960.00940.081440.0009504.75
SL1a_020.662180.00970.082880.0009513.35
SL1a_030.656260.00980.082460.0009510.85
SL1a_040.658400.01000.082420.0009510.55
SL1a_050.671280.01060.083050.0009514.35
SL1a_060.690620.01250.086080.0010532.36
SL1a_070.687260.01230.086670.0010535.86
SL1a_080.671320.01170.084310.0010521.86
SL1a_090.683690.01240.085820.0010530.86
SL1a_100.679360.01260.084660.0010523.96
SL1a_110.689500.01330.086370.00105346
Annealed half (a) weighted mean age and 95% confidence
520 ± 7 Ma (n = 11, MSWD = 16)
Both SL1 halves weighted mean age and 95% confidence
524 ± 6 Ma (n = 22, MSWD = 20)
Table A4. U-Pb isotope ratios and 206Pb/238U ages of the large, heterogeneous Sri Lankan zircon crystal SL2 with multiple analyses spots on both, the radiation damage-affected unannealed half (u) and the heat-treated, annealed half (a) with its crystallinity recovered.
Table A4. U-Pb isotope ratios and 206Pb/238U ages of the large, heterogeneous Sri Lankan zircon crystal SL2 with multiple analyses spots on both, the radiation damage-affected unannealed half (u) and the heat-treated, annealed half (a) with its crystallinity recovered.
Analysis SpotIsotope RatiosAges [Ma]
207Pb/235U±2σ206Pb/238U±2σ206Pb/238U±2σ
SL2u_11.145800.02760.091450.0012564.17
SL2u_20.749780.01870.086780.0012536.57
SL2u_30.712830.01900.085560.0012529.27
SL2u_40.665280.01790.085740.0012530.37
SL2u_50.638460.02240.084310.0014521.88
SL2u_60.644430.01940.084520.0012523.17
SL2u_70.648260.01890.084410.0012522.47
SL2u_80.640510.01990.083710.0012518.27
Unannealed half (u) weighted mean age and 95% confidence
531 ± 12 Ma (n = 8, MSWD = 17)
SL2a_10.660470.01630.083270.0011515.67
SL2a_20.683700.01780.082620.0011511.77
SL2a_30.662190.01730.083240.0011515.57
SL2a_40.642470.01700.083630.0011517.87
SL2a_50.649470.01900.082840.0012513.17
SL2a_60.634620.01810.083400.0012516.47
SL2a_70.626650.01930.082110.0012508.77
SL2a_80.629650.01960.081970.0012507.97
Annealed half (a) weighted mean age and 95% confidence
514 ± 3 Ma (n = 8, MSWD = 1.1)
Both SL2 halves weighted mean age and 95% confidence
522 ± 7 Ma (n = 16, MSWD = 15)
Table A5. U-Pb isotope ratios and 206Pb/238U ages of the large, heterogeneous zircon crystal MalawiX from the Mt. Malosa Massif, Malawi with multiple analyses spots on both, the radiation damage-affected unannealed half (u) and the heat-treated, annealed half (a).
Table A5. U-Pb isotope ratios and 206Pb/238U ages of the large, heterogeneous zircon crystal MalawiX from the Mt. Malosa Massif, Malawi with multiple analyses spots on both, the radiation damage-affected unannealed half (u) and the heat-treated, annealed half (a).
Analysis SpotIsotope RatiosAges [Ma]
207Pb/235U±2σ206Pb/238U±2σ206Pb/238U±2σ
MalawiXu_010.119310.00640.017240.0003110.22
MalawiXu_021.281050.03320.027910.0005177.53
MalawiXu_030.120400.00570.017860.0003114.12
MalawiXu_040.120480.00390.017460.0003111.62
MalawiXu_050.115120.00260.017230.0002110.11
MalawiXu_060.136470.01250.018470.00061184
MalawiXu_070.132010.01290.018670.0005119.23
MalawiXu_080.643190.02990.023120.0006147.34
MalawiXu_090.394920.01550.020270.0004129.43
MalawiXu_100.110540.00330.017500.0003111.92
MalawiXu_110.115050.01060.017550.0005112.23
MalawiXu_120.110670.00420.017680.00031132
MalawiXu_130.107880.00400.017520.00031122
Unannealed half (u) weighted mean age and 95% confidence
112 ± 1.3 Ma (n = 9, MSWD = 3.4)
MalawiXa_010.117860.00480.017190.0003109.92
MalawiXa_020.564360.01870.021050.0004134.33
MalawiXa_030.123540.00590.017680.00031132
MalawiXa_040.116480.00340.017180.0002109.81
MalawiXa_050.113600.00250.016740.00021071
MalawiXa_060.117610.00620.017620.0003112.62
MalawiXa_070.961490.04290.024960.0008158.95
MalawiXa_080.476040.01810.020460.0004130.63
MalawiXa_090.207880.01050.018690.0004119.43
MalawiXa_100.111530.00340.017360.00031112
MalawiXa_110.106860.00910.017740.0005113.33
MalawiXa_120.113350.00430.017210.00031102
MalawiXa_130.105920.00380.017190.0003109.82
Annealed half (a) weighted mean age and 95% confidence
110 ± 1.6 Ma (n = 9, MSWD = 5.4)
Both halves weighted mean age and 95% confidence
111± 1 Ma (n = 18, MSWD = 5.2)
Notes: Measurement spots with discordant isotopic ratios within ±2σ are written in grey-colour and were not considered for calculation of the weighted mean 206Pb/238U age.
Table A6. U-Pb isotope ratios and 206Pb/238U ages of the large, heterogeneous zircon crystal MalawiA from the Mt. Malosa Massif, Malawi with multiple analyses spots on both, the radiation damage-affected unannealed half (u) and the heat-treated annealed half (a).
Table A6. U-Pb isotope ratios and 206Pb/238U ages of the large, heterogeneous zircon crystal MalawiA from the Mt. Malosa Massif, Malawi with multiple analyses spots on both, the radiation damage-affected unannealed half (u) and the heat-treated annealed half (a).
Analysis SpotIsotope RatiosAges [Ma]
207Pb/235U±2σ206Pb/238U±2σ206Pb/238U±2σ
MalawiAu_10.293150.01710.019150.0005122.33
MalawiAu_20.121010.01080.017150.0005109.63
MalawiAu_30.256110.01640.018570.0005118.63
MalawiAu_40.203700.01430.018850.0005120.43
MalawiAu_51.701620.05430.031810.0008201.85
MalawiAu_60.132620.01020.017970.0005114.83
MalawiAu_70.157080.01130.017840.00051143
MalawiAa_10.224440.01260.019000.0005121.43
MalawiAa_20.196240.01330.017910.0005114.43
MalawiAa_30.177840.01290.017560.0005112.23
MalawiAa_40.257480.01560.018980.0005121.23
MalawiAa_51.971380.05920.033900.0008214.95
MalawiAa_60.287090.01530.019260.00051233
MalawiAa_70.552530.02450.021510.0006137.24
Notes: Isotopic ratios obtained from a number of measurement spots were found to be discordant within ±2σ (written in grey-colour). Calculation of mean weighted 206Pb/238U ages is therefore considered misleading for this sample. The age of this sample was inferred from the lower intercept of the discordia 112.3 ± 1.6 Ma (see Figure 7b).
Table A7. Trace-element concentrations from multiple spot measurements on zircon sample MalawiX from Mt. Malosa, Malawi as obtained by LA-ICP-MS (given in ppm, except for Hf and REE totals that are given in wt. %; totals of REEs comprise Y + La to Lu). Note, that measurement spots for trace-elements differ from those chosen for U-Pb isotopic dating (Table A5 and Figure 7a). Spot locations are indicated in a PL hyperspectral map with the integrated intensity of the Dy3+ signal given color-coded.
Table A7. Trace-element concentrations from multiple spot measurements on zircon sample MalawiX from Mt. Malosa, Malawi as obtained by LA-ICP-MS (given in ppm, except for Hf and REE totals that are given in wt. %; totals of REEs comprise Y + La to Lu). Note, that measurement spots for trace-elements differ from those chosen for U-Pb isotopic dating (Table A5 and Figure 7a). Spot locations are indicated in a PL hyperspectral map with the integrated intensity of the Dy3+ signal given color-coded.
Element Concentrations [ppm] Minerals 10 00083 i001
Spot:1234567891011121314
Hf1.221.291.320.981.001.071.060.450.440.441.021.000.990.97
P165415281533208311701002942377395388992838982933
Ca<51<43<43823259<45<41<4659<46160<4649<46
Ti221632465565111046
Fe27549632156127886711369333282229109382016154172
Nb750857881231111879249584134243994251014734835
Ta2733553984444737363626317565450
Pb2729301821151430729433022191720
Th55228648994816,94615,49615,40813,449543852025811722215,23017,87321,653
U382412434220252205205449243304842161253255303
Y20,04519,66719,561967010,00222,85322,73713,30013,27613,30814,08423,18722,64923,861
La10137667100007500
Ce318285280181152524550226232241237621615698
Pr0016214321123536
Nd548240676755262531308375120
Sm192225107762522309693106104288287379
Eu4551312434117161819494962
Gd17219220236337213561298627621665578143714751724
Tb939699104112372363196194204177386391429
Dy14431434144110551125328732662007198920651746336133843540
Ho567552550299314774769514509522469781769783
Er28792744271611421207244724741587157915841759247423772403
Tm618586582195205361368211209205306365347343
Yb47714488452813151409219522591130112110682165222820992051
Lu53050651214215521722310010094244223201196
REEs [wt %]3.153.063.051.531.533.473.462.041.992.012.193.553.473.66
Table A8. Trace-element concentrations from multiple spot measurements on zircon sample MalawiA from Mt. Malosa, Malawi as obtained by LA-ICP-MS (given in ppm, except for Hf and REE totals that are given in wt.-%; totals of REEs comprise Y + La to Lu). Note, that measurement spots for trace-elements differ from those chosen for U-Pb isotopic dating (Table A6 and Figure 7a). Spot locations are indicated in a PL hyperspectral map with the integrated intensity of the Dy3+ signal given color-coded.
Table A8. Trace-element concentrations from multiple spot measurements on zircon sample MalawiA from Mt. Malosa, Malawi as obtained by LA-ICP-MS (given in ppm, except for Hf and REE totals that are given in wt.-%; totals of REEs comprise Y + La to Lu). Note, that measurement spots for trace-elements differ from those chosen for U-Pb isotopic dating (Table A6 and Figure 7a). Spot locations are indicated in a PL hyperspectral map with the integrated intensity of the Dy3+ signal given color-coded.
Element Concentrations [ppm] Minerals 10 00083 i002
Spot:1234567891011121314
Hf1.261.121.121.181.161.071.061.071.081.121.151.181.171.19
P74656675686068555087686489139
Ca<72<67526558<63<71<68<65<62<62<61<64<62
Ti974564211215116316634707356
Fe911045813020113184119828531142252378
Nb41712419525565291680137113761395682754127221071592
Ta71814771826272513162086
Pb9101010101199998101113
Th300411401190534516691387349361303255313013921521
U125152147146151155133130128136128148159181
Y1067984769171155281577795518779795966321506311650901380749777
La112511392110599100613351771245141217898
Ce25097105607917021454669578098226724820732155987871306424
Pr1579317361249590096595794213315414911438895
Nd13254756435310338803782355633373056106010847507115529947
Sm71511071216689657911717604517659557172535575460
Eu5868776262584437305843101185326
Gd105487499810761074789599502409941696139920503525
Tb2151481682422441471139578201146225264434
Dy15781039115519601972109890779865516521228146414422123
Ho350239255471482271239218182422330303273361
Er1122891935176417821026980919778170813699858461033
Tm193169175332339198199189160335278169146165
Yb14951390142827112755165216981625137828042377130311181211
Lu174167172314320195207200172332287151132141
REEs [wt %]2.162.762.712.802.792.792.782.702.402.752.243.613.844.19
Table A9. Trace-element concentrations from multiple spot measurements on zircon sample SL2 from the Ratnapura district, Sri Lanka, as obtained by LA-ICP-MS. Concentrations given in ppm, except for Hf that are given in wt.-%. Totals of REEs comprise Y + La to Lu. Note, that measurement spots for trace-elements differ from those chosen for U-Pb isotopic dating (Table A4, Figure 4).
Table A9. Trace-element concentrations from multiple spot measurements on zircon sample SL2 from the Ratnapura district, Sri Lanka, as obtained by LA-ICP-MS. Concentrations given in ppm, except for Hf that are given in wt.-%. Totals of REEs comprise Y + La to Lu. Note, that measurement spots for trace-elements differ from those chosen for U-Pb isotopic dating (Table A4, Figure 4).
Element Concentrations [ppm]
Spot:12345678910111213
Hf0.790.780.800.630.620.640.640.780.750.720.980.700.80
P30332542414447477634464327
Ca<55<58<60<57<59<56<55<59<66<57<66<70<60
Ti4.74.34.56.36.15.15.27.824.06.74.85.65.5
Fe<2<2<2<2<2<2<2<2<247477<3
Nb0.60.60.40.70.70.80.80.81.70.50.40.50.6
Ta0.10.10.10.60.50.50.40.40.40.10.00.20.2
Pb16116016117816922722522015621566273106
Th1581711321641522042042431763253745266
U498492500541515691675667483663206843323
Y4594754491281171501501564905148553087
Lab.d.b.d.0.0030.09b.d.0.0050.0050.0020.380.040.010.030.005
Ce2.42.52.21.71.72.22.22.2157.42.33.72.32.1
Prb.d.0.020.030.02b.d.0.020.010.010.240.08b.d.0.09b.d.
Nd0.30.40.40.20.20.40.10.31.91.4b.d.1.50.1
Sm1.11.10.90.40.30.50.60.61.62.70.22.90.3
Eu0.30.40.30.10.10.10.20.21.10.80.10.60.1
Gd6.86.86.71.81.82.72.32.99.110.40.910.91.1
Tb2.32.42.30.70.60.80.80.82.83.10.43.20.4
Dy3233309810101137405416
Ho1314134344515152163
Er73747120192323247580148214
Tm1818185566619194204
Yb18819318655526461661951994420440
Lu333432101011111234359367
REEs [ppm]8298548112362182762712871038922168950164

References

  1. Pabst, A. The metamict state. Am. Mineral. J. Earth Planet. Mater. 1952, 37, 137–157. [Google Scholar]
  2. Ewing, R.C. The metamict state: 1993—The centennial. Nucl. Instrum. Methods Phys. Res. Sect. B Beam Interact. Mater. At. 1994, 91, 22–29. [Google Scholar] [CrossRef]
  3. Nasdala, L.; Grambole, D.; Ruschel, K. Review of effects of radiation damage on the luminescence emission of minerals, and the example of He-irradiated CePO4. Mineral. Petrol. 2013, 107, 441–454. [Google Scholar] [CrossRef]
  4. Weber, W.J. Radiation-induced defects and amorphization in zircon. J. Mater. Res. 1990, 5, 2687–2697. [Google Scholar] [CrossRef]
  5. Trachenko, K.O.; Dove, M.T.; Salje, E.K. Atomistic modelling of radiation damage in zircon. J. Phys. Condens. Matter 2001, 13, 1947. [Google Scholar] [CrossRef]
  6. Holland, H.D.; Gottfried, D. The effect of nuclear radiation on the structure of zircon. Acta Crystallogr. 1955, 8, 291–300. [Google Scholar] [CrossRef]
  7. Chakoumakos, B.C.; Oliver, W.C.; Lumpkin, G.R.; Ewing, R.C. Hardness and elastic modulus of zircon as a function of heavy-particle irradiation dose: I. In situ α-decay event damage. Radiat. Eff. Defects Solids 1991, 118, 393–403. [Google Scholar] [CrossRef]
  8. Palenik, C.S.; Nasdala, L.; Ewing, R.C. Radiation damage in zircon. Am. Mineral. 2003, 88, 770–781. [Google Scholar] [CrossRef]
  9. Beirau, T.; Nix, W.D.; Pöllmann, H.; Ewing, R.C. Radiation-induced effects on the mechanical properties of natural ZrSiO4: Double cascade-overlap damage accumulation. Phys. Chem. Miner. 2018, 45, 435–442. [Google Scholar] [CrossRef]
  10. Balan, E.; Neuville, D.R.; Trocellier, P.; Fritsch, E.; Muller, J.-P.; Calas, G. Metamictization and chemical durability of detrital zircon. Am. Mineral. 2001, 86, 1025–1033. [Google Scholar] [CrossRef]
  11. Mathieu, R.; Zetterström, L.; Cuney, M.; Gauthier-Lafaye, F.; Hidaka, H. Alteration of monazite and zircon and lead migration as geochemical tracers of fluid paleocirculations around the Oklo–Okélobondo and Bangombé natural nuclear reaction zones (Franceville basin, Gabon). Chem. Geol. 2001, 171, 147–171. [Google Scholar] [CrossRef]
  12. Geisler, T.; Trachenko, K.; Ríos, S.; Dove, M.T.; Salje, E.K. Impact of self-irradiation damage on the aqueous durability of zircon (ZrSiO4): Implications for its suitability as a nuclear waste form. J. Phys. Condens. Matter 2003, 15, L597–L605. [Google Scholar] [CrossRef]
  13. Geisler, T.; Pidgeon, R.T.; Kurtz, R.; Van Bronswijk, W.; Schleicher, H. Experimental hydrothermal alteration of partially metamict zircon. Am. Mineral. 2003, 88, 1496–1513. [Google Scholar] [CrossRef]
  14. Kuiper, Y.D. Isotopic age constraints from electron microprobe U-Th-Pb dates, using a three-dimensional concordia diagram. Am. Mineral. 2005, 90, 586–591. [Google Scholar] [CrossRef]
  15. Zamyatin, D.A.; Shchapova, Y.V.; Votyakov, S.L.; Nasdala, L.; Lenz, C. Alteration and chemical U-Th-total Pb dating of heterogeneous high-uranium zircon from a pegmatite from the Aduiskii Massif, Middle Urals, Russia. Mineral. Petrol. 2017, 111, 475–497. [Google Scholar] [CrossRef]
  16. Lumpkin, G.R. Ceramic waste forms for actinides. Elements 2006, 2, 365–372. [Google Scholar] [CrossRef]
  17. Lumpkin, G.R. Ceramic host phases for nuclear waste remediation. In Experimental and Theoretical Approaches to Actinide Chemistry; Gibson, J.K., de Jong, W.A., Eds.; John Wiley & Sons: Hoboken, NJ, USA, 2018. [Google Scholar]
  18. Weber, W.J.; Ewing, R.C.; Catlow, C.R.A.; De La Rubia, T.D.; Hobbs, L.W.; Kinoshita, C.; Matzke, H.; Motta, A.T.; Nastasi, M.; Salje, E.K.H. Radiation effects in crystalline ceramics for the immobilization of high-level nuclear waste and plutonium. J. Mater. Res. 1998, 13, 1434–1484. [Google Scholar] [CrossRef] [Green Version]
  19. Ewing, R.C. The design and evaluation of nuclear-waste forms: Clues from mineralogy. Can. Mineral. 2001, 39, 697–715. [Google Scholar] [CrossRef] [Green Version]
  20. Chakoumakos, B.C.; Murakami, T.; Lumpkin, G.R.; Ewing, R.C. Alpha-decay-induced fracturing in zircon: The transition from the crystalline to the metamict state. Science 1987, 236, 1556–1559. [Google Scholar] [CrossRef]
  21. Sláma, J.; Košler, J.; Condon, D.J.; Crowley, J.L.; Gerdes, A.; Hanchar, J.M.; Horstwood, M.S.; Morris, G.A.; Nasdala, L.; Norberg, N. Plešovice zircon—A new natural reference material for U–Pb and Hf isotopic microanalysis. Chem. Geol. 2008, 249, 1–35. [Google Scholar] [CrossRef]
  22. Nasdala, L.; Wenzel, M.; Vavra, G.; Irmer, G.; Wenzel, T.; Kober, B. Metamictisation of natural zircon: Accumulation versus thermal annealing of radioactivity-induced damage. Contrib. Mineral. Petrol. 2001, 141, 125–144. [Google Scholar] [CrossRef]
  23. Nasdala, L.; Hanchar, J.M.; Rhede, D.; Kennedy, A.K.; Váczi, T. Retention of uranium in complexly altered zircon: An example from Bancroft, Ontario. Chem. Geol. 2010, 269, 290–300. [Google Scholar] [CrossRef] [Green Version]
  24. Yada, K.; Tanji, T.; Sunagawa, I. Application of lattice imagery to radiation damage investigation in natural zircon. Phys. Chem. Miner. 1981, 7, 47–52. [Google Scholar] [CrossRef]
  25. Murakami, T.; Chakoumakos, B.C.; Ewing, R.C.; Lumpkin, G.R.; Weber, W.J. Alpha-decay event damage in zircon. Am. Mineral. 1991, 76, 1510–1532. [Google Scholar]
  26. Weber, W.J.; Ewing, R.C.; Wang, L.-M. The radiation-induced crystalline-to-amorphous transition in zircon. J. Mater. Res. 1994, 9, 688–698. [Google Scholar] [CrossRef]
  27. Lumpkin, G.R.; Smith, K.L.; Gieré, R. Application of analytical electron microscopy to the study of radiation damage in the complex oxide mineral zirconolite. Micron 1997, 28, 57–68. [Google Scholar] [CrossRef]
  28. Meldrum, A.; Boatner, L.A.; Weber, W.J.; Ewing, R.C. Radiation damage in zircon and monazite. Geochim. Cosmochim. Acta 1998, 62, 2509–2520. [Google Scholar] [CrossRef]
  29. Ewing, R.C. Nuclear waste forms for actinides. Proc. Natl. Acad. Sci. USA 1999, 96, 3432–3439. [Google Scholar] [CrossRef] [Green Version]
  30. Capitani, G.C.; Leroux, H.; Doukhan, J.C.; Ríos, S.; Zhang, M.; Salje, E.K.H. A TEM investigation of natural metamict zircons: Structure and recovery of amorphous domains. Phys. Chem. Miner. 2000, 27, 545–556. [Google Scholar] [CrossRef]
  31. Nasdala, L.; Lengauer, C.L.; Hanchar, J.M.; Kronz, A.; Wirth, R.; Blanc, P.; Kennedy, A.K.; Seydoux-Guillaume, A.-M. Annealing radiation damage and the recovery of cathodoluminescence. Chem. Geol. 2002, 191, 121–140. [Google Scholar] [CrossRef] [Green Version]
  32. Nasdala, L.; Reiners, P.W.; Garver, J.I.; Kennedy, A.K.; Stern, R.A.; Balan, E.; Wirth, R. Incomplete retention of radiation damage in zircon from Sri Lanka. Am. Mineral. 2004, 89, 219–231. [Google Scholar] [CrossRef]
  33. Ríos, S.; Salje, E.K.; Zhang, M.; Ewing, R.C. Amorphization in zircon: Evidence for direct impact damage. J. Phys.: Condens. Matter 2000, 12, 2401. [Google Scholar] [CrossRef] [Green Version]
  34. Park, S.; Lang, M.; Tracy, C.L.; Zhang, J.; Zhang, F.; Trautmann, C.; Kluth, P.; Rodriguez, M.D.; Ewing, R.C. Swift heavy ion irradiation-induced amorphization of La2Ti2O7. Nucl. Instrum. Methods Phys. Res. Sect. B Beam Interact. Mater. At. 2014, 326, 145–149. [Google Scholar] [CrossRef]
  35. Zhang, M.; Salje, E.K. Infrared spectroscopic analysis of zircon: Radiation damage and the metamict state. J. Phys. Condens. Matter 2001, 13, 3057. [Google Scholar] [CrossRef]
  36. Farnan, I.; Salje, E.K.H. The degree and nature of radiation damage in zircon observed by 29Si nuclear magnetic resonance. J. Appl. Phys. 2001, 89, 2084–2090. [Google Scholar] [CrossRef]
  37. Kinny, P.D.; Compston, W.; Williams, I.S. A reconnaissance ion-probe study of hafnium isotopes in zircons. Geochim. Cosmochim. Acta 1991, 55, 849–859. [Google Scholar] [CrossRef]
  38. Pidgeon, R.T. Calibration of zircon standards for the Curtin SHRIMP II. In Proceedings of the Abstract of Eighth International Conference on Geochronology; Cosmochronology and Isotope Geology, U.S. Geological Survey: Reston, VA, USA, 1994. [Google Scholar]
  39. Nasdala, L.; Hofmeister, W.; Norberg, N.; Martinson, J.M.; Corfu, F.; Dörr, W.; Kamo, S.L.; Kennedy, A.K.; Kronz, A.; Reiners, P.W. Zircon M257—A homogeneous natural reference material for the ion microprobe U-Pb analysis of zircon. Geostand. Geoanal. Res. 2008, 32, 247–265. [Google Scholar] [CrossRef]
  40. Nasdala, L.; Corfu, F.; Valley, J.W.; Spicuzza, M.J.; Wu, F.-Y.; Li, Q.-L.; Yang, Y.-H.; Fisher, C.; Münker, C.; Kennedy, A.K. Zircon M127–A homogeneous reference material for SIMS U–Pb geochronology combined with hafnium, oxygen and, potentially, lithium isotope analysis. Geostand. Geoanal. Res. 2016, 40, 457–475. [Google Scholar] [CrossRef] [Green Version]
  41. Stern, R.A. A New Isotopic and Trace-Element Standard for the Ion Microprobe: Preliminary Thermal Ionization Mass Spectrometry (TIMS) U-Pb and Electron-Microprobe Data; Natural Resources Canada, Geological Survey of Canada: Ottawa, ON, Canada, 2001.
  42. Santos, M.M.; Lana, C.; Scholz, R.; Buick, I.; Schmitz, M.D.; Kamo, S.L.; Gerdes, A.; Corfu, F.; Tapster, S.; Lancaster, P. A new appraisal of Sri Lankan BB zircon as a reference material for LA-ICP-MS U-Pb geochronology and Lu-Hf isotope tracing. Geostand. Geoanal. Res. 2017, 41, 335–358. [Google Scholar] [CrossRef] [Green Version]
  43. Nasdala, L.; Irmer, G.; Wolf, D. The degree of metamictization in zircons: A Raman spectroscopic study. Eur. J. Mineral. 1995, 7, 471–478. [Google Scholar] [CrossRef] [Green Version]
  44. Lenz, C.; Nasdala, L. A photoluminescence study of REE3+ emissions in radiation-damaged zircon. Am. Mineral. 2015, 100, 1123–1133. [Google Scholar] [CrossRef] [Green Version]
  45. Nasdala, L.; Akhmadaliev, S.; Artac, A.; Chanmuang, C.; Habler, G.; Lenz, C. Irradiation effects in monazite–(Ce) and zircon: Raman and photoluminescence study of Au-irradiated FIB foils. Phys. Chem. Miner. 2018, 45, 855–871. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  46. Bau, M. Controls on the fractionation of isovalent trace elements in magmatic and aqueous systems: Evidence from Y/Ho, Zr/Hf, and lanthanide tetrad effect. Contrib. Mineral. Petrol. 1996, 123, 323–333. [Google Scholar] [CrossRef]
  47. Nicholas, J.V. Origin of the luminescence in natural zircon. Nature 1967, 215, 1476. [Google Scholar] [CrossRef]
  48. Gaft, M.; Panczer, G.; Reisfeld, R.; Shinno, I. Laser-induced luminescence of rare-earth elements in natural zircon. J. Alloy. Compd. 2000, 300, 267–274. [Google Scholar] [CrossRef]
  49. Panczer, G.; De Ligny, D.; Mendoza, C.; Gaft, M.; Seydoux-Guillaume, A.M.; Wang, X. Raman and fluorescence. In Applications of Raman Spectroscopy to Earth Sciences and Cultural Heritage; Dubessy, J., Caumon, M.C., Rull, F., Eds.; EMU Notes in Mineralogy; European Mineralogical Union and Mineralogical Society of Great Britain and Ireland: London, UK, 2012; Volume 12, pp. 61–82. [Google Scholar]
  50. Friis, H.; Finch, A.A.; Williams, C.T.; Hanchar, J.M. Photoluminescence of zircon (ZrSiO4) doped with REE3+ (REE= Pr, Sm, Eu, Gd, Dy, Ho, Er). Phys. Chem. Miner. 2010, 37, 333–342. [Google Scholar] [CrossRef]
  51. Lenz, C.; Nasdala, L.; Talla, D.; Hauzenberger, C.; Seitz, R.; Kolitsch, U. Laser-induced REE3+ photoluminescence of selected accessory minerals—An “advantageous artefact” in Raman spectroscopy. Chem. Geol. 2015, 415, 1–16. [Google Scholar] [CrossRef] [Green Version]
  52. Lenz, C.; Talla, D.; Ruschel, K.; Škoda, R.; Götze, J.; Nasdala, L. Factors affecting the Nd3+ (REE3+) luminescence of minerals. Mineral. Petrol. 2013, 107, 415–428. [Google Scholar] [CrossRef] [Green Version]
  53. Geisler, T.; Burakov, B.E.; Zirlin, V.; Nikolaeva, L.; Pöml, P. A Raman spectroscopic study of high-uranium zircon from the Chernobyl lava. Eur. J. Mineral. 2005, 17, 883–894. [Google Scholar] [CrossRef]
  54. Kempe, U.; Seltmann, R.; Graupner, T.; Rodionov, N.; Sergeev, S.A.; Matukov, D.I.; Kremenetsky, A.A. Concordant U–Pb SHRIMP ages of U-rich zircon in granitoids from the Muruntau gold district (Uzbekistan): Timing of intrusion, alteration ages, or meaningless numbers. Ore Geol. Rev. 2015, 65, 308–326. [Google Scholar] [CrossRef]
  55. Kempe, U.; Trinkler, M.; Pöppl, A.; Himcinschi, C. Coloration of natural zircon. Can. Mineral. 2016, 54, 635–660. [Google Scholar] [CrossRef]
  56. Kempe, U.; Trullenque, G.; Thomas, R.; Sergeev, S.; Presnyakov, S.; Rodionov, N.; Himcinschi, C. Substitution-induced internal strain and high disorder in weakly radiation damaged hydrothermal zircon from Mt. Malosa, Malawi. Eur. J. Mineral. 2018, 30, 659–679. [Google Scholar] [CrossRef]
  57. Ruschel, K.; Nasdala, L.; Kronz, A.; Hanchar, J.M.; Többens, D.M.; Škoda, R.; Finger, F.; Möller, A. A Raman spectroscopic study on the structural disorder of monazite–(Ce). Mineral. Petrol. 2012, 105, 41–55. [Google Scholar] [CrossRef]
  58. Lenz, C.; Thorogood, G.; Aughterson, R.; Ionescu, M.; Gregg, D.J.; Davis, J.; Lumpkin, G.R. The quantification of radiation damage in orthophosphates using confocal μ-luminescence spectroscopy of Nd3+. Front. Chem. 2019, 7, 13. [Google Scholar] [CrossRef] [PubMed]
  59. Soman, A.; Geisler, T.; Tomaschek, F.; Grange, M.; Berndt, J. Alteration of crystalline zircon solid solutions: A case study on zircon from an alkaline pegmatite from Zomba–Malosa, Malawi. Contrib. Mineral. Petrol. 2010, 160, 909–930. [Google Scholar] [CrossRef]
  60. Nasdala, L.; Kronz, A.; Hanchar, J.M.; Tichomirowa, M.; Davis, D.W.; Hofmeister, W. Effects of natural radiation damage on back-scattered electron images of single crystals of minerals. Am. Mineral. 2006, 91, 1739–1746. [Google Scholar] [CrossRef]
  61. Váczi, T. A new, simple approximation for the deconvolution of instrumental broadening in spectroscopic band profiles. Appl. Spectrosc. 2014, 68, 1274–1278. [Google Scholar] [CrossRef]
  62. Elhlou, S.; Belousova, E.; Griffin, W.L.; Pearson, N.J.; O’Reilly, S.Y. Trace element and isotopic composition of GJ-red zircon standard by laser ablation. Geochim. Cosmochim. Acta Suppl. 2006, 70, A158. [Google Scholar] [CrossRef]
  63. Jackson, S.E.; Pearson, N.J.; Griffin, W.L.; Belousova, E.A. The application of laser ablation-inductively coupled plasma-mass spectrometry to in situ U–Pb zircon geochronology. Chem. Geol. 2004, 211, 47–69. [Google Scholar] [CrossRef]
  64. Griffin, W.L.; Powell, W.J.; Pearson, N.J.; O’Reilly, S.Y. GLITTER: Data reduction software for laser ablation ICP-MS. In Laser Ablation ICP-MS in the Earth Sciences: Current Practices and Outstanding Issues; Sylvester, P., Ed.; MAC Short-Course; Mineralogical Association of Canada: Quebec City, QC, Canada, 2008; pp. 308–311. [Google Scholar]
  65. Belousova, E.; Griffin, W.L.; O’Reilly, S.Y.; Fisher, N.L. Igneous zircon: Trace element composition as an indicator of source rock type. Contrib. Mineral. Petrol. 2002, 143, 602–622. [Google Scholar] [CrossRef]
  66. Zeug, M.; Nasdala, L.; Wanthanachaisaeng, B.; Balmer, W.A.; Corfu, F.; Wildner, M. Blue Zircon from Ratanakiri, Cambodia. J. Gemmol. 2018, 36, 112–132. [Google Scholar] [CrossRef] [Green Version]
  67. Nagli, L.; Bunimovich, D.; Katzir, A.; Gorodetsky, O.; Molev, V. The luminescence properties of Dy-doped high silicate glass. J. Non-Cryst. Solids 1997, 217, 208–214. [Google Scholar] [CrossRef]
  68. Nasdala, L.; Götze, J.; Hanchar, J.M.; Gaft, M.; Krbetschek, M.R. Luminescence techniques in Earth sciences. In Spectroscopic Methods in Mineralogy; Beran, A., Libowitzky, E., Eds.; EMU Notes in Mineralogy; European Mineralogical Union and Mineralogical Society of Great Britain and Ireland: London, UK, 2004; Volume 6, pp. 43–91. [Google Scholar]
  69. Trachenko, K.; Dove, M.T.; Salje, E.K. Large swelling and percolation in irradiated zircon. J. Phys. Condens. Matter 2002, 15, L1. [Google Scholar] [CrossRef] [Green Version]
  70. Carnall, W.T.; Goodman, G.L.; Rajnak, K.; Rana, R.S. A systematic analysis of the spectra of the lanthanides doped into single crystal LaF3. J. Chem. Phys. 1989, 90, 3443–3457. [Google Scholar] [CrossRef]
  71. Finch, A.A.; Garcia-Guinea, J.; Hole, D.E.; Townsend, P.D.; Hanchar, J.M. Ionoluminescence of zircon: Rare earth emissions and radiation damage. J. Phys. D Appl. Phys. 2004, 37, 2795. [Google Scholar] [CrossRef]
  72. Steely, A.N.; Hourigan, J.K.; Juel, E. Discrete multi-pulse laser ablation depth profiling with a single-collector ICP-MS: Sub-micron U–Pb geochronology of zircon and the effect of radiation damage on depth-dependent fractionation. Chem. Geol. 2014, 372, 92–108. [Google Scholar] [CrossRef]
  73. Marillo-Sialer, E.; Woodhead, J.; Hanchar, J.M.; Reddy, S.M.; Greig, A.; Hergt, J.; Kohn, B. An investigation of the laser-induced zircon ‘matrix effect’. Chem. Geol. 2016, 438, 11–24. [Google Scholar] [CrossRef]
  74. Wiemer, D.; Allen, C.M.; Murphy, D.T.; Kinaev, I. Effects of thermal annealing and chemical abrasion on ca. 3.5 Ga metamict zircon and evidence for natural reverse discordance: Insights for UPb LA-ICP-MS dating. Chem. Geol. 2017, 466, 285–302. [Google Scholar] [CrossRef] [Green Version]
  75. Allen, C.M.; Campbell, I.H. Identification and elimination of a matrix-induced systematic error in LA–ICP–MS 206Pb/238U dating of zircon. Chem. Geol. 2012, 332, 157–165. [Google Scholar] [CrossRef]
  76. Solari, L.A.; Ortega-Obregón, C.; Bernal, J.P. U–Pb zircon geochronology by LAICPMS combined with thermal annealing: Achievements in precision and accuracy on dating standard and unknown samples. Chem. Geol. 2015, 414, 109–123. [Google Scholar] [CrossRef]
  77. Shnukov, S.E.; Andreev, A.V.; Savenok, S.P. Admixture elements in zircons and apatites: A tool for provenance studies of terrigenous sedimentary rocks. In Proceedings of the European Union of Geosciences 1997, Strasbourg, France, 23–27 March 1997. [Google Scholar]
  78. Hoskin, P.W.O. Trace-element composition of hydrothermal zircon and the alteration of Hadean zircon from the Jack Hills, Australia. Geochim. Cosmochim. Acta 2005, 69, 637–648. [Google Scholar] [CrossRef]
  79. De Hoog, J.C.M.; Lissenberg, C.J.; Brooker, R.A.; Hinton, R.; Trail, D.; Hellebrand, E. Hydrogen incorporation and charge balance in natural zircon. Geochim. Cosmochim. Acta 2014, 141, 472–486. [Google Scholar] [CrossRef] [Green Version]
  80. Nasdala, L.; Beran, A.; Libowitzky, E.; Wolf, D. The incorporation of hydroxyl groups and molecular water in natural zircon (ZrSiO4). Am. J. Sci. 2001, 301, 831–857. [Google Scholar] [CrossRef] [Green Version]
Figure 1. (a) Typical laser-induced photoluminescence spectra (λexc = 473 nm) of zircon from Ratanakiri, Cambodia (Rata) and Ratnapura, Sri Lanka (M257, G3, N17). Emissions of Sm3+ and Dy3+ were identified in the yellow spectral range between 550 and 620 nm. Crystal-field split sub-level bands broaden and decrease in intensity with progressive accumulation of radiation damage from damage-free (Rata) to completely amorphised sample N17. (b) The Dy3+ emissions of partly radiation-damaged zircon samples M146 and G3 are characterised by broadened crystal-field split sub-level bands which superimpose a broad background component that resembles the shape of sample’s N17 emission of Dy3+.
Figure 1. (a) Typical laser-induced photoluminescence spectra (λexc = 473 nm) of zircon from Ratanakiri, Cambodia (Rata) and Ratnapura, Sri Lanka (M257, G3, N17). Emissions of Sm3+ and Dy3+ were identified in the yellow spectral range between 550 and 620 nm. Crystal-field split sub-level bands broaden and decrease in intensity with progressive accumulation of radiation damage from damage-free (Rata) to completely amorphised sample N17. (b) The Dy3+ emissions of partly radiation-damaged zircon samples M146 and G3 are characterised by broadened crystal-field split sub-level bands which superimpose a broad background component that resembles the shape of sample’s N17 emission of Dy3+.
Minerals 10 00083 g001
Figure 2. (a) A schematic energy-level diagram of Dy3+ (4f9 electron configuration) with its major luminescence transitions after excitation with a blue 473 nm laser (electronic level energies after Carnall et al. [70]). The free-ion levels are fine-split into so called “Stark levels” in dependence of the local crystal field. A stressed, altered, or completely destroyed structural environment has a great impact on crystal-field fine splitting, i.e., the luminescence emission. (b) Graphical representation of a molecular-dynamics simulation of a single 30 keV recoil from a decay event in zircon (modified after Trachenko et al. [69]). The arrow indicates the moving direction of the decayed daughter nuclei. Blue balls represent random substitution of REE3+ ions (i.e., Dy3+) on Zr-sites. Radiation damage alters the local structural environment of Dy3+ ions in affected areas.
Figure 2. (a) A schematic energy-level diagram of Dy3+ (4f9 electron configuration) with its major luminescence transitions after excitation with a blue 473 nm laser (electronic level energies after Carnall et al. [70]). The free-ion levels are fine-split into so called “Stark levels” in dependence of the local crystal field. A stressed, altered, or completely destroyed structural environment has a great impact on crystal-field fine splitting, i.e., the luminescence emission. (b) Graphical representation of a molecular-dynamics simulation of a single 30 keV recoil from a decay event in zircon (modified after Trachenko et al. [69]). The arrow indicates the moving direction of the decayed daughter nuclei. Blue balls represent random substitution of REE3+ ions (i.e., Dy3+) on Zr-sites. Radiation damage alters the local structural environment of Dy3+ ions in affected areas.
Minerals 10 00083 g002
Figure 3. Graphical fitting results that estimate the amorphous fraction fa from the laser-induced PL emission of Dy3+ (4F9/26H13/2) using two different least-square deconvolution procedures on the example of zircon samples Rata, M146 and G3. (a) Individual crystal-field split sub-level bands are modulated with individual Gauss-Lorentz functions with spectral position, FWHM and amplitude treated as free variables that are optimised to fit best with experimental data. The amorphous component is modelled and best described using a Gaussian function with a variable amplitude, but with fixed FWHM of 540 cm−1 and peak position at 17,360 cm−1 (~576 nm). (b) Classic least-square (CLS) fitting procedure that takes advantage of multiple references, so-called “loadings”. Those comprise spectra that are considered to represent non-radiation damaged, crystalline samples (e.g., sample Rata or other annealed zircon samples) and a reference from a completely amorphised sample (N17, cf., Figure 1a). Note that band-broadening is accounted for by using several residue spectra of partly-radiation damaged samples with the amorphous component subtracted according to the fitting presented in (a).
Figure 3. Graphical fitting results that estimate the amorphous fraction fa from the laser-induced PL emission of Dy3+ (4F9/26H13/2) using two different least-square deconvolution procedures on the example of zircon samples Rata, M146 and G3. (a) Individual crystal-field split sub-level bands are modulated with individual Gauss-Lorentz functions with spectral position, FWHM and amplitude treated as free variables that are optimised to fit best with experimental data. The amorphous component is modelled and best described using a Gaussian function with a variable amplitude, but with fixed FWHM of 540 cm−1 and peak position at 17,360 cm−1 (~576 nm). (b) Classic least-square (CLS) fitting procedure that takes advantage of multiple references, so-called “loadings”. Those comprise spectra that are considered to represent non-radiation damaged, crystalline samples (e.g., sample Rata or other annealed zircon samples) and a reference from a completely amorphised sample (N17, cf., Figure 1a). Note that band-broadening is accounted for by using several residue spectra of partly-radiation damaged samples with the amorphous component subtracted according to the fitting presented in (a).
Minerals 10 00083 g003
Figure 4. Photoluminescence hyperspectral maps of two large, heterogeneous Sri Lankan zircon single-crystals, SL1 (a) and SL2 (b). The crystals were cut in halves, with one half left in the original, radiation-damaged state and the other half annealed at 1300 °C for 3 days (see inset sketch). Mapping parameter is the amorphous fraction fa as determined with the PL method (cf., Figure 3b) using a software-based least-squares fitting applied to a point-by-point scanned data-set. Radiation damage is detected in the original halves (left) with high fa in zones with elevated U and Th concentrations (cf., Table A2 in Appendix A). The annealed halves (right) fully recovered crystallinity and show no amorphous component. Pairs of LA-ICP-MS analyses-spots on each half are indicated.
Figure 4. Photoluminescence hyperspectral maps of two large, heterogeneous Sri Lankan zircon single-crystals, SL1 (a) and SL2 (b). The crystals were cut in halves, with one half left in the original, radiation-damaged state and the other half annealed at 1300 °C for 3 days (see inset sketch). Mapping parameter is the amorphous fraction fa as determined with the PL method (cf., Figure 3b) using a software-based least-squares fitting applied to a point-by-point scanned data-set. Radiation damage is detected in the original halves (left) with high fa in zones with elevated U and Th concentrations (cf., Table A2 in Appendix A). The annealed halves (right) fully recovered crystallinity and show no amorphous component. Pairs of LA-ICP-MS analyses-spots on each half are indicated.
Minerals 10 00083 g004
Figure 5. (a) The calculated α-dose of various Sri Lankan zircon samples (cf., Table A1 and Table A2 in Appendix A) plotted against the amorphous fraction fa determined using the PL method (cf., Figure 3b). Within analytical uncertainty (2σ), data points lie on the regression curve (red) given by Zhang and Salje [35], who applied the direct-impact model to data of Sri Lankan zircon obtained from various alternative, bulk methods. (b) Obtained concordant U-Pb ages plotted against U + Th concentrations from multiple measurement spots on the annealed and unannealed halves of zircon samples SL1, SL2 and MalawiX. Linear regressions of data points from radiation-damaged, unannealed halves are indicated. Note, that bias of 206Pb/238U ages is present irrespective of the radiation damage accumulated (fa indicated for selected measurement spots).
Figure 5. (a) The calculated α-dose of various Sri Lankan zircon samples (cf., Table A1 and Table A2 in Appendix A) plotted against the amorphous fraction fa determined using the PL method (cf., Figure 3b). Within analytical uncertainty (2σ), data points lie on the regression curve (red) given by Zhang and Salje [35], who applied the direct-impact model to data of Sri Lankan zircon obtained from various alternative, bulk methods. (b) Obtained concordant U-Pb ages plotted against U + Th concentrations from multiple measurement spots on the annealed and unannealed halves of zircon samples SL1, SL2 and MalawiX. Linear regressions of data points from radiation-damaged, unannealed halves are indicated. Note, that bias of 206Pb/238U ages is present irrespective of the radiation damage accumulated (fa indicated for selected measurement spots).
Minerals 10 00083 g005
Figure 6. Comparison of spectral parameters of the Raman ν3(SiO4) and laser-induced Dy3+ photoluminescence sublevel-bands of Sri Lankan zircon samples that are characterised by different degrees of radiation-damage accumulation. Measurement points from the unannealed halves of sample SL1 and SL2 are indicated in Figure 4. Photoluminescence sublevel-bands considered here are indicated as Dy1 and Dy2 in Figure 3a. Note, that Dy3+ sublevel-band parameters were fitted with Gauss-Lorentz functions after subtraction of the ‘amorphous’ broad-band component (cf., Figure 3). (a) The band-width of PL sublevel Dy1 closely correlates with that obtained from the Raman ν3(SiO4) mode. Both, increasingly broaden with progressive accumulation of structural radiation damage. (b) Spectral band positions are inversely correlated. Whereas Raman band ν3 shifts to lower wavenumbers, PL bands Dy1 and Dy2 shift to higher peak positions with progressive damage accumulation.
Figure 6. Comparison of spectral parameters of the Raman ν3(SiO4) and laser-induced Dy3+ photoluminescence sublevel-bands of Sri Lankan zircon samples that are characterised by different degrees of radiation-damage accumulation. Measurement points from the unannealed halves of sample SL1 and SL2 are indicated in Figure 4. Photoluminescence sublevel-bands considered here are indicated as Dy1 and Dy2 in Figure 3a. Note, that Dy3+ sublevel-band parameters were fitted with Gauss-Lorentz functions after subtraction of the ‘amorphous’ broad-band component (cf., Figure 3). (a) The band-width of PL sublevel Dy1 closely correlates with that obtained from the Raman ν3(SiO4) mode. Both, increasingly broaden with progressive accumulation of structural radiation damage. (b) Spectral band positions are inversely correlated. Whereas Raman band ν3 shifts to lower wavenumbers, PL bands Dy1 and Dy2 shift to higher peak positions with progressive damage accumulation.
Minerals 10 00083 g006
Figure 7. (a) PL hyperspectral maps of two zircon crystals from pegmatites of the Mt. Malosa massif, Malawi (see inset photographs) that were cut in half, with one half annealed at 1300 °C for three days (cf., Figure 4). The amorphous fraction fa determined according to the deconvolution procedure of the Dy3+ photoluminescence as described above is given colour-coded. Note, that the scale may vary, e.g., to highlight subtle differences in radiation damage accumulation within sample MalawiA. Spots of LA-ICP-MS, detailed Raman and PL measurements are indicated. (b) U-Pb concordia and results of LA-ICP-MS isotope ratios obtained from measurement spots indicated in (a).
Figure 7. (a) PL hyperspectral maps of two zircon crystals from pegmatites of the Mt. Malosa massif, Malawi (see inset photographs) that were cut in half, with one half annealed at 1300 °C for three days (cf., Figure 4). The amorphous fraction fa determined according to the deconvolution procedure of the Dy3+ photoluminescence as described above is given colour-coded. Note, that the scale may vary, e.g., to highlight subtle differences in radiation damage accumulation within sample MalawiA. Spots of LA-ICP-MS, detailed Raman and PL measurements are indicated. (b) U-Pb concordia and results of LA-ICP-MS isotope ratios obtained from measurement spots indicated in (a).
Minerals 10 00083 g007
Figure 8. Results of LA-ICP-MS trace-element analyses of zircon samples from Sri Lanka (SL2) and Mt. Malosa pegmatites, Malawi obtained from their annealed crystal halves. Detailed results are given in Table A7, Table A8 and Table A9. (a) Plots of mean chondrite-normalised REE concentrations of sample SL2 (n = 13), MalawiX (n = 14), and MalawiA (n = 14) with their limits given colour-shaded, respectively. (b) Zircon provenance discrimination diagram after Shnukov [77] using Hf and Y concentrations. Expected fields of source rock type are indicated (details see [65]).
Figure 8. Results of LA-ICP-MS trace-element analyses of zircon samples from Sri Lanka (SL2) and Mt. Malosa pegmatites, Malawi obtained from their annealed crystal halves. Detailed results are given in Table A7, Table A8 and Table A9. (a) Plots of mean chondrite-normalised REE concentrations of sample SL2 (n = 13), MalawiX (n = 14), and MalawiA (n = 14) with their limits given colour-shaded, respectively. (b) Zircon provenance discrimination diagram after Shnukov [77] using Hf and Y concentrations. Expected fields of source rock type are indicated (details see [65]).
Minerals 10 00083 g008
Figure 9. (a) Raman ν3(SiO4) band widths of the Malawian samples (X, A) plotted against the α-dose calculated from the age of the samples and the U and Th concentrations (Table A2) obtained from multiple LA-ICP-MS measurement spots (indicated in Figure 7). (b) The amorphous fraction fa as determined from the PL of Dy3+ from same measurement spots. The correlation fa vs. α-dose of Sri Lankan zircon according to the direct impact model after Zhang and Salje [35] is given as red curve (cf., Figure 5a). An additional (green) curve represents this correlation with Sri Lankan α-doses corrected with a factor of ×0.55 (Ba = 4.96 × 10−19 g) that accounts for geothermal annealing and represents the ’effective’ dose Sri Lankan zircon samples have experienced.
Figure 9. (a) Raman ν3(SiO4) band widths of the Malawian samples (X, A) plotted against the α-dose calculated from the age of the samples and the U and Th concentrations (Table A2) obtained from multiple LA-ICP-MS measurement spots (indicated in Figure 7). (b) The amorphous fraction fa as determined from the PL of Dy3+ from same measurement spots. The correlation fa vs. α-dose of Sri Lankan zircon according to the direct impact model after Zhang and Salje [35] is given as red curve (cf., Figure 5a). An additional (green) curve represents this correlation with Sri Lankan α-doses corrected with a factor of ×0.55 (Ba = 4.96 × 10−19 g) that accounts for geothermal annealing and represents the ’effective’ dose Sri Lankan zircon samples have experienced.
Minerals 10 00083 g009

Share and Cite

MDPI and ACS Style

Lenz, C.; Belousova, E.; Lumpkin, G.R. The In-Situ Quantification of Structural Radiation Damage in Zircon Using Laser-Induced Confocal Photoluminescence Spectroscopy. Minerals 2020, 10, 83. https://doi.org/10.3390/min10010083

AMA Style

Lenz C, Belousova E, Lumpkin GR. The In-Situ Quantification of Structural Radiation Damage in Zircon Using Laser-Induced Confocal Photoluminescence Spectroscopy. Minerals. 2020; 10(1):83. https://doi.org/10.3390/min10010083

Chicago/Turabian Style

Lenz, Christoph, Elena Belousova, and Gregory R. Lumpkin. 2020. "The In-Situ Quantification of Structural Radiation Damage in Zircon Using Laser-Induced Confocal Photoluminescence Spectroscopy" Minerals 10, no. 1: 83. https://doi.org/10.3390/min10010083

APA Style

Lenz, C., Belousova, E., & Lumpkin, G. R. (2020). The In-Situ Quantification of Structural Radiation Damage in Zircon Using Laser-Induced Confocal Photoluminescence Spectroscopy. Minerals, 10(1), 83. https://doi.org/10.3390/min10010083

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop