Micro-Raman—A Tool for the Heavy Mineral Analysis of Gold Placer-Type Deposits (Pianu Valley, Romania)

In the current study, different heavy minerals typical of gold placer deposits were identified by means of micro-Raman spectroscopy, and their chemical composition analyzed and discussed (garnet, kyanite, staurolite, zircon, allanite, monazite, xenotime, rutile, anatase, cassiterite, titanite, barite). Even complex solid solution series, such as those of garnets, can be deciphered with the aid of systematic trends observed in Raman line frequencies. The ν1 mode in garnets will shift from high to low frequencies as a function of the ionic radius of the X2+ cation, from Mg2+, to Fe2+ and Mn2+, while the presence of Ca2+ will make the band to be shifted strongly to even lower wavenumbers. This approach has successfully been taken to differentiate between polymorph triplets such as kyanite-sillimanite-andalusite and rutile-anatase-brookite. Minerals under consideration with high contents of REE, U and Th are affected by intensive metamictization, particularly zircon and titanite. Raman peak features, such as shape, symmetry and intensity, respond to this radiation damage of the lattice and enable fine-tuning of these heavy minerals, such as in the case of fluorite (fetid fluorite).


Introduction
Detrital heavy minerals provide insights into the mineral chemistry of the entire source area, and when cast in the role of provenance markers, environment indicators and ore guides, can be successfully used in applied and genetic geosciences [1][2][3][4][5].
Placers are clastic sedimentary deposits formed by mechanical and chemical processes operative in aeolian, fluvial and coastal-marine environments. In particular, weathering of metamorphic, igneous or reworked sedimentary source rocks, sorting on transport and in situ alteration at the trap site play paramount roles during formation of placer-type deposits [6][7][8]. When it comes to an economic assessment of placer deposits, it is especially gold, platinum group elements (PGE), platinum group minerals (PGM), Ti minerals, zircon and rare-earth elements (REE) minerals that attract the attention of economic geologists and mineralogists [9,10]. There is an ever-increasing demand for heavy minerals worldwide in various branches of the industry, mainly those dealing with new green technologies and The Pianu Valley is renowned in Romania for its important alluvial gold deposits. The area has a long-lasting history of gold exploration and exploitation reaching back to antiquity. Traces of panning gold by rudimentary and traditional methods can still be observed today in the environs of Pianu de Sus [31]. Pianu Valley takes an outstanding position as the only site where native platinum has been reported in Romania [32]. The following heavy minerals were identified as major components within the Pianu Valley: garnet, zircon, titanite, magnetite, tourmaline, kyanite, hornblende, Cr-spinel, monazite [32][33][34].
The Pianu Valley is located in the north-eastern part of the Sebes , Mountains and is made up of Upper Cretaceous and Quaternary sedimentary rocks that rest on the Getic Domain or Median Dacides, according to Săndulescu [35]. In terms of the lithofacies, it belongs to the Sebes , -Lotru variegated litho-group, which is characterized by gneisses, amphibolites, ultrabasic magmatic rocks, quartzites, migmatites, kyanite-bearing micaschists, meso-and epimetamorphic schists, granites, pegmatites and quartziferrous porphyries [31,36]. The Getic Realm is intersected by several shear zones [37,38]. These structures are considered as the primary source of the Pianu Valley gold and rare-elements accumulations [32,36,39].
The alluvial-fluvial drainage system spread along the Pianu Valley is Cretaceous in age, and its rocks are denominated as conglomerate, breccia, sandstone, shale, marl, and clay, while the Pleistocene to Holocene deposits include sand, gravel with quartzite boulders [40].

Materials and Methods
Heavy minerals were sampled in the Pianu Valley area at different locations at which gold occurrences had been identified. The samples were pre-concentrated in the field by manual panning with a 35 cm diameter plastic pan in order to remove as much as possible the light minerals (quartz, feldspars, carbonates). The samples were subsequently enriched in the laboratory using diiodomethane heavy liquid (CH 2 I 2 ), density of 3.325 g/mL. After the separation procedure with heavy liquid, the resulting fractions consisted of mineral grains with high specific gravity. These heavy mineral fractions were embedded in epoxy resin and prepared as polished sections.
All samples were examined by means of micro-Raman spectroscopy and Electron Probe Microanalysis (EMPA) at Slovak Academy of Sciences from Banská Bystrica, Slovakia. The microprobe measurements were performed on carbon-coated polished sections using the JEOL JXA-8530F. The Raman spectra were obtained on polished sections using the micro-Raman Spectrometer LabRAM HR800 (Horiba Jobin-Yvon, Kyoto, Japan) coupled with an Olympus BX51 optical petrographic microscope (100× /0.80 magnifying objective) and a CCD detector (1024 × 256 pixels). The samples were irradiated using the 632.8 nm excitation line of a He-Ne laser with a power of 17mW. Several density filters were used to prevent from heating, deterioration and fluorescence effect of the samples. The micro-Raman spectra were recorded in the range of 50-1400 cm −1 with an exposure time of 2-20 s per frame, 10-20 acquisitions, in order to improve the signal-noise ratio. Spectral accuracy was calibrated with Rayleigh line and the 520.7 cm −1 line of a Si standard. The system resolution was 2 cm −1 , and the wavenumber accuracy was ±1 cm −1 . The band fitting and was done by combined Gaussian/Lorentzian amplitude function. The Raman spectra were reduced in intensity for proper display in stacked graphs and better comparison.

Results and Discussion
Raman investigations of the heavy minerals from the Pianu Valley placer deposit revealed the presence of the following minerals: garnet, kyanite, staurolite, zircon, allanite, monazite, xenotime, rutile, anatase, cassiterite, titanite and barite. All minerals were confirmed by chemical analyses through electron microprobe measurements (EPMA). Their vibrational features are discussed in detail in the next sections.
The chemical analyses show that most of the studied garnets have an almandine-rich composition with variable contents of pyrope, spessartine and grossular (Table S1). The main cation in Y site is Al 3+ Representative Raman spectra are shown in Figure 1, and the Raman bands are listed in Table 1 together with their assignment. The spectra were selected from a total of 27 acquired on different garnet compositions.
A shifting behavior was observed for all Raman peaks based on the chemical composition. These variations are quite complex, since we are dealing with hexary systems such as garnet groups. Furthermore, we will discuss in detail the most intense peak dominating the garnet spectra, which was attributed to ν 1 vibrational modes of SiO 4 , as being mostly representative in behavior also of the other bands.
The For a better understanding of these peak variations, it has to be considered that line shifting in pyralspites is caused by atomic mass, atomic structure, ionic radius and polarization of X 2+ cations, which strongly affect the unit-cell size and bonds force constants in the crystal structure [43,46]. Out of all these factors, the ionic radius has the strongest influence on the position of the Raman peaks. In pyrope-almandine-spessartine-grossular series, the ionic radii increase from Mg 2+ (0.89) to Fe 2+ (0.92) , Mn 2+ (0.96) and Ca 2+ (1.12) Å. Therefore, the ν 1 mode will shift from high to low frequencies from Mg 2+ , to Fe 2+ and Mn 2+ , while the presence of Ca 2+ will cause the band to be shifted strongly to even lower wavenumbers because of the much larger ionic radius [17]. This decrease of the wavenumbers with increasing ionic radius is also supported by the correlations with the chemical results discussed above.
Since in the case of garnets solid solutions we are dealing with four different cations that can substitute each other and form mixed intermediate members, it is important to consider how much each cation contributes with its ionic radius to the X 2+ site occupancy. By calculating a weighted average of the ionic radii and molar ratios for each garnet composition, it can be observed that it strongly correlates with the measured Raman ν 1 mode wavenumbers (r = −0.86). This value can be calculated as (r Mg ·x Pyr + r Fe ·x Alm + r Mn ·x Spe + r Ca ·x Gro ), where r is the ionic radius and x represents the mol% of each end-member. Another similar approach is to use the ν 1 mode frequencies of pure end-members as standard values weighted by their molar ratio. The ν 1 band appears at 926 cm −1 in pyrope, 915 cm −1 in almandine, 907 cm −1 in spessartine, and 879 cm −1 for grossular [16]. Within this study, by using these values together with the end-members' mol%, a very good correlation with the real frequencies of the ν 1 mode was observed (r = 0.88). This confirms that by applying this procedure backwards and taking into account up to six peaks with unimodal behavior, one can roughly estimate the chemical concentration of any garnet intermediate member using only the Raman analyses, as previously discussed in more detail by Bersani et al. [16] and Smith [17].
As for the other vibrational modes that are present in the Raman spectra of the studied garnets, it was noticed that ν 3 (SiO 4 ) modes at 852-867, 1028-1047 cm −1 , ν 2 (SiO 4 ) modes at 466-482, 553-559 cm −1 , and T(X 2+ ) mode at 162-170 cm −1 show a similar trend as the already discussed ν 1 (SiO 4 ) vibrations. An inverse behavior was noticed for the peaks at 211-251, 343-355 and 369-379 cm −1 , which are assigned to R,T(SiO 4 ). The other peaks show no direct correlation with the chemical data, most likely having a multi-modal character.
The Raman spectra of kyanite ( Figure 2) show a large number of clear and sharp peaks, especially in the low wavenumber region. The bands from the 892-1001 cm −1 interval can be attributed to the ν 1 and ν 3 stretching modes of SiO 4 (symmetric and antisymmetric). Two extra peaks were observed in the studied samples at 1113 and 1185 cm −1 (Figure 2-kyanite 2). As for the ν 2 and ν 4 bending modes, these appear in the spectra between 607 and 669 cm −1 [19,47]. The strong peak at 485-487 cm −1 can be assigned to bending vibrations of Al-O-Si bonds, as suggested by Zhai et al. [19]. In the low region of the spectra, below 450 cm −1 there are several Raman bands that probably also overlap with the external modes. Zhai et al. [19] assigned these bands to Al-O bending modes and explained it as a cause of poor symmetry and degeneracy in the kyanite structure. The Raman spectra of kyanite ( Figure 2) show a large number of clear and sharp peaks, especially in the low wavenumber region. The bands from the 892-1001 cm −1 interval can be attributed to the ν1 and ν3 stretching modes of SiO4 (symmetric and antisymmetric). Two extra peaks were observed in the studied samples at 1113 and 1185 cm −1 (Figure 2-kyanite 2). As for the ν2 and ν4 bending modes, these appear in the spectra between 607 and 669 cm −1 [19,47]. The strong peak at 485-487 cm −1 can be assigned to bending vibrations of Al-O-Si bonds, as suggested by Zhai et al. [19]. In the low region of the spectra, below 450 cm −1 there are several Raman bands that probably also overlap with the external modes. Zhai et al. [19] assigned these bands to Al-O bending modes and explained it as a cause of poor symmetry and degeneracy in the kyanite structure.
Some differences were observed between kyanite Raman spectra from different analyzed samples. For instance, all peaks from the 200-564 cm −1 region seem to be shifted ~3 cm −1 (Figure 2kyanite 1 and 2). Instead, the stretching modes showed an opposite trend: from 892 to 900 cm −1 and 952 to 962 cm −1 . Although the last-mentioned peak might be in fact two overlapping bands with varying intensities (sample orientation), the clear shifting of the other lines might be explained by a slight change in chemical composition. As the EPMA results showed (Table S2), the kyanite 2 sample has a higher content of FeO (0.29 wt.%) and TiO2 (0.11 wt.%), which might cause structure deformations. Nevertheless, the Raman spectroscopy proved again to be very useful in identifying kyanite out of the other polymorphs (sillimanite and andalusite) based on structural features, since the chemical analyses were not enough to clearly differentiate between them.  Some differences were observed between kyanite Raman spectra from different analyzed samples. For instance, all peaks from the 200-564 cm −1 region seem to be shifted~3 cm −1 (Figure 2-kyanite  1 and 2). Instead, the stretching modes showed an opposite trend: from 892 to 900 cm −1 and 952 to 962 cm −1 . Although the last-mentioned peak might be in fact two overlapping bands with varying intensities (sample orientation), the clear shifting of the other lines might be explained by a slight change in chemical composition. As the EPMA results showed (Table S2), the kyanite 2 sample has a higher content of FeO (0.29 wt.%) and TiO 2 (0.11 wt.%), which might cause structure deformations. Nevertheless, the Raman spectroscopy proved again to be very useful in identifying kyanite out of the other polymorphs (sillimanite and andalusite) based on structural features, since the chemical analyses were not enough to clearly differentiate between them.
The Raman spectra of the staurolite samples show less peaks than kyanite and have a lower quality with rather broad bands and background noise (Figure 2). This is due to the overlapping modes and the low degree of crystallinity of the analyzed staurolite samples. A similar spectrum was reported previously by Makreski et al. [47], but with some differences in peak positions. According to the assignment of those authors, the peaks observed at 895, 939 and 1000 cm −1 can be attributed to the stretching modes of SiO 4 tetrahedra, while the bending vibrations occur at 639 and 796 cm −1 . The peaks from the low region of the spectra appear between 304 and 500 cm −1 and might be assigned to Al-O vibrational modes.
The analyzed staurolites suggest a shifting of the 939 cm −1 band to lower wavenumbers in samples with increasing TiO 2 concentrations correlated with a decreased ZnO, FeO and MgO content, while the other peaks maintain the same frequencies.

REE Minerals
Among the minerals enriched in rare earth elements (REE), zircon, allanite, monazite and xenotime were identified in the heavy mineral fractions of Valea Pianu samples.

Zircon
The Raman spectra acquired on the studied zircons show some differences, as indicated in Figure 3. The most intense peaks appear at 352-358, 437-440 and 1002-1009 cm −1 . The 1002-1009 cm −1 band is attributed to the ν 3 stretching mode of SiO 4 tetrahedra, while the ν 1 mode was observed at 970-974 cm −1 . Another internal mode was identified at 437-440 cm −1 that was assigned to the ν 2 bending vibration [18,48,49,53,54]. The 352-358 cm −1 line was previously assigned to the ν 4 bending mode [55], but more detailed studies based on group-theory analysis of the total number of Raman-active modes and structural investigations concluded that this line is due to external rotational modes of SiO 4 [18,48,49,53]. The bands from the low region of the spectra, 202, 214 and 224 cm −1 are also attributed to the external modes [49]. Additional bands were observed (Figure 3) which are not consistent with those reported in the literature and might be caused by sample impurities. less damaged. This might suggest that the line shifting may appear even in non-metamict zircons, but it's actually due to substitutions in cationic site which are causing structural deformations. This radiation damage of the crystal structure can also be observed in the Raman spectra of common ore mineral called fetid fluorite, a heavy mineral scattered in the clastic apron around U-bearing fluorite deposits [56]. The name has been derived from native fluorine in the lattice which spreads a garliclike smell.  There are some differences to be observed in the Raman spectra of the zircons under study that are considered a result of metamictization. Since natural zircon contains significant amounts of U and Th, the structure can strongly be damaged in the course the bombardment of alpha particles during the radioactive decay [18,48,49]. The metamict zircons have a high degree of structural disorder and consist of a mixture of defect-rich clusters with amorphous regions on a submicrometric scale [18]. These structural changes occur as point defects, SiO 4 tetrahedra polymerization, appearance of ZrO 2 , and even a complete transformation to an amorphous state. The metamict zircons were extensively studied by means of Raman spectroscopy [18,48,49], and it was observed that the degree of damage is reflected by a broadening of the Raman bands, lower peak intensities and lines shifting to smaller wavenumbers. Fresh zircons display sharp and well resolved peaks as those represented in Figure 3-zircon 1 and 2. Therefore, zircon 3 crystal (Figure 3) might be affected by a certain degree of metamictization due to its broader Raman bands observed at lower frequencies. Looking at its chemical composition (Table S3) the U and Th contents are relatively low (0.07 wt.% UO 2 and 0.01 wt.% ThO 2 ). Out of all analyzed zircons, the one with the highest U and Th values (0.16 wt.% UO 2 and 0.10 wt.% ThO 2 ) shows shifted peaks at lower frequencies as well (i.e., ν 3 mode-1005 cm −1 ). Judging from its very clear spectrum and sharp Raman lines, it can be considered that its structure is less damaged. This might suggest that the line shifting may appear even in non-metamict zircons, but it's actually due to substitutions in cationic site which are causing structural deformations. This radiation damage of the crystal structure can also be observed in the Raman spectra of common ore mineral called fetid fluorite, a heavy mineral scattered in the clastic apron around U-bearing fluorite deposits [56]. The name has been derived from native fluorine in the lattice which spreads a garlic-like smell.

Allanite
Allanite was identified in the present study as a single crystal with three different compositional zones. The mineral is of allanite-(Ce) type, with Ce being partially substituted by La and Nd. The literature data is very limited concerning the Raman studies of allanite. Lopez and Frost [51] reported the Raman spectrum of an "allanite" sample, but their results are considerably different from those obtained in the present study. Most likely, the authors misidentified their sample, since no chemical or X-ray analyses can confirm their assumption. The Raman spectra of allanite identified in Valea Pianu samples (Figure 3) better resemble a Cr-rich allanite-group mineral reported by Varlamov et al. [57], and the Raman frequencies fit perfectly with those reported by Chukanov and Vigasina [58] for allanite-(Ce). A tentative assignment of the Raman modes can be made by comparing it with epidote spectrum [59]. Therefore, the observed peaks in the high-frequency region (872, 961 and 1046 cm  (Figure 3). No peak shifting was observed in the spectra that might be due to different chemical composition.

Monazite
The Raman spectrum of monazite shows typical features of phosphate minerals (Figure 3). The spectra are dominated by the intense line at 972-975 cm −1 that is assigned to symmetric stretching ν 1 modes of PO 4 units while the antisymmetric ν 3 stretch is observed at 1055-1060 cm −1 . The bending vibrations are attributed to the Raman bands from 620 and 570 cm −1 (ν 4 modes) and from 464-472 cm −1 (ν 2 modes). The remaining lines that appear below 450 cm −1 are due to the external or lattice modes that involve M-O movements: 171, 222, 396 and 415 cm −1 [20,21,60,61].
Some differences were observed in the Raman lines position for the measured monazites, especially for the symmetric modes ν 1 (972-975 cm −1 ) and ν 2 (464-472 cm −1 ). The studied samples are monazite-(Ce) with additional La, Nd and Th. It seems that there is a systematic shifting caused by the dominant cation that occupies M site. The stretching bands tend to shift to higher wavenumbers with the decrease of ionic radius of the main cation or with its increasing atomic number in the case of REE [20]. For the monazite with the highest frequency of the ν 1 mode (975 cm −1 ), it was measured to have the highest content of Y (2.04 wt.% Y 2 O 3 ), which has a smaller ionic radius (Table S4).
The Raman spectrum of xenotime 2 type (Figure 3) shows the ν 1 mode at a decreased intensity and lower frequency (997 cm −1 ), while the ν 3 mode is completely diminished. Similar behavior was reported by Švecová et al. [52] for an altered xenotime affected by radiation damage with a high sub-micrometer porosity. An alteration state of xenotime 2 might explain the different shape of the spectrum. Moreover, the xenotime 2 sample has lower concentrations of U, Th and Y, most likely being leached due to its advanced degree of alteration.

Oxides
The oxides identified in the Valea Pianu gold deposit by means of Raman spectroscopy are rutile, anatase and cassiterite (Table S5). Titanium dioxides occur in nature as three polymorphs: the stable phase rutile, and metastable anatase and brookite (TiO 2 ). Rutile and anatase are both tetragonal with different space groups (Figures S8-S9). Tin dioxide has a rutile-like crystalline structure ( Figure S10) and it is known as cassiterite (SnO 2 ).
The Raman spectra of the studied oxides show well-resolved and intense peaks (Figure 4). These bands can be assigned in terms of factor group analysis. Therefore, rutile has four Raman-active modes with symmetries B 1g , E g , A 1g and B 2g , out of which the B 2g mode has a very low intensity, and it can appear as a weak band at 827 cm −1 [28,29,63]. These modes are associated with symmetric stretching (E g ), symmetric bending (B 1g ) and antisymmetric bending (A 1g ) vibrations of O-Ti-O bonds [28]. The B 1g mode is observed in rutile at 144 cm −1 , E g mode at 447 cm −1 , while the 612 cm −1 peak is assigned to the A 1g mode (Figure 4). An additional broad band appears in rutile at around 243 cm −1 , which might be attributed to disorder-induced scattering or the second-order effect [28]. The spectrum of anatase shows six Raman-active modes as 1A 1g , 2B 1g and 3E g symmetries. The bands at 142, 195 and 637 cm −1 are assigned to the E g modes, while the 394 cm −1 peak is attributed to the B 1g mode. The Raman line at 514 cm −1 in anatase is a doublet of the A 1g and B 1g vibrational modes [29,63,64].
Minerals 2020, 10, x FOR PEER REVIEW 12 of 18 different space groups (Figures S8-S9). Tin dioxide has a rutile-like crystalline structure ( Figure S10) and it is known as cassiterite (SnO2). The Raman spectra of the studied oxides show well-resolved and intense peaks (Figure 4). These bands can be assigned in terms of factor group analysis. Therefore, rutile has four Raman-active modes with symmetries B1g, Eg, A1g and B2g, out of which the B2g mode has a very low intensity, and it can appear as a weak band at 827 cm −1 [28,29,63]. These modes are associated with symmetric stretching (Eg), symmetric bending (B1g) and antisymmetric bending (A1g) vibrations of O-Ti-O bonds [28]. The B1g mode is observed in rutile at 144 cm −1 , Eg mode at 447 cm −1 , while the 612 cm −1 peak is assigned to the A1g mode ( Figure 4). An additional broad band appears in rutile at around 243 cm −1 , which might be attributed to disorder-induced scattering or the second-order effect [28]. The spectrum of anatase shows six Raman-active modes as 1A1g, 2B1g and 3Eg symmetries. The bands at 142, 195 and 637 cm −1 are assigned to the Eg modes, while the 394 cm −1 peak is attributed to the B1g mode. The Raman line at 514 cm −1 in anatase is a doublet of the A1g and B1g vibrational modes [29,63,64].
In the case of cassiterite, group theory predicts four Raman-active modes. Several Raman bands are observed (Figure 4), out of which the most intense one at 634 cm −1 is assigned to the A1g mode. The other modes appear at 125 cm −1 (B1g mode), 452 cm −1 (Eg mode), and 775 cm −1 (B2g mode) [27,65,66]. Two additional bands are observed at 499 and 833 cm −1 . The Raman lines at 243 and 286 cm −1 are attributed to the transverse (TO) and longitudinal optical (LO) modes, respectively [27].  In the case of cassiterite, group theory predicts four Raman-active modes. Several Raman bands are observed (Figure 4), out of which the most intense one at 634 cm −1 is assigned to the A 1g mode. The other modes appear at 125 cm −1 (B 1g mode), 452 cm −1 (E g mode), and 775 cm −1 (B 2g mode) [27,65,66]. Two additional bands are observed at 499 and 833 cm −1 . The Raman lines at 243 and 286 cm −1 are attributed to the transverse (TO) and longitudinal optical (LO) modes, respectively [27].

Titanite and Barite
Titanite (CaTiSiO 5 ) is a nesosilicate with a structure composed of SiO 4 tetrahedra, TiO 6 octahedra and CaO 7 polyhedra ( Figure S11). The Ti atoms are in an off-centered position within the octahedra. TiO 6 octahedra share corners, forming zig-zag chains that are linked by SiO 4 units. The entire framework contains large cavities that enclose the Ca atoms [67,68].
The Raman spectrum of titanite ( Figure 5) shows multiple bands that are rather broad and sometimes asymmetric. The band assignment considers the strong influence of the Ti-O bonds, which show a stretching mode at 608 cm −1 . The vibrational modes of SiO 4 tetrahedra are visible at 859 and 871 cm −1 for the symmetric stretching ν 1 and at 913 cm −1 for the ν 3 antisymmetric stretch. Three other intense lines appear at 423, 468 and 541 cm −1 that can be assigned to the ν 4 bending modes [69,70]. Up to 10 more bands appear in the low region of the spectrum, which are attributed to the external modes ( Figure 5).

Titanite and Barite
Titanite (CaTiSiO5) is a nesosilicate with a structure composed of SiO4 tetrahedra, TiO6 octahedra and CaO7 polyhedra ( Figure S11). The Ti atoms are in an off-centered position within the octahedra. TiO6 octahedra share corners, forming zig-zag chains that are linked by SiO4 units. The entire framework contains large cavities that enclose the Ca atoms [67,68].
The Raman spectrum of titanite ( Figure 5) shows multiple bands that are rather broad and sometimes asymmetric. The band assignment considers the strong influence of the Ti-O bonds, which show a stretching mode at 608 cm −1 . The vibrational modes of SiO4 tetrahedra are visible at 859 and 871 cm −1 for the symmetric stretching ν1 and at 913 cm −1 for the ν3 antisymmetric stretch. Three other intense lines appear at 423, 468 and 541 cm −1 that can be assigned to the ν4 bending modes [69,70]. Up to 10 more bands appear in the low region of the spectrum, which are attributed to the external modes ( Figure 5).
As in the case of zircon, titanite can be affected by metamictization due to its relatively high contents of REE, U, and Th. Zhang et al. [68] showed that the Raman spectrum of titanite is strongly dependent on the degree of metamictization and it changes by broadening peaks, shifting frequencies, decreasing overall intensity and new additional bands. Barite (BaSO4) is a sulfate mineral with SO4 tetrahedra as fundamental unit that surround the Ba 2+ cations ( Figure S12). The Raman spectrum of barite ( Figure 5) shows very clear and sharp bands As in the case of zircon, titanite can be affected by metamictization due to its relatively high contents of REE, U, and Th. Zhang et al. [68] showed that the Raman spectrum of titanite is strongly dependent on the degree of metamictization and it changes by broadening peaks, shifting frequencies, decreasing overall intensity and new additional bands.
Barite (BaSO 4 ) is a sulfate mineral with SO 4 tetrahedra as fundamental unit that surround the Ba 2+ cations ( Figure S12). The Raman spectrum of barite ( Figure 5) shows very clear and sharp bands that can be attributed to the four fundamental vibrational modes of SO 4 tetrahedra. Therefore, the most intense line at 987 cm −1 is assigned to the ν 1 stretching mode and the very weak one at 1142 cm −1 is due to the ν 3 mode. The bending vibrations occur at 461 cm −1 for the ν 2 symmetric type and at 616 and 628 cm −1 for the ν 4 antisymmetric modes [22,23].

Conclusions
The micro-Raman technique proved to be a powerful tool for identifying a large number of heavy minerals in a typical gold placer deposit from Pianu Valley, Romania. The following minerals were observed: garnet, kyanite, staurolite, zircon, allanite, monazite, xenotime, rutile, anatase, cassiterite, titanite and barite. The results were confirmed by EPMA chemical analysis. Even if Raman spectroscopy can be used as a stand-alone technique for most of the minerals, in some cases it is recommended to use additional techniques for extra caution in structural and chemical interpretation of the minerals.
Systematic shifting of the Raman bands is very sensitive to structural deformations, even at a low degree of cation substitutions, especially for solid solution series. Such shifting behavior was observed in garnets, kyanite, staurolite and zircon. Several Raman bands in garnets spectra are strongly correlated with the average ionic radius in the X 2+ site and therefore peaks positions are dependent to the molar ratios of each end-member. The chemical composition of intermediate members can be estimated by considering up to six peaks with unimodal behavior only by means of Raman spectroscopy.
The distinction between polymorphic minerals can be realized only by structural investigations, micro-Raman being the fastest way to obtain such information. This approach was successfully used complementary to chemical analyses in order to differentiate between kyanite-sillimanite-andalusite and titanium dioxide polymorphs (rutile-anatase-brookite).
Metamictization processes can affect some important accessory minerals due to their relatively high contents of U and Th radioactive elements. The degree of radiation damage can be assessed by evaluating the Raman spectra. Features like broad peaks, frequency shifting and lower intensities are directly correlated with the metamictization degree in zircon, titanite and monazite, even if the latter is more resistant to radioactive damage.