About this Journal Submit a Manuscript Table of Contents
Advances in Condensed Matter Physics
Volume 2012 (2012), Article ID 439617, 10 pages
Research Article

Analysis of Erbium and Vanadium Diffusion in Porous Silicon Carbide

1Ioffe Physical-Technical Institute of the Russian Academy of Sciences, 26 Polytechnicheskaya Street St., Petersburg 194021, Russia
2Mathematical Department, Faculty of Radiophysics, Nizhny Novgorod State University, 65 Il’insky Street, Nizhny Novgorod 603950, Russia

Received 5 March 2012; Revised 1 June 2012; Accepted 15 June 2012

Academic Editor: Jörg Fink

Copyright © 2012 Marina G. Mynbaeva et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.


Experimental data on diffusion of erbium and vanadium in porous and nonporous silicon carbide at 1700 and 2200°C have been used for modelling diffusion in porous SiC. It is shown that the consideration of pore structure modification under annealing via vacancy redistribution allows for satisfactory description of dopant diffusion. As expected, important contribution to the diffusion in the porous medium is found to be made by the walls of the pores: in SiC, the vacancy surface diffusion coefficient on the walls appears to exceed that in the bulk of the material by an order of magnitude. When thermal treatment transforms pore channels into closed voids, pathways for accelerated diffusion cease to exist and diffusion rates in porous and nonporous SiC become similar.

1. Introduction

Silicon carbide (SiC) is a semiconductor that possesses properties highly suitable for high-power, high-frequency, and high-temperature microelectronic applications. For a long time, selective doping in SiC technology has been implemented mostly with ion implantation. Diffusion doping was considered to be impractical due to low diffusion coefficients of impurities at temperatures below 1800°C, and the interest in diffusion in SiC was generally related to applications requiring “drive-in” annealing after implantation (see, e.g., [13]). Papers devoted to “classic” thermal diffusion in SiC were mostly limited to early years [4, 5]. However, recent achievements in fabrication of SiC-based microelectronic devices using diffusion doping [1, 2, 6] stirred renewed interest in diffusion in SiC.

Diffusion in SiC is believed to proceed mostly via kick-out or vacancy-related mechanisms [15]. Vacancies define many physical properties of SiC [7, 8], and one of the current subjects of research is their contribution to the formation and properties of porous SiC, a material with many prospective applications [9, 10]. The very formation of porous structure in SiC is believed to proceed via vacancy migration and clustering [11], and there are all grounds to believe that experimentally observed modification of the structure under thermal annealing [12, 13] is also due to vacancy redistribution. In this respect, diffusion in porous SiC represents a subject of special interest, as pores serve as the sources of vacancies, which should contribute to diffusion. In this paper, we consider diffusion in SiC with porous structure formed by the vacancy mechanism, and model a dynamic process, where diffusion of dopants proceeds against the background of fluxes of vacancies responsible for pore ripening during high-temperature annealing.

2. Experimental Details

To form the porous structure, 6H–SiC wafers were subjected to surface anodization in HF. The anodization mode was chosen such as to ensure formation of pores via vacancy-related mechanism, with no chemical etching of the surface involved [11]. As a result of anodization, porous layers were formed with thickness significantly greater than the expected diffusion depth. The dopants were erbium and vanadium. Erbium doping of SiC is attractive for infrared applications relying on radiative intrashell transitions in Er ions placed in semiconductor lattice [14], while vanadium is a dopant of choice for fabricating semi-insulating SiC [15]. In each experiment, the dopant was simultaneously introduced into SiC layer with a “nanoporous” structure, and into a reference sample of nonporous SiC. The diffusion was carried out in the so-called sandwich-type cell. Conditions in the cell prevent the crystal surface from decomposition in the course of high-temperature annealing, thus making possible quasi-equilibrium diffusion doping of nonporous SiC at temperatures °C [5]. The diffusion was carried out at and 2200°C for a period of time of 30 and 120 min. The dopant profiles were investigated using secondary ion mass spectroscopy (SIMS).

3. Experimental Results

To follow changes in the porous structure induced by the high temperature, a control set of samples was first subjected to thermal anneal under the same conditions, which were later used for diffusion. Figure 1 shows a typical example of the modification of porous structure in 6H-SiC as a result of anneal at °C for 120 and 240 min. As can be seen, initial porous structure, which appeared in a form of branched thin pores with the pore stem diameter of 30 to 40 nm (a), after 120 min of anneal evolved into a set of “cylinders” with the diameter of 100 to 200 nm (b). Further annealing transformed the “cylinders” into closed voids, with some voids showing tendency to faceting (c). The pattern of porous structure after anneal at 2200°C for 120 min was similar to that shown in Figure 1(c), as with temperature increasing the modification of porous structure proceeds faster [13]. According to Cheremskoy et al. [16], the closing and faceting of the voids confirm that pore modification proceeds according to ripening mechanism, and it is due to vacancy redistribution.

Figure 1: Transmission (a) and scanning ((b) and (c)) electron microscopy images of the cleft edge of porous layer in SiC: straight after anodization (a), after 1700°C anneal for 120 min (b), and after 1700°C anneal for 240 min (c).

SIMS measurements with depth profiling, performed both after the above anneals and dopant diffusion, did not reveal any serious deviation from SiC stoichiometry in the processed samples. The examples of dopant diffusion profiles are presented in Figure 2. As was shown elsewhere [17], only a small part of such profiles related to the area close to the surface could be fitted with complementary error function, which testified to the concentration-dependent diffusion. Figure 2(a) presents erbium profiles in porous and nonporous SiC after diffusion at 2200°C for 30 and 120 min. As is seen, the erbium surface concentration in the samples was close to 3 × 1018 cm3. For samples diffused for 30 min, we observed a difference in the diffusion depth in porous and nonporous specimens. For samples diffused for 120 min, this difference became somewhat larger, but considering the difference in the surface concentration, not as large as could be expected if the “classic” inverse square root diffusion law was applicable.

Figure 2: Experimental profiles showing erbium (a) and vanadium (b) concentration versus depth after diffusion at 2200°C for 30 min (curves 1 and 2) and 120 min (curves 1′ and 2′), in: 1: nonporous SiC, 2: porous SiC.

A similar situation was observed with vanadium diffusion at 2200°C. For all the diffused samples, the vanadium surface concentration was ~2 × 1017 cm−3 (Figure 2(b)). The short-time (30 min) diffusion yielded diffusion depth in porous SiC larger than that in nonporous material. Increasing the anneal time to 120 min, the difference in diffusion depth in porous and nonporous specimens remained almost the same.

4. Modelling

For modelling, let us consider diffusion of impurity atoms into a porous layer. The diffusion starts at the moment of time and proceeds through an interface with . It can be described with the second Fick law written in the following form: with the following boundary and initial conditions: where is the Boltzmann constant, is molar volume; is the chemical potential, and are the initial and final pore volumes, respectively, is the molar gas constant, , , and are the dimensions of doped sample, is the impurity concentration, and and are the diffusion coefficients in the bulk and on the surface (pore walls) of the material, respectively (see, e.g., [18]). As follows from the experimental data, we should consider the dependence of and on the vacancy and dopant concentration. These dependences can be described using the functions suggested by Gotra [19], Zorin et al. [20], and Ryssel and Ruge [21], respectively, Here, and (where index relates to diffusion in bulk and relates to surface diffusion) are linear (i.e., independent on dopant concentration, the linearity being indicated by index ) coefficients of the approximations. These (as well as solubility limits) have been extracted from the experimental data presented earlier [17]. The second multipliers in (3) describe concentration dependence of diffusion coefficients [19]. is the dopant solubility limit. Value of the parameter (which is equal to an integer value in the range of 1 to 3 [19]) depends on the interaction between dopant atoms and defects [19, 22]. The third multipliers in (3) describe vacancy dependences of the diffusion coefficients [20, 21]. is the equilibrium vacancy distribution. Four fitting parameters, , , , can be considered as scale multipliers of the limit of solubility of dopant and equilibrium vacancy concentration. The vacancy concentration can be described as follows: with the initial and boundary conditions: where designates initial pore surface, stands for the surface of the outward boundary of the porous area, is the absolute value of the velocity of the movement of the surface of the growing or shrinking pore [18], and are volumetric and surface diffusion coefficients of vacancies; is the specific surface energy [16], and , where is interatomic distance. The appropriateness of using (4) for describing the evolution of a pore under annealing was considered in detail by Kitayama et al. [18]. The initial and boundary conditions for this equation represent three-dimensional generalization of those given by Cheremskoy et al. [16], with the consideration of pore size changes.

A typical vacancy distribution in the vicinity of a pore, which was obtained as a result of modelling using (4) is shown in Figure 3. Here, “x” is the coordinate across the pore: the maximum value of vacancy concentration corresponds to the center of the pore. The initial pore was regarded as a cylinder with the radius of 40 nm. The parameters of modelling were chosen so that the simulated annealing would lead to the growth of pore “cylinder” with the diameter of ~200 nm, which corresponds to the experimental data presented in Figure 1(b). The figure shows how annealing leads to the increase in pore diameter, and confirms qualitative adequacy of our model for consideration of pore modification under heat treatment on the microscopic scale.

Figure 3: Modelled vacancy distribution in the vicinity of a pore in porous SiC before (curve 1) and after (curve 2) anneal at 1700°C for 120 min.

To consider the diffusion, let us transform the diffusion coefficients into the following form: and , where and are average values of diffusion coefficients, , , , , (see, e.g., [2224]). We determined solutions of (1) and (4) as the following power series: Substitution of the series (6) into (1), (4), and appropriate boundary and initial conditions gives us possibility to obtain equations for zero-order approximation of dopant and vacancy concentration and , their corrections and , and the boundary and initial conditions for all functions and : Solutions of (7) with conditions (8) are presented in the appendix. Though the values of the vacancy and dopant concentration can be determined with more accuracy, second-order approximation usually gives satisfactory qualitative and some quantitative results (see, e.g., [2224]).

5. Discussion

Let us analyze the experimental diffusion profiles using the relations obtained in the previous section. We will restrict our consideration to one dimension only, which is inward the porous layer (direction x). Considering the experimental results, we may assume that the results of 30 min-long anneals were somehow affected by transient processes, yet 120 min-long diffusion proceeded under “steady-state” conditions. This assumption is supported by the fact that the propagation of the diffusion fronts as presented in Figure 2 did not obey the inverse square root law. We fitted experimental profiles shown in Figure 2 by minimization of mean-squared error between experimental and theoretical results. The solid lines correspond to erbium profiles in porous and nonporous SiC after diffusion at 1700°C according to the procedure developed above for 120 min, that is, long diffusion anneals. The parameters obtained are listed in Table 1 and the fitting curves are shown in Figure 4. Values of fitting parameters are presented in the Table 2. The parameters , , , and give us corrections to experimental values of dopant solubility limit and equilibrium vacancy distribution. Their usefulness is determined by the fact that one can obtain quite a large difference between experimental data taken from different references.

Table 1: Diffusion parameters obtained from the fitting.
Table 2: Fitting parameters for Figure 4.
Figure 4: Distribution in non-porous SiC (curves 1) and porous SiC (curves 2) of: erbium after diffusion at 1700°C for 120 min (a), vanadium after diffusion at 2200°C for 120 min (b), and erbium after diffusion at 2200°C for 120 min (c). Points show experimental data, solid curves show results of fitting.

It follows from the results obtained that the vacancy diffusion coefficient on the surface of the pore walls exceeds that in the bulk of the material approximately by an order of magnitude. Erbium appears to diffuse in SiC much faster than vanadium, seemingly due to stronger concentration dependence of its diffusion coefficient, that is, concentration dependence of diffusion coefficient leads to higher diffusion coefficient for erbium in comparison with that of vanadium. Summarizing the whole set of diffusion data obtained (including our earlier data [17]), it appears that when a short-term diffusion anneal is used, diffusion in porous SiC proceeds faster than that in nonporous SiC simply because porous structure does not have enough time to evolve into closed “cylinders.” For long-term anneals, pathways for accelerated diffusion (pore walls) disappear as pore channels get transformed into closed voids, and diffusion coefficients in porous and nonporous SiC become close. This effect becomes more pronounced with increasing diffusion temperature. At low temperatures °C, the diffusion in porous SiC always proceeds much faster than that in nonporous SiC due to weak pore structure modification.

6. Conclusion

Modelling based on experimental data on the diffusion of erbium and vanadium in porous silicon carbide at 1700 and 2200°C shows that this process is greatly affected by the diffusion of vacancies, which is a typical process in the material with vacancy-related porosity under thermal treatment. Consideration of the vacancy redistribution under annealing allows for satisfactory describing the dopant diffusion in porous SiC and yields basic diffusion parameters. The vacancy diffusion coefficient on the surface of the pore walls appears to exceed that in the bulk of the material by an order of magnitude. The values of “bulk” diffusion coefficients obtained in this work agree with the data obtained for typical dopants in SiC (see, e.g., [25]). Generally, it appears that the consideration of vacancy redistribution during high-temperature annealing allows for correct description of diffusion in a porous semiconductor.


Equations for the functions and have been solved by the standard Fourier approach (see, e.g., the work by Carslaw and Jaeger [26]). Solutions of (7) with conditions (8) can be written in the form: where , where ,


This work was partly supported by the grant of the President of Russia (Project no. MK-548.2010.2).


  1. A. V. Bolotnikov, P. G. Muzykov, and T. S. Sudarshan, “Investigation of two-branch boron diffusion from vapor phase in n-type 4H-SiC,” Applied Physics Letters, vol. 93, no. 5, Article ID 052101, 3 pages, 2008. View at Publisher · View at Google Scholar
  2. A. V. Bolotnikov, P. G. Muzykov, A. E. Grekov, and T. S. Sudarshan, “Improvement of 4H-SiC power p-i-n diode switching performance through local lifetime control using boron diffusion,” IEEE Transactions on Electron Devices, vol. 54, no. 6, pp. 1540–1544, 2007. View at Publisher · View at Google Scholar · View at Scopus
  3. E. Friedland, J. B. Malherbe, N. G. van der Berg et al., “Study of silver diffusion in silicon carbide,” Journal of Nuclear Materials, vol. 389, no. 2, pp. 326–331, 2009. View at Publisher · View at Google Scholar · View at Scopus
  4. A. Y. Vodakov and E. N. Mokhov, “Diffusion and solubility of impurities in silicon carbide,” in Proceedings of the 3rd International Conference on Silicon Carbide 1973, R. C. Marshall, J. W. Faust Jr., and C. E. Ryan, Eds., pp. 508–520, South Carolina University Press, Florida, Fla, USA, 1974.
  5. A. Y. Vodakov and E. N. Mokhov, “Point defects in silicon carbide,” Institute of Physics Conference Series, vol. 137, pp. 197–206, 1994.
  6. Y. Gao, S. I. Soloviev, and T. S. Sudarshan, “Planar 4H- and 6H-SiC p-n diodes fabricated by selective diffusion of boron,” Solid-State Electronics, vol. 45, no. 12, pp. 1987–1990, 2001. View at Publisher · View at Google Scholar · View at Scopus
  7. M. G. Mynbaeva, K. D. Mynbaev, V. A. Ivantsov, A. A. Lavrent'ev, B. A. Grayson, and J. T. Wolan, “Semi-insulating porous SiC substrates,” Semiconductor Science and Technology, vol. 18, no. 6, pp. 602–606, 2003. View at Publisher · View at Google Scholar · View at Scopus
  8. R. Rurali, E. Hernández, P. Godignon, J. Rebollo, and P. Ordejón, “First-principles studies of the diffusion of B impurities and vacancies in SiC,” Physical Review B, vol. 69, no. 12, Article ID 125203, 2004. View at Scopus
  9. R. M. Feenstra and C. E. C. Wood, Eds., Porous Silicon Carbide and Gallium Nitride: Epitaxy, Catalysis, and Biotechnology Applications, John Wiley & Sons, Chichester, UK, 2008.
  10. M. G. Mynbaeva and K. D. Mynbaev, “Technological Applications of Porous SiC,” in Nanocrystals and Quantum Dots of Group IV Semiconductors, T. V. Torchinska and V. Yu. Vorobiev, Eds., pp. 255–276, American Scientific, New York, NY, USA, 2010.
  11. M. G. Mynbaeva, D. A. Bauman, and K. D. Mynbaev, “On the role of vacancies in pore formation in the course of anodizing of silicon carbide,” Physics of the Solid State, vol. 47, no. 9, pp. 1630–1636, 2005. View at Publisher · View at Google Scholar · View at Scopus
  12. J. Bai, G. Dhanaraj, P. Gouma, M. Dudley, and M. Mynbaeva, “Porous SiC for HT chemical sensing devices: an assessment of its thermal stability,” Materials Science Forum, vol. 457–460, no. 2, pp. 1479–1482, 2004. View at Scopus
  13. M. Mynbaeva, A. Lavrent'ev, I. Kotousova, A. Volkova, K. Mynbaev, and A. Lebedev, “On current limitations in porous SiC applications,” Materials Science Forum, vol. 483–485, pp. 269–272, 2005. View at Scopus
  14. V. Glukhanyuk and A. Kozanecki, “Site selective studies of Er3+ emission centers in Er-implanted 6H-SiC,” Applied Physics Letters, vol. 89, no. 21, Article ID 211114, 2006. View at Publisher · View at Google Scholar · View at Scopus
  15. N. T. Son, P. Carlsson, A. Gällström, B. Magnusson, and E. Janźn, “Deep levels and carrier compensation in V-doped semi-insulating 4H-SiC,” Applied Physics Letters, vol. 91, no. 20, Article ID 202111, 2007. View at Publisher · View at Google Scholar · View at Scopus
  16. P. G. Cheremskoy, V. V. Slezov, and V. I. Betekhtin, Pores in Solids, Energoatomizdat, Moscow, Russia, 1990.
  17. M. G. Mynbaeva, E. N. Mokhov, A. A. Lavrent'ev, and K. D. Mynbaev, “High-temperature diffusion doping of porous silicon carbide,” Technical Physics Letters, vol. 34, no. 9, pp. 731–733, 2008. View at Publisher · View at Google Scholar · View at Scopus
  18. M. Kitayama, T. Narushima, and A. M. Glaeser, “Wulff shape of alumina: II, experimental measurements of pore shape evolution rates,” Journal of the American Ceramic Society, vol. 83, no. 10, pp. 2572–2583, 2000. View at Scopus
  19. Z. Yu. Gotra, Technology of Microelectronic Devices, Radio and Communications, Moscow, Russia, 1991.
  20. E. I. Zorin, P. V. Pavlov, and D. I. Tetelbaum, Ion Doping of Semiconductors, Energiya, Moscow, Russia, 1975.
  21. H. Ryssel and I. Ruge, Ion Implantation, Teubner, Stuttgart, Germany, 1978.
  22. E. L. Pankratov, “Influence of the spatial, temporal, and concentrational dependence of the diffusion coefficient on dopant dynamics: optimization of annealing time,” Physical Review B, vol. 72, no. 7, Article ID 075201, 2005. View at Publisher · View at Google Scholar · View at Scopus
  23. E. L. Pankratov and B. Spagnolo, “Optimization of impurity profile for p-n-junction in heterostructures,” European Physical Journal B, vol. 46, no. 1, pp. 15–19, 2005. View at Publisher · View at Google Scholar · View at Scopus
  24. E. L. Pankratov, “Dopant diffusion dynamics and optimal diffusion time as influenced by diffusion-coefficient nonuniformity,” Russian Microelectronics, vol. 36, no. 1, pp. 33–39, 2007. View at Publisher · View at Google Scholar · View at Scopus
  25. V. A. Dmitriev and M. G. Spencer, “SiC fabrication technology: growth and doping,” Semiconductors and Semimetals, vol. 52, pp. 21–75, 1998. View at Publisher · View at Google Scholar · View at Scopus
  26. H. S. Carslaw and J. C. Jaeger, Conduction of Heat in Solids, Oxford University Press, Oxford, UK, 1964.