Seebeck coefficient of silicon nanowire forests doped by thermal diffusion

  1. Shaimaa Elyamny1,2ORCID Logo,
  2. Elisabetta Dimaggio1ORCID Logo and
  3. Giovanni Pennelli1

1Dipartimento di Ingegneria della Informazione, Università di Pisa, Via G.Caruso, I-56122 Pisa, Italy
2Electronic Materials Research Department, Advanced Technology and New Materials Research Institute, City of Scientific Research and Technological Applications (SRTA-City), New Borg El-Arab City, 21934, Alexandria, Egypt

  1. Corresponding author email

Associate Editor: E. Meyer
Beilstein J. Nanotechnol. 2020, 11, 1707–1713.
Received 07 Aug 2020, Accepted 14 Oct 2020, Published 11 Nov 2020


Thermoelectric generators made by large arrays of nanowires perpendicular to a silicon substrate, that is, so-called silicon nanowire forests are fabricated on large areas by an inexpensive metal-assisted etching technique. After fabrication, a thermal diffusion process is used for doping the nanowire forest with phosphorous. A suitable experimental technique has been developed for the measurement of the Seebeck coefficient under static conditions, and results are reported for different doping parameters. These results are in good agreement with numerical simulations of the doping process applied to silicon nanowires. These devices, based on doped nanowire forests, offer a possible route for the exploitation of the high power factor of silicon, which, combined with the very low thermal conductivity of nanostructures, will yield a high efficiency of the conversion of thermal to electrical energy.

Keywords: nanowires; Seebeck coefficient; thermal conductivity; thermoelectricity


Thermoelectric generators for direct conversion of heat into electrical power will certainly play a decisive role in the next generation of energy harvesting and energy scavenging systems. However, a large-scale application of thermoelectric devices requires the development of materials that have good thermoelectric features and are, at the same time, of low cost, technologically affordable and sustainable. Silicon has a very high power factor S2σ [1-4] (S is the Seebeck coefficient and σ is the electrical conductivity). This, combined with the reduced thermal conductivity when nanostructured [5-10], makes it very suitable for thermoelectric applications. As added value, silicon is cheap and abundant. Also, there are established technologies for processing Si and Si is biocompatible. Furthermore, the use of silicon for thermoelectric generator devices will make them technologically compatible with standard CMOS devices. The main requirement for the use of silicon as thermoelectric material is the development of techniques for the low-cost fabrication and interconnection of a large number of nanostructures to generate a significant amount of power. Metal-assisted chemical etching (MACE) [11-14] of silicon is very promising because it gives the opportunity to fabricate large numbers of nanowires with high aspect ratio, perpendicular to a silicon substrate, that is, so-called silicon nanowire (SiNW) forests. The process is very suitable for the large-scale fabrication of nanostructured devices useful for several applications, such as sensing, photovoltaics, energy storage (supercapacitors), and, in particular, thermoelectric applications [15-17]. There are two main requirements to the fabrication of a leg for a thermoelectric generator that is based on a large number of silicon nanowires perpendicular to a substrate: 1) Electrical contacts need to fabricated on the top of a silicon nanowire forest, which can be achieved by copper electrodeposition [18]. 2) The optimum doping concentration of the nanowires for the exploitation of the maximum power factor of silicon [3] needs to be found. Both the Seebeck coefficient and the electrical conductivity depend on the doping concentration. In particular, S decreases with increasing doping concentrations. At the same time, high doping concentrations yield high values of electrical conductivity σ. Hence, a trade-off between S and σ needs to be found, such that the product S2σ is maximized. Nanowires with an average diameter of 80 nm and a length of several hundreds of micrometers can easily be obtained from MACE on n- and p-doped substrates, with concentrations up to 1018 cm−3. However, the maximum of the power factor is achieved for doping concentrations greater than 1019 cm−3 [2,3]. Unfortunately MACE yields very unreliable results at such high doping concentrations. In the case of p+ doping, parameters for a satisfactory fabrication of monocrystalline SiNW forests have been found [17]. However, the reliable fabrication of n+ SiNWs by MACE is still an open issue [19,20]. Hence, it is currently not possible to fabricate an optimized generator module based on two legs with opposite heavy doping.

A possible solution is to dope the silicon nanowires by thermal diffusion [21], after their fabrication by MACE. In this work, we present the measurement of the Seebeck coefficient of SiNWs doped through a diffusion process, based on a solid source (see Methods section). The main point is that the diffusion process must be performed in a single step (predeposition step), because a drive-in step, typical of the standard diffusion processes currently applied in the semiconductor industry, would require an oxidation of the surface for the trapping of the doping species. However, the reduction of the thermal conductivity, which is the aim of the nanoscale structuring, relies on the roughness of the nanowire surfaces, and the smoothing produced by the oxide growth would heavily reduce this effect. Therefore, it is mandatory to perform a single-step diffusion process, which results in a nonuniform doping concentration in the nanowire. Here, we report the measurement of the Seebeck coefficient after different doping processes, and compare the measurements with numerical simulations that take into consideration the nonuniform doping of the silicon nanowires. We found a very good agreement between experimental measurements and simulations of the doping process.


Fabrication and doping of silicon nanowire forests

Silicon nanowire forests have been fabricated by a simple and inexpensive process based on one-pot metal-assisted chemical etching (MACE) [22] (Figure 1).


Figure 1: (Left) Schematics of the fabrication process and (right) SEM images of silicon nanowire forests. The forests consist of nanowires with a length of several tens of micrometers, placed perpendicularly to the silicon substrate.

Silicon chips of roughly 1 × 1 cm2 have been cut from n-doped (phosphorous) commercial silicon ⟨100⟩ wafers with a nominal resistivity of 10 Ω·cm (nominal doping concentration 1015 cm−3). The chips, mounted on a custom-made apparatus to provide the stirring during the etch, have been soaked in a HF/AgNO3 solution (HF (48%)/AgNO3 (0.1 M)/H2O, volume ratio 16:5:60, see Figure 1). The stirring is crucial for the uniformity of the etching. Hence, the custom-made apparatus that holds the sample changes at random (every few seconds) the direction and the speed of the rotation. During the etch, a temperature-controlled bath allowed to maintain a stable temperature of 18 ± 0.5 °C. The etching time determines the final length of the nanowires. Several SiNW forests with nanowire lengths between 6.5 μm (30 min etching) and 41 μm (3 h etching) have been fabricated for this work. At the end of the MACE process step, the chips were covered with granular silver, which has been removed by etching in a HNO3/H2O 1:1 solution for 2 min. SEM inspection confirmed the complete removal of silver from the SiNW forests, which are covered by a thin layer of oxide as a result of the HNO3 etching.

The doping of the silicon nanowires has been carried out by thermal diffusion from a solid source. At first, the chips with the SiNW forests, with a surface of roughly 1 × 1 cm2 have been cleaned in buffered HF (BHF) for 1 min, to remove the SiO2 grown during the HNO3 etching. The chips have then been placed in a tubular quartz furnace together with the solid source. As a solid source, we used ceramic wafers provided by Techneglass (PhosPlus TP-250). The face of the chips with the nanowires has been placed in contact with the ceramic wafer. The sealed tube has been cleaned with a flux of pure nitrogen for several minutes; the nitrogen flux has been maintained during the whole thermal doping process. The temperature has been raised to the target temperature, with a ramp of 20 °C/min. Once the target temperature had been reached, it has been maintained for the chosen doping time (ten minutes, typically). Then, the chips have been allowed to cool maintaining the nitrogen flux.

After doping by thermal diffusion, a contact has been provided on top of the silicon nanowire forest exploiting the copper electrodeposition method described in a previous publication [18]. At first, a Cr (for adhesion)/Cu double metal layer has been deposited by thermal evaporation both on top of the nanowires and on the other side of the silicon wafer, that is, at the bottom of the substrate. The Cu layer on the top has been used as a seed for the electrodeposition of copper, performed at a constant current of 800 A/m2 for 2 min.

At the end of the process, we obtained a SiNW forest placed between a top copper contact and the silicon substrate, contacted through the bottom metal layer. The doping of the SiNW forest depends on the diffusion process (temperature and time), while the substrate remains slightly doped, with a nominal resistivity of 10 Ω·cm. The result is a leg of a silicon-based thermoelectric generator, as shown in Figure 1.

Measurement of the thermal conductivity and of the Seebeck coefficient

Single-leg thermoelectric generators have been characterized with a measurement apparatus based on the guarded hot-plate technique, schematically shown in Figure 2.


Figure 2: Left: sketch of the measurement apparatus. Right: thermal resistance of several samples as a function of the nanowire length (doping at 800 °C for 10 min).

This apparatus has already been used in a previous publication for the measurement of the thermal conductivity of undoped silicon nanowire forests [23]. Essentially, the Si chip with the contacted nanowire forest is squeezed between two aluminum blocks maintained at different temperatures TH = THot and TC = TCold. Given the surface A of the sample, the squeezing pressure has been set to 20 MPa for all measured samples. A third block THRef is maintained at the same temperature as TH (within ±0.1°), so that the electrical power needed to heat TH through Joule heating is the thermal power [Graphic 1] = V × I that crosses the sample in the absence of heat losses. The system is in vacuum and the irradiation loss is limited by a thermal shield maintained at TH. Hence, the thermal resistance ℜ can be evaluated as ℜ = (THTC)/[Graphic 2] Finite-element modeling has been used to refine the temperature measurement (see our previous publication [23] for more details). Samples of different lengths, obtained after different etching times, and with different doping temperatures, have been fabricated and measured. The plot of Figure 2 shows ℜ × A for several samples, where A is the area of the sample, as a function of the nanowire length.

Each point is the average of multiple measurements on the same sample, achieved with different thermal power values, [Graphic 3] = V × I, between 0.5 and 2 W. Each measurement has been made under static conditions. The electrical power V × I has been applied to the TH heater, and an automatic system controlled the heater of the THRef block, so that THref was within ±0.1° from TH. Once the electrical power V × I had been applied, a sufficiently long time (roughly 1 h) has been waited for the thermal stabilization of the system, before recording TH, TC and the thermal power [Graphic 4] through the sample. The linear fit of ℜ × A as a function of the nanowire length L is also reported in the plot of Figure 2. The reciprocal of the slope is the thermal conductivity kt, multiplied by the coverage factor (or filling factor) ν, where ν is the ratio between the total cross-section surface of the nanowires and the surface of the sample. From the plot, a value of kt = 4.2 ± 0.4 results for the samples doped at 800 °C.

For the estimation of the filling factor, several SEM images have been taken of different locations of each sample. The area of the samples is roughly 1 × 1 cm2, each image covers an surface of at least 50 × 50 μm2. A software for image reduction (ImageJ) has been used to determine the filling factor in each image, and an average has been calculated repeating the procedure for at least ten images. The resulting average filling factor was ν = 0.3 ± 0.02. These values are comparable, within the experimental errors, with that measured on undoped samples, reported in our previous work [8]. The intercept with the vertical axis of the linear fits, shown in the plots of Figure 2, is the thermal resistance of the contacts ℜCA = (1.8 ± 0.6) × 10−5 m2·K/W. It has been assumed that the contact thermal resistance ℜCA was the same for all the samples, since it depends on the mechanical pressure, which has been set to 20 MPa during all measurements.

The thermal resistance of the contacts ℜCA is fundamental for the evaluation of the Seebeck coefficient S. Simultaneously to the measurement of the thermal resistance of each sample, the voltage between the aluminum blocks TH and TC has been measured (Seebeck voltage Vseebeck for THTC). However, the temperature difference ΔT = THTC includes also the temperature drop due to the contact thermal resistance, that is, ΔT = ΔTNW + ΔTCONTACTS. Hence, the evaluation of S simply as the ratio between the measured values of VS and ΔT (S = VST) would give an underestimated value, because ΔTNW is smaller than the measured value of ΔT = THTC. The measured value of the contact thermal resistance ℜCA, together with the value of the heat flux through the sample [Graphic 5] = V × I, can be used for a correct evaluation of ΔTNW:


For each sample, several Seebeck voltage values Vseebeck have been recorded for different heating power values [Graphic 6], that is, for different values of ΔT. As explained in the description of the thermal resistance, each measurement was carried out under static conditions, allowing for a thermal stabilization of the system before recording the data. The temperature drop across the nanowires, ΔTNW, has then been evaluated exploiting the value of the thermal contact resistance.

Results and Discussion

Figure 3 shows the measured voltage drop as a function of the temperature difference between the ends of the nanowires, calculated as described above. The slope of the linear fit gives the Seebeck coefficient, reported in the legend of the figure for the different values of the doping process parameters (doping temperature and time). The results are S = −0.88 mV/K for the undoped sample (nominal resistivity 10 Ω·cm, estimated doping concentration 1015 cm−3), S = −0.41 mV/K for the sample doped at 700 °C for 10 min, and S = −0.20 mV/K for the sample doped at 800 °C for 10 min.


Figure 3: Measured voltage drop (Seebeck voltage ΔVseebeck) as a function of the temperature difference between the ends of the nanowires, ΔTNW. The slope is the Seebeck coefficient of the nanowires.

In the case of undoped samples, it is very difficult to establish the final charge carrier concentration, because it is strongly affected by the surface states of the nanostructures. In the case of doped samples, it is presumable that the doping concentration inside the nanowires is not uniform, because phosporous diffuses from the surfaces generating an exponential (error function-like) profile inside the nanowire. It must be noted that the standard diffusion doping process, used in the fabrication of integrated devices, consists of two steps. First, there is a predeposition step, in presence of a phosphorous-rich atmosphere, which results in an exponential doping profile (such as the doping performed here). Second, there is a so-called drive-in process in oxidizing atmosphere. This second step, performed at temperatures in excess of 1000 °C, is necessary to allow for the penetration of the doping species (phosphorous, in our case) into the silicon. It requires an oxidizing environment, at least in a first preliminary phase, to grow a thin SiO2 layer at the surface as a barrier for the doping species, forcing the diffusion into silicon. In the specific case of thermoelectric applications, we cannot afford an oxidization process, because it would smooth the surfaces. The reduction of the thermal conductivity, which is essential for good thermoelectric features, relies on the roughness of the surfaces, which would be compromised by an oxidation step. Hence, the best solution is to perform the doping using only the predeposition step.

Given the doping parameters, we gave an estimation of the doping distribution into the nanowires by simulating the diffusion process. Phosphorous diffusion in silicon is a well-known process, widely characterized for its importance in the fabrication of integrated electronic devices and circuits [24].

The diameter of the nanowires fabricated by MACE is not uniform, but it is distributed around an average value of 80 nm, which depends on the etching parameters. The transport properties (thermal conductivity and Seebeck coefficient) have been measured on macroscopic samples (several square millimeters of surface). Therefore, their value is averaged over a large number of nanowires with different diameters. For the simulations of the doping diffusion process, we considered the average diameter. Figure 4 shows finite-element simulations of doping diffusion in nanowires with a diameter of 80 nm and length of several micrometers.


Figure 4: Numerical simulations of phosphorous diffusion in silicon nanowires. The doping concentration is maximum on the surfaces and decreases exponentially along the radius toward the center.

The simulations were carried out solving the 2D diffusion equation [Graphic 7] = −DPND(x,y,t) in the circular cross section of a typical nanowire, where ND(x,y,t) is the doping concentration as a function of the position and of the time, t is the time and DP is the diffusion coefficient of phosphorous at the given temperature T, evaluated as [Graphic 8], with D0 = 0.79 cm2/s and EP = 3.29 eV, following [24]. The radial doping profiles are also reported in Figure 4. The Seebeck coefficient of the doped nanowires is difficult to evaluate because the doping is not uniform. There are experimental papers [2,25-27] reporting the Seebeck coefficient S = S(n) measured on bulk silicon with uniform doping concentration n. As shown in our previous work [3], the best logarithmic fit of the reported experimental data is S(n) = 0.0077 – 1.26 × 10−4 ln(n) (n is the doping concentration in m−3, S as absolute value). This relationship between S = S(n) and the doping concentration can be used to evaluate the Seebeck coefficient in our nanowires with nonuniform doping. The doping profile has been assumed as n(x,y) = ND(x,y), where ND(x,y) is the result of the simulation of the doping process. Hence, the Seebeck coefficient depends on the position in the nanowire. That is, S = S(n(x,y)), where (x,y) is a generic point in the cross section. Also the electrical conductivity depends on doping, that is, σ(x,y) = σ(n(x,y)) It has been estimated using the formula by Arora [28], which is widely used in silicon device simulations. The Seebeck coefficient of the nanowires has then been calculated as:


This formula can be easily derived by considering the nanowires as many parallel thermoelectric generators SiΔT, each with its resistance Ri (conductance Gi = 1/Ri). It is straightforward to obtain [Graphic 9]. Taking into account the distributed parameters we derive the formula above. The simulations yield (as absolute values) S = 0.448 mV/K for the nanowires doped at 700 °C for 10 min, and S = 0.207 mV/K for the nanowires doped at 800 °C for 10 min. This values are in very good agreement with the experimental results.


The development of efficient silicon-based thermoelectric generators require the improvement of low-cost processes for the fabrication of large assemblies of nanostructures, which also require an optimized doping concentration to achieve the maximum power factor. Metal-assisted chemical etching is a very advantageous technique, because it allows for the fabrication of silicon nanowires with high aspect ratio on a very large scale. However, MACE is incompatible with the high doping concentrations necessary for the optimization of the power factor. Our work is devoted to the investigation of the Seebeck coefficient of large-area silicon nanowire forests, doped by thermal diffusion after their fabrication. At first, we presented a reliable measurement procedure for the measurement of the Seebeck coefficient of macroscopic samples, made of large collections of nanometric structures (vertical silicon nanowires). Then, we demonstrated that the measured Seebeck coefficient is compatible with simulations of the diffusion process, which results in a nonuniform radial doping concentration.

Hence, the fabrication by MACE of slightly doped SiNW forests, and their doping by diffusion after fabrication, is a possible route for the exploitation of nanostructured silicon for thermoelectric purposes. Procedures for the fabrication of macroscopic nanostructured-Si generators, based on interconnected p- and n-doped legs, are under development.


  1. Bux, S. K.; Blair, R. G.; Gogna, P. K.; Lee, H.; Chen, G.; Dresselhaus, M. S.; Kaner, R. B.; Fleurial, J.-P. Adv. Funct. Mater. 2009, 19, 2445–2452. doi:10.1002/adfm.200900250
    Return to citation in text: [1]
  2. Bennett, N. S. Phys. Status Solidi A 2017, 214, 1700307. doi:10.1002/pssa.201700307
    Return to citation in text: [1] [2] [3]
  3. Dimaggio, E.; Pennelli, G. Nanotechnology 2018, 29, 135401. doi:10.1088/1361-6528/aaa9a2
    Return to citation in text: [1] [2] [3] [4]
  4. Gadea Díez, G.; Sojo Gordillo, J. M.; Pacios Pujadó, M.; Salleras, M.; Fonseca, L.; Morata, A.; Tarancón Rubio, A. Nano Energy 2020, 67, 104191. doi:10.1016/j.nanoen.2019.104191
    Return to citation in text: [1]
  5. Melosh, N. A.; Boukay, A.; Diana, F.; Gerardot, B.; Badolato, A.; Petroff, P. M.; Heath, J. R. Science 2003, 300, 112–115. doi:10.1126/science.1081940
    Return to citation in text: [1]
  6. Hochbaum, A. I.; Chen, R.; Delgado, R. D.; Liang, W.; Garnett, E. C.; Najarian, M.; Majumdar, A.; Yang, P. Nature 2008, 451, 163–167. doi:10.1038/nature06381
    Return to citation in text: [1]
  7. Feser, J. P.; Sadhu, J. S.; Azeredo, B. P.; Hsu, K. H.; Ma, J.; Kim, J.; Seong, M.; Fang, N. X.; Li, X.; Ferreira, P. M.; Sinha, S.; Cahill, D. G. J. Appl. Phys. 2012, 112, 114306. doi:10.1063/1.4767456
    Return to citation in text: [1]
  8. Pennelli, G.; Nannini, A.; Macucci, M. J. Appl. Phys. 2014, 115, 084507. doi:10.1063/1.4866994
    Return to citation in text: [1] [2]
  9. Lim, J.; Hippalgaonkar, K.; Andrews, S. C.; Majumdar, A.; Yang, P. Nano Lett. 2012, 12, 2475–2482. doi:10.1021/nl3005868
    Return to citation in text: [1]
  10. Pennelli, G.; Dimaggio, E.; Macucci, M. IEEE Trans. Nanotechnol. 2018, 17, 500–505. doi:10.1109/tnano.2018.2816119
    Return to citation in text: [1]
  11. Huang, Z.; Geyer, N.; Werner, P.; de Boor, J.; Gösele, U. Adv. Mater. (Weinheim, Ger.) 2011, 23, 285–308. doi:10.1002/adma.201001784
    Return to citation in text: [1]
  12. Kim, J.; Han, H.; Kim, Y. H.; Choi, S.-H.; Kim, J.-C.; Lee, W. ACS Nano 2011, 5, 3222–3229. doi:10.1021/nn2003458
    Return to citation in text: [1]
  13. Peng, K.; Yan, Y.; Gao, S.; Zhu, J. Adv. Funct. Mater. 2003, 13, 127–132. doi:10.1002/adfm.200390018
    Return to citation in text: [1]
  14. To, W.-K.; Tsang, C.-H.; Li, H.-H.; Huang, Z. Nano Lett. 2011, 11, 5252–5258. doi:10.1021/nl202674t
    Return to citation in text: [1]
  15. Curtin, B. M.; Fang, E. W.; Bowers, J. E. J. Electron. Mater. 2012, 41, 887–894. doi:10.1007/s11664-012-1904-1
    Return to citation in text: [1]
  16. Xu, B.; Fobelets, K. J. Appl. Phys. 2014, 115, 214306. doi:10.1063/1.4881781
    Return to citation in text: [1]
  17. Elyamny, S.; Dimaggio, E.; Magagna, S.; Narducci, D.; Pennelli, G. Nano Lett. 2020, 20, 4748–4753. doi:10.1021/acs.nanolett.0c00227
    Return to citation in text: [1] [2]
  18. Dimaggio, E.; Pennelli, G. Nano Lett. 2016, 16, 4348–4354. doi:10.1021/acs.nanolett.6b01440
    Return to citation in text: [1] [2]
  19. Qi, Y.; Wang, Z.; Zhang, M.; Yang, F.; Wang, X. J. Phys. Chem. C 2013, 117, 25090–25096. doi:10.1021/jp407720e
    Return to citation in text: [1]
  20. Dimaggio, E.; Narducci, D.; Pennelli, G. J. Mater. Eng. Perform. 2018, 27, 6279–6285. doi:10.1007/s11665-018-3605-2
    Return to citation in text: [1]
  21. Krali, E.; Durrani, Z. A. K. Appl. Phys. Lett. 2013, 102, 143102. doi:10.1063/1.4800778
    Return to citation in text: [1]
  22. Magagna, S.; Narducci, D.; Alfonso, C.; Dimaggio, E.; Pennelli, G.; Charaï, A. Nanotechnology 2020, 31, 404002. doi:10.1088/1361-6528/ab9b47
    Return to citation in text: [1]
  23. Pennelli, G.; Elyamny, S.; Dimaggio, E. Nanotechnology 2018, 29, 505402. doi:10.1088/1361-6528/aae5b6
    Return to citation in text: [1] [2]
  24. Fisher, D. Diffusion in Silicon; Scitec Publications, 1998.
    Return to citation in text: [1] [2]
  25. Geballe, T. H.; Hull, G. W. Phys. Rev. 1955, 98, 940–947. doi:10.1103/physrev.98.940
    Return to citation in text: [1]
  26. Brinson, M. E.; Dunstant, W. J. Phys. C: Solid State Phys. 1970, 3, 483–491. doi:10.1088/0022-3719/3/3/001
    Return to citation in text: [1]
  27. Stranz, A.; Kähler, J.; Merzsch, S.; Waag, A.; Peiner, E. Microsyst. Technol. 2012, 18, 857–862. doi:10.1007/s00542-011-1390-z
    Return to citation in text: [1]
  28. Arora, N. D.; Hauser, J. R.; Roulston, D. J. IEEE Trans. Electron Devices 1982, 29, 292–295. doi:10.1109/t-ed.1982.20698
    Return to citation in text: [1]

© 2020 Elyamny et al.; licensee Beilstein-Institut.
This is an Open Access article under the terms of the Creative Commons Attribution License ( Please note that the reuse, redistribution and reproduction in particular requires that the authors and source are credited.
The license is subject to the Beilstein Journal of Nanotechnology terms and conditions: (

Back to Article List

Other Beilstein-Institut Open Science Activities

Keep Informed

RSS Feed

Subscribe to our Latest Articles RSS Feed.


Follow the Beilstein-Institut


Twitter: @BeilsteinInst