Mechanical properties of commercial copper current-collector foils

Jianyu Zhua, Jiemin Fengbc and Zhansheng Guo*ab
aShanghai Institute of Applied Mathematics and Mechanics, Shanghai University, Shanghai 200072, China. E-mail: davidzsguo@shu.edu.cn; Fax: +86-21-36033287; Tel: +86-21-56331451
bShanghai Key Laboratory of Mechanics in Energy Engineering, Shanghai University, Shanghai 200072, China
cDepartment of Mechanics, College of Science, Shanghai University, Shanghai 200072, China

Received 31st July 2014 , Accepted 30th October 2014

First published on 30th October 2014


Abstract

The functionality and reliability of the current collector (CC) are crucial to design and fabricate electrodes for Li-ion batteries because the CC serves as the bridge between external electronic and internal Li-ion transports. Therefore, understanding the mechanical behavior of CCs is of great importance for battery design and manufacturing. In this paper, we report the measured values of the elastic moduli of six commercial copper current-collector (CCC) foils. Measurements were performed using three techniques: a standard microtensile testing machine equipped with a laser sensor, dynamic mechanical analysis (DMA), and nanoindentation. For electrolytic copper (E-Cu) foils, we find elastic moduli of approximately 70 GPa, and for rolled copper (R-Cu) foils, we find elastic moduli of approximately 50 GPa. Values for yield strength and fracture strength of the foils were determined from load-deflection curves; the results are consistent with values recommended by the manufacturer. Crystalline structures, which influence values for the elastic moduli of the foils, were investigated by X-ray diffraction. Surface morphologies of the foils before testing and the fracture morphologies after testing were studied by scanning electron microscopy.


1. Introduction

Lithium-ion batteries (LIBs) play an important role in storing electrical energy. However, most studies on LIB functionality have focused on performance of the cathode, anode, and electrolyte materials.1 To better understand the degradation mechanisms in LIBs, all components in the battery, including the separator2 and current collector (CC), must be tested under real-life conditions. In LIBs, copper is used as the CC for the anode and aluminum is used as the CC for the cathode.1 In recent years, many researchers have studied such degradation mechanisms as diffusion-induced stress and structural volume changes when Li ions are inserted into or extracted from active materials of electrodes.3–10 In contrast, less attention has been paid to CC foils, which play an important role in the cycle performance of LIBs. Mechanical failure of CC foils in LIBs can result in external electronic transport failures followed by battery failure.

Only a few studies have considered the influence of copper CCs (CCCs) on Li-ion diffusion and diffusion-induced stress. Zhang et al.3–5 used simulations to study the role of the CC in diffusion-induced stress of active materials of electrodes; however, in their simulation, they just used bulk properties of copper. Gao et al.11,12 showed that the CC might affect LIB performance. Nagpure et al.13 studied the influence of copper CCs on LIBs. They found that lithium impurities in the CCC lead to degradation of its thermal and electrical behavior; therefore, such impurities cannot be ignored in efforts to understand aging mechanisms that impact the life and performance of LIBs. Sa et al.14 found that the CC also plays an extremely important role in rechargeable LIBs. The physical and chemical properties of the CC can impact the performance of LIBs, and use of different CCs can result in significant differences in the performance of LIBs. Kim et al.15 demonstrated that the cycle life of a silicon–graphite composite electrode was significantly improved by using a CCC with a modified surface morphology. Cho et al.16 studied the effects of CCC substrates on electrochemical properties of Li/Si thin-film cells. They discovered that the initial capacity depends on microstructures in the copper foil and that cycle performance depends on surface roughness.

Characterizing of the mechanical properties for CCC foils, including elastic moduli, yield strengths, and fracture strengths, should enable more reasonable and reliable LIB designs and manufacturing processes. Although bulk copper is well understood, the same is not true for copper foils. This is because of the intricacies and complexities involved in fabricating the specimens and subsequent handling during testing. Since CCs affect the performance of LIBs, we report experimental results for elastic moduli of CCC foils measured by a microtensile testing machine, dynamic mechanical analysis (DMA), and nanoindentation. These tests were performed on CCC foil samples having six different thicknesses and fabricated via two different processes. Yield strengths and fracture strengths of the CCC foils were obtained from stress–strain curves obtained from the microtensile tests. Crystalline structures and surface morphologies were investigated by X-ray diffraction (XRD) and scanning electron microscopy (SEM), respectively.

2. Experiments

2.1 Materials

All Cu foils used in this investigation had purities of 99.95% and were produced industrially by rolling (R-Cu) (Lingbao Wason Copper Foil Co., Ltd., China and Dongguan Wah Wei Copper Foil Technology Ltd., China) or by electrolytic extraction (E-Cu) (Lingbao Wason Copper Foil Co., Ltd., China and Dongguan Wah Wei Copper Foil Technology Ltd., China). They were tested in their as-received conditions. Sample thickness, manufacturing process, and mechanical properties are summarized in Table 1.
Table 1 Thicknesses and mechanical properties of copper foils tested in this study
Thickness (μm) Process Surface Tensile strength (MPa) Elongation (%)
8 E-Cu Double side shiny ≥295 ≥4.0
9 R-Cu Single side shiny ≥245 ≥2.5
12 E-Cu Double side shiny ≥295 ≥5
15 E-Cu Double side shiny ≥295 ≥6
18 R-Cu Single side shiny ≥245 ≥5
35 R-Cu Single side shiny ≥245 ≥5


2.2 Microtensile tests

Tensile properties were determined using a microtensile testing machine (Zwick/roll Germany, load capacity between 10 and 2000 N) equipped with a laser sensor. To avoid the influence of detrimental vibrations, the entire testing machine was mounted on an optical table. CCC foils were glued to specially designed grips and were aligned using an xy die. A double-bladed cutter was used to prepare specimens according to the ASTM standard E345-93 (2008).17 If necessary, cutting edges of the double-bladed cutter (D-tx003, Ningbo Donrie Metal Manufacture Co., Ltd, China) should be lubricated with a material such as stearic acid in alcohol or a similar alternative. The finished specimens were examined under approximately 20× magnification to determine whether edges were smooth and that there were no surface scratches or creases. Specimens showing discernible surface scratches, creases, or edge discontinuities were rejected. Detailed descriptions of the noncontact optical method can be found in numerous other researches.18–20

2.3 Dynamic mechanical analysis tests

Elastic moduli can also be determined using a dynamic mechanical analysis machine (TA Instruments DMA Q800). Uniaxial testing clamps were used. A double-bladed cutter (D-tx003, Ningbo Donrie Metal Manufacture Co., Ltd, China) was used to cut specimens from the foils to the desired dimensions: 60 mm × 3 mm × varying thickness. A very small preload was applied to the samples to ensure that the specimens were always in tension. All samples were tested at a frequency of 1 Hz.

2.4 Nanoindentation tests

Nanoindentation techniques have been developed for probing mechanical properties of materials and have become popular experimental techniques. Our nanoindentation tests were conducted using a Hysitron TriboIndenter (Hysitron Inc., Minneapolis, MN, USA) with the continuous stiffness measurement option under constant strain-rate control. The tip on the nanoindenter system was a Berkovich diamond tip, which is preferred for low-load indentation testing. Measurements were performed by loading to the maximum load, then unloading to 95% and holding for 50 s to allow for the thermal-drift correction of the sample, and then unloading completely. The elastic modulus of each specimen was determined from the load–displacement curve using Fischer-Cripps method.21

2.5 Crystalline structure and surface morphology characterizations

Crystalline structures and crystallographic orientations of the CCC foils were determined using grazing incident-angle XRD (18 kW D/MAX2550, Rigaku, Japan). The power of the XRD (Cu Kα radiation) was fixed at 40 kV and 25 mA. The angle of the incident X-ray beam was fixed at 0.5°. To compare diffraction intensities of different sputtering samples, the XRD diffraction angle (2θ) was spanned from 30° to 100°. Surface morphologies of the CCC foils were determined by a desktop scanning electron microscope (Sh3000, Hirox, Korea) before and after tensile testing.

3. Results and discussion

3.1 Elastic modulus

Fig. 1–3 show typical stress–strain curves (load–displacement curves) for E-Cu and R-Cu CC foils obtained by microtensile, DMA, and nanoindentation tests, respectively. Copper foils of all six thicknesses had similar stress–strain curves (not shown here). For each specimen, six samples were subjected to tensile testing, three were subjected to DMA, and three indents were used for nanoindentation. Elastic moduli of the CCC foils were determined in the linear elastic region of the stress–strain curves from the microtensile test machine. The resulting values for elastic moduli are shown in Table 2. Fig. 1 shows that stress–strain curves were all reproducible, especially when the strain was below 1%. Fig. 2(a) and (b) show that elongation of the E-Cu foils was larger than those of R-Cu foils. Table 2 shows that the elastic moduli of all CCC foils were much different from that of bulk copper (for bulk Cu, a value in the range 110–126 GPa is usually as assumed in ref. 19 and 22). According to Mechant et al.,19,20 differences in texture explain variations in experimentally determined elastic moduli from 68 GPa up to 132 GPa. The texture of the copper foils in our study was studied by XRD and will be discussed below. Our results indicate that the elastic moduli for CCC foils cannot be replaced with the value for bulk Cu in calculations that predict diffusion-induced stress or fracture of electrodes in LIBs. Table 2 shows that the elastic moduli from different test methods agree reasonably well; i.e., elastic moduli obtained from stress–strain curves of microtensile tests have excellent reliability. A comparison of the elastic moduli of E-Cu (8 μm, 12 μm, and 15 μm) with those of R-Cu (9 μm, 18 μm, and 35 μm) shows that E-Cu foils have higher elastic moduli (approximately 70 GPa) than R-Cu foils (approximately 51 GPa).
image file: c4ra07675c-f1.tif
Fig. 1 Typical stress–strain curves of CCC foils by microtensile test (a) 8 μm E-Cu foils and (b) 9 μm R-Cu foils.
Table 2 Test results for elastic moduli of CCC foils
Thickness (μm) Elastic modulus (GPa)
Microtensile test DMA test Nanoindentation test
Six sample values average Three sample values average Three indents values average
8 68 71 70
9 53 51 52
12 69 70 68
15 71 67 72
18 50 49 48
35 56 52 54



image file: c4ra07675c-f2.tif
Fig. 2 Typical stress–strain curves of CCC foils by DMA test (a) 15 μm E-Cu foils and (b) 35 μm R-Cu foils.

image file: c4ra07675c-f3.tif
Fig. 3 Load–displacement curves of CCC foils by nanoindentation test (a) 12 μm E-Cu foils and (b) 35 μm R-Cu foils.

3.2 Yield strength and fracture strength

Fig. 4 shows the yield strengths for CCC foils having six different thicknesses, and Fig. 5 does the same for the fracture strength. Fig. 4 shows that decreases in foil thickness significantly increased yield strength, while Fig. 5 shows that decreases in thickness had little influence on fracture strength. Since the grain size of thin foil is smaller than the corresponding foil thickness, the foil in the thickness direction may contain few grains. Thus, the increase in the yield strength can be attributed to a grain size effect, as described by the Hall–Petch relationship. The measured fracture strengths were consistent with the values provided by the manufacturer.
image file: c4ra07675c-f4.tif
Fig. 4 Yield strength of CCC foils as a function of thickness.

image file: c4ra07675c-f5.tif
Fig. 5 Fracture strength of CCC foils as a function of thickness. The inset is the real picture of sample after microtensile test.

Fig. 6 shows the elongations of E-Cu and R-Cu CCC foils with different thicknesses. For the same fabrication process, elongation decreased remarkably with decreasing foil thickness, indicating that the plasticities of thicker foils were better than those of the thinner ones. The larger plasticities of thicker foils can be attributed to the larger grains in those foils, and to more movable dislocations stored in the grains compared to those in the thinner foils. However, the thinner foils contained relatively fewer dislocation sources, which resulted in higher breaking stresses. The measured elongations of all samples were in good agreement with those provided by the manufacturer.


image file: c4ra07675c-f6.tif
Fig. 6 Elongation of CCC foils as a function of thickness.

3.3 Crystalline structure characterization

Fig. 7 shows the XRD patterns obtained from CCC foils. Five peaks have been observed and compared with the standard powder diffraction card of JCPDS, copper file no. 04-0836. The XRD study confirms that these foils are FCC coppers. Fig. 7(a) shows the XRD results of the three different thicknesses of E-Cu CC foils. All E-Cu CC foils were found to be highly textured in the (111) and (200) planes while the intensity peak value of (220) plane was weak. Fig. 7(b) shows that all R-Cu CC foils were also highly textured in the (111) and (200) planes with a strong (220) plane. Note that the 35 μm-thick R-Cu CC foil shows a texture-like structure formed along a certain direction because it only has a strong peak in the (220) plane. According to the results of Sarada et al.,23 the foils with a strong peak in the (111) plane had high electrical conductivity, which is the key factor for use as CC in LIBs. This may be a reason why most commercial industries often select E-Cu foils and not R-Cu foils for anode CCs. Of course, R-Cu foils are usually used in devices in which the electrical conductivity requirement is a little lower than that of a battery.
image file: c4ra07675c-f7.tif
Fig. 7 XRD patterns of CCC foils(a) E-Cu foils XRD data and (b) XRD spectra for R-Cu foils XRD data.

Many microstructure factors, such as crystallographic texture, grain size, grain orientation, grain boundaries, and structural defects, have important influence on mechanical properties, such as hardness, Young's modulus, yielding strength, and fracture toughness.24–31 Here, we just provide a simple qualitative analysis of crystallographic texture on Young's modulus on the basis of the XRD results. The theoretical calculation and quantitative analysis of microstructure factors affecting mechanical properties which need more additional experiment methods are being studied. Crystallographic textures have a significant effect because the elastic modulus of copper is direction specific.26 For example, the Young's modulus (the stiffness of a material in tension) of a copper single crystal varies from 67 GPa (low value) (〈100〉 direction) to 192 GPa (high value) (〈111〉 direction) depending on orientation. In contrast polycrystalline copper has the same modulus of 111 GPa in all directions.32 The elastic modulus of a (1**) plane is higher than that of a (2**) plane. i.e., the (111) plane, which has the lowest surface energy when an fcc structure dominates the crystal,25 has the highest elastic modulus. From the XRD patterns of our CCC foils (Fig. 7), the intensity of (111) or (200) plane of E-Cu CC foils is found to be much higher than that of R-Cu foils. For example, the minimum peak value of E-copper is 106[thin space (1/6-em)]647 cps (8 μm) while the maximum peak value in (111) of R-copper is 53[thin space (1/6-em)]920 cps (18 μm); this is likely an important reason why the elastic moduli of E-copper are higher than those of R-Cu foils. Notably, the elastic modulus measured by the microtensile and DMA tests is a mixed modulus over all crystal planes. In contrast, elastic moduli determined by the nanoindentation test are more discrete than macroscopic test methods because the indenter size is at nanoscale.

3.4 Surface morphology characterization

The Surface morphologies of the CCC foils are shown in Fig. 8. The figure shows that surface morphologies were very different depending on whether the surface was shiny or matte. These differences mean that the bonding mechanisms of CCC surfaces with active materials will differ depending on the method of surface treatment. The influence of different surface treatments on electrochemical performance and interfacial bonding strength of electrodes needs further study.
image file: c4ra07675c-f8.tif
Fig. 8 Typical SEM photographs of CCC foils (a) SEM photo of shining-treated side of CCC foils and (b) SEM photo of matte-treated side of CCC foils.

Fig. 9(a) and (b) present SEM images of surface morphologies close to the fracture zones of the 15 μm-thick E-Cu and 9 μm-thick R-Cu CC foils, respectively. Fig. 9(a) is for an E-Cu CC foil and reveals typical knife-edge fracture morphology, while Fig. 9(b) is for an R-Cu CC foil and is characterized by a strong crystallographic appearance. Extensive slip bands appear in the 9 μm-thick foil (Fig. 9(b)), but only faint slip traces appear in the 15 μm-thick foil (Fig. 9(a)). This type of failure was observed for all samples, irrespective of their thicknesses. More detailed analyses of fracture mechanisms are under study.


image file: c4ra07675c-f9.tif
Fig. 9 SEM fractographs of CCC foils tensile tested (a) 15 μm E-Cu current collector foil and (b) 9 μm R-Cu current collector foil.

Conclusions

Elastic moduli of commercial CCC foils were determined using microtensile, DMA, and nanoindentation techniques. Yield strengths, fracture strengths, and elongations were determined using a standard microtensile testing machine equipped with a laser sensor. For the E-Cu foils, the measured elastic moduli were approximately 70 GPa, while those for the R-Cu foils were approximately 50 GPa. Crystalline structures and surface morphologies of the CCC foils were investigated by XRD and SEM, respectively. These results for mechanical properties (elastic modulus, yield strength, fracture strength, and elongation), microstructures, and surface morphologies should contribute to building models for studying battery performance. Only by using such measured values can accurate estimates be established for the performance of LIBs. Even though specific mechanisms for some of the studied phenomenon remain elusive, the measured values and behaviors should contribute to predictions of the cycle life of LIBs.

Acknowledgements

The authors gratefully acknowledge the financial supports by the Science and Technology Commission of Shanghai Municipality (no. 12ZR1410200), the National Science Foundation of China (no. 11472165, 11332005), Shanghai Municipal Education Commission (no. 13ZZ070) and Shanghai leading Academic Discipline Project (no. S30106).

Notes and references

  1. C. Daniel and J. O. Besenhard, Handbook of Battery Materials, Wiley-VCH Verlag & Co. KGaA, Boschstr., Weinheim, Germany, 2nd edn, 2011, vol. 12, p. 69469 Search PubMed.
  2. J. Chen, Y. Yan, T. Sun, Y. Qi and X. Li, Deformation and fracture behaviors of microporous polymer separators for lithium ion batteries, RSC Adv., 2014, 4, 14904 RSC.
  3. J. Zhang, B. Lu, Y. Song and X. Ji, Diffusion induced stress in layered Li-ion battery electrode plates, J. Power Sources, 2012, 209, 220–227 CrossRef CAS PubMed.
  4. B. Lu, Y. C. Song, Z. S. Guo and J. Q. Zhang, Analysis of delamination in thin film electrodes under galvanostatic and potentiostatic operations with Li-ion diffusion from edge, Acta Mech. Sin., 2013, 29(3), 348–356 CrossRef PubMed.
  5. Y. Song, X. Shao, Z. Guo and J. Zhang, Role of material properties and mechanical constraint on stress-assisted diffusion in plate electrodes of lithium ion batteries, J. Phys. D: Appl. Phys., 2013, 46(10), 105307 CrossRef.
  6. I. Ryu, J. W. Choi, Y. Cui and W. D. Nix, Size-dependent fracture of Si nanowire battery anodes, J. Mech. Phys. Solids, 2011, 59, 1717 CrossRef CAS PubMed.
  7. K. Zhao, M. Pharr, J. J. Vlassak and Z. Suo, Fracture of electrodes in lithium-ion batteries caused by fast charging, J. Appl. Phys., 2010, 108, 073517 CrossRef PubMed.
  8. M. Pharr, Z. Suo and J. J. Vlassak, Measurements of the Fracture Energy of Lithiated Silicon Electrodes of Li-ion Batteries, Nano Lett., 2013, 13(11), 5570–5577 CrossRef CAS PubMed.
  9. G. Bucci, S. P. V. Nadimpalli, V. A. Sethuraman, A. F. Bower and P. R. Guduru, Measurement and modeling of the mechanical and electrochemical response of amorphous Si thin film electrodes during cyclic lithiation, J. Mech. Phys. Solids, 2014, 62, 276–294 CrossRef CAS PubMed.
  10. H. Haftbaradaran, X. Xiao, M. W. Verbrugge and H. Gao, Method to deduce the critical size for interfacial delamination of patterned electrode structures and application to lithiation of thin-film silicon islands, J. Power Sources, 2012, 206, 357–366 CrossRef CAS PubMed.
  11. H. Haftbaradaran, X. Xiao and H. Gao, Critical film thickness for fracture in thin-film electrodes on substrates in the presence of interfacial sliding, Modell. Simul. Mater. Sci. Eng., 2013, 21, 074008 CrossRef.
  12. M. E. Stournara, X. Xiao, Y. Qi, P. Johari, P. Lu, B. W. Sheldon, H. Gao and V. B. Shenoy, Li segregation induces structure and strength changes at the amorphous Si/Cu interface, Nano Lett., 2013, 13(10), 4759–4768 CrossRef CAS PubMed.
  13. S. C. Nagpure, R. G. Downing, B. Bhushan and S. S. Babu, Discovery of lithium in copper current collectors used in batteries, Scr. Mater., 2012, 67, 669–672 CrossRef CAS PubMed.
  14. Q. N. Sa and Y. Wang, Ni foam as the current collector for high capacity C–Si composite electrode, J. Power Sources, 2012, 208, 46–51 CrossRef CAS PubMed.
  15. Y. L. Kim, Y. K. Sun and S. M. Lee, Enhanced electrochemical performance of silicon-based anode material by using current collector with modified surface morphology, Electrochim. Acta, 2008, 53, 4500–4504 CrossRef CAS PubMed.
  16. G. B. Cho, B. K. Lee, W. C. Sin, K. K. Cho, H. J. Ahn, T. H. Tae and K. W. Kim, Effects of Cu current collector as a substrate on electrochemical properties of Li/Si thin film cells, J. Mater. Sci., 2006, 41, 313–315 CrossRef PubMed.
  17. Standard Test Methods of Tension Testing of Metallic Foil. E345-93 (Reapproved 2008).
  18. B. Weiss, V. Gröger, G. Khatibi, A. Kotas, P. Zimprich, R. Stickler and B. Zagar, Characterization of mechanical and thermal properties of thin Cu foils and wires, Sens. Actuators, A, 2002, 99, 172–182 CrossRef CAS.
  19. H. D. Merchant, G. Khatibi and B. Weiss, Elastic and elastoplastic response of thin copper foil, J. Mater. Sci., 2004, 39, 4157–4170 CrossRef CAS.
  20. H. D. Merchant, B. Weiss and G. Khatibi, Fatigue and fracture properties of thin metallic foils, Int. J. Fract., 2001, 109, 69–89 CrossRef.
  21. A. C. Fischer-Cripps, Illustrative Analysis of Load–Displacement Curves in Nanoindentation, J. Mater. Res., 2007, 22(11), 3075–3086 CrossRef CAS.
  22. D. Ma, K. Xu, J. He and J. Lu, Evaluation of the mechanical properties of thin metal films, Surf. Coat. Technol., 1999, 116–119, 128–132 CrossRef CAS.
  23. B. V. Sarada, Ch. L. P. Pavithra, M. Ramakrishma, T. N. Rao and G. Sundararajan, Highly (111) textured copper foils with high hardness and high electrical conductivity by pulse reverse electrodeposition, Electrochem. Solid-State Lett., 2010, 13, D40–D42 CrossRef CAS PubMed.
  24. D. J. Dunstan, J. U. Galle, X. D. Hou, K. M. Y. P’ng, A. J. Bushby, B. Yang and D. Kiener, Yield and plastic flow of soft metals in small volumes loaded in tension and flexure, Philos. Mag., 2012, 92, 3199–3215 CrossRef CAS.
  25. D. A. Porter and K. E. Easterling, in Phase Transformation in Metals and Alloys, Chapman & Hall, London, 1992 Search PubMed.
  26. J. H. Choi, S. Y. Kang and D. N. Lee, Relationship between deposition and recrystallization textures of copper and chromium electrodeposits, J. Mater. Sci., 2000, 35(16), 4055–4066 CrossRef CAS.
  27. W. C. Overton Jr and J. Gafney, Temperature Variation of the Elastic Constants of Cubic Elements. I. Copper, Phys. Rev., 1955, 98, 969–977 CrossRef.
  28. B. Merle, E. W. Schweitzer and M. Göken, Thickness and grain size dependence of the strength of copper thin films as investigated with bulge tests and nanoindentations, Philos. Mag., 2012, 92(25–27), 3172–3187 CrossRef CAS.
  29. J. D. Giallonardo, U. Erb, K. T. Aust and G. Palumbo, The influence of grain size and texture on the Young's modulus of nanocrystalline nickel and nickel–iron alloys, Philos. Mag., 2011, 91(36), 4594–4605 CrossRef CAS.
  30. S. J. Skrzypek, W. ratuszek and A. Bunsch, et al., Crystallographic texture and anisotropy of electrolytic deposited copper coating analysis, Journal of Achievements of Materials and Manufacturing Engineering, 2010, 43, 264–268 Search PubMed.
  31. D. E. J. Armstrong, A. J. Wilkinson and S. G. Roberts, Measuring anisotropy in Young’s modulus of copper using microcantilever testing, J. Mater. Res., 2009, 24, 3268–3276 CrossRef CAS.
  32. http://www.aem.umn.edu/people/faculty/shield/c1b1/pcb4.htm.

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