Stark sublevels of Er3+–Yb3+ codoped Gd2(WO4)3 phosphor for enhancing the sensitivity of a luminescent thermometer

Hongyu Lua, Ran Menga, Haoyue Haoa, Yunfeng Baia, Yachen Gaob, Yinglin Songa, Yuxiao Wang*a and Xueru Zhang*a
aDepartment of Physics, Harbin Institute of Technology, Harbin 150001, China. E-mail: wangyx@hit.edu.cn; xrzhang@hit.edu.cn
bCollege of Electronic Engineering, Heilongjiang University, Harbin 150080, China

Received 19th April 2016 , Accepted 6th June 2016

First published on 7th June 2016


Abstract

We report a strategy for enhancing the sensitivity of an optical thermometer by utilizing the Stark sublevels. Under 980 nm excitation, the upconversion emission originating from the Stark levels of Er3+ (4S3/2(2)/4S3/2(1) and 4F9/2/4F9/2(1)) is observed in a monoclinic phase Gd2(WO4)3:Er3+/Yb3+ phosphor, which is synthesized through the co-precipitation method. The temperature sensing behavior is studied over the range of 296–620 K based on thermal coupling levels 2H11/2/4S3/2(2), 2H11/2/4S3/2, 2H11/2/4S3/2(1), 4S3/2(2)/4S3/2(1) and 4F9/2(2)/4F9/2(1) using fluorescence intensity ratio (FIR) technology. The sensitivity based on Stark sublevels 2H11/2/4S3/2(2) or 2H11/2/4S3/2(1) is almost twice more than that based on the traditional 2H11/2/4S3/2 levels. The obtained maximum sensitivity of the optical thermometer is 16.5 × 10−3 K−1 at 395 K (2H11/2/4S3/2(2)). These results suggest that the use of Stark levels is a promising approach for enhancing the sensitivity of optical thermometers.


Introduction

Highly sensitive and highly accurate temperature sensing is a challenging task in both scientific research and industrial production. Recently, rare earth (RE) doped materials have attracted great attention as temperature probes, especially for their potential application in intracellular environments and strong electromagnetic fields (without interferences).1–8 As we know, the FIR technique is an effective approach to achieve accurate measurements by exploiting the temperature-dependent luminescence intensities derived from two thermal coupling levels, which are unaffected by spectrum losses and fluctuations in excitation intensity.9–11 Photon upconversion (UC) is basically a nonlinear optical effect, which can convert low energy photons to high energy photons via multiphoton processes.12,13 As one of the most widely used activators, erbium ions (Er3+) have abundant ladder-like electronic energy levels.14–17 Under near-infrared (NIR) excitation, green UC emissions originating from 2H11/2 and 4S3/2 to ground state transitions are highly efficient and the energy gap between them matches well with the requirement for the thermal coupling levels (200 cm−1 < ΔE < 2000 cm−1).18–20

According to previous literature, the Stark sublevels are considered as thermal coupling levels due to the small energy gap.21,22 Jaque et al. investigated the subtissue thermal sensing behavior based on the Stark sublevels (4F5(2) and 4F5(1)) of Nd3+ ions.23 Dong et al. reported highly penetrating bio-imaging and temperature sensing using the Stark sublevels (3H4(2) and 3H4(1)) of Tm3+ ions.24 Yin et al. achieved high sensitivity sensing in the low temperature range using the Stark sublevels (1G4(2) and 1G4(2)) of Tm3+ ions.25 Although the Stark energy levels of RE ions, such as Tm3+ (1G4(2) and 1G4(1)), Ho3+ (5F5(2) and 5F5(1)), Nd3+ (4F3/2(2) and 4F3/2(1)), and so on have been widely used in optical thermometers,26–29 few reports focused on the Stark energy levels of Er3+ ions. Generally, the larger the energy gap between the thermal coupling levels is, the higher the sensor sensitivity. Because the sensitivity is restricted to the energy gap of the thermal coupling levels, utilizing thermal coupling levels with a large energy gap may be a valid strategy to improve sensitivity. In addition, due to its remarkable thermal stability, distinct crystal structure and low phonon energy, Gd2(WO4)3 is regarded as a promising host material for temperature sensing.30–32

In the present work, we prepared the monoclinic phase Gd2(WO4)3:10%Er3+/10%Yb3+ phosphor via the co-precipitation method. Emission bands arising due to the transitions from Stark sublevels were observed. Using FIR technology, the temperature sensing behavior was studied over the range of 296–620 K based on five pairs of thermal coupling levels.

Experimental

The Gd2(WO4)3:10%Er3+/10%Yb3+ phosphor was synthesized via a co-precipitation method. In a typical procedure, 10 mL of solution containing Gd(NO3)3 (1.8 mmol), Er(NO3)3 (0.1 mmol) and Yb(NO3)3 (0.1 mmol) was slowly added into 30 mL of Na2WO4 (3 mmol) solution, and then the mixed solution was magnetically stirred for 30 min. The product was collected by centrifugation (8000 rpm, 10 min) and washed twice with deionized water. The precipitate was dried at 70 °C for 2 h. Finally, the Gd2(WO4)3:Yb3+/Er3+ phosphor was obtained after being calcined at 900 °C for 4 h. The pH of the Na2WO4 solution was adjusted to be lower than 10 using ammonia. Many repeated experiments indicated that a 10%Er3+ and 10%Yb3+ codoped Gd2(WO4)3 phosphor is an optimal result for relatively strong red emission, which is more suitable for measurement and application.

Upconversion luminescence spectra were measured using HORIBA Jobin Yvon iHR550 spectrometers which employ a photon-counting detection system using a Hamamatsu R928 photomultiplier tube. Powder X-ray diffraction (XRD) patterns were recorded using an X-ray diffractometer (Panalytical Empyrean).

Results and discussion

The XRD pattern of the Gd2(WO4)3:10%Er3+/10%Yb3+ phosphor is shown in Fig. 1. All the diffraction peaks of the present material match well with the monoclinic phase Gd2(WO4)3 (JCPDS 23-1076) and there is an absence of impurity peaks, indicating that high purity monoclinic phase Gd2(WO4)3 nanoparticles have been successfully synthesized. Under 980 nm excitation, the UC spectra of Gd2(WO4)3:Er3+/Yb3+ at various temperatures are shown in Fig. 2. The intensity of each emission band shows an obvious decrease with increasing temperature. One red and two green emission bands of the Er3+ ions are observed. Fig. 3 gives the schematic for the UC luminescence mechanism for the Gd2(WO4)3:Er3+/Yb3+ phosphor under 980 nm excitation. It is interesting to find that both the longer wavelength green emission band and the red emission band are strongly divided into two parts, which are attributed to the 4S3/2 and 4F9/2 levels being split into the 4S3/2(2)/4S3/2(1) and 4F9/2(2)/4F9/2(1) levels, respectively. For a typical upconversion mechanism of an Er3+ ion, the population at the 2H11/2 and 4F9/2 (4F9/2(2)/4F9/2(1)) levels can be realized through two-photon processes.33,34 The 2H11/2 and 4S3/2 (4S3/2(2)/4S3/2(1)) levels are mainly populated via multiphonon relaxation of the 4F7/2 level.35 The luminescence bands centered at 525, 544, 551, 657 and 666 nm are derived from the 2H11/24I15/2, 4S3/2(2)4I15/2, 4S3/2(1)4I15/2, 4F9/2(2)4I15/2 and 4F9/2(1)4I15/2 transitions, respectively. Owing to the small energy gap, these sublevels can be thermally populated from each other indicating that these two levels are thermally coupled. In addition, it is well known that the 2H11/2 and 4S3/2 levels are thermal coupling levels. Therefore, any pair of the 2H11/2, 4S3/2(2) and 4S3/2(1) levels is a thermal coupling level. According to the emission spectra (Fig. 2), the energy gaps of the thermal coupling levels are estimated as about 485, 600, 719, 234 and 205 cm−1 for the 2H11/2/4S3/2(2), 2H11/2/4S3/2, 2H11/2/4S3/2(1), 4S3/2(2)/4S3/2(1) and 4F9/2(2)/4F9/2(1) levels, respectively.
image file: c6ra10138k-f1.tif
Fig. 1 XRD pattern of Gd2(WO4)3:10%Er3+/10%Yb3+ and the standard data for monoclinic phase Gd2(WO4)3 (JCPDS 23-1076).

image file: c6ra10138k-f2.tif
Fig. 2 Temperature-dependent upconversion emission spectra of the Gd2(WO4)3:Er3+/Yb3+ phosphor with 980 nm excitation. The inset shows the emission spectra in the wavelength range 600–725 nm.

image file: c6ra10138k-f3.tif
Fig. 3 Schematic energy level diagram for the Er3+ and Yb3+ system and the proposed UC emission mechanism under 980 nm excitation.

Since the relative population of thermal coupling levels follows the Boltzmann distribution, the FIR of the thermal coupling levels can be represented as:36–38

 
image file: c6ra10138k-t1.tif(1)
where I1 (lower level) and I2 (upper level) are the integrated intensities of the two thermal coupling levels. Ai, gi and ωi are the radiative transition probability, the degeneracy and the angular frequency of corresponding transitions, respectively. T, k and C (C = g2A2ω2/g1A1ω1) are the absolute temperature, the Boltzmann constant and the proportionality factor, respectively. ΔE is the energy gap of the thermal coupling levels.

According to FIR technique, appropriate thermal coupling levels are a vital and decisive factor in the temperature sensing behavior. Thus, the temperature sensing behavior based on three pairs of thermal coupling levels (2H11/2/4S3/2(2), 2H11/2/4S3/2 and 2H11/2/4S3/2(1)) are investigated. The FIRs as a function of temperature are plotted in Fig. 4 over the temperature range 296–620 K. The corresponding data can be fitted with exponential functions, which can be expressed as: I525/I544 = 24.2[thin space (1/6-em)]exp(−790/T), I525/(I544 + I551) = 13.36[thin space (1/6-em)]exp(−890/T), and I525/I551 = 28.5[thin space (1/6-em)]exp(−968/T).


image file: c6ra10138k-f4.tif
Fig. 4 (a) The FIRs based on thermal coupling levels 2H11/2/4S3/2(2), 2H11/2/4S3/2(1) and 2H11/2/4S3/2 as a function of temperature. (b) The sensitivity based on thermal coupling levels 2H11/2/4S3/2(2), 2H11/2/4S3/2(1) and 2H11/2/4S3/2 as a function of temperature.

For sensing applications, the absolute sensitivity (S) is a very important parameter, which has been calculated using the formula:39–41

 
image file: c6ra10138k-t2.tif(2)

The corresponding sensitivity curves are calculated and plotted in Fig. 4(b). Compared with ordinary thermal coupling levels (2H11/2/4S3/2), the sensitivities of optical thermometers based on both 2H11/2/4S3/2(2) and 2H11/2/4S3/2(1) levels increase dramatically. Obviously, the sensitivity based on the 2H11/2 and 4S3/2(2) levels is the highest over the temperature range of 296–485 K. Upon further increasing the temperature, the sensitivity based on the 2H11/2 and 4S3/2(1) levels is dominant over the temperature range of 485–620 K. Based on the three pairs of thermal coupling levels, the sensitivities first increase, and then after a certain temperature they gradually decrease as the temperature increases. The corresponding maximal sensitivity values are 16.5 × 10−3 K−1 at 395 K (2H11/2/4S3/2(2)), 15.9 × 10−3 K−1 at 485 K (2H11/2/4S3/2(1)) and 8.1 × 10−3 K−1 at 446 K (2H11/2/4S3/2). These significantly enhanced sensitivities are attributed to the energy gaps of the thermal coupling levels and the corresponding radiative transition probabilities (Ai). The above demonstrations reveal that the high sensitivity of an optical thermometer based on 2H11/2/4S3/2(2) (296–485 K) and 2H11/2/4S3/2(1) (485–620 K) was perfectly realized in the Gd2(WO4)3:Er3+/Yb3+ phosphor, whose sensitivities are about twice those based on the traditional 2H11/2 and 4S3/2 levels. Table 1 lists the sensitivities of typical Er3+/Yb3+ doped optical thermometers.

Table 1 Values of the maximum sensitivity of typical Er3+/Yb3+ doped matrices, including the temperature of the maximum relative thermal sensitivity (Smax) and the temperature range
Material Smax (× 10−3 K−1) (temperature K) Temperature range (K) References
Er, Yb:Ba5Gd8Zn4O21 3.2 (490 K) 260–490 K 42
Er, Yb:NaYF4 3.68 (363 K) 223–403 K 43
Er, Yb:YVO4 11.69 (380 K) 300–485 K 44
Er, Yb:SrWO4 14.98 (403 K) 299–518 K 45
Er, Yb:NaGd(WO4)2 11.49 (453 K) 293–573 K 46
Er, Yb:Gd2O3 8.4 (570 K) 298–723 K 47
Er, Yb:YAG 1.7 (404 K) 298–573 K 48
Er, Yb:NaGdTiO4 4.5 (480 K) 300–510 K 49
Er, Yb:Gd2(WO4)3 16.5 (395 K) 296–620 K This work


In addition, the temperature sensing behavior based on the Stark sublevels of 4S3/2 and 4F9/2 (4S3/2(2)/4S3/2(1) and 4F9/2(2)/4F9/2(1)) is investigated over the range 296–620 K. Fig. 5 show the FIRs (a) and the sensitivities (b) based on the Stark sublevels as a function of temperature. It is worth noting that the sensitivities based on the Stark levels exhibit a similar trend whereby they keep decreasing as the temperature increases over the whole temperature range. The sensitivity based on the 4S3/2(2) and 4S3/2(1) levels is higher than that based on 4F9/2(2) and 4F9/2(1), which is mainly attributed to the difference in the energy gap. The maximum sensitivities are observed to be 1.31 × 10−3 K−1 and 1.27 × 10−3 K−1 at 296 K based on the Stark sublevels of the Er3+ ions, 4S3/2(2)/4S3/2(1) and 4F9/2(2)/4F9/2(1), respectively.


image file: c6ra10138k-f5.tif
Fig. 5 (a) The FIRs based on the Stark levels 4S3/2(2)/4S3/2(1) and 4F9/2(2)/4F9/2(1) as a function of temperature. (b) The sensitivities based on the Stark levels 4S3/2(2)/4S3/2(1) and 4F9/2(2)/4F9/2(1) as a function of temperature.

Conclusion

In summary, the monoclinic phase 10%Er3+/10%Yb3+ codoped Gd2(WO4)3 phosphor is synthesized successfully via the co-precipitation method and is characterized using XRD. Under 980 nm excitation, visible luminescence arising from Stark sublevels to ground state transitions (4S3/2(2)4I15/2 (544 nm), 4S3/2(1)4I15/2 (551 nm), 4F9/2(2)4I15/2 (657 nm) and 4F9/2(1)4I15/2 (666 nm)) is observed. Using the FIR technique, the temperature sensing behavior is investigated based on thermal coupling levels 2H11/2/4S3/2(2), 2H11/2/4S3/2, 2H11/2/4S3/2(1), 4S3/2(2)/4S3/2(1) and 4F9/2(2)/4F9/2(1). The high sensitivity of the optical thermometer based on the 2H11/2/4S3/2(2) (296–485 K) and 2H11/2/4S3/2(1) (485–620 K) levels are perfectly realized, whose sensitivities are about twice more than those based on the traditional 2H11/2 and 4S3/2 levels. The obtained maximum sensitivity of the optical thermometer is 16.5 × 10−3 K−1 at 395 K (2H11/2/4S3/2(2)). These merits indicate that the Gd2(WO4)3:Er3+/Yb3+ phosphor is robust in high sensitivity temperature sensing.

Acknowledgements

This work has been supported by the grant of National Natural Science Foundation of China under project No. 11374079, 11474078 and 61275117.

References

  1. J. Zhou, Q. Liu, W. Feng, Y. Sun and F. Y. Li, Chem. Rev., 2015, 115, 395–465 CrossRef CAS PubMed.
  2. X. D. Wang, O. S. Wolfbeis and R. J. Meier, Chem. Soc. Rev., 2013, 42, 7834–7869 RSC.
  3. D. Yang, Z. Hou, Z. Cheng, C. Li and J. Lin, Chem. Soc. Rev., 2015, 44, 1416–1448 RSC.
  4. H. Dong, L.-D. Sun and C.-H. Yan, Nanoscale, 2013, 5, 5703–5714 RSC.
  5. D. Jaque and F. Vetrone, Nanoscale, 2012, 4, 4301–4326 RSC.
  6. F. Vetrone, R. Naccache, A. Zamarron, A. J. de la Fuente, F. Sanz-Rodriguez, L. M. Maestro, E. M. Rodriguez, D. Jaque, J. G. Sole and J. A. Capobianco, ACS Nano, 2010, 4, 3254–3258 CrossRef CAS PubMed.
  7. G. Chen, H. Qju, P. N. Prasad and X. Chen, Chem. Rev., 2014, 114, 5161–5214 CrossRef CAS PubMed.
  8. Y. Sun, W. Feng, P. Yang, C. Huang and F. Li, Chem. Soc. Rev., 2015, 44, 1509–1525 RSC.
  9. K. Zheng, Z. Liu, C. Lv and W. Qin, J. Mater. Chem. C, 2013, 1, 5502–5507 RSC.
  10. S. Zheng, W. Chen, D. Tan, J. Zhou, Q. Guo, W. Jiang, C. Xu, X. Liu and J. Qiu, Nanoscale, 2014, 6, 5675–5679 RSC.
  11. X. Xu, Z. Wang, P. Lei, Y. Yu, S. Yao, S. Song, X. Liu, Y. Su, L. Dong and J. Feng, ACS Appl. Mater. Interfaces, 2015, 7, 20813–20819 CAS.
  12. W. Zheng, P. Huang, D. Tu, E. Ma, H. Zhu and X. Chen, Chem. Soc. Rev., 2015, 44, 1379–1415 RSC.
  13. X. Liu, R. Deng, Y. Zhang, Y. Wang, H. Chang, L. Huang and X. Liu, Chem. Soc. Rev., 2015, 44, 1479–1508 RSC.
  14. D. Chen, Z. Wan, Y. Zhou, X. Zhou, Y. Yu, J. Zhong, M. Ding and Z. Ji, ACS Appl. Mater. Interfaces, 2015, 7, 19484–19493 CAS.
  15. B. Chen, B. Dong, J. Wang, S. Zhang, L. Xu, W. Yu and H. Song, Nanoscale, 2013, 5, 8541–8549 RSC.
  16. K. Zheng, W. Song, G. He, Z. Yuan and W. Qin, Opt. Express, 2015, 23, 7653–7658 CrossRef PubMed.
  17. E. C. Ximendes, U. Rocha, C. Jacinto, K. U. Kumar, D. Bravo, F. J. López, E. M. Rodríguez, J. García-Solé and D. Jaque, Nanoscale, 2016, 8, 3057–3066 RSC.
  18. R. Dey and V. K. Rai, Dalton Trans., 2014, 43, 111–118 RSC.
  19. B. Singh, A. Parchur, R. Ningthoujam, P. Ramakrishna, S. Singh, P. Singh, S. Rai and R. Maalej, Phys. Chem. Chem. Phys., 2014, 16, 22665–22676 RSC.
  20. L. Li, C. F. Guo, S. Jiang, D. K. Agrawal and T. Li, RSC Adv., 2014, 4, 6391–6396 RSC.
  21. O. A. Savchuk, J. J. Carvajal, M. C. Pujol, E. W. Barrera, J. Massons, M. Aguilo and F. Diaz, J. Phys. Chem. C, 2015, 119, 18546–18558 CAS.
  22. A. K. Soni, R. Dey and V. K. Rai, RSC Adv., 2015, 5, 34999–35009 RSC.
  23. U. s. Rocha, C. Jacinto da Silva, W. Ferreira Silva, I. Guedes, A. Benayas, L. Martinez Maestro, M. Acosta Elias, E. Bovero, F. C. van Veggel, J. A. Garcia Sole and D. Jaque, ACS Nano, 2013, 7, 1188–1199 CrossRef CAS PubMed.
  24. N.-N. Dong, M. Pedroni, F. Piccinelli, G. Conti, A. Sbarbati, J. E. Ramírez-Hernández, L. M. Maestro, M. C. Iglesias-de la Cruz, F. Sanz-Rodriguez and A. Juarranz, ACS Nano, 2011, 5, 8665–8671 CrossRef CAS PubMed.
  25. H. Suo, C. Guo, Z. Yang, S. Zhou, C. Duan and M. Yin, J. Mater. Chem. C, 2015, 3, 7379–7385 RSC.
  26. L. Li, H. Lin, X. Zhao, Y. Wang, X. Zhou, C. Ma and X. Wei, J. Alloys Compd., 2014, 586, 555–560 CrossRef CAS.
  27. P. Haro-González, L. M. Maestro, M. Trevisani, S. Polizzi, D. Jaque, J. G. Sole and M. Bettinelli, J. Appl. Phys., 2012, 112, 054702 CrossRef.
  28. A. Dwivedi, K. Mishra and S. Rai, J. Phys. D: Appl. Phys., 2015, 48, 435103 CrossRef.
  29. L. Li, L. Zheng, W. Xu, Z. Liang, Y. Zhou, Z. Zhang and W. Cao, Opt. Lett., 2016, 41, 1458–1461 CrossRef PubMed.
  30. M. Yin, Y. Liu, L. Mei, M. S. Molokeev, Z. Huang and M. Fang, RSC Adv., 2015, 5, 73077–73082 RSC.
  31. Y. Zeng, Z. Li, L. Wang and Y. Xiong, CrystEngComm, 2012, 14, 7043–7048 RSC.
  32. J. Li, J. Sun, J. Liu, X. Li, J. Zhang, Y. Tian, S. Fu, L. Cheng, H. Zhong and H. Xia, Mater. Res. Bull., 2013, 48, 2159–2165 CrossRef CAS.
  33. N. Niu, P. Yang, F. He, X. Zhang, S. Gai, C. Li and J. Lin, J. Mater. Chem., 2012, 22, 10889–10899 RSC.
  34. N. Niu, F. He, S. Gai, C. Li, X. Zhang, S. Huang and P. Yang, J. Mater. Chem., 2012, 22, 21613–21623 RSC.
  35. L. Marciniak, K. Prorok, L. Francés-Soriano, J. Pérez-Prieto and A. Bednarkiewicz, Nanoscale, 2016, 8, 5037–5042 RSC.
  36. X. Wang, Q. Liu, Y. Bu, C.-S. Liu, T. Liu and X. Yan, RSC Adv., 2015, 5, 86219–86236 RSC.
  37. G. Jiang, S. Zhou, X. Wei, Y. Chen, C. Duan, M. Yin, B. Yang and W. Cao, RSC Adv., 2016, 6, 11795–11801 RSC.
  38. N. Rakov and G. S. Maciel, Opt. Lett., 2014, 39, 3767–3769 CrossRef CAS PubMed.
  39. W. Tang, S. Wang, Z. Li, Y. Sun, L. Zheng, R. Zhang, B. Yang, W. Cao and M. Yu, Appl. Phys. Lett., 2016, 108, 061902 CrossRef.
  40. S. León-Luis, V. Monteseguro, U. Rodríguez-Mendoza, M. Rathaiah, V. Venkatramu, A. Lozano-Gorrín, R. Valiente, A. Muñoz and V. Lavín, RSC Adv., 2014, 4, 57691–57701 RSC.
  41. W. Yu, W. Xu, H. Song and S. Zhang, Dalton Trans., 2014, 43, 6139–6147 RSC.
  42. H. Suo, C. Guo and T. Li, J. Phys. Chem. C, 2016, 120, 2914–2924 CAS.
  43. B. Dong, R. N. Hua, B. S. Cao, Z. P. Li, Y. Y. He, Z. Y. Zhang and O. S. Wolfbeis, Phys. Chem. Chem. Phys., 2014, 16, 20009–20012 RSC.
  44. M. K. Mahata, K. Kumar and V. K. Rai, Sens. Actuators, B, 2015, 209, 775–780 CrossRef CAS.
  45. A. Pandey, V. K. Rai, V. Kumar, V. Kumar and H. Swart, Sens. Actuators, B, 2015, 209, 352–358 CrossRef CAS.
  46. J. Liao, L. Nie, Q. Wang, S. Liu, H.-R. Wen and J. Wu, RSC Adv., 2016, 6, 35152–35159 RSC.
  47. Y. Tian, B. Tian, P. Huang, L. Wang and B. Chen, RSC Adv., 2015, 5, 14123–14128 RSC.
  48. G. Liu, L. Fu, Z. Gao, X. Yang, Z. Fu, Z. Wang and Y. Yang, RSC Adv., 2015, 5, 51820–51827 RSC.
  49. D. He, C. Guo, S. Jiang, N. Zhang, C. Duan, M. Yin and T. Li, RSC Adv., 2015, 5, 1385–1390 RSC.

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