Geochemistry and Genesis of Beryl Crystals in the LCT Pegmatite Type, Ebrahim-Attar Mountain, Western Iran

: Ebrahim-Attar granitic pegmatite, which is distributed in southwest Ghorveh, western Iran, is strongly peraluminous and contains minor beryl crystals. Pale-green to white beryl grains are crystallized in the rim and central parts of the granite body. The beryl grains are characterized by low contents of alkali oxides (Na 2 O = 0.24–0.41 wt.%, K 2 O = 0.05–0.17 wt.%, Li 2 O = 0.03–0.04 wt.%, and Cs 2 O = 0.01–0.03 wt.%) and high contents of Be 2 O oxide (10.0 to 11.9 wt.%). The low contents of alkali elements (oxides), low Na/Li (apfu) ratios (2.94 to 5.75), and variations in iron oxide (FeO= 0.28–1.18 wt.%) reveal a poorly evolved magmatic source of the beryl grains. Low abundances of rare earth elements ( Σ REE = 0.8–4.9 ppm) with high 87 Sr /86 Sr (i) ratios of 0.739 ± 0.036 for the beryl grains and 0.7081 for the host granites infer that the primary magma was directly produced by partial melting of the upper continental crust (UCC). The crystallization temperature of the Ebrahim-Attar granitic pegmatite changes from 586 to 755 ◦ C (average = 629 ◦ C), as calculated based on the zircon saturation index. Furthermore, the quartz geobarometer calculation shows that crystallization occurred at pressures of approximately 233–246 MPa. This pressure range is a promising condition for saturation of Be in magma. During granitic magma crystallization, the melt was gradually saturated with Be, and then beryl crystallized in the assemblage of the main minerals such as quartz and feldspar. Likewise, the host granite is characterized by high ratios of Nb/Ta (4.79–16.3) and Zr/Hf (12.2–19.1), and peraluminous signatures are compatible with Be-bearing LCT (Li-Ce and Ta) pegmatites.

The chemical composition of beryl could be useful to understand the evolution of granitic pegmatite [14].The concentrations of alkali elements (Na, Li, and Cs) and divalent elements (Fe, Mg) are important for understanding the geochemical signatures of the host granitic pegmatite [11,14,15].The temperature and composition of magma generally control the solubility of Be in the melt [16].Be has low solubility in magma with high alumina at low temperatures [16].Consequently, volatility has a critical role in increasing the solubility of Be, which has high solubility in melts that are rich in H 2 O, Li, P, F, and B [17].Therefore, a high concentration of Be shows magmatic [18,19] or hydrothermal sources [20,21].
As mentioned above, there are several geological settings, but herein, we focus on granitic pegmatite.Beryl is the main beryllium-bearing phase and is formed during the late stage of granite melts with high abundances of Li, Cs, and Ta, which is called LCT pegmatite [22,23].Most LCT pegmatites have peraluminous signatures that are determined by the melting of pelitic rocks related to the thickness of the continental crust during collision [24] to a postcollisional tectonic regime [25].Some beryl-bearing pegmatites have been reported in Iran, including Khajeh Morad in south Mashhad [26,27], Kamari-Zaman Abad-Dehno in south Hamedan [28,29], and Ebrahim-Attar in southwest Ghorveh [26,28,[30][31][32] (Figure 1).The Ebrahim-Attar granitic pegmatite of Ghorveh district in the Sanandaj-Sirjan zone (SaSZ) is one of the most important semigemstone sources in Iran (Figure 2a).The granitic pegmatite is widely exposed in two nearby district areas called Ebrahim-Attar 1 and Ebrahim-Attar 2. It is an active feldspar mine and, to a lesser extent, mica and beryl [30].
Herein, we focus on the Ebrahim-Attar beryl-bearing granitic pegmatite because there is a lack of information on the source and geochemistry of beryl.Overall, a few studies have been published on host granite without focusing on the beryl chemistry, genesis, and economic value of this semiprecious mineral.Therefore, major oxides, trace elements including rare earth elements, and Sr isotopes are measured in the separated mineral beryl to provide detailed information about the genesis and mineralization of the mineral and its host rock.Detailed mapping, field observation, petrography, and chemical analysis of beryl grains were performed to characterize the mineralization properties.The goal of this study is to illustrate the crystal chemistry of beryl grains.We investigate the main process (magmatic/hydrothermal) that leads to the formation of semigemstone beryl mineralization and the geochemical signature of source melt based on beryl chemistry.Trace element abundances of the beryl grains can track the mineralization pathway and provide valuable information about the geochemical features of the parental granite melt.In this research, we report the chemical signatures of beryl crystals by measuring the major elements by applying XRF, trace elements including rare earth elements (REEs), and Sr isotope ratios using ICP-MS and TIMS for the first time and provide more accurate and precise data to evaluate the geochemical features of host granitic pegmatite and beryl mineralization, respectively.[33] with positions of beryl-bearing pegmatites including (1) Khajeh Morad [26,27], (2) Kamari and Zaman Abad [28,29], and (3) Ebrahim-Attar granitic pegmatite (this study).

Geological Setting and Mineralogy
The Sanandaj-Sirjan zone (SaSZ) is a metamorphic-magmatic belt in western Iran and has a width of 50-150 km and a length of 1200 km [33][34][35].In the northern SaSZ, the Ghorveh district consists of a Jurassic metamorphic complex, a suite of greenschist, marble, amphibolite, and quartzite with several stages of deformation [30], and some igneous rock, aplite, and pegmatite with ages of 180-140 Ma intruded into the Jurassic metamorphic complex [36-39] (Figure 2a).The largest pegmatite outcrop is found in the Ebrahim-Attar Mountains in southwest Ghorveh, west of Iran (Figure 2b).
The Ebrahim-Attar Mountains contain a sequence of marble, schist, and metabasite, and granitic pegmatite occurs in the dolomitic layer.A small tungsten-bearing skarn zone occurs at their contact with a thickness of 1-2 m (Figures 1 and 2).This mountain is surrounded by the 157 Ma Moshir Abad granite [38] in the north and the 157.9 (±1.6)Ma Ghalaylan granite in the south [36] as determined by zircon U-Pb dating.Rb-Sr dating yielded an age of 102.5 ± 6.1 Ma for the Ebrahim-Attar pegmatite [30] (Figure 2b).The predominant signatures of this pluton are highly peraluminous leucogranite with megacrysts of K-feldspar.The Ebrahim-Attar granite has been well documented mineralogically, geochemically, and isotopically by Azizi et al. [30].[40][41][42].This shows that the Jurassic metamorphic complex is cut by various granites.(b) The geological map of the study area indicates that the Ebrahim-Attar granitic pegmatite cuts the Jurassic metamorphic rocks [41].The Ebrahim-Attar pegmatite has two parts-1 and 2.
The Ebrahim-Attar granite cuts the Jurassic marble in the study area (Figure 3a).Minor skarn mineralization was observed in the contact.The granitic intrusion from the rim to the center varies from fine-grained leucogranite to coarse-grained pegmatite (Figure 3b,c).In the center, the crystal size ranges from millimeters to tens of centimeters (Figure 3d-f).Although the Ebrahim-Attar rock is not a typically zoned pegmatite, based on our field observations, the mineral sequence and size of the crystals, it is divided into two subzones: (1) The thick border zone, which includes equigranular leucogranite.The mineralogy is mainly feldspar, quartz, and muscovite (Figure 3d).( 2) The inner zone is marked with mineral sequences of large feldspar+quartz+muscovite with less biotite and tourmaline (Figure 3e).A graphic texture is observed in the inner pegmatite, and intergrowth of quartz and feldspar is identified in a hand specimen (Figure 3f).Overall, the mineral assemblage of both subzones 1 and 2 consists of quartz, alkali feldspar (orthoclase and microcline), plagioclase (albite), muscovite, tourmaline, and beryl [30,32] and the beryl grains appear in pale green color (Figure 4).5).The chondrite-normalized REE pattern shows a flat LREE trend with slightly enriched HREEs (Figure 5) and pronounced negative Eu anomalies (Eu/Eu* = 0.01-0.14 and Figure 5).The Rb-Sr and Sm-Nd analyses for the Ebrahim-Attar pegmatites are also presented by Azizi et al. [30].The 87 Sr/ 86 Sr(i) ratios (=0.7081), 143 Nd/ 144 Nd (i) (=0.5122), and ENd(t) (=−5.8 to −1.6) were calculated based on 102 Ma.
The existence of beryl shows that the Ebrahim-Attar pegmatite was fertile [43].The chemical compositions of the host rocks and the color of beryl are proper indicators of the fractionation degree of magma [43].For example, less-evolved pegmatite has a low Ce concentration [31], and its beryl is green.Highly evolved pegmatite has a higher Cs content [32], and its beryl is white to pink.The Ebrahim-Attar pegmatite has low Cs content (Cs = 0.6-4.7 ppm), and its beryl is pale green, indicating that the beryl crystallized at a low rate of magma fractionation.
Large beryl crystals (10-20 cm) with euhedral and prismatic shapes and pale-green to white colors appear as mineral parageneses of quartz and feldspar in mega-crystals (Figure 4).

Analytical Method
In this research, ten homogeneous beryl grains from the pegmatite zone were separated, and seven fresh and unaltered crystals were selected for petrological and geochemical studies.The concentrations of major oxides and trace elements, including rare earth elements (REEs) and Sr isotope ratios, were measured in the selected grains in this research.Here, our analytical methods are briefly summarized below.
Approximately 10-20 g of each sample was powdered with an agate mill.First, 0.5 g of each sample powder was mixed with 5.0 g of lithium tetraborate, and a glass bead was prepared for XRF analysis.Major elements except Be were measured by WD-XRF (Rigaku ZSX Primus II) at Nagoya University (NU), Japan.For measurement of loss on ignition (LOI), 0.4 g of powder samples were taken in already-weighed quartz crucibles and heated at 1000 • C in a furnace for 6 h.Based on the mass difference before and after heating, the values in weight percent of LOI were calculated for each sample.
For the quantitative analysis of beryllium (Be) and trace elements, including REEs, approximately 30 mg of each sample was fused with 200 mg of potassium fluoride (KF) in a platinum crucible at 1000 • C for 10 min at NU.A blank sample of 200 mg KF was also fused and prepared.After cooling, it was dissolved in H 2 O, moved into an open PTFE beaker, and then dried.Next, several drops of concentrated HNO 3 (68%) were added to the dried samples and dried again.To exclude the remaining silicon (Si) completely from the dried samples, 1 mL of concentrated hydrofluoric acid (HF, 38%) and 0.5 mL of concentrated HNO 3 were added and dried.Then, 1 to 2 mL of 2% HNO 3 and 2 drops of concentrated HF were added to the dried sample until completely dissolved.These sample solutions were diluted with 2% HNO 3 for inductively coupled plasma mass spectrometry (ICP-MS) analysis.
For Sr isotope analysis with quantitative analysis of REEs, approximately 40 mg of each sample was digested by HF and perchloric acid (HClO 4 ) in an open PTFE beaker on a hotplate at a temperature of approximately 140 • C for 72 h at NU.After drying the digested solution, it was dissolved in 6 M HCl.White precipitates appeared in 6 M HCl, although the samples were once dissolved during HF digestion.The supernatants (dissolved parts) were separated and divided into two portions for 87 Sr/ 86 Sr analysis and quantitative analysis of REEs.Both portions were dried again.The isotope fraction was dried and then dissolved in 3 mL of 2.4 M HCl and loaded on a cation exchange column (BioRad AG50W-X8, 200-400 mesh) with HCl eluent to isolate Sr.The REE fraction was dissolved in 2% HNO 3 solution.The REE quantitative analysis data were compared with those obtained for the KF-fused samples for a cross-check.
Beryllium and trace elements were quantified by ICP-MS (Agilent 7700x) at NU. Since the KF-fused sample solution had high matrix element concentrations, isobaric interferences and blanks of alkaline earth elements were carefully checked for ICP-MS analysis.The internal standard of indium (In) was used to reduce the effect of ionization efficiencies for the high-matrix sample solution.The helium (He) collision mode and nogas mode were both used to check the isobaric interferences by molecular ions, including K + and F − .In addition, to check the accuracy of the trace element concentrations in the high-matrix sample solutions, two sample solutions with different dilution ratios were prepared and cross-checked.The blanks of some elements, such as Rb, were high (10 µg, Table 1) because of the impurities in the KF reagent.All data were corrected for blanks.Analytical errors for each element by ICP-MS were less than 10%.
Approximately 150 ng of the Sr fraction with 1 µL of 2 M H 3 PO 4 was loaded onto a Ta filament, and the isotope ratio of Sr was measured by VG Sector 54-30 thermal ionization mass spectrometry (TIMS) at NU. 87 Sr/ 86 Sr ratios were normalized to 86 Sr/ 88 Sr = 0.1194.NIST SRM 987 was used as the Sr isotope standard, and repeated analysis of the standard during the sample measurement yielded 0.710240 ± 0.000010 (1 SD, n = 17).
The REE abundances (Table 1) are low, and the total REE contents (ΣREEs) range from 1.21 to 5.41 ppm (average = 2.46 ppm).The total contents of light REEs (ΣLREEs) vary from 0.43 to 4.33 ppm (average = 1.72 ppm), and the HREE contents range from 0.131 to 0.244 ppm (average = 0.193 ppm).Yttrium contents range from 0.35 to 0.83 ppm (average = 0.56 ppm) (Table 1).The REE concentrations normalized by chondrite values [44] show flat LREE trends and strong depletions in HREEs, especially in Tb-Lu, and all of the patterns are similar to each other (Figure 5).All samples strongly show negative Eu anomalies (Eu/Eu* = 0.019-0.072,average= 0.033).The REE patterns of beryl grains are more fractionated than those of the host granitic pegmatite.The ΣREE values of beryl grains are lower than those of the host pegmatite (3.88-30.3ppm; [30]).In the host pegmatite, REE patterns are different from those of the beryl grains.The host pegmatites show slightly HREE-enriched patterns, although the beryl grains show HREE-depleted patterns.Meanwhile, both the beryl grains and their host rocks have negative Eu anomalies, which are clearly related to feldspar fractionation [45].

Crytallochemistry of Beryl
The general formula of beryl is well documented by Bragg and West (1926) [46].Beryl involves rings with 6 Si-O tetrahedra ( T1 ), which are associated with Be-O tetrahedra ( T2 ) at both vertical and lateral positions.Al-O is situated in an octahedral ( O ) site and is linked to six rings of Si-O tetrahedral position.The center of Si rings involves channels ( ch ) that prefer to place alkali elements such as Na + , Cs + , and K + and water [9,47].Al in the O site can be substituted by Fe 3+ , Cr, V, and Fe 2+ according to Equations ( 1) and ( 2): where O A 1+ is an alkali such as Na + , Cs + , and K + , O and C denote the site of substitution, ions of comparable size are shown as R 2+ and R 3+ , and charge balance is achieved by incorporation of alkalis in the channels.
Belov [48] suggested the substitution of Li + for Be 2+ .Aurisicchio et al. [47] re-evaluated the geochemistry of beryl based on a negative correlation between Be and Li and defined the mechanism of substitution as follows: where C A1 + is alkalis, and superscripts T2 and C denote the site of replacement.The beryl samples from the Ebrahim-Attar area show a negative correlation between Be and Li, supporting the substitution of Li at the tetrahedral position ( T2 ) (Figure 6a).Furthermore, a weak positive correlation between Fe+Mg and Al is evidence of minor substitution to the absence of substitution of Fe+Mg at the octahedral site (Figure 6b).According to the chemical composition of Ebrahim-Attar beryls, the apfu of the main component was calculated based on 18 oxygen atoms (Table 2).Thus, the chemical formula is calculated based on 18 oxygen atoms, and the average value of the main component of Ebrahim-Attar beryls is: ch (Na 0.054 ,K 0.014 ,Cs 0.001 ) T2 (Be The abovementioned replacement in beryl causes the colored crystal.Hence, the presence of Fe 2+ causes a green color in the Ebrahim-Attar beryls.

87 Sr/ 86 Sr Isotope Ratios in the Beryl Grains
The 87 Sr/ 86 Sr ratios of the seven beryl samples are listed in Table 3.The 87 Sr/ 86 Sr ratios show a wide range from 0.7381 to 0.7734 with an average of 0.7642 (n = 6), except for sample BRYL3, which has an extremely high value ( 87 Sr/ 86 Sr = 0.8060), and the separated beryl grains do not have a homogenous composition.Alkali elements such K, Rb, and Cs could occupy a large space in a silicon tetrahedral ring in the beryl structure, and charge balance generally occurs by substitution of Al 3+ in Si 4+ .More than 90% of various types of beryl have high concentrations of Rb (5-50 ppm) and low Sr (usually less than 1 ppm), resulting in a high 87 Rb/ 86 Sr ratio in beryl [49].Based on the concentrations of Rb and Sr in beryl, 87 Rb/ 86 Sr ratios were calculated for the seven samples (Table 3).The 87 Sr/ 86 Sr vs. 87 Rb/ 86 Sr diagram for the six beryl samples, except BRYL3, shows a weak, positive trend.The 87 Sr/ 86 Sr (i) ratio, 0.739 ± 0.036 (Figure 7), was calculated based on ISOPLOT software version 4.15.

Beryl and Its Parental Magma
The chemical composition of beryl grains is associated with the geochemical features of their host rocks [50,51].Barton 8b).The Ebrahim-Attar pegmatite has a high content of SiO 2 (71.4-81.0wt.%) and low content of Fe 2 O 3 (0.19-0.69 wt.%), MgO (0.04-0.12 wt.%), and CaO (0.33-1.73 wt.%), which are similar to peraluminous fertile granite [52].The Ebrahim-Attar pegmatites have low contents of ΣREEs (3.88-30.3ppm).Mostly, the ΣREE abundances in fertile granites vary from 1 to 20 times higher than those of chondrite [52] and are similar to those of the Ebrahim-Attar pegmatites (Figure 5).The Ebrahim-Attar granite contains green muscovite and tourmaline with coarse-grained graphic texture, which is a common feature in fertile granite [43,53].Fertile granite is defined as the parental granite of more evolved rare-metal pegmatite dikes [43].The granitic pegmatite is related to LCT pegmatite, which has high contents of Li, Cs, Ta, Nb, Rb, B, and Sn [16,43,55].Černý et al. [56] reported that in a fertilized system, lithophile REEs increase gradually during the progress of crystal fractionation from granite to pegmatite.Rare alkali elements such as Rb and Cs, high-field-strength elements (HFSE) such as Nb and Ta, and some ratios such as Rb/Sr, Nb/Ta, and Zr/Hf are key indicators of granite fractionation [57].
The Ebrahim-Attar granitic pegmatites have high ratios of Nb/Ta (4.79-16.3)and Zr/Hf (12.2-19.1)compared with the highly evolved Altay system [58] with Nb/Ta < 5 and Zr/Hf < 10 [58].Highly evolved granite is distributed with low Nb/Ta and Zr/Hf ratios [59].The K/Rb ratio in the more fractionated leucogranite is less than 50 [59] and varies from 69.6 to 194 in the Ebrahim-Attar granitic pegmatites.This is consistent with less-evolved pegmatite.The most evolved pegmatite is located in the distal position to the parental granite, and the less-evolved pegmatite is in the proximal position.Therefore, the distance of parental magma could control the degree of evolution [60].Considering these geochemical features, the Ebrahim-Attar granitic pegmatite can be classified as a less evolved, proximal pegmatite.
Herein, we compare the Ebrahim-Attar granitic pegmatite with Be-bearing LCT pegmatites, such as Zealand station dikes, pegmatitic-aplite beryl-rich dikes, and New Brunswick, Canada [57].The Ebrahim-Attar pegmatite has almost the same contents of Rb, Nb, Cs, and Ta and the same ratios of Nb/Ta, Zr/Hf, and Rb/Sr as the Zealand station dikes (Table 4).
Therefore, the crystallization pressure of Ebrahim-Attar granitic pegmatite ranges from 233-246 MPa at a depth near 8.0-9.5 km.
It is well documented that early alkali-poor beryl forms prismatic crystals.The prismatic habit appears in the late stage of pegmatite evolution.Additionally, alkali substitution in Be-tetrahedra causes lateral growth [69].Thus, the transition of a prismatic to a tabular pattern with alkali-rich beryl reflects the hydrothermal process in the late stage of pegmatite evolution [69,70].The accommodation of alkali elements such as Li, Na, K, Rb, and Cs and substitution of Be 2+ and Li + in the beryl structure along the c-axis cause the tabular beryl to form [71].The Ebrahim-Attar beryl has a prismatic euhedral shape and is interpreted as the magmatic source, consistent with the geochemistry of beryl (Figure 4).
With increasing alkali element contents in the Ebrahim-Attar beryls, the BeO value shows a weak decreasing trend (Figure 9), and the alkali element contents probably control the shape and Be concentration [72].This could explain the high BeO content (10.0 to 11.9 wt.%) in the alkali-poor beryl in the Ebrahim-Attar granite.
The contents of alkali elements (Na, Cs, and Li) and divalent metals such as Fe and Mg in beryl are also used to determine the geochemical signatures of host pegmatite [11,15,18].The Cs content of beryl increases during magma evolution (crystallization and differentiation) and is convenient to distinguish primitive (undifferentiated) magma from evolved magma [69,70].Beryl in evolved pegmatite becomes Cs-dominant [69].The Cs content is 106-260 ppm in the Ebrahim-Attar pegmatite beryl (Table 1).The evolution of pegmatitic melts can be inferred from the Na/Li vs. Cs relationship in beryl grains [9,73].Trueman and Cerny [73] show a progressive evolution trend in pegmatite based on the variation of Na/Li (apfu) and Cs content (wt.%) in beryl.Figure 10 shows that the Ebrahim-Attar beryl samples plot in field A of the magmatic trend from Trueman and Cerny [73], and it is concluded that it is a less-evolved pegmatite.The content of Na 2 O in the Ebrahim-Attar beryls varies from 0.24 to 0.41 (wt.%), and the FeO content ranges from 0.28 to 1.18 (wt.%); that is, the Ebrahim-Attar beryls plot in the magmatic field on the Na 2 O vs. FeO diagram [74] (Figure 11).One sample shows a high content of FeO (1.18 wt.%), which is probably because of fine biotite inclusions in the beryl grain.

Beryl Mineralization
Four main factors controlling the solubility and saturation of Be in granitic magma are (1) melt composition, (2) fractionation, (3) temperature [12,75], and (4) possible reaction of fluxing elements such as B, P, and F [75,76].The main source for LCT-type pegmatite is S-type granites [77], and the S-type granites have a high value of alumina, a high ratio of K/(K+Na), and high contents of Li and Cs [77], which are mainly produced from the anatexis of the metapelites rocks [77].
The Ebrahim-Attar granitic pegmatite samples show 87 Sr/ 86 Sr(i) = 0.7081 and 143 Nd/ 144 Nd(i) = 0.5122, calculated based on 102 Ma [30].The Ebrahim-Attar beryl samples show similar 87 Sr/ 86 Sr(i) = 0.739 ± 0.036.This feature clarifies a contribution of the older continental crust and/or metasediment [30] and sufficiently supports the source of Be-bearing LCT pegmatite.Clay minerals have a high amount of rare alkali elements [78].During metamorphism, white and black mica is produced due to the increase in temperature.These minerals provide a large amount of storage for alkali-earth elements, including Be [79].During partial melting, muscovite breaks down, and some trace elements, such as Be, enter the melt, leading to Be being sufficiently enriched and saturated in the pegmatiteforming magma to produce beryl [62].The Ebrahim-Attar granitic pegmatite is highly peraluminous with a high content of SiO 2 .The BASH (BeO-Al 2 O 3 -SiO 2 -H 2 O) mineral assemblage is stable when silica activity is high [4] and controls the diffusion of nonalkali elements [80].The degree of fractionation is an alternative process that controls the concentration of incompatible elements such as Be in the melt [60].Crystallization of quartz and feldspar increases the content of Be in the peraluminous granitic melt [74].
Selway et al. [43] suggested schematic complex zonation from parental peraluminous granite to rare-metal granite based on Černý [81] (Figure 12a).The contents of incompatible elements increase with the distance from the source of granitic magma and the degree of fractionation [43].The solubility of lithophile rare elements increases with increasing fluxing phases such as F, B, Li, and H 2 O and decreases with the crystallization temperature of pegmatite [82].These volatile components cause the melt to migrate far from the source [43], and successively, the incompatible elements can precipitate far from the original melt.Due to the decreasing temperature, the incompatible elemental solubility gradually decreases, and finally, some zonation occurs for these elements [43] (Figure 12).The elemental zonation from the source is as follows: Be pegmatite, Be-Nb-Ta pegmatite, Li-Be-Ta-Nb pegmatite, and Li-Be-Cs-Ta-Nb pegmatite.Higher contents of incompatible elements are observed in the most distal pegmatite [43] (Figure 12a).Schematic overview of regional zoning in a cogenetic parent + pegmatite group pegmatite increasing in the degree of evolution with increasing distance from the parent granite.Incompatible elements are more enriched in the most distal pegmatite [12,55].(a) By increasing the distance from peraluminous S-type magma, the fractionation degree increased, the volatile value could be enriched, the zoning became more complex, and rare metals precipitated in the distal part.(b) During the undercooling of the pegmatite dike, in the proximal part, barren pegmatite formed with the mineral sequence of plagioclase (An 10-15 ), muscovite, and K-feldspar, beryl was saturated at approximately 600 • C, and Be-pegmatite was produced.In the distal pegmatite, rare metal mineralization occurred.Modified from Selway et al. [43].
In general, the complexity of zonation and the degree of crystal fractionation in pegmatite increase with distance from the granite source [79].Unzonation pegmatites are the main sources for feldspar, mica, and quartz deposits [79].Beryl is the first rare-metal mineral crystallized in the proximal, muscovite-rich LCT pegmatite [18,73].Therefore, based on beryl dominance, Ebrahim-Attar granitic pegmatite could be classified as Bepegmatite, consistent with less-evolved granites.
Based on the model suggested by Selway et al. [43], Sn and W mineralization occurred at the contact of barren rocks with Be-pegmatite (Figure 12a).This is also consistent with tungsten mineralization as a skarn type in the contact of granitic hosts with carbonate rocks and near the Ebrahim-Attar Be-pegmatite [30].All lines of evidence show that the Ebrahim-Attar pegmatite was generated proximal to the parental granite.According to field observations, beryl occurs as megacrystals at the granite and inner pegmatite border, showing a higher degree of fractionation due to the host granite.Be-bearing pegmatite would be saturated at temperatures below 600 • C in peraluminous granite [52] (Figure 12b).Overall, the existence of Be over the saturation level and the lower temperature are two essential factors to initiate the growth of beryl [12].The crystallization temperature in the Ebrahim-Attar granitic pegmatite varies from 586 to 755 • C with an average of 629 • C and is lower than that in the inner part of the pegmatite, which is favorable for beryl crystallization.The temperature of the rare element pegmatite is lower (approximately 400-450 • C) because the margin of the pegmatite dyke is quickly cooled [77].London et al. [12] suggested that barren pegmatite occurs near the main silicate melt during the undercooling of silicate dike, with mineral sequences including plagioclase (An 10-15 ), muscovite, and K-feldspar (Figure 7b).When the temperature reaches approximately 600 • C, Be must be saturated to form Be-pegmatite (Figure 12b).The temperature for the Ebrahim-Attar pegmatite is consistent with Be-pegmatite on the models for LCT pegmatite of London et al. [12] and Selway et al. [43] (Figure 12).Be has low solubility in the melt because its solubility is related to Al 2 O 3 activity [16].The solubility of Be could be increased by the existence of a fluxing component [75].Fluxing components such as F, B, P, and H 2 O in metamorphosed sedimentary rocks lead to peraluminous feathers, which are essential for the genesis of LCT pegmatite [16].Furthermore, 87 Sr/ 86 Sr(i) (0.7081) and 143 Nd/ 144 Nd(i) (0.5122) provide evidence of continental crust origin for the Ebrahim-Attar granitic pegmatite [30].The peraluminous signature supports metasedimentary and Be sources in the Ebrahim-Attar granitic pegmatite.In the transition of peraluminous host rock to pegmatite, the Be content increases in the melt.In addition, the composition of plagioclase shifts from oligoclase in granite to albite in pegmatite, leading to some Be extraction.The Ebrahim-Attar pegmatite is not well zoned, with a low degree of evolution.These lines of evidence are consistent with proximal pegmatites with low rare metals.The proximal pegmatite contains low Be, indicating that the beryl could be saturated and precipitated in mineral sequences during cooling.

Conclusions
The highlighted conclusions on the geochemistry of the beryl and its host rock in the Ebrahim-Attar area are as follows: 1.
According to the mineral sequence, beryl is the only Be-bearing mineral, and its host rock is classified as Be-pegmatite in the Ebrahim-Attar area.

2.
The host rock has an LCT pegmatite and is a less-evolved pegmatite with high values of Nb/Ta (4.79-16.3)and Zr/Hf (12.2-19.1).

3.
Pale-green to white-colored beryl in the host granite appears as euhedral crystals and is alkali-poor beryl with low amounts of Na 2 O, K 2 O, Li 2 O, and Cs 2 O. 4.
Concentrations of Cs and Na/Li values in the beryl grains imply a magmatic signature and suggest a relation to less-evolved granite. 5.
High 87 Sr/ 86 Sr ratios of both the beryl and its host granite are consistent with continental crust contributions from their sources.

Figure 2 .
Figure 2. (a) Simplified geological map of the Ghorveh area based on 1:100,000 geological maps of Ghorveh, Songhor, and Sanandaj[40][41][42].This shows that the Jurassic metamorphic complex is cut by various granites.(b) The geological map of the study area indicates that the Ebrahim-Attar granitic pegmatite cuts the Jurassic metamorphic rocks[41].The Ebrahim-Attar pegmatite has two parts-1 and 2.

Figure 4 .
Figure 4. (a-d) Euhedral hexagonal prismatic crystal of beryl grains with pale-green to white colors in the Ebrahim-Attar granitic pegmatite.

Figure 6 .
Figure 6.(a) Li vs. Be diagram for the Ebrahim-Attar beryl samples, showing a negative correlation, which supports that Li is substituted with Be.(b) Fe + Mg vs. Al diagram, showing a weak positive correlation, which supports the minor substitution of Fe+Mg at the octahedral site.Li, Be, Fe, Mg, and Al (apfu) are calculated based on 18 oxygen atoms.

Figure 9 .
Figure 9.The plot of total alkali element content (Na 2 O + K 2 O + Li 2 O + Cs 2 O) vs. BeO of the Ebrahim-Attar beryl shows a weak negative correlation.

Figure 10 .
Figure 10.Diagram of Na/Li vs. Cs to compare the geochemical compositions of Ebrahim-Attar beryl with magmatic trend proposed by Trueman and Cerny [73] and divided pegmatite into four groups: (A) barren and geochemically primitive beryl type and less-evolved pegmatite; (B) geochemically evolved beryl-columbite and beryl-columbite-phosphate pegmatite; (C) albite-spodumene and complex pegmatite; (D) highly fractionated Li, Cs, and Ta-rich complex pegmatite.The beryl samples are plotted in less-evolved pegmatite.

Figure 12 .
Figure12.Schematic overview of regional zoning in a cogenetic parent + pegmatite group pegmatite increasing in the degree of evolution with increasing distance from the parent granite.Incompatible elements are more enriched in the most distal pegmatite[12,55].(a) By increasing the distance from peraluminous S-type magma, the fractionation degree increased, the volatile value could be enriched, the zoning became more complex, and rare metals precipitated in the distal part.(b) During the undercooling of the pegmatite dike, in the proximal part, barren pegmatite formed with the mineral sequence of plagioclase (An 10-15 ), muscovite, and K-feldspar, beryl was saturated at approximately 600 • C, and Be-pegmatite was produced.In the distal pegmatite, rare metal mineralization occurred.Modified from Selway et al.[43].

Table 3 .
Rb-Sr abundances and Sr isotope ratios of the beryl grains.