Geochemical Transactions

Fluid inclusions and geological relationships indicate that rodingite formation in the Asbestos ophiolite, Québec, occurred in two, or possibly three, separate episodes during thrusting of the ophiolite onto the Laurentian margin, and that it involved three fluids. The first episode of rodingitization, which affected diorite, occurred at temperatures of between 290 and 360°C and pressures of 2.5 to 4.5 kbar, and the second episode, which affected granite and slate, occurred at temperatures of between 325 and 400°C and pressures less than 3 kbar. The fluids responsible for these episodes of alteration were moderately to strongly saline (~1.5 to 6.3 m eq. NaCl), rich in divalent cations and contained appreciable methane. A possible third episode of alteration is suggested by primary fluid inclusions in vesuvianite-rich bodies and secondary inclusions in other types of rodingite, with significantly lower trapping temperatures, salinity and methane content. The association of the aqueous fluids with hydrocarbon-rich fluids containing CH4 and higher order alkanes, but no CO2, suggests strongly that the former originated from the serpentinites. The similarities in the composition of the fluids in all rock types indicate that the ophiolite had already been thrust onto the slates when rodingitization occurred.

Structural studies using SEM and TEM indicate the local-disordering of the structure as a result of charge-transfer process between the inserted lithium and the molybdenite layer. Selected area electron diffraction measurements indicate the large variations in the diffusivity of lithium confirming that the charge-transfer is different along and perpendicular to the layers in molybdenite. Thermal heating of molybenite surface in air at 400°C induces surface oxidation, which is slow during the first hour of heating and then increases significantly. The SEM results indicate that the crystals formed on the molybdenite surface as a result of thermal oxidation exhibit regular thin-elongated shape. The average size and density of the crystals on the surface is dependent on the time of annealing; smaller size and high density during the first one-hour and significant increase in size associated with a decrease in density with further annealing.

Background
Sulfide minerals and the associated geological/physical/ chemical processes are an active research topic for mineralogists, geochemists, and geotechnical/environmental engineers. Sulfide ores constitute a major source of metals, especially noble metals. In addition, the rich diversity in crystal chemistry, surface reactivity, phase transformations, stability, thermodynamics, and electronic properties makes the sulfide minerals attractive for a wide variety of industrial applications, such as lubricants and catalysts [1,2]. Molybdenite (MoS 2 ) belongs to the family of the transition-metal dichalcogenide (TMD) minerals with the formula MX 2 (where M = Cd, Ti, Mo, Sn and X = I, S, Se). Due to their layered structure, TMDs are often referred to as two-dimensional (2D) solids. The 2D structure of these minerals is due to the strong covalent or ionic bonding within a layer while individual layers are held together by weaker van-der-Waals forces. Even though the latter are often referred to as "van der Waals" type of interactions, some contributions from covalent and ionic interactions are also possible, particularly in the metal inserted complexes [3,4]. These compounds exhibit anisotropic physical properties, such as different conductivity parallel and perpendicular to the layers, ranking MoS 2 the most anisotropic 2D material after graphite [3].
Molybdenite is the essential ore mineral of the molybdenum industry for production of Mo metal and Mo-based compounds, such as sodium and calcium molybdates, ammonium paramolybdate, and molybdenum trioxide [5]. Molybdenite-based formulations are extensively used in industrial machinery and weapons for lubrication [6]. Molybdenite is a widely used catalyst in CO hydrogenation and hydrodesulfurization processes for the production of cleaner fuels [7,8]. Mo is a redox-sensitive trace metal and becomes enriched in sulfidic, reducing, and organic rich sediments. Molybdenite is found in deepsealed veins associated with scheelite, wolframite, topaz, and fluorite [9].
The mineral molybdenite belongs to the group VI TMD compounds that adopt the 2H b polytype, in which the metal atoms are staggered. The perspective view of the structure of MoS 2 is shown in Fig. 1 [9]. Naturally occurring molybdenite has been found to exist mostly in 2H stacking polytypes but 3R ones are also available.
The objective of the present work is to understand the effects of thermally-induced oxidation and alkali metal (Li) ion insertion in molybdenite surfaces. Understanding the oxidation behavior of MoS 2 is crucial because several applications of MoS 2 are influenced by its degree of oxidation. Similarly, information on the mechanism of foreign metal-atom insertion into the molybdenite structure is important to understand the associated electronics structure changes. Therefore, the present investigation was made using a wide variety of analytical techniques such as atomic force microscopy (AFM), scanning electron microscopy (SEM), transmission electron microscopy (TEM), Raman spectroscopy (RS), Rutherford backscattering spectrometry (RBS), and nuclear reaction analysis (NRA) to elucidate how the oxidation proceeds in these minerals at higher temperatures. In addition, the structural effects associated with metal insertion at room temperature were evaluated. Our choice of lithium (Li) was due to the following reasons: (1) The chemistry of Li, being a simple metal, is well known, (2) the metal insertion in molybdenite is only possible with small, strongly reducing guests, such as the alkali metals [3,10] and (3) comparison of the results will be made easy in view of several existing reports on Li in synthetic molybdenite [3,10]. Using the familiar surface analytical techniques, namely AFM, Perspective view of the structure of 2H b -MoS 2 Figure 1 Perspective view of the structure of 2H b -MoS 2 . The crystallographic axes for the H structure are those for standard hexagonal systems with a(x,y) parallel and c(z) perpendicular to the layers. The coordination of Mo atoms are trigonal prismatic with the stacking sequence AbA-BaB along the c direction.

S
Mo 1120 _ c a b SEM and TEM, is to probe the changes in surface structure and morphological features due to thermal oxidation and metal-ion insertion in molybdenite. The combined use of RS and TEM measurements allow probing the information on the local chemical structure and bonding and are quite useful to evaluate the impact of metal-ions on the local structural environment in molybdenite. Using the less familiar ion-beam analytical techniques, such and NRA, is mainly to understand the thermal oxidation behavior of molybdenite. High-energy ion-beam analysis using nuclear reactions, particularly NRA measurements, provide an important tool to measure the absolute concentration of lighter elements, such as hydrogen (H), carbon (C), nitrogen (N), and oxygen (O), and makes it quite useful to understand the science and processes at geochemical media [11][12][13]. Therefore, we have used the NRA to probe the surface oxide growth on molybdenite surface.

Experimental
Two sets of naturally-occurring molybdenite samples were used in this work. The minerals, found in the mineral storage/collection at the Department of Geological Sciences at the University of Michigan, of unknown origin were employed for thermal oxidation experiments. Molybdenite samples obtained from Japan were employed for lithium-insertion measurements. Atomic force microscopy measurements were performed using a Digital Instruments AFM (NanoScope IV). The measurements were made in tapping mode. Scanning electron microscopy observations were made using two different instruments. A JEOL (model 6150) electron microscope was used to image the oxidation behavior of the molybdenite surfaces as a function of heating time. SEM imaging experiments on lithium-reacted molybdenite samples were made using a high-resolution electron microscope (Hitachi S-4700). Raman scattering (RS) spectra were measured using a Jobin-Yvon U1000 double-pass spectrometer equipped with a cooled, low-noise photomultiplier tube (ITT FW130). The incident light used for the experiments was the 515 nm Ar ion laser. TEM analysis was performed using a JEOL TEM 2010F at a 200 kV acceleration voltage. Phase and structure of the material were monitored using selected area electron diffraction (SAED

Pristine molybdenite surfaces
The AFM image of a representative pristine natural MoS 2 sample surface is shown in Fig. 2. Molybdenite surfaces are often atomically flat over relatively large areas as seen in this work (Fig. 2) and other reports [14][15][16]. AFM images, in the area of scanning, show no specific surface features such as steps or kinks. The electron diffraction pattern of the molybdenite crystal is shown in (Fig. 3).
The hexagonal lattice of the 2H polytype is evident from the image (Fig. 3).
The long-wavelength lattice vibrations of molybdenite are interesting and often provide a model system to understand the linear-chain model of the layered-type chalcogenide minerals. Group theory predicts two infrared (IR) and four Raman-active modes for 2H-MoS 2 [16][17][18]. The Raman scattering spectrum of the pristine molybdenite surface recorded at ambient temperature is shown in Fig.  4 along with the associated mode assignments. The A 1g mode at 407 cm -1 is an intralayer mode involving the motion of S atoms along the c axis. The E g mode at 383 cm -1 is an intralayer vibrational mode involving motion of Mo+S atoms in the basal plane. The peak at 286 cm -1 with weak intensity and E 1g symmetry is observed, which involves S atoms in the basal plane. The spectrum of the pure 2H-MoS 2 exhibits a rigid-layer (RL) mode E g at 32 cm -1 . This mode is of interlayer type involving rigid motion of neighboring sandwiches in anti-phase. In layered crystals, the frequency of the RL mode provides direct information on the strength of the interlayer forces in these crystals, since for such RL motions the restoring forces are provided entirely by layer-layer interactions [19]. Figure 5 shows the SEM and electron diffraction micrographs of Li-inserted molybdenite produced by inserting the samples in n-butyl-lithium. It can be seen that some defects are created near the edges and in the steps of the specimen, after 0.2 hours of insertion (arrows in Fig. 5a). After 2 hours of insertion reaction, superlattice spots appear (label s in Fig. 5b). Notice also the splitting of the main spots denoted by the letter m. The micrograph taken from the same area reveals that the specimen is heavily disordered or contains a significant number of defects as a result of metal-insertion reaction in the molybdenite structure (Fig. 5c). The distribution of lithium vs. distance from the edges of the specimen has been determined by tracer studies using solid-state nuclear track detectors (SSNTD). A typical SSNTD image of Li 0.3 MoS 2 is also shown in Fig. 5 (Fig. 5d). Continued insertion reaction for AFM image of the molybdenite surface 20 hours results in the appearance of superlattice spots in areas far from the edges. Splitting of the ( ) spots still occurs, the inner spots being very intense while the outer ones are faint as shown in Fig. 5 (Fig. 5e). This is an indication that the intercalated system did not reach saturation, even locally, although the diffusion of lithium has progressed, implying large variations in the diffusivity of lithium in molybdenite.

Lithium insertion
The results of the lithium insertion reaction and associated structural modifications are discussed as follows. Insertion of foreign metal and/or ionic species into a host structure implies two main conditions: (1) the structural ability to accept ions in empty sites, and (2) the presence of acceptor levels in the electronic structure. From a crystal structure point of view, there are several possible locations to insert lithium into molybdenite (Fig. 6) as the unit cell of 2H b -MX 2 possesses tetrahedral and octahedral intersandwich sites (this is equally applicable for inserting any other metal ions into the structure). Among the group-VI

RS intensity (cps)
Raman shift (cm -1 ) Electron micrographs of Li-inserted MoS 2 by the n-butyllithium method TMDs, MoS 2 is a typical example where insertion reactions induce a local structure modification. In this particular case, Mo presents a trigonal prismatic (TP) S coordination which changes to the octahedral (Oh) one, i.e. the TP → Oh transition [3,[20][21][22]. The structure modification is accompanied by an increase in the Mo-S bond ionicity in agreement with the respective stability of the new atomic arrangement, the Coulomb repulsion between partially charged ligands favoring the octahedral form. Also, comparison of the d-band density of states for 2H-MoS 2 and hypothetical 1T-LiMoS 2 shows that the occupied bands which contain six states are lower in the case of the Oh phase corresponding to the glide process between Mo and S atoms. This is a good example of destabilization through lithium reduction. The transformation from TP to Oh coordination is attributed to a process which is driven by a lowering of the electronic energy for the octahedral structure when the charge transfer occurs from Li ions to the MoS 2 layer (electrons are donated by Li) upon insertion reaction [3,10].

Thermal oxidation
The NRA results are shown in Fig. 7. The curves are shown as a function of heating time of molybdenite surfaces in air. The molybdenite surface, introduced into the highvacuum chamber immediately after cleavage, did not show any peak corresponding to oxygen. This observation indicates that there is no presence of oxygen on the pristine molybdenite surface before the thermal oxidation experiments. We claim that we start our thermal oxidation experiments with atomically clean molybdenite surfaces since the detection limit of NRA to determine the oxygen concentration (using the 16  heating. This observation indicates that the oxidation during the first hour of heating is slow leading to oxidation only on the surface. A significant increase in peak intensity and change in shape (becoming more broad) with further heating after the first hour indicates that oxidation is progressively proceeding from the surface layers into the deeper layers.
SEM micrographs of the molybdenite surfaces are shown in Fig. 8 as a function of heating time. The thermallyinduced oxidation completely changes the morphology of the molybdenite surface. The SEM images indicate the formation of small oxide crystals on the surface within the first hour of heating. More significant changes are observed after 2 to 3 hours of heating time. The crystals formed are light greenish indicating the oxide phases. The number of the oxide crystals formed on the surface increases significantly as a result of thermal oxidation. The crystals formed on the molybdenite surface as a result of thermal oxidation are assuming a regular shape with an elongated thin structure, which is in good agreement with the shapes reported for Mo oxides [23]. Oxide crystals formed on the surface as a result thermally-induced oxidation during the first one hour of heating are, however, longer in size extending over 4-5 microns in length and about 2 microns in width. The SEM image shown for samples after 3 hours of heating clearly shows the density of such crystals formed on the surface increases significantly. The increase in the number of crystals at 3 hours of heating increases the surface roughness with random distribution of the crystals in addition to the irregular shapes. However, an interesting phenomenon is the size of the crystals. It can be seen that the crystals formed at 3 hours of heating are smaller in size when compared to the ones formed in the first hour of heating. It may be that the tendency of all crystals to grow within the limited surface area may inhibit the growth of larger crystals in the neighboring medium. The accumulation of several crystals contributing to the rough morphology is also observed in some areas, which supports the idea that the higher density of oxide crystals may inhibit the growth of large and uniform crystals as seen at the initial stages of oxidation. In other words, the formation of continuous oxide layers, initially formed on the top-most surface layers, may influence the lateral diffusion of the particles to join the already formed crystals and to grow into larger crystals.

Conclusion
The thermally-induced oxidation behavior and lithium metal-insertion reaction effects on molybdenite surfaces are studied using a combination of structural imaging and spectroscopic measurements. Starting from pristine and well-ordered molybdenite layers, local disordering of the structure occurs as a result of charge-transfer process between the inserted metal ions and the molybdenite layer. Heating of the molybdenite surface in air at 400°C induces oxidation. Oxidation in the first hour of heating is initially slow and then increases significantly.

Authors' contributions
CVR carried out sample collection, performed the measurements of SEM, AFM and NRA on the prestine and oxidized MoS 2 samples, interpreted data, and drafted the manuscript, UB is involved in planning, coordination and support for all the measurements and assisted with data interpretation and drafting the manuscript, VS assisted with the preparation of samples, experiments and data collection, and interpretation of NRA results obtained, CMJ assisted with sample collection and performed the lithium insertion, Raman spectroscopy and TEM measurements and participated in drafting the manuscript. All authors have approved the significance of the work, interpretation of results, and contents of the final manuscript.
SEM images of molybdenite surfaces as a function of heating time Figure 8 SEM images of molybdenite surfaces as a function of heating time.