Investigating the nano-tribological properties of chemical vapor deposition-grown single layer graphene on SiO2 substrates annealed in ambient air

Quan Wang*ab, Bing Baia, Yun Lia, Yan Jiangc, Laipeng Mad and Naifei Rena
aSchool of Mechanical Engineering, Jiangsu University, Zhenjiang 212013, P.R. China. E-mail: wangq@mail.ujs.edu.cn
bState Key Laboratory of Transducer Technology, Chinese Academy of Sciences, Shanghai 200050, P.R. China
cSchool of Material Science & Engineering, Jiangsu University, Zhenjiang 212013, P.R. China
dShenyang National Laboratory for Materials Science, Institute of Metal Research, Chinese Academy of Sciences, Shenyang 110016, P.R. China

Received 15th October 2014 , Accepted 6th January 2015

First published on 6th January 2015


Abstract

Nano-tribological properties of graphene have attracted a lot of research interest in the last few years. In this work, X-ray photoelectron spectroscopy was used to study the distribution of chemical groups in chemical vapor deposition (CVD) of single layer graphene transferred onto a 90 nm SiO2/Si substrate. It was demonstrated that the graphene was oxidized after thermal treatment at 520 °C in ambient air, as indicated by the formation of C[double bond, length as m-dash]O and C–OH bonds. Significantly enhanced D-band and decreased 2D-band were found in the Raman spectrum. Blue shift occurred for the G-band and 2D-band after thermal oxidation. The nano-tribological properties of graphene before and after thermal oxidation were studied with atomic force microscopy. A remarkable increase of friction was found on the surface of graphene after thermal oxidation. This was contributed to by the increased adhesion and decreased Young's modulus of the functionalized graphene, inducing the increase of contact stiffness. In addition, the adhesion force between the tip and the samples was discussed as an important factor affecting tribological properties on the nano-scale.


1 Introduction

Graphene is a two-dimensional, zero bandgap material with carbon atoms arranged in a honeycomb lattice.1–3 Since graphene first became experimentally accessible in 2004,4 it has attracted considerable research attention owing to its unique properties such as high specific surface area, ultra-high intrinsic mobility and great mechanical strength.5–7

Many recent works have concentrated on the mechanical and tribological performance of graphene.8–10 Lee et al. found that there was a decreased friction force between SiN tip and graphene flake prepared on silicon oxide with the increasing number of graphene layers. The “puckering effect” was proposed to explain this phenomenon.9 Filleter et al. found that friction force was higher on graphite compared to epitaxial bilayer graphene on SiC substrate surface, and “electron-phonon coupling” theory was presented in their work.10 However, the large area graphene with controlled shape and layers is needed in the practical applications of graphene. Currently, the increasingly mature chemical vapor deposition (CVD) method was widely applied to the preparation of graphene which is utilized in graphene-based devices and the flexible display devices. Graphene prepared by CVD is sensitive to the impurities and ion in air ambient, such as the adsorption of H2O, oxygen molecules and the chemical functionalization of graphene surface in air ambient. It is not clear how doping and functionalization influence the nano-tribological properties of graphene. Therefore, it must be better understood to design the stability of graphene-based devices. Nan et al. studied the thermal stability of graphene in air ambient by Raman spectroscopy, and found that single layer graphene (SLG) prepared by CVD would show defects at ∼500 °C. These defects were initially SP3 type and became vacancy-like at higher temperature.11 However, there are few research on the change of bonding energy structure and mechanics characteristics of graphene after thermal treatment in ambient. In this paper, we present the research on the nano-tribological of graphene before and after thermal treatment in ambient.

In the term of device integration, the preparation of devices with thin gate dielectric layer has become preferred with the rapid development of CMOS process. We use the high quality SLG prepared by CVD method transferred onto Si substrate with 90 nm thermal oxidation SiO2 in this study. X-ray photoelectron spectroscopy (XPS) was used to analyze the distribution of chemical groups in the CVD single layer graphene. It revealed that the graphene was oxidized, which was indicated by the formation of C[double bond, length as m-dash]O, C–OH bond after annealing in ambient atmosphere at 520 °C. This can be manifested as the appearance of a disorder-induced D-band in Raman spectrum. Then the surface topography and nano-tribological properties of graphene before and after thermal oxidation were investigated with atomic force microscopy (AFM). These results would provide valuable information for the application of graphene-based devices.

2 Experimental

Graphene samples are prepared by CVD method. Graphene is grown on thin copper foils about 25 μm at 1000 °C with a flow of H2 (30 sccm) and CH4 (50 sccm), the growth pressure of ∼100 Pa, and the growth time of 20 min. Graphene was transferred to the Si substrate with 90 nm thermal oxidation SiO2. Optical microscopy is used to estimate the thickness, which is further confirmed by Raman spectroscopy.

As the physical or chemical adsorption of H2O molecules on sample surface has great influence on surface properties of two-dimensional material, the samples were annealed in rapid thermal annealing furnace in N2 atmosphere with the temperature at ∼380 °C for 2 hours. Thermal oxidation process of graphene was carried out in STF-1200X tube furnace and the graphene samples were kept in it for 5 min at 520 °C in air ambient.

Raman spectra of graphene before and after thermal oxidation was recorded at room temperature by using WiTec alpha 300R system with 532 nm laser wavelength. The laser power was kept below 2 mW to avoid laser-induced heating or damage for graphene.12 The 5 μm × 5 μm Raman images was obtained by recording Raman spectrum at every point when the sample moves along the x or y axis with step size of 200 nm. X-ray photoelectron spectroscopy (XPS, K-Alpha, Thermo Scientific) was performed using amonochromated Al Kα irradiation. The chamber pressure was about 3 × 10−8 Torr under testing conditions. Peak split and fitting of C 1s of graphene samples before and after thermal oxidation were accomplished using XPSPEAK 4.1.

The surface morphology and nano-tribological properties of different samples were measured using a MFP3D-SA atomic force microscope (AFM, Asylum Research) in contact mode. A non-conductive silicon nitride probe (Veeco NP-20, spring constant k = 0.37 N m−1) was employed for AFM measurements with scan rate of 1 Hz and scan angle of 90°. There are various influence factors in the measurement of nano-tribological, such as roughness of substrate,13 scanning speed14 etc. In order to improve the reliability of the experimental results, parameter settings (scan size, scan rate, scan angle et al.), type of probe and test environment (room environment) were kept the same in all measurements. The values of friction presented below are average value of Root-Mean-Square (RMS) obtained from scan area about 1 μm × 1 μm. To obtain the adhesive force between the AFM tip and the samples surface, the force-distance curves were recorded, and the pull-off force was considered to be the adhesive force. The value of adhesion force was the average of RMS obtained from the 10 pt × 10 pt force-distance map. Owing to the complexity of the calibration of the torsional force constant, and the same cantilever applied in all measurements, the output voltages were directly used as the relative frictional force.

3 Results and discussion

In order to avoid the influence of the number of graphene layers on the measurements of nano-tribological properties, Raman spectrum was used to be an effective means for the precise representation of the number of layers of the test area. Fig. 1(b) shows the Raman spectra of pristine graphene which is the randomly selected test point in Fig. 1(a), and two representative peaks of graphene observed at ∼2680 cm−1 (2D-band) and ∼1580 cm−1 (G-band). The intensity ratio of 2D and G peaks (I2D/IG) is >1.5, indicating that the graphene is single layer graphene (I2D/IG of bilayer and multilayer would be less than 1.5). The insert shows a single Lorentzian 2D peak with a full width at half maximum of ∼30 cm−1 which confirms that the graphene is a SLG.15 Fig. 1(a) presents the Raman image of the intensity ratio of 2D and G peaks for 5 μm × 5 μm test area of CVD graphene. Most of the region is SLG with I2D/IG >1.5. There is no D-band observed in Fig. 1(b), indicating that there is a small amount of intrinsic defects in these pristine CVD graphene.
image file: c4ra12437e-f1.tif
Fig. 1 Raman spectrum of pristine SLG with small amount of intrinsic defects. (a) is the Raman image of the intensity ratio of 2D and G peaks (I2D/IG). (b) is the Raman spectrum fitted by Lorentz. The insert gives the Lorentz fitting of 2D peak.

The Raman spectra of SLG before and after thermal treatment at 520 °C in air are shown in Fig. 2(b). The significant enhancement of D-band can be observed clearly due to the reaction of graphene with O2 in the open furnace, which creates defects such as carbon–oxygen SP3 bonds and vacancies after thermal treatment. Significant degradation and blue shift of 2D and G peaks are observed after thermal treatment in air, which is caused by the increasing defects and hole doping effect resulting from the charge transfer between adsorbed molecules (H2O and O2) and graphene. Fig. 2(a) shows the Raman image of intensity ratio of D and G peaks (ID/IG) is 0.2–0.4, indicating the induced defects are uniformly distributed after thermal treatment.16,17


image file: c4ra12437e-f2.tif
Fig. 2 (a) is the Raman image of the intensity ratio of D and G band (ID/IG) of thermal treated graphene. (b) is the Raman spectra of SLG by CVD before and after thermal treatment at 520 °C in air ambient.

In order to determine the defect types of graphene and distribution of bond energy before and after thermal treatment in air, XPS analysis was carried out to substantiate directly whether the graphene was oxidized or not after thermal treatment in air ambient. Fig. 3(a) shows the C 1s of pristine graphene. A typical signal at C–C (284.5 eV) was obtained through Gaussian fitting. Fig. 3(b) shows XPS spectrum of C 1s for graphene after thermal treatment in air can be deconvoluted into three Gaussian peaks, which are typical signals of C–C (284.5 eV), C–OH (285.6 eV), and C[double bond, length as m-dash]O (288.4 eV). XPS analysis fully demonstrates that graphene was oxidized faintly with generating of SP3 hybridization after thermal treatment at 520 °C in air ambient.


image file: c4ra12437e-f3.tif
Fig. 3 XPS spectra of C 1s for graphene (a) before and (b) after thermal treatment at 520 °C in ambient atmosphere, the ratio between C[double bond, length as m-dash]C, C–OH, and C[double bond, length as m-dash]O is 10[thin space (1/6-em)]:[thin space (1/6-em)]3[thin space (1/6-em)]:[thin space (1/6-em)]1.

Fig. 4(a) presents the morphology of pristine graphene on 90 nm SiO2/Si. The ripples, originating from the cool-down and transfer process of the sample preparation, are shown in the morphology. Moreover, a small quantity of residuals of the PMMA resist appear as highly protruding objects. Fig. 4(b) and (c) present the trace and retrace image of friction. To identify graphene on substrate, friction image is more effective compared with morphology image. Shown in Fig. 4(d), the friction loop is consisted of trace and retrace friction force profiles of the line. The width of the loop can be used to roughly represent the practical friction of different area, indicating that friction on SiO2 is about four times larger than that on graphene.


image file: c4ra12437e-f4.tif
Fig. 4 (a) is the morphology image of pristine graphene on 90 nm SiO2/Si. (b) and (c) are the trace and retrace image of friction. (d) is friction value line profiles obtained from forward and backward scans across SiO2 and SLG. The position of scan line is shown in (a) (green line), (b) (black line) and (c) (red line), respectively. The scan size is 10 μm × 10 μm and the color bar used in (b) and (c) is the same one.

Compared with Fig. 4(a) and 5(a) gives a more micromesh morphology image of the wrinkles of pristine graphene clung gently on the substrate. The wrinkles would be preserved if they existed before the graphene was transferred to the target substrate.18 In addition, SLG, clung on the substrate by the weak van der Waals force interaction with substrate, could also exist folded structure.19 Moreover, residuals between graphene and substrate in transferred process is a main factor for the influence of topography of graphene on substrate. The presence of the winkles and bubbles also caused the change of friction between tip and sample with fluctuation of the surface of sample, indicating nano-tribological test is associated with the roughness of the surface, as is shown in Fig. 5(b) and (c). The surface become smoother after thermal treatment in air, supported as the values of mean square roughness values (Ra) from 0.85 nm to 0.44 nm. The flexible PMMA on graphene become stiffness after thermal treatment, resulting in a larger friction, shown as the light and black spot in Fig. 5(e) and (f), respectively. In order to minimize the influence of changing morphology on the test of nano-tribological, all of the friction data presented below are average values of Root-Mean-Square (RMS) obtained from scan area about 1 μm × 1 μm.


image file: c4ra12437e-f5.tif
Fig. 5 (a) and (d) are the morphology images of graphene on SiO2/Si substrate before and after thermal oxidation respectively. Followings are the friction images of (b) (e) trace and (c) (f) retrace direction corresponding to (a) and (b), respectively.

Fig. 6 shows the changes of friction between tip and graphene versus varied load. The relative friction coefficients (RFC) of friction for graphene before and after oxidized can be extracted from the slopes of the fitting lines in Fig. 6. There is an obvious enhancement of RFC of the graphene after oxidized, observed in the Fig. 6. It is noted that the friction of both samples could not be zero despite without normal load. According to the Amonton's law, the friction between tip and sample can be expressed as follow in nanoscale:20

 
Lf = μ(Lload + Ladh), (1)
where μ is friction coefficient. Lload and Ladh are the normal load and the adhesion force between the two objects. Due to the biggish load LloadLadh, the effect of adhesion force on friction can be neglected in the measurements of macroscopical friction. But Lload (about 0–30 nN in this research) is equal to Ladh (about 0–100 nN with our measurements) in the measurements in nano-tribological. So Ladh is crucial for the characterization of nano-tribological, especially the two-dimensional material with excellent elastic properties. In this work, when Lload = 0, Ladh plays the role of Lload, leading to the so-called “initial friction”.


image file: c4ra12437e-f6.tif
Fig. 6 A graph showing the effect of varying load on the friction force (measured as the width of the lateral deflection loop). Error bars represent the standard deviation of collected data.

Fig. 7 shows the statistics histogram of friction signals and adhesion force for pristine graphene, oxygenated graphene and 90 nm SiO2/Si substrate. There is a significant increase of friction and adhesion on graphene after oxidation. But it is much smaller compared to the substrate, indicating that graphene can be an excellent solid lubricant whether it is oxidized or not. The corresponding relationship between the change of friction and adhesion suggests that adhesion is a decisive factor in the nano-tribological.


image file: c4ra12437e-f7.tif
Fig. 7 Comparison of friction signals and adhesion force for pristine graphene, oxygenated graphene and 90 nm SiO2/Si substrate (normal load = 7 nN). Error bars represent the standard deviation of collected data.

There are several possible mechanisms proposed for the nano-tribological behavior of graphene, such as the energy dissipation mechanism,21 the mechanism of electron phonon coupling9 and “puckering effect”.8,22 Due to the inseparable relationship between friction in nano-scale and the elastic property of material,23 the exchange of mechanical property with graphene after chemical modification cannot be neglected on the impact of its tribological property.

The nanoscale friction force in the lateral force microscope (LFM) stick-slip motion is determined by the slope of the LFM profile (the total lateral stiffness) on the measurements of nano-tribological property of graphene clung on the substrate with AFM.24 The relationship between total lateral stiffness is and the friction signal (Torsion of the probe) can be expressed as25

 
image file: c4ra12437e-t1.tif(2)

kgraphene-substrate is large enough to be neglected in eqn (2) when graphene and substrate interact with each other tightly, and klever is not be considered because the same lever was used in this experiment. The contact stiffness determined by the lateral deformation in the tip-to-graphene contact as the sole variable in eqn (2) can be expressed as

 
ktip-graphene = 8G*a, (3)
where G* can be calculated by G* = [(2 − ν1)2/G1 + (2 − ν2)2/G2]−1. G1 and G2, ν1 and ν2 are shear modulus and Poisson ratio of silicon nitride tip and graphene, respectively. a is the contact radius between tip and sample which can be expressed with Hertz theory:
 
image file: c4ra12437e-t2.tif(4)
where E* = [(1 − ν1)2/E1 + (1 − ν2)2/E2]−1. E1 and E2 are the Young's modulus (elasticity modulus) of silicon nitride tip and graphene respectively. R is the radii of curvature (ROC) for AFM probe tip, Lload is the normal load applied by the tip and Ladh is the adhesion force between tip and sample. Young's modulus of pristine SLG prepared by mechanical exfoliation is ∼1.0 TPa6 and SLG prepared by chemical reduction of graphene oxide is ∼0.25 ± 0.15 TPa.26 Paci et al. have proved that the oxidation of graphene will lead to a greatly reduce of its Young's modulus and breaking strength.27 Due to the decreasing of Young's modulus of graphene and the increasing of adhesion force between tip and sample, with the consideration of eqn (2)(4), there would be a larger total lateral stiffness ktotal when the tip slip on the graphene after oxidized. Therefore, an increasing friction signal can be observed by AFM.

As discussed above, the adhesion between tip and sample is a significant factor in nano-tribological tests with AFM. In other words, it also can be considered as an additional load. It is well-known that environment change for a short period of time is not easy to cause the functionalization of surface. But there would be a strong impact on the chemical bond on the surface of graphene because of the interaction with the gas and ion in the air. Therefore, it is imperative to know the influence of environment in the adhesion force on the graphene surface. Q. Li et al. found that the adhesion force between tip and mechanical exfoliated multilayer graphene was increased with the increase of the time exposed to the air, resulting a larger friction on the surface of graphene.28

The adhesion force is generally consists of capillary force LC, van der Waals force LvdW, electrostatic force LE and bonding force LB which can be expressed as

 
Ladh = LC + LvdW + LE + LB. (5)

Because the samples and the probes have placed in air for a long time, there is no additional charge on the surface. The surfaces of the tip and samples are saturated with chemical bonds, LE and LB can be neglected in this research. Herein, LC and LvdW are the main factors in the influence of adhesion between tip and graphene before and after thermal oxidation.

LC is a main contributor to adhesion and closely related to surface wettability. SP3 hybridization was introduced in pristine graphene after thermal oxidation and C–OH bonds were formatted on the surface of graphene. Due to the function of hydrogen bond, water vapor in the air is easy to be adsorbed on the surface, leading to the enhancement of capillarity between tip and sample. In addition, the contact area between graphene and substrate was increased after thermal oxidation (winkles and bubbles were decreased). The van der Waals force is increased because the interaction between tip and substrate is enhanced, as the Lennard-Jones potential strengthen with λ−6, where λ is the distance between tip and substrate. The van der Waals interaction between tip and SiO2 is much stronger than that of graphene. And hence, large LC and LvdW may be the key factor in the relatively high adhesion as is shown in Fig. 7.

4 Conclusions

In this work, high-quality single layer graphene prepared by CVD were transferred onto 90 nm SiO2/Si substrate. Thermal treatment for the transferred graphene was taken at 520 °C in air ambient. The microstructure and tribological properties of the samples were investigated. These results show that the graphene were oxidized, as indicated by the formation of C[double bond, length as m-dash]O, C–OH bond after annealed in air atmosphere. A significant increase of friction was found on the surface of graphene which is attributed to the decreased Young's modulus and increased adhesion after thermal oxidation. The adhesion force between tip and samples was discussed as an important factor on affecting of tribological properties in nano-scale.

Acknowledgements

This work was supported by the National Natural Science Foundation of China (no. 51175238 and 11102075), the Industrial Supporting Plan of Zhenjiang City (no. GY2013011), "Six talent peaks" of high level talent selection and training project of Jiangsu Province (no. 2013-ZBZZ-031), the Priority Academic Program Development of Jiangsu Higher Education Institutions and Graduate research and innovation projects in Jiangsu University (no. KYXX_0029). The authors would like to thank Rong Ding, from Taizhou Sunano New Energy Co., Ltd for his technical direction and support in Raman spectrometry. Laipeng Ma, from Institute of Metal Research (IMR), CAS, is appreciated for the graphene sample transfer to our substrates and he acknowledges the support from NSFC (no. 51102241) and MOST (no. 2012AA030303).

References

  1. K. S. Novoselov, A. K. Geim, S. V. Morozov, D. Jiang, M. I. Katsnelson, I. V. Grigorieva, S. V. Dubonos and A. A. Firsov, Nature, 2005, 438, 197 CrossRef CAS PubMed.
  2. A. K. Geim, Science, 2009, 324, 1530 CrossRef CAS PubMed.
  3. Y. Zhang, Y. W. Tan, H. L. Stormer and P. Kim, Nature, 2005, 438, 201 CrossRef CAS PubMed.
  4. K. S. Novoselov1, A. K. Geim, S. V. Morozov, D. Jiang, Y. Zhang, S. V. Dubonos, I. V. Grigorieva and A. A. Firsov, Science, 2004, 306, 666 CrossRef PubMed.
  5. J. H. Chen, C. Jang, S. Xiao, M. Ishigami and M. S. Fuhrer, Nat. Nanotechnol., 2008, 3, 206 CrossRef CAS PubMed.
  6. C. Lee, X. Wei, J. W. Kysar and J. Hone, Science, 2008, 321, 385 CrossRef CAS PubMed.
  7. G. Xie, M. Forslund and J. Pan, ACS Appl. Mater. Interfaces, 2014, 6, 7444 CAS.
  8. Y. Jiang, Y. Li, B. Liang, X. Yang, T. Han and Z. Wang, Nanotechnology, 2012, 23, 495703 CrossRef PubMed.
  9. C. Lee, Q. Li, W. Kalb, X. Z. Liu, H. Berger, R. W. Carpick and J. Hone, Science, 2010, 328, 76 CrossRef CAS PubMed.
  10. T. Filleter, J. L. McChesney, A. Bostwick, E. Rotenberg, K. V. Emtsev, T. Seyller, K. Horn and R. Bennewitz, Phys. Rev. Lett., 2009, 102, 086102 CrossRef CAS.
  11. H. Nan, Z. Ni, J. Wang, Z. Zafar, Z. Shi and Y. Wang, J. Raman Spectrosc., 2013, 44, 1018 CrossRef CAS.
  12. D. Teweldebrhan, V. Goyal and A. A. Balandin, Nano Lett., 2010, 10, 1209 CrossRef CAS PubMed.
  13. E. SuáKim and S. YoonáLee, Nanoscale, 2013, 5, 3063 RSC.
  14. J. Chen, I. Ratera, J. Y. Park and M. Salmeron, Phys. Rev. Lett., 2006, 96, 236102 CrossRef.
  15. A. C. Ferrari, J. C. Meyer, V. Scardaci, C. Casiraghi, M. Lazzeri, F. Mauri, S. Piscanec, D. Jiang, K. S. Novoselov, S. Roth and A. K. Geim, Phys. Rev. Lett., 2006, 97, 187401 CrossRef CAS.
  16. H. Xu, Y. Chen, J. Zhang and H. Zhang, Small, 2012, 8, 2833 CrossRef CAS PubMed.
  17. Z. H. Ni, H. M. Wang, Z. Q. Luo, Y. Y. Wang, T. Yu, Y. H. Wu and Z. X. Shen, J. Raman Spectrosc., 2010, 41, 479 CrossRef CAS.
  18. G. H. Han, F. Güneş, J. J. Bae, E. S. Kim, S. J. Chae, H. J. Shin, J. Y. Choi, D. Pribat and Y. H. Lee, Nano Lett., 2011, 11, 4144 CrossRef CAS PubMed.
  19. M. Ishigami, J. H. Chen, W. G. Cullen, M. S. Fuhrer and E. D. Williams, Nano Lett., 2007, 7, 1643 CrossRef CAS PubMed.
  20. G. Xie, J. Ding, B. Zheng and W. Xue, Tribol. Int., 2009, 42, 115–121 CrossRef PubMed.
  21. A. Smolyanitsky, J. P. Killgore and V. K. Tewary, Phys. Rev. B: Condens. Matter Mater. Phys., 2012, 85, 035412 CrossRef.
  22. J. S. Choi, J. S. Kim, I. S. Byun, D. H. Lee, M. J. Lee, B. H. Park, C. Lee, D. Yoon, H. Cheong, K. H. Lee, Y. W. Son, J. Y. Park and M. Salmeron, Science, 2011, 333, 607 CrossRef CAS PubMed.
  23. S. Kwon, J. H. Ko, K. J. Jeon, Y. H. Kim and J. Y. Park, Nano Lett., 2012, 12, 6043 CrossRef CAS PubMed.
  24. R. W. Carpick, D. F. Ogletree and M. Salmeron, Appl. Phys. Lett., 1997, 70, 1548 CrossRef CAS PubMed.
  25. J. Colchero, M. Luna and A. M. Baro, Appl. Phys. Lett., 1996, 68, 2896 CrossRef CAS PubMed.
  26. C. Gómez-Navarro, M. Burghard and K. Kern, Nano Lett., 2008, 8, 2045 CrossRef PubMed.
  27. J. T. Paci, T. Belytschko and G. C. Schatz, J. Phys. Chem. C, 2007, 111, 18099 CAS.
  28. Z. Deng, A. Smolyanitsky, Q. Li, X. Q. Feng and R. J. Cannara, Nat. Mater., 2012, 11, 1032 CAS.

This journal is © The Royal Society of Chemistry 2015
Click here to see how this site uses Cookies. View our privacy policy here.