Introduction
The family of hydrated sodium-calcium sulfates with the general formula NaxCay(SO4)0.5x+y·nH2O includes a dozen natural and synthetic phases. This family demonstrates a truly enigmatic case: despite chemical simplicity, the majority of these compounds are poorly characterised. Until recently, four valid mineral species belonging to this family were known: wattevilleite, Na2Ca(SO4)2·4H2O (?) (Palache et al., Reference Palache, Berman and Frondel1951; Sabelli and Trosti-Ferroni, Reference Sabelli and Trosti-Ferroni1985; Leverett and Williams, Reference Leverett and Williams2007), hydroglauberite, Na10Ca3(SO4)8·6H2O (Slyusareva, Reference Slyusareva1969), eugsterite, Na4Ca(SO4)3·2H2O (Vergouwen, Reference Vergouwen1981), and omongwaite, Na2Ca5(SO4)6·3H2O (Mees et al., Reference Mees, Hatert and Rowe2008). None of these minerals were studied using a single-crystal diffraction method due to the imperfectness of their crystals. Moreover, unit-cell data are still unknown for wattevilleite, hydroglauberite and eugsterite. Among synthetic hydrated Na–Ca sulfates, the crystal structure has only been solved for NaCa3(SO4)3.5·1.7H2O which was found to be structurally similar to bassanite CaSO4·0.5H2O (Freyer et al., Reference Freyer, Reck, Bremer and Voigt1999). This synthetic Na–Ca sulfate and omongwaite have close powder X-ray diffraction patterns that allowed us to determine unit-cell parameters of the latter and consider its crystal structure as bassanite-related (Mees et al., Reference Mees, Hatert and Rowe2008).
The present paper is devoted to the new mineral species antofagastaite (Cyrillic: антофагастаит), ideally Na2Ca(SO4)2·1.5H2O, which belongs to the same chemical family but significantly differs from all its other representatives as well as from known synthetic Na–Ca sulfates. The mineral is named after the Antofagasta Province (Northern Chile) in which its type locality, the Coronel Manuel Rodríguez mine is situated. We note that the similar name ‘antofagastite’ was proposed in 1938 for a mineral with the formula CuCl2·2H2O, considered initially as a new species (Palache and Foshag, Reference Palache and Foshag1938). However, it was later shown to be identical to the earlier described eriochalcite and the name ‘antofagastite’ was therefore discredited (Palache et al., Reference Palache, Berman and Frondel1951).
Both the new mineral antofagastaite and its name have been approved by the IMA Commission on New Minerals, Nomenclature and Classification (IMA-CNMNC), IMA2018-049. The type specimen is deposited in the collection of the Fersman Mineralogical Museum of the Russian Academy of Sciences, Moscow, with the catalogue number 96263.
After completion of the investigation of antofagastaite from the Coronel Manuel Rodríguez mine (the holotype) and its approval as a new species by the IMA-CNMNC, the same mineral was identified by us in the material from the Arsenatnaya fumarole at the Tolbachik volcano, Kamchatka, Russia, and studied including the crystal structure determination. We consider it of interest to include a brief description of this, the second find of antofagastaite in the present paper.
Occurrence and general appearance
Specimens which became the type material of the new mineral were collected by us in January 2016 in dumps of the abandoned Coronel Manuel Rodríguez copper mine near the city of Antofagasta, Mejillones Peninsula, Mejillones, Antofagasta Province, Antofagasta Region, Chile (Mills et al., Reference Mills, Kampf, Dini and Molina2012).
Antofagastaite from the type locality is a supergene mineral formed in the oxidation zone of sulfide–quartz veins under conditions of extremely arid climate, typical for the Atacama Desert. It occurs in cavernous aggregates mainly composed of sideronatrite (partly altered to metasideronatrite), aubertite and gypsum. Other associated minerals are ferrinatrite, glauberite, amarillite and an unidentified Fe phosphate. Quartz is a relict gangue mineral found in this assemblage.
At the Coronel Manuel Rodríguez mine antofagastaite occurs as prismatic crystals up to 0.5 mm × 1 mm × 5 mm (typically less, up to 1 mm long) elongated along [010] and commonly flattened on [100]. They are usually coarse, sometimes blocky, slightly divergent and spear-like. The major form of antofagastaite crystals is the pinacoid {100} (Figs 1a and b). The lateral faces were not indexed. Faces of the pinacoid {010} are observed on well-terminated crystals (Fig. 1c). Parallel or near-parallel intergrowths of crystals are observed (Fig. 1c) as well as chaotic, open-work aggregates up to 1 cm across occurring in cavities. The new mineral overgrows sideronatrite/metasideronatrite and aubertite or forms intimate intergrowths with these hydrous Na–Fe sulfates (Figs 1d and 2). Some antofagastaite crystals are sprinkled by tiny tabular crystals of an unidentified Fe phosphate and contain inclusions of glauberite and amarillite (Figs 1b and 3).
Morphology of crystals and aggregates of antofagastaite from the Coronel Manuel Rodríguez mine: (a) cluster of flattened prismatic crystals; (b) board-shaped crystal with overgrowing tiny crystals of an unidentified Fe phosphate; (c) parallel intergrowth of prismatic crystals; (d) open-work group of slightly divergent prismatic crystals with metasideronatrite crystals (lighter). Scanning electron microscopy (SEM) images: (a)-(c) – in secondary electron (SE) mode, (d) – in back-scattered electron (BSE) mode, specimen no. CMR-5052.

Colourless prismatic crystals of antofagastaite from the Coronel Manuel Rodríguez mine: (a) on yellowish metasideronatrite; (b) on orange–yellow sideronatrite (partially altered to metasideronatrite) with blue aubertite. Field of view width: (a) 3.8 mm and (b) 8.9 mm. Photo: I.V. Pekov and A.V. Kasatkin, specimen no. CMR-5052.
Polished section of a cluster of antofagastaite crystals (1; porosity is very visible) with glauberite (2), amarillite (3) and an unidentified Fe phosphate (4). Coronel Manuel Rodríguez mine. SEM (BSE) image, specimen no. CMR-5052.

At the Tolbachik volcano, antofagastaite was found in the Arsenatnaya fumarole located on the summit of the Second scoria cone of the Northern Breakthrough of the Great Tolbachik Fissure Eruption of 1975–1976. The description of this active and richly mineralised fumarole has been given by Pekov et al. (Reference Pekov, Koshlyakova, Zubkova, Lykova, Britvin, Yapaskurt, Agakhanov, Shchipalkina, Turchkova and Sidorov2018). Antofagastaite occurs in the upper, moderately hot zone (temperatures measured by us in different areas of this zone are 70–150°C) of the fumarole in which the complex interactions between volcanic gas, atmospheric agents (water, water vapour and CO2) and earlier formed sublimate minerals take place. Antofagastaite forms in cavities long prismatic, lath-shaped to acicular crystals up to 0.03 mm × 0.3 mm associated closely with syngenite (Fig. 4) and blödite. An earlier, high-temperature assemblage of sublimate minerals include: hematite, anhydrite, aphthitalite, fluorophlogopite, sanidine (As-bearing), johillerite, svabite, tilasite and zubkovaite Ca3Cu3(AsO4)4 (IMA2018-008, Pekov et al., Reference Pekov, Lykova, Agakhanov, Belakovskiy, Vigasina, Britvin, Turchkova, Sidorov and Scheidl2019).
Abundant lath-shaped to acicular crystals of antofagastaite on massive syngenite crystals. Arsenatnaya fumarole, Tolbachik volcano. SEM (SE) image, specimen no. Tolb-6225.

Physical properties and optical data
Antofagastaite is transparent and colourless with a white streak. It has a vitreous lustre and does not fluoresce under ultraviolet light. It is brittle, the Mohs’ hardness is ca 3. There are two directions of distinct cleavage observed under the microscope, one of which is definitely (001); the fracture is uneven. Density measured by flotation in heavy liquids (bromoform + hexane) is 2.42(1) g cm−3. Density calculated using the empirical formula is 2.465 g cm−3. Some lower values of the measured density in comparison with the calculated one is caused by a slight porosity, observable in the crystals (Fig. 3).
Antofagastaite is optically biaxial (–), α = 1.489(2), β = 1.508(2), γ = 1.510(2) (589 nm), 2Vmeas. = 40(10)° and 2Vcalc. = 36°. Dispersion of optical axes is strong with r > v. Optical orientation is Z = a. Under the microscope the crystals are colourless and nonpleochroic.
Infrared spectroscopy
In order to obtain an infrared (IR) absorption spectrum of antofagastaite (holotype: Fig. 5), the powdered sample of the mineral was mixed with dried KBr, pelletised, and analysed using an ALPHA FTIR spectrometer (Bruker Optics) with a resolution of 4 cm−1 and 16 scans. The IR spectrum of an analogous pellet of pure KBr was used as a reference.
Powder infrared absorption spectra of (a) antofagastaite from the Coronel Manuel Rodríguez mine and (b) syngenite from Kalush, Precarpathian region, Ukraine [drawn using data from (Chukanov, Reference Chukanov2014)]. The weak band at ~1420 cm−1 in the IR spectrum of syngenite corresponds to an admixed carbonate.

Wavenumbers of absorption bands in the IR spectrum of antofagastaite and their assignments (cm−1; s – strong band, w – weak band, sh – shoulder) are: 3598 (O–H-stretching vibrations of water molecules, a weak Hbond), 3373 and 3300sh (O–H-stretching vibrations of water molecules, medium-strengths Hbonds) 1685w, 1629w (bending vibrations of water molecules), 1220, 1153s, 1128s, 1090s and 1060sh [ν3(F2) – asymmetric stretching vibrations of SO42− anions], 1015sh and 982 [ν1(A1) – symmetric stretching vibrations of SO42− anions), 758 (librational vibrations of water molecules: see Gerakines et al., Reference Gerakines, Schutte, Greenberg and van Dishoeck1995), 656 and 613 [ν4(F2) – bending vibrations of SO42− anions], 480 and 458 [ν2(E) bending mode of SO42− anions]. Very weak absorptions in the range 2100–2250 cm−1 correspond to first overtones and combination modes of the SO42− stretching modes.
Two bands of nondegenerate modes of H2O molecules observed at 1685 and 1629 cm−1 indicate the presence of two kinds of significantly non-equivalent water molecules. The ν1 symmetric stretching mode is inactive in the IR spectrum of MO4 groups having the configuration of a regular tetrahedron, and the corresponding band intensity increases with the increase of distortion degree. Consequently, the strongest band in the range of symmetric stretching vibrations of SO42– anions (i.e. in the range 960–1025 cm−1) observed at 982 cm−1 is to be assigned to the most distorted S2B-centred tetrahedron S2BO4 (see Table 4). Bands of stretching vibrations of BO33–, CO32– and NO3− anions (in the range from 1250 to 1550 cm−1) are not observed.
The following empirical correlations between O–H stretching frequencies in IR spectra of minerals and O···O and H···O distances (from structural data) were established by Libowitzky (Reference Libowitzky1999):
It should be noted that at high frequencies (above 3500 cm−1) substantial deviations from the correlations (1) and (2) are possible because O–H stretching frequencies depend not only on O···O and H···O distances, but also on the nature of cations coordinating O–H groups and H2O molecules, as well as on the angle O–H···O, and the influence of these factors becomes most evident in the case of weak hydrogen bonds. The equations 1 and 2 predict that maximum possible values of O–H stretching frequencies for minerals are 3592 and 3632 cm−1, respectively. However, e.g. for magnesium serpentines, brucite and kaolinite observed frequencies are close to 3700 cm−1.
According to the correlation (1), the band at 3598 cm−1 in the IR spectrum of antofagastaite formally corresponds to the absence of hydrogen bond: d(O···O) = ∞. According to the correlation (2), the band at 3598 cm−1 should be assigned to a very weak hydrogen bond: d(H···O) = 2.333 Å. Consequently, the peak at 3598 cm−1 nearly conforms to the overlapping bands which correspond to the hydrogen bonds with the OW2···O1 and OW2–OW1 distances of 3.16 and 3.09 Å, respectively (see description of the crystal structure below). The absorption maximum at 3373 cm−1 and the shoulder at 3300 cm−1 correspond to d(O···O) distances of 2.78 and 2.74 Å, respectively, which is close to the OW1···O4 distance of 2.773 Å. Consequently, these bands can be assigned to symmetric and antisymmetric vibrations of OW1H2 molecules.
The IR spectrum of antofagastaite shows a remote resemblance to that of the related mineral syngenite (curve (b) in Fig. 5). The absence of the band at ~3600 cm−1 (corresponding to weak hydrogen bands) in the IR spectrum of syngenite indicates the absence of additional H2O molecules like those occurring at the site OW2 in the structure of antofagastaite. The singlet of bending vibrations of water molecules observed at 1680 cm−1 in the spectrum of syngenite is in agreement with the presence of a single site of H2O in this K–Ca sulfate. Consequently, the additional (as compared to syngenite) bands at 1629 and 3598 cm−1 can be used for the evaluation of the occupancy of the OW2 site in antofagastaite.
Most bands, especially those of librational and stretching vibrations of OW1, in the IR spectrum of antofagastaite (at 758 and in the range 3300–3400 cm−1, respectively) are somewhat broadened compared to the similar bands of syngenite. This may be due to a disordering of the structure (including disordering of hydrogen bonds formed by OW1) due to statistical alteration of vacancies and H2O molecules at the OW2 site.
Chemical composition
Contents of Na, Ca and S in antofagastaite were determined using a Jeol JSM-6480LV scanning electron microscope equipped with an INCA-Wave 500 wavelength-dispersive spectrometer (Laboratory of Analytical Techniques of High Spatial Resolution, Department of Petrology, Moscow State University), with an acceleration voltage of 20 kV, a beam current of 20 nA, and a beam diameter of 3 µm. Contents of other elements with atomic numbers higher than oxygen are below detection limits. For the holotype, H2O was determined by gas chromatography of products of the mineral ignition at 1200°C with a Vario Micro cube analyser (Elementar GmbH, Germany). Analytical data are given in Table 1. The absence of B, C and N in detectable amounts is confirmed by both crystal structure and IR spectroscopy data.
Chemical data for antofagastaite.

*Averaged for eight spot analyses; **calculated by total difference; S.D. – standard deviation.
The empirical formulae of two samples of the new mineral (based on 8 O atoms belonging to sulfate anions per formula unit (pfu), with all H belonging to H2O molecules – see structure data below) are: Na2.06Ca0.95S2.01O8·1.35H2O for the holotype and Na2.03Ca0.98S2.00O8·1.38H2O for the specimen from Tolbachik. The simplified formula is Na2Ca(SO4)2·1.5H2O, which requires Na2O 20.31, CaO 18.37, SO3 52.46, H2O 8.86, total 100.00 wt.%.
X-ray crystallography
Powder X-ray diffraction data of antofagastaite (Table 2) were collected with a Rigaku R-AXIS Rapid II single-crystal diffractometer equipped with a cylindrical image plate detector (radius 127.4 mm) using Debye-Scherrer geometry, CoKα radiation (rotating anode with VariMAX microfocus optics), 40 kV, 15 mA and exposure 15 min. Angular resolution of the detector is 0.045° in 2θ units (pixel size is 0.1 mm). The data were integrated using the software package Osc2Tab (Britvin et al., Reference Britvin, Dolivo-Dobrovolsky and Krzhizhanovskaya2017). Unit-cell parameters refined from the powder data for the holotype specimen are as follows: a = 6.463(3), b = 6.867(2), c = 9.467(4) Å, β = 104.47(4)° and V = 406.8(4) Å3.
Powder X-ray diffraction data (d in Å) of the holotype specimen of antofagastaite.

* For the calculated pattern, only reflections with intensities ≥ 0.5 are given; ** for the unit-cell parameters calculated from single-crystal data.
The strongest lines are given in bold.
Single-crystal X-ray diffraction studies were carried out using a Bruker D8 Venture diffractometer with a micro-focus X-ray tube operated with MoKα radiation (λ = 0.71073 Å) at 50 kV and 1 mA. More than a half of the Ewald sphere was collected with frame widths of 0.5° in ω, and a 50 s count time for each frame.
The diffraction data were integrated and corrected for absorption using a multi-scan type model using the APEX and SADABS Bruker programs. The unit-cell parameters of antofagastaite [a = 6.4596(4), b = 6.8703(5), c = 9.4685(7) Å, β = 104.580(4)°, V = 406.67(5) Å3and Z = 2 in P21/m] were determined and refined by least-squares techniques on the basis of 2751 reflections. The crystal structure of antofagastaite was solved in space group P21/m by the charge-flipping method, using the SUPERFLIP program (Palatinus and Chapuis, Reference Palatinus and Chapuis2007) and refined to R 1 = 0.057 by means of the SHELX program (Sheldrick, Reference Sheldrick2015). All atoms sites were refined anisotropically. Positions of hydrogen atoms of the fully occupied water OW1 site were taken from the difference-Fourier maps. Hydrogen positions of the one-half occupied water OW2 site were not determined. The main crystallographic information is summarised in Table 3.
Crystallographic data, data collection information and structure refinement details for antofagastaite Na2Ca(SO4)2·1.5H2O.

One relatively large residual electron density peak (2.04 e− Å−3) was observed in the difference-Fourier maps on the final stages of refinement. However, its attribution to an additional O position of a water molecule or to metal cations in the structure of antofagastaite remains ambiguous due to its very low occupancy as refined. Fractional atomic coordinates, thermal displacement parameters, selected interatomic distances, and bond-valence sums for the structure of antofagastaite are listed in Tables 4–5. The crystallographic information file has been deposited with the Principal Editor of Mineralogical Magazine and is available as Supplementary material (see below).
One small (0.03 mm × 0.03 mm × 0.1 mm) prismatic crystal of antofagastaite from the Arsenatnaya fumarole, Tolbachik volcano, was also studied by the means of single-crystal X-ray diffraction. The obtained space group P21/m is the same and unit-cell parameters [а = 6.427(6), b = 6.785(8), c = 9.44(1) Å, β = 105.03(2)° and V = 397.6(8) Å3] are nearly identical to those of the holotype specimen of antofagastaite from the Coronel Manuel Rodríguez mine described above. The crystal structure of anthofagastaite from Tolbachik was solved (R 1 = 7.1%) and proved to be identical to that of the holotype, including the disorder of one of the sulfate groups and H2O molecules.
Fractional atomic coordinates, site occupancy factors (s.o.f.), equivalent/isotropic and anisotropic displacement parameters (U, Å2) of atoms in antofagastaite.

*U iso.
Selected interatomic distances (in Å) and bond-valence sums (BVS, in vu) in the structure of antofagastaite. Distances not considered as bonds are given in italics.

* BVS values divided by two, according to site occupancy factors (= 0.5) of the O5A/B, O6A/B and OW2 sites.
Crystal structure: description and discussion
There are two symmetrically independent S sites (Fig. 6a) in antofagastaite. The S1 site is ordered and demonstrates a regular tetrahedral coordination environment (<S1–O> = 1.470 Å, 6.07 valence units, vu), whereas the S2 site is split over S2A (<S2A–O> = 1.464 Å, 6.18 vu) and S2B (<S2B–O> = 1.467 Å, 6.13 vu) sites with s.o.f. of 0.5 per each. There is a positional disorder also observed for O5 and O6 sites coordinating atom S2.
Coordination of S6+ in the crystal structure of antofagastaite (a). Ball-and-stick and polyhedral representations of the chains formed by CaO7(H2O) polyhedra are shown in (b) and (c), respectively. Ordered S1O4 and disordered S2O4 tetrahedra (yellow) are added in (c). One-half occupied sites are shown by striped balls.

One symmetrically independent Ca site is coordinated by seven oxygen atoms and one water molecule (Fig. 6b). Ca–O distances lie in the range of 2.380–2.555 Å (2.07 vu). CaO7(H2O) polyhedra share common edges, thus, forming the chains depicted in Fig. 6b. Sulfate tetrahedra decorate Ca-centred chains in bidentate and monodentate fashions (Fig. 6c).
The Na site coordination environment is irregular due to the presence of disordered O5 and O6 sites and partially occupied water molecules. Na–O bond lengths vary from 2.242 Å to 2.994 Å.
The structure of antofagastaite consists of ordered and disordered blocks (Fig. 7a). Generally, the structure is similar to that described previously for syngenite K2Ca(SO4)2·H2O (Corazza and Sabelli, Reference Corazza and Sabelli1967; Bokiy et al., Reference Bokiy, Pal'chik and Antipin1978; Ballirano et al., Reference Ballirano, Belardi and Maras2005) (Fig. 7). Antofagastaite and syngenite also demonstrate a close relationship in symmetry, unit-cell parameters, powder X-ray patterns and optical characteristics (Table 6). However, no solid solution between these minerals is observed. The difference between antofagastaite and syngenite is clearly demonstrated by the paragenetic association of them found in the Arsenatnaya fumarole (Fig. 4) in which both minerals are chemically very close to the end-members.
Crystal structures of (a) antofagastaite Na2Ca(SO4)2·1.5H2O; (b) syngenite K2Ca(SO4)2·H2O; and (c) glauberite Na2Ca(SO4)2 projected along the b axis. One-half occupied sites are shown by stripped balls in antofagastaite. The unit cells are outlined.

Comparative data for antofagastaite and syngenite.

*Data for the holotype are given.
The crystal structure of syngenite (Fig. 7b) is characterised by a lower content of water molecules. Intrusion of additional water molecules in the syngenite structure type results in disorder of the one of the two tetrahedral sulfate groups (Fig. 8) observed in antofagastaite (Fig. 7a). Hence, the split into two half-occupied conformations results in a row of sulfate tetrahedra oriented towards each other by their apices along [100] with a water molecule located between them with a 50% probability.
General projection of the arrangement of the (a) ordered S1O4 tetrahedra and OW1 water molecules and (b) disordered S2O4 tetrahedra with OW1 and OW2 water molecules in the crystal structure of antofagastaite. One-half occupied sites are shown by striped balls.

Note also that the antofagastaite crystal structure is significantly different from the chemically closely related but anhydrous mineral glauberite Na2Ca(SO4)2 (Araki and Zoltai, Reference Araki and Zoltai1967) (Fig. 7c).
The crystal structure of antofagastaite contains 1.5 H2O molecules pfu, which are located in one fully occupied OW1 site and in one half-occupied OW2 site. Hydrogen atoms of the OW1 molecule localised from the difference-Fourier electron density map result in OW1–H1 and H1···O4 distances of 0.95(4) Å and 1.87(5) Å, respectively (Fig. 8b). A similar atomic surrounding of the one symmetrically independent water molecule was reported in the syngenite structure: OW–H = 0.82(5) Å and H···O = 1.90(5) Å (Bokiy et al., Reference Bokiy, Pal'chik and Antipin1978).
The presence of a 0.5 vacancy and the occurrence of several possible H-bond acceptor-anions shown in Fig. 9 make the localisation of hydrogen atoms of the half-occupied OW2 water molecule impossible.
Atomic situation surrounding the OW2 site in the crystal structure of antofagastaite. The highest maximum of the residual electron density localised from the difference-Fourier map is represented by a dark green ball. One-half occupied sites are shown by striped balls.

We note the presence of the maximum positive difference-Fourier peak of 2.04 e − Å−3 located at the distances of only 1.96 Å from the OW1 site and of 1.53 Å from its hydrogen atoms (Fig. 9). Attribution of this peak to an O atom of an additional water molecule may provide 0.12 additional H2O. However, we believe that the structure of antofagastaite contains exactly 1.5 H2O, which is constrained to this amount as a function of the neighbouring disordered S2O4 tetrahedra, and, thus, the peak cannot be assigned to a water molecule. As shown in Fig. 10, the residual peak is comparable with the overall noise level within the difference-Fourier map and, thus, it may be ignored in the final structure model of antofagastaite.
Fragment of the antofagastaite structure and isosurfaces corresponding to the maxima of the difference-Fourier electron density shown with a boundary condition of 0.7 e− Å−3 (light green). The highest maximum is localised and displayed by a black spot (boundary condition of 2.0 e− Å−3). Displacement ellipsoids are drawn at the 50% probability level.

Acknowledgements
We thank Gerald Giester, Peter Leverett, František Laufek and anonymous referee for valuable comments. This study was supported by the Russian Science Foundation, grant no. 16-17-10085 (X-ray diffraction and structural studies). The technical support by the SPbSU X-Ray Diffraction Resource Center is acknowledged.
Supplementary material
To view supplementary material for this article, please visit https://doi.org/10.1180/mgm.2019.17.














