The crucial importance of tungsten in our modern technologies explains why this element is included in the list of critical raw materials published by the European Union (European Union, 2020). In most geological contexts, only two mineral substances are carrying this rare element: scheelite, CaWO₄, and “wolframite”, which designate members of a solid solution between the two species hubnerite, MnWO₄, and ferberite, FeWO₄.

In the Central African tungsten province, which extends from southwestern Uganda to central Rwanda and eastern Congo, most tungsten-bearing minerals of economic value are constituted by varieties of ferberite (Sahama, 1981). This species generally occurs as black platy crystals, showing a characteristic perfect cleavage plane and crystallizing as a primary phase in mineralized quartz veins. A remarkable variety of ferberite, called “reinite”, forms large pseudo-octahedral crystals without any cleavage plane and shows a dull black colour and an earthy lustre. This unusual variety mainly occurs in the “tungsten belt” of central Rwanda (Goldmann et al., 2013) and is produced by a pseudomorphic replacement of scheelite by ferberite.

In such “reinite” pseudomorphs, numerous late alteration phases occur, mainly consisting of tungsten and iron oxides (Sahama, 1981; Goldmann et al., 2013). The centre of the pseudomorphs is occupied by white earthy masses of anthoinite [Al(WO3)(OH)₃], the latest alteration product of the sequence, but several tungstate species are formed between anthoinite and the rim of ferberite. These species include tungstite [WO3⋅ H₂O], hydrotungstite and meymacite [polymorphs of WO3⋅ 2H₂O], hydrokenoelsmoreite [□2W₂O₆(H₂O)] (previously “ferritungstite” and “alumotungstite”), raspite [PbWO₄], mpororoite [WAlO₃(OH)3⋅ 2H₂O], yttrotungstite-(Ce) [CeW₂O₆(OH)₃] (previously “cerotungstite”), and yttrotungstite-(Y) [YW₂O₆(OH)₃] (previously “yttrotungstite”).

During a field campaign in Rwanda, the Luxembourgish geologist Paul Antun collected, in 1977, several samples from the famous Nyakabingo tungsten mine, close to Kigali. Those samples were donated to the Natural History Museum of Luxembourg, and its curator, Simon Philippo, discovered the yttrotungstite-(Nd) samples that form yellow-orange bladed crystals associated with anthoinite and hydrokenoelsmoreite.

The mineral, which is the Nd-equivalent of yttrotungstite-(Ce) and yttrotungstite-(Y), has been named in agreement with the nomenclature of the yttrotungstite group. It has been approved by the Commission on New Minerals, Nomenclature, and Classification of the International Mineralogical Association (IMA-CNMNC) under number IMA 2023-064 (Bosi et al., 2023); its abbreviation is “Ytgs-Nd” (Warr, 2021). Holotype samples used for the photographs, for the EMP analyses, and for the X-ray powder diffraction measurements are deposited in the collections of the Natural History Museum of Luxembourg, catalogue number WPA504. Part of the holotype used for the single-crystal X-ray diffraction measurements and for the optical measurements is deposited in the collections of the Laboratory of Mineralogy, University of Liège, Belgium (catalogue number ULG 21981).

The sample used for this study was collected in the Nyakabingo tungsten mine, Kigali, Rwanda (latitude 1°52^′33^′′ S, longitude 30°02^′36^′′ E). The region is well-known as an important source of tungsten and was historically called the “Tungsten belt”. The Nyakabingo deposit, located at 10 km NNW of Kigali, is part of the Mezoproterozoic Karagwe-Ankole Belt of Central Africa, a metallogenic province characterized by W- and Sn-bearing quartz veins, as well as by Nb–Ta–Sn-bearing granitic pegmatites (Hulsbosch et al., 2016).

The Nyakabingo deposit is hosted in a sequence showing an alternation of pyritiferous black shales, quartzites, and psammites that were affected by a low-grade greenschist metamorphism (De Clercq et al., 2008; Dewaele et al., 2016). The W-bearing quartz veins were dated by the Ar–Ar method on muscovite, giving an age of 992.4 ± 1.5 Ma, corresponding to the age of the G4 granites that are certainly the source of the W-bearing vein systems (Dewaele et al., 2016). Two different stages of W mineralization were observed, the first one producing scheelite and massive ferberite and the second one producing the ferberite pseudomorphs after scheelite (“reinite”) (Dewaele et al., 2016).

The quartz veins at Nyakabingo contain the major phases of muscovite, feldspar, pyrite, arsenopyrite, and xenotime-group minerals. The economically valuable minerals are ferberite, cassiterite, and molybdenite. During the last meteoric alteration stage, ferberite, feldspar, and pyrite were altered with the formation of clay minerals, iron oxides, and tungsten secondary phases such as tungstite and anthoinite (Dewaele et al., 2016).

Yttrotungstite-(Nd) was found as a secondary mineral occurring in the oxidized zones of the Nyakabingo tungsten deposit. Dewaele et al. (2016) describe, for the type locality, the following minerals in these oxidation zones: ferberite, anthoinite, hydrokenoelsmoreite, meymacite, raspite, yttrotungstite-(Ce), and iron oxides. In our samples that show a honeycomb texture, we identified yttrotungstite-(Nd) as yellow-orange platy crystals, associated with anthoinite and with microcrystalline yellowish hydrokenoelsmoreite octahedra (Fig. 1).

Yttrorungstite-(Nd) forms tiny yellow-orange tabular crystals, in association with anthoinite and hydrokenoelsmoreite (Fig. 1A). The crystals are transparent and reach a length of maximum of 100 µm; they frequently form radiated groups (Fig. 1B). The dominant form is {010}, and the β angle of ca. 105°, between a and c, is clearly visible on the platy crystals (Fig. 1A). Twinning was observed optically, parallel to (001), in comparison with yttrotungstite-(Ce) (Sahama et al., 1970).

The streak is white, the lustre is vitreous, and the mineral is non-fluorescent. The Mohs hardness is 1 by analogy with the yttrotungstite-(Y) and yttrotungstite-(Ce). Observed cleavage planes are (010) perfect and (100) and (001) excellent; the tenacity is brittle. The density has not been directly determined due to the small size of the crystals and the small amount of available material, but the calculated density is 6.25 g cm⁻³, obtained from the EMPA data and the single-crystal unit-cell parameters.

Optically, the mineral is biaxial, but the optical sign was impossible to measure due to the orientation of the lamellae parallel to the optical plane. Refraction indices, determined under sodium light (λ= 589 nm), are α= 1.90(2) and γ= 2.10(2). The optical plane is (010), parallel to the lamellae and to the perfect cleavage plane. X ∧ c = 6°, Z ∧ a = 9°. Pleochroism is from light yellow (X) to yellow (Z).

The chemical composition of yttrotungstite-(Nd) (Table 1) was determined with an electron microprobe at the Department of Earth Sciences “Ardito Desio” in Milano, Italy. The instrument was a Jeol JXA-8200, WDS electron microprobe, using an acceleration voltage of 15 kV, a beam current of 5 nA, and a beam diameter of 3 µm. A total of 17 point analyses were realized, but the water content was not determined due to the small amount of available material. The absence of fluorine in the mineral was confirmed by energy-dispersive X-ray spectrometry (EDS) measurements.

The empirical formula of yttrotungstite-(Nd) was calculated on the basis of 15 positive charges, assuming 1 (OH) group and 1 H₂O molecule per formula unit, according to the structural data (see below). This formula corresponds to the following: (Nd_0.36Sm_0.13Ce_0.12Pr_0.08La_0.06Gd_0.06Y_0.05Dy_0.03Yb_0.02 Ca_0.01)_Σ0.92(W_2.02P_0.02)_Σ2.04O₇(OH)(H₂O). The simplified formula is (Nd,Sm,Ce,Pr,La,Gd,Y)W₂O₇(OH)(H₂O), and the ideal formula is NdW₂O₇(OH)(H₂O), which requires WO₃ 72.99, Nd₂O₃ 22.72, and H₂O 4.29, for a total of 100.00 wt %.

X-ray powder diffraction data were collected at the University of Liège with a Rigaku Xcalibur four-circle diffractometer used in Debye–Scherrer geometry and equipped with a MoKα radiation source (λ= 0.71073 Å) and with an EOS charge-coupled device (CCD) detector (Table 2). Unit-cell parameters, refined from the powder data with the LCLSQ program (Burnham, 1991), are as follows: a= 5.84(1), b= 8.64(1), c= 7.045(8) Å, β= 106.1(1)°, V= 341.3(6) ų, Z = 2, space group P21/m.

Single-crystal X-ray diffraction studies were carried out with the same instrument on a crystal of yttrotungstite-(Nd) measuring 0.057×0.088×0.098 mm. A total of 350 frames with a spatial resolution of 1° were collected by the φ/ω scan technique, with a counting time of 20 s per frame, in the range of 6.00°<2θ<57.70°. A total of 870 reflections were extracted from these frames, corresponding to 808 unique reflections. The unit-cell parameters refined from these reflections, a= 5.8534(3), b= 8.6881(3), c= 7.0488(4) Å, β= 105.336(5)°, and V= 345.70(3) ų, are in very good agreement with those refined from the X-ray powder diffraction data. Data were corrected for Lorentz, polarization, and absorption effects, the latter with an empirical method using the SCALE3 ABSPACK scaling algorithm included in the CrysAlisRED package (Oxford Diffraction, 2007). Further details on the structure solution and refinement are given in Table 3.

The crystal structure of yttrotungstite-(Nd) (Figs. 2 and 3) was refined in space group P21/m. Scattering curves for neutral atoms, together with anomalous dispersion corrections, were taken from the International Tables for X-ray Crystallography, Vol. C (Wilson, 1992). In the final refinement cycle, all atoms were refined anisotropically, leading to the final R1 value 0.0362. Atomic coordinates and site occupancies, as well as equivalent isotropic displacement parameters, are given in Table 4.

The crystal structure is characterized by edge-sharing WO₆ octahedra, forming kinked chains parallel to the b axis (Fig. 2). The NdO₈ complex polyhedra connect three chains together by sharing edges or corners with adjacent WO₆ octahedra. A view perpendicular to the b axis shows that the NdO₈ and WO₆ polyhedra occur in two different types of planes, stacked parallel to the [-101] direction (Fig. 3).

Neodymium occurs in a very large eight-coordinated site, which can be described as two tetragonal pyramids sharing one common basal edge and sharing the opposite basal edge with a square (Fig. 4). The average Nd–O bond length is 2.446 Å, and the average W–O bond length is 1.955 Å (Table 5). The WO₆ octahedron is strongly distorted, with three short bonds (1.764–1.881 Å) and three long bonds (1.961–2.204 Å; Table 5).

Bond-valence sums were calculated for yttrotungstite-(Nd) according to the empirical parameters of Gagné and Hawthorne (2015) (Table 6). Those sums are in fairly good agreement with the ideal values expected for W (5.94) and Nd (3.31), as well as with those expected for the four oxygen atoms O1 to O4 (1.91–2.06). OH5 shows a bond-valence sum very close to the ideal value for an OH group (0.84), and OW6 shows a bond-valence sum of 0.33, indicating that this site is occupied by a water molecule. The end-member formula of yttrotungstite-(Nd) consequently contains one OH group and one water molecule per formula unit and can be written as NdW₂O₇(OH)(H₂O).

The Raman spectrum of yttrotungstite-(Nd) (Fig. 5) was recorded at the University of Luxembourg on a Renishaw inVia micro-Raman spectrometer. The microscope connected to the spectrometer is a Leica equipped with a 50× objective; the working distance was 8 mm. Four wavelengths were used at 785, 633, 532, and 442 nm. The spectra given in Fig. 5 were collected at 633 nm, which offers the best spectral resolution, on four different mineral grains. The laser power was kept at ∼ 0.5 mW to avoid any heating effect. The data were collected in three cycles of 30 s exposure with a 1200 g mm⁻¹ grating, in a stepwise manner in order to cover the full range up to 4000 cm⁻¹. The Renishaw software Wire 5.1 was used for spectra acquisition. The spectra displayed in Fig. 5 are raw data without any processing, particularly without any baseline correction.

The spectrum is characterized by stretching vibrations of the W–O bonds, occurring between ca. 600 and 1000 cm⁻¹, as well as by stretching vibrations of the OH groups around 3300–3600 cm⁻¹. The absence of peaks at ca. 1600 cm⁻¹, which would correspond to water bending vibrations, is certainly due to polarization effects. By using the empirical correlation between the wavenumbers of the W–O stretching vibrations and the W–O bond lengths established by Hardcastle and Wachs (1995), it was possible to assign the intense peak at 930 cm⁻¹ to the vibrations of the W1–O1 and W1–O3 bonds (calculated wavelengths 900 and 901 cm⁻¹, respectively) and the peaks at 834 and 684 cm⁻¹ to the vibrations of the W1–O2^′ and W1–O4 bonds (calculated wavelengths 722 and 620 cm⁻¹, respectively). The two long W1–OH5 and W1–O2 bonds (Table 5) produce vibrations around 400 cm⁻¹, but they cannot be accurately attributed to any peak since such long W–O bonds have lost most of their covalency and consequently give rise to weak peaks difficult to assign (Hardcastle and Wachs, 1995).

In the OH-stretching region, five peaks occur that are produced by the vibrations of the OH5 group and of the OW6 water molecule. By using the empirical correlation between the OH-stretching wavenumbers and the O–H…O distances, established by Libowitsky (1999), it is possible to accurately calculate the theoretical wavenumbers corresponding to the O–H…O distances measured in the crystal structure model (Table 5). A comparison of these calculated values with the wavenumbers observed on the Raman spectrum (Fig. 5) indicates that the vibrations of the OH5 groups correspond to the bands at 3350, 3384, and 3434 cm⁻¹, while the vibrations of the OW6 water molecule produce the bands at 3522 and 3600 cm⁻¹ (Table 5).

The new mineral species described in the present paper, yttrotungstite-(Nd), belongs to the yttrotungstite group and is isostructural with yttrotungstite-(Ce) (Sahama et al., 1970) and with yttrotungstite-(Y) (Davis and Smith, 1971). These two species, originally described as “cerotungstite” and “yttrotungstite”, were renamed by Nickel and Mandarino (1987). Yttrotungstite-(Nd) can be numbered in the Strunz classification as 04.FD.20 (Strunz and Nickel, 2002), and its Dana number is 49.03.03.04 (Gaines et al., 1997). It corresponds to the unnamed mineral UM2004-//-O:CeHNdSmW, described in the Clara Mine, Wolfach, Germany, by Kolitsch et al. (2004). A comparison of the physical properties of minerals belonging to the yttrotungstite group is given in Table 7. The compatibility index of yttrotungstite-(Nd), according to Mandarino (1981), is 1-(KP/KC) = 0.0583, good (calculated from the density of 6.25 g cm⁻³ obtained from single-crystal unit-cell parameters and electron microprobe analyses).

According to the structural data, the end-member formula of yttrotungstite-(Nd) is written as NdW₂O₇(OH)(H₂O) since the mineral contains one OH group and one water molecule per formula unit. In the yttrotungstite group, however, the formulae of yttrotungstite-(Ce) and yttrotungstite-(Y) were still expressed as CeW₂O₆(OH)₃ and Y(W,Fe,Si,Al,Ti)₂(O,OH,H₂O)₉ in the IMA-CNMNC list, respectively. Structural data are available for yttrotungstite-(Y) (Davis and Smith, 1971) and by using the atomic coordinates of this species, bond-valence calculations indicate bond-valence sums of 1.13 on the O1 site (which corresponds to our OH5 site) and of 0.31 on the Wat6 site (which corresponds to our OW6 site). These data therefore confirm the presence of one water molecule and one OH group in the structure of yttrotungstite-(Y), thus allowing its formula change to YW₂O₇(OH)(H₂O). The same procedure is applied to the formula of yttrotungstite-(Ce), which becomes CeW₂O₇(OH)(H₂O), in agreement with the uniformization of the yttrotungstite group accepted by the CNMNC (Bosi et al., 2026).