<<

Article from the Somma-Vesuvius Complex: New Data and Revised Formula

Taras L. Panikorovskii 1,2,* ID , Nikita V. Chukanov 3, Vyacheslav S. Rusakov 4 ID , Vladimir V. Shilovskikh 5, Anton S. Mazur 6, Giuseppina Balassone 7 ID , Gregory Yu. Ivanyuk 2 and Sergey V. Krivovichev 1,2 1 Department of Crystallography, Institute of Earth Sciences, St. Petersburg State University, University Embankment, 7/9, St. Petersburg 199034, Russia; [email protected] 2 Kola Science Centre, Russian Academy of Sciences, 14 Fersman Street, Apatity 184200, Russia; [email protected] 3 Institute of Problems of Chemical Physics, Russian Academy of Sciences, Chernogolovka 142432, Russia; [email protected] 4 Faculty of Physics, Moscow State University, Vorobievy Gory, Moscow 119991, Russia; [email protected] 5 Department of Colloid Chemistry, Insitute of Chemistry, St. Petersburg State University, University Av. 26, St. Petersburg 198504, Russia; [email protected] 6 Center for Magnetic Resonance, St. Petersburg State University, University Av. 26, St. Petersburg 198504, Russia; [email protected] 7 Dipartimento di Scienze della Terra, dell’Ambiente e delle Risorse, Università degli Studi di Napoli Federico II, Compless Universitario Monte S. Angelo, Via Cinthia 26, 80126 Napoli, Italy; [email protected] * Correspondence: [email protected]; Tel.: +7-812-363-69-14

Received: 2 December 2017; Accepted: 14 December 2017; Published: 16 December 2017

Abstract: At present, the vesuvianite group of minerals consists of eight members, six of which are distinguished by the dominant cation in the Y1(A,B) five-coordinated site. We investigated two vesuvianite samples from the type locality by electron microprobe analysis, Mössbauer and infrared spectroscopy, TGA/DSC, MAS NMR, single- and powder X-ray diffraction. The crystal structures of these samples (# 27844 and 51062 from the Vesuvius collection, Fersman Mineralogical Museum, Moscow) have been refined to R1 = 0.027 and R1 = 0.035, respectively. Both samples have the space group P4/nnc; a = 15.5720(3) and 15.5459(3), c = 11.8158(5) and 11.7988(4), respectively. In both samples low-occupied T1 and T2 sites are populated by minor B and Al, which agrees with their high-temperature origin. According to our experimental results, the general revised crystal-chemical VII−IX V VI formula of vesuvianite can be written as X19 Y1 Y12(Z2O7)4(ZO4)10(W)10, where X are seven- to nine-coordinated sites of Ca with minor Na, K, Fe2+ and REE impurities; VY has a square pyramidal coordination and is occupied predominantly by Fe3+ with subordinate Mg, Al, Fe2+ and Cu2+; VIY has octahedral coordination and is predominantly occupied by Al with subordinate Mg, Fe2+, Fe3+, Mn2+, 3+ Mn , Ti, Cr and Zn; ZO4 = SiO4, sometimes with subordinate AlO4 and/or (OH)4, and W = OH, F, with minor O and Cl. The idealized charge-balanced formula of the vesuvianite end-member without 3+ 2+ 2+ 2+ 2+ subordinate cations is Ca19Fe (Al10Me 2)(Si2O7)4(SiO4)10O(OH)9, where Me = Fe , Mg , Mn .

Keywords: vesuvianite-group minerals; Somma-Vesuvius volcanic complex; vesuvianite; nomenclature;

1. Introduction The vesuvianite-group minerals are widespread in different contact rocks (including skarns formed during contact or regional metamorphism of limestones; in garnetized gabbros, mafic and ultramafic rocks, and serpentinites) of metamorphic, volcanic and hydrothermal origin [1]. They crystallize in

Minerals 2017, 7, 248; doi:10.3390/min7120248 www.mdpi.com/journal/minerals Minerals 2017, 7, 248 2 of 16 a wide range of PT conditions (0–8 kbar and 200–800 ◦C) at the greenschist up to granulite facies of metamorphism and can be considered as a geothermometer [2]. Because of their crystal structure flexibility, vesuvianite-group minerals can contain variable amounts of di- and trivalent cations, and they are stable under reducing and oxidizing conditions [3]. First found by Kappeler [4] at the Somma-Vesuvius volcanic complex, vesuvianite was initially confused with , schorl, and even obsidian [5]. It was established later as a species with its present name by Werner [6]. Many authors tried to determine its formula on the basis of the chemical data [7–10], but only structural investigations made it possible to propose the first rational formula for vesuvianite, Ca10Al4(Mg,Fe)2Si9O34(OH)4 [11]. Another formula provided by Machatschki [12] is X19Y13Z18(O,OH)76, where X = Ca (partly with minor Na, K, Mn), Y = Al, Fe3+, Fe2+, Mg, Ti, Zn, Mn, and Z = Si, which is comparable with the improved formula, Ca19(Mg,Fe,Al,Ti,Mn)5Al8(O,OH)10(SiO4)10(Si2O7)4 [13], determined on the basis of accurate single-crystal XRD studies. incorporation into additional T sites (T1 with tetrahedral coordination and T2 with triangular coordination) of the vesuvianite structure is reflected in the following formula: 3+ 2+ 3+ 3+ 2+ 4+ X19Y13Z18T0–5O68W10 [14], where X = Ca, Na, REE , Pb and Sb ; Y = Al, Mg, Fe , Fe , Ti , Mn, Cu and Zn; Z = Si; T = B; W = (OH, F, O). Over the past 50 years, numerous substitution schemes have been established that significantly expand the crystal-chemical diversity of the vesuvianite group:

− − − T1B3+ + YMg2+ + 2O10,11O2 ↔T1 + YAl3+ + 2O10,11(OH) and T2B3+ + 2O10,11O2 ↔T2 + O10,11(OH)− in [15] Y1Mn3+↔Y1Fe3+ in manganvesuvianite [16] Z1,2 4− Z1,2 4− (SiO4) ↔ (H4O4) [17] O10OH− + O11OH−↔O10F− + O11F− in fluorvesuvianite [18] O10OH−↔O10Cl− [19] X3Bi3+ + Y2,3Mg2+↔X3Ca2+ + Y2,3Al3+ [20] Y2Mg2+ + Y3Ti4+↔Y2Al3+ + Y3Al3+ [21] X1,4Ca2+ + Y2,3Al3+↔X1,4Na+ + Y2,3Ti4+ in “natrovesuvianite” [22] Y1Cu2+ + Y2,3Mn3+↔Y1Fe3+ + Y2,3Mg2+ in cyprine [23] Y1Mg2+ + Y3Al3+↔Y1Fe3+ + Y3Mg2+ in magnesiovesuvianite [24] Y1Fe3+↔Y1Al3+ in alumovesuvianite [25]

At present, the vesuvianite group consists of seven mineral species (Table1), distinguished by the dominant component at five-coordinated Y1 site, as well as T1 and T2 and anionic W positions. Taking into account cation ordering in the octahedral positions and the incorporation of additional 4− B, Na and (H4O4) into the structure, the general formula of the vesuvianite-group minerals can be written as follows (Z = 2): X16X12X4Y1Y24Y38T0–5(Z2O7)4[(ZO4)10−x(H4O4)x](W)9O1–3, where x < 3, X are seven- to nine-coordinated sites (Ca, Na, K, Fe2+, REE), X4 has square antiprismatic coordination (Ca, Na), Y1 has square pyramidal coordination (Fe3+, Mg, Al, Fe2+, Cu2+), Y2 and Y3 have octahedral coordination (Al, Mg, Zn, Fe2+, Fe3+, Mn2+, Mn3+, Ti, Cr, Zn), T is the additional site with triangular or tetrahedral coordination (B, Fe), ZO4 [SiO4, , (OH)4, AlO4] and Z2O7 (Si2O7) and the additional anionic position W can be occupied by OH, F, or minor O, Cl [12,13,15,17,19,21,26–29]. The crystal structure of vesuvianite-group minerals contains half-populated cation sites arranged along the fourfold axis in the Y1–X4–X4–Y1 sequence with the Y1–X4 and X4–X4 distances less than 1.3 Å and 2.5 Å, respectively [28]. Cation ordering at the Y1 and X4 sites produces different ordering schemes and leads to the lower symmetry P4/nnc → P4/n and P4/nnc → P4nc [13,30], which results in the appearance of subsites (for example Y1A and Y1B instead Y1, X4A and X4B instead X4, etc.). Minerals 2017, 7, 248 3 of 16

Table 1. Dominant components in crystallographic sites of vesuvianite-group minerals 1.

Mineral/Formula X1, X2, X3 X4 Y1 Y2 Y3 T1 T2 O10 O11 O12 Reference Vesuvianite s.s. Ca Ca Fe3+ Al Al   OH OH  [30] Fluorvesuvianite Ca Ca Fe2+ Al Al   FF  [18] Manganvesuvianite Ca Ca Mn3+ Al Al   OH OH  [16] Cyprine Ca Ca Cu Al Al   OH OH  [23] Magnesiovesuvianite Ca Ca Mg Al Al   OH OH  [24] Alumovesuvianite Ca Ca Al Al Al   OH OH  [25] Wiluite Ca Ca Mg Al Al B B O O O [14,29] Hongheite Ca Ca Fe2+ Al Fe3+  BOO  [31] 1 Vesuvianite-group minerals are characterized by domain structure, and their symmetry is connected with the kind of cation ordering along the four-fold axis [space groups are either tetragonal (P4/nnc, P4/n, P4nc, or P-4) or monoclinic (P2/n or Pn)] [12,32–34]. Some contain growth sectors showing triclinic distortion [35].

However, the general formula of vesuvianite given in the International Mineralogical Association (IMA) (http://ima-cnmnc.nrm.se/imalist.htm), (Ca,Na)19(Al,Mg,Fe)13(SiO4)10(Si2O7)4(OH,F,O)10 [36], does not explicitly show the occupancies of the Y1 site. This simplified formula is now out of date since it refers to five vesuvianite-group minerals: vesuvianite sensu stricto, manganvesuvianite, cyprine, magnesiovesuvianite and alumovesuvianite. For this reason, we report the results of a multimethodological study of two vesuvianite samples from the type locality of this mineral, i.e., the Somma-Vesuvius complex, Campania, Italy. The studied samples nos. 27844 and 51062 originate from the Vesuvius collection of the Fersman Mineralogical Museum, Moscow. The aim of this work is to establish a correct formula of vesuvianite, which does not contradict the IMA Commission on New Minerals, Nomenclature and Classification (CNMNC) guidelines [37].

2. Occurrence Vesuvianite from the Somma-Vesuvius volcanic complex (Campania region, Southern Italy), typically occurs in various kinds of xenolithic rocks, commonly associated with various Plinian eruptives and found as ejecta in pyroclastic deposits [38–41]. The typical occurrence of vesuvianite is related to xenoliths of skarn, a largely mineralogically zoned, coarse-grained Ca-Mg-Fe(Mn)-rich silicate rocks [42]. At Somma-Vesuvius, these rocks are very complex in mineralogical and textural aspects. Most of skarns at Somma-Vesuvius complex show a sequence of mono- or bimineral zones dominated by Ca- and/or Mg-silicates, such as clinopyroxene, phlogopite, , , olivine, or , and other minerals (perovskite, spinel, calcite, etc.) commonly with a central cavity containing euhedral crystals of vesuvianite, , wollastonite, gehlenite, or anorthite. Mono- or bimineral skarns, as well as composite ejecta (consisting of two or more rock types, such as skarn-marble, skarn-, cumulate-skarn or cumulate skarn-marble are common and document the close spatial association of at least some of these rock types at depth) with sharp contacts of skarns to marbles or are common [42–44]. The mineralogical composition of skarn ejecta generally includes highly variable amounts of vesuvianite, wollastonite, clinopyroxene (diopside, hedenbergite, low- augite, i.e., “fassaite”), anorthite, phlogopite, clinoamphibole, garnet (mainly of the - series), forsterite, and humite-group minerals. Other minerals locally observed as minor to accessory phases include melilite (gehlenite), feldspathoids (leucite, nepheline, sodalite, haüyne, nosean, , meionite, davyne, balliranoite), fluorapatite, fluorite, cuspidine, zircon, perovskite, calcite, baddeleyite, spinels (spinel sensu stricto, magnetite), REE-minerals, sulfides, etc. [41,42]. The vesuvianite-bearing calc-silicate ejecta are most likely related to the Plinian eruptions of Avellino (3800 BC) and Pompeii (79 AD) [44]. The investigations of fluid inclusions as well as the isotope compositions of Pb, Nd and Sr have been carried out, among others, on vesuvianite crystals as reported in [42]. According to [42], the formation of skarns can be related to assimilation of carbonate-wall rocks by the alkaline magma at moderate depths (<5 km), and consequent exsolution of CO2-rich vapor and complex saline melts from the contaminated magma that reacted with the carbonate rocks to form Minerals 2017, 7, 248 4 of 16 Minerals 2017, 7, 248 4 of 16 carbonate rocks to form skarns. According to fluid inclusion studies [42], the estimated temperature skarns.of the formation According skarn to fluid minerals inclusion at Somma-Vesuvius studies [42], the complex estimated is in temperature the range 800–900 of the formation °C. skarn mineralsTogether at Somma-Vesuvius with the skarn complex occurrence, is in the vesuvian range 800–900ite has◦C. been also observed in thermally metamorphosedTogether with carbonate the skarn nodules occurrence, (“marbles” vesuvianite [5]), has in beenterpreted also observed as the result in thermally of recrystallization metamorphosed and carbonatein part, devolatilization nodules (“marbles” without [5]), interpretedsignificant asmetaso the resultmatic of interaction recrystallization with andaqueous in part, fluids devolatilization or silicate withoutmelts [42]. significant In these metasomatic calcitic marbles, interaction vesuvianite with aqueous is variably fluids or associated silicate melts with [42 ].(i) In wollastonite, these calcitic marbles,clinopyroxene, vesuvianite sodalite is variablyand humite; associated (ii) clinopyroxen with (i) wollastonite,e, sodalite, clinopyroxene,cuspidine, forsterite, sodalite nepheline and humite; and (ii)davyne; clinopyroxene, or (iii) sanidine, sodalite, fluorite, cuspidine, fluorapatite forsterite, and sulfides nepheline [5]. and davyne; or (iii) sanidine, fluorite, fluorapatiteIn addition, and sulfides vesuvianite [5]. occurs as an accessory mineral in syenitic K-feldspar-rich nodulesIn addition, [5,39]. In vesuvianite this rock, it occurs is locally as an observed accessory in mineral vugs and in syeniticcavities K-feldspar-richas grains (<1 mm) nodules among [5 ,39the]. Intabular this rock, sanidine it is locallycrystals, observed in association in vugs with and garnet, cavities clinoamphibole, as grains (<1 mm) clinopyroxene, among the tabulardavyne, sanidine titanite, crystals,or with influorite, association zircon with and garnet, wöhlerite-group clinoamphibole, minerals clinopyroxene, [5]. In [39], davyne, accessory titanite, vesuvianite or with fluorite,crystals zircon(~200–500 and µm) wöhlerite-group were observed minerals in the [K-feldspar5]. In [39], accessorysyenitic ejecta vesuvianite together crystals with other (~200–500 minorµ m)to weretrace observedphases that in can the be K-feldspar alternatively syenitic represented ejecta together by clinoamphibole, with other minorclinopyroxene, to trace phasesbiotite, plagioclase, that can be alternativelysodalite, nosean, represented nepheline, by hematite, clinoamphibole, magnetite, clinopyroxene, fluorite and biotite, titanite. plagioclase, Vesuvianite sodalite, also occurs nosean, in nepheline,sanidinite hematite,ejecta found magnetite, at the fluoriteSan Sebastiano and titanite. valley Vesuvianite at Vesuvius also in occurs association in sanidinite with ejectapyrochlore, found atmagnetite, the San Sebastianowöhlerite group valley minerals, at Vesuvius baddeleyite, in association zircon, with fluorite, pyrochlore, ferrohornblende, magnetite, wöhlerite phlogopite group and minerals,nepheline baddeleyite,[45]. zircon, fluorite, ferrohornblende, phlogopite and nepheline [45]. In samples 51062 and 27844 (Figure 11),), euhedraleuhedral vesuvianitevesuvianite occursoccurs inin associationassociation withwith grossulargrossular and clinozoisite or as a monomineralic aggregate, respectively;respectively; the large sizes of the crystals indicate their skarnskarn provenance. provenance. In In particular, particular, the associationthe association of vesuvianite of vesuvianite with garnet, with wollastonite,garnet, wollastonite, gehlenite orgehlenite anorthite or inanorthite the first in sample the first is typical sample for is zonedtypical skarns for zo ofned the skarns Somma-Vesuvius of the Somma-Vesuvius complex [42]. Itcomplex is worth [42]. nothing It is worth that the nothing occurrence that the of clinozoisiteoccurrence of in theclinozoisite sample 51062in therepresents, sample 51062 tothe represents, authors’ knowledge,to the authors’ the knowledge, first record ofthe this first mineral record atof Somma-Vesuvius.this mineral at Somma-Vesuvius.

Figure 1. Short-prismaticShort-prismatic vesuvianite vesuvianite crystals crystals (1) (1) associated associated with with grossular grossular (2), (2), and and clinozoisite clinozoisite (3), (3), in insamples samples (A ()A 51062) 51062 andand (B) ( B27844) 27844. .

3. Materials and Methods Both vesuvianitevesuvianite samples samples studied studied in thisin this work work (27844 (27844and 51062 and) 51062 form well-shaped) form well-shaped short-prismatic short- crystalsprismatic up crystals to 6 mm up long to (Figure6 mm 1long) with (Figure the main 1) prismaticwith the formsmain {100}prismatic and {110}forms terminated {100} and by {110} the pyramidalterminated {101} by the and pyramidal pinacoidal {101} {001} and faces. pinacoidal The crystals {001} are faces. lustrous The and crystals semi-transparent. are lustrous Theirand semi- color variestransparent. from greenish Their color brown varies (27844 from) to greenish brownish brown green (27844 and orange-brown) to brownish ingreen the marginal and orange-brown zone (51062 in). the marginalChemical zone analyses (51062 of). both samples (Table2, each analysis is the mean from 4 to 7 points) were obtainedChemical using analyses a HITACHI of both S-3400N samples scanning (Table electron 2, each microscope analysis is (Tokyo, the mean Japan) from equipped 4 to 7 points) with INCA were Waveobtained 500 WDSusing spectrometer.a HITACHI S-3400N The system scanning was operated electron at 20microscope kV and 10 (Tokyo nA, and, Japan) the electron equipped beam with was focusedINCA Wave to a 5500µm WDS spot. spectrometer. The following The standards system were was used:operated rutile at (Ti),20 kV wollastonite and 10 nA, (Si, and Ca), the forsterite electron (Mg),beam fluoritewas focused (F), cryolite to a 5 (Na),µm spot. halite The (Cl), following spessartine stan (Mn,dards Al), were hematite used: (Fe),rutile danburite (Ti), wollastonite (B), synthetic (Si, Ce-glassCa), forsterite (Ce),synthetic (Mg), fluorite Nd-glass (F), (Nd). cryolite The content(Na), halite of boron (Cl), was spessartine measured (Mn, at the Al), Faculty hematite of Geology, (Fe), Universitydanburite (B), of Warsaw,synthetic byCe-glass a CAMECA (Ce), synthetic SX100 instrument Nd-glass (Nd). operated The atcontent 6 kV andof boron 80 nA was for measured 50–100 s countingat the Faculty at each of pointGeology, (Kα line,University PC2 crystal of Warsaw, (Ni/C-LSM), by a CAMECA standard deviationSX100 instrument 0.25 wt % operated B2O3). No at other6 kV elementsand 80 nA having for 50–100 atomic s counting numbers at above each point eight (K wereα line, found PC2 above crystal their (Ni/C-LSM), detection standard limits. The deviation content 0.25 of Hwt2 O% inB227844O3). Nowas other determined elements byhaving means atomic of thermogravimetric numbers above eight analysis were using found NETZSCH above their STA detection 449 F3

Minerals 2017, 7, 248 5 of 16

Jupiter thermoanalyzer (Exton, PA, USA) during 10 ◦C/min heating of a 68.4 mg sample from room temperature to 1100 ◦C in a dynamic argon atmosphere, with aluminum oxide standard.

Table 2. Single-crystal data, data collection and structure refinement parameters.

Sample 27844 51062 Temperature/K 293(2) 293(3) Tetragonal Tetragonal Space group P4/nnc P4/nnc a = b (Å) 15.5720(3) 15.5459(3) c (Å) 11.8158(5) 11.7988(4) Volume (Å3) 2865.16(16) 2851.48(16) Z 2 2 −3 ρcalc (g/cm ) 3.395 3.405 µ (mm−1) 2.839 2.837 F(000) 2941.0 2901.0 Crystal size (mm3) 0.21 × 0.15 × 0.14 0.22 × 0.19 × 0.17 Radiation MoKα (λ = 0.71073) MoKα (λ = 0.71073) 2θ range for data collection (◦) 5.232–54.982 6.802–54.974 −20 ≤ h ≤ 19 −13 ≤ h ≤ 20 Index ranges −15 ≤ k ≤ 20 −14 ≤ k ≤ 15 −15 ≤ l ≤ 14 −15 ≤ l ≤ 7 Reflections collected 14,428 5687 1653 [R = 0.0440 1622 [R = 0.0274 Independent reflections int int Rsigma = 0.0178] Rsigma = 0.0267] Data/restraints/parameters 1653/0/162 1622/0/162 Goodness of fit on F2 1.184 1.114 R = 0.027 R = 0.030 Final R indices [I > 2σ (I)] 1 1 wR2 = 0.072 wR2 = 0.078 R = 0.028 R = 0.035 Final R indices [all data] 1 1 wR2 = 0.072 wR2 = 0.080 Largest diff. peak/hole (e Å−3) 0.56/−1.05 0.53/−1.08

The 27Al nuclear magnetic resonance (NMR) spectrum was obtained at room temperature by means of a Bruker Avance III 400 WB spectrometer (Billerica, MA, USA) at 104.24 MHz. The magic- angle spinning (MAS) spin rate of the rotor was 20 kHz. For these investigations, a single-pulse sequence with a pulse length of 86 kHz was used, with the recycle delay of 1 s and 4096 scans. The 1H NMR spectrum was obtained at room temperature using a Bruker Advance III 400 WB spectrometer operating at 400.23 MHz. The MAS spin rate of rotor was 20 kHz. A single-pulse sequence was used with a pulse length of 100 kHz, recycle delay of 20 s and 32 scans. The Fe2+/Fe3+ ratio is given in accordance with the Mössbauer spectrometry data. Gamma-resonance studies of 27844 and 51062 were carried out using a commercial Mössbauer spectrometer (WissEl Wissenschaftliche Elektronik GmbH, Starnberg, Germany) in order to clarify the valence state of iron and the coordination environment(s) of the Fe. The spectrometer was equipped with a standard radioactive source, 57Co, in a rhodium matrix (Ritverc). The absorption-mode measurements have been performed at room temperature. The samples were crushed in an agate mortar, distributed homogeneously within the punch-holder and packed in the form of a tablet with a density of 112 mg/cm2 in order to optimize the conditions for measurement. The spectral fitting of the experimental Mössbauer spectrum was carried out by the method of least squares under the assumption of Lorentz-shaped spectral lines. Isomer shifts were measured relative to metallic iron. The mathematical analysis was performed with the commercial Mössbauer software SpectRelax Version 2.1 (Moscow State University, Moscow, Russia). Minerals 2017, 7, 248 6 of 16

A FTIR spectrum of 27844 was collected using a Bruker ALPHA spectrometer at room temperature in the range 360–3800 cm−1 (KBr pellet, resolution 4 cm−1, 16 scans). A pure KBr pellet was used as a reference. Powder X-ray diffraction data (Table A1) were collected with a Bruker Phaser D2 diffractometer in the 2θ range of 10–80◦ (CuKα; 1.5418 Å), with a scanning step of 0.02◦ (1 s per step) in 2θ. Single-crystal X-ray diffraction experiments were carried out using an Agilent Technologies Xcalibur Eos diffractometer (Agilent Technologies, Santa Clara, CA, USA) operated at 50 kV and 40 mA. A hemisphere of three-dimensional data was collected at 293 K for each sample using monochromatic Mo Kα radiation, with frame widths of 1◦ and 20 s count for each frame. The crystal-to-detector distance was 45 mm. For the investigated samples, P4/nnc was chosen as the most probable space group. Only 79 reflections violating the absence conditions in the range 2σ(I) < I < 188σ(I) with 2 maximum Fo = 881 were found for 27844, whereas only 15 such reflections with 2σ(I) < I < 382σ(I) 2 and maximum Fo = 2021 were observed for 51062. Refinements of both structures in space groups P4/n and P4nc did not indicate any reasons for the symmetry lowering from the P4/nnc space group. The crystal structures of 27844 and 51062 were refined to R1 = 0.027 and R1 = 0.035 for 1653 and 1622 unique observed reflections with |Fo| ≥ 4σF, respectively, using the SHELX program [46]. Empirical absorption corrections were applied in the CrysAlisPro program complex [47] using spherical harmonics, implemented in the SCALE3 ABSPACK scaling algorithm. The final difference Fourier map showed no features higher than 0.53/0.56 e− for 27844 and 51062, respectively. The experimental details and crystallographic parameters are given in Table2. Crystallographic information files (CIFs) are deposited in the electronic Supplementary Materials.

4. Results

4.1. Chemical Composition The crystal-chemical formulae of the investigated samples (Table3) were calculated on the basis of 19 X-cations (Ca + Na + REE). Both compositions are in good agreement with the previous data for vesuvianite from the Monte-Somma-Vesuvius complex [38,48].

Table 3. Chemical composition of 27844 and 51062 1.

wt % apfu Component Component 27844 51062 27844 51062

SiO2 36.41 35.89 Si 17.81 17.99 Al2O3 15.82 15.93 Al 9.12 9.41 TiO2 bd 0.41 Ti bd 0.16 1 3+ Fe2O3 3.48 2.97 Fe 1.28 1.12 FeO 1 0.98 1.29 Fe2+ 0.40 0.54 MnO 0.27 0.38 Mn 0.11 0.16 MgO 2.51 3.21 Mg 1.83 2.40 CaO 36.19 35.19 Ca 18.97 18.90 Ce2O3 bd 0.33 Ce bd 0.06 Nd2O3 bd 0.22 Nd bd 0.04 Na2O 0.03 bd Na 0.03 bd 2 B2O3 0.83 0.79 B 0.70 0.68 Cl 0.06 0.05 Cl 0.05 0.04 F 1.29 1.22 F 2.00 1.93 3 H2O 1.60 1.76 OH 5.22 5.88 –O=F,Cl −0.56 −0.52 O 73.04 74.75 Total 98.90 98.57 1 Although Fe was measured as Fe2O3, it is reported as Fe2O3 and FeO based on the results from Mössbauer spectroscopy data; 2 Boron content measured using a CAMECA SX-100 instrument at the same grains (PC2 crystal, 3 6 kV and 80 nA for 50–100 s at each point, mean of 7 point analyses); The content of H2O measured by thermogravimetry (TGA) method for 27844 and calculated from the structural data for 51062. Minerals 2017, 7, 248 7 of 16

Minerals 2017, 7, 248 7 of 16 4.2.Minerals Thermogravimetric 2017, 7, 248 Analysis and Differential Scanning Calorimetry (TGA/DSC) 7 of 16 4.2. Thermogravimetric Analysis and Differential Scanning Calorimetry (TGA/DSC) The TGA curve (Figure2) of 27844 contains one major step of weight loss, accompanied by 4.2. ThermogravimetricThe TGA curve (Figure Analysis 2) andof 27844 Differential contains Scanning one major Calorimetry step of weight(TGA/DSC) loss, accompanied by two two endothermic effects at 980 ◦C and 1030 ◦C in the DSC-curve, which can be assigned to separated endothermic effects at 980 °C and 1030 °C in the DSC-curve, which can be assigned to separated dehydroxylationThe TGA curve processes (Figure for 2) the of 27844 O(11)H contains and O(10)H one major groups step [ 23of ,weight49]. According loss, accompanied to [49,50] by one-step two dehydroxylationendothermic effects processes at 980 for°C theand O(11)H 1030 °C and in O(10)Hthe DSC-curve, groups [23,49].which canAccording be assigned to [49,50] to separated one-step weight loss in the temperature range of 820–1090 ◦C is characteristic for low-temperature vesuvianite. weightdehydroxylation loss in the processestemperature for rangethe O(11)H of 820–1090 and O(10)H °C is characteristic groups [23,49]. for low-temperatureAccording to [49,50] vesuvianite. one-step Weight loss in the range of 990–1150 ◦C responds to the high-symmetry vesuvianite. The total weight Weightweight lossloss inin thethe temperature range of 990–1150 range °Cof 820–1090responds °C to isthe characteristic high-symmetry for low-temperature vesuvianite. The vesuvianite.total weight loss for 27844 is 1.60 wt %, which corresponds to 5.22 OH-groups per formula unit. As a result lossWeight for loss27844 in isthe 1.60 range wt of%, 990–1150 which corresponds °C responds to to 5.22 the high-symmetryOH-groups per vesuvianite. formula unit. The As total a result weight of of dehydroxylation, vesuvianite transforms into the mixture of wollastonite, grossular, quartz and dehydroxylation,loss for 27844 is 1.60 vesuvianite wt %, which transforms corresponds into the to mi 5.22xture OH-groups of wollastonite, per formula grossular, unit. quartzAs a result and anof anamorphousdehydroxylation, amorphous glass-like glass-like vesuvianite substance. substance. transforms into the mixture of wollastonite, grossular, quartz and an amorphous glass-like substance.

FigureFigure 2. Thermogravimetric2. Thermogravimetric (wt (wt %, green%, green line) line) and and differential differential scanning scanning calorimetry calorimetry (blue (blue line) line) curves forcurvesFigure27844 .for2. Thermogravimetric27844. (wt %, green line) and differential scanning calorimetry (blue line) curves for 27844. 4.3. Solid State Magic-Angle Spinning Nuclear Magnetic Resonance (MAS NMR) 4.3. Solid State Magic-Angle Spinning Nuclear Magnetic Resonance (MAS NMR) 4.3. SolidThe 27StateAl MAS Magic-Angle NMR spectrum Spinning (Figure Nuclear 3) Magnetic of 27844 Resonance contains one(MAS broad NMR) peak, which is centered at The 27Al MAS NMR spectrum (Figure3) of 27844 contains one broad peak, which is centered −2.5 ppmThe 27 andAl MAS can NMRbe assigned spectrum to the(Figure octahedrally 3) of 27844 coordinated contains one Al broad at the peak, Y2 and which Y3 issites centered [41–52]. at at −2.5 ppm and can be assigned to the octahedrally coordinated Al at the Y2 and Y3 sites [41–52]. Similarly−2.5 ppm toand wiluite, can be broadening assigned toof the peakoctahedrally at 60 ppm coordinated can be assigned Al at tothe a traceY2 and amount Y3 sites of Al [41–52]. in the Similarly to wiluite, broadening of the peak at 60 ppm can be assigned to a trace amount of Al in the 4-coordinatedSimilarly to wiluite, T1 site broadening (<1 apfu) [29]. of the The peak weak at 60symmetric ppm can peaksbe assigned observed to a neartrace ±190 amount ppm of belong Al in the to 4-coordinated T1 site (<1 apfu)[29]. The weak symmetric peaks observed near ±190 ppm belong to spinning4-coordinated sidebands. T1 site (<1 apfu) [29]. The weak symmetric peaks observed near ±190 ppm belong to spinning sidebands.1 spinningThe sidebands.H MAS NMR spectrum (Figure 3) of 27844 contains two peaks centered at 6.9 and 1.5 ppm, The 1H MAS NMR spectrum (Figure3) of 27844 contains two peaks centered at 6.9 and 1.5 ppm, respectively.The 1H MAS Calculations NMR spectrum of the O···O(Figure distances 3) of 27844 using contains the equation two peaks δiso centered(ppm) = at79.05 6.9 and− 0.255 1.5 dppm, (O‒ respectively. Calculations of the O··· O distances using the equation δ (ppm) = 79.05 − 0.255d H···O)(pm)respectively. [53] Calculations provided ofthe the values O···O ofdistances 2.83 Å usingand 3.04 the Å.equation The most δiso(ppm) intenseiso = 79.05peak −at 0.255 6.9 dppm (O‒ (O–HcorrespondsH···O)(pm)··· O)(pm) [53] to [the 53provided] O11 provided‒H1a···O7 the thevalues and values O10of 2.83 of‒H2···O10 2.83 Å and Å andhydrogen 3.04 3.04 Å. Å.The interactions, The most most intense intensewhich peak is peak in at agreement at6.9 6.9 ppm ppm correspondswithcorresponds the distances to to the the O11–H1a 2.784(4) O11‒H1a···O7 ···andO7 2.74–2.82 andand O10–H2O10 Å‒ reportedH2···O10··· O10 byhydrogen other authors interactions, interactions, [54–59]. which A which low is intense isin inagreement agreement line at with1.5with theppm the distances candistances be assigned 2.784(4) 2.784(4) to and and the 2.74–2.82 2.74–2.82H1b site, ÅsinceÅ reportedreported its chemical by other shift authors authors agrees [54–59]. [well54–59 with]. A A low the low intenseO11 intense‒H1···O11 line line at at 1.5distances1.5 ppm ppm can can of be 3.01–3.04 be assigned assigned Å to[54,57]. to the the H1b H1b site, site, since since its its chemical chemical shiftshift agreesagrees wellwell with the O11 O11–H1‒H1···O11··· O11 distancesdistances of of 3.01–3.04 3.01–3.04 Å Å [54 [54,57].,57].

Figure 3. (A) 1H MAS NMR and (B) 27Al MAS NMR spectra of 27844. FigureFigure 3. 3.(A (A) 1) H1H MAS MAS NMRNMR andand ((B)) 27AlAl MAS MAS NMR NMR spectra spectra of of 27844.27844.

Minerals 2017, 7, 248 8 of 16

Minerals 2017, 7, 248 8 of 16 4.4. Mössbauer Measurenments 4.4. Mössbauer Measurenments The Mössbauer spectra of 27844 and 51062 can be described as a superposition of five and four symmetricThe Mössbauer doublets, respectively spectra of 27844 (Figure and4). 51062 The assignments can be described based onas a the supe crystal-structurerposition of five refinement and four andsymmetric Mössbauer doublets, parameters respectively of these doublets(Figure (isomer4). The shiftsassignmentsIS, quadrupole based splittingson the crystal-structureQS, line widths andrefinement areas under and Mössbauer resonance parameters doublets S) of according these doublets to [21 ,(isomer60,61] are shifts listed IS, inquadrupole Table4. Trivalent splittings iron QS, dominatesline widths in and both areas samples: underC doublets resonance with doubletsIS = 0.38, S 0.38,) accordingQS = 0.49, to 0.52[21,60,61] and S =are 54, listed 60% are in assignedTable 4. toTrivalent Fe3+ at ironY3 octahedra, dominates whereas in both samples:E doublet C withdoubletsIS = with 0.48, IS 0.41, = 0.38,QS 0.38,= 1.08, QS 1.31 = 0.49, and 0.52S = and 12.8, S = 16.3 54, corresponds60% are assigned to Fe 3+to Feat3+ the at five-coordinatedY3 octahedra, whereasY1 site. E doubletD doublets with with IS = IS0.48,= 0.82,0.41, 0.84,QS = QS1.08,= 1.31 0.40, and 0.43 S and= 12.8,S = 16.3 18, 15% corresponds correspond to toFe divalent3+ at the ironfive-coordinated at the octahedral Y1 Ysite.3 site. D doublets The B doublets with IS with = 0.82,IS = 1.02,0.84, 1.03,QS =QS 0.40,= 2.92,0.43 and 2.92 S and = 18,S = 15% 18, 15%,correspond similar to to divalent doublets iron reported at the inoctahedral [21,61–63 ],Y3 correspond site. The B todoublets minor amountswith IS = of 1.02, Fe2+ 1.03,at the QSY1 site.= 2.92, The 2.92G doublet and S with= 18,IS 15%,= 0.92 similar and QS to =doublets 2.40 observed reported in the in Mössbauer[21,61–63], spectracorrespond of 27844 to minorand wiluite amounts [29 ]of can Fe be2+ at assigned the Y1 tosite. the The octahedrally G doublet coordinated with IS = 0.92 Fe2+ andincorporated QS = 2.40 intoobserved the Y 2in site. the Mössbauer spectra of 27844 and wiluite [29] can be assigned to the octahedrally coordinated Fe2+ incorporated into the Y2 site. Table 4. Parameters of the Mössbauer spectra of 27844 and 51062 at 293 K. Table 4. Parameters of the Mössbauer spectra of 27844 and 51062 at 293 K. Designation of the Isomer Shift Quadrupole Line Width Relative Assignment 1 QuadrupoleDesignation Doublet of the (mm/s)Isomer Shift SplittingQuadrupole (mm/s) Splitting (mm/s)Line Width AreaRelative (%) Area Assignment 1 QuadrupoleSample Doublet (mm/s) (mm/s) 27844 (mm/s) (%) Sample 27844 VI 3+ C 0.3770.377± 0.002 ± 0.002 0.492 ± 0.4920.004 ± 0.004 54.4 ±54.40.6 ± 0.6 Fe VIinFeY3+ 3in Y3 V 3+ E 0.4810.481± 0.006 ± 0.006 1.076 ± 1.0760.011 ± 0.011 12.8 ± 12.80.7 ± 0.7 Fe VinFeY3+ 1in Y1 VI 2+ D 0.8210.821± 0.003 ± 0.003 0.396 ± 0.3960.004 ± 0.004 0.460 0.460± 0.004 ± 0.004 17.9 ± 17.90.5 ± 0.5 Fe VIinFeY2+ 3in Y3 V 2+ B 1.0391.039± 0.012 ± 0.012 2.920 ± 2.9200.031 ± 0.031 6.1 ± 6.10.5 ± 0.5 Fe VinFeY2+ 1in Y1 VI 2+ G 0.9210.921± 0.009 ± 0.009 2.400 ± 2.4000.022 ± 0.022 8.8 ± 8.80.5 ± 0.5 Fe VIinFeY2+ 2in Y2 Sample 5106251062 C 0.377 ± 0.007 0.528 ± 0.014 59.9 ± 0.4 VIFe3+ in Y3 C 0.377 ± 0.007 0.528 ± 0.014 59.9 ± 0.4 VIFe3+ in Y3 E 0.405 ± 0.004 1.312 ± 0.010 16.3 ± 0.5 VFe3+ in Y1 E 0.405 ± 0.004 1.312 ± 0.0100.432 ± 0.004 16.3 ± 0.5 VFe3+ in Y1 D 0.842 ± 0.030 0.438 ± 0.060 0.432 ± 0.004 14.8 ± 0.4 VIFe2+ in Y3 D 0.842 ± 0.030 0.438 ± 0.060 14.8 ± 0.4 VIFe2+ in Y3 B 1.030 ± 0.007 2.920 ± 0.014 9.0 ± 0.3 VFe2+ in Y1 B 1.030 ± 0.007 2.920 ± 0.014 9.0 ± 0.3 VFe2+ in Y1 1 The symbols and the assignment of the quadrupole doublets are given in accordance with the data 1 The symbols and the assignment of the quadrupole doublets are given in accordance with the data from [21,29,61] basedfrom on[21,29,61] the analysis based of Mössbauer on the analysis spectra of 17Mössbauer vesuvianite spectra samples of from 17 vesuvianite different localities. samples Taking from into different account Fe–Olocalities. distances Taking and site into multiplicities, account Fe–O quadrupole distances doublet andd withsite anomalouslymultiplicities, low quadrupole quadrupole splitting doublet could d with not Y beanomalously assigned to Fe low at distortedquadrupole1 sites. splitting could not be assigned to Fe at distorted Y1 sites.

FigureFigure 4.4. MössbauerMössbauer spectraspectra ofof ((AA)) 2784427844 andand ((BB)) 5106251062..

Minerals 2017, 7, 248 9 of 16 Minerals 2017, 7, 248 9 of 16 Minerals 2017, 7, 248 9 of 16 4.5.4.5. X-ray X-ray Crystallography Crystallography 4.5. X-ray Crystallography TheThe unit-cell unit-cell dimensions dimensions ofof27844 27844 werewere determineddetermined from from the the X-ray X-ray powder-diffraction powder-diffraction pattern pattern (Figure(Figure5The) by5) unit-cellby Rietveld Rietveld dimensions refinement refinement of using using27844 the thewere program program determined TopasTopas from 4.24.2 [[64] the64] and X-ray and are are powder-diffraction in in a agood good agreement agreement pattern with with thethe(Figure single-crystal single-crystal 5) by Rietveld XRD XRD data: refinementdata:a a= = 15.5779(6)15.5779(6) using the Å, program c = 11.8156(6) 11.8156(6) Topas Å, 4.2 Å, V V[64] = =2867.3(4) 2867.3(4)and are Åin3Å .a 3good. agreement with the single-crystal XRD data: a = 15.5779(6) Å, c = 11.8156(6) Å, V = 2867.3(4) Å3.

FigureFigure 5. 5.Powder Powder diffraction diffraction patternpattern ofof 27844.. All All reflections reflections refer refer to to vesuvianite vesuvianite (Table (Table A1). A1 ). Figure 5. Powder diffraction pattern of 27844. All reflections refer to vesuvianite (Table A1). Structural models of 27844 and 51062 are close to that reported in [48], with site designations Structural models of 27844 and 51062 are close to that reported in [48], with site designations accordingStructural to [15]. models We pay of particular 27844 and attention 51062 are to closethe occupancy to that reported of the Yin1 site[48], that with plays site adesignations key role in according to [15]. We pay particular attention to the occupancy of the Y1 site that plays a key role in theaccording diversity to of[15]. the We vesu payvianite-group particular attention minerals. to the occupancy of the Y1 site that plays a key role in the diversity of the vesuvianite-group minerals. the diversityIn the crystal of the structure vesuvianite-group of vesuvianite, minerals. the five-coordinated Y1 site (Table 1 and Figure 6) can be occupiedIn theIn the crystalby crystal Fe3+ structure, Festructure2+, Mn of3+ ,of vesuvianite,Mn ve2+suvianite,, Cu, Mg the the and five-coordinatedfive-coordinated Al [16,23,30,51,52,55,65]. Y11 site site (Table (Table In most 1 andand cases, Figure Figure including 6)6) can can be be 3+ 3+ 2+2+ 3+3+ 2+2+ occupiedpreviouslyoccupied by by Festudied Fe, Fe, Fevesuvianite, Mn, Mn ,, MnMn from,, Cu, the Mg Somma-Vesuvius and Al Al [16,23,30,51,52,55,65]. [16,23 ,volcanic30,51,52 ,complex,55,65]. In In Femost most3+ dominates cases, cases, including including in the 3+3+ previouslyYpreviously1 site [48]. studied Instudied accordance vesuvianite vesuvianite with from the from 27 theAl the MAS Somma-Vesuvius Somma-Vesuvius NMR data, five-coordinated volcanic volcanic complex, complex, Al3+ Fe in Fe 27844dominates dominates is absent, in in and the the Y 1 27 27 3+3+ sitetheY [148 finalsite]. In [48]. occupancy accordance In accordance of withthe Ywith the1 site the Alis determined MASAl MAS NMR NMR by data, Fedata, and five-coordinated five-coordinated Mg. In the structures AlAl inin of 2784427844 27844 isis absent,and absent, 51062 and and, the final occupancy of the Y1 site is determined by Fe and Mg. In the structures of 27844 and 51062, theiron final is occupancypredominant of theat theY1 siteY1 site: is determined (Fe0.72Mg0.28 by)1.00Fe and and (Fe Mg.0.73Mg In0.27 the), structuresrespectively, of or,27844 takingand into51062 , iron is predominant at the Y1 site: (Fe0.72Mg0.28)1.00 and (Fe0.73Mg0.27), respectively, or, taking into ironaccount is predominant Mössbauer atdata, the FeY3+10.50 site:Mg0.28 (FeFe0.722+0.22Mg and0.28 Fe)1.003+0.47andMg0.27 (FeFe0.732+0.26Mg, respectively.0.27), respectively, Scattering or, takingof the Y into2 account Mössbauer data, Fe3+3+0.50Mg0.28Fe2+2+0.22 and Fe3+0.473+Mg0.27Fe2+0.26, respectively.2+ Scattering of the Y2 accountsite in Mössbauerboth samples data, is close Fe 0.50to 13Mg ē,0.28 andFe the0.22 averageand Fe 0.27bondFe lengths0.26, respectively. are equal to Scattering 1.896 and of the1.895siteY2 sitein Å both for in both 27844samples samples and is 51062 close is close, torespectively, 13 to ē 13, ande¯, and the which theaverage average agrees O> with bond bond the lengths site lengths occupancy are are equal equal by to toAl1.896 1.896atoms and and 1.895 Å for 27844 and 51062, respectively, which agrees well with the site occupancy by Al atoms 1.895only. Å forThe27844 octahedraland 51062 Y3 site, respectively, has minor admixtures which agrees of welliron with(mainly the siteFe3+), occupancy which is confirmed by Al atoms by only.a only. The octahedral Y3 site has minor admixtures of iron (mainly Fe3+), which is confirmed by a Theslightly octahedral higherY mean3 site length has minor of the admixtures bonds: of iron 1.966 (mainly and 1.961 Fe3+), Å, which respectively. is confirmed The total by refined a slightly slightly higher mean length of the bonds: 1.966 and 1.961 Å, respectively. The total refined higherY3 site mean occupancies length of are the (Al <0.90Y3–O>Fe0.10)1.00 bonds: and (Al 1.9660.91Fe and0.09), 1.961respectively. Å, respectively. Both the contents The total and refined oxidationY3 site Y3 site occupancies are (Al0.90Fe0.10)1.00 and (Al0.91Fe0.09), respectively. Both the contents and oxidation occupanciesstates of iron are have (Al beenFe confirmed) and by (Al the MössbauerFe ), respectively. data. Both the contents and oxidation states states of iron have0.90 been0.10 confirmed1.00 by the0.91 Mössbauer0.09 data. of iron have been confirmed by the Mössbauer data.

FigureFigure 6. 6.The The crystal crystal structure structure ofof 2784427844 projectedprojected along ( (AA)) cc axisaxis and and (B (B) )b baxis.axis. T-sitesT-sites have have low low Figure 6. The crystal structure of 27844 projected along (A) c axis and (B) b axis. T-sites have low occupation: T1 ≈ 10% and T2 ≈ 25%. The O11 site with mixed (O0.75F0.25) occupancy represent by blue occupation: T1 ≈ 10% and T2 ≈ 25%. The O11 site with mixed (O0.75F0.25) occupancy represent by circles,occupation: T1 atoms ≈ 10% are and red T2 circles. ≈ 25%. The O11 site with mixed (O0.75F0.25) occupancy represent by blue bluecircles, circles, oxygen oxygen atoms atoms areare red red circles. circles. The mean bond lengths and scattering factors of the tetrahedral Z1, Z2 and Z3 sites are in TheagreementThe mean mean O> full bond bond occupancieslengths lengths andand of these scatteringscattering sites byfactors Si atoms of of the the only. tetrahedral tetrahedral ZZ1,1, ZZ2 2and and Z3Z 3sites sites are are in agreementin agreement with with the the full full occupancies occupancies of of these these sites sites by by SiSi atomsatoms only.only.

Minerals 2017, 7, 248 10 of 16

Minerals 2017, 7, 248 10 of 16 In the refinement model, the 7- to 9-coordinated X1, X2, X3 and X4 sites are completely populated by Ca.In the The refinementX sites have model, full occupancies the 7- to 9-coordinated except for the X eightfold-coordinated1, X2, X3 and X4 sitesX4 positionare completely that is populatedhalf-populated by Ca. and The is situatedX sites have in the full structure occupancies channels. except According for the eightfold-coordinated to [66], REEs incorporate X4 position into the thatX3 site, is half-populated but their low contentand is situated (1.5% of in the theX 3structure site occupancy) channels. cannot According be precisely to [66], determined REEs incorporate by the intostructural the X3 analysis. site, but Similartheir low reasoning content (1.5% applies of alsothe X to3 thesite admixtureoccupancy) of cannot Na in thebe preciselyX4 site [22 determined]. by theIncorporation structural analysis. of small Similar amounts reasoning of boron applies into also the to vesuvianite the admixture structure of Na does in the not X4 leadsite [22]. to the splittingIncorporation of the O7 siteof small and the amounts appearance of boron of the into O12 th sitee vesuvianite as it does in structure wiluite [67 does]. The not low-occupied lead to the Tsplitting2 site (Figure of the 7O7) is site populated and the byappearance about 25% of boron,the O12 which site as prevented it does in wiluite us from [67]. the The determination low-occupied of Tits2 coordination.site (Figure 7) However,is populated the by presence about 25% of 3-coordinated boron, which boronprevented is confirmed us from the by FTIR.determination The observed of its − coordination.scattering factors However, of the Tthe1 sitepresence of 1.4 andof 3-coordi 1.0 e correspondnated boron to is theconfirmed refined occupanciesby FTIR. The Al observed0.106 and ‒ 0.106 0.072 Alscattering0.072, respectively. factors of the Due T1 to site the of low 1.4 occupanciesand 1.0 e correspond of the T1 to site, the the refined mean occupancies bond Al lengths, and Al 1.827, respectively.and 1.830 Å, areDue significantly to the low longeroccupancies than the of the theoretical T1 site, onethe formean tetrahedrally bond coordinated lengths, Al 1.827 with and full 1.830occupancy. Å, are At significantly the same time, longer the than presence the theoreti of smallcal amounts one for (~0.45 tetrahedrallyapfu) of tetrahedrallycoordinated Al coordinated with full occupancy.boron at the AtT1 the site same cannot time, be the excluded. presence of small amounts (~0.45 apfu) of tetrahedrally coordinated boronAll at oxygenthe T1 site sites cannot have beenbe excluded. refined as fully occupied by oxygen atoms, excluding the O7 site that possessesAll oxygen anomalously sites have high been displacement refined as fully parameters. occupied by In theoxygen final atoms, model, excluding the O7 site the wasO7 site refined that possessesas (O0.75F 0.25anomalously)1.00 with plausiblehigh displacement displacement parameters. parameters. In the Small final amountsmodel, the of O7 Cl determinedsite was refined by the as (Oelectron0.75F0.25) microprobe1.00 with plausible (0.05 anddisplacement 0.04 apfu for parameters.27844 and Small51062 amounts, respectively) of Cl aredetermined probably by situated the electron in the microprobeO10 site [19 (0.05]. and 0.04 apfu for 27844 and 51062, respectively) are probably situated in the O10 site [19]. Taking into into account account all all of of the the data data given given above, above, the final the finalcrystal-chemical crystal-chemical formulae formulae of the of27844 the X1 X1 X22.00X2 X8.003 X3 8.00XX44 0.97 0.03 1.00Y1 and27844 51062and 51062 samplessamples can can be be written,written, as as follows: follows:(Ca) 2.00(Ca) (Ca)(Ca)8.00 (Ca)(Ca)8.00 (Ca(Ca0.97Na0.03))1.00 Y1 3+0.503+ 0.28 2+0.22 2+Σ1.00Y2 3.85Y2 2+0.15 Σ4.00Y2+3 5.26 Y1.833 3+0.54 2+0.26 3+0.11 Σ8.002+Z1 2.00Z2 8.00Z3 Z1 8.00 (Fe(Fe Mg0.50MgFe0.28Fe) 0.22)(AlΣ1.00 Fe(Al3.85) Fe (Al0.15)Σ4.00Mg (AlFe 5.26MgFe 1.83MnFe 0.54) Fe (Si)0.26Mn0.11(Si))Σ8.00(Si)(Si)(O2.00 )Z68.002 T1+T2(AlZ3 0.44B0.25 4.31)Σ5.00TW1+(OHT2 5.65F2.00O1.30Cl0.05)Σ9.00W and X1(Ca)2.00X2(Ca)8.00X3(CaX1 7.90Ce0.06XNd2 0.04)8.00 (Si)8.00 (Si)8.00(O)68.00 (Al0.44B0.254.31)Σ5.00 (OH5.65F2.00O1.30Cl0.05)Σ9.00 and (Ca)2.00 (Ca)8.00 X43 1.00Y1 3+0.47 0.27 X2+40.26 Σ1.00YY2 1 Σ3+4.00Y3 5.01 2+2.10 3+0.53 Y22+0.13 0.12 Y3 0.11 Σ8.00Z1 2.003+Z2 8.002+Z3 (Ca)(Ca7.90Ce(Fe0.06NdMg0.04)8.00Fe (Ca)) 1.00 (Al)(Fe 0.47(AlMg0.27MgFe 0.26Fe)Σ1.00Fe (Al)TiΣ4.00Mn(Al5.01) Mg(Si)2.10Fe (Si)0.53Fe (Si0.13 )8.00(O)68.00T1+T2(AlZ0.291 B0.25 4.46Z2)Σ5.00W(OHZ3 5.88F2.00O1.08ClT0.041+T)Σ29.00. The Fe content inW each Y site was assigned Ti0.12Mn0.11)Σ8.00 (Si)2.00 (Si)8.00 (Si)8.00(O)68.00 (Al0.29B0.254.46)Σ5.00 (OH5.88F2.00O1.08Cl0.04)Σ9.00. accordingThe Fe content to the incrystal-structure each Y site was data, assigned whereas according the Fe2+ to/Fe the3+ ratio crystal-structure is given in accordance data, whereas with the MössbauerFe2+/Fe3+ ratio spectroscopy is given in data. accordance with the Mössbauer spectroscopy data.

Figure 7. Graphs of observed electron density along along the the cation cation rods rods (fourfold (fourfold axis axis is is centered centered at at 0.75; 0.75; 0.75; z) in the crystal structures of:of: (A) low-symmetry P4/nn vesuvianitevesuvianite from from Ak Akhmatovskayahmatovskaya Pit and −3 (B) 51062. Projections are ontoonto thethe (010)(010) plane,plane, contourcontour intervalsintervals areare 0.1 0.1e eÅÅ−3.. Displacement Displacement ellipsoids ellipsoids are drawn at the 50% probability level.

4.6. Infrared Spectroscopy The IR spectrum of 27844 (Figure 8) is typical for high-symmetry vesuvianite-group minerals. It contains strong bands of Si–O stretching in the Si2O7 (at 1018 and 976 cm−1) and SiO4 groups (at

Minerals 2017, 7, 248 11 of 16

4.6. Infrared Spectroscopy

MineralsThe 2017 IR, 7, spectrum 248 of 27844 (Figure8) is typical for high-symmetry vesuvianite-group minerals.11 of 16 −1 It contains strong bands of Si–O stretching in the Si2O7 (at 1018 and 976 cm ) and SiO4 groups 914(at 914cm−1 cm with−1 shoulderswith shoulders at 900 atcm 900−1 and cm 870−1 andcm−1 870), as cmwell−1 as), asSi–O–Si well asand Si–O–Si O–Si–O and bending O–Si–O vibrations, bending partlyvibrations, mixed partly with the mixed M–O with stretching the M–O vibrations, stretching where vibrations, M = Al, where Mg, Fe,M Ti,= Al,Mn Mg, (below Fe, 700 Ti,Mn cm− (below1). The band700 cm at− 7991). Thecm−1 band is assigned at 799 cmto bending−1 is assigned vibrations to bending of the M vibrations···O–H groups. of the M··· O–H groups.

FigureFigure 8. 8. InfraredInfrared spectrum spectrum of of 2784427844..

TheThe region region of of O–H O–H stretching stretching vibrations is dominated byby thethe peakpeak atat 35653565 cmcm−−11 (D(D band band according according toto [[68],68], whichwhich corresponds corresponds to theto the vibrations vibrations of the of O11–H1 the O11–H1 bond in bond the presence in the ofpresence F in the of neighboring F in the neighboringO11 site [69 ].O11 The site weak [69]. bands The atweak 3450 bands cm−1 atand 3450 3629 cm cm−1 and−1 are 3629 related cm−1 toare the related OH groups to the coordinatedOH groups coordinatedby Ti and to by the Ti and O11–H1 to the··· O11–H1·O7 fragment··O7 fragment with thewith angle the angle <120 <120°,◦, respectively respectively [69 [69].]. The The weakweak band atat 3210 3210 cm cm−1 −(J1 band)(J band) corresponds corresponds to the to OH the groups OH groups oriented oriented along the along c-axis the [68]c-axis forming [68] strong forming hydrogen strong bonds.hydrogen bonds. −−11 33−− TwoTwo weak bands at 1568 andand 14781478 cmcm areare assigned assigned to the BO 3 anionsanions (the (the TT2-centered2-centered triangle) triangle) withwith shortened (as compared toto wiluite)wiluite) B–OB–O bonds.bonds. TheThe shouldersshoulders atat 1170 1170 cm cm−−11 and 10751075 cmcm−−11 areare tentativelytentatively assigned assigned to to boron boron atoms atoms having having 4-fold 4-fold coordination. coordination.

5.5. Discussion Discussion GeneralizationGeneralization of of the the vesuvianite vesuvianite formula formula is is an an important important step step for for the the development development of of the the nomenclaturenomenclature forfor this groupthis of minerals.group Theof tentativeminerals. vesuvianite The formula,tentative (Ca,Na) vesuvianite19(Al,Mg,Fe) 13formula,(SiO4)10 (Ca,Na)(Si2O7)419(OH,F,O)(Al,Mg,Fe)1013,(SiO does4)10 not(Si2 specifyO7)4(OH,F,O) the occupancies10, does not specify of the Y the1(A,B) occupancies site [36], of which the Y are1(A,B) responsible site [36], whichfor the are observed responsible species for diversity the observed of the species vesuvianite-group diversity of the minerals vesuvianite-group [16,21,29]. According minerals to[16,21,29]. the IMA AccordingCNMNC guidelinesto the IMA for CNMNC compositional guidelines criteria for compositional [37], at least one criteria structural [37], at site least in one the structural potential newsite inmineral the potential should benew predominantly mineral should occupied be predominan by a differenttly occupied chemical by component a different than chemical that which component occurs thanin the that equivalent which occurs site in anin existingthe equivalent mineral site species in an [16 existing]. In this mineral context, species recently [16]. approved In this memberscontext, recentlyof the vesuvianite approved groupmembers with of Mn the3+ vesuvianite, Al3+, Mg2+ group, Cu2+ withas predominant Mn3+, Al3+, cationsMg2+, Cu at2+ the as Ypredominant1 site can be cationsconfused at the with Y1 vesuvianite site can be confused itself, according with vesuvianite to the present itself, formulaaccording of to this the mineral. present formula Therefore, of this it is mineral.crucially Therefore, important it to is determine crucially important a new vesuvianite to determine formula, a new which vesuvianite would formula, be in accordance which would with thebe inabove-mentioned accordance with IMA the above-mentioned rule. IMA rule. InIn most cases,cases, vesuvianitevesuvianite with with the theP 4/P4/nncnncsymmetry symmetry has has only only Fe andFe and Mg atMg the atY the1 site Y1 [55 site]. There[55]. Thereare only are three only publicationsthree publications available availa thatble report that report vesuvianite vesuvianite with with the Fe the2+ /FeFe2+3+/Feratio3+ ratio determined determined by bydirect direct chemical chemical or spectroscopic or spectroscopic methods methods (Figure9 )[(Figure21,29, 639)]. [21,29,63]. Our data demonstrate Our data thedemonstrate predominance the predominance of Fe3+ at the Y1 site for both vesuvianite samples from the type locality, which is consistent with previous conclusions [38,48]. According to our experimental results, the new formula of vesuvianite can be written as VII−IXX19VY1VIY12(ZO7)4(ZO4)10(W)10, where X are seven- to nine- coordinated sites of Ca with minor Na, K, Fe2+ and REE impurities; VY has square pyramidal coordination of dominant Fe3+ and subordinate Mg, Al, Fe2+ and Cu2+; VIY has octahedral coordination

Minerals 2017, 7, 248 12 of 16 of Fe3+ at the Y1 site for both vesuvianite samples from the type locality, which is consistent with previous conclusions [38,48]. According to our experimental results, the new formula of vesuvianite Minerals 2017, 7, 248 VII−IX V VI 12 of 16 can be written as X19 Y1 Y12(ZO7)4(ZO4)10(W)10, where X are seven- to nine-coordinated sites of Ca with minor Na, K, Fe2+ and REE impurities; VY has square pyramidal coordination of dominant and is dominantly occupied by Al with subordinate Mg, Fe2+, Fe3+, Mn2+, Mn3+, Ti, Cr and Zn; ZO4 = Fe3+ and subordinate Mg, Al, Fe2+ and Cu2+; VIY has octahedral coordination and is dominantly SiO4, sometimes with subordinate AlO4, (O4H4)4−, and W = OH, F, Cl and minor O [12,13,15,17,19,21‒ occupied by Al with subordinate Mg, Fe2+, Fe3+, Mn2+, Mn3+, Ti, Cr and Zn; ZO = SiO , sometimes 29]. Since the sum of occupancies of the T1 and T2 sites is below 1 apfu, their populations4 4 are not with subordinate AlO , (O H )4−, and W = OH, F, Cl and minor O [12,13,15,17,19,21–29]. Since the taken into account.4 4The4 simplified vesuvianite formula may be written as sum of occupancies of the T1 and T2 sites is below 1 apfu, their populations are not taken into account. 3+ Ca19Fe (Al,Mg,Fe)12(SiO4)10(Si2O7)4(OH,F,O)10. The idealized3+ charge-balanced formula of the The simplified vesuvianite formula may be written as Ca19Fe (Al,Mg,Fe)12(SiO4)10(Si2O7)4(OH,F,O)10. vesuvianite end-member without subordinate cations is Ca19Fe3+(Al10Me2+2)(Si2O7)4(SiO4)10O(OH)9, The idealized charge-balanced formula of the vesuvianite end-member without subordinate cations is 2+ 2+ 2+ where Ме3+ = Fe , Mg2+ , Mn . 2+ 2+ 2+ Ca19Fe (Al10Me 2)(Si2O7)4(SiO4)10O(OH)9, where Me = Fe , Mg , Mn .

Figure 9. Triangular compositional diagram of the Y1 site in the P4/nncnnc vesuvianitevesuvianite based based on on single- single- 2+ 3+ crystalcrystal X-ray diffraction withwith FeFe2+/Fe3+ ratioratio obtained obtained from from Mössbauer Mössbauer data data [63] [63] (black (black circle), circle), [21] [21] (red (red circle),circle), [29] [29] (orange circle) andand presentpresent studystudy (blue(blue circles)).circles)).

6. Conclusions Conclusions By revising vesuvianite-group of of minerals minerals [16,21,29] [16,21,29] two two vesuvianite samples samples 2784427844 and 51062 fromfrom type locality were investigated by electronelectron microprobe analysis, Mössbauer and infrared spectroscopy, TGA/DSC, MAS MAS NMR, NMR, single-c single-crystalrystal and and powder powder X-ray X-ray diffraction. diffraction. Present studies studies reveal reveal square-pyramidal square-pyramidal Y1 positionY1 position as a key as site a in key the site diversity in the of diversityvesuvianite- of groupvesuvianite-group minerals. On minerals. the basis On theof basisthis ofappr thisoach approach we wedetermine determine vesuvianite vesuvianite formula formula as 3+ 2+ Ca19FeFe3+(Al(Al1010MeMe2+2)(Si2)(Si2O27)O4(SiO7)4(SiO4)10O(OH)4)10O(OH)9, which9, which would would not not contradict contradict present present CNMNC CNMNC guidelines guidelines forfor compositional compositional criteria. The vesuvianite from from the Somma-Vesuvius comple complexx demonstrate disorder inside the structural channels and overall P4/nncnnc symmetry,symmetry, which which partially partially connected connected with with the the incorporation incorporation of of boron boron at at ◦ additional T1, T2 sites and agrees with the estimated temperature of crystallization at 800–900 °C.C.

Supplementary Materials: The following are available online at http://www.mdpi.com/2075-163X/7/12/248/s1. Supplementary Materials: The following are available online at http://www.mdpi.com/2075-163X/7/12/248/s1. Acknowledgments: This work was supported by the President of Russian Federation grant for leading scientific Acknowledgments:schools (to S.V.K.). Authors This work are was grateful supported to Sergey by the Aksenov, President Piotr of DzierRussianzanowski˙ Federation and grant Evgeny for Galuskin leading scientific for their schoolscontribution (to S.V.K.). to the experimentalAuthors are grateful studies andto Sergey useful Aksenov, discussions. Piotr Anonymous Dzierżanowski referees, and handling Evgeny Galuskin editor are for thanked their contributionfor critical reviews to the and experimental helpful comments. studies and useful discussions. Anonymous referees, handling editor are thanked for critical reviews and helpful comments.

Author Contributions: Taras L. Panikorovskii, Nikita V. Chukanov, Gregory Yu. Ivanyuk, Giuseppina Balassone and Sergey V. Krivovichev wrote the paper. Taras L. Panikorovskii conceived and designed the single- crystal, powder diffraction, TGA/DSC experiments. Nikita V. Chukanov performed and analyzed infrared data. Vyacheslav S. Rusakov performed and analyzed Mössbauer data. Vladimir V. Shilovskikh performed microprobe measurements. Anton S. Mazur performed MAS NMR measurements. Gregory Yu. Ivanyuk prepared the photo of vesuvianite species.

Minerals 2017, 7, 248 13 of 16

Author Contributions: Taras L. Panikorovskii, Nikita V. Chukanov, Gregory Yu. Ivanyuk, Giuseppina Balassone and Sergey V. Krivovichev wrote the paper. Taras L. Panikorovskii conceived and designed the single-crystal, powder diffraction, TGA/DSC experiments. Nikita V. Chukanov performed and analyzed infrared data. Vyacheslav S. Rusakov performed and analyzed Mössbauer data. Vladimir V. Shilovskikh performed microprobe measurements. Anton S. Mazur performed MAS NMR measurements. Gregory Yu. Ivanyuk prepared the photo of vesuvianite species. Conflicts of Interest: The authors declare no conflict of interest.

Appendix A

Table A1. Powder X-ray diffraction data for 27844 1.

dmeas Å Imeas dcalc Å Icalc hkl dmeas Å Imeas dcalc Å Icalc hkl 10.83 4 10.84 2 110 2.081 5 2.076 8 623 5.85 4 5.85 2 002 2.060 4 2.060 3 543 4.673 4 4.667 6 202 2.038 9 2.047 13 730 4.000 6 4.027 4 222 2.022 4 2.023 6 642 3.857 2 3.845 2 400 2.008 2 2.001 1 731 3.463 15 3.460 15 322 1.9941 8 1.9957 12 633 3.234 4 3.229 6 402 1.9602 6 1.6712 6 651 3.062 5 3.065 5 313 1.9281 2 1.9265 5 116 3.038 8 3.031 6 510 1.9171 2 1.9098 2 713 2.995 9 2.995 14 501 1.9029 4 1.9030 2 206 2.939 18 2.933 45 511 1.8823 7 1.8898 11 216 2.897 5 2.895 6 323 1.8647 2 1.8635 2 515 2.744 100 2.748 100 440 1.7960 1 1.7950 2 831 2.659 3 2.654 2 530 1.7614 9 1.7587 18 714 2.586 58 2.587 60 531 1.7158 2 1.7149 3 910 2.522 6 2.523 5 314 1.6768 17 1.6756 11 734 2.492 1 2.483 1 611 1.6603 22 1.6631 26 436 2.454 45 2.451 48 620 1.6215 30 1.6226 31 526 2.433 7 2.432 2 324 1.5853 2 1.5854 1 824 2.367 9 2.368 7 541 1.5712 5 1.5712 4 327 2.345 5 2.345 3 404 1.5580 7 1.5579 6 616 2.320 8 2.319 7 105 1.5397 2 1.5362 2 735 2.293 5 2.296 5 334 1.5318 1 1.5312 2 626 2.269 2 2.268 2 631 1.5235 4 1.5237 5 942 2.193 5 2.194 9 710 1.5089 2 1.5084 1 1021 2.177 5 2.169 5 701 1.5028 1 1.5029 1 1002 2.158 3 2.164 3 711 1.4979 8 1.4979 5 636 2.141 2 2.143 2 305 1.4753 1 1.4752 2 1022 2.121 11 2.122 14 315 1 The eight strongest lines are highlighted bold.

References

1. Deer, W.A.; Howie, R.A.; Zussman, J. Rock-Forming Minerals: Orthosilicates, Vol 1A, 1st ed.; Geological Society: London, UK, 1982; pp. 701–719. 2. Gnos, E.; Armbruster, T. Relationship among metamorphic grade, vesuvianite “rod polytypism”, and vesuvianite composition. Am. Mineral. 2006, 91, 862–870. [CrossRef] 3. Olesch, M. Natürliche und synthetische Fe-haltige Vesuviane. Fortschr. Mineral. 1979, 57, 114–115. (In German) 4. Kappeler, M.A. Prodromus Crystallographiae de Crystallis Improprie sic Dictis Commentarium; Heinrich, R.W., Ed.; HR Wyssing: Lucerna, Switzerland, 1723; pp. 1–43. 5. Russo, M.; Punzo, I. I Minerali del Somma-Vesuvio, 1st ed.; AMI: Cremona, Italy, 2004; pp. 1–317. 6. Werner, A.G. Klaproth’s Beiträge: Über Vesuvian; Klaproth, M.H., Ed.; Becker und Compagnie: Berlin, Germany, 1795; pp. 1–34. 7. Clarke, F. The constitution of the silicates. Bull. U. S. Geol. Surv. 1895, 125, 109–110. Minerals 2017, 7, 248 14 of 16

8. Sjögren, H. Analyser pa tvenne vesuvian-varieter och vesuvianens kemiska constitution I Allmänhet. Geol. I Stockh. Forh. 1895, 17, 267–271. [CrossRef] 9. Weibull, M. Studien über vesuvian. Z. Krystallogr. 1896, 35, 26–35. [CrossRef] 10. Jannasch, P.; Weingarten, P. Uber die chemische Zusammensetzung und Konstitution des Vesuvians und des Wiluits. Z. Anorg. Allg. Chem. 1896, 11, 40–48. [CrossRef]

11. Warren, B.E.; Modell, D.I. The structure of vesuvianite Ca10Al4(Mg,Fe)2Si9O34(OH)4. Z. Kristallogr. 1931, 78, 422–432. [CrossRef] 12. Machatschki, F. Zur Formel des Vesuvian. Z. Kristallogr. 1932, 81, 148–152. [CrossRef] 13. Coda, A.; Giusta, D.A.; Isetti, G.; Mazzi, F. On the structure of vesuvianite. Atti Accad. Sci. Torino 1970, 105, 1–22. 14. Groat, L.A.; Hawthorne, F.C.; Ercit, T.S.; Grice, J.D. Wiluite, Ca19(Al,Mg,Fe,Ti)13(B,Al,)5Si18O68(O,OH)10, a new mineral species isostructural with vesuvianite, from the Sakha Republic, Russian Federation. Can. Mineral. 1998, 36, 1301–1304. [CrossRef] 15. Groat, L.A.; Hawthorne, F.C.; Ercit, T.S. The chemistry of vesuvianite. Can. Mineral. 1992, 33, 19–48. 16. Armbruster, T.; Gnos, E.; Dixon, R.; Gutzmer, J.; Hejny, C.; Döbelin, N.; Medenbach, O. Manganvesuvianite and tweddillite, two new Mn3+-silicate minerals from the Kalahari manganese fields, South Africa. Mineral. Mag. 2002, 66, 137–150. [CrossRef] 17. Galuskin, E.V.; Galuskina, I.O.; Sitarz, M.; Stadnicka, K. Si-deficient, OH-substituted, boron-bearing vesuvianite from the Wiluy River, Yakutia, Russia. Can. Mineral. 2003, 41, 833–842. [CrossRef] 18. Britvin, S.N.; Antonov, A.A; Krivovichev, S.V.; Armbruster, T.; Burns, P.C.; Chukanov, N.V. Fluorvesuvianite, 2+ Ca19(Al,Mg,Fe )13[SiO4]10[Si2O7]4O(F,OH)9, a new mineral species from Pitkäranta, Karelia, Russia: Description and crystal structure. Can. Mineral. 2003, 41, 1371–1380. [CrossRef] 19. Galuskin, E.V.; Galuskina, I.O.; Dzierzanowski,˙ P. Chlorine in vesuvianites. Miner. Pol. 2005, 36, 51–61. 20. Hålenius, U.; Bosi, F.; Gatedal, K. Crystal structure and chemistry of skarn-associated bismuthian vesuvianite. Am. Mineral. 2013, 98, 566–573. [CrossRef] 21. Aksenov, S.M.; Chukanov, N.V.; Rusakov, V.S.; Panikorovskii, T.L.; Gainov, R.R.; Vagizov, F.G.; Rastsvetaeva, R.K.; Lyssenko, K.A.; Belakovskiy, D.I. Towards a revisitation of vesuvianite-group nomenclature: The crystal structure of Ti-rich vesuvianite from Alchuri, Shigar valley, Pakistan. Acta Crystallogr. B 2016, 72, 744–752. [CrossRef] [PubMed] 22. Panikorovskii, T.L.; Krivovichev, S.V.; Yakovenchuk, V.N.; Shilovkikh, V.V.; Mazur, A.S. Crystal chemistry of Na-bearing vesuvianite from fenitized gabbroid of the Western Keivy (Kola peninsula, Russia). Zap. Ross. Mineral. Obsh. 2016, 145, 83–95. (In Russian) 23. Panikorovskii, T.L.; Shilovskikh, V.V.; Avdontseva, E.Y.; Zolotarev, A.A.; Pekov, I.V.; Britvin, S.N.; 2+ Hålenius, U.; Krivovichev, S.V. Cyprine, Ca19Cu (Al,Mg)12Si18O69(OH)9, a new vesuvianite-group mineral from the Wessels mine, South Africa. Eur. J. Mineral. 2017, 29, 295–307. [CrossRef] 24. Panikorovskii, T.L.; Shilovskikh, V.V.; Avdontseva, E.Y.; Zolotarev, A.A.; Karpenko, V.Y.; Mazur, A.S.;

Yakovenchuk, V.N.; Krivovichev, S.V.; Pekov, I.V. Magnesiovesuvianite, Ca19Mg(Al,Mg)12Si18O69(OH)9, a new vesuvianite-group mineral. J. Geosci. 2017, 1, 25–36. [CrossRef] 25. Panikorovskii, T.L.; Chukanov, N.V.; Aksenov, S.M.; Mazur, A.S.; Avdontseva, E.Y.; Shilovskikh, V.V.;

Krivovichev, S.V. Alumovesuvianite, Ca19Al(Al,Mg)12Si18O69(OH)9, a new vesuvianite-group member from the Jeffrey mine, Asbestos, Estrie Region, Québec, Canada. Mineral. Petrol. 2017, 111, 833–842. [CrossRef] 26. Dyrek, K.; Platonov, A.N.; Sojka, Z.; Zabinski,˙ W. Optical absorption and EPR study of Cu2+ ions in vesuvianite (“cyprine”) from Sauland, Telemark, Norway. Eur. J. Mineral. 1992, 4, 1285–1289. [CrossRef] 27. Panikorovskii, T.L.; Krivovichev, S.V.; Zolotarev, A.A.; Antonov, A.A. Crystal chemistry of low-symmetry (P4nc) vesuvianite from the Kharmankul’ Cordon (South Urals, Russia). Zap. Ross. Mineral. Obsh. 2016, 145, 94–104. (In Russian) 28. Panikorovskii, T.L.; Krivovichev, S.V; Galuskin, E.V.; Shilovskikh, V.V.; Mazur, A.S.; Bazai, A.V. Si-deficient, OH-substituted, boron-bearing vesuvianite from Sakha-Yakutia, Russia: A combined single-crystal, 1H MAS-NMR and IR spectroscopic study. Eur. J. Mineral. 2016, 28, 931–941. [CrossRef] 29. Panikorovskii, T.L.; Mazur, A.S.; Bazai, A.V.; Shilovskikh, V.V.; Galuskin, E.V.; Chukanov, N.V.; Rusakov, V.S.; Zhukov, Y.M.; Avdontseva, E.Y.; Aksenov, S.M.; et al. X-ray diffraction and spectroscopic study of wiluite: Implications for the vesuvianite-group nomenclature. Phys. Chem. Mineral. 2017, 44, 577–593. [CrossRef] Minerals 2017, 7, 248 15 of 16

30. Giuseppetti, G.; Mazzi, F. The crystal structure of a vesuvianite with P4/n symmetry. Tsch. Mineral. Petrol. 1983, 31, 277–288. [CrossRef] 31. Xu, J.; Li, G.; Fan, G.; Ge, X.; Zhu, X.; Shen, G. Hongheite, IMA 2017-027. CNMNC Newsletter No. 39, October 2017, page 1283. Mineral. Mag. 2017, 81, 1279–1286. 32. Veblen, D.R.; Wiechmann, M.J. Domain structure of low-symmetry vesuvianite from Crestmore, California. Am. Mineral. 1991, 76, 397–404. 33. Allen, F.M.; Burnham, C.W. A comprehensive structure-model for vesuvianite: Symmetry variation and crystal growth. Can. Mineral. 1992, 30, 1–18. 34. Groat, L.A.; Hawthorne, F.C.; Ercit, T.S.; Putnis, A. The symmetry of vesuvianite. Can. Mineral. 1993, 31, 617–635. 35. Tanaka, T.; Akizuki, M.; Hudoh, Y. Optical properties and crystal structure of triclinic growth sectors in vesuvianite. Mineral. Mag. 2002, 66, 261–274. [CrossRef] 36. Fitzgerald, S.; Leavens, P.B.; Rossman, G.R.; Yap, G.P.A.; Rose, T. Vesuvianite from Pajsberg, Sweden, and the role of Be in the vesuvianite structure. Can. Mineral. 2016, 54, 1525–1537. [CrossRef] 37. Nickel, E.H.; Grice, J.D. The IMA Commission on New Minerals and Mineral Names: Procedures and guidelines on mineral nomenclature, 1998. Can. Mineral. 1998, 36, 913–926. [CrossRef] 38. Balassone, G.; Talla, D.; Beran, A.; Mormone, A.; Altomare, A.; Moliterni, A.; Mondillo, N.; Saviano, M.; Petti, C. Vesuvianite from Somma-Vesuvius volcano (Southern Italy): Chemical, X-ray diffraction and single-crystal polarized FTIR investigations. Period. Mineral. 2011, 80, 369–384. [CrossRef] 39. Balassone, G.; Bellatreccia, F.; Mormone, A.; Biagioni, C.; Pasero, M.; Petti, C.; Mondillo, N.; Fameli, G. Sodalite-group minerals from the Somma-Vesuvius volcanic complex, Italy: A case study of K-feldspar-rich xenoliths. Mineral. Mag. 2012, 76, 191–212. [CrossRef] 40. Balassone, G.; Kahlenberg, V.; Altomare, A.; Mormone, A.; Rizzi, R.; Saviano, M.; Mondillo, M. Nephelines from the Somma-Vesuvius volcanic complex (Southern Italy): Crystal chemical, structural and genetic investigations. Mineral. Petrol. 2014, 108, 71–90. [CrossRef] 41. Balassone, G.; Bellatreccia, F.; Ottolini, L.; Mormone, A.; Petti, C.; Ghiara, M.R.; Altomare, A.; Saviano, M.; Rizzim, R.; D’oraziom, L. Sodalite-group minerals from Somma-Vesuvius volcano (Naples, Italy): A combined EPMA, SIMS and FTIR crystal chemical study. Can. Mineral. 2016, 54, 583–604. [CrossRef] 42. Gilg, A.H.; Lima, A.; Somma, R.; Belkin, H.E.; De Vivo, B.; Ayuso, R.A. Isotope geochemistry and fluid inclusion study of skarns from Vesuvius. Mineral. Petrol. 2001, 73, 145–176. [CrossRef] 43. Hermes, O.D.; Cornell, W.C. Petrochemical significance of xenolithic nodules associated with potash-rich lavas of Somma-Vesuvius volcano. In NSF Final Technical Report; University Rhode Island: Kingston, RI, USA, 1978; p. 58. 44. Joron, J.L.; Metrich, N.; Rosi, M.; Santacroce, R.; Sbrana, A. Chemistry and petrography. CNR Quad. Ric. Sci. 1987, 114, 105–174. 45. Gruppo Mineralogico Geologico Napoletano. Available online: http://www.gmgn.it/index.html (accessed on 1 November 2017). 46. Sheldrick, G.M. A short history of SHELX. Acta Crystallogr. A 2008, 64, 112–116. [CrossRef][PubMed] 47. Agilent Technologies. CrysAlis CCD and CrysAlis RED; Oxford Diffraction Ltd.: Oxfordshire, UK, 2014. 48. Elmi, C.; Brigattim, M.F.; Pasqualim, L.; Montecchi, M.; Laurora, A.; Malferrari, D.; Nannarone, S. High-temperature vesuvianite: Crystal chemistry and surface considerations. Phys. Chem. Miner. 2011, 38, 459–468. [CrossRef] 49. Zabi´nski,˙ W.; Wactawska, Z.; Paluszkiewicz, C. Thermal decomposition of vesuvianite. J. Therm. Anal. 1996, 46, 1437–1447. [CrossRef] 50. Foldvari, M. Handbook of Thermogravimetric System of Minerals and Its Use in Geological Practice; Geological Institute of Hungary: Budapest, Hungary, 2011; pp. 1–180. 51. Phillips, B.L.; Allen, F.M.; Kirkpatrick, R.J. High-resolution solid-state 27Al NMR spectroscopy of Mg-rich vesuvianite. Am. Mineral. 1987, 72, 1190–1194. 52. Olejniczak, Z.; Zabi´nski,W.˙ 27Al NMR study of white vesuvianite from Piz Lunghin, Switzerland. Miner. Pol. 1996, 27, 41–45. 53. Yesinowski, J.P.; Eckert, H.; Rossman, G.R. Characterization of hydrous species in minerals by high-speed 1H MAS-NMR. J. Am. Chem. Soc. 1988, 110, 1367–1375. [CrossRef] Minerals 2017, 7, 248 16 of 16

54. Lager, G.A.; Xie, Q.; Ross, F.K.; Rossman, G.R.; Armbruster, T.; Rotella, F.J.; Schultz, A.J. Hydrogen-atom position in P4/nnc vesuvianite. Can. Mineral. 1999, 37, 763–768. 55. Ohkawa, M.; Yoshiasa, A.; Takeno, S. Crystal chemistry of vesuvianite: Site preferences of square-pyramidal coordinated sites. Am. Mineral. 1992, 77, 945–953. 56. Ohkawa, M.; Armbruster, T.; Galuskin, E. Structural investigation of low symmetry vesuvianite collected from Tojyo, Hiroshima, Japan: Implications for hydrogarnet-like substitution. J. Mineral. Petrol. Sci. 2009, 104, 69–76. [CrossRef] 57. Pavese, A.; Prencipe, M.; Tribaudino, M.; Aagaard, S.S. X-ray and neutron single-crystal study of P4/n vesuvianite. Can. Mineral. 1998, 36, 1029–1037. 58. Armbruster, T.; Gnos, E. “Rod” polytypism in vesuvianite: Crystal structure of a low-temperature P4nc with pronounced octahedral cation ordering. Schweiz. Mineral. Petrogr. 2000, 80, 109–116. 59. Galuskin, E.V. Vesuvianite-Group Minerals from Achtarandite Rocks (Wiluy River, Yakutia), 1st ed.; University of Silesia: Katowice, Poland, 2005; pp. 123–137. (In Polish) 60. Kraczka, J.; Zabi´nski,W.˙ Mössbauer study of iron in some vesuvianites. Mineral. Pol. 2003, 34, 37–44. 61. Rusakov, V.S.; Kovalchuk, R.V.; Borovikova, E.Y.; Kurazhkovskaya, V.S. State of iron atoms in high vesuvianites according to Mössbauer spectroscopy data. Zap. Ross. Mineral. Obsh. 2006, 135, 91–100. (In Russian) 62. Manning, P.G.; Tricker, M.J. Optical absorption and Mössbauer spectral studies of iron and site-populations in vesuvianites. Can. Mineral. 1975, 13, 259–265. 63. Groat, L.A.; Evans, R.J.; Cempírek, J.; McCammon, C.; Houzar, S. Fe-rich and As-bearing vesuvianite and wiluite from Kozlov, Czech Republic. Am. Mineral. 2013, 98, 1330–1337. [CrossRef] 64. Bruker AXS. Topas V4.2: General Profile and Structure Analysis Software for Powder Diffraction Data; Bruker AXS GmbH: Karlsruhe, Germany, 2009. 65. Groat, L.A.; Evans, R.J. Crystal chemistry of Bi- and Mn-bearing vesuvianite from Langban, Sweden. Am. Mineral. 2012, 97, 1627–1634. [CrossRef] 66. Fitzgerald, S.; Leavens, P.B.; Rheingold, A.L.; Nelen, J.A. Crystal structure of a REE-bearing vesuvianite from San Benito County, California. Am. Mineral. 1987, 72, 625–628. 67. Galuskin, E.V.; Galuskina, I.O.; Stadnicka, K.; Armbruster, T.; Kozanecki, M. The crystal structure of Si-deficient, OH-substituted, boron-bearing vesuvianite from the Wiluy River, Sakha-Yakutia, Russia. Can. Mineral. 2007, 45, 239–248. [CrossRef] 68. Groat, L.A.; Hawthorne, F.C.; Rossman, G.R.; Scott, T.E. The infrared spectroscopy of vesuvianite in the OH region. Can. Mineral. 1995, 33, 609–626. 69. Chukanov, N.V.; Panikorovskii, T.L.; Chervonnyi, A.D. On the relationships between crystal-chemical characteristics of vesuvianite-group minerals and their IR spectra. Zap. Ross. Mineral. Obsh. 2017, in print. (In Russian)

© 2017 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (http://creativecommons.org/licenses/by/4.0/).