Fluoralforsite, Ba5(PO4)3F – a new apatite-group mineral from the Hatrurim Basin, Negev Desert, Israel | Mineralogical Magazine | Cambridge Core (original) (raw)

Introduction

The apatite-supergroup minerals can be described by the general formula IX_M_12VII_M_23(IV_T_O4)3_X_ (Z = 2), where the Roman numbers indicate the ideal coordination number; M and T are cationic sites and X is the anionic site. They crystallise in the hexagonal or pseudo-hexagonal system. The structure consists of columns parallel to the c axis made of face-sharing _M_1-centred polyhedra. The coordination number of _M_1 can be divided into six closer and three more distant ligands (oxygens). The closer ones form a metaprism, whereas taking into account the more distant oxygens leads to a tricapped trigonal prism (Henderson et al., Reference Henderson, Bell, Charnock, Knight, Wendlandt, Plant and Harrison2009; Pasero et al., Reference Pasero, Kampf, Ferraris, Pekov, Rakovan and White2010). The _M_1 columns are connected by corner-sharing tetrahedra (T sites). The six columns of tricapped trigonal prisms combined with tetrahedra form a channel with hexagonal symmetry along the _c_-axis (Henderson et al., Reference Henderson, Bell, Charnock, Knight, Wendlandt, Plant and Harrison2009). The inter-column space (channels) is occupied by _M_2 cations and anionic sites (X = Cl, F and OH). The coordination number of _M_2 sites varies depending on the type of cations and anions at the X site. In the archetype structure (space group _P_63/m), _M_2 is seven-fold coordinated – a distorted pentagonal bipyramid (Pasero et al., Reference Pasero, Kampf, Ferraris, Pekov, Rakovan and White2010). Anionic sites lie in the centre of the channel. However, their position along the _c_-axis depends on the type of adjacent anions and cations (Hughes et al., Reference Hughes, Cameron and Crowley1989, Reference Hughes, Cameron and Crowley1990).

The high flexibility of the chemical composition characterises the apatite structure-type. Currently, the following elements are mentioned as the species-defining ones: M = Ca2+, Sr2+, Ba2+, Pb2+, Mn2+, Cd2+, Na+, Y3+, La3+, Ce3+, Nd3+and Bi3+; T = P5+, As5+, V5+, Si4+, S6+ and B3+; X = F–, (OH)–, Cl– and O2– (Pasero et al., Reference Pasero, Kampf, Ferraris, Pekov, Rakovan and White2010). Based on chemical and crystal-chemical grounds, the following five mineral groups have been distinguished within the apatite supergroup: apatite group, hedyphane group, belovite group, britholite group and ellestadite group. Among all these groups, the minerals belonging to the apatite group are the most widespread in Nature. This group includes minerals with pentavalent elements at the T site (P, V and As) and containing the same prevailing (species-defining) bivalent cation at the _M_1 and _M_2 sites. Therefore, the simplified formula, M 5(T_O4)3_X, is used in literature for this group (Pasero et al., Reference Pasero, Kampf, Ferraris, Pekov, Rakovan and White2010). Most of the apatite-group members crystallise in space group _P_63/m. In contrast, the hedyphane group can host different cations at the _M_1 and _M_2 sites, resulting in various hexagonal and pseudo-hexagonal symmetries. On the other hand, belovite-type minerals crystallise in lower symmetry due to the ordered arrangement of cations at two independent _M_1 sites, _M_1 and _M_1’. _P_63 symmetry is characteristic of minerals containing divalent cations at all M sites, and P bar3\bar{3}bar3 is characteristic symmetry of minerals containing REE3+ and monovalent cations at split _M_1 sites and a divalent cation at the _M_2 site (Pasero et al., Reference Pasero, Kampf, Ferraris, Pekov, Rakovan and White2010; Pekov et al., Reference Pekov, Zubkova, Husdal, Kononkova, Agakhanov, Zadov and Pushcharovsky2012).

Fluoralforsite, ideally Ba5(PO4)3F is the Ba-analogue of fluorapatite and the F-analogue of alforsite Ba5(PO4)3Cl (Newberry et al., Reference Newberry, Essene and Peacor1981) belonging to the apatite group. It was found in coarse-grained rankinite paralava in the Hatrurim Basin (Negev Desert), Hatrurim Complex in Israel (Fig. 1), which is a well-known source of new minerals, especially Ba species such as zadovite, aradite, gurimite, hexacelsian, bennesherite and mazorite (Galuskin et al., Reference Galuskin, Gfeller, Armbruster, Galuskina, Vapnik, Murashko, Włodyka and Dzierżanowski2015a; Galuskina et al., Reference Galuskina, Galuskin, Vapnik, Prusik, Stasiak, Dzierżanowski and Murashko2017a; Juroszek et al., Reference Juroszek, Galuskina, Krüger, Krüger, Vapnik and Galuskin2022; Krzątała et al., Reference Krzątała, Krüger, Galuskina, Vapnik and Galuskin2022). It is the fourth barium phosphate belonging to the apatite supergroup, after alforsite Ba5(PO4)3Cl (apatite group), kuannersuite-(Ce) NaCeBa3(PO4)3F0.5Cl0.5, (belovite group) and miyahisaite (Sr,Ca)2Ba3(PO4)3F (hedyphane group; Table 1). The synthetic analogue of fluoralforsite was reported by Mathew et al. (Reference Mathew, Mayer, Dickens and Schroeder1979). A comparison of the main crystallographic properties of barium apatite-supergroup minerals and synthetic analogues is presented in Table 1.

Figure 1. Samples of rankinite paralava from (a) wadi Zohar tributary and (b), Gurim Anticline Hatrurim Basin, in which fluoralforsite was found. The areas framed by red rectangles in (a) and (b) are magnified in Figs 3a,b and 3e,f, respectively. The main minerals of this paralava are distinguished by colour: Ti-rich andradite is brown to black, melilite of the gehlenite–alumoåkermanite series is brown to yellow-light brown, rankinite, fluorapatite, kalsilite and cuspidine are pale pink, light-grey or colourless, secondary Ca-hydrosilicates are chalk white or beige.

Table 1. Comparison of the structural and physical properties of fluoralforsite with selected minerals of the apatite supergroup and synthetic Ba5(PO4)3F.

Fluoralforsite, Ba5(PO4)3F, was named according to the Commission on New Minerals, Nomenclature and Classification of the International Mineralogical Association (IMA-CNMNC) recommendations for the new minerals of the apatite supergroup (Pasero et al., Reference Pasero, Kampf, Ferraris, Pekov, Rakovan and White2010; Hatert et al., Reference Hatert, Mills, Pasero and Williams2013). The mineral was named by adding the adjectival prefix for the anion, fluor, to the root name, alforsite, giving fluoralforsite. Consequently, the mineral symbol is Fafr.

A small polished fragment of paralava from the tributary of wadi Zohar with grains of fluoralforsite has been deposited in the mineralogical collection of the Natural History Museum of Bern, Bernstrasse 15, CH-3005 Bern, Switzerland, catalogue number 44951. This work provides data about the new mineral and discusses the crystallographic aspects and ordering of cations and anions in detail.

Occurrence and description of fluoralforsite

Fluoralforsite Ba5(PO4)3F was found in coarse-grained rankinite paralava in the Hatrurim Basin (Negev Desert), Hatrurim Complex in Israel (Fig. 1). The unique Hatrurim Complex is composed of pyrometamorphic rocks and their alteration products, which are distributed on both sides of the Dead Sea Transform Fault in Israel, the West Bank, and Jordan (Fig. 2; Gross, Reference Gross1977; Burg et al., Reference Burg, Starinsky, Bartov and Kolodny1991; Geller et al., Reference Geller, Burg, Halicz and Kolodny2012). The pyrometamorphic rocks were formed as a result of the high-temperature metamorphism of a sedimentary protolith at sanidinite facies (>600°C and low pressure). They are represented mainly by spurrite marble, larnite pseudoconglomerates, gehlenite hornfelses and different types of paralavas (Vapnik et al., Reference Vapnik, Sharygin, Sokol and Shagam2007; Novikov et al., Reference Novikov, Vapnik and Safonova2013; Galuskina et al., Reference Galuskina, Vapnik, Lazic, Armbruster, Murashko and Galuskin2014). In the Hatrurim Basin, the common type of paralava occurs as lenses within the hornfelses mainly composed of gehlenite, andradite and larnite. The lenses range in size from about ten centimetres to a few millimetres. The occurrence of the paralava in hornfels suggests that the high temperature of the pyrometamorphic processes caused partial protolith melting on a local scale. Moreover, rock-forming minerals of coarse-grained paralava contain melt inclusions with 1200–1250°C homogenisation temperature (Sharygin et al., Reference Sharygin, Vapnik, Sokol, Kamenetsky and Shagam2006). Andradite–flamite, gehlenite–flamite and walstromite–kalsilite intergrowths are present in the paralava, which also indicate crystallisation of rocks from melt (Gfeller et al., Reference Gfeller, Widmer, Krüger, Galuskin, Galuskina and Armbruster2015; Krzątała et al., Reference Krzątała, Krüger, Galuskina, Vapnik and Galuskin2020). The characteristic feature of rankinite paralavas is the occurrence of the Ba, V-bearing enclaves, which are a source of new minerals species such as zadovite, BaCa6[(SiO4)(PO4)](PO4)2F, aradite, BaCa6[(SiO4)(VO4)](VO4)2F, hexacelsian, BaAl2Si2O8, gurimite Ba3(VO4)2, bennesherite Ba2FeSi2O7, pliniusite Ca5(VO4)3F and mazorite Ba3(PO4)2.

Figure 2. (a) Schematic map of the Middle East with outcrops of pyrometamorphic rocks of the Hatrurim Complex (dark yellow). The framed part is magnified in b. (b) Schematic map of the Hatrurim Basin. The sampling areas of rankinite paralava are marked by a red pin (an unnamed tributary of wadi Zohar) and a blue pin (Gurim Anticline) (modified after Britvin et al., Reference Britvin, Murashko, Vapnik, Polekhovsky and Krivovichev2015).

Fluoralforsite has been found in the rankinite-rich paralava (the holotype locality 31°11.3'N, 35°16.6'E) from the upper reaches of the tributary of the wadi Zohar. Additionally, it was found in a paralava sample from Gurim Anticline (Fig. 2). It occurs in small enclaves in association with the following Ba-minerals: walstromite, zadovite–aradite series, gurimite, bennesherite, barioferrite and baryte (Fig. 3). The rock-forming minerals of those paralava are mainly represented by rankinite, Ti-rich andradite, members of the gehlenite–alumoåkermanite series, members of the fluorapatite–fluorellestadite series, kalsilite and wollastonite (Fig. 3). Additionally, the following minerals occur in the association: pliniusite, fresnoite, celsian, members of the magnetite–magnesioferrite–trevorite series, U–Nb-bearing cuspidine, nepheline, native copper and hashemite. The paralava voids are filled by low-temperature mineralisation represented by zeolites, afwillite and tobermorite-like Ca-hydrosilicates.

Figure 3. Back-scattered electron (BSE) images of fluoralforsite in paralava near the wadi Zohar tributary (a–d) and Gurim Anticline (e,f). Crystals of fluoralforsite occur between rock-forming minerals in small enclaves enriched in barium. The white frames correspond to the magnified fragments shown on the right side. Mineral abbreviations after Warr (Reference Warr2021): Aåk – alumoåkermanite, Adr – andradite, Baf – barioferrite, Bnh – bennesherite, Brt – baryte, Csp – cuspidine, Fafr – fluoralforsite, Fap – fluorapatite, Fmt – flamite, Gh – gehlenite, Gur – gurimite, HSi-Ca – hydrosilicate of calcium, Kls – kalsilite, Mag – magnetite, Mzo – mazorite, Rnk – rankinite, Wo – wollastonite, Ws – walstromite, Zad – zadovite.

The crystals of fluoralforsite are colourless and transparent up to 50 μm in length and characterised by a white streak and a vitreous lustre. The cleavage is not observed. The tenacity is brittle and the fracture is uneven. Fluoralforsite is uniaxial (–) with refractive indices (589 nm) n ω = 1.689(3) and n ɛ = 1.687(3). Based on the average chemical composition and unit cell parameters, the calculated density is 4.57 g/cm3. Micro-hardness measurements gave a mean VHN25 value of 320.4 kg/mm2 (ranging from 301 to 340 kg/mm2) based on 12 measurements. This value corresponds to a Mohs hardness of ca. 4–4½.

Experimental methods

The morphology and composition of the minerals were studied using optical microscopy, a scanning electron microscope (Phenom XL and FEI Quanta 250, Institute of Earth Science, Faculty of Natural Sciences, University of Silesia in Katowice, Poland) and an electron microprobe analyser (EMPA, Cameca SX100, Institute of Geochemistry, Mineralogy and Petrology, University of Warsaw). Chemical analyses were carried out (wavelength dispersive spectroscopy mode, 15 kV, 10 nA and ~1 μm beam diameter) using the following lines and standard materials: F_K_α – synthetic fluorophlogopite; Cl_K_α – sodalite; Na_K_α – albite; Si_K_α and Ca_K_α – diopside; K_K_α – orthoclase; Fe_K_α – Fe2O3; Ba_L_α and S_K_α– baryte; Sr_L_α – celestine; V_K_α – V2O5; and P_K_α – YPO4.

The Raman spectra of fluoralforsite were recorded on a WITec alpha 300R Confocal Raman Microscope (Institute of Earth Science, Faculty of Natural Sciences, University of Silesia in Katowice, Poland) equipped with an air-cooled solid-state laser (488 nm) and a CCD camera operating at –61°C. The laser radiation was coupled with a microscope through a single-mode optical fibre with a diameter of 3.5 μm. An air Zeiss LD EC Epiplan-Neofluan DIC–100/0.75NA objective was used. Raman scattered light was focused on a broadband single-mode fibre with an effective pinhole size of ~30 μm and a monochromator with a 600 mm–1 grating. The power of the laser at the sample position was ~40 mW. An integration time of 5 s, an accumulation of 20 scans and a resolution of 3 cm–1 were chosen. The monochromator was calibrated using the Raman scattering line of a silicon plate (520.7 cm–1). Spectra processing, such as baseline correction and smoothing, was performed using the SpectraCalc software package GRAMS (Galactic Industries Corporation, NH, USA). Band fitting was performed using a Gauss–Lorentz cross-product function, with a minimum number of component bands used for the fitting process.

Single-crystal X-ray diffraction data were collected with ω scans at different phi-kappa-theta settings using a four-circle Rigaku Synergy R diffractometer (University of Bern, Switzerland), equipped with a PhotonJet-R source, featuring a Mo-rotating anode (Mo_K_α = 0.71073 Å) and a Hypix Arc 100° detector. The measurement was carried out at room temperature. The investigated sample consisted of a grain (~0.05 × 0.04 × 0.03 mm) from the wadi Zohar tributary. The analysis of the collected reflections showed that the selected grain consisted of several phases. Approximately 40% of the total reflections were indexed with a hexagonal unit cell: a = 10.00 Å, c =7.54 Å; V = 654 Å3, corresponding to fluoralforsite. The second set of reflections (~30%) was indexed by the rhombohedral cell a = 7.083 Å, c = 25.498 Å; and V =1107.5 Å3, ascribed to the mineral zadovite (BaCa6[(SiO4)(PO4)](PO4)2F; R bar3\bar{3}bar3m; a = 7.097 Å and c = 25.73 Å) (Galuskin et al. Reference Galuskin, Gfeller, Galuskina, Pakhomova, Armbruster, Vapnik, Włodyka, Dzierżanowski and Murashko2015b). No specific phase could be identified in the remaining diffraction spots. Experimental details are summarised in Table 2. The data processing was performed using the software package CrysalisPro 42.89a. The data reduction was carried out by using an empirical absorption correction. The fluoralforsite structure was solved in the _P_63/m space group by direct methods using the software SHELXS (Sheldrick, Reference Sheldrick2008) and then refined using SHELXL (Sheldrick, Reference Sheldrick2015) in the WinGX package (Farrugia, Reference Farrugia2012). The crystallographic information file has been deposited with the Principal Editor of Mineralogical Magazine and is available as Supplementary material (see below). Nevertheless, the integration, solution, and refinement attempts were carried out for test purposes in the following space groups _P_63, P bar3\bar{3}bar3, and _P_21/m (Supplementary Tables S1–S6). All of them converged to inferior to _P_63/m refinement. For better comparison with literature data, the final structural model was refined by using the atomic coordinates reported by Hata et al. (Reference Hata, Marumo, Iwai and Aoki1979) for alforsite (_P_63/m).

Table 2. X-ray data collection and crystal structure refinement parameters of fluoralforsite.

Chemical composition

The empirical formulas of fluoralforsite were calculated on the basis of 8 cations (Table 3). The following crystal chemical formulas have been obtained on the crystals from the wadi Zohar tributary and Gurim Anticline locality, respectively: $ ( {{\rm B}{\rm a}_{3.81}{\rm C}{\rm a}_{0.97}}$ $ {{\rm N}{\rm a}_{0.07} {\rm K}_{0.05}{\rm S}{\rm r}_{0.05}{\rm F}{\rm e}_{0.05}} ) _{\Sigma 5} ( {{\rm P}^{5 + }_{2.32} {\rm V}^{5 + }_{0.29} {\rm S}^{6 + }_{0.22} {\rm S}{\rm i}^{4 + }_{0.17} } ) _{\Sigma 3}{\rm O}_{12}( {{\rm F}_{0.85}{\rm C}{\rm l}_{0.13}} ) _{\Sigma 0.98} $ $ ( {{\rm B}{\rm a}_{3.92}{\rm C}{\rm a}_{0.88}{\rm N}{\rm a}_{0.08}{\rm K}_{0.02}{\rm S}{\rm r}_{0.11} ) _{\Sigma 5.01}} ( {{\rm P}^{5 + }_{2.17} {\rm V}^{5 + }_{0.41} {\rm S}^{6 + }_{0.23} {\rm S}{\rm i}^{4 + }_{0.11}{\rm A}{\rm l}_{0.07}} ) _{\Sigma 2.99} ![](/binary/version/id/urn:cambridge.org:id:binary:20240209120312010−0659:S0026461X23000580:S0026461X23000580inline1c.png) ![](/binary/version/id/urn:cambridge.org:id:binary:20240209120312010−0659:S0026461X23000580:S0026461X23000580inline1c.png) {{\rm O}_{12.02}( {{\rm F}_{0.92}{\rm C}{\rm l}_{0.02}} ) _{\Sigma 0.94.}}$ The analyses of fluoralforsite crystals reveal a significant substitution of Ca2+, nearly equal to one apfu. Thus, the simplified and ideal formula can be written as (Ba,Ca)2Ba3(PO4)3F and Ba5(PO4)3F, respectively.

Table 3. Chemical data (in wt. %) of fluoralforsite from wadi Zohar tributary and Gurim Anticline*.

Raman investigation

The Raman spectra of phosphate minerals of the apatite group are characterised by strong ν1 symmetric stretching vibrations of (PO4)3− tetrahedra at 950–990 cm–1. The weak bands of ν2 symmetric bending vibrations, ν3 antisymmetric stretching vibrations and ν4 antisymmetric bending vibration are observed in 400–500 cm–1, 1000–1150 cm–1 and 550 and 660 cm–1 regions, respectively. Bands below 300 cm–1 are related to the lattice vibrations (Litasov and Podgornykh, Reference Litasov and Podgornykh2017). In the Raman spectra of fluoralforsite (Fig. 4), the main band of (PO4)3− groups vibrations occur at ~942/941 cm–1. The bands at ~410–430 cm–1 are related to ν2 symmetric bending vibrations. The weak bands at ~1008–1048 cm–1 can be assigned to ν3(PO4)3– antisymmetric stretching vibrations, whereby ν4 antisymmetric bending vibration occurs at ~550–590 cm–1 (Levitt and Condrate, Reference Levitt and Condrate1970; Meegoda et al., Reference Meegoda, Bonner, Loutts, Stefanos and Miller1999). The bands at 840/830 cm−1 can be attributed to asymmetric stretching vibrations (ν3) of VO43– or symmetric stretching vibration (ν1) of SiO44– (Schulte et al., Reference Schulte, Buchter, Chai, Scheps and Kokta1995; Comodi et al., Reference Comodi, Liu, Stoppa and Woolley1999; Frost et al., Reference Frost, Palmer, Čejka, Sejkora, Plášil, Bahfenne and Keeffe2011; Galuskin et al., Reference Galuskin, Gfeller, Galuskina, Pakhomova, Armbruster, Vapnik, Włodyka, Dzierżanowski and Murashko2015b). Low-intensity bands at ~995/991 cm–1 can be assigned to the vibrations of (SO4)2− groups. The bands from moderate to weak intensity in the 80–250 cm–1 region can be connected with lattice or Ba(Ca)–O modes. The remaining bands in the Raman spectrum of fluoralforsite from both localities and their assignments are [cm−1]: 349, 365/350 (ν2VO43–), 379/385 (ν4VO43–), 625/626 (ν4SO42–), 840/830 (ν3VO43– or ν1SiO44–), 860/857 (ν1VO43–), 995/991 (ν1SO42–) (Fig. 4). Accurate band assignment in the Raman spectra of fluoralforsite is hampered by overlapping bands corresponding to vanadate, sulfate and silicate anionic-group vibrations. The spectra analysis from 3600 cm–1 to 3700 cm–1 indicates the absence or very low content of hydroxyl groups in the fluoralforsite channel.

Figure 4. The Raman spectrum of fluoralforsite from (a) the type locality and (b) Gurim Anticline Hatrurim Basin, Hatrurim Complex, Israel.

The structure of fluoralforsite was solved and refined in space group _P_63/m. During the first refinement cycles, _M_1, _M_2 and tetrahedral (_P_1) sites were located (Table 4). Initially, _M_1 and _M_2 sites and tetrahedral position (P_1) were refined using barium and phosphorus scattering curves, respectively. The next refinement cycle allowed for the location of the anion position (F site, Wyckoff position 4_e), which was refined with fluorine scattering factors. However, a strong peak in the difference-Fourier maps was observed along the c_-axis at 1.0 Å from the F site. This residual electron density was modelled as a low-occupied site (2_b position), corresponding to the Cl position in alforsite (Table 4; Hata et al. Reference Hata, Marumo, Iwai and Aoki1979; Newberry et al., Reference Newberry, Essene and Peacor1981). In order to model the residual electron density close to the F site, a new low-occupied site F1 (at ~0.6 Å from F) was inserted in the model. Similarly, the residual peak close to M2, was modelled by an additional site, _M_2A, refined with barium scattering factors (Tables 4, 5). A calculation of the mean electron number of the _M_1 site gave values of 38.56 e –, indicating the significant substitution of a lighter cation. Owing to the complex chemistry of the crystal investigated, _M_1 and _P_1 sites were described with mixed scattering curves according to the chemical microanalyses. _M_1 occupancy was refined with Ba vs. Ca curves and finally converged to 0.501(3) Ba and 0.498(3) Ca. The tetrahedral position was refined using a constrained mixed-scattering curve (2.8 P + 0.2 V). The refined site-scattering values are compatible with the chemical formula calculated from the EMPA data (Tables 2, 3). The location of Ca at the _M_1 position is in agreement with the suggestion of Pasero and co-authors guidelines for the apatite supergroup, according to which smaller cations fill the M_1 site and larger cations the M_2 site (Pasero et al., Reference Pasero, Kampf, Ferraris, Pekov, Rakovan and White2010). All sites, except low-occupied Cl, F, F1 and M_2A, were refined using anisotropic displacement parameters (Table 5). The pronounced displacement of O3 was modelled by splitting the site over O3 and O3A (12_i positions) positions, 0.5 Å apart (Table 4, Table 5). The O2 site was located at the 12_i position, whereas O1 was at the 6_h special position. The final _R_-value converged to 1.92% for 1009 unique reflections and 64 refined parameters.

Table 4. Atom coordinates, equivalent isotropic or isotropic (*) displacement parameters and site occupancies of fluoralforsite.

Table 5. Anisotropic displacement parameters (Å2).

Structure description

Fluoralforsite has an apatite structure type and crystallises in the space group P_63/m_ with the following unit cell parameters: a =10.0031(2), c = 7.5382(2) Å and V = 653.23(3) Å3 (Table 2). The _M_1 site, building columns along the _c_-axis, is occupied equally by Ba and Ca (Figs 5, 6). The _M_1–O bond distances vary from 2.60Å (_M_1–O1) to 2.90Å (_M_1–O3). It is worth noting that two different _M_1–O2 bond distances (2.40Å and 2.79 Å) result as a consequence of the O2 position, which is shifted with respect to the mirror plane. However, the two equivalent O2 positions (at z and – z + ½) cannot be simultaneously present in the coordination sphere of _M_1 (Table 6, Fig. 6). Hence, the nine-fold coordination is maintained, similar to the archetype structure (Table 6). As mentioned above, the O2 site lies slightly above the mirror plane (z = 0.21, Table 4), resulting in two possibilities of tetrahedra tilting (Fig. 5b). The tetrahedron connecting _M_1 columns, is occupied mainly by phosphorus with slight impurities of vanadium. The average <_T_–O> bond distance is 1.55Å, as for the pure phosphorus tetrahedron (Mathew et al., Reference Mathew, Mayer, Dickens and Schroeder1979; Newberry et al., Reference Newberry, Essene and Peacor1981). This value also agrees with the estimated bond distance (1.56 Å) considering the V, S and Si substitutions at the _P_1 site, as shown by chemical analyses (Nyfeler and Armbruster, Reference Nyfeler and Armbruster1998; Hawthorne et al., Reference Hawthorne, Krivovichev, Burns, Alpers, Jambor and Nordstrom2000; Baikie et al., Reference Baikie, Schreyer, Wei, Herrin, Ferraris, Brink, Topolska, Piltz, Price and White2014). The cation position in the inter-column space is partially disordered along the _c_-axis and split over two barium sites, _M_2 (Occupancy 0.93) and _M_2A (Occ. 0.037) (Table 4; Fig. 7). The channel contains three partially occupied anionic sites (Fig. 5c). Fluorine occupies two 4e positions, F and F1, respectively. In turn, the 2b position is a low-occupied (Occ. = 0.112(10)) Cl site corresponding to the Cl position in the alforsite (Fig. 7; Hata et al., Reference Hata, Marumo, Iwai and Aoki1979; Newberry et al., Reference Newberry, Essene and Peacor1981). The anionic arrangement is described in detail in the discussion session.

Figure 5. (a) Arrangement of _M_1 and P polyhedra in the crystal structure of fluoralforsite projected along [110]. _M_2 and anionic sites are not shown for the sake of clarity. The circle marks the part magnified in b; (b) phosphorus tetrahedron with split O2 and O3 sites; (c) crystal structure of fluoralforsite projected along [001]. The triangle highlights the part of the structure magnified in Fig. 7. Turquoise and light violet polyhedra represent the _M_1 prism and P tetrahedra, respectively. _M_2 sites are shown as a turquoise sphere, red spheres correspond to oxygens and dark blue and green spheres represent fluorine and chlorine, respectively. The dotted line indicates the unit cell. Drawn using _CrystalMaker_® software.

Figure 6. (a) Two variants of an _M_1 metaprism coordinated by O1 and O2, displaced off the mirror plane; (b) ordered _M_1 column arrangement along [001]. The radii proportions of Ca/Ba have been preserved. Drawn using _CrystalMaker_® software.

Table 6. Selected interatomic distances (Å).

Figure 7. Anion positions (_c_-axis vertical) with respect to the mirror plane in fluoralforsite. Atoms are represented as thermal ellipsoids. Colour code as in Fig. 5. The dark grey triangle represents the mirror plane. The dotted line indicates _c_-axis.

Discussion

The a unit-cell parameter of fluoralforsite varies considerably compared to that of alforsite, though only slightly from the synthetic phase Ba5(PO4)3F (Table 1; Mathew et al., Reference Mathew, Mayer, Dickens and Schroeder1979; Newberry et al., Reference Newberry, Essene and Peacor1981). The difference in the a parameter between fluoralforsite and alforsite may be ascribed to the different radii of the cations (i.e. Ba and Ca) at the _M_1 site. Similarly, the average bond distances of the tricapped trigonal prism for the investigated fluoralforsite (<_M_1–O> = 2.59 Å; Table 6) are significantly shorter than those in alforsite (<_M_1–O> = 2.83 Å), as a result of Ca2+ substitution at the _M_1 position. Due to the displacement of the O2 position off the mirror plane, two variants of _M_1 coordination are possible (Fig. 6a). In the first version, the _M_1–O2 bond distance is equal to 2.79 Å. In that case, the _M_1 position is preferentially occupied by Ba2+. In the second variant, the _M_1–O2 distance is equal to 2.40 Å and the cationic position is occupied by Ca2+. As a result, an ordered arrangement within individual columns should take place (Fig. 6b). Assuming the cations are ordering at _M_1, the symmetry would lower to P_63, similar to fluorcaphite and fluorstrophite (belovite group) or even to P bar3\bar{3}bar3 such as in kuannersuite-(Ce) and hydroxylhedyphane (Table 1; Friis et al., Reference Friis, Balić-Žunić, Pekov and Petersen2004; Pekov et al., Reference Pekov, Britvin, Zubkova, Pushcharovsky, Pasero and Merlino2010; Biagioni et al., Reference Biagioni, Hålenius, Pasero, Karlsson and Bosi2019). Monoclinic symmetry is also probable, such as in pure synthetic hydroxyl- and chlorapatite, where anionic column ordering results in a double b unit cell parameter (Hughes and Rakovan, Reference Hughes, Rakovan, Kohn, Rakovan and Hughes2002; Pasero et al., Reference Pasero, Kampf, Ferraris, Pekov, Rakovan and White2010). However, in fluoralforsite the P_63/m symmetry is maintained due to different ordering of Ca/Ba among adjacent columns. The mirror plane at ¼ and ¾ z would be lost if there were likewise ordering in adjacent columns. The possibility of different space groups (_P_63, P bar3\bar{3}bar3, and _P_21/m) were tested (Supplementary Tables S1–S6); however, no indication of ordering at the two symmetry-independent _M_1 and _M_1’ sites was observed.

In fluoralforsite, tilting of _M_1 polyhedra is caused by the displacement of the O1–O2 edge off the mirror plane (Fig. 5b), whereas in kuannersuite-(Ce) (P bar3\bar{3}bar3, belovite group) and hydroxylhedyphane (P bar3\bar{3}bar3, hedyphane group) a displacement off the polyhedra edge parallel to c was observed. We consider that the tilting is related to anions or _M_1-cations ordering (Friis et al., Reference Friis, Balić-Žunić, Pekov and Petersen2004; Biagioni et al., Reference Biagioni, Hålenius, Pasero, Karlsson and Bosi2019). In fluoralforsite, the O2 sites are shifted from the unique position (.m) in order to obtain a favourable bond distance with the large Ba cation and the smaller Ca cation. Nevertheless, the chemical analyses indicate a slight dominance of barium at the _M_1 sites, which may induce local disorder.

In addition, the hypothetical cations ordering at _M_1 sites would also affect the O1 site, which, similarly to O2, lies on the mirror plane. The refined values of O1 atomic displacement parameters (Table 5) indicate an elongation along the _c_-axis, even though the amount of the displacement is less pronounced than that observed for O2 (Table 5). Apart from O2 and O1, additional sites show a displacement along the _c_-axis: _M_2, which is refined as a split position between _M_2 and _M_2A; and fluorine, which is also modelled by two sites, F and F1, 0.61 Å apart along c. (Table 4). The splitting of _M_2 was not observed in alforsite or in the synthetic counterpart, though it has been noticed in the apatite structure, where the Cl anion is shifted towards the mirror plane (Sudarsanan and Young, Reference Sudarsanan and Young1978; Hughes et al., Reference Hughes, Cameron and Crowley1990; Hughes and Rakovan, Reference Hughes, Rakovan, Kohn, Rakovan and Hughes2002). The _M_2 displacement in fluoralforsite might be triggered by the presence of Cl. This would also justify greater displacement of O2 with respect to O1: the _M_2–O2 distance is shorter (2.78 Å) than _M_2–O1 (3.04 Å). The re-arrangement of _M_2 has a stronger effect on the oxygen at O2 and O3 with respect to that at O1.

It is also conceivable that the splitting of the _M_2 position results from two overlapping effects, that is: (1) the local order at _M_1 sites and the subsequent split of O2, O3 and O1; and (2) the disorder of anions along [001].

Compared to fluorapatite, F– is displaced off the mirror plane in fluoralforsite, forming two symmetrical partially-occupied equivalent sites, corresponding to the hydroxyl group in hydroxylapatite (Mathew et al., Reference Mathew, Mayer, Dickens and Schroeder1979; Hughes et al., Reference Hughes, Harlov and Rakovan2018; Hazrah and Antao, Reference Hazrah and Antao2022). In the hexagonal apatite structure, the position of the anionic sites depends on the size of the anions and _M_2 cations (Hata et al., Reference Hata, Marumo, Iwai and Aoki1979; Hughes et al., Reference Hughes, Cameron and Crowley1990). Therefore, the splitting of the F site with respect to fluorapatite can be considered a result of the larger Ba size compared to Ca. However, F anions are displaced ~0.74 Å (considering an average position between F and F1 sites) below or above the mirror plane in fluoralforsite, whereas in synthetic Ba5(PO4)3F, the displacement is significantly smaller (~0.24 Å). The greater distance of F anions from the mirror plane is ascribed to the minor presence of the hydroxyl group.

Previous studies show that complete disordering along anion columns is unrealistic due to contacts being too short with neighbouring anions (Hughes et al., Reference Hughes, Cameron and Crowley1989, Reference Hughes, Cameron and Crowley1990; Biagioni et al., Reference Biagioni, Bosi, Hålenius and Pasero2017; Pieczka et al., Reference Pieczka, Biagioni, Gołębiowska, Jeleń, Pasero and Sitarz2018). The anion position in hexagonal ternary apatites depends on the adjacent anionic sites in the channel (Hughes and Rakovan, Reference Hughes, Rakovan, Kohn, Rakovan and Hughes2002). Therefore, the ordering arrangement of anions is required within individual columns, but the disordering of columns throughout the apatite structure allows the hexagonal symmetry to be maintained (Hughes et al., Reference Hughes, Cameron and Crowley1989; Hughes and Rakovan, Reference Hughes, Rakovan, Kohn, Rakovan and Hughes2002; Biagioni et al., Reference Biagioni, Bosi, Hålenius and Pasero2017). In fluoralforsite, the F sites, located above and below the mirror plane and Cl site cannot be occupied simultaneously (Table 7, Fig. 8). The arrangement of F and Cl can be as follows: F (above the mirror plane) – Cl – F (below mirror plane; Fig 7). The adjacency of two chlorines is implausible due to the low content in the structure.

Table 7. Distances (Å) between anions in the channels of fluoralforsite structure. Grey colour indicates too-short distances for simultaneous occupancy of both sites.*

Figure 8. (a) Unit cell (dotted line) of fluoralforsite projected along [110] with the mirror planes at ¼ z and ¾ z marked by solid grey lines. (b) Possible F– and Cl– arrangement along the channel parallel to [001]. Colour code as in Fig. 5. ‘b’ subscript – below mirror plane, ‘a’ subscript – above mirror plane. The dotted line indicates the _c_-axis. Drawn using _CrystalMaker_® software.

Considering the chlorine and hydroxyl group content, the sample studied can be regarded as a ternary barium phosphate. An ordered arrangement of Ba/Ca at the _M_1 site is very likely to occur within individual columns, however adjacent columns are disordered with respect to each other, similar to the anionic columns of hexagonal apatites (Hughes et al., Reference Hughes, Cameron and Crowley1989; Hughes and Rakovan, Reference Hughes, Rakovan, Kohn, Rakovan and Hughes2002). The ordering of cations at the _M_1 position is characteristic of belovite-group minerals, where the _M_1 site is split into _M_1 and _M_1’ and the symmetry is lowered to _P_63 or P bar3\bar{3}bar3 (Pasero et al., Reference Pasero, Kampf, Ferraris, Pekov, Rakovan and White2010). In our study, the chemical analyses and single-crystal X-ray diffraction data indicate the dominance of barium at _M_2 sites and slight occupancy at the _M_1 site, enabling the sample investigated to be included in the apatite group. A similar issue has been investigated by Biagioni et al. (Reference Biagioni, Hålenius, Pasero, Karlsson and Bosi2019) for hydroxylhedyphane Ca2Pb3(AsO4)3(OH) (P bar3\bar{3}bar3), in which the _M_1 site splits into _M_1 and _M_1’ but with Ca as the dominant cation at both sites. On that basis the authors classified the mineral as belonging to the hedyphane group and suggested the possible occurrence of intermediate phases between hedyphane and belovite groups (Biagioni et al., Reference Biagioni, Hålenius, Pasero, Karlsson and Bosi2019). Likewise, the data provided here about the new mineral species – fluoralforsite indicate the possibility of the existence of intermediate phases between belovite and apatite groups and suggest a possible new member of the belovite group with end-member formula BaCaBa3(PO4)3F and symmetry _P_63 or P bar3\bar{3}bar3.

Genesis of fluoralforsite

Fluoralforsite is a novel barium mineral from rankinite-rich paralavas that occurs in intergranular space and enclaves. Studies on the pyrometamorphic rocks of the Hatrurim Basin have distinguished two main types of enclaves, ferrite and silicate–phosphate. The origin of them is related to the crystallisation of residual melt between large crystals of paralava rock-forming minerals. The first type is dominated by minerals such as khesinite, Fe-spinels (magnetite, magnesioferrite, trevorite and cuprospinel), harmunite, barioferrite, perovskite and a nataliakulikite-like phase (Galuskina et al., Reference Galuskina, Galuskin, Pakhomova, Widmer, Armbruster, Krüger, Grew, Vapnik, Dzierażanowski and Murashko2017b; Krzątała et al., Reference Krzątała, Panikorovskii, Galuskina and Galuskin2018). Walstromite or zadovite minerals can be dominant in the second type – silicate–phosphate enclaves enriched in elements such as Ba, K, V, Nb, U and Sr. The walstromite-type enclaves, which co-occur with wollastonite and cuspidine have been described by Krzątała et al. (Reference Krzątała, Krüger, Galuskina, Vapnik and Galuskin2020). Fluoralforsite was found in zadovite enclaves, which are usually associated with fluorapatite (Fig. 3). The formation of a Ba-bearing enclave is a result of non-equilibrium between the primary minerals and residual melt, indicated by the melting of fluorapatite rims (Fig. 3d). It led to the enrichment of the silica-rich residual melt in phosphorus. Consequently, zadovite and fluoralforsite were formed at the final crystallisation stage of the paralava.

Fluoralforsite is a novel barium mineral from rankinite-rich paralavas that occurs in intergranular space and enclaves. Studies on the pyrometamorphic rocks of the Hatrurim Basin have distinguished two main types of enclaves, ferrite and silicate–phosphate. The origin of them is related to the crystallisation of residual melt between large crystals of paralava rock-forming minerals. The first type is dominated by minerals such as khesinite, Fe-spinels (magnetite, magnesioferrite, trevorite and cuprospinel), harmunite, bario ferrite, perovskite and a nataliakulikite-like phase (Galuskina et al., Reference Galuskina, Galuskin, Pakhomova, Widmer, Armbruster, Krüger, Grew, Vapnik, Dzierażanowski and Murashko2017b; Krzątała et al., Reference Krzątała, Panikorovskii, Galuskina and Galuskin2018). Walstromite or zadovite minerals can be dominant in the second type – silicate–phosphate enclaves enriched in elements such as Ba, K, V, Nb, U and Sr. The walstromite-type enclaves, which co-occur with wollastonite and cuspidine have been described by Krzątała et al. (Reference Krzątała, Krüger, Galuskina, Vapnik and Galuskin2020). Fluoralforsite was found in zadovite enclaves, which are usually associated with fluorapatite (Fig. 3). The formation of a Ba-bearing enclave is a result of non-equilibrium between the primary minerals and residual melt, indicated by the melting of fluorapatite rims (Fig. 3d). It led to the enrichment of the silica-rich residual melt in phosphorus. Consequently, zadovite and fluoralforsite were formed at the final crystallisation stage of the paralava.

In summary, rock-forming minerals such as melilites of the alumoåkermanite–gehlenite series, Ti-rich andradite, fluorapatite and calcium silicates such as rankinite and/or wollastonite crystallised from the primary melt. Subsequently, residual fluid/melt rich in halides, mainly fluorine and chlorine and incompatible elements reacted with early-formed minerals. These changes enhanced the mineral variety and contributed to the formation of new minerals such as fluoralforsite. It is noteworthy that the mineralisation of both ferrite and silicate–phosphate enclave types originate from similar processes – high-temperature alteration between pre-existing pyrometamorphic minerals and residual melt (Galuskina et al., Reference Galuskina, Galuskin, Pakhomova, Widmer, Armbruster, Krüger, Grew, Vapnik, Dzierażanowski and Murashko2017b; Galuskin et al., Reference Galuskin, Galuskina, Gfeller, Krüger, Kusz, Vapnik, Dulski and Dzierżanowski2016, Reference Galuskin, Gfeller, Galuskina, Armbruster, Krzątała, Vapnik, Kusz, Dulski, Gardocki, Gurbanov and Dzierżanowski2017).