Introduction
Niobium and tantalum are critical metals that have been used increasingly in advanced technologies (e.g. information technology, energy, aerospace, national defence and military, Linnen et al., Reference Linnen, Van Lichtervelde and Černý2012). Members of the columbite–tantalite solid-solution series are a major mineral resource of Nb and Ta, and have the formula, AB2O6, where A is occupied by Fe2+ and Mn2+ (Mg2+ in rare cases) and B is dominated by Nb5+ and Ta5+ (Černý and Ercit, Reference Černý, Ercit, Möller, Černý and Saupé1989; Ercit, Reference Ercit1994; Romer et al., Reference Romer, Smeds and Ern1996; Melcher et al., Reference Melcher, Graupner, Gäbler, Sitnikova and Dewaele2015), and belong to the orthorhombic crystal system. There are three approved end-members in the columbite group: columbite-(Fe) [Fe2+Nb2O6]; columbite-(Mg) [Mg2+Nb2O6]; and columbite-(Mn) [Mn2+Nb2O6]. In this paper the shortened form ‘columbite’ refers to all columbite end-members. This mineral group typically occurs in highly fractionated granites and granitic pegmatites (Linnen et al., Reference Linnen, Van Lichtervelde and Černý2012; Feng et al., Reference Feng, Liang, Zhang, Wang, Zhou, Yang, Gao, Wang and Ding2019, Reference Feng, Liang, Linnen, Zhang, Zhou, Zhang and Gao2020).
In recent years, columbite has become an increasingly important mineral for U–Pb geochronology, and can provide robust age constraints for highly evolved granites and rare-element pegmatites (Romer and Smeds, Reference Romer and Smeds1994; Smith et al., Reference Smith, Foster, Romer, Tindle, Kelley, Noble, Horstood and Breaks2004; Baumgartner et al., Reference Baumgartner, Romer, Moritz, Sallet and Chiaradia2006; Melleton et al., Reference Melleton, Gloaguen, Frei, Novák and Breiter2012; Deng et al., Reference Deng, Li, Zhao, Hu, Hu, Selby and Souza2013; Che et al., Reference Che, Wu, Wang, Gerdes, Ji, Zhao, Yang and Zhu2015; Melcher et al., Reference Melcher, Graupner, Gäbler, Sitnikova, Oberthür, Gerdes, Badanina and Chudy2017; Lupulescu et al., Reference Lupulescu, Chiarenzelli, Pecha, Singer and Regan2018; Yan et al., Reference Yan, Qiu, Wang, Wang, Wei, Li, Zhang, Li and Liu2018; Zhou et al., Reference Zhou, Qin, Tang, Wang and Sakyi2018; Feng et al., Reference Feng, Liang, Zhang, Wang, Zhou, Yang, Gao, Wang and Ding2019, Reference Feng, Liang, Linnen, Zhang, Zhou, Zhang and Gao2020). However, analogous to zircon, the U–Pb isotopic ages of columbite are also discordant, and the principal reason is metamictisation. Although the metamictisation of zircon, the conventional mineral for age determination, has been studied intensely (Nasdala et al., Reference Nasdala, Hanchar, Kronz and Whitehouse2005, and references therein), that of columbite has received less interest (e.g. Lumpkin, Reference Lumpkin1998; Feng et al., Reference Feng, Liang, Linnen, Zhang, Zhou, Zhang and Gao2020).
Raman spectra reflect the inelastic Raman scattering of light that causes shifts in the frequency of the excitation laser. These frequency shifts depend solely on the inherent vibrational and rotational energy-level structure of the molecules interacting with the excitation laser. Different substances, therefore, have characteristic Raman spectra reflecting their unique molecular structures. Raman spectroscopy has many advantages, including small detection limits, simple sample preparation, and fast, economic, and non-destructive data acquisition (Nasdala et al., Reference Nasdala, Smith, Kaindl and Zieman2004). It has been used widely for quantifying the degree of metamictisation of zircon in studies such as those of Nasdala et al. (Reference Nasdala, Irmer and Wolf1995), Zhang et al. (Reference Zhang, Salje, Farnan, Graeme-Barber, Daniel, Ewing, Clark and Leroux2000), Balan et al. (Reference Balan, Neuville, Trocellier, Fritsch, Muller and Calas2001) and Nasdala et al. (Reference Nasdala, Beran, Libowitzky and Wolf2001b, Reference Nasdala, Smith, Kaindl and Zieman2004, Reference Nasdala, Hanchar, Kronz and Whitehouse2005). These investigations demonstrated that the band position and full width at half maximum (FWHM) of a characteristic band of zircon [such as υ3(SiO4)] changes systematically with the degree of metamictisation, and that this can be used to indicate quantitatively the degree of zircon metamictisation. At present, there are few Raman spectroscopy studies on columbite. Husson et al. (Reference Husson, Repelin, Dao and Brusset1977) analysed the Nb–O and M–O bonds of synthetic columbite-type niobates MNb2O6 (where M = Mg, Ca, Mn, Fe, Co, Ni, Cu, Zn and Cd). Moreira et al. (Reference Moreira, Rubinger, Krambrock and Dias2010) conducted a systematic Raman spectroscopic study on a single, highly crystalline columbite-(Mn) crystal. Their results demonstrated that the Raman bands below 150 cm–1 and those at 250–380 cm–1 were related to A-site cations, whereas the other bands were more likely to be controlled by B–O bonds. The positions of strong Raman bands were influenced mainly by the presence of Nb and Ta in the B site (Husson et al., Reference Husson, Repelin, Dao and Brusset1977). Moreover, Moreira et al. (Reference Moreira, Rubinger, Krambrock and Dias2010) demonstrated that the strongest band (A 1g vibration mode) occurred at ~882 cm–1. Feng et al. (Reference Feng, Liang, Linnen, Zhang, Zhou, Zhang and Gao2020) conducted an integrated Raman spectroscopic, compositional, and U–Pb geochronological study on two columbite crystals from No. 1 pegmatite dyke at Dakalasu and discussed the relationship between metamictisation and age concordance. The relationship between metamictisation and Raman spectra, however, was not discussed.
In this contribution, we systematically investigated the Raman spectra of a representative suite of columbite samples from rare-element pegmatites in northwest China. The major–minor-element concentrations and U–Pb isotopic ratios of these samples were determined in a previous study (Feng et al., Reference Feng, Liang, Zhang, Wang, Zhou, Yang, Gao, Wang and Ding2019, Reference Feng, Liang, Linnen, Zhang, Zhou, Zhang and Gao2020, Reference Feng, Liang, Lei, Ju, Zhang, Gao, Zhou and Wu2021). The factors influencing the key parameters [i.e. position (p) and FWHM] of the strongest Raman band of columbites are explored, and the relationship between these parameters and the degree of metamictisation of columbite is discussed.
Previous studies have shown that, in common with zircon U–Pb isotopic ages, the U–Pb isotopic ages of columbite can also become discordant with metamictisation being one of the principal reasons (Romer, Reference Romer2003; Feng et al., Reference Feng, Liang, Zhang, Wang, Zhou, Yang, Gao, Wang and Ding2019). Therefore, an in-depth understanding of the metamictisation of columbite is of considerable importance for the rational selection of analytical points, avoidance of invalid analyses and interpretation of columbite U–Pb isotopic ages obtained using in situ radiometric age determination techniques.
Sample descriptions
Six columbite samples were investigated in this study. Sample JR-1 was collected from the intermediate zone of the Jing'erquan No. 1 spodumene pegmatite located in the Eastern Tianshan, Xinjiang, China. Details of the regional and deposit geology are described by Yao et al. (Reference Yao, Xu, Yang, Wu and Geng2020), Liu et al. (Reference Liu, Wang, Jeon, Hou, Xue, Zhou, Chen, Zhang and Xi2020) and Feng et al. (Reference Feng, Liang, Lei, Ju, Zhang, Gao, Zhou and Wu2021); information pertinent to this study are summarised below. This Li–Nb–Ta-mineralised pegmatite dyke occurs in the early Carboniferous Xiaorequanzi–Wutongwozi belts in the eastern section of the Quoltag Arc (Liu et al., Reference Liu, Wang, Jeon, Hou, Xue, Zhou, Chen, Zhang and Xi2020). The ore minerals hosted by the pegmatite dyke mainly consist of spodumene, columbite-group minerals and lepidolite. Aggregates of green, scaly mica occur in the middle of the pegmatite and are associated with coarse-grained microcline. The grain size of columbite in sample JR-1 (Fig. 1) is in the range of 1–2 cm; it is associated mainly with spodumene, albite, quartz and phosphate minerals. The columbite-(Mn) is compositionally homogeneous, with an average #Mn [= molar Mn/(Mn + Fe)] of 0.84 and #Ta [= molar Ta/(Nb + Ta)] of 0.08 (Feng et al., Reference Feng, Liang, Lei, Ju, Zhang, Gao, Zhou and Wu2021). The U–Pb Concordia age obtained for this sample is 250.8 ± 1.0 Ma (Feng et al., Reference Feng, Liang, Lei, Ju, Zhang, Gao, Zhou and Wu2021).
Photos of columbite samples from the investigated pegmatites: JR-1; DKLS 107; DKLS 108; K802 (the first two rows of grains); K803 (the two lower rows) and K650. In samples JR-1 to DKLS 108, ‘a’, ‘b’ and ‘c’ represent the long, middle and short axes of the crystals, respectively, which were determined using electron back-scatter diffraction [EBSD]).
Samples DKLS 107 and DKLS 108 were collected from the intermediate zone of the No. 1 pegmatite dyke at Dakalasu, which is located in the Altai orogenic belt, Xinjiang, China (Feng et al., Reference Feng, Liang, Linnen, Zhang, Zhou, Zhang and Gao2020). To date, the No. 1 pegmatite is the largest pegmatite dyke in the area and is characterised by significant Be, Nb, and Ta mineralisation. The pegmatite consists of a medium- to coarse-grained graphic texture zone and a megacrystic zone with miarolitic cavities (Wang et al., Reference Wang, Chen and Xu2003; Zou and Li, Reference Zou and Li2006). Both samples are single crystals several centimetres in size (Fig. 1). They are compositionally columbite-(Mn), with #Mn ratios of 0.79 and 0.72, and #Ta of 0.45 and 0.19, for DKLS 107 and DKLS 108, respectively (Feng et al., Reference Feng, Liang, Linnen, Zhang, Zhou, Zhang and Gao2020). The two samples have comparable #Mn ratios, but distinctly different #Ta, indicative of significant Nb–Ta fractionation and limited Fe–Mn fractionation (Feng et al., Reference Feng, Liang, Linnen, Zhang, Zhou, Zhang and Gao2020).
Samples K802 and K803 were collected from the intermediate zones of the 802 and 803 spodumene pegmatite dykes, respectively, of Kalu'an, which is located in the Kalu'an–Azubai pegmatite field in the Altai orogenic belt (Feng et al., Reference Feng, Liang, Zhang, Wang, Zhou, Yang, Gao, Wang and Ding2019). Exposed strata in the area are dominated by the Middle–Upper Silurian Kulumuti Group (Ma et al., Reference Ma, Zhang, Tang, Lv, Zhang and Zhao2015). Both samples have been prepared as mounts of fine-grained mineral separates (>200 grains, Fig. 1) with the composition of columbite-(Mn). The #Mn ratios of the core, mantle and rim of individual grains in samples K802 and K803 are relatively constant, ranging from 0.90 to 0.95. The #Ta ratios of individual mineral grains, however, is highly variable. For sample K802, the average #Ta of the mantle of columbite grains (0.43) is notably higher than that of the core (0.14) and rim (0.20). Similarly, for sample K803, the average #Ta ratio of the Ta-rich mantle (0.49) is notably higher than that of the core (0.18), the Nb-rich mantle (0.18), and the rim (0.18) (Feng et al., Reference Feng, Liang, Zhang, Wang, Zhou, Yang, Gao, Wang and Ding2019).
Sample K650 is a mount comprising >200 columbite grains (Fig. 1) separated from the intermediate zone of the No. 650 pegmatite dyke of Kalu'an. The major-element compositions of these columbite grains were determined in the present investigation.
All of the separated columbite grains from each of the samples are columbite-(Mn). The compositional (electron microprobe) data of the columbite grains reported in previous studies (Feng et al., Reference Feng, Liang, Zhang, Wang, Zhou, Yang, Gao, Wang and Ding2019, Reference Feng, Liang, Linnen, Zhang, Zhou, Zhang and Gao2020, Reference Feng, Liang, Lei, Ju, Zhang, Gao, Zhou and Wu2021) are compiled in Supplementary Table S1.
Analytical methods
Back-scattered electron imaging and electron back-scatter diffraction
Back-scattered electron images were obtained using a FEI Quanta 650 Environmental Scanning Electron Microscope at the Xi'an Key Laboratory for Mineralization and Efficient Utilization of Critical Metals, Chang'an University to characterise the internal structure of columbite grains. The accelerating voltage was adjusted between 15–20 kV to optimise image quality. The spot size of the electron beam was 3.9 μm.
Electron back-scatter diffraction (EBSD) analysis of single-crystal samples JR-1, DKLS 107 and DKLS 108 was completed at the State Key Laboratory of Continental Dynamics, Northwestern University using a FEI Quanta FEG 450 Field Emission Environmental Scanning Electron Microscope equipped with a Nordlys Nano electron back-scatter diffractometer produced by Oxford Instruments. The surface of the samples was polished prior to analysis. During analyses, the samples were inclined at 70° in the sample chamber, which was maintained under low vacuum (~40 Pa). The accelerating voltage was 20 kV and the beam spot size was 7.0 μm. Following the collection of a background signal, measurements were conducted on a series of spots from core to rim on individual columbite crystals to determine the orientation of the crystal lattice and obtain diffraction patterns. Calibration of the Kikuchi pattern was accomplished automatically using the instrument.
Major-element analysis
The major-element composition of columbite grains in sample K650 was determined at the Xi'an Key Laboratory for Mineralization and Efficient Utilization of Critical Metals, Chang'an University using a JEOL JXA-8100 electron microprobe (EPMA). The accelerating voltage and electron beam current were 15 kV and 20 nA, respectively. During EPMA, the peak counting time was 30 s for Sn and F, 10 s for W, Sc, Nb and Ta, and 15 s for Mn, Fe, Ca, Ti and Al. The compositions and formulae of the analysed columbite grains were calculated on the basis of six oxygen atoms. Details of the analytical parameters and the standards used for instrument calibration are summarised in Supplementary Table S2.
Raman spectroscopy
To characterise the relationship between metamictisation of columbite and their key Raman bands, the Raman spectra of the samples were determined at the Xi'an Key Laboratory for Mineralization and Efficient Utilization of Critical Metals, Chang'an University, using a LabRam HR evolution 800 mm confocal laser Raman spectrometer manufactured by HORIBA. The wavelength and energy of the excitation laser were set to 532 nm and 100 mW, respectively. The laser beam was focused on the sample surface through a 50× objective lens. A 100 μm pinhole and a 50% attenuation filter were used for controlling the signal intensity. A grating with a groove density of 600 grooves/mm was used for light dispersion. Spectra were acquired over a wavenumber range of 0–1200 cm–1 and represent the average of three accumulations of 28 s each. The Labspec 6 software provided by HORIBA was used for baseline correction and smoothing of Raman spectra, calculation of band positions and FWHM values. Spectral resolution was ~1.2 cm–1. Spectral fitting and correction for experimental band broadening used the method of Váczi and Tamás (2014).
Laser ablation inductively coupled plasma mass spectrometry
Laser ablation inductively coupled plasma mass spectrometry (LA–ICP–MS) U–Pb age determination of sample JR-1 was performed using a Thermo Fisher Scientific iCAP-Q model ICP–MS coupled to a RESOlution S155 193nm ArF excimer laser ablation system at the State Key Laboratory for Mineral Deposits Research, Nanjing University. The diameter of the laser beam was set to 43 μm. The laser output energy and repetition rate were 105 mJ and 4 Hz, respectively. Coltan 139 (with a reference age of 506.2 ± 5.0 Ma, Che et al., Reference Che, Wu, Wang, Gerdes, Ji, Zhao, Yang and Zhu2015) was used as the external standard for calibrating U, Th and Pb isotopic ratios. Zircon standard samples 91500 and PL (Plešovice) were used as secondary standards to monitor U, Th and Pb isotope signals (Feng et al., Reference Feng, Liang, Lei, Ju, Zhang, Gao, Zhou and Wu2021). For each spot analysis, the signal of the gas background was collected for 20 s, followed by 50 s signal collection time on the standards or the unknown samples. The dwell time for 204Pb, 206Pb, and 208Pb was 15 ms, for 207Pb was 30 ms, for 232Th and 238U was 10 ms, and for 202Hg was 6 ms. The concentrations of the rare earth elements (REE) were determined simultaneously with U, Th and Pb. Each experiment consisted of the following sequence of analyses: one analysis of NIST SRM 610, two analyses of 91500, two analyses of Coltan 139, and one analysis of PL, followed by six sample analyses. For a detailed description of the analytical procedure, see Che et al. (Reference Che, Wu, Wang, Gerdes, Ji, Zhao, Yang and Zhu2015). The 206Pb/238U, 207Pb/235U, 207Pb/206Pb and 208Pb/232Th ratios and U–Pb radiometric age were calibrated using Iolite v3.7 (Hellstrom et al., Reference Hellstrom, Paton, Woodhead, Hergt and Sylvestor2008) and Isoplot 3.0 (Ludwig, Reference Ludwig2003). Due to the low common Pb (204Pb) concentrations in the columbite samples, no correction was applied for common Pb.
Results
Electron back-scattered diffraction
Back-scattered electron (BSE) images of representative columbite grains in samples K802 and K803 are shown in Fig. 2. The BSE images illustrate that the columbite grains in sample K802 are compositionally zoned, with the core and rim generally coloured grey, and the mantle being relatively bright (e.g. K802-1, K802-2, K802-13 and K802-9 in Fig. 2). Compared to sample K802, columbite grains in sample K803 exhibit oscillatory zonation (e.g. K803-3, K803-15 and K803-16 in Fig. 2).
Representative back-scattered electron images of columbite in samples K802 and K803 from the Kalu'an pegmatite (modified after Feng et al., Reference Feng, Liang, Zhang, Wang, Zhou, Yang, Gao, Wang and Ding2019). The red dashed circles represent the location of laser ablation analyses, whereas the red solid dots indicate the position of electron microprobe and Raman spectroscopic analyses.
The crystallographic orientations of the coarse-grained columbite crystals in samples JR-1, DKLS 107, and DKLS 108 (Fig. 1) are illustrated in Fig. 3. The a, b and c axes represent the long, medium and short axis directions, respectively, corresponding to the a, b and c axes in Fig. 1.
Crystallographic orientations of columbite grains in samples JR-1, DKLS 107 and DKLS 108 obtained using EBSD. The brown balls represent divalent cations of Fe and Mn in the A-site; the blue balls denote pentavalent cations Nb and Ta in the B-site; and the red balls represent O anions. The crystal structures of the columbite are modified after Tarantino and Zema (Reference Tarantino and Zema2005).
Electron probe microanalysis
The compositional data from previous work on samples K802, K803, DKLS 107, DKLS 108 and JR-1 by Feng et al. (Reference Feng, Liang, Zhang, Wang, Zhou, Yang, Gao, Wang and Ding2019, Reference Feng, Liang, Linnen, Zhang, Zhou, Zhang and Gao2020, Reference Feng, Liang, Lei, Ju, Zhang, Gao, Zhou and Wu2021) are compiled in Supplementary Table S1. Twenty-eight analyses were acquired here for further analysis of columbite grains in sample K650. The analytical results are provided in Table 1. Columbite grains in this sample are compositionally columbite-(Mn). The grains are compositionally homogeneous, with no obvious change in element concentrations from core to rim.
Major-element content of columbite in sample K650.*
*Notes: n.d. = not detected; S.D. = standard deviation; #Mn = molar Mn/(Mn + Fe) and #Ta =molar Ta/(Nb + Ta).
Raman spectroscopy
Totals of 92, 20, 20, 38, 32 and 49 Raman spectra were acquired for columbite grains in samples JR-1, DKLS 107, DKLS 108, K802, K803 and K650, respectively. The strongest Raman band (i.e. the A 1g vibration mode related to the B–O bond) for columbite-(Mn) from the No. 1 pegmatite at Jing'erquan (JR-1) occurs at 866.5–872.7 cm–1, with corresponding FWHM values of 21.1–30.2 cm–1. The strongest Raman band for columbite-(Mn) from the No. 107 and No. 108 pegmatites at Dakalasu (DKLS 107 and DKLS 108) occur at 852.3–866.2 cm–1 and 854.1–864.8 cm–1, respectively, with corresponding FWHM values of 26.5–54.1cm–1 and 32.7–50.7cm–1. The strongest Raman band of columbite-(Mn) in the Kalu'an pegmatite (K802, K803 and K650) are 873.1–878.2 cm–1, 870.1–877.9 cm–1 and 871.5–878.3 cm–1, respectively, with corresponding FWHM values of 13.4–20.4 cm–1, 13.4–26.3 cm–1 and 10.1–26.9 cm–1, respectively. Details of the key parameters of the characteristic Raman bands in these samples are provided in Supplementary Table S3. Representative Raman spectra of columbite-(Mn) from each pegmatite are shown in Fig. 4, with information on their Raman characteristic band (band position and FWHM) provided in Table 2. The position and FWHM of the characteristic Raman band for the columbite-(Mn) grains have negative correlation (Fig. 5).
Representative Raman spectra for the columbite samples. All spectra show the strongest Raman band occurring at ~850–880 cm–1.
Bivariate diagram illustrating the correlation between the position and FWHM of the characteristic Raman band for columbite (A 1g vibration mode related to B–O bond) from all of the pegmatite samples.
Representative positions and FWHM values of the characteristic Raman band of columbite (A 1g vibration mode related to the B–O bond) in each pegmatite.*
*Notes: #Mn = molar Mn/(Mn + Fe) and #Ta = molar Ta/(Nb + Ta); p = band position.
U–Pb geochronology
Feng et al. (Reference Feng, Liang, Zhang, Wang, Zhou, Yang, Gao, Wang and Ding2019, Reference Feng, Liang, Linnen, Zhang, Zhou, Zhang and Gao2020, Reference Feng, Liang, Lei, Ju, Zhang, Gao, Zhou and Wu2021) acquired the U–Pb radiometric ages of columbite grains for samples JR-1, DKLS 107, DKLS 108, K802, and K803; the U–Pb isotopic data are compiled in Supplementary Table S4. An additional 39 analyses were collected here for further analysis for columbite-(Mn) sample JR-1 to establish a complete dataset including Raman spectra and columbite composition data. The spots analysed were distributed approximately evenly between the core, mantle and rim of the columbite-(Mn) grain in regions devoid of inclusions and fractures. Uranium–Th concentrations, isotope ratios, and ages are provided in Table 3. Average U and Th concentrations of 39 points are 1454.3 ± 279.0 ppm and 17.3 ± 5.6 ppm, respectively. The 206Pb/238U and 207Pb/235U ratios are 0.0398 ± 0.0018 and 0.2902 ± 0.0428, respectively. The apparent 206Pb/238U ages range from 230.9 to 304.5 Ma, with the majority of data points having a 206Pb/238U age of ~250 Ma, which is consistent with the age reported by Feng et al. (Reference Feng, Liang, Lei, Ju, Zhang, Gao, Zhou and Wu2021).
U–Pb geochronological results for sample JR-1 obtained using LA–ICP–MS.
Alpha-decay dose and displacement per atom calculations
Alpha-decay dose (D) and displacement per atom (dpa) are two important parameters used to characterise the degree of metamictisation of columbite (Lumpkin and Ewing, Reference Lumpkin and Ewing1988; Murakami et al., Reference Murakami, Chakoumakos, Ewing, Lumpkin and Weber1991; Davis and Krogh, Reference Davis and Krogh2000; Lumpkin et al., Reference Lumpkin, Leung and Ferenczy2012; Feng et al., Reference Feng, Liang, Linnen, Zhang, Zhou, Zhang and Gao2020). The alpha-decay dose (D) can be calculated using the equation defined by Holland and Gottfried (Reference Holland and Gottfried1955), Lumpkin and Ewing (Reference Lumpkin and Ewing1988), and Lumpkin et al. (Reference Lumpkin, Leung and Ferenczy2012):
where D is the α-decay dose in α-events/mg, N238, N235 and N232 are the number of atoms of 238U, 235U (equal to 238U/137.88) and 232Th in 1 mg of columbite, respectively, and τ238, τ235 and τ232 are the mean lives (equal to 1/λ) of 238U, 235U and 232Th, respectively. The age (t in Ma) utilised here is based on the U–Pb geochronological results. In per α-decay event, ~1500 displaced atoms are produced and the corresponding dpa values are calculated using the following formula (Matzke, Reference Matzke1982; Weber et al., Reference Weber, Turcotte and Roberts1982; Weber and Roberts, Reference Weber and Roberts1983; Van Konynenburg and Guinan, Reference Van Konynenburg and Guinan1983; Vance et al., Reference Vance, Kariorsis, Cartz, Wong, Wicks and Ross1984; Eyal and Fleischer, Reference Eyal and Fleischer1985; Lumpkin and Ewing, Reference Lumpkin and Ewing1988; Feng et al., Reference Feng, Liang, Linnen, Zhang, Zhou, Zhang and Gao2020):
where D is the α-decay dose obtained from equation 1, M is the molecular weight (in mg) calculated using the EPMA data, and Nf and NA are the number of atoms per formula unit and Avogadro's number, respectively. Metamictisation of minerals represents the combined effects of radiation damage accumulation and annealing (Nasdala et al., Reference Nasdala, Wenzel, Vavra, Irmer, Wenzel and Kober2001a, Reference Nasdala, Beran, Libowitzky and Wolf2001b). Therefore, the annealing rate needs to be considered to correct the α-decay dose of columbite. Alpha-decay dose can be corrected for the effects of annealing using the average lifetime (200 Ma) of α-recoil trajectories of columbite (Lumpkin, Reference Lumpkin1992, Reference Lumpkin1998; Feng et al., Reference Feng, Liang, Linnen, Zhang, Zhou, Zhang and Gao2020). The correction is done as follows:
The degrees of metamictisation for samples JR-1, DKLS 107, DKLS 108, K802 and K803 were calculated using equations 3 and 4. Some of the U–Pb isotope ratios and trace-element concentrations required for the calculations of sample JR-1 are provided in Table 3. The full dataset is provided in Supplementary Table S1 and S4. The average degree of metamictisation calculated for each sample is provided in Table 4.
Corrected α-decay doses, displacements per atom, and key parameters of the strongest Raman bands of the analysed columbite grains.*
*The unit for aD is α-events/mg; adpa represents the amount of atomic displacement caused by each α decay; p – band position; S.D. – standard deviation.
Feng et al. (Reference Feng, Liang, Linnen, Zhang, Zhou, Zhang and Gao2020) demonstrated that the aD values of fully crystalline and fully metamict columbite are <0.2 × 1016 α-events/mg and >1.0 × 1016α-events/mg, respectively, with columbite in a transitional state having aD values between these. According to Table 4, columbite in samples JR-1, K802 and K803 have relatively high degrees of crystallinity and did not undergo metamictisation. Columbite in samples DKLS 107 and DKLS 108, however, are in a state transitional to fully crystalline and fully metamict.
Discussion
Factors influencing the position and FWHM of the strongest columbite Raman band
Composition
Moreira et al. (Reference Moreira, Rubinger, Krambrock and Dias2010) demonstrated that the strongest Raman band of columbite occurs at ~882 cm–1 and was related to the symmetric stretching vibration (A 1g) of the BO6 chemical bond. As there is a difference in bond length between Nb–O and Ta–O, the average bond length of B–O should be related to the amount of Nb and Ta in this site. Therefore, the contents of Nb and Ta in columbite probably influences the position of its strongest Raman band (Moreira et al., Reference Moreira, Rubinger, Krambrock and Dias2010). In contrast, variation in the amount of Fe and Mn has no significant impact on the bond length of the B–O chemical bond in columbite (Tarantino and Zema, Reference Tarantino and Zema2005), and should have no effect on the position and FWHM of its strongest Raman band (Moreira et al., Reference Moreira, Rubinger, Krambrock and Dias2010). There are no significant variations in the Fe and Mn contents of columbite grains in samples K802, K803 and K650, however there are significant variations in their Nb and Ta contents (Table 2). The Raman spectra of columbite in samples K802, K803 and K650, which were acquired in the same locations as the EPMA, demonstrate that the positions of the strongest Raman bands occur within a narrow range of wavenumbers of 870.1–878.3 cm–1 and FWHM values of 13.2–26.3 cm–1. Additionally, no correlation exists between band position, FWHM values, and Nb/Ta ratio of columbite (Table 4). Therefore, the Nb and Ta contents of columbite do not significantly influence the position and FWHM value of the strongest Raman band of columbite in the samples studied here.
Crystal orientation
Moreira et al. (Reference Moreira, Rubinger, Krambrock and Dias2010) determined the Raman spectra of the same columbite-(Mn) crystal at different orientations. Their results demonstrated that the strongest band (A 1g vibration mode) was always present at ~882 cm–1, but the FWHM value varied over a narrow range of 12–19 cm–1. Different A 1g, B 1g, B 2g and B 3g vibration modes were recognised in the Raman spectra, however, when a, b and c were parallel or perpendicular to the polarisation direction of the incident laser. This indicates that changing the orientation of columbite crystals has only a limited effect on the position and FWHM value of the strongest band. In this study, crystal orientations were obtained on columbite crystals in samples JR-1 and DKLS 108 using EBSD (see Fig. 1 and Fig. 3). Raman spectra were obtained when the c axes of the crystals were parallel and perpendicular to the vibration direction of the excitation laser; the results are illustrated in Fig. 6. For sample JR-1, when the c axis is perpendicular to the vibration direction of the excitation laser, the position of the strongest band is 872–876 cm–1 and the FWHM value is 17–24 cm–1. When the c axis is parallel to the vibration direction of the excitation laser, the position (868–869 cm–1) and FWHM (24–25 cm–1) of the strongest band are similar. For sample DKLS 108, when the c axis is perpendicular to the vibration direction of the excitation laser, the position of the strongest band is 862–864 cm–1 and the FWHM value is 31–34 cm–1. When the c axis is parallel to the vibration direction of the excitation laser, the position (861–866 cm–1) and FWHM (27–33 cm–1) of the strongest Raman band are, again, similar, indicating that crystal orientation has limited effect on the position and FWHM values of the strongest Raman band of columbite.
Raman spectra of columbite grains in samples: (a) JR-1; and (b) DKLS 108. Analyses JR-1.1 and JR-1.2 represent spectra of crystals where the c axis is perpendicular to the vibration direction of the excitation laser, whereas JR-1.3 and JR-1.4 represent spectra of crystals where the c axis is parallel to the vibration direction of the excitation laser. Analyses DKLS 108-1 and DKLS 108-2 represent spectra of crystals where the c axis is parallel to the vibration direction of the excitation laser, whereas DKLS 108-3 and DKLS 108-4 represent spectra of crystals where the c axis is perpendicular to the vibration direction of the excitation laser.
Radiation damage
Metamictisation of U- and Th-rich minerals is attributed to the destruction of the crystal lattice by radiation damage (mainly α-decay). When the cumulative radiation damage exceeds the self-healing rate of the mineral, metamictisation will occur (Ewing, Reference Ewing1975; Nasdala et al., Reference Nasdala, Wenzel, Vavra, Irmer, Wenzel and Kober2001a; Romer, Reference Romer2003). Previous studies have shown that the band position and FWHM of the (υ3(SiO4)) Raman band of zircon changes systematically with degree of metamictisation (Nasdala et al., Reference Nasdala, Irmer and Wolf1995; Zhang et al., Reference Zhang, Salje, Farnan, Graeme-Barber, Daniel, Ewing, Clark and Leroux2000; Balan et al., Reference Balan, Neuville, Trocellier, Fritsch, Muller and Calas2001; Nasdala et al., Reference Nasdala, Beran, Libowitzky and Wolf2001b; Nasdala et al., Reference Nasdala, Smith, Kaindl and Zieman2004; Nasdala et al., Reference Nasdala, Hanchar, Kronz and Whitehouse2005). Similarly, with increasing degree of metamictisation, the position of the strongest Raman band of columbite in each sample decreases gradually and the FWHM increases gradually (see Table 4). This indicates that radiation damage has a significant influence on the Raman characteristic band of columbite.
Murakami et al. (Reference Murakami, Chakoumakos, Ewing, Lumpkin and Weber1991) noted that the α-decay dose threshold for fully metamict zircon is 8.0 × 1015 α-events/mg. According to Feng et al. (Reference Feng, Liang, Linnen, Zhang, Zhou, Zhang and Gao2020), the α-decay dose threshold for fully metamict columbite is 1 × 1016 α-events/mg (after correcting for the effects of annealing), which is much larger than that of zircon. This indicates that columbite is probably more resistant to radiation damage than zircon in highly fractionated granites and granitic pegmatites. Therefore, the U–Pb age of columbite is more suitable than that of zircon to constrain the age of highly fractionated granites and granite pegmatites.
Correlation between position and FWHM of the strongest Raman band (A1g vibration mode)
The position and FWHM values of the Raman characteristic band of the analysed columbites are correlated negatively (Fig. 5). This negative correlation is evident in columbite-(Mn) grains collected from the same pegmatite, and also in those characterised by different degrees of metamictisation. From Table 4, we know that these samples have different degrees of metamictisation, however the positions of the Raman characteristic bands are all negatively correlated with FWHM values according to Fig. 5. The correlation between the two parameters can be quantified using the following equation (R 2 = 0.9431):
This correlation equation can be used to quantitatively characterise the relationship between the position and FWHM of the strongest Raman band of columbite. This negative correlation among variably metamict columbite is similar to that for zircon that has undergone variable degrees of metamictisation (Nasdala et al., Reference Nasdala, Hanchar, Kronz and Whitehouse2005). This supports the suggestion that the variation in position of the characteristic Raman band of columbite is not controlled by its composition or crystal orientations, but rather is primarily controlled by its degree of metamictisation.
Quantifying the degree of metamictisation of columbite using Raman spectroscopy
The results that are summarised in Table 4 demonstrate that, with increasing degree of metamictisation, the position of the strongest Raman band of columbite shifts towards lower wavenumbers and its FWHM value increases. A relationship can, therefore, be established to characterise quantitatively the degree of metamictisation of columbite using the position and FWHM values of the strongest Raman band (Fig. 7), in a manner similar to that proposed for zircon (Nasdala et al., Reference Nasdala, Hanchar, Kronz and Whitehouse2005). In contrast to the relationship for zircon, however, our results suggest that band position and FWHM are related to degree of metamictisation of columbite via a logarithmic relationship (Fig. 7). For columbite, the quantitative relationship between α-decay dose aD (× 1015 α-events/mg) and the FWHM (cm–1), displacements per atom (adpa) and the FWHM (cm–1), α-decay dose aD (× 1015 α-events/mg) and the band position p (cm–1), and displacements per atom (adpa) and the band position p (cm–1) can be quantified, respectively, via following relationships:
Bivariate diagrams illustrating the relationships between: (a) alpha-decay dose (aD) vs. FWHM; (b) alpha-decay dose (aD) vs. band position (p); (c) displacements per atom (adpa) vs. FWHM; and (d) displacements per atom (adpa) vs. band position (p) for the characteristic Raman band of columbite (A 1g vibration mode related to B–O bond) from all of the pegmatite samples.
On the basis of these relationships, once the FWHM value of the strongest band of columbite reaches ~50 cm–1 (Fig. 7), the crystal is likely to be completely metamict (the associated alpha-decay dose [aD] corrected for annealing is 1.0 × 1016 α-events/mg) (Lumpkin, Reference Lumpkin1998; Feng et al., Reference Feng, Liang, Linnen, Zhang, Zhou, Zhang and Gao2020) and not be suitable for geochronological studies.
Conclusions
On the basis of systematic Raman spectroscopic analysis of columbite grains from six samples (JR-1, DKLS 107, DKLS 108, K802, K803 and K650) from five pegmatites (the Jing'erquan No. 1 spodumene pegmatite, the Dakalasu No. 1 pegmatite dyke, the No. 802 and No. 803 spodumene pegmatites, and the No. 650 pegmatite dyke of Kalu'an), the following conclusions can be drawn regarding controls on the parameters of their strongest Raman bands: (1) the main factor that influences the position and FWHM value of the Raman characteristic band (A 1g vibration mode related to B–O bond) of columbite-(Mn) is metamictisation caused by radiation damage; composition and crystal orientations have limited influence; (2) the position of the strongest Raman band of columbite is correlated negatively with FWHM values; and (3) with increasing degree of metamictisation, the band position of the characteristic Raman band shifts towards lower wavenumbers and the FWHM value systematically increases.
The results presented here demonstrate that the position and FWHM value of the characteristic Raman band for columbite can be used to quantitatively assess their degree of metamictisation.
Acknowledgements
This research was supported financially by The National Natural Science Foundation of China (grant No. 41902073 and 91962214), Department of Science and Technology of Shaanxi Province, China (grant No. 2020JM-215), and the Fundamental Research Funds for the Central Universities (grant No. 300102272504). The authors would like to thank Dr. Huan Hu at Nanjing University for their assistance with EPMA and LA–ICP–MS analyses. We appreciate the constructive comments from the two anonymous reviewers that greatly helped to improve the manuscript.
Supplementary material
The supplementary material for this article can be found at https://doi.org/10.1180/mgm.2023.18.
Competing interests
The authors declare none.
Open access
