Introduction
The Bukov mine (Rožná–Olší uranium district, Czech Republic) represents the most significant selenide mineralisation in the Bohemian Massif in terms of size of the selenide (especially berzelianite and umangite) aggregates and their abundance (Kvaček, Reference Kvaček1973, Reference Kvaček1979). At present, seventeen selenide species are reported from the Bukov mine (e.g. https://www.mindat.org/loc-757.html), including descriptions of the new minerals bukovite (Johan and Kvaček, Reference Johan and Kvaček1971) and sabatierite (Johan et al., Reference Johan, Kvaček and Picot1978). During the systematic investigation of selenides from the various occurrences of the Bohemian Massif, we identified a possible Se analogue of ullmannite in a Ni–Co rich selenide association from the Bukov mine. Thus, a crystal-chemical investigation was undertaken, which confirmed its identity with synthetic NiSbSe and led to the proposal of the new mineral species, kvačekite.
This new mineral and its name were approved by the Commission on New Minerals, Nomenclature and Classification of the International Mineralogical Association (IMA2023–095, Pauliš et al., Reference Pauliš, Dolníček, Sejkora, Pour, Laufek, Ulmanová and Vymazalová2024). Kvačekite is named after Milan Kvaček (13. 2. 1930 to 21. 5. 1993), a Czech mineralogist and the author of the first comprehensive descriptions of selenides from western Moravia (Kvaček, Reference Kvaček1973, Reference Kvaček1979). Milan Kvaček is the author of over sixty publications in the fields of mineralogy, chemistry and geochemistry. He also co-operated closely with Zdeněk Johan from BRGM (Orléans) on the description of six new selenides from the Bohemian Massif – bukovite, fischesserite, hakite, permingeatite, petrovicite and sabatierite. Its mineral symbol, in accord with Warr (Reference Warr2021), is Kvč. The holotype material (polished section) is deposited in the mineralogical collection of the Department of Mineralogy and Petrology of the National Museum, Prague, Czech Republic (catalogue number P1P 26/2023).
Occurrence and mineral description
Occurrence
The new mineral was discovered in material from Milan Kvačeḱs research collection originating from the now abandoned Bukov mine, located on the northeastern edge of the village of Bukov, 35 km southeast of Nové Město na Moravě, Vysočina Region, Czech Republic. The GPS coordinates of the type occurrence are 49°27'27.333“N, 16°13'43.224”E.
The Bukov mine is a part of the Rožná–Olší uranium ore district with a total mine production of 23,000 t U in the years 1957–2017 (Kříbek et al., Reference Kříbek, Žák, Dobeš, Leichmann, Pudilová, René, Scharm, Scharmová, Hájek, Holeczy, Hein and Lehmann2009, Reference Kříbek, Knésl, Dobeš, Veselovský, Pořádek, Škoda, Čopjaková, Leichmann and Košek2022). This uranium ore district is an example of shear-zone hosted, late Variscan and post-Variscan hydrothermal mineralisation (Kříbek et al., Reference Kříbek, Žák, Dobeš, Leichmann, Pudilová, René, Scharm, Scharmová, Hájek, Holeczy, Hein and Lehmann2009). Uranium mineralisation is bound to a system of NNW–SSE trending faults and zones of cataclasites, which are hosted in high-grade metamorphic rocks with a Moldanubian Zone of the Bohemian Massif. Mineralised cataclasites are 4–15 m thick and are traceable over a horizontal distance of 15 km. The mineralised zones strike 340–355° and dip WSW at angles of 45–70° (Kříbek et al., Reference Kříbek, Knésl, Dobeš, Veselovský, Pořádek, Škoda, Čopjaková, Leichmann and Košek2022). The host rocks mainly comprise gneisses and amphibolites, with small bodies of calc-silicate gneiss, marble, serpentinite and pyroxenite. The vein and disseminated mineralisation at Rožná comprises (1) pre-uranium quartz–sulfide and carbonate–sulfide mineralisation, (2) uranium, and (3) post-uranium carbonate–quartz–sulfide mineralisation (Kříbek et al., Reference Kříbek, Žák, Dobeš, Leichmann, Pudilová, René, Scharm, Scharmová, Hájek, Holeczy, Hein and Lehmann2009, Reference Kříbek, Knésl, Dobeš, Veselovský, Pořádek, Škoda, Čopjaková, Leichmann and Košek2022). Uranium mineralisation in some carbonate veins is accompanied by selenides, which are generally slightly younger than the main uranium mineralisation (Kříbek et al., Reference Kříbek, Žák, Dobeš, Leichmann, Pudilová, René, Scharm, Scharmová, Hájek, Holeczy, Hein and Lehmann2009). The most significant occurrences of selenides appear in calcite–uraninite veins in the upper levels (5th to 9th) of the Bukov mine in the underlying rocks of the uranium-bearing shear zone R1. Selenide mineralisation occurs in two distinct types. The first is represented by massive aggregates of predominant berzelianite up to several tens of cm in size, usually in association with Cu sulfides (chalcopyrite, bornite and chalcocite). The second type represents disseminated selenide grains, more rarely tiny veinlets, formed by diverse selenides of often microscopic dimensions (Kvaček, Reference Kvaček1973, Reference Kvaček1979). The full list of minerals (more than 110 species including 20 selenides) recorded from this area is given on https://www.mindat.org/loc-27262.html.
Kvačekite was identified in samples of white calcite gangue in association with nickeltyrrellite–tyrrellite, berzelianite, hakite-(Zn), hakite-(Cd), eucairite, clausthalite and gold.
Experimental procedures
The small size of kvačekite grains and its intergrowths with other minerals prevented its extraction and isolation in an amount sufficient for the relevant crystallographic and structural investigations. Therefore, these measurements were performed on synthetic NiSbSe.
The synthetic NiSbSe phase was prepared using an evacuated silica glass tube method in the Laboratory of Experimental Mineralogy of the Czech Geological Survey in Prague. Nickel (99.999%), antimony (99.999%) and selenium (99.9999%) were used as starting materials for synthesis. The evacuated tube with its charge was sealed and annealed at 800°C. After cooling by a cold-water bath, the charge was ground into powder in acetone using an agate mortar and pestle, and thoroughly mixed to homogenise. The pulverised charge was sealed in an evacuated silica-glass tube again, and heated at 300°C for 3 months. The experimental product was then rapidly quenched in cold water. The size of obtained NiSbSe grains ranges from 2 to 20 μm.
Physical and optical properties
Kvačekite occurs as anhedral grains up to 15 μm in size (Fig. 1), but usually only 2–5 μm across, in nickeltyrrellite, usually intergrown with hakite-(Zn) or clausthalite. Synthetic kvačekite is greyish white in colour with a grey streak and opaque in transmitted light; it has a metallic lustre. A cleavage was not observed; it is brittle with a conchoidal fracture. The calculated density for the empirical formula (Z = 4) is 7.607 g.cm–3. Mohs hardness is assumed to be 5–6, in agreement with other members of the cobaltite group. In reflected light, kvačekite is white with a faint yellowish shade; bireflectance, pleochroism and anisotropy are absent. Internal reflections were not observed. Reflectance spectra for kvačekite and synthetic NiSbSe were measured in air with a TIDAS MSP400 spectrophotometer attached to a Leica microscope (100× objective) using a WTiC (Zeiss No. 370) standard, with a square sample measurement field of ca. 4 × 4 μm. The results for the 400–700 nm range are given in Table 1 and plotted in Fig. 2.
a The reference wavelengths required by the Commission on Ore Mineralogy (COM) are given in bold.
Chemical composition
Chemical analyses of kvačekite from Bukov and synthetic NiSbSe were performed using a Cameca SX100 electron microprobe operating in wavelength-dispersive mode (25 kV, 20 nA and 0.7 μm beam size). The following standards and X-ray lines were used to minimise line overlaps: Ag (AgLα), Au (AuMα), Bi2Se3 (BiMβ), CdTe (CdLα), chalcopyrite (CuKα and SKα), Co (CoKα), GaAs (GaLα), Ge (GeLα), HgTe (HgMα), Mo (MoLα), Ni (NiKα), NiAs (AsLβ), PbS (PbMα), PbSe (SeLβ), PbTe (TeLα), pyrite (FeKα), Sb2S3 (SbLα), Sn (SnLα), Tl(Br,I) (TlLα) and ZnS (ZnKα). Peak counting times were 20 s for all elements, and 10 s for each background. Arsenic, Ag, Au, Bi, Cd, Fe, Ga, Ge, Hg, Mo, Pb, Sn, Te, Tl and Zn were all found to be below the detection limits (0.02–0.15 wt.%). Raw intensities were converted to the concentrations of elements using the automatic ‘PAP’ (Pouchou and Pichoir, Reference Pouchou, Pichoir and Armstrong1985) matrix-correction procedure.
Analytical data for kvačekite from Bukov (17 analyses) and synthetic NiSbSe (6 analyses) are given in Table 2. On the basis of 3 atoms per formula unit (apfu), the empirical chemical formula of kvačekite is (Ni0.95Cu0.04Co0.03)Σ1.02Sb1.00(Se0.97S0.01)Σ0.98 and Ni1.00Sb1.02Se0.98 for synthetic NiSbSe. The ideal formula is NiSbSe, which requires Ni 22.63, Sb 46.93 and Se 30.44, a total of 100.00 wt.%.
(σ) – estimated standard deviation; n = number of spot analyses.
Raman spectroscopy
The Raman spectra of kvačekite from Bukov and synthetic NiSbSe were collected in the range 40–4000 cm–1 using a DXR dispersive Raman Spectrometer (Thermo Scientific) mounted on a confocal Olympus microscope. The Raman signal was excited by an unpolarised green 532 nm solid state, diode pumped laser and detected by a CCD detector. The experimental parameters were: 100× objective, 10 s exposure time, 100 exposures, 50 μm slit spectrograph aperture and 1.5 mW laser power level. The eventual thermal damage of the measured points was excluded by visual inspection of the excited surface after measurement, by observation of possible decay of spectral features in the start of excitation and checking for thermal downshift of Raman lines. The instrument was set up by a software-controlled calibration procedure using multiple neon emission lines (wavelength calibration), multiple polystyrene Raman bands (laser frequency calibration) and standardised white-light sources (intensity calibration). Spectral manipulations were performed using the Omnic 9 software (Thermo Scientific). The Raman bands of kvačekite appear only in the area 250–40 cm–1, the observed main bands at 199, 179 and 146 cm–1, with weak to moderate intensity bands at 127, 102, 87, 68, 57 and 45 cm–1 corresponding very well with Raman bands of synthetic NiSbSe (Fig. 3). Strong bands in the area 250–120 cm–1 and medium intensity bands in the range 120–70 cm–1 can be assigned to stretching and bending vibrations of Sb(Ni)–Se bonds. Weak Raman bands below 70 cm–1 should be connected with external modes – translations and librations (Škácha et al., Reference Škácha, Buixaderas, Plášil, Sejkora, Goliáš and Vlček2014; Sejkora et al., Reference Sejkora, Buixaderas, Škácha and Plášil2018).
X-ray diffraction data
The small size of natural kvačekite prevented its extraction and direct investigation by means of X-ray diffraction. Therefore, powder X-ray diffraction data were collected on the synthetic phase NiSbSe at room temperature using a Bruker D8 Advance diffractometer equipped with a solid-state LynxEye detector and secondary monochromator producing CuKα radiation, housed at the Department of Mineralogy and Petrology, National Museum, Prague, Czech Republic. The instrument was operating at 40 kV and 40 mA. In order to minimise the background, the powder was placed on the surface of a flat silicon wafer. The powder pattern was collected in the Bragg–Brentano geometry in the range 3–75°2θ, step size 0.01° and counting time of 20 s per step. The positions and intensities of diffractions were found and refined using the Pearson VII profile-shape function of the ZDS program package (Ondruš, Reference Ondruš1993). The obtained experimental data correspond well with the powder X-ray diffraction pattern (Table 3) calculated using the software PowderCell2.3 (Kraus and Nolze, Reference Kraus and Nolze1996) on the basis of our structural model of NiSbSe. The unit-cell parameters for cubic space group P213 (# 198) were refined by the least-squares program of Burnham (Reference Burnham1962) as follows: a = 6.09013(13) Å, V = 225.881(15) Å3 and Z = 4. The refined unit-cell parameters agree very well with data a = 6.0868(6) Å (Foecker and Jeitschko, Reference Foecker and Jeitschko2001) and a = 6.090(3) Å (Hulliger, Reference Hulliger1963) given for synthetic NiSbSe.
Intensity and d hkl (in Å) were calculated using the software PowderCell2.3 (Kraus and Nolze, Reference Kraus and Nolze1996) on the basis of our refined structural model. Only reflections with I calc > 0.5 are listed. The six strongest reflections are given in bold.
Crystal structure
The source data for Rietveld refinement of the synthetic phase NiSbSe were collected at room temperature using a Bruker D8 Advance diffractometer equipped with LynxEye XE detector and CuKα radiation (Czech Geological Survey, Prague). The instrument was operated at 40 kV and 30 mA. In order to minimise the background, the powder was placed on the surface of a flat silicon wafer. The powder pattern was collected in the Bragg–Brentano geometry in the range 10–140°2θ, step 0.015° and counting time of 1.4 s per step. The details of data collection and basic crystallographic facts are given in Table 4.
Rietveld refinement was carried out using the Topas 5 program and involved refinement of unit-cell parameter, atomic coordinates, isotropic size and strain parameters. Isotropic displacement parameters were fixed during the refinement (B iso = 0.2 Å2). The structure model proposed by Hulliger (Reference Hulliger1963) and subsequently by Foecker and Jeitschko (Reference Foecker and Jeitschko2001) for synthetic NiSbSe was confirmed, the refined coordinates did not deviate more than 0.002 from their starting positions. The refined unit-cell parameter a = 6.0902(2) Å agrees very well with the data a = 6.0868(6) Å (Foecker and Jeitschko, Reference Foecker and Jeitschko2001) and a = 6.090(3) Å (Hulliger, Reference Hulliger1963) given for synthetic NiSbSe. The refinement converged to R wp = 0.056, R B = 0.067 and the final Rietveld plot is shown in Fig. 4. The atomic positions are summarised in Table 5. The crystallographic information file has been deposited with the Principal Editor of Mineralogical Magazine and is available as Supplementary material (see below).
a Isotropic displacement parameters were fixed during the refinement.
Synthetic NiSbSe crystallises in the cubic space group P213 and has the ullmannite type structure. It can be considered as a ternary ordered variant of the pyrite-type structure MX2 (Pa $\bar{3}$), where Sb and Se atoms replace S atoms and are ordered at two independent crystallographic positions X and Y. Consequently, the symmetry is lowered from Pa $\bar{3}$ (pyrite type) to P213 (ullmannite type). Derivation of the kvačekite crystal structure from the pyrite-type and corresponding group–subgroup relation is shown in Fig. 5. The disordered (i.e. pyrite-type) structure model with the Sb/Se mixed site has been also tested in the Rietveld refinement. As indicated in Fig. 6, this structured model yields significantly worse profile agreement factors and cannot fit all observed reflections in the powder diffraction pattern of synthetic NiSbSe. Therefore, the disordered pyrite-type structure can be excluded.
Another ordering of anion atoms occurs in the cobaltite-type structure (CoAsS), which shows orthorhombic symmetry (Pca21). No indication of peak splitting indicating orthorhombic symmetry have been observed for the synthetic analogue of kvačekite.
The crystal structure of kvačekite is formed by corner-sharing [NiSb3Se3] octahedra (Fig. 7). There are three Ni octahedra sharing every anion. The Sb and Se atoms form covalent pairs (i.e. dumbells), with bonding distance of 2.665 Å. A similar bonding distance of 2.636 Å in Sb–Se covalent pairs was observed in milotaite, PdSbSe (Paar et al., Reference Paar, Topa, Makovicky and Culetto2005). This is slightly longer than observed Ni–Se and Ni–Sb distances of 2.469(3) and 2.579(3) Å, respectively. Considering Ni–Se distances, slightly longer NiSe separations of 2.503 Å were observed in synthetic NiSe (Unoki et al., Reference Unoki, Yoshiasa, Kitahara, Nishiayama, Tokuda, Sugiyama and Nakatsuka2021), an analogue of sederholmite, which adopts a nickeline-type structure. Ullmannite shows a Ni–Sb distance of 2.563 Å (Foecker and Jeitschko, Reference Foecker and Jeitschko2001) which can be compared with the Ni–Sb separations in the kvačekite structure (2.579(3) Å). In addition to ullmannite, the most closely related mineral species is milotaite, PdSbSe (Paar et al., Reference Paar, Topa, Makovicky and Culetto2005).
Electron back-scattered diffraction
The structural correlation between the synthetic NiSbSe and the natural kvačekite was confirmed by electron back-scattered diffraction (EBSD). For that purpose, we used a TESCAN Mira 3GMU scanning electron microscope combined with EBSD system Symmetry (Oxford Instruments). The natural sample surface was re-polished with colloidal silica (OP-U) for 20 min to remove an amorphised layer formed due to standard diamond polishing. The EBSD patterns were collected and processed using the proprietary software AZtec 6.0 (Oxford Instruments). The Kikuchi patterns obtained from the natural material (35 measurements on different spots on five grains of natural kvačekite) were found to match the patterns generated from our refined structural model for synthetic NiSbSe (Fig. 8). The Kikuchi bands (11, containing 48 reflectors) were detected using the optimised accuracy indexing mode, using band edges, with a Hough resolution of 65. The values of the mean angular deviation (MAD; goodness-of-fit in the solution) between the calculated and measured Kikuchi bands range between 0.11° and 0.36°. These values hence reveal a very good match; mean angular deviations <1° are considered as indicators of an acceptable fit.
Mineralogy of the associated selenides
A detailed investigation of the selenides associated with kvačekite was undertaken using the same electron microprobe and under the same analytical conditions and correction procedures as previously given for kvačekite.
Berzelianite
Berzelianite belongs among the most abundant selenides within the samples studied. It occurs as anhedral grains up to 1 cm in size in calcite gangue in association with other selenides. In reflected light it is greenish blue and isotropic.
The chemical composition of berzelianite (Table 6) is relatively uniform and is close to the ideal formula Cu2–xSe (with x ≈ 0.18). In addition to the main cation (Cu), minor contents of Ag (up to 0.004 apfu), Co (up to 0.031 apfu) and Ni (up to 0.020 apfu) were found, Co and Ni may be an effect of sub-microscopic inclusions of nickeltyrrellite. In the anion part, the range of SSe–1 substitution is limited to 0.26 apfu S. The mean composition of berzelianite (calculated from 69 point analyses on the basis of 1 Se+S apfu) is (Cu1.81Co0.01)Σ1.82(Se0.89S0.11)Σ1.00.
Nickeltyrrellite–tyrrellite solid solution
Minerals of the nickeltyrrellite–tyrrellite solid solution are abundant in the association studied. They form euhedral grains (metacrysts) up to 300 μm in size (Fig. 9a) and groups of crystals up to 2 mm in length (Fig. 9b), both enclosed in berzelianite. More rarely, they occur as anhedral aggregates up to 800 μm in size in calcite gangue without other selenides (Fig. 9c). In reflected light they are yellowish brown and isotropic.
The chemical composition of the minerals studied corresponds to the formulae Cu(Ni,Co)2Se4–Cu(Co,Ni)2Se4 (Förster et al., Reference Förster, Ma, Grundmann, Bindi and Stanley2019). The determined representation of CoNi–1 substitution (Fig. 10) indicates the presence nickeltyrrellite with Ni contents in the wide range of 1.08–2.43 apfu and on the other hand tyrrellite, where Co contents (1.04–1.29 apfu) only slightly prevail over Ni (0.83–1.29 apfu). A wide extent of CoNi–1 substitution was observed even within individual grains, e.g. from 0.96 to 2.43 apfu Ni in a single grain of 40 μm size. For both minerals, the majority of measured Cu contents (0.55–1.06 apfu) is too low to completely fill the Cu position. It is likely that other cations (Co and/or Ni) have entered this position; a similar situation was described for part of nickeltyrrellite grains from El Dragón by Förster et al. (Reference Förster, Ma, Grundmann, Bindi and Stanley2019). The range of SSe–1 substitution in nickeltyrrellite (0.01–0.25 apfu S) is close to that of tyrrellite (0.02–0.34 apfu). Representative compositions of minerals of the nickeltyrrellite–tyrrellite solid solutions and the corresponding empirical formulae are given in Table 7.
Clausthalite
Clausthalite is a rarer phase, it forms anhedral grains up to 300 μm in size in nickeltyrrellite–tyrrellite and partly replaces earlier berzelianite (Fig. 9d). In reflected light, it is isotropic and has a white colour.
Clausthalite contains, in addition to Pb and Se (Table 6), minor contents of Cu (up to 0.05 apfu) and Co, Ni (both up to 0.02 apfu) and is S-free. Its empirical formula (mean of 23 point analyses on the basis of 2 apfu) is (Pb0.97Cu0.02Ni0.01Co0.01)Σ1.01Se0.99.
Eucairite
Rare eucairite occurs as anhedral grains up to 100 μm in size and partly replaces earlier berzelianite (Fig. 9e), usually in association with clausthalite and nickeltyrrellite–tyrrellite.
Tiny veinlets of gold (with up to 17 wt.% Ag and 5 wt.% Hg) in eucairite grains were observed only occasionally. In reflected light, eucairite is brown and strongly anisotropic in blue green – blue grey rotation tints.
Its chemical composition (Table 6) is close to the ideal formula CuAgSe; only minor contents of Ni, Co (both up to 0.016 apfu) and Cd (up to 0.004 apfu) were determined. The mean composition of eucairite (calculated from 17 point analyses on the basis of 3 apfu) is Ag1.00Cu1.02Co0.01Ni0.01Se0.96.
Hakite-(Zn) and hakite-(Cd)
Hakite forms anhedral grains up to 100 μm in size in berzelianite (Fig. 9f) and nickeltyrrellite–tyrrellite, locally in intergrowths with kvačekite and clausthalite. Hakite replaces earlier berzelianite and clausthalite and can contain relics of these minerals and of nickeltyrrellite. In reflected light it is brownish and isotropic.
According to the nomenclature scheme of minerals of the tetrahedrite group (Biagioni et al., Reference Biagioni, George, Cook, Makovicky, Moëlo, Pasero, Sejkora, Stanley, Welch and Bosi2020) the studied minerals correspond to hakite-(Zn) and two point analyses to hakite-(Cd), both with regular minor Hg contents (Fig. 11). Both minerals were recently defined as new minerals on samples from the Příbram uranium ore district (Czech Republic) by Sejkora et al. (Reference Sejkora, Biagioni, Škácha, Musetti and Dolníček2024). The trigonal position of both studied members from Bukov is predominantly occupied by Cu, with minor amounts of Ag (up to 0.60 apfu). Zinc, and in the case of two point analyses Cd, are dominant Me 2+ elements; there are two trends in minor contents of Hg (Fig. 11) – approximately 0.13 apfu and in the range 0.19–0.39 apfu. The calculated Cu2+ contents reach up to 0.60 apfu. Locally minor contents of Ni and Co were observed, reaching up to 0.19 and 0.16 apfu, respectively. The studied hakites are pure Sb members and the range of SSe–1 substitution is limited to 1.70 apfu S. Representative chemical analyses of hakite-series minerals and the corresponding empirical formulae are given in Table 8.
Conclusions
Kvačekite is a new ternary phase in the Ni–Sb–Se system and a member of the cobaltite group. It does not correspond to any valid or invalid unnamed mineral (Smith and Nickel, Reference Smith and Nickel2007; http://cnmnc.units.it/). Its identity with synthetic cubic NiSbSe (Hulliger, Reference Hulliger1963; Foecker and Jeitschko, Reference Foecker and Jeitschko2001) was confirmed by a study of their chemical composition, reflectance measurements, Raman spectroscopy and electron back-scattered diffraction (EBSD) measurements on the natural sample.
Kvačekite occurs in a Ni–Co–Cu rich mineral association, which also represents the second world occurrence for nickeltyrrellite (Förster et al., Reference Förster, Ma, Grundmann, Bindi and Stanley2019), hakite-(Zn) and hakite-(Cd) (Sejkora et al., Reference Sejkora, Biagioni, Škácha, Musetti and Dolníček2024). This association hosted by a calcite gangue exhibits a limited range of SSe–1 substitution (Fig. 12). The selenide mineralisation studied originated from Ni–Co–Cu-bearing fluids at temperatures above 112°C (absence of umangite, Chakrabarti and Laughlin, Reference Chakrabarti and Laughlin1981) with neutral-to-weakly-alkaline character, high oxygen fugacity and a Se2/S2 fugacity ratio greater than unity.
Acknowledgements
The helpful comments of anonymous reviewer, Panagiotis Voudouris, Peter Leverett, Associate Editor Ian Graham and Principal Editor Stuart Mills are greatly appreciated. The study was financially supported by the Ministry of Culture of the Czech Republic (long-term project DKRVO 2024-2028/1.II.a; National Museum, 00023272) and by the Grant Agency of the Czech Republic (project No 22-26485S).
Supplementary material
The supplementary material for this article can be found at https://doi.org/10.1180/mgm.2024.32.
Competing interests
The authors declare none.