Introduction
Of the 60 known uranyl–sulfate minerals, nearly half were first discovered during the last 10 years in the mines of Red Canyon in southeast Utah, USA. The Blue Lizard mine, in particular, has been a prolific source, and is now the type locality for 22 uranyl–sulfate minerals (see Plášil et al., Reference Plášil, Kampf, Ma and Desor2023; Kampf et al., Reference Kampf, Olds, Plášil and Marty2023a), with more awaiting characterisation. The new mineral libbyite, (NH4)2(Na2□)[(UO2)2(SO4)3(H2O)]2⋅7H2O, described herein, is the latest to be described from the Blue Lizard mine. Like several of the other new uranyl sulfates from here, libbyite contains a uranyl–sulfate structural unit that has not previously been found in Nature.
Libbyite is named in honour of American nuclear chemist Willard F. Libby (1908–1980) for his work on nuclear and radiochemistry. Dr. Libby's long and illustrious career, following a Ph.D. in chemistry from the University of California at Berkeley in 1933, included chemistry professorships at UC Berkeley, the University of Chicago (Institute for Nuclear Studies) and the University of California at Los Angeles (UCLA), where he was Director of the Institute of Geophysics and Planetary Physics (IGPP). During his tenure at the University of Chicago, Dr. Libby developed the method of radiocarbon dating (published in 1952) for which he was awarded the Nobel Prize in Chemistry for 1960.
The new mineral and name (symbol Ly) were approved by the Commission on New Minerals, Nomenclature and Classification of the International Mineralogical Association (IMA2022-091, Kampf et al., Reference Kampf, Olds, Plášil, Nash and Marty2023b). The description is based on four cotype specimens, all micromounts, deposited in the collections of the Natural History Museum of Los Angeles County, 900 Exposition Boulevard, Los Angeles, CA 90007, USA, catalogue numbers 76267, 76268, 76269 and 76270. Specimen 76267 is also a cotype for zincorietveldite (Kampf et al., Reference Kampf, Olds, Plášil and Marty2023a).
Occurrence
Libbyite was found by two of the authors (ARK and JM) in efflorescent crusts on mine walls underground in the Blue Lizard mine (37°33'26"N 110°17'44"W), Red Canyon, White Canyon District, San Juan County, Utah, USA. The mine is ~72 km west of the town of Blanding, Utah, and ~22 km southeast of Good Hope Bay on Lake Powell. Detailed historical and geological information on the Blue Lizard mine is described elsewhere (e.g. Kampf et al., Reference Kampf, Plášil, Kasatkin and Marty2015), and is derived primarily from a report by Chenoweth (Reference Chenoweth1993). Abundant secondary uranium mineralisation in Red Canyon is associated with post-mining oxidation of asphaltite-rich sandstone beds laced with uraninite and sulfides in the damp underground environment. Libbyite is a rare mineral found in association with bobcookite, coquimbite, halotrichite, metavoltine, rhomboclase, römerite, tamarugite, voltaite, zincorietveldite and other potentially new minerals on matrix comprised mostly of subhedral to euhedral, equant quartz crystals that are recrystallised counterparts of the original grains of the sandstone.
Morphology, physical properties and optical properties
Libbyite occurs as tightly intergrown aggregates of equant, somewhat rounded, light green–yellow crystals (Fig. 1). No crystal forms could be measured, but {001}, {011} and {111} appear likely. Merohedral twinning is likely because of the noncentrosymmetric space group, but was not observed. The streak is very pale green yellow. The mineral fluoresces strong green under 405 nm ultraviolet illumination. Crystals are transparent with vitreous lustre. The tenacity is brittle and the fracture is curved. The Mohs hardness is ~2½ based on scratch tests. Cleavage is excellent on {001}. The density could not be measured because the mineral is soluble in Clerici solution and there is insufficient material available for physical measurement. The calculated density based upon the empirical formula is 3.465 g⋅cm–3. The mineral is easily soluble in room-temperature H2O. Libbyite is optically uniaxial (–) with ω = 1.581(2) and ɛ = 1.540(2) measured in white light. The pleochroism is O = yellow, E = pale yellow; O > E. The Gladstone–Dale compatibility (Mandarino, Reference Mandarino2007) 1 – (K p/K c) is 0.001 (superior) based on the empirical formula using k(UO3) = 0.118, as provided by Mandarino (Reference Mandarino1976).
Raman spectroscopy
Raman spectroscopy was conducted on a Horiba XploRA PLUS using a 100× (0.9 NA) objective. Libbyite was very sensitive to the 532 nm diode laser and exhibited strong fluorescence. Consequently, the spectrum from 2000 to 60 cm–1 was recorded with a 785 nm laser (100 μm slit and 1800 gr/mm diffraction grating). The spectrum was featureless between 2000 and 1300 cm–1 so only the spectrum from 1300 and 100 cm–1 is shown in Fig. 2. The approximate wavenumbers and tentative band assignments are labelled in the expanded portion of Fig. 2. The band assignments are based primarily upon those for uranyl–sulfate minerals provided by Čejka (Reference Čejka, Burns and Finch1999) and Plášil et al. (Reference Plášil, Buixaderas, Cejka, Sejkora, Jehlicka and Novak2010).
The presence of three symmetrically distinct SO4 tetrahedra in the structure of libbyite leads to the multiple split bands for the SO4 modes. According to the empirical relationship of Bartlett and Cooney (Reference Bartlett and Cooney1989), the very strong ν1 (UO2)2+ symmetric stretching vibration at 864 cm–1 corresponds to an approximate U–OUr bond length of 1.75 Å, in good agreement with the average U1–OUr bond length from the X-ray data: 1.765 Å.
Chemical composition
Electron probe microanalyses (EPMA) were performed at the University of Utah on a Cameca SX-50 electron microprobe with four wavelength dispersive spectrometers and using Probe for EPMA software. Analytical conditions were 15 kV accelerating voltage, 10 nA beam current and 10 μm beam diameter and 3 points were analysed. Raw X-ray intensities were corrected for matrix effects with a ϕρ(z) algorithm (Pouchou and Pichoir, Reference Pouchou and Pichoir1991). No other elements were detected. There was major beam damage and, despite efforts to apply time-dependent corrections, it was obvious that they could not account for major rapid losses of N and Na. Consequently, we calculated (NH4)2O and Na2O based on the ideal formula (see below). The significant H2O loss under vacuum and during analyses resulted in higher concentrations for the remaining constituents than are to be expected for the fully hydrated phase; therefore, the other analysed constituents have been normalised to provide a total of 100% when combined with the calculated H2O content. Analytical data are given in Table 1.
* Based on N + K = 1, Na = 2 and O = 41 apfu.
S.D. = standard deviation.
The empirical formula (calculated on the basis of 41 O atoms per formula unit) is (NH4)1.92K0.08Na2.00U4.00S6.00O41H18.00. The simplified formula is (NH4,K)2(Na,□)3[(UO2)2(SO4)3(H2O)]2⋅7H2O and the ideal formula is (NH4)2(Na2□)[(UO2)2(SO4)3(H2O)]2⋅7H2O, which requires (NH4)2O 2.74, Na2O 3.26, UO3 60.20, SO3 25.27, H2O 8.53, total 100 wt.%.
X-ray crystallography
Powder X-ray diffraction (PXRD) data were recorded using a Rigaku R-Axis Rapid II curved imaging plate microdiffractometer with monochromatised MoKα radiation. A Gandolfi-like motion on the φ and ω axes was used to randomise the sample. Observed d values and intensities were derived by profile fitting using JADE Pro software (Materials Data, Inc.). The powder data are presented in Supplementary Table S1. The unit-cell parameters refined from the powder data using JADE Pro with whole pattern fitting are a = 10.7032(18), c = 31.864(6) Å and V = 3650.3(1) Å3.
The single-crystal structure data were collected at room temperature using the same diffractometer and radiation noted above. The best crystal found exhibited relatively high mosaicity. The mosaicity coupled with the close spacing of reflections along c caused problems in integration, which, in turn, forced us to drop some frames. This is the cause of the rather low completeness value of 92.7%. The structure data for libbyite were processed using the Rigaku CrystalClear software package, including the application of an empirical multi-scan absorption correction using ABSCOR (Higashi, Reference Higashi2001). The structure was solved using SHELXT (Sheldrick, Reference Sheldrick2015a). Refinement proceeded by full-matrix least-squares on F 2 using SHELXL-2016 (Sheldrick, Reference Sheldrick2015b). The large b parameter in the final weighting function reflects the relatively low quality of the data. Data collection and refinement details are given in Table 2, atom coordinates and displacement parameters in Table 3, selected bond distances in Table 4 and a bond-valence analysis 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).
R int = Σ|F o2–F o2(mean)|/Σ[F o2]. GoF = S = {Σ[w(F o2–F c2)2]/(n–p)}½. R 1 = Σ||F o|–|F c||/Σ|F o|. wR 2 = {Σ[w(F o2–F c2)2]/Σ[w(F o2)2]}½; w = 1/[σ2(F o2)+(aP)2+bP] where a is 0.0739, b is 169.7151 and P is [2F c2+Max(F o2,0)]/3.
* Refined occupancies: N = N0.88(6)/K0.12(6); Na2 = 0.90(6)
*NH4+–O bond valence parameters are from García-Rodríguez et al. (Reference García-Rodríguez, Rute-Pérez, Piñero and González-Silgo2000); U+6–O and S+6–O bond-valence parameters are from Gagné and Hawthorne (Reference Gagné and Hawthorne2015). Hydrogen-bond strengths based on O–O bond lengths are from Ferraris and Ivaldi (Reference Ferraris and Ivaldi1988). Negative values indicate donated hydrogen-bond contributions.
Description of the structure
Two U sites (U1 and U2) in the structure of libbyite are each surrounded by seven O atoms forming squat UO7 pentagonal bipyramids. This is a typical coordination for U6+ in which the two short apical bonds of the bipyramid constitute the uranyl group (see Burns, Reference Burns2005). The two apical O atoms of the bipyramids (OUr) form short bonds with the U, and this unit comprises the UO22+ uranyl group. Five equatorial O atoms (Oeq) complete the U coordinations.
There are three S sites (S1, S2 and S3) each centring an SO4 tetrahedron. The SO4 tetrahedra share corners with the equatorial O atoms of the UO7 bipyramids to form a uranyl–sulfate sheet with the composition [(UO2)2(SO4)3(H2O)]2– (Fig. 3). Within this sheet, the U1 bipyramid shares four of its Oeq corners with SO4 groups and the U2 bipyramid shares all five of its Oeq corners with SO4 groups. The uranyl–sulfate sheet in the structure of libbyite is unique in the mineral kingdom, however it has the same topology as the sheet in the synthetic phase K(H5O2)[(UO2)2(SeO4)3(H2O)] (Gurzhiy et al., Reference Gurzhiy, Tyumentseva, Krivovichev, Tananaev and Myasoedov2012) and as those in seven other synthetic uranyl selenates listed by Krivovichev (Reference Krivovichev2009) with graph cc2–2:3–4 and ring symbol 4361. The charge-deficiency-per-anion (CDA), which is one of the useful measures to quantify the bond-valence characteristics of the structural units of minerals, is 0.14 valence units for libbyite. It is noteworthy that the value found for such an ‘exotic’ sheet is not far from the CDA values observed for other sheet uranyl sulfates, for example, for some of the zippeite group of minerals (see Plášil et al., Reference Plášil, Kampf, Ma and Desor2023).
The region between the uranyl–sulfate sheets (Fig. 4) contains one NH4 site (N), two Na sites (Na1 and Na2) and four H2O sites (OW2, OW3, OW4 and OW5). The NH4 site is seven coordinated (for N–O < 3.3 Å). The Na1 site is six coordinated and the Na2 site is seven coordinated. The sheets are linked to each other in the [001] direction via NH4–O, Na–O and hydrogen bonds.
The structural formula, [(NH4)1.75K0.25](Na2.81□0.19)[(UO2)2(SO4)3(H2O)]2⋅7H2O, based on our refinement has an excess positive charge of 0.81. In the final stages of refinement, we allowed interlayer cation occupancies to refine freely for an indication of their preferred site occupancies. The Na1 site refined to full occupancy and the Na2 site to 0.90 occupancy; however, both Na sites, as well as the NH4 site, exhibited large anisotropic displacement parameters. This suggests that the total of the refined Na and NH4 site occupancies is larger than reality. To test this, a refinement was done with the site occupancies adjusted to Na1: 0.62, Na2: 0.69, N: N0.96K0.04 corresponding to the structural formula [(NH4)1.92K0.08](Na2.00□1.00)[(UO2)2(SO4)3(H2O)]2⋅7H2O, which is the same as the EPMA empirical formula. The refinement was well behaved, converging to R 1 = 0.0590, with reasonable U eq values of 0.022, 0.027 and 0.055 for the Na1, Na2 and N sites, respectively.
Electron probe microanalysis of libbyite was very challenging because the mineral is very sensitive to the electron beam. Preliminary energy dispersive spectroscopy (EDS) indicated sufficient N to account for full occupancy of the NH4 site in the structure and the structure refinement indicated full occupancy of the site together with greater K than indicated by the EPMA. The preliminary EDS also indicated ~2.5 Na atoms per formula unit. Considering the results of the structure refinement noted in the previous paragraph, we believe that the most realistic ideal formula for libbyite is (NH4)2(Na2□)[(UO2)2(SO4)3(H2O)]2⋅7H2O.
New structure types such as those found in libbyite provide important metrics when surveying for potential relationships between the observed crystal chemistry and the, usually unmeasured, conditions of formation. The diverse and densely co-mingled associations of uranyl sulfates have so far made the determination of paragenesis and the measurement of conditions of formation impossible for most species. However such an understanding is valuable when evaluating the long-term disposal of nuclear waste for certain proposed repositories, where highly soluble uranyl–sulfate minerals could form and enhance radionuclide mobility. The rarest Red Canyon uranyl sulfates occur in just one or two specimens that were confined to a very small footprint of efflorescence underground, and have afforded only a few crystals for analyses, limiting our understanding of chemical and substitutional variability. Most species have yet to be reproduced synthetically and recent efforts to do so proceed essentially by luck, with little guidance from Nature. These difficulties necessitate continued underground collecting of uranyl sulfates, together with solution data that may complement synthetic work. Sadly, recent closures of mines in Red Canyon prevent this research.
Acknowledgements
Structures Editor Peter Leverett and two anonymous reviewers are thanked for their constructive comments on the manuscript. We are grateful to retired miner Dan Shumway of Blanding, Utah, for advice and assistance in our collecting efforts in Red Canyon. This study was funded, in part, by the John Jago Trelawney Endowment to the Mineral Sciences Department of the Natural History Museum of Los Angeles County. JP acknowledges the support of the Czech Science Foundation (GACR 20-11949S).
Supplementary material
The supplementary material for this article can be found at https://doi.org/10.1180/mgm.2023.26.
Competing interests
The authors declare none.