Introduction
The monazite group consists of monoclinic orthophosphates of the APO4 type of light rare-earth elements (LREE: La–Gd), Th, and Ca, crystallising with a P21/n space-group symmetry with nine-fold coordinated A cations. This group includes monazite-(La), monazite-(Ce), monazite-(Nd), monazite-(Sm), monazite-(Gd) and cheralite (Levinson, Reference Levinson1966; Graeser and Schwander, Reference Graeser and Schwander1987; Masau et al., Reference Masau, Černý, Cooper and Chapman2002; Linthout, Reference Linthout2007; Ondrejka et al., Reference Ondrejka, Uher, Ferenc, Milovská, Mikuš, Molnárová, Škoda, Kopáčik and Bačík2023). Monazite-(Ce) is a common accessory mineral in granitic igneous and gneissic metamorphic rocks, and also occurs in detrital sands derived from such crystalline rocks. Monazite-(La) and monazite-(Nd) are much rarer, whereas monazite-(Sm) and especially monazite-(Gd), are extraordinarily rare, growing only in specific magmatic and post-magmatic environments.
Among the rarer varieties of monazite, monazite-(Nd) is more common and its occurrence has been described in gneisses, carbonatites, sandstones, metamorphic magnetite ore deposits, and veins associated with hydrothermal mineralisation (Graeser and Schwander, Reference Graeser and Schwander1987; Del Blanco et al., Reference Del Blanco, Costas Ulbrich, Echeveste and Vlach1998; Kusiak and González-Álvarez, Reference Kusiak and González-Álvarez2006; Pršek et al., Reference Pršek, Ondrejka, Bačík, Budzyń and Uher2010; Števko et al., Reference Števko, Uher, Ondrejka, Ozdín and Bačík2014; Mitchell and Smith, Reference Mitchell and Smith2017; Ondrejka et al., Reference Ondrejka, Uher, Ferenc, Milovská, Mikuš, Molnárová, Škoda, Kopáčik and Bačík2023). Monazite-(Sm) has only been described in a pod of lepidolite-subtype granitic pegmatite from the Annie Claim #3 pegmatitic leucogranite, which is part of the Greer Lake intrusion in southeastern Manitoba (Masau et al., Reference Masau, Černý, Cooper and Chapman2002), and as a phase accompanying monazite-(Gd) in a REE–U–Au hydrothermal quartz vein in the Western Carpathians, Slovakia (Ondrejka et al., Reference Ondrejka, Uher, Ferenc, Milovská, Mikuš, Molnárová, Škoda, Kopáčik and Bačík2023). Monazite-(Nd) occurs mostly as crystals with an average of >12 to up to ~25 wt.% Nd2O3. Two known exceptions include monazite-(Nd) containing 30.32 wt.% Nd2O3 (0.44 Nd atoms per formula unit; apfu) (Graeser and Schwander, Reference Graeser and Schwander1987) and monazite from sandstones belonging to the Grinnell Formation at Red Rock Canyon in the Waterton–Glacier International Peace Park, southwestern Alberta, Canada, mentioned by Kusiak and González-Álvarez (Reference Kusiak and González-Álvarez2006) as containing up to 31 wt.% Nd2O3 (0.45 Nd apfu). The Sm2O3 contents of monazite-(Sm) are known to be relatively low. Masau et al. (Reference Masau, Černý, Cooper and Chapman2002) reported only up to 14.29 wt.% Sm2O3 (0.197 Sm apfu) in the holotype monazite-(Sm) while Ondrejka et al. (Reference Ondrejka, Uher, Ferenc, Milovská, Mikuš, Molnárová, Škoda, Kopáčik and Bačík2023) measured Sm2O3 contents of 17.32 wt.% (0.236 Sm apfu) in monazite-(Sm) and up to 17.68 wt.% Sm2O3 (0.241 Sm apfu) in monazite-(Gd).
Monazite-group minerals and xenotime-(Y) are common accessory minerals in the anatectic pegmatite bodies of the Julianna system at Piława Górna in Lower Silesia, Poland (Szuszkiewicz et al., Reference Szuszkiewicz, Szełęg, Pieczka, Ilnicki, Nejbert, Turniak, Banach, Łodziński, Różniak and Michałowski2013). The monazite occurrences are almost exclusively represented by monazite-(Ce) with varying compositions, however individual crystals of monazite-(Nd) and monazite-(Sm) have been encountered in one vein, described as the Blue Beryl Dyke. The main aim of this contribution was to document their occurrence and clarify the reasons behind their atypically high Nd and Sm enrichment.
Geology and mineralogy of the Blue Beryl Dyke
The Julianna pegmatite system crops out in an active amphibolite–migmatite quarry (50.70327, 16.73677) at Piława Górna in the Góry Sowie Block (Szuszkiewicz et al., Reference Szuszkiewicz, Szełęg, Pieczka, Ilnicki, Nejbert, Turniak, Banach, Łodziński, Różniak and Michałowski2013). The unit constitutes a central part of the Sudetes at the NE margin of the Bohemian Massif in SW Poland. The pegmatites crystallised from anatectic melts produced by partial melting of metasedimentary–metavolcanic rocks during retrograde amphibolite-facies metamorphism ca. 370–380 Ma (Van Breemen et al., Reference Van Breemen, Bowes, Aftalion and Żelaźniewicz1988; Timmermann et al., Reference Timmermann, Parrish, Noble and Kryza2000). Most of the anatectic melts were trapped in neosomes in migmatites, whereas the rest occur in pegmatite bodies located along tectonic zones in the Góry Sowie metamorphic rocks (Szuszkiewicz et al., Reference Szuszkiewicz, Szełęg, Pieczka, Ilnicki, Nejbert, Turniak, Banach, Łodziński, Różniak and Michałowski2013; Pieczka et al., Reference Pieczka, Szuszkiewicz, Szełęg, Janeczek and Nejbert2015). The pegmatite system consists of a network of co-genetic dykes and vein- and nest-type apophyses that vary in size and degree of geochemical evolution and represent barren and rare-element pegmatites. They individually represent various subclasses that correspond to the niobium–yttrium–fluorine (NYF) and lithium–caesium–tantalum (LCT) petrogenetic families in the pegmatite classification of Černỳ and Ercit (Reference Černý and Ercit2005), however overall, the system exhibits hybrid NYF + LCT characteristics. Wise et al. (Reference Wise, Müller and Simmons2022) classified the Góry Sowie pegmatites as an example of type 1 DPA (direct product of anatexis) pegmatites, showing a typical Be–Nb–Ta–P–Li–B geochemical signature.
The Blue Beryl Dyke, which was exposed in 2010, extended over a few tens of metres horizontally with a maximum thickness of ~3–4 m. The zoning was typical up to the blocky unit with interstices among feldspars filled with quartz, but without a massive core. The characteristic features of the dyke were the presence of short prismatic crystals of pale-bluish low-Na,Cs beryl, abundant cassiterite and columbite-group minerals, and subordinate, few millimetre-sized, greenish crystals of gahnite associated with ferro- and zinconigerites-2N1S, (Al,Fe,Zn)2(Al,Sn)6O11(OH) and (Zn,Al,Fe)3(Al,Fe,Ti)8O15(OH), and -6N6S, (Al,Fe,Zn)3(Al,Sn,Fe)8O15(OH) and (Zn,Al)7(Al,Fe3+,Ti,Mg)16O31(OH), as well as genthelvite, Be3Zn4(SiO4)3S (Pieczka et al., Reference Pieczka, Szełęg and Szuszkiewicz2019). Tourmaline-supergroup minerals were represented by Zn-enriched foitite, □(Fe2+2Al)Al6(Si6O18)(BO3)3(OH)3(OH), evolving into schorl, NaFe2+3Al6(Si6O18)(BO3)3(OH)3(OH); Li-bearing tourmalines were unique, occurring only as small dark greenish crystals of Fe-bearing elbaite, Na(Li1.5Al1.5)Al6(Si6O18)(BO3)3(OH)3OH, and fluor-elbaite, Na(Li1.5Al1.5)Al6(Si6O18)(BO3)3(OH)3F, in a few-centimetre-sized muscovite flakes, and as tiny inclusions in quartz and albite (Pieczka et al., Reference Pieczka, Gołębiowska, Jeleń, Włodek, Szełęg and Szuszkiewicz2018). Lithium-bearing mineralisation was also represented by minute aggregates of spodumene, LiAlSi2O6, veinlets of a Cs-dominant dark mica, and single aggregates of partly altered lithiophilite, LiMn2+(PO4), associated with ferrowodginite, Fe2+ Sn4+Ta2O8, and wodginite, Mn2+Sn4+Ta2O8. The REE-bearing phases were represented by rare to accessory monazite-(Ce), monazite-(Nd), monazite-(Sm), xenotime-(Y), allanite-(Ce) and traces of rhabdophane-group phosphates.
Methods
Electron-probe microanalyses of the monazite-group minerals from the Blue Beryl Dyke were performed at the Inter-Institute Analytical Complex for Minerals and Synthetic Substances at the University of Warsaw, Poland, with a CAMECA SX 100 electron microprobe. Wavelength-dispersive X-ray spectrometry was conducted under the following conditions: accelerating voltage of 15 kV, beam current of 70 nA, beam diameter of 2 μm, peak count-time of 20 s, and background times of 10 s before and 10 s after the peak. The following calibration materials (analytical lines, diffracting crystals, and mean detection limits in wt.% element) were used: diopside – Ca (Kα, PET, 0.01) and Si (Kα, TAP, 0.01); YAG – Y (Lα, PET, 0.03); LaP5O14 – La (Lα, PET, 0.02); CeP5O14 – Ce (Lα, PET, 0.03); PrP5O15 – Pr (Lβ, LIF, 0.06); GdP5O14 – Gd (Lα, LIF, 0.03); REE1 – Dy (Lβ, LIF, 0.16) and Er (Lα, LIF, 0.09); REE2 – Eu (Lβ, LIF, 0.15), Ho (Lβ, LIF, 0.16) and Tm (Lα, LIF, 0.13); REE3 – Sm (Lα, LIF, 0.09) and Yb (Lα, LIF, 0.07); REE4 – Nd (Lβ, LIF, 0.09), Tb (Lα, LIF, 0.10) and Lu (Lα, LIF, 0.15); galena – Pb (Mα, PET, 0.04); ThO2 – Th (Mα, PET, 0.04); and UO2 – U (Mβ, PET, 0.04), where REE1–REE4 are silica glass standards containing the respective lanthanides. The raw data were reduced with the PAP routine of Pouchou and Pichoir (Reference Pouchou, Pichoir, Heinrich and Newbury DE1991). Corrections for interferences of REE lines were made on the basis of factors derived for the used microprobe considering the interferences indicated by Reed and Buckley (Reference Reed and Buckley1998). The chemical formulae of monazite were normalised in relation to 4 O apfu. The REE-distribution patterns were normalised to the REE contents in CI chondrite (O'Neill, Reference O'Neill2016). The decomposition of chondrite-normalised REE patterns was performed with the online interactive polynomial fitting using the BlambdaR application (Anenburg and Williams, Reference Anenburg and Williams2022). This application allowed quantitative fitting on the estimated uncertainty for individual ln(REE) in percent (s) and the description of REE patterns through their decomposition into orthogonal polynomial functions (shape components) with coefficients known as the lambda coefficients (λ0, λ1, λ2, λ3, λ4) and the tetrad effect coefficients (τ1, τ2, τ3, τ4), in addition to the calculation of statistically predicted REEλτ concentrations.
Results
The single crystals of monazite-(Nd) (~10 μm in size) and monazite-(Sm) (~30 μm in size) found in the Blue Beryl Dyke were associated with columbite-group minerals that underwent alteration into fluor- to oxy/hydroxy-calciomicrolite and oxy/hydroxy-plumbomicrolite free of fluorine and were accompanied by more common monazite-(Ce). However, monazite-(Ce) was much less abundant than in other dykes of the pegmatite system, forming rare and disseminated minute grains, not more than several tens to 100 μm in size. Monazite-(Nd) occurred as an inclusion in columbite-(Fe) associated with more common inclusions of cassiterite, scheelite, uraninite and BiVO4 (pucherite or clinobisvanite). Monazite-(Sm) formed an intergrowth with columbite-(Mn). Back-scattered electron (BSE) images of the monazites and the respective X-ray element maps depicting Ce, Nd, and Sm distributions in monazite-(Sm) are shown in Fig. 1. The grain of monazite-(Sm) is clearly heterogeneous with a visible correlation between Ce and Nd that were concentrated in different parts of the examined crystal to Sm, although the crystal was Sm-dominated overall. Representative electron-probe microanalyses of the three monazite species are presented in Table 1. The complete set of microprobe analyses of the dyke monazites is available as Supplementary Material (Table S1). The analysed monazite-(Ce) grains show diverse compositions evolving towards more enriched middle rare earth element (MREE, Sm–Dy) contents, including up to 0.272 Nd apfu (19.19 wt.% Nd2O3), 0.152 Sm apfu (11.07 wt.% Sm2O3), and 0.075 Gd apfu (5.61 wt.% Gd2O3). The increasing MREE contents reach their maximum concentration in monazite-(Nd) containing 0.486 Nd apfu (34.12 wt.% Nd2O3), 0.095 Sm apfu (6.90 Sm2O3 wt.%) and 0.029 Gd apfu (2.18 Gd2O3 wt.%), as well as in a more compositionally varied monazite-(Sm), containing up to 0.261 Nd apfu (18.52 wt.% Nd2O3), 0.462 Sm apfu (33.22 wt.% Sm2O3) and 0.148 Gd apfu (11.05 wt.% Gd2O3). All monazite species exhibit significant Eu, Dy and Er concentrations, and monazite-(Ce) is also enriched in Y. The monazite-(Nd) had a composition consisting of (Nd0.486Ce0.257Sm0.095Pr0.077Gd0.0.029La0.026Eu0.019Dy0.002Er0.002U0.001) (P0.999Si0.005)O4 (one spot analysis), while monazite-(Sm) was assigned a composition of (Sm0.462Nd0.191Gd0.148Ce0.113Pr0.028La0.018Eu0.005Dy0.005Er0.005Yb0.002Ca0.012Th0.010)(P0.996Si0.007)O4 in the portion richest in Sm.
Notes: b.d.l. – below detection limit; Σan. (cat.) – total of anions (cations); N = number of spot analysis.
Chondrite-normalised REE patterns (REECN) show continuously decreasing concentrations of La and Ce from monazite-(Ce), through monazite-(Nd), to monazite-(Sm); similar concentrations of Pr, Nd and Eu were observed in monazite-(Ce) and monazite-(Sm), with a distinct enrichment of these elements in monazite-(Nd) and increasing concentrations of Sm and Gd from monazite-(Ce) to monazite-(Sm) (Fig. 2). Monazite-(Sm) and monazite-(Ce) are characterised by a negative Eu anomaly, whereas the REE pattern of monazite-(Nd) reveals no Eu anomaly. The decomposition of the REECN patterns of monazite-(Sm) and monazite-(Ce) was made with the BlambdaR application (Anenburg and Williams, Reference Anenburg and Williams2022) for two populations of λ and τ coefficients: (1) (λ0, λ1, λ2, λ3, λ4) and (τ1, τ2, τ3, 0) and (2) (λ0, λ1, λ2, λ3, 0) and (τ1, τ2, τ3, 0) with the considered Eu anomaly, leading to optimal values of r 2 > 0.99 and χ2 in the range of 1.6–18.7 for s = 0.05 (Table S2). The symbol ‘0’ in place of λ4 and τ4 denotes that the coefficients were excluded from the fitting procedure, e.g. the tetrad effect τ4 was not fitted due to the lack or scarcity of analytical data for elements ranging from Er to Lu. The REECN pattern of the monazite-(Nd) was fitted with the (λ0, λ1, λ2, λ3, λ4) and (0, 0, τ3, 0) coefficients without a considered Eu anomaly, which led to the best statistical parameters of r 2 > 0.998 and χ2 = 2.46 for the assumed s value (Table S2). The negative Eu anomaly, EuCN/Euλτ, in the REECN patterns of monazite-(Sm) varied from 0.03–0.09 [average 0.06(2)], those of monazite-(Ce) from 0.22–0.49 [average 0.33(12)], and in the case of monazite-(Nd), which did not reveal a clear Eu anomaly, the value was estimated at 0.97 on the basis of fitting with the BlambdaR application.
Genetic implications
To the best of our knowledge, the monazite-(Nd) and monazite-(Sm) from Piława Górna are the two most Nd- and Sm-enriched crystals worldwide. This could imply that these elements, with abundances in the Earth's continental crust of ~20 and 3.9 ppm, respectively (Rudnick, Reference Rudnick and White2018), have been concentrated by ~1.5×104 times in monazite-(Nd) and from up to 8×103 times (Nd) to above 70×103 times (Sm) in monazite-(Sm). Thus, it is clear that geochemical fractionation of REE itself within a bulk volume of the pegmatite-forming melt was not the only crucial factor for the chemical composition of these crystals. When explaining the crystal-chemistry of monazite-(Gd), Ondrejka et al. (Reference Ondrejka, Uher, Ferenc, Milovská, Mikuš, Molnárová, Škoda, Kopáčik and Bačík2023) suggested the role of a localised process in a chemically near-closed system leading to the progressive concentration of MREE in fluids as a result of the promoted crystallisation of early monazite-(Ce) and xenotime-(Y), which selectively decreased the LREE and HREE contents in the crystallisation environment. The abnormally high Sm- and Nb-enrichment of the Blue Beryl Dyke monazite-(Sm) and monazite-(Nd), and exceptional rarity of these monazite species in the dyke also suggest local REE geochemical anomalies in the pegmatite-forming melt–fluid system related to variations of its chemical and physical parameters, even at the micro scale. These led to the rapid growth of accessory phases under supercooling conditions which scavenged selected elements leading to their local depletion that was not balanced by diffusion processes. We suggest that the outstripping crystallisation of xenotime-(Y) progressively increased the local abundances of LREE (Gd to La), first enabling the occasional crystallisation of monazite-(Sm) with a large negative Eu anomaly; this then led to the crystallisation of (Sm,Nd)-bearing monazite-(Ce) with greater Eu3+ concentrations than in the monazite-(Sm), and finally to the formation of monazite-(Nd) with a smooth REECN pattern. As is known, the main reason for the Eu anomaly is in the presence of Eu2+ or Eu3+ species in the crystallisation environment (Weill and Drake, Reference Weill and Drake1973; Bau, Reference Bau1991; Deer et al., Reference Deer, Howie and Zussman2001). At higher temperatures and in more reducing magmas, Eu occurs as Eu2+ which is compatible with Ca-bearing melt components, mainly becoming incorporated into plagioclase crystals. At lower temperatures and under more oxidising conditions, the element will occur as Eu3+, becoming incompatible with the melt components but compatible with other trivalent REE. Thus, the decreasing negative Eu anomaly in monazite-(Sm) through monazite-(Ce) to monazite-(Nd) [from 0.06(2) to 0.33(12) and 0.97, respectively] indicates increasing amounts of Eu3+ in probably more oxidising fluids coexisting with the pegmatite-forming melt, which were subsequently remobilised by the coeval local hydrothermal alteration of plagioclase. The fluid also affected the REE fractionation, decreasing the tetrad effect on the third, second, and first tetrads. These observations provide additional insights into the evolution of the crystallisation environment of the monazites.
Monazite-(Sm) exhibiting a negative Eu anomaly has been found as an intergrowth with partly digested columbite-(Mn) undergoing alteration to fluor- and oxy/hydroxycalciomicrolite and finally to F-free plumbomicrolite associated with ferrowodginite, wodginite, F-bearing tourmaline, fluorapatite and Cs-bearing mica. This indicates relatively early crystallisation of the mineral from a pegmatite-forming melt evolving into a residual melt coexisting with aqueous fluids. High-temperature hydrothermal fluids are commonly associated with the presence of strong REE-complexing ligands such as F–, Cl–, OH–, and PO43–, which are reported to be the most effective agents in the mobilisation, transport, fractionation, and deposition of REE species (Gieré, Reference Gieré1990; Jiang et al., Reference Jiang, Wang, Xu and Zhao2005; Migdisov and Williams-Jones, Reference Migdisov and Williams-Jones2014; Williams-Jones, Reference Williams-Jones, Simandl and Neetz2015; Migdisov et al., Reference Migdisov, Williams-Jones, Brugger and Caporuscio2016). Among these anions, Cl– forms the weakest complexes with REE and F– the strongest ones (Migdisov et al., Reference Migdisov, Williams-Jones, Brugger and Caporuscio2016). However, as chloride is a dominant ligand in natural hydrothermal systems, it is mainly responsible for hydrothermal transport of the REE. Thus, Nd and Sm may have undergone additional fractionation within the Blue Beryl Dyke during the hydrothermal stage due to the different stabilities of LREE-chloride complexes, and differences in their mobility and solubility. This fractionation should decrease with increasing atomic number of the REE, which was corroborated in the Blue Beryl Dyke monazites by the decreasing tetrad effect on the third, second, and first tetrads from monazite-(Sm) through monazite-(Ce) to monazite-(Nd). Due to significant differences between the stability/mobility of LREE-chloride complexes and those involving HREE, this process may have led to Sm- and Nd-enrichment resulting in the local saturation in SmPO4 and NdPO4 at different stages of the system's evolution. The destabilisation of Nd–Cl and Sm–Cl complexes could also have been related to the neutralisation of the solution's acidity, e.g. through F consumption during the formation of F-bearing tourmaline, fluorapatite, or fluor-calciomicrolite.
In melt–fluid systems, Th and U, the two elements that are commonly associated with REE in monazite-group minerals, are highly compatible. Among the studied monazite crystals, monazite-(Ce) was the most enriched in Th and U, containing up to 14.91 wt.% ThO2 and 1.76 wt.% UO2. The amounts of ThO2 in monazite-(Sm) were <2 wt.% and UO2 was absent altogether; in turn, small amounts of UO2 were measured in monazite-(Nd) (0.15 wt.%), whereas ThO2 was absent. According to Keppler and Wyllie (Reference Keppler and Wyllie1990, Reference Keppler and Wyllie1991), the Th and U partitioning coefficients between fluid and melt are very low if water is the only volatile component present. However, this partitioning increases significantly as the F– concentration increases, indicating the extraction of certain amounts of Th and U from the melt through the formation of fluoride complexes in the fluid. Chloride anions may form complexes with U, but not with Th; the latter will not be extracted under such conditions. Consequently, fluids with higher Cl– concentrations coexisting with more peraluminous silicate melts will reduce the initial U contents of the melts as the latter will partition into the fluid phase. These data on the behaviours of Th and U in melt–fluid systems corroborate the crystallisation of monazite-(Sm) from a melt which became progressively enriched in volatiles, including F and Cl. The system evolved progressively into a hydrated F-bearing residual melt from which the (Sm,Nd)-bearing monazite-(Ce) crystallised and which lost U through the formation of (U–Cl)-complexes in the fluid. Monazite-(Nd), occurring as ~10 μm inclusions in columbite-(Fe), had a different REECN pattern shaped by λ coefficients that differed from those characteristic of monazite-(Sm) and monazite-(Ce); this monazite also had the weakest tetrad effect only on the third tetrad, and a negligible Eu anomaly. The small UO2 contents and the absence of ThO2 as well as its coexistence with inclusions of other phases crystallised from fluids (cassiterite, scheelite, uraninite and a BiVO4 phase) indicate late crystallisation from fluids devoid of fluorine, which potentially could be coeval with the formation of F-free oxy/hydroxy-plumbomicrolite, as well as possible formation through secondary remobilisation and in situ crystallisation. These data agree with the observations of Gysi et al. (Reference Gysi, Harlov and Miron2018) who showed that the REECN profiles of monazites could be classed into two groups based on their LREE contents: one group showing the limited fractionation of La, Ce, Pr and Nd, and the other group showing a strong geochemical fractionation and including enhanced concentrations of Sm, Eu and Gd. These authors explained the differentiation as resulting from the incorporation of REE with varying ionic radii to match the geometric parameters of the monazite structure.
Conclusions
The presence of highly fractionated accessory mineral phases such as monazite-(Sm) and monazite-(Nd) in the pegmatites described here resulted from the interaction of several factors. The most important of these appear to be:
The presence of highly fractionated pegmatite melts and hydrothermal fluids during the final stages of pegmatite solidification, which were enriched in REE and other metals in the form of complexes mainly with Cl–, F–, and OH– ligands.
A very rapid growth of all accessory phases under supercooling conditions which scavenged selected elements leading to their local depletion in the melt. On the scale of the bulk pegmatite body, this process strongly differentiated the local concentrations of these elements, which were not subsequently balanced by diffusion processes.
The crystallisation of pegmatite-forming melts under supercooling conditions forced the growth of accessory phases that consumed elements that were available in the local environment, and which could occupy acceptable structural positions within the structures of the crystallising phases.
Acknowledgements
We thank Monika A. Kusiak and Martin Ondrejka for helpful and constructive comments and suggestions, which greatly improved the manuscript. We are also very grateful to Sam Broom-Fendley for language corrections. This study was supported by AGH University of Science and Technology grant 16.16.140.315.
Supplementary material
The supplementary material for this article can be found at https://doi.org/10.1180/mgm.2023.38.
Competing interests
The authors declare none.