Gadolinite-( Nd ) , a new member of the gadolinite supergroup from Fe-REE deposits of Bastnäs-type , Sweden

A newmember of the gadolinite supergroup, gadolinite-(Nd), IMA2016-013, ideally Nd2FeBe2Si2O10, was found in the Malmkärra mine, ∼3.5 km WSW of Norberg, south-central Sweden, where it occurs in association with fluorbritholite-(Ce), västmanlandite-(Ce), dollaseite-(Ce), bastnäsite-(Ce) and tremolite. Gadolinite-(Nd) forms anhedral grains up to 150 μm in size, commonly occurring as aggregates of olive green colour. The crystals are transparent with vitreous to adamantine lustre. Gadolinite-(Nd) is brittle with conchoidal fracture, no cleavage or parting was observed. It has a white streak, the Mohs hardness is 6.5–7 and the calculated density is 4.86 g cm. Optically, the mineral is weakly pleochroic in shades of olive green, biaxial (–), nα = 1.78(1), nβ(calc.) = 1.80, nγ = 1.81(1) measured in white light, 2V(meas.) = 62(3)°. Electron-microprobe and laser ablation inductively coupled plasma mass spectrometry analysis [in wt.%] provided SiO2 21.77, Y2O3 5.49, La2O3 2.78, Ce2 O3 14.04, Pr2O3 3.28, Nd2O3 19.27, Sm2O3 5.30, Eu2O3 0.24, Gd2O3 4.10, Tb2O3 0.36, Dy2O3 1.32, Ho2O3 0.18, Er2O3 0.38,MgO 0.51, CaO 0.14,MnO 0.10, FeO 10.62, B2O3 0.10, BeOcalc. 8.99, H2Ocalc. 0.55 and total 99.52 giving the following empirical formula (based on 2 Si): (Nd0.632 Ce0.472Y0.268Sm0.168Gd0.125Pr0.110La0.094Dy0.039Ca0.014Er0.011Tb0.011Eu0.008 Ho0.008)∑1.957(Fe0.816Mg0.070Mn0.008)∑0.894(Be1.984B0.016)∑2.000Si2O9.660OH0.337. A weak Raman vibration band at ∼3525 cm confirms the presence of water in the structure. Gadolinite-(Nd) is monoclinic, P21/c, with a = 4.8216(3) Å, b = 7.6985(4) Å, c = 10.1362(6) Å, β = 90.234(4)°, V= 376.24(6) Å 3 and Z = 2. The strongest X-ray diffraction lines are [dobs in Å (hkl) Irel]: 4.830 (100) 72, 3.603 (021) 37, 3.191( 112) 52, 3.097 (013) 35, 2.888 (121) 100, 2.607(113) 49, 2.412 (200) 24. Along with the Malmkärra mine, gadolinite(Nd) was also recorded also at Johannagruvan and Nya Bastnäs. The minerals of the gadolinite subgroup together with fluorbritholite-(Ce) incorporate the highest fraction of medium-to-heavy rare-earth elements among associated rare-earth elementminerals in theMalmkärramine and possibly in all Bastnäs-type deposits.


Introduction
THE mineral gadolinite has a special place in the history of science, and specifically so in chemistry and mineralogy. It was the study of what was then called "ytterbite" from the Ytterby pegmatite north of Stockholm, Sweden, by the Finnish-Swedish chemist Johan Gadolin (1760Gadolin ( -1852, that led to the first discovery of a rare-earth element (REE) in the form of "yttria" (Gadolin, 1794). Gadolin's yttria was later proven to be a mixture of REEs, dominated by yttrium, while the "ytterbite" was finally named gadolinite (Ekeberg, 1802), now gadolinite-(Y). Presently, the REEs are very high on the European Union's lists of critical metals (http:// europa.eu/rapid/press-release_IP-14-599_en.htm) and neodymium is the most in demand. Systematic research towards the understanding of REE oreforming processes, a better knowledge of the distribution of these metals among the individual ore minerals, and the origin of these metals is of increasing interest.
The newly-defined gadolinite supergroup (Bačík et al., 2017) includes silicates, phosphates and arsenates with the general formula A 2 MQT 2 O 8 j 2 , where the individual structural sites are occupied as follows A: Ca, REE (Y + lanthanoids) > U, Th, Pb, Mn 2+ and Bi; M: Fe 2+ , □ (vacancy) > Mg, Mn 2+ , Zn, Cu, Al; Q: Be, B > Li; T: Si > P, As, B, Be, S; j: O, OH > F (Bačík et al., 2017). Based on the T-site occupancy, the gadolinite supergroup is divided into the gadolinite [Si 4+ > (P 5+ + As 5+ )] and herderite [(P 5+ + As 5+ ) > Si 4+ ] groups. Moreover, there are two subgroups within the gadolinite group, of which the gadolinite subgroup is dominated by REE 3+ cations, and the datolite subgroup is dominated by divalent cations, in the A site. Accordingly, the Q-site is occupied predominantly by Be 2+ in the gadolinite subgroup and by B 3+ in the datolite subgroup. The herderite group is also subdivided into the herderite and drugmanite subgroups (cf. Table 1).
The members of the gadolinite subgroup with a magmatic origin are commonly metamict due to significant amounts of actinides substituting for REE, which complicates their structural investigation. Well-crystalline samples are scarce and usually of a hydrothermal origin (Demartin et al., 1993).
A neodymium-dominant analogue of gadolinite-(Y) or gadolinite-(Ce) was suggested by Holtstam and Andersson (2007) from Bastnäs-type skarn deposits (Bergslagen ore province, south-central Sweden), but due to the lack of a sufficient amount of material, they did not provide any data apart from the chemical microanalyses. This paper provides a full characterization, including chemical composition, structural, physical and optical properties of the new mineral gadolinite-(Nd), which was approved recently by the Commission on New Minerals, Nomenclature and Classification (CNMNC) of the International Mineralogical Association (IMA), under the number 2016-013. The name is given in analogy with the other members of the gadolinite subgroup (Bačík et al., 2017), as the neodymiumdominant species. Type material is deposited under catalogue number B 11298 at the Moravian Museum, Brno, Czech Republic.

S134
Geological setting and the type locality Gadolinite-(Nd) was found in samples from the Malmkärra, Johanna and Nya Bastnäs mines, all belonging to the Fe-REE 'Bastnäs-type' skarn deposits in the Bergslagen ore province (Geijer, 1961;Holtstam and Andersson, 2007;Jonsson et al., 2014;Holtstam et al., 2014). The present mineral characterization was only undertaken on material from the Malmkärra mine (60°4′N, 15°51′ E), ∼3.5 km WSWof Norberg due to the very small size of gadolinite-(Nd) crystals (<30 µm) at the other localities. The Bastnäs-type deposits occur in the northwestern part of the Bergslagen ore province, south central Sweden, in the so-called 'REE-line' (Jonsson and Högdahl, 2013;Goodenough et al., 2016). This is an over 100 km long, north-easttrending, narrow belt within Svecofennian (1.91-1.88 Ga), mainly felsic metavolcanic rocks with discontinuous, interlayered carbonate (marble) horizons. The volcanic sequences were intensively metasomatized at the synvolcanic stage, and much of the carbonates have been replaced by skarn, hosting the localized Fe-REE mineralizations. The skarn-hosted Fe-oxide-REE-polymetallic mineralization at Bastnäs, in the central part of the REEline, has been dated from the Re-Os contents of molybdenite to give ages between 1.90 and 1.84 Ga. The mineralization is suggested to be the products of reactions between carbonate interlayers with Si, F, Cl, S, CO 2 and REE-bearing, hightemperature fluids of a magmatic origin (Holtstam et al., 2014;Sahlström et al., 2015). The term 'Bastnäs-type deposits' was introduced by Geijer (1961) to characterize a group of similar, Fe + REE ± Cu ± Mo ± Bi ± Au mineralized skarns that supposedly have a common origin to the Nya Bastnäs deposit. Based mainly on differences in mineral assemblages, Holtstam and Andersson (2007) divided them into two subtypes: deposits of subtype 1 are essentially only enriched in lightrare-earth elements (LREE), whereas those of subtype 2 are enriched in both light and heavyrare-earth elements (LREE and HREE, respectively).
The Malmkärra deposit (subtype 2) is hosted by a partly dolomitic marble that forms a conformable interlayer within the surrounding, typically extensively altered metavolcanic rocks. The skarn ore body, formed through reactions and replacement of parts of the carbonate interlayer, occurs along the contact with the more extensively altered metavolcanic rocks. Magnetite is the sole ore mineral that was mined, and tremolitic amphibole, humitegroup minerals, commonly serpentine-altered, biotite and phlogopite make up the main skarn silicates; sulfide minerals ( pyrite, chalcopyrite and molybdenite) are locally present as well (Geijer, 1936). Allanites sensu lato represent the most widespread REE minerals scattered in the skarn, as well as in the iron ore. During mining, a ∼0.5 m thick zone of massive REE silicates [dominated by fluorbritholite-(Ce), västmanlandite-(Ce) and dollaseite-(Ce)] was encountered in the deeper levels of the mine, along the contact between the marble and the skarn iron ore (Geijer, 1936;Andersson, 2004). The REE minerals reported from the Malmkärra mine are summarized in Table 2.  (2007), Holtstam et al. (2005).
The Malmkärra iron mine has a long history. The oldest preserved record is from 1664, with major mining operations taking place from 1885 up to its final closure in 1936 (Geijer and Magnusson, 1944;Andersson, 2004). Malmkärra is also the type locality of västmanlandite-(Ce), see Holtstam et al. (2005).

Mineral association, appearance and physical properties
The skarn rock sample (∼10 cm × 10 cm × 5 cm) containing gadolinite-(Nd) was found on the dumps of the Malmkärra mine. Gadolinite-(Nd) forms anhedral grains up to 150 µm in size, occurring as aggregates associated with tremolite, dollaseite-(Ce), västmanlandite-(Ce), fluorbritholite-(Ce) and bastnäsite-(Ce) (Fig. 1). Gadolinite-(Nd) is transparent and, under the optical microscope in thin section, it is pale olive green, while larger fragments are olive green. It has a white streak and a vitreous to adamantine lustre. No visible fluorescence under shortwave or longwave ultraviolet light was observed. The Mohs hardness is between 6.5 and 7. Gadolinite-(Nd) is brittle with a conchoidal fracture; no cleavage or parting was observed. The density (D calc ) is 4.86 g cm -3 , and was calculated using the empirical formula and unit-cell parameters obtained from the single-crystal X-ray diffraction (XRD) data. A density measurement could not be performed directly due to the lack of suitable material. Optically, gadolinite-(Nd) is weakly pleochroic in shades of olive green, biaxial (-), n α = 1.78(1), n β(calc.) = 1.80, n γ = 1.81(1) measured in white light. The poor optical quality of the measured fragments did not allow more precise determination of the refractive indices. The angle of the optical axes, 2V (meas.) : 62 (3)°was FIG. 1. Back-scattered electron images of gadolinite-(Nd) and associated REE minerals from (a-c) Malmkärra and (d ) Johannagruvan. gd-gadolinite-(Nd), ds-dollaseite-(Ce), fb-fluorbritholite-(Ce), tr-tremolite, vl-västmanlandite-(Ce), bsbastnäsite-(Ce). The aggregate from which the crystals for single-crystal X-ray diffraction and optical study were extracted is indicated by an ellipse. The zoning of the gadolinite-(Nd) from Johannagruvan (d ) reflects Y-lanthanoids variability; the Y-rich zones are darker. S136 measured using extinction by means of a spindle stage and computed by the EXCALIBR II software. Dispersion is strong (r < v).

Composition
The compositions of gadolinite-(Nd) and associated minerals were determined by means of electron-probe microanalysis (EPMA) with a CAMECA SX100 in wavelength-dispersive mode at an accelerating voltage of 15 kV, a beam current 20 nA and a beam diameter of 5 µm. The standards, X-ray lines and monochromators used for analysis are listed in Table 2. Aluminium, Pb, Tm, Yb and Lu were also analysed for, but their contents were below the detection limit. The peak counting times varied from 10 s for the main elements to 60 s for the minor elements and high and low energy background was counted for ½ of the peak counting time of the relevant analytical line. The raw data were processed using the X-phi matrix correction routine (Merlet, 1994). Based on the counting statistics, the measurement error expressed as 2σ is approximately <1 rel.% for concentrations of ∼20 wt.% and ∼8 rel.% for concentrations of ∼1=wt.%. The contents of Be, B and other trace elements below the detection limit of the EPMAwere sought for by means of an Agilent 7500ce quadrupole Inductively Coupled Plasma Mass Spectrometer (ICP-MS; Agilent 7500ce, Santa Clara, CA, USA) with an attached UP 213 laser ablation (LA) system (New Wave Research, Inc., Fremont, CA, USA). The LA-ICP-MS system consists of a nanosecond laser Nd:YAG operating at 213 nm ( pulse duration of 4.2 ns) and a SuperCell ablation chamber. Ablated material was transported from the sample chamber using helium carrier gas (1 L min -1 ) and mixed with argon (0.6 L min -1 ) prior to the torch. Potential polyatomic interferences were minimized by a collision reaction cell in He mode (1 mL min -1 ). The contents of major to trace elements were determined after laser ablation of individual spots at the following conditions: diameter 55 µm, strength of laser beam 4 J cm -2 , frequency 10 Hz and spot Si for internal calibration and the calculations were performed via peak area of isotopes.
The empirical formula has been calculated on the basis of 2 atoms per formula unit of Si and the amount of OH has been charge balanced to keep the empirical formula electroneutral.
The composition of gadolinite-(Nd) was examined from three thin sections prepared from a single piece of the massive REE-silicate-rich sample. The gadolinite shows quite variable chemical composition, particularly in the REE pattern and the Y/ (REE) ratio, but all analyses (n = 11) are Nddominant and correspond to gadolinite-(Nd), see Table 3. The calcium content is low (0.14-0.24 wt.% CaO). Holtstam and Andersson (2007) performed Mössbauer spectroscopy on the gadolinite-(Ce/Y) from Malmkärra, which revealed the presence of exclusively divalent Fe. Therefore, all Fe is assumed as divalent in gadolinite-(Nd). Apart from Fe, the M site is occupied by Mg (0.51-0.66 wt.% MgO) and Mn (0.10-0.14 wt.% MnO). All these cations do not fill the M site completely, thus a minor vacancy (0.093-0.143 □ pfu) is present. Seven grains, previously analysed by EPMA, were subsequently analysed by LA-ICP-MS to determine the Be and B content as well as the whole suite of the REEs and some trace elements (Table 4). A chondritenormalized REE pattern (using a combination of EPMA and LA-ICP-MS data) shows the maxima at Nd and depletion of the lightest and heavy REE as well as a moderate Eu anomaly (Eu N */Eu N = 0.16), see Fig. 2. The content of BeO obtained from the LA-ICP-MS analysis ranges from 8.69 to 11.02 wt.% with a mean of 10.11 wt.%. The theoretical BeO content, calculated from stoichiometry, is ∼9.0 wt.%. This discrepancy, ∼10 rel.%, is still within the analytical error of LA-ICP-MS for Be (∼10 rel.%). The average content of B 2 O 3 is 0.10 wt.%. Laser ablation ICP-MS also revealed very low U and Th contents in the gadolinite (∼80 and ∼5 ppm, respectively). Such a low radioactive element content does not affect the crystallinity of the mineral. A small amount of OH (0.34 pfu) resulted from the calculation of an empirical formula. The empirical formula obtained from the mean of three analyses of the crystal used for structural determination is: A part of the largest grain (150 µm × 100 µm) of gadolinite-(Nd), carefully measured by EPMA, was extracted subsequently from the thin section for the single-crystal X-ray diffraction study. The other parts of this grain were used for determination of the optical properties.

Raman spectroscopy
Raman analysis of gadolinite-(Nd), performed on uncoated polished sections, was undertaken by means of a Horiba LabRAM HR Evolution spectrometer. This dispersive, edge-filter-based system was equipped with an Olympus BX 41 optical microscope, a diffraction grating with 600 grooves per millimetre, and a Peltier-cooled, Sibased charge-coupled device (CCD) detector. Raman spectra of REE-bearing minerals are commonly obscured by laser-induced emissions of the REE (e.g. Nasdala et al., 2012;Lenz et al., 2015). After careful tests with different lasers (473, 523 and 633 nm), the 633 nm He-Ne laser (10 mW S138 at the sample surface) was selected for spectral acquisition to minimize analytical artefacts. With the Olympus 100× objective (numerical aperture 0.9) and the system being operated in the confocal mode, the lateral resolution was ∼1 μm. Wavenumber calibration was done using the Rayleigh line and Ne-lamp emissions. The wavenumber accuracy was >0.5 cm -1 and the spectral resolution was ∼2 cm -1 . Band fitting was done after appropriate background correction, assuming combined Lorentzian-Gaussian band shapes using PeakFit (Jandel Scientific Software). The Raman spectrum of gadolinite-(Nd) is shown in Fig. 3. Following the spectroscopic  S140 study of isostructural datolite and herderite (Frost et al., 2013(Frost et al., , 2014Goryainov et al., 2015)  The local maximum at 3370 cm -1 is caused by luminescence as proven by excitation using a 532 nm laser.

X-ray crystallography and structure determination
Powder X-ray diffraction data for gadolinite-(Nd), carefully extracted from the polished section by a diamond microgrinder and powdered under isopropyl alcohol in an agate mortar, were obtained using an Empyrean (PANalytical) powder diffractometer with the solid-state PIXcel 3D detector operating in Bragg-Brentano geometry. Data were collected over the range of 5-80°2θ with a step size of 0.013°and counting time of 2 s per step. Positions and intensities of diffraction maxima were obtained from a profile fitting performed with High-Score3 software (PANalytical). Unit-cell parameters were refined using least-squares based on the positions of 40 reflections utilizing Celref software (Laugier and Bochu, 2004). The experimental powder pattern was indexed in accordance with the calculated values of intensities obtained from the single-crystal structure refinement. Powder X-ray diffraction data are given in Table 5. Unit-cell parameters refined from the powder data are a = 4.826(1), b = 7.717(2), c = 10.147(2) Å, β = 90.16°and V = 377.9(2) Å 3 . A platy crystal fragment (0.07 mm × 0.06 mm × 0.01 mm) was examined on an Oxford Diffraction Gemini single-crystal diffractometer with Atlas S2 CCD detector, using monochromatized MoKα 1.14, 1.05 Weighting scheme, weights σ, 1/(σ 2 (F) + 0.0004F 2 ) Δρ min , Δρ max (e Å -3 ) -1.21, 1.56 S141 radiation (λ = 0.71073 Å) from a conventional sealed X-ray tube and collimated with fibre optics and Mo-Enhance collimator. The unit cell was refined from 1323 reflections using least-squares with the Crysalis software (2014, Agilent Technologies, Oxford, UK); ω rotational scans (frame width 1°, counting time 170 s per frame) were adopted to collect a set of three-dimensional intensity data. From a total of 2955 reflections, 875 were independent and 607 classified as observed, with [I obs > 3σ(I)]. Corrections for background, Lorentz effects and polarization were applied; Gaussian correction for absorption (shape + empirical scaling) was done in JANA2006 (Petríček et al., 2014) giving a data set with R int of ∼10%.
Reduction of the data was performed using the Crysalis package. The crystal structure of gadolinite-(Nd) was refined from single-crystal data using the structure model for gadolinite-(Y) published in Demartin et al. (1993) by full-matrix least-squares using JANA2006 (Petríček et al., 2014) based on F 2 . All non-oxygen atoms were refined with anisotropic displacement parameters. Refinement for 57 parameters and 18 constraints converged to R = 0.0371, wR = 0.0780 for 607 observed reflections (Goof = 1.14). Details of data collection, crystallographic data and refinement details are given in Table 6. Atom coordinates and displacement parameters are listed in Table 7; selected interatomic distances are provided in Table 8. The structure refinement of gadolinite-(Nd) provides no significant departure from the published structure results of the isostructural gadolinite-(Y) (Miyawaki et al., 1984;Demartin et al., 1993;Cámara et al., 2008). Minerals of the gadolinite group possess a layer structure composed of two distinct sheets, which are alternating in the c direction (Bačík et al., 2014). One sheet is represented by TO 4 and QO 4 tetrahedra and it alternates with a second sheet of AO 8 polyhedra and MO 6 octahedra. In the structure of gadolinite-(Nd), the T-site is occupied dominantly by Si 4+ , the Q-site by Be, the A-site by Nd 3+ (and other REE) and the M-site by Fe 2+ . During the refinement we treated A-sites as only occupied by Nd. We did not model substitutions of a large number of REE that have very similar scattering curves for X-rays as we are convinced that such a fit would not make much sense. The refined site occupancies for Nd and Be turned out they be slightly lower than unity (due to the effects of mixing with lighter elements at particular sites). The Mg substitution for Fe, as documented from the wavelength-dispersive spectroscopy analyses, was not successfully confirmed by the structure refinement. The structural formula for gadolinite-   (7) 0.003 (6) 0.013 (7) 0.003 (5) 0.003(5) −0.001(5) 1.96(4)/ 3.79 (7) O1 0.6835(14) 0.3964(9) 0.3597 (7)

S142
(Nd) is therefore Nd 1.942 FeBe 1.931 Si 2.069 O 10 , Z = 2. This formula is not electro-neutral due to the fact that the refined occupancies are only a proxy to the real situation, the refinement of all (excluding REE) substituting components is beyond the quality of current dataset. However, the refinement of the data suggested that a small amount of Si can enter the Q-sites and vice-versa.

Concluding remarks
Rare-earth element ores from the Bastnäs-type deposits are in general characterized by LREE enriched and HREE depleted patterns (e.g. Holtstam et al., 2014). Taking into consideration the low activity of P and elevated activity of CO 2 , the primary REE minerals are represented as silicates and carbonates. Among the REE phases of the Bastnäs-type mineralization at the Malmkärra mine, the gadolinite-group minerals together with fluorbritholite represent the phases which accommodate the highest amount of middlerare-earth elements (MREE) and HREE into their structure. On the other hand, epidote-supergroup minerals, including västmanlandite, are strongly depleted in MREE and HREE (Fig. 2).
Notably, the new mineral may represent a locally important sink for the critical metal neodymium in deposits of the Bastnäs-type. Geijer (1936) reported törnebohmite as occurring locally in "not insignificant amounts" at Malmkärra, while Holtstam and Andersson (2007) make no mention of having observed this mineral. Considering their, in part similar, optical properties, it seems possible that Geijer's törnebohmite could represent (a misidentified) gadolinite. Additionally, new observations also suggest that it may be more abundant than hitherto known in other Bastnäs-type mineralizations, specifically in the Johanna and Nya Bastnäs mines.