Научная статья на тему 'Size-dependent concentration of N0 paramagnetic centres in HPHT nanodiamonds'

Size-dependent concentration of N0 paramagnetic centres in HPHT nanodiamonds Текст научной статьи по специальности «Физика»

i Надоели баннеры? Вы всегда можете отключить рекламу.
Область наук
Ключевые слова

Аннотация научной статьи по физике, автор научной работы — Yavkin B.V., Mamin G.V., Gafurov M.R., Orlinskii S.B.

Size-calibrated commercial nanodiamonds synthesized by high-pressure high-temperature (HPHT) technique were studied by high-frequency Wand conventional X-band electron paramagnetic resonance (EPR) spectroscopy. The numbers of spins in the studied samples were estimated. The coreshell model of the HPHT nanodiamonds was proposed to explain the observed dependence of the concentration of the N 0 paramagnetic centers. Two other observed paramagnetic centers are attributed to the two types of structures in the nanodiamond shell.

i Надоели баннеры? Вы всегда можете отключить рекламу.
iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.
i Надоели баннеры? Вы всегда можете отключить рекламу.

Текст научной работы на тему «Size-dependent concentration of N0 paramagnetic centres in HPHT nanodiamonds»

ISSN 2072-5981

Volume 17, Issue 1 Paper No 15101, 1-7 pages 2015

http://mrsej.kpfu.ru http://mrsej.ksu.ru

Established and published by Kazan University Sponsored by International Society of Magnetic

Resonance (ISMAR) Registered by Russian Federation Committee on Press,

August 2, 1996 First Issue was appeared at July 25, 1997

© Kazan Federal University (KFU)*

"Magnetic Resonance in Solids. Electronic Journal" (MRSey) is a

peer-reviewed, all electronic journal, publishing articles which meet the highest standards of scientific quality in the field of basic research of a magnetic resonance in solids and related phenomena. MRSey is free for the authors (no page charges) as well as for the readers (no subscription fee). The language of MRSey is English. All exchanges of information will take place via Internet. Articles are submitted in electronic form and the refereeing process uses electronic mail. All accepted articles are immediately published by being made publicly available by Internet (http://MRSez.kpfu.ru).

Editors-in-Chief Jean Jeener (Universite Libre de Bruxelles, Brussels) Boris Kochelaev (KFU, Kazan) Raymond Orbach (University of California, Riverside)

Executive Editor Yurii Proshin (KFU, Kazan) mrsej@kpfu.ru



Vadim Atsarkin (Institute of Radio Engineering and Electronics, Moscow) Yurij Bunkov (CNRS, Grenoble) Mikhail Eremin (KFU, Kazan) David Fushman (University of Maryland,

College Park)

Hugo Keller (University of Zürich, Zürich) Yoshio Kitaoka (Osaka University, Osaka) Boris Malkin (KFU, Kazan) Alexander Shengelaya (Tbilisi State University, Tbilisi) Jörg Sichelschmidt (Max Planck Institute for Chemical Physics of Solids, Dresden) Haruhiko Suzuki (Kanazawa University,

Kanazava) Murat Tagirov (KFU, Kazan) Dmitrii Tayurskii (KFU, Kazan) Valentin Zhikharev (KNRTU, Kazan)

In Kazan University the Electron Paramagnetic Resonance (EPR) was discovered by Zavoisky E.K. in 1944.


Size-dependent concentration of N0 paramagnetic centres in HPHT nanodiamonds

B.V. Yavkin *, G.V. Mamin, MR. Gafurov, S.B. Orlinskii Kazan Federal University, Kremlevskaya 18, 420008 Kazan, Russia * E-mail: boris.yavkin@gmail.com (Received November 23, 2015; revised December 8, 2015; accepted December 11, 2015)

Size-calibrated commercial nanodiamonds synthesized by high-pressure high-temperature (HPHT) technique were studied by high-frequency W- and conventional X-band electron paramagnetic resonance (EPR) spectroscopy. The numbers of spins in the studied samples were estimated. The core-shell model of the HPHT nanodiamonds was proposed to explain the observed dependence of the concentration of the N0 paramagnetic centers. Two other observed paramagnetic centers are attributed to the two types of structures in the nanodiamond shell.

PACS: 75.20.-g, 76.30.Mi, 81.05.ug

Keywords: diamond, nanodaimond, EPR, ESR, crystal defect 1. Introduction

Nanodiamonds (NDs) have attracted increasing attention worldwide after a number of breakthroughs in synthesis, purification, isolation and their surface modification techniques achieved in the late 1990s [1]. At present, NDs can be synthesized at a relatively low cost by various synthesis techniques: detonation, laser ablation, high-energy ball milling of high-pressure high-temperature (HPHT) diamond microcrystals, plasma-assisted chemical vapour deposition (CVD), autoclave synthesis from supercritical fluids, chlorination of carbides, ion irradiation of graphite, electron irradiation of carbon onions [2]. Various applications of nanodiamonds have been suggested in recent years, particularly in the biological and medical fields, due to their excellent biocompatibility and superior optical properties [3-5]. The importance of nanodiamond surface study is also pumping by few vibrant research fields, like in situ magnetometry, thermometry and electric field measurements, quantum computing [6-8].

In as-produced detonation nanodiamonds the primary diamond particles are tightly bounded to each other forming large aggregates which are difficult to separate using conventional treatments [9]. New disaggregation approaches have been developed to obtain mono-dispersed nanodiamonds required by most biomedical applications. To date single-digital nanodiamonds with the size of about 4-5 nm in diameter can be achieved routinely and even smaller diamond particles are of interest [10]. The cause of strong aggregation of nanodiamond and the mechanism for the interaction between nanodiamonds and adsorbant are still challenges [11].

Samples produced using different synthesis techniques often have distinct sizes and surfaces which lead to significant differences in chemical reactivity and affinity to specific adsorbates. Purification treatments always introduce additional chemical terminations to the nanodiamond surface. Therefore, the characterisation of the ND nanoparticles and ND surfaces reported in the literature is often inconsistent [12].

Numerous and solid contributions to the study of nanodiamonds of various origin, size and surface termination were made by microscopic techniques such as transmission electron (TEM) [13], atomic force (AFM) [14] and confocal microscopy [15], NMR techniques [16], Raman spectroscopy [17], optically detected magnetic resonance (ODMR) [18, 19] technique. Unfortunately, the EPR experiments over the long period of time were on the backstage of nanodiamond research due to its limited sensitivity and lack of direct and unambiguous assignment of the observed signals [13]. In addition, despite large nitrogen content in nanodiamond, paramagnetic centers of nitrogen and

nitrogen-vacancy centers that one could routinely observe in bulk diamond crystals [20-22] have not been observed in nanodiamond materials until recently [23-28].

In this work we investigate the evolutions of X-(9.6 GHz) and W-(94 GHz) bands EPR spectra in a series of the size-calibrated (10-1000 nm) commercial nanodiamonds prepared by the HPHT technique by exploiting conventional continuous wave (CW) methods of detection for qualitative and quantitative characterization of the core-shell structures of nanodiamonds.

2. Materials and Methods

A series of the size-calibrated nanodiamonds obtained from Tomei Diamond, Japan was studied: seven samples with the average size from 1 |im down to 10 nm, produced by HPHT synthesis from the diamond crystals of type Ib with the high nitrogen content in form of single substitutional nitrogen. Procedure of small nanodiamond particle production could be found elsewhere [9]. The sizes (d) of the investigated nanodiamonds and their notations throughout the article are presented in Table 1.

The EPR experimental data were acquired at room temperature by using the capabilities of Bruker ElexSys 580/680 spectrometer. CW EPR spectra in the W-band microwave frequency were used to identify the contributions from several paramagnetic centers to the overall EPR spectrum. To attain quantitative information about the paramagnetic centers of different nature, conventional X-band spectroscopy was employed and MgO:Mn crystal powder with the concentration of paramagnetic species of 5-1014 spins per sample was used as a reference. To avoid saturation effects, the microwave power was chosen to be 0.2 |W for both X- and W-band measurements. The modeling of experimental EPR spectra was realized using EasySpin toolbox for Matlab [29].

3. Results

The W-band EPR spectra and their simulations under the assumption of the presence of three paramagnetic centers with S = 1/2 (denoted as N0, SC1 and SC2) are shown in Figure 1. The g-factors of two of them (SC1 and SC2) were calculated with respect to the reference of g = 2.0024 for N0 center. The changes of g-factors with the MD size-were not monitored. The values of g-factors and linewidths for all of the samples and obtained paramagnetic centers are listed in Table 2.

The N0 paramagnetic center shows up as a five lines powder EPR spectrum due to the anisotropic hyperfine interaction with 14N of 99.6 % natural abundance and nuclear spin I = 1, described by the spin Hamiltonian H = gfi BS+ ASZIZ + A± (SXIX + SYIY ) with g-factor of 2.0024 and hyperfine

interaction constants A± = 2.92 mT; A\ = 4 08 mT. This description agrees both qualitatively and

quantitatively with the previously published data on the nitrogen donors in the monocrystalline and powder diamond samples [20, 30]. Two other centers named SC1 and SC2 have g-factors of 2.0028(2) and 2.0029(2), respectively, and described by the lorentz and gauss lineshapes. The g-factors of SC1 and SC2 centers were determined relatively to the g-factor of N0 center that was kept constant in our fittings for the whole size series; hyperfine interaction parameters and linewidth of nitrogen center were kept constant as well. It is interesting to note that the lineshape and width of SC1 component are practically coincide with those reported in one of the first study of deformation-induced paramagnetic centers in diamond [31, 32]. Separation of SC1 and SC2 centers became possible only by using high-frequency (high-field) EPR. One might try to describe both SC1 and SC2 as a single

Table 1. Size distribution of MD nanodiamond samples

Sample Median size d, nm

MD10 5-15

MD20 15-25

MD50 40-64

MD100 90-129

MD200 180-229

MD500 480-529

MD1000 980-1029



MD1000 MD500 MD200 MD100

MD200 experiment
















magnetic field, mT

magnetic field, mT

Figure 1.

(Left panel) CW EPR spectra of nanodiamond samples (solid lines) and their corresponding simulations (red dashed lines). (Right panel) example of the spectral decomposition of the EPR spectrum for MD200 assuming the presence of three contributions. The spectra are detected at room temperature at 93.966 GHz.

paramagnetic centre with the small anisotropy of g-factor, although this model would fail to explain the difference in the behaviour of signal intensities that correspond to the particular g-factor components with the nanoparticle size. Anisotropic g-factor model will lead to a constant intensities ratio between the signal components which is not the case for our experiment.

Quantitative study of number of paramagnetic centers was implemented by exploiting X- and W-band EPR measurements as follows. First, in X-band we have determined a total spin concentration in the sufficient amount of MD1000 (ca. 15 mg) and compared that with the spin concentration in the reference sample. Comparison was done by exploiting the double integration of the CW spectra. In order to keep the microwave tuning parameters unchanged we used only MD1000: the samples with the smaller sizes shifted the tuning conditions away from reference sample conditions. Second, we have re-calculated the extracted in this way the total spin concentration to the mass of the MD1000 used for the W-band experiments, which was determined to be about 3.2-1018 spin per gram. No changes in the tuning-matching conditions with the MD sample sizes were obtained in the W-band experiments that probably could be explained by the small amount of substances within the W-band tubes (less than 1 mg) and, therefore, the influence of the species on the microwave field distribution in the cavity is negligible. Third, we have extracted the relative concentrations of the observed paramagnetic centers in a series of size-calibrated samples using MD1000 as a reference for the W-band. The results of the numerical measurements are presented in Figures 2 and 3.

Table 2. Spectroscopic parameters of SC1, SC2 and N0 centers obtained from the simulation of the W-band EPR spectra

Sample gSC1 Lorentzian linewidth, mT gSC2 Gaussian linewidth, mT gN0 Gaussian linewidth, mT A± ,mT A ,mT

MD10 2.0028(2) 2.0(1) 2.0029(2) 0.8(1) - - 2.92 4.08

MD20 2.0028(2) 1.9(1) 2.0029(2) 0.9(1) - - 2.92 4.08

MD50 2.0028(2) 2.0(1) 2.0029(2) 0.7(1) 2.0024 0.6 2.92 4.08

MD100 2.0028(2) 1.8(1) 2.0029(2) 0.7(1) 2.0024 0.6 2.92 4.08

MD200 2.0028(2) 1.8(1) 2.0029(2) 0.7(1) 2.0024 0.6 2.92 4.08

MD500 2.0028(2) 1.8(1) 2.0029(2) 0.7(1) 2.0024 0.6 2.92 4.08

MD1000 2.0028(2) 1.8(1) 2.0029(2) 0.7(1) 2.0024 0.6 2.92 4.08

O experiment

linear approximation

<n 1 -





1E-3 0.01 0.1

1/d, inverse particle size, nm-1

O experiment — — linear approximation











iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.


co0.1 7 '





i i 11 1E-3



—I—I—I I I I-





1/d, inverse particle size, nm"

Figure 2. Relative concentrations of SC1 (left panel) and SC2 (right panel) paramagnetic centers as a function of inverse particle size along with their linear approximations for particle sizes larger than 20 nm.

Figure 3. Relative concentrations of N0 paramagnetic centers as a function of inverse particle size along with approximation for particles larger than 50 nm. Solid blue line is drawn according to the core-shell model approximation with h = 20 nm; red vertical dashed line corresponds to d = 50 nm.

4. Discussion

In order to explain the dependencies observed in Figure 2 and 3, one can propose a core-shell model of nanoparticle structure that is schematically shown in Figure 4. This is a relatively standard model for nanoparticle spectroscopy interpretation. For instance, it was used in joint XRD and EPR spectroscopy measurements of rare-earth doped fluorite nanoparticles studied in the work [33]. We suppose that N0 centers are mainly concentrated in the core limited by the radius RC of the nanodiamond particle with the cross-section of diameter d while in the surface layer of thickness h other types of paramagnetic species prevail. The plots (figure 3, mainly) allow estimating the thickness as h « 20(5) nm.

The lack of EPR signal from N0 centers for MD20 and MD10 samples could be explained by either conversion of the neutrally charged nitrogen N0 to paramagnetic-silent charge state or due to the nitrogen flushing outside of the diamond lattice region of the particle due to its energy-induced

preferences [34]. Indeed, there are some reports for the nitrogen-vacancy (NV) charge conversion with the nanoparticle size [35, 36] or induced by surface chemical treatment [37] as well as by optical illumination [38]. Despite different electron spin density distribution for N0 and NV, the certain similarity in their electronic structure and energy level positions within the band gap makes it possible to allow the analogous behavior with the nanoparticle size.

The quantitative analysis of paramagnetic signals ascribed to surface layer(s) of particle is also of special interest. As stated above, the SC1 contribution was observed earlier in the study of deformed diamond crystals, but the SC2 component which is observed in this work was not detected directly in EPR spectra of nanodiamond before, except for instance the work of Nadolynny [24]. The EPR lineshape and width of SC1 center, along with almost perfect numerical agreement between total spin number of SC1 center and previously reported deformation-induced center in diamond [31] is probably a sufficient argument to relate this center to the mechanically-deformed paramagnetic species that were observed in diamonds previously. Linear dependence of the SC1 concentration with the inverse particle size at d > 50 nm (Figure 2) confirms homogeneous distribution of damage-induced paramagnetic centers over the geometrically distinct shell with the width h. The same behavior is to observe for SC2 centre (Figure 2). However, the MD10 and MD20 samples exhibit deviation of the SC1 and SC2 size-dependent concentration from the linear law, namely the concentration of SC2 is independent on the particle size while the EPR intensity for SC1 centre shows a sharp decrease. The origin of this decrement is not discussed in this work, but it is clear that the presented core-shell model is not an adequate one for the particles below 50 nm.

5. Conclusion

For the first time, based on the combination of conventional X-band and high-frequency (W-band) EPR measurements, the core-shell model of the HPHT-synthesized nanodiamond particle is proposed and shown to be in consistency with the numerous Raman, optical and NMR studies. Apparently, the nanoparticle consists of three distinct structural layers characterized by the specific paramagnetic centers. Nitrogen N0 centers are located within the diamond core of nanoparticle, SC1 paramagnetic center are filled the inner surface region with the thickness of about 20 nm of the deformed diamond lattice while SC2 centers are localized in the outer shell of mainly non-diamond carbon atoms.

The spin-Hamiltonian parameters for the SC1 center are in qualitative agreement with previously published papers devoted to the studies of the deformed synthetic diamond with the sizes from hundreds of micrometers down to 1 ^m. SC2 center exhibits purely surface behavior and its concentration is inversely proportional, in large extent, to nanoparticle size. It is shown that the amount of the paramagnetic nitrogen centers becomes negligible at HPHT diamonds dimensions of less than approx. 40-50 nm probably due to their conversion to the nonparamagnetic state.


Authors are grateful to Igor Vlasov from GPI RAS for providing MD samples for research. This work was partially funded by the subsidy allocated to Kazan Federal University for the project part of the state assignment in the sphere of scientific activities.

Figure 4. Schematic model of nanoparticle, where d is a cross section diameter (full size) of particle, Rc is the diamond core radius, and h is the width of deformed diamond lattice region


1. Osawa E., in Synthesis, Properties and Applications of Ultrananocrystalline Diamond, edited by D. Gruen, O. Shenderova and A. Vul', Springer, Netherlands (2005) p. 231.

2. Ösawa E. Diamond and Related Materials 16, 2018 (2007)

3. Fu C.-C., Lee H.-Y., Chen K., Lim T-S., Wu H.-Y., Lin P.K., Wei P.-K., Tsao P.-H., Chang H.-C., Fann W. Proceedings of the National Academy of Sciences, 104, 727 (2007)

4. Schrand A.M., Huang H., Carlson C., Schlager J.J., Ösawa E., Hussain S.M., Dai L., Journal of Physical Chemistry B 111, 2 (2007)

5. Chow E.K., Zhang X.-Q., Chen M., Lam R., Robinson E., Huang H., Schaffer D., Osawa E., Goga A., Ho D. Science TranslationalMedicine 3, 73ra21 (2011)

6. Dolde F., Fedder H., Doherty M.W., Nobauer T., Rempp F., Balasubramanian G., Wolf T., Reinhard F., Hollenberg L.C.L., Jelezko F., Wrachtrup J. Nature Physics 7, 459 (2011)

7. Kucsko G., Maurer P.C., Yao N.Y., Kubo M., Noh H.J., Lo P.K., Park H., Lukin M D. Nature, 500, 54 (2013)

8. Acosta V.M., Bauch E., Ledbetter M.P., Santori C., Fu K.-M.C., Barclay P.E., Beausoleil R.G., Linget H., Roch J.F., Treussart F., Chemerisov S., Gawlik W., Budker D. Physical Review B 80, 115202 (2009)

9. Osawa E. Pure and Applied Chemistry 80, 1365 (2008)

10. Vlasov I.I., Shiryaev A.A., Rendler T., Steinert S., Lee S.-Y., Antonov D., Voros M., Jelezko F., Fisenko A.V., Semjonova L.F., Biskupek J., Kaiser U., Lebedev O.I., Sildos I., Hemmer P.R., Konov V.I., Gali A., Wrachtrup J. Nature Nanotechnology 9, 54 (2014)

11. Stehlik S., Varga M., Ledinsky M., Jirasek V., Artemenko A., Kozak H., Ondic L., Skakalova V., Argentero G., Pennycook T., Meyer J.C., Fejfar A., Kromka A., Rezek B. Journal of Physical Chemistry C 119, 27708 (2015)

12. Shenderova O.A., McGuire G.E. Biointerphases 10, 030802 (2015)

13. Shames A.I., Panich A.M., Kempinski W., Alexenskii A.E., Baidakova M.V., Dideikin A.T., Osipov V.Y., Siklitski V.I., Osawa E., Ozawa M., Vul A.Y. Journal of Physics and Chemistry of Solids 63, 1993 (2002)

14. Balasubramanian G., Chan I.Y., Kolesov R., Al-Hmoud M., Tisler J., Shin C., Kim C., Wojcik A., Hemmer P.R., Krueger A., Hanke T., Leitenstorfer A., Bratschitsch R., Jelezko F., Wrachtrup J.,

Nature 455, 648 (2008)

15. Gruber A., Dräbenstedt A., Tietz C., Fleury L., Wrachtrup J., von Borczyskowski C. Science 276, 2012 (1997)

16. Panich A.M. Critical Reviews in Solid State and Materials Sciences 37, 276 (2012)

17. Ferrari A.C., Robertson J. Philosophical Transactions of the Royal Society of London A: Mathematical, Physical and Engineering Sciences 362, 2269 (2004)

18. Bradac C., Gaebel T., Naidoo N., Sellars M.J., Twamley J., Brown L.J., Barnard A.S., Plakhotnik T., Zvyagin A.V., Rabeau J.R. Nature Nanotechnology 5, 345 (2010)

19. Yavkin B.V, Soltamov V.A., Babunts R.A., Anisimov A.N., Baranov P.G., Shakhov F.M., Kidalov S.V., Vul A.Ya, Mamin G.V., Orlinskii S.B. Applied Magnetic Resonance 45, 1035 (2014)

20. Orlinskii S.B., Bogomolov R.S, Kiyamova A.M., Yavkin B.V., Mamin G.V., Turner S., Van Tendeloo G, Shiryaev A.A., Vlasov I.I., Shenderova O. Nanoscience and Nanotechnology Letters 3, 63 (2011)

21. Smith W.V., Sorokin P.P., Gelles I.L., Lasher G.J. Physical Review 115, 1546 (1959)

22. Loubser J.H.N., von Wyk J.A. Reports on Progress in Physics 41, 1201 (1978)

23. Baranov P.G., Soltamova A.A., Tolmachev D.O., Romanov N.G., Babunts R.A., Shakhov F.M., Kidalov S.V., Vul' A.Y., Mamin G.V., Orlinskii S.B., Silkin N.I. Small 7, 1533 (2011)

24. Nadolinny V.A., Golushko V.V., Palyanov Y.N., Shatsky V.S., Moseenkov S.I., Yuryeva O.P., Volodin A.M. Applied Magnetic Resonance 39, 295 (2010)

25. Baranov P.G., Il'in I.V., Soltamova A.A., Vul' A.Y., Kidalov S.V., Shakhov F.M., Mamin G.V., Orlinskii S.B., Salakhov M.K. JETP Letters 89, (2009)

26. Yavkin B.V., Gafurov M.R., Kharintsev S.S., Mamin G.V., Goovaerts E., Salakhov M.K., Osin Yu.N., Orlinskii S.B. Journal of Physics: Conference Series 478, 012001 (2013)

27. Yavkin B.V., Mamin G.V., Orlinskii S.B., Kidalov S.V., Shakhov F.M., Vul A.Ya., Soltamova A.A., Soltamov V.A., Baranov P.G. Applied Magnetic Resonance 44, 1235 (2013)

28. Fionov A.V., Lund A., Chen W.M., Rozhkova N.N., Buyanova I.A., Emel'yanova G.I., Gorlenko L.E., Golubina E.V., Lokteva E.S., Osawa E., Lunin V.V. Chemical Physics Letters 493, 319 (2010)

29. Stoll S., Schweiger A. Journal of Magnetic Resonance 178, 42 (2006)

30. Danil'chuck G.S., Ganyuk L.N., Koval'skii A.E., Pogoretskii P.P., Podzyarei G.A., Shul'man L.A.

Theoretical and Experimental Chemistry 1, 237 (1965)

31. Samsonenko N.D., Sobolev E.V. Journal of Experimental and Theoretical Physics Letters 5, 250 (1967)

32. Walters G.K., Estle T.L. Journal of Applied Physics 32, 1854 (1961)

33. Irisova I., Kiiamov A., Korableva S., Rodionov A., Tayurskii D., Yusupov R. Applied Magnetic Resonance 46, 515 (2015)

34. Gamarnik M Y. Physical Review B 54, 2150 (1996)

35. Petrakova V., Taylor A., Kratochvilova I., Fendrych F., Vacik J., Kucka J., Stursa J., Cigler P., Ledvina M., Fiserova A., Kneppo P., Nesladek M. Advanced Functional Materials 22, 812 (2012)

36. Rondin L., Dantelle G., Slablab A., Grosshans F., Treussart F., Bergonzo P., Perruchas S., Gacoin T., Chaigneau M., Chang H.C., Jacques V., Roch J.F. Physical Review B 82, 115449 (2010)

37. Hauf M.V., Grotz B., Naydenov B., Dankerl M., Pezzagna S., Meijer J., Jelezko F., Wrachtrup J. Stutzmann M., Reinhard F., Garrido J.A. Physical Review B 83, 081304(R) (2011)

38. Iakoubovskii K., Adriaenssens G.J., Nesladek M. Journal of Physics: Condensed Matter 12, 7843 (2000)

i Надоели баннеры? Вы всегда можете отключить рекламу.