Bobtraillite from Gejiu hyperagpaitic nepheline syenite, southwestern China: new occurrence and crystal structure

. A second occurrence of bobtraillite is described from the Gejiu nepheline syenite, southwestern China. The extremely rare and complex boron-bearing zirconium silicate is associated with albite, orthoclase, jadeite, ﬂuorite, andradite, titanite, as well as other REE and zirconium-bearing minerals, catapleiite, moxuanx-ueite, låvenite, eudialyte, britholite-(Ce)


Introduction
Bobtraillite is an extremely rare cyclosilicate with a unique composition and complex structure that has previously only been reported from a single locality.Only two boroncontaining zirconium silicates have been found in nature so far: bobtraillite and rogermitchellite.According to the total structural information content, they are both classified into complex mineral structure (> 1000 bits per unit cell) (Grew et al., 2016); especially, rogermitchellite has the most complex crystal structure of all the boron minerals with I G, total = 3019 bits and is the 20th most complex of the 4443 structures reported to data (Krivovichev, 2013;Krivovichev et al., 2022).The two chemically and structurally related minerals are both found in hyperagpaitic nepheline syenites at Mont Saint-Hilaire, Quebec, Canada.The Canadian bobtraillite was originally described by McDonald and Chao (2005), the electron microprobe analyses of bobtraillite yielding an empirical formula (Na 11.20 Ca 1.22 ) 12.42 (Sr 10.59 Ba 0.16 ) 10.75 (Zr 12.69 Y 0.63 Nb 0.61 Hf 0.14 ) 14.07 Si 41.64 B 6 O 132 (OH) 12 q 12H 2 O, and the ideal formula was given as (Na,Ca) 13 Sr 11 (Zr,Y,Nb) 14 Si 42 B 6 O 132 (OH) 12 q 12H 2 O. Bobtraillite was considered as a structural hybrid of benitoite-group minerals due to its similarities in chemical composition and crystal structure.However, it is also unique among all those members for the BO 4 tetrahedra, which are linked to isolated ZrO 6 octahedra.
In this paper, we describe a second occurrence of bobtraillite from Gejiu sodalite-nepheline syenite in Yunnan Province, southwestern China.We present a detailed description of the mineral paragenesis, Raman spectroscopy, complete chemical analysis, and crystal structure refinements for the extremely rare complex zirconium silicate.

Occurrences and paragenesis
Bobtraillite was found in the Gejiu alkaline intrusive complex (23 • 29 40 N, 103 • 4 41 E).This alkaline complex is associated with gabbro, granitic, and alkaline intrusions in the western part of the Cathaysia block, southwestern China (Fig. 1a).The alkaline intrusions consist mainly of early-formed silica-saturated alkali syenites from the northeastern part (Baiyunshan intrusions) and later-formed silicaundersaturated feldspathoid syenites from the northern part of the complex (Changlinggang intrusions) (Fig. 1b).The studied material is found in the feldspathoid (nephelinesodalite) syenite.The studied material was found in hand specimen 18CL18 collected by Wang et al. (2021) near the location for their sample 18CL07.Wang et al. (2021) gave a detailed account of the geology of the Changlinggang intrusions.

Raman spectroscopy
Raman spectra were collected using polished grains, and analyses were performed at Tianjin Center, China Geological Survey using a Renishaw in Via Raman microscope with a 532 nm laser (40 mW, 1 µm) and a measurement range of 4000-150 cm −1 (Fig. 3).In general, the medium-intensity band between 3467-3585 cm −1 is a superposition of O-H stretching bands of OH and H 2 O.The strong and broad bands at 1522, 1766, 1989, and 2374 cm −1 , which are similar to the Raman spectra of other hydrous materials containing acid groups (Chukanov et al., 2022), are assigned to BO-H stretching vibrations.These bands correspond to extremely strong hydrogen bonds in hydrated proton complexes involving proton formed as a result of partial dissociation of an acid group.The strong sharp band at 1006 cm −1 is assigned to v 3 (SiO 4 , BO 4 ) antisymmetric stretching, and the weak shoulders at 983 and 1091 cm −1 arise from the symmetric v 1 (SiO 4 , BO 4 ) stretching vibrations, and relatively weak peaks at 630 cm −1 can be identified as v 2 (SiO 4 , BO 4 ) symmetric bending.The band at 517 cm −1 can be due to Zr-O stretching vibrations.The weak peak at 366 cm −1 , which is commonly found in the zirconium silicate, may arise from v 4 (SiO 4 , BO 4 ) antisymmetric bending or external lattice vibrations, whereas the weak band at 201 cm −1 can be assigned to lattice vibrations.

Chemical composition
Bobtraillite was identified in the hand specimen 18CL18.The polished section 18CL18-5 was analyzed at the Beijing Research Institute of Uranium Geology with a JXA-8100 electron microprobe under operating conditions of 15 kV accelerating voltage, 10 nA beam current, and a 5 µm beam diameter.The standards used were albite (Na, Al), phlogopite (K, Fe, Mg, Si), monazite-(Ce) (P, REE, Th), plagioclase (Ca), celestite (Sr), zircon (Zr), rutile (Ti), yttrium Al garnet (Y), pure metallic Nb (Nb), and pure metallic U (U). was determined by laser ablation inductively coupled plasma mass spectrometry (LA-ICP-MS) analysis.The analysis was carried out at Tianjin Center, China Geological Survey, using an Agilent 7900 equipped with a RESOlution LR laser ablation system.The operating conditions were as follows: beam diameter = 29 µm, alongside a laser pulse rate of 6 Hz with an energy density of approximately 3 J cm −2 .NIST SRM 610 was used for the standard.The measured content of B 2 O 3 is consistent with the result (B 2 O 3 3.13 %) calculated as the difference of Si + B = 48 apfu.H 2 O (5.22 %) was calculated by stoichiometry from the results of the crystal structure analysis, and the hydrous nature was confirmed by means of OH stretching vibration absorption in the Raman spectrum.Complete analytical results are given in Table 1.

X-ray crystallography
Single-crystal X-ray studies were carried out at the School of Geosciences and Info-Physics, Central South University, with a Rigaku XtaLAB Synergy diffractometer and CuKα radiation at 50 kV and 1 mA.The intensity data were corrected for X-ray absorption using multi-scan method, and empirical absorption correction was performed us-    ing CrysAlisPro software spherical harmonics, which was implemented in SCALE3 ABSPACK scaling algorithm.
The systematic absences suggest that for the space group P 3c1, the refined cell parameters are a = 19.6977(6)Å, c = 9.9770(3) Å, and V = 3352.4(2)Å 3 .The structure was solved with Olex2 (Dolomanov et al., 2009) and refined with the SHELXL-program (Sheldrick, 2015).Scattering factors for neutral atoms were used initially: Si and Zr sites were fully occupied by Si and Zr, respectively, whereas the B site was refined with B vs. Si, the Na(1) and Na(2) sites were refined with Na vs. , and the Sr site was refined with Sr vs. Ca.After several cycles of anisotropic refinement for all the atoms, the final R 1 was 0.0586 for 1908 independent reflections with F o > 4σ (F o ).Crystal data collection and refinement details are reported in Table 2. Atom coordinates and isotropic displacement parameters are given in Table 3, anisotropic displacement parameters in Table 4, and selected bond distances in Table 5.The bond-valence calculations are obtained using the bond-valence parameters of Brese and O'Keeffe (1991), in Table 6.
In the crystal studied, there was a mixed occupancy by B (0.92 apfu) and Si (0.08 apfu) at the B site, while other four Si sites are fully occupied by Si.The average Si-O bond lengths range from 1.609 to 1.630 Å, and average B-O/OH  bond distance is 1.499 Å.The BVS calculated for Si and B are within error the same as the ideal values (Table 6).14)] (Fig. 4b).The BVS calculated for Na(1) and Na(2) site is 0.889 and 0.157 v.u., respectively, in perfect agreement with their refined occupancies.Interestingly, when we quantitatively analyzed the B content in bobtraillite, we unexpectedly discovered that Li is also present, which was not analyzed in the bobtraillite from Mont Saint-Hilaire.It successfully explains the insufficient cations in the Na site of bobtraillite from Gejiu syenite.

Ideal formula
The empirical formula of bobtraillite from Gejiu syenite was calculated on the basis of 156 O atoms per formula unit (apfu); B content was calculated as the difference of Si + B = 48 apfu; H content was estimated by stoichiometry in order to achieve the electrostatic balance.

Genesis of bobtraillite
Although only two occurrences of bobtraillite have been reported so far, it is almost certainly a secondary mineral.This is strongly demonstrated in the original description of bobtraillite from the type locality by McDonald and Chao (2005); there it is very clearly shown that fine-grained bobtraillite, associated with other secondary carbonates and other hydrous minerals, occurs in cavities in the early-formed alkaline breccias.In our case, bobtraillite occurs in the fineto coarse-grained nepheline-sodalite syenites.The assemblages where bobtraillite occurs usually form pseudomorphs of titanite or moxuanxueite, indicating it is the decomposition product of the primary minerals.It should be noted that jadeite usually occurs in the central part of the aggregates (Fig. 2), indicating bobtraillite formed at the same stage, and thus the jadeite should play a key role in discussing the alteration process.In silica-undersaturated rocks, where jadeite could form through the following Reaction (R1) given by Boettcher and Wyllie (1968): The Reaction (R1) between nepheline and albite to produce jadeite is in the presence of excess H 2 O (Gupta, 2015).This is consistent with the observed albite always occur-ring around mineral pseudomorph assemblage in the laterformed silica-undersaturated feldspathoid syenites from the Changlinggang intrusions.Previous studies have experimentally confirmed that this jadeite could form at a lower pressure in silica-undersaturated rocks than in silica-saturated rocks (Popp and Gilbert, 1972).Andersen et al. (2010) suggested that catapleiite could replace primary zirconium minerals under a relatively high H 2 O activity environment.The increasing activity of water can also be substantiated by the replacement of primary anhydrous moxuanxueite by late-formed hydrous bobtraillite and catapleiite.This process could occur at a relatively low temperature (≤ 348 • C) (Chakrabarty et al., 2016), which is significantly lower than the magma crystallization temperature of Gejiu feldspathoid syenites (898 • C) (Wang et al., 2021).
In addition, jadeite could also form from the breakdown of analcime according to the Reaction (R2) given by Kennedy (1961).The reaction could explain the paucity of hydrous analcime in the Gejiu feldspathoid syenites within the hydrothermal assemblage.The release of H 2 O, with the Ca, Sr, and Zr elements released from the breakdown of moxuanxueite to form the secondary hydrous zirconium silicates, for example, bobtraillite and catapleiite.It is worth noting that analcime is considered a decomposition product of nepheline (Henderson and Gibb, 1977).That means there are at least two stages of hydrothermal process.The first stage is https://doi.org/10.5194/ejm-35-65-20233.279(6) Mean [8]  2.671(6) Mean [10]  2.792( 6 the replacement of primary nepheline and alkali feldspar by analcime, and then the secondary analcime continues to alter to jadeite. The source of the B is still an open question.However, another boron silicate (stillwellite) has been found in both the Gejiu and the Mont Saint-Hilaire hyperagpaitic syenites (McDonald and Chao, 2005), which is considered to crystallize during the late pegmatitic and succeeding hydrothermal stages (Bailey, 2006;Chakrabarty et al., 2016).Furthermore, Li may be a more abundant constituent in hyperagpaitic rocks than is commonly realized.The Gejiu nepheline-sodalite syenite contains 107 ppm Li 2 O (Wang et al., 2021), and several lithium minerals have been found associated with the hyperagpaitic intrusions at Mont Saint-Hilaire (Van Velthuizen and Chao, 1989;Khomyakov et al., 1990;Chao et al., 1991;Chao and Ercit, 1991;McDonald et al., 2013), where the mineral assemblages resemble those in the Gejiu alkaline intrusive complex.The possibility of Li enrichment in hyperagpaitic intrusive complexes reminds us that cation-deficient minerals in such rocks, which have only been analyzed by the electron microprobe, need to be checked for Li, Be, and B, constituents that either cannot be analyzed with the electron microprobe or analyzed only with great difficulty.These under-reported light elements often play a crucial role in oreforming processes during the late stages of magmatic crystallization (Alferyeva et al., 2011).
Data availability.The Crystallographic Information File data of bobtraillite are available in the Supplement.
Author contributions.YW wrote the manuscript.FN improved the structure results and manuscript.ZH oversaw the research.XG carried out single-crystal X-ray diffraction and the structural refinement.GD and ZY supplied the specimens, images, and geological background information.GF conducted the microprobe analyses.ZX measured the Li 2 O (%) and B 2 O 3 (%) content using LA-ICP-MS.KQ carried out the scanning electron microscope (SEM) studies and assisted in all data processing.
Competing interests.The contact author has declared that none of the authors has any competing interests.
Disclaimer.Publisher's note: Copernicus Publications remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Acknowledgements.Kimberly Tait, Nikita V. Chukanov, and Associate Editor Edward Grew are thanked for their constructive comments on the manuscript.The authors thank Zongyao Tai for help with collecting electron microprobe (EMP) data and Lijuan Ye for Raman spectrum studies.Yufei Wang is thanked for the field work and the sample collecting.Review statement.This paper was edited by Edward Grew and reviewed by Nikita V. Chukanov and Kimberly Tait.

Figure 4 .
Figure 4.The crystal structure for bobtraillite (Na ions shown in yellow, Sr in gray, O in red, OH in pink, H 2 O in cherry red, and vacancy in white).

Financial support .
This research was supported by National Key R&D Programmes (grant nos.92062217 and 92062105), National Natural Science Foundation of China (NSFC) (grant no.42072054), and China Scholarship Council (CSC) (grant nos.202106400047 and 202108575009).
Notes: bond valence sums were calculated with the site-occupancy factors given in Table4.Calculations were done using the equation and constants ofBrown (1977),S = exp[(R 0 − d 0 )/b].

Figure 5 .
Figure 5.The giant pinwheel motif view of bobtraillite along [001] (four-unit cells, the dotted line is a single unit cell).The Zr(1)O 6 octahedra are shown in green, the Zr(2)O 6 octahedra in olive green, the SiO 4 tetrahedra in indigo, and BO 4 tetrahedra in green-blue.Other legends as in Fig. 4.
Notes: a B 2 O 3 was calculated by stoichiometry based on Si + B = 48; b H 2 O was calculated by stoichiometry.
Chao, 2005) material corresponds to the chemical formula (Na 9.70 Li 0.42 K 0.08 ) 10.20 (Sr 10.61 Ca 1.14 Fe 0.07 ) 11.82 (Zr 12.87 Ti 0.53 Nb 0.31 RE 0.14 U 0.02 Th 0.01 ) 13.88 (Si 42.41 B 5.59 Al 0.02 ) 48.02 O 132 (OH) 12 q 12H 2 O. Differently from previous studies (McDonald andChao, 2005), if we introduce a small number of divalent calcium ions into the two Na sites of our studied material, the R value increases.In addition, we are confident in stating that our BVS calculation is in agreement with the absence of Ca at the two Na sites.Instead, all the Ca is substituting for Sr at the Sr site.Taking into account the results of the crystal structure refinement, the ideal end-member formula of bobtraillite could be written as (Na, ) 12 ( , Na) 12 Sr 12 Zr 14 (Si 3 O 9 ) 10 [Si 2 BO 7(OH) 2 ] 6 q 12H 2 O.In order to keep the charge balance of the end-member formula, the total site population of Na is set to 10 apfu, ideally corresponding to (wt %) Na 2 O 4.89, SrO 19.63, ZrO 2 29.40, SiO 2 25.16, B 2 O 3 3.30, H 2 O 5.12, and a total of 100.