Introduction
Pure bismuth tellurides are a group of five minerals (ehrigite, hedleyite, pilsenite, tsumoite and tellurobismuthite) with varying characteristics, including their compositions and crystallinity (Anthony et al., Reference Anthony, Bideaux, Bladh and Nichols1990). One of the ways to differentiate them from each other is the Bi/Te ratio of their structural formulas (e.g. Cook et al., Reference Cook, Ciobanu, Wagner and Stanley2007; Gu et al., Reference Gu, Watanabe, Hoshino and Shibata2001), which varies from 0.667 for tellurobismuthite (Bi2Te3) to 2.667 for the newly discovered ehrigite (Bi8Te3; Ciobanu et al., Reference Ciobanu, Cook, Yao, Slattery, Wade and Dima2024). Most bismuth tellurides belong to the trigonal/hexagonal crystal system and have a layered structure (Anthony et al., Reference Anthony, Bideaux, Bladh and Nichols1990). Their diversity is controlled by different factors, specifically the order of layer alternation, the tolerance of their crystalline structure to incorporating substitution elements, and the possibility of an ordered distribution of substitution elements in their structure (Belogub et al., Reference Belogub, Moloshag, Novoselov and Kotlyarov2011; Cook et al., Reference Cook, Ciobanu, Wagner and Stanley2007).
In terms of the composition of these minerals, Bi is commonly substituted by Pb, whereas Te is substituted by S and less commonly by Se, leading to the existence of several minerals, many of which are rather rare, with varying complexity, including a few solid solutions (e.g. Anthony et al., Reference Anthony, Bideaux, Bladh and Nichols1990; Cook et al., Reference Cook, Ciobanu, Wagner and Stanley2007; Gu et al., Reference Gu, Watanabe, Hoshino and Shibata2001; Alexandre et al., Reference Alexandre, Heine, Fayek, Potter and Sharpe2019). However, when only Pb and Se substitutions in bismuth tellurides are considered (Fig. 1), it becomes apparent that it is much more common for one or the other of these substitutions to exist to the point of defining a distinct mineral, but it is rather uncommon for both Pb and Se substitutions to define a specific discrete mineral, with poubaite (PbBi2Te2Se2) being the unique example (Fig. 1). Other cations, such as Fe, Cu, Sb, Sn and Ag, can also substitute for Bi in pure bismuth tellurides to varying degrees, from very low to sufficiently high to form distinct minerals (Anthony et al., Reference Anthony, Bideaux, Bladh and Nichols1990). In general terms, the compositions of bismuth tellurides and its substitutions are relatively well studied, even though some knowledge gaps remain; the most notable of these involve the compositional variability of these minerals and the extent of the cation and anion substitutions.

Figure 1. Theoretical positions of minerals in Bi–Pb–Te–Se space. The pure bismuth tellurides, tsumoite, tellurobismuthite and the Pb-bearing kochkarite and rucklidgeite were subject of a literature data compilation.
As a part of a large-scale project on the comparative ore mineralogy of several volcanogenic massive sulfide deposits in the Flin Flon – Snow Lake Mineral Belt (Manitoba, Canada; Fig. 2), a bismuth telluride mineral, or rather a range of compositions, with significant and highly variable degrees of both Pb and Se substitutions was observed at the Stall Lake VMS deposit. It was soon appreciated that the degree of Pb substitution was comparable to that of other Pb-rich bismuth tellurides such as rucklidgeite ((Bi,Pb)3Te4) and kochkarite (PbBi4Te7). However, the degree of Se substitution in this phase was significantly higher than that of any of these minerals, without reaching the levels of kawazulite (Bi2Te2Se), skippenite (Bi2TeSe2), or telluronevskite (Bi3TeSe2).

Figure 2. Simplified geological context of the Stall Lake and other VMS deposits within North America (a), the Trans-Hudson Orogen of Canada (b), and the Snow Lake area (c). The approximate extent of the Flin Flon – Snow Lake Mineral Belt (FF-SL MB) is indicated in (b). SK=Saskatchewan, MB=Manitoba.
It was the unusually high degrees of Pb and Se substitutions in this phase that prompted a comprehensive comparison with the all literature data available for similar minerals, specifically the almost pure bismuth tellurides tsumoite and tellurobismuthite, and the Pb-rich bismuth tellurides rucklidgeite and kochkarite. Tsumoite and tellurobismuthite were preferred for this comparison rather than other pure bismuth tellurides, as their Bi/Te ratios are very similar to those of the Stall Lake phase. Rucklidgeite and kochkarite were selected because their Pb contents are very similar to that of the Stall Lake phase. The conclusions of this comparison are presented in this paper, together with observations regarding Pb and Se substitutions in bismuth tellurides.
Geological background of the Stall Lake deposit
The mined-out Stall Lake deposit is one of the 12 volcanogenic massive sulfide (VMS) deposits from the Snow Lake arc assemblage in Manitoba, Canada, itself the easternmost part of the Flin Flon – Snow Lake Mineral Belt (Fig. 2). The Snow Lake assemblage was formed at ca. 1.87–1.85 Ga as the result of intraoceanic accretion (Lewry and Collerson Reference Lewry and Collerson1990; Lucas et al., Reference Lucas, Stern, Syme, Reilly and Thomas1996).
The Stall Lake deposit was the third largest in the Snow Lake area, with ∼6.4 million tonnes total ore containing 4.41 wt.% Cu, 0.5 wt.% Zn, 0.01 wt.% Pb, 12.34 g/t Ag and 1.41 g/t Au (Pehrsson et al., Reference Pehrsson, Gibson and Gilmore2016). It was discovered in 1956 by ground electromagnetics followed by drilling and mined between 1964 and 1982. The host lithology is dominated by gneisses (quartz-eye gneiss, staurolite-garnet gneiss and sillimanite gneiss) and schists (chlorite-garnet schist and chlorite schist), with minor andesitic tuffs and amphibolite sills (Studer, Reference Studer1982).
The ore at the Stall Lake deposit was contained in seven sulfide lenses, each dipping between ∼30 and 60° to the SW, with a maximum length between ∼50 and 300 m, and thickness varying between ∼10 and 50 m (Studer, Reference Studer1982). These lenses were zoned and heterogeneous, dominated by massive coarse-grained pyrrhotite with lesser pyrite, chalcopyrite and sphalerite, minor galena, magnetite and gahnite, and traces of hessite (Ag2Te) and the Bi–Te phase that is the subject of this paper (Studer, Reference Studer1982; our observations; Fig. 3). The main ore mineral, chalcopyrite, is typical, with minor Pb (0.36 wt.%) and Bi (0.09 wt.%); sphalerite contains 6.38 wt.% Fe, 0.37 wt.% Cu and 0.35 wt.% Pb (EMPA, our unpublished data).

Figure 3. Examples of the Stall Lake VMS deposit bismuth tellurides (Bi–Te) in reflected light (a and b) and in back-scattered electron images (c–f), together with other minerals present (Po: pyrrhotite; Py, pyrite; Cpy, chalcopyrite; Mag, magnetite; Hes, hessite, and Amph, amphibole).
Sampling and analytical methodology
Approximately 50 representative drill core samples from the Stall Lake deposit main ore bodies were collected from the Manitoba Geological Survey’s Centennial core storage facility at Baker Lake, ∼20 km south of Flin Flon (Manitoba, Canada), as a part of a larger study of the Flin Flon – Snow Lake Mineral Belt VMS deposits ore mineralogy. Of these, 20 high-grade ore samples, composed mostly of sulfide minerals, were selected for further analyses. Polished thin sections were prepared from these samples and were used for petrographic observations and electron microprobe analyses (EMPA).
Optical microscopy (reflected light) and scanning electron microscopy (SEM) were used to identify the ore minerals present at Stall Lake. The SEM instrument used was a Jeol JSM-6390 at Brandon University, Canada.
Compositions of the Bi–Te phase were obtained using wavelength dispersive spectroscopy (WDS) at the University of Manitoba with a Cameca SX100 electron microprobe. The elements analysed were: Se (X-ray line used La, analytical crystal used TAP); Pb (Ma, LPET); S (Ka, LPET); Sb (La, PET); Ag (La, PET); Bi (Mb, PET); Cu (Ka, LLIF); Ni (Ka, LLIF); Co (Ka, LLIF); Fe (Ka, LLIF); and Te (La, LLIF). The analytical conditions were: 15 keV, 20 nA and 1 mm beam; peak counting times were 20 s (except for Te, Se and S: 30 s), with the background counted for half of that time, once on each side of the peak. The ZAF corrections followed the PAP procedure (Pouchou and Pichoir, Reference Pouchou and Pichoir1984; Pouchou and Pichoir, Reference Pouchou, Pichoir and Armstrong1985). Reference materials are natural and synthetic CdSe (for Se), PbTe (for Pb and Te), pyrite (for S and Fe), stibnite (for Sb), silver (for Ag), bismuthinite (for Bi), chalcopyrite (for Cu), pentlandite (for Ni) and cobalt (for Co).
It was impossible to collect any structural data for the Bi–Te phase from the Stall Lake deposit, due to its small grain size (Fig. 3) and very low quantities. If X-ray diffraction analysis had been performed on a bulk sample, the signal emanating from the small Bi–Te phase grains would be too low to be distinguished from the background noise. Further, micro-Raman is impossible to apply, because the method is not applicable to telluride minerals, which oxidise under the laser beam and thus do not produce a recognisable spectrum or signature (Alexandre and Aisida, Reference Alexandre and Aisida2023). In the absence of structural data, this work relies only on the composition of the Bi–Te phase observed at the Stall Lake deposit.
A significant component of this work is the literature data compilation of the composition of Bi–Te minerals comparable to the Stall Lake phase, used to situate these new data within a relevant compositional context. To that end, all publicly available literature on tsumoite, tellurobismuthite, rucklidgeite and kochkarite are summarised here. These minerals were selected as they are comparable to the Stall Lake phase, on the basis of the Bi/Te (or total cation/total anion) ratio of the structural formula (0.73 to 1.05) and Pb contents. The results of this compilation are provided, together with the Stall Lake phase data, in Table 1 and will be discussed in the Results section; the references used are Rucklidge (Reference Rucklidge1969); Zav’yalov and Begizov (Reference Zav’yalov and Begizov1977); Haranczyk (Reference Haranczyk1978; Spiridonov et al. (Reference Spiridonov, Ershova and Tananaeva1989); McQueen (Reference McQueen1990); McQueen and Solberg (Reference McQueen and Solberg1990); Kase et al. (Reference Kase, Kusachi and Kishi1993); Gu et al. (Reference Gu, Watanabe, Hoshino and Shibata2001); Cepedal et al. (Reference Cepedal, Fuertes-Fuente, Martın-Izard, Gonzalez-Nistal and Rodrıguez-Pevida2006); Cook et al. (Reference Cook, Ciobanu, Wagner and Stanley2007); Ivashchenko et al. (Reference Ivashchenko, Sundblad, Toritsin, Golubev and Lavrov2008); Sejkora et al. (Reference Sejkora, Litochleb, Plasil and Bures2009); Belogub et al. (Reference Belogub, Moloshag, Novoselov and Kotlyarov2011); Voloshin et al. (Reference Voloshin, Cherniavskiy and Voitehobskiy2012); Maslennikov et al. (Reference Maslennikov, Maslennikova, Large, Danyushevsky, Herrington and Stanley2013); Jian et al. (Reference Jian, Lehmann, Mao, Ye, Li, Zhang, Feng and Ye2014); Sabau et al. (Reference Sabau, Berbeleac, Negulescu and Costin2014); Xue and Campbell (Reference Xue and Campbell2014); Stenco et al. (Reference Stenco, Sejkora, Macek, Tucek and Zitnn2014); Roberts (Reference Roberts2017); Haing (Reference Haing2019); Kondratieva et al. (Reference Kondratieva, Anisimova and Kardashevskaia2021); Kasatkin et al. (Reference Kasatkin, Kuznetsov and Arzamastev2022); Alfonso et al. (Reference Alfonso, Colque, Garcia-Valles, Martínez, Yubero, Anticoi and Sidki-Rius2023); Kasprowicz (Reference Kasprowicz2023); and Ma et al. (Reference Ma, Chen, Zhang, Yang, Liu, Lai and Gu2024).
Table 1. Stall Lake Bi–Te minerals. Summary (min, mean, max) of the analytical results for the Stall Lake bismuth telluride compositional groups seen in Figs 4 and 5 (1 to 5) and the literature data compilation for tsumoite, tellurobismuthite, rucklidgeite, and kochkarite (6 to 12 in Fig. 5; based on the references in the text) Cobalt, Ni, Zn, Mn, As, Sn, and Hg were also analysed in the Stall Lake phase, but were always just about or below the quantification limit; Zn, As and Hg were occasionally reported in the literature data, however the values were always very near the quantification limit for the respective method (predominantly electron microprobe), and are not reported here.

Notes: bdl = below detection limit
Analytical results and literature data compilation findings
In this section, the Stall Lake Bi–Te phase is described (Fig. 3) together with its composition, followed by the results of the compilation of all available literature data on the composition of tsumoite, tellurobismuthite, kochkarite and rucklidgeite. The summary data for all minerals, using the minimum, maximum and mean values (both in terms of wt.% and apfu), are provided in Table 1; all data are also represented in Figs 4 and 5.

Figure 4. Se and Pb substitutions in the Stall Lake bismuth tellurides (a and b, respectively), compared to those of the pure bismuth tellurides tsumoite and tellurobismuthite and the Pb-bearing kochkarite and rucklidgeite, based on the references given in the text (Sampling and analytical methodology). The empirically defined Stall Lake compositional groups are numbered 1 to 5. Note that weight percent (wt.%) are used rather than atoms per formula unit (apfu), as minerals with different numbers of atoms in their structural formula (2, 5, 7 and 12) are being compared.

Figure 5. The empirically defined Stall Lake bismuth tellurides compositional groups (numbered 1 to 5), compared to the literature data for tsumoite, tellurobismuthite, kochkarite, and rucklidgeite (numbered 6 to 12), in a Pb substitution vs. Se substitution diagram, based on the references given in the text (Sampling and analytical methodology). The extent of substitution was based on the weight percent (wt.%), even though fairly similar results would have been obtained if atoms per formula unit (apfu) were used, given that the substitutions are expressed as ratios.
Mineral description
On the basis of optical microscope (reflected light; Fig. 3a, b) and SEM observations (Fig. 3c–f), the Stall Lake Bi–Te phase has the following characteristics. The colour is pale yellow to white (Fig. 3a, b). Due to the small grain size, anisotropy was difficult to estimate, though appears to be very low. Its reflectivity is very high, though not as high as galena and is estimated to be ∼80% (Fig. 3a, b). The phase varies from a few μm to ∼30 μm across (Fig. 3). Its morphology varies between mostly isomorphic (Fig. 3e, f) to moderately elongated (Fig. 3d); it is mostly anhedral and occasionally subhedral (Fig. 3c, f).
The Stall Lake Bi–Te phase is associated with the major minerals present at this deposit, mostly with pyrite and pyrrhotite, though also chalcopyrite and magnetite (Fig. 3). It is also sometimes closely associated with hessite (Ag2Te; Fig. 3c, e). On the basis of its textural relationships with other minerals, it probably formed later than the main stage (coarse grained) minerals and together with other (minor and rather small grained) sulfides, possibly during the lower to middle almandine–amphibolite-facies metamorphism that affected the deposit area (maximum conditions 5 kb and 535°C; Froese and Moore, Reference Froese and Moore1980), in imitation of other deposits in the near vicinity (Alexandre et al., Reference Alexandre, Heine, Fayek, Potter and Sharpe2019; Alexandre and Aisida, Reference Alexandre and Aisida2023).
The composition of the Stall Lake deposit Bi–Te phase
At a first glance, this phase is dominated by Bi and Te, with concentrations anywhere between ∼41 and 56 wt.% for the former and between ∼35 and 41 wt.% for the latter, when all analyses are considered. The major cation substitutions are Pb, varying between 3.0 and 11.0 wt.%, Fe (0.2 to 1.4 wt.%), Cu (up to 0.9 wt.%), Ag (up to 3.2 wt.%), and traces of Sb and Au (Table 1). Copper, Ni, Zn, Mn, As, Sn and Hg were also analysed, but were only just at the quantification limit or below. The main anion substitution is Se, varying overall from 0.3 to 7.9 wt.%; S never exceeds ∼0.3 wt.% (Table 1).
When these data are considered from the point of view of the two main substitutions, Pb for Bi and Se for Te (Fig. 4), five clearly discrete compositional groups can be empirically defined (Fig. 5), as described below.
Compositional groups 1 and 2 are fairly similar to each other and are defined by the two lowest amounts of Pb substitutions (∼3.3 wt.% and 5.1 wt.% on average) and the lowest Se substitution (∼0.3 wt.% for both groups, on average; Figs 4 and 5). Iron and Cu are present in both groups (∼1.2 wt.% for Fe, ∼0.7 and ∼0.6 wt.% for Cu); Au, Ag and Sb are present in trace amounts (Table 1). Their average structural formulae, on the basis of two atoms in imitation of tsumoite, can be written as (Bi0.841Fe0.068Pb0.051Cu0.034)Σ0.996(Te0.981Se0.012S0.011)Σ1.004 for group 1 and (Bi0.841Pb0.080Fe0.069Cu0.031)Σ1.024(Te0.952Se0.012 S0.012)Σ0.976 for group 2.
Compositional group 3 is notable for increased Pb and Se substitutions, with average Pb of ∼7.2 wt.% and average Se of ∼0.75 wt.% (Table 1). There is less Fe (∼0.9 wt.%), Ag (∼0.3 wt.%) and Cu (∼0.1 wt.%); Sb and Au are present at trace level (Table 1). Sulfur is relatively low, at ∼0.1 wt.%. The average structural formula for this compositional group can be written, on the basis of 12 atoms as in kochkarite, as (Bi4.034Pb0.654Fe0.307Cu0.039)Σ5.089(Te6.674Se0.179S0.056)Σ6.909.
Compositional groups 4 and 5 contain higher Pb (∼10.9 wt.% on average) than group 3, however they are remarkable by the largest degree of Se substitutions, with, on average, ∼4.9 wt.% for group 4 and ∼7.2 wt.% for group 5 (Table 1; Fig. 4). Though Se in these two groups is indeed high, it does not reach the levels of poubaite (∼13.6 wt.%; Anthony et al., Reference Anthony, Bideaux, Bladh and Nichols1990; Johan et al., Reference Johan, Picot and Ruhlmann1987; Kalinin, Reference Kalinin2021). There is less Fe (∼0.4 wt.%) and Cu (∼0.2 and ∼0.3 wt.%, respectively) than the other groups, but more Ag (∼0.6 and ∼0.7 wt.%, on average; Table 1); Sb and Au are present at trace level. Sulfur is also the lowest among all Stall Lake analyses, at ∼0.1 wt.%; Table 1). The average structural formula for these two compositional groups can be written, also based on 12 atoms, as (Bi3.819Pb0.973Fe0.140Cu0.045)Σ5.089(Te5.722Se1.154S0.034)Σ6.911 for group 4 and (Bi3.775Pb0.958Fe0.138Cu0.076)Σ5.064 (Te5.217Se1.655S0.064)Σ6.936 for group 5.
The main point of these observations is the continuous increase in both Pb and Se substitution from group 1 to group 5, without Se reaching poubaite levels (Fig. 5). Iron and Cu are always present in measurable amounts but tend to decrease with increasing degree of Pb and Se substitution; Ag is the highest in groups 4 and 5, whereas Au and Sb are always present at trace amounts (Table 1). It should also be noted that the average Bi/Te (or total cations/total anions) ratio for groups 1 and 2 is close to 1, whereas in the other three groups this ratio is closer to 0.74 (Table 1), justifying using the tsumoite and kochkarite structural formula patterns (2 and 12 atoms, respectively), as indicated above.
Literature data compilation
Tsumoite literature data fall into two compositional groups. The first corresponds to a fairly typical tsumoite with 62.3 wt.% Bi and 36.4 wt.% Te, and relatively low degrees of cation and anion substitutions (∼1.6 wt.% and ∼0.33 wt.% in total, respectively; Table 1). The main substitution cations are Pb, Fe, Sb and Ag, all at ∼0.4 wt.%; Au is also present (∼0.1 wt.%). Its structural formula can be written as (Bi0.984Fe0.027Sb0.010Sb0.009)Σ1.038(Te0.943S0.010Se0.009)Σ0.962.
The second compositional group for tsumoite, observed in two deposits (from the Ural mountains and from the Karelian Peninsula in Russia (Belogub et al., Reference Belogub, Moloshag, Novoselov and Kotlyarov2011; Ivashchenko et al., Reference Ivashchenko, Sundblad, Toritsin, Golubev and Lavrov2008, is characterised by the presence of Pb substitution (∼4.6 wt.% of Pb on average, ranging from 3.2 to 5.9 wt.%) and was labelled Pb-rich tsumoite. No other substitution cations or anions were observed (Table 1); its structural formula can be written as (Bi0.927Pb0.075)Σ1.002Te0.998.
Kochkarite literature data also define two discrete two compositional groups, distinguished from each other by their Pb concentrations (Figs 4 and 5, Table 1). The major group corresponds to typical kochkarite, with, on average, ∼43.7 wt.% Bi, 11.8 wt.% Pb, 1.3 wt.% Ag, and 1.0 wt.% Fe; Sb and Cu are also present (Table 1). The anions are dominated by Te (∼45.0 wt.%); the amount of Se is relatively low (0.2 wt.%; Table 1). Significantly, the data show a high degree of variation of all elements present, with Pb, for instance, varying from ∼9.5 to ∼16.0 wt.% (Table 1). The average structural formula for this group is (Bi3.645Pb1.055Fe0.346Ag0.231)Σ5.365(Te6.539S0.049Se0.047)Σ6.635.
The second kochkarite group, observed at the Deer Horn intrusion-related polymetallic deposit in British Columbia, Canada (Roberts, Reference Roberts2017), contains more Bi (∼43.7 wt.% on average), less Pb (∼7.0 wt.% on average), some Ag (∼1.0 wt.%), Sb (∼0.5 wt.%) and traces of Fe, Au and Cu (Table 1); its anion substitutions are very low (Se ≈ 0.03 wt.% on average, no S). This compositional group was labelled ‘low-Pb kochkarite’ and its average structural formula is (Bi4.031Pb0.648Ag0.170Sb0.086)Σ4.962(Te7.029Se0.007)Σ7.038.
Rucklidgeite is another mineral for which two distinct sub-populations were observed in the compilation dataset; it is also the one with the highest variation in the concentrations of most elements (Fig. 4b; Table 1). Most analyses form a group of more or less typical rucklidgeite, with ∼40.0 wt.% Bi on average, ∼16.1 wt.% Pb (varying from ∼10.8 to ∼20.9 wt.%), ∼1.6 wt.% Ag (∼0.4 to ∼4.7 wt.%), and some Fe and Sb (∼0.5 wt.% for both; Table 1). This group contains ∼44.4 wt.% Te on average, and limited Se (∼0.3 wt.% on average, ranging from below detection to ∼1.1 wt.%); S is minimal. Its average structural formula is (Bi2.098Pb0.709Ag0.159Fe0.104)Σ3.124(Te3.820Se0.037S0.007)Σ3.876.
A minor subpopulation of rucklidgeite, from a VMS deposit in Quebec, Canada (Kasprowicz, Reference Kasprowicz2023), is characterised by high Se (∼3.1 wt.% on average: Table 1) and was labelled high-Se rucklidgeite. The concentrations of all other elements in this group are similar to those in the main rucklidgeite population (Table 1). Its average structural formula is (Bi1.997Pb0.829Ag0.140Fe0.097)Σ3.114(Te3.444Se0.414S0.027)Σ3.886.
Tellurobismuthite also shows some compositional variation, in particular in its Se concentrations, which vary from below detection to ∼2.7 wt.% (∼0.45 on average; Table 1). The average Bi concentrations are ∼51.9 wt.%; Te is present at ∼47.1 wt.% on average (Table 1). Tellurobismuthite has the lowest Pb concentrations of all the minerals considered (∼0.3 wt.% on average); Fe, Sb, Cu and Ag are also present (∼1.0 wt.%, ∼0.4, wt.%, ∼0.4 wt.%, and ∼0.2 wt.%, on average, respectively). The average structural formula of tellurobismuthite, based on this compilation, is (Bi1.899Fe0.130Cu0.046Ag0.011)Σ2.120(Te2.822Se0.043S0.015)Σ2.880.
Discussion
The literature compilation of all publicly available compositional data for tsumoite, tellurobismuthite, kochkarite and rucklidgeite (based on the references at the end of the section on Sampling and analytical methodology), provide some interesting insights and some general comments can be made about their chemical purity and their compositional variability; these questions have never been addressed. The composition of the Stall Lake deposit phases will also be discussed, in the context of the other bismuth tellurides (e.g. Cook et al., Reference Cook, Ciobanu, Wagner and Stanley2007; Gu et al. Reference Gu, Watanabe, Hoshino and Shibata2001).
Extent of substitutions in bismuth tellurides
All minerals considered and documented by both the literature compilation and the new data from Stall Lake contain at least some Pb, Fe, Cu, Sb, Ag and Au (Table 1). Silver is most abundant in rucklidgeite (∼1.6 wt.%) and kochkarite (∼1.4 wt.%); Fe is most abundant in kochkarite and tellurobismuthite (∼1.0 wt.% on average). Tellurobismuthite also has the most Sb and Cu, both at ∼0.4 wt.% on average. The average concentrations of Pb vary from ∼0.3 wt.% in tellurobismuthite to ∼16.1 wt.% in high-Se rucklidgeite. Overall, the cations observed in these minerals, other than Bi and Pb, typically combine to ∼2 wt.%.
All minerals considered show some degree of Se substitution, with Se varying from below detection limit to ∼3.6 wt.% for high-Se rucklidgeite (VMS deposit in Quebec; Kasprowicz, Reference Kasprowicz2023). The amounts of Se are low for tsumoite and kochkarite (both at ∼0.2 wt.% on average), higher for tellurobismuthite (∼0.5 wt.% on average, up to 2.7 wt.%; Cook et al., Reference Cook, Ciobanu, Wagner and Stanley2007), and highest for high-Se rucklidgeite (∼3.1 wt.% on average; Kasprowicz, Reference Kasprowicz2023; Table 1). Overall, Se and S typically combine to just under 1 wt.% in these minerals.
From this perspective, it could be argued that there is no such thing as pure bismuth telluride (in clear contradiction of the first three words of this paper): all and any of them will always contain a certain amount of both cation and anion substitutions, as demonstrated by the literature data compilation and the new data obtained at Stall Lake. The cation substitutions are mostly Pb, but also, in decreasing amounts, Ag, Fe, Sb, Cu and Au; the anion substitutions are predominantly Se and much less S (Table 1). As a practical corollary, these are the bare minimum of elements that should always be analysed in any bismuth telluride.
Compositional variability of bismuth tellurides
Another important observation originating from the literature data compilation is the significant compositional variability observed for all the minerals concerned, to the point of defining distinct sub-groups (high-Pb tsumoite, low-Pb kochkarite, high-Se rucklidgeite; Fig. 5). That variability, in particular for the Pb concentrations, is the highest for kochkarite and rucklidgeite: in fact, their Pb concentrations overlap to a very high extent (Fig. 5). This observation raises questions about the chemical definition of these two minerals (and also about labelling them in any particular investigation, as discussed below; see also Cook et al., Reference Cook, Ciobanu, Wagner and Stanley2007 and Gu et al., Reference Gu, Watanabe, Hoshino and Shibata2001). For instance, it should be possible to distinguish kochkarite and rucklidgeite from each other by the Bi/Pb ratio in their structural formulae as this is ideally 0.750 for rucklidgeite ((Bi,Pb)3Te4) and 0.714 for kochkarite (PbBi4Te7) (Anthony et al., Reference Anthony, Bideaux, Bladh and Nichols1990). However, the average total cation/total anion ratios resulting from this compilation are on the order of 0.80 for kochkarite, rucklidgeite and high-Se rucklidgeite, and ∼0.71 for low-Pb kochkarite (Table 1): these are significantly different from the theoretical values.
Labelling, or naming, of bismuth tellurides
A direct and practical consequence of the above comments is that sometimes it might be problematic, in the absence of structural information, to know which bismuth telluride is observed when new data are collected. The present investigation is a pertinent example of this difficulty.
The Stall Lake deposit bismuth tellurides fall into five empirically defined compositional groups, as described above (Mineral description; Table 1; Fig. 5), and the question that emerged immediately was to what minerals these compositions correspond (again, in the absence of structural data). This question is most easily answered for compositional groups 1 and 2, which seem to closely resemble the high-Pb tsumoite visible in the literature data compilation (Fig. 5; Table 1). Specifically, (1) the cation/anion ratios in their structural formulas are very similar, at ∼1 (no other bismuth tellurides have the same ratio; telluronevskite is not an option because of its high Se; Anthony et al., Reference Anthony, Bideaux, Bladh and Nichols1990); and (2) the amounts of Pb in the Stall Lake groups 1 and 2 (∼3.3 and ∼5.1 wt.%) are fairly comparable to those of literature data high-Pb tsumoite (∼4.6 wt.%; Table 1; VMS deposits in the Ural, Russia: Belogub et al., Reference Belogub, Moloshag, Novoselov and Kotlyarov2011). It seems that these two compositional groups can be regarded as high-Pb tsumoite (e.g. Cook et al., Reference Cook, Ciobanu, Wagner and Stanley2007).
The other three compositional groups from Stall Lake are more difficult to label. The total cation/total anion ratio of their average structural formulae are ∼0.74, which is closest to theoretical rucklidgeite (0.75; Anthony et al., Reference Anthony, Bideaux, Bladh and Nichols1990). However, the literature compilation (Table 1; references at the end of Sampling and analytical methodology) reveals that kochkarite and rucklidgeite have cation/anion ratios between 0.71 (low-Pb kochkarite) and 0.81 (kochkarite, rucklidgeite), indicating that it may be indeed challenging to decide what the Stall Lake compositions are.
Based on its amount of Pb, the Stall Lake deposit compositional group 3 (Pb ≈ 7.2 wt.%) is closest to low-Pb kochkarite (compilation group 8: Pb ≈ 7.0 wt.%; Roberts, Reference Roberts2017; Table 1); the amount of Pb is significantly lower than that of rucklidgeite (ranging from ∼10.8 to ∼20.9 wt.%; references at the end of Sampling and analytical methodology). Hence, this bismuth telluride compositional group can probably be named low-Pb kochkarite.
The greatest difficulty is deciding what the Stall Lake compositional groups 4 and 5 are, the high-Se bismuth tellurides. The average cation/anion ratio of their structural formulas (∼0.73) are not sufficiently close to either literature data kochkarite or rucklidgeite (∼0.8) to be of any help. Based uniquely on the Pb contents of these two groups (at ∼10.9 wt.%; Table 1), kochkarite was preferred, as its Pb content is ∼11.8 wt.% on average (Table 1). This is not ideal and serves to illustrate the difficulties in naming unknown bismuth tellurides.
The high-Se bismuth telluride phases from Stall Lake
The difficulty described above is compounded by the very high degree of Se substitution observed in the Stall Lake deposit compositional groups 4 and 5 (average Se ≈ 4.9 and ≈ 7.2 wt.%, respectively; Figs 4 and 5), which had never been observed in either kochkarite or rucklidgeite. These concentrations are lower than in poubaite (Se ∼13.6 wt.% on average; Anthony et al., Reference Anthony, Bideaux, Bladh and Nichols1990; Kalinin, Reference Kalinin2021; Johan et al., Reference Johan, Picot and Ruhlmann1987); poubaite also has more Pb (∼20.1 wt.% on average). The high amounts of Se in the Stall Lake compositional groups 4 and 5 clearly demonstrate the potential for anion substitutions in bismuth tellurides, as noted above, however this also raises the question of what mineral(s) these two groups correspond to.
The two most likely hypotheses are (1) the existence of new minerals or (2) the existence of a high-Se kochkarite or rucklidgeite. New minerals with these compositions would have the ideal structural formulas of Bi4PbTe6Se (group 4) and Bi4PbTe5Se2 (group 5). However, such a hypothesis cannot be addressed here because of the absence of structural data, highlighting in the process the absolute necessity to collect, when at all possible, crystal structure data for any mineral studied. The second hypothesis seems more likely and, on the basis of the Pb contents of these compositional groups (both are, on average, ∼10.9 wt.%; Table 1) the term ‘high-Se kochkarite’ was preferred for them. This labelling is far from satisfactory or conclusive and should be taken with a pinch of scepticism.
One corollary of the high Se observed in the Stall Lake bismuth tellurides is speculation regarding a possible solid solution between kochkarite and poubaite or between rucklidgeite and poubaite (e.g. Kase et al., Reference Kase, Kusachi and Kishi1993). The former possibility can be dismissed immediately, as kochkarite is cubic (point group 4/m 3̅ 2/m), whereas poubaite is hexagonal (point group 3̅ 2/m; Anthony et al., Reference Anthony, Bideaux, Bladh and Nichols1990). The latter possibility remains theoretically possible, as both poubaite and rucklidgeite are hexagonal, with the same point group (3̅ 2/m). They also have comparable amounts of Pb (∼20.1 for poubaite and ∼16.1 for high-Se rucklidgeite) and the same cation/anion ratio of 0.75 (Anthony et al., Reference Anthony, Bideaux, Bladh and Nichols1990). In such a scenario, the high-Se phases from Stall Lake (compositional groups 4 and 5) could be considered as belonging to this hypothetical solid solution if they are in fact (low-Pb) rucklidgeite. The high-Se literature data for rucklidgeite (compositional group 11; Kasprowicz, Reference Kasprowicz2023) would also belong to such a solid solution (Fig. 6). However, the lesser amount of Pb in Stall Lake compositional groups 4 and 5 compared to poubaite, and the far from consistent or continuous Te–Se anti-correlation (Fig. 6) are arguments against such a possibility. It can be concluded that while a poubaite–rucklidgeite solid solution might be theoretically possible, the present work, specifically the existence of a high-Se bismuth telluride compositions at Stall Lake, is not a sufficient proof.

Figure 6. A speculative rucklidgeite–poubaite solid solution, visualised by the Se substitution of Te in rucklidgeite, high-Se rucklidgeite, the two high-Se Stall Lake compositional groups, and poubaite. The atoms per formula unit (apfu) were calculated on the basis of 7 atoms, for a generic formula (Bi,Pb)3(Te,Se)4, which describes both rucklidgeite and poubaite. The Se–Te anti-correlation is far from perfect in terms of consistency and connectivity, which, combined with the low Pb concentrations in the Stall Lake bismuth telluride and the lack of crystal structure data, renders the existence of such a solid solution currently unsubstantiated.
Summary
A few interesting generalisations can be made from the literature data compilation of four bismuth tellurides and new data obtained at the Stall Lake VMS deposit, the first is that these minerals have a highly variable composition (Fig. 4, Table 1). They all contain at least a modest degree of cation and anion substitutions, with (in decreasing concentrations) Pb, Ag, Fe, Cu, Sb and Au substituting for Bi, and Se and S substituting for Te. All these elements should be therefore routinely analysed for in bismuth tellurides. The amount of substitution cations other than Bi and Pb is in the vicinity of 2 wt.%, Se and S combine to under 1 wt.%.
Bismuth tellurides are not only more chemically complex than previously appreciated, but demonstrate fairly high compositional variability, to the point of defining specific mineral varieties such as high-Pb tsumoite, low-Pb kochkarite and high-Se rucklidgeite (Fig. 5). Their empirical compositions are commonly different from their ideal structural formulae (Fig. 5). The highest variability is observed in the amount of Pb in kochkarite and rucklidgeite, with the former varying between ∼9.5 and ∼16.0 wt.%, and the latter between ∼10.8 and ∼20.9 wt.% (Table 1). These two minerals’ compositions overlap to such an extent that distinguishing them from each other conclusively might, in some cases, only be possible using crystal structural information. In practical terms, this compositional variability can sometimes represent a difficulty with the classification of an unknown bismuth telluride, as this investigation demonstrates.
Two new high-Se bismuth telluride compositions were discovered at the Stall Lake VMS deposit, containing, on average, ∼43.4 wt.% Bi, 10.9 wt.% Pb, 39.6 and 36.8 wt.% Te, and 4.9 and 7.2 wt.% Se (Table 1; Fig. 5).These two compositional groups were tentatively labelled high-Se kochkarite. Such a high extent of Se substitution had never been observed for kochkarite or rucklidgeite and opens a few possibilities, one of which is the existence of two new minerals, with the ideal structural formulae of Bi4PbTe6Se and Bi4PbTe5Se2. It can be speculated that the Stall Lake high-Se bismuth telluride belongs to a hypothetical rucklidgeite–poubaite solid solution (Fig. 6), even though the available data, from literature or the newly obtained, do not conclusively support such a possibility.
As defining a new mineral or a solid solution would strongly rely on crystal structure data, which were not available for the Stall Lake phases, these speculations cannot be fully addressed at present. However, these new data, showing high and variable Pb and Se substitutions, emphasise the high compositional variability and the high extent of both cation and anion substitutions in bismuth tellurides, which is the main significance of this contribution.
Acknowledgements
T. Heine, from Flin Flon in Canada, assisted with sample collection. A. Baig, from Brandon University, assisted with sample preparation and SEM microscopy. P. Yang, from the University of Manitoba, assisted with EMPA. S. Kienas and R. Schira, from Brandon University’s library reference team, assisted with finding references for bismuth tellurides. The reviewers, Dr. Belogub and an anonymous one, as well as the Principal Editor Dr. Mitchell and the Associate Editor Dr. Good, significantly contributed to the quality of this paper through constructive and helpful comments and suggestions.
Competing interests
The authors declare none.