Open Access Nano Commentary

Generation of high photocurrent in three-dimensional silicon quantum dot superlattice fabricated by combining bio-template and neutral beam etching for quantum dot solar cells

Makoto Igarashi12, Weiguo Hu12, Mohammad Maksudur Rahman12, Noritaka Usami3 and Seiji Samukawa124*

Author Affiliations

1 Institute of Fluid Science, Tohoku University, 2-1-1 Katahira, Aoba, Sendai 9808577, Japan

2 Japan Science and Technology Agency, CREST, 5 Sanbancho, Chiyoda, Tokyo 1020075, Japan

3 Institute for Materials Research, Tohoku University, 2-1-1 Katahira, Aoba, Sendai, 9808577, Japan

4 WPI Advanced Institute for Materials Research, Tohoku University, 2-1-1 Katahira, Aoba, Sendai, 9808577, Japan

For all author emails, please log on.

Nanoscale Research Letters 2013, 8:228  doi:10.1186/1556-276X-8-228

The electronic version of this article is the complete one and can be found online at:

Received:7 March 2013
Accepted:1 May 2013
Published:15 May 2013

© 2013 Igarashi et al.; licensee Springer.

This is an Open Access article distributed under the terms of the Creative Commons Attribution License (, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.


We fabricated a three-dimensional (3D) stacked Si nanodisk (Si-ND) array with a high aspect ratio and uniform size by using our advanced top-down technology consisting of bio-template and neutral beam etching processes. We found from conductive atomic microscope measurements that conductivity became higher as the arrangement was changed from a single Si-ND to two-dimensional (2D) and 3D arrays with the same matrix of SiC, i.e., the coupling of wave functions was changed. Moreover, our theoretical calculations suggested that the formation of minibands enhanced tunneling current, which well supported our experimental results. Further analysis indicated that four or more Si-NDs basically maximized the advantage of minibands in our structure. However, it appeared that differences in miniband widths between 2D and 3D Si-ND arrays did not affect the enhancement of the optical absorption coefficient. Hence, high photocurrent could be observed in our Si-ND array with high photoabsorption and carrier conductivity due to the formation of 3D minibands.

Si nanodisk; Aspect ratio; Photocurrent; Miniband


Quantum dot superlattices (QDSLs) have attracted a great deal of interest from both physical scientists and device researchers. Electron wave functions diffuse and overlap, which merge discrete quantum levels into minibands, with quantum dots approaching and forming a quasi-crystal structure. This band rearrangement has significant applications for many novel optoelectronic/electronic devices [1-15]. For example, quantum dot solar cells, the most exciting photovoltaic device with more than 63% conversion efficiency, have to utilize minibands for carrier transport and additional optical transitions.

Ideal QDSLs present a great challenge to current nanotechnologies. Several technologies (e.g., chemical solution methods and molecular beam epitaxy (MBE)) have convincingly been used to fabricate relatively uniform quantum dots; however, very few technologies can finitely arrange QDs to form a quasi-crystal structure. The well-developed MBE technology can only achieve very limited control on the direction of growth, which induces a mixed state with the wetting layer. The most direct idea is to develop a top-down nanotechnology. However, nanometer-order sizes exceed most light/electron beam limitations, and suitable masks seem impossible to create. The neutral beam (NB) etching and ferritin bio-template we developed have recently brought about a great breakthrough in that we successfully fabricated two-dimensional (2D) array Si nanodisks (Si-NDs) with sub-10 nm, high density (>1011 cm-2), and quasi-hexagonal crystallization [16-20].

Photovoltaic conversion efficiency was determined by light absorbance and carrier collection efficiency. Our previous work has proven that wave function coupling relaxes the selection rule to induce additional optical transitions [21,22]. We first observed enhanced conductivity in 2D and three-dimensional (3D) array Si-NDs with a SiC matrix in this study. Moreover, we calculated electronic structures and current transport, which theoretically suggested that minibands enhanced conductivity, within envelope function theory and the Anderson Hamiltonian method. These enhanced optical and electrical properties indicated a potential application for the highly efficient quantum dot solar cells.


The fabrication of the 3D Si-ND array was based on bio-template and NB processes. Figure 1 schematically illustrates the fabrication flow, which started with (Figure 1a) a 2-nm-thick SiC film and 4-nm-thick poly-Si being deposited alternately four times on the n-doped Si substrate using a high-vacuum sputtering system and electron beam evaporation. Then a 3-nm-thick SiO2 layer was fabricated as a surface oxide (called NBO-SiO2 after this) by the NB oxidation process we developed at a low temperature of 300°C [16]. Figure 1b has a 2D array of bio-template molecules (Listeria-Dps) that was deposited on the surface of the NBO-SiO2. Figure 1c shows the bio-template protein shell that was removed by annealing it in an oxygen atmosphere to obtain a 2D array of iron cores as a uniform mask for the etching process. Figure 1d shows the etching process that was carried out with nitrogen trifluoride gas/hydrogen radical treatment (NF3 treatment) to remove the surface SiO2, which was carried out with NB etching to remove the poly-Si. Here we performed a one-step etching and found a well-aligned vertical etching profile due to high etching selectivity between the iron cores and etched material and the low selectivity of 1.3 between Si and SiC. The etching process has been detailed elsewhere [17-19]. Figure 1e shows that the iron cores were then removed by HCl wet cleaning, and then the remaining surface SiO2 was removed by NF3 treatment. Figure 1f shows that the SiC was deposited between pillars, which were stacked Si-NDs, by the sputtering system. The diameter, space between NDs, and average ND center-to-ND center distance corresponded to 6.4, 2.3, and 8.7 nm in the structure. The size distribution of the Si-NDs was less than 10% for all samples [19,21]. We prepared three types of Si-ND arrangements, as seen in Figure 2: separated Si-NDs as a single QD, a 2D array of Si-NDs as a 2D QDSL, and a 3D array of Si-NDs as a 3D QDSL. The electrical conductivity and optical absorption in QDSLs were methodically, experimentally, and theoretically investigated with these samples to study the effect of wave function coupling between QDs.

thumbnailFigure 1. Schematic of the fabrication flow for 3D array of Si-NDs with SiC interlayer. (a) Deposition of 2-nm-thick SiC, 4-nm-thick poly-Si, and 3-nm-thick SiO2 layers. (b) Arrangement of 2D array of bio-template molecules on the surface. (c) Removal of bio-template protein shell by annealing in oxygen atmosphere. (d) NF3 treatment to remove surface SiO2 and NB etching to remove surface multilayers of poly-Si and SiC. (e) Removal of iron cores with HCl and NF3 treatment to etch remaining surface SiO2. (f) SiC deposition on Si-NDs.

thumbnailFigure 2. Schematics of the three types of Si-ND arrangements. (a) Separated Si-NDs as single QD, (b) 2D array of Si-NDs as 2D QDSL, and (c) 3D array of Si-NDs as 3D QDSL.

Results and discussion

Conductive atomic force microscopy (c-AFM) has been used to investigate conductivity, as seen in Figure 3. Changing the matrix from SiO2 to SiC greatly increases current (I) and decreases threshold voltage (V), according to comparisons of the 2D arrays of Si-NDs. Although a primary factor should be macro-conductivity differences between SiC and SiO2, one cause is minibands that enhance conductivity, which was revealed in a later theoretical simulation. More significantly, conductivity became higher as the arrangement was changed from a single Si-ND to 2D and 3D arrays with the same matrix of SiC, i.e., the coupling of wave functions was changed. Note that conductivity in the 3D array was higher than that in the 2D array, even though the total thickness of the QDSL expanded. These results indicate that the formation of minibands both in-plane and out-of-plane (vertically) might enhance carrier conductivity in QDSLs.

thumbnailFigure 3. I-V curves of single Si-ND, 2D, and 3D arrays of Si-NDs measured by c-AFM. Red, blue, and green lines plot results for the 3D array, 2D array, and single Si-ND with SiC matrix. Black line plots the results for 2D array Si-NDs with SiO2 matrix.

We considered resonant tunneling to be a theoretical mechanism that could explain our experimental results on the basis of these results. Therefore, we theoretically investigated enhanced conductivity due to the formation of minibands. Our developed top-down nanotechnology achieved great flexibility in designing parts for the quantum structure, such as the independently controllable diameter and thickness, high aspect ratio, and different matrix materials. The finite element method duly described the complex quantum structures. The electronic structure and wave function within envelope function theory are presented as.

<a onClick="popup('','MathML',630,470);return false;" target="_blank" href="">View MathML</a>


Here we mainly took into consideration the matrix material, realistic geometry structure, and number of stacking layers. The results are presented in Figure 4. A distinct feature is that electron wave functions are more strongly confined in the Si-NDs in the SiO2 matrix due to the higher band offset of the Si/SiO2 interface. Thus, they resulted in higher quantum levels. In addition, stronger confinement means weaker coupling of the wave function and narrower minibands in the same geometry alignment. By stacking our NDs from one layer to ten layers, the miniband in Figure 5 gradually broadens, and at around four to six layers, the miniband width seems to saturate. The probability of the wave function diffusing into the barrier exponentially reduces with distance, which indicates that wave function coupling exponentially saturates as the number of layers increases. Perhaps four- or six-layer NDs are sufficient to maximize the advantage of minibands.

thumbnailFigure 4. Calculated results for electron spatial possibilities. In three lateral coupled NDs and miniband width in 2D array of Si-NDs. Square wave functions and quantum levels of coupled NDs with (a) SiC matrix and (b) SiO2 matrix.

thumbnailFigure 5. Calculated results for miniband width in 3D array of Si-NDs. Thickness, diameter, and space between NDs were assumed to correspond to 4.0, 6.4 and 2.0 nm.

Chang et al. [23] considered interdot coupling with the Anderson Hamiltonian model to deduce tunneling current density as

<a onClick="popup('','MathML',630,470);return false;" target="_blank" href="">View MathML</a>


Here E(kxy) is related to the energy discrepancy, t, due to in-plane ND coupling E(kxy) = 2t[cos(kxR) + cos(kyR)]. We simulated the I-V properties of our structures with this. The results are in Figure 6. The calculated results also revealed that the wider minibands in the SiC matrix resulted in better transport properties than those in the SiO2 matrix. A simplified, but not too obscure, explanation is that the formation of minibands broadens the resonance levels to increase joint-state density. Carrier transport in this two-barrier structure mainly depends on resonant tunneling. Moreover, if the Coulomb blockade effect is neglected, the tunneling joint-state density in Equation 2 can be simplified as a parabola function with a resonant peak at ~E0– E(kxy). The formation of minibands broadens the resonant peak to allow more states to approach maximum, which results in enhanced current. Thus, wider minibands mean a higher current density and lower threshold voltage, as can be seen in the Si-NDs in the SiC matrix. In addition, the 2D array of Si-NDs in the SiC matrix has a lower miniband level, E0, which also shifts the I-V curves to a lower threshold voltage. This tendency closely matches that in our experimental results, and due to the larger tunneling resistance in the SiO2 interlayer (Ct), the threshold voltage (V) is further increased in realistic I-V curves. Moreover, conductivity in the 2D and 3D arrays of Si-NDs was enhanced due to the same mechanism that broadened the wave functions and formed wider minibands. As these were also very consistent with the trend in our experimental results, they clarified that the formation of minibands both in-plane and out-of-plane could enhance carrier transport in QDSLs. Enhanced conductivity is very important for electronic/optoelectronic devices, which indicates high charge injection efficiency in lasers and carrier collection efficiency in solar cells.

thumbnailFigure 6. Simulation results for I-V properties of our sample structures. Red, blue, and green lines plot calculated results for 3D array, 2D array, and single Si-ND with SiC matrix. Black line plots results for 2D array Si-NDs with SiO2 matrix.

Optical absorption was then investigated by measuring the transmittance of samples using ultraviolet-visible-near-infrared spectroscopy. Our previous work demonstrated that the formation of minibands perpendicular to incident light could enhance photon absorption, i.e., 2D minibands could improve the absorption coefficient in the 2D array of Si-NDs [21,22]. Therefore, we investigated what effect 3D minibands had on optical absorption in this study. Figure 7 shows the absorption coefficients in the 2D and 3D Si-ND array samples prepared on transparent quartz substrates. The absorption coefficient in the 3D array was almost the same as that in the 2D array, and the calculated bandgap energy of both samples was 2.2 eV. Moreover, the change in the miniband width between the samples should be 3.85 meV, as shown in Figure 5 (0.95 meV in single layer and 4.80 meV in four layers). Therefore, it seems that the change of 3.85 meV in the miniband width is not sufficiently large to affect photon absorption.

thumbnailFigure 7. Absorption coefficients of 2D and 3D arrays of Si-NDs with SiC matrix. Blue and red lines correspond to 2D and 3D arrays of Si-NDs.

Finally, we fabricated a p++-i-n Si solar cell with a 3D array of Si-NDs as an absorption layer, as shown in Figure 8, and measured the amount of possible photocurrent generated from the Si-ND layers where the high doping density (>1020 cm-3) of the p++-Si substrate prevented photocurrent from being generated inside the substrate itself. Here we found that the generated short-circuit current density from the p++-i-n solar cell was 2 mA/cm2, where the largest possible photocurrent generated in the Si-ND layers and n-Si emitter was estimated to be 3.5 mA/cm2 for the former and 1.0 mA/cm2 for the latter [22]. Since 1 mA/cm2 is the highest possible value for photocurrent from the n-Si emitter according to this estimate, the actual value should be lower than the calculated value. Therefore, we found that out of the total photocurrent of 2 mA/cm2, much more of it (>1 mA/cm2) was contributed to by Si-ND. This confirms that most of the observed photocurrent originated from the carrier generated at the Si-ND itself because of high photoabsorption and carrier conductivity due to the formation of 3D minibands in our Si-ND array.

thumbnailFigure 8. I-V characteristics of p++-i-n solar cell. Current-voltage characteristics in dark (blue line) and under sunlight (red line).


We developed an advanced top-down technology to fabricate a stacked Si-ND array that had a high aspect ratio and was of uniform size. We found from c-AFM measurements that conductivity increased as the arrangement was changed from a single Si-ND to 2D and 3D arrays with the same matrix of SiC. This enhancement was most likely due to the formation of minibands, as suggested by our theoretical calculations. Moreover, the change in out-of-plane minibands did not affect the absorption coefficient. This enhanced transport should work in the collection efficiency of high carriers in solar cells.


c-AFM: Conductive atomic force microscopy; I-V: Current-voltage; MBE: Molecular beam epitaxy; ND: Nanodisk; QDSL: Quantum dot superlattices.

Competing interests

The authors declare that they have no competing interests.

Authors' contributions

MI and SS conceived and designed the experiment, fabricated the silicon nanodisk samples, performed electrical and optical measurements, analyzed these data, and wrote the paper. MMR and NU fabricated the solar cell structures and analyzed the I-V data. WH performed the theoretical calculations. All authors discussed the results, commented on the manuscript, and read and approved the final version.


This work is supported by the Japan Science and Technology Agency (JST CREST) and the Grant-in-Aid for Japan Society for the Promotion of Science (JSPS) Fellows.


  1. Luque A, Marti A: Increasing the efficiency of ideal solar cells by photon induced transitions at intermediate levels.

    Phys Rev Lett 1997, 78:5014. Publisher Full Text OpenURL

  2. Konakov AA, Burdov VA: Optical gap of silicon crystallites embedded in various wide-band amorphous matrices: role of environment.

    J Phys Condens Matter 2010, 22:215301. PubMed Abstract | Publisher Full Text OpenURL

  3. Thamdrup LH, Persson F, Bruus H, Kristensen A, Flyvbjerg H: Experimental investigation of bubble formation during capillary filling of SiO2 nanoslits.

    Appl Phys Lett 2007, 91:163505. Publisher Full Text OpenURL

  4. Conibeer G, Green MA, Corkish R, Cho Y, Cho EC, Jiang CW, Fangsuwannarak T, Pink E, Huang Y, Puzzer T, Trupke T, Richards B, Shalav A, Lin KL: Silicon nanostructures for third generation photovoltaic solar cells.

    Thin Solid Films 2006, 511–512:654. OpenURL

  5. Cho EC, Park S, Hao X, Song D, Conibeer G, Park SC, Green MA: Silicon quantum dot/crystalline silicon solar cells.

    Nanotechnol 2008, 19:245201. Publisher Full Text OpenURL

  6. Conibeer G, Green MA, Cho EC, König D, Cho YH, Fangsuwannarak T, Scardera G, Pink E, Huang Y, Puzzer T, Huang S, Song D, Flynn C, Park S, Hao X, Mansfield D: Silicon quantum dot nanostructures for tandem photovoltaic cells.

    Thin Solid Films 2008, 516:6748. Publisher Full Text OpenURL

  7. Nuryadi R, Ikeda H, Ishikawa Y, Tabe M: Ambipolar Coulomb blockade characteristics in a two-dimensional Si multidot device.

    IEEE Trans Nanotechnol 2003, 2:231. Publisher Full Text OpenURL

  8. Cordan AS, Leroy Y, Goltzene A, Pepin A, View C, Mejias M, Launois H: Temperature behavior of multiple tunnel junction devices based on disordered dot arrays.

    J Appl Phys 2000, 87:345. Publisher Full Text OpenURL

  9. Uchida K, Koga J, Ohba R, Takagi SI, Toriumi A: Silicon single-electron tunneling device fabricated in an undulated ultrathin silicon-on-insulator film.

    J Appl Phys 2001, 90:3551. Publisher Full Text OpenURL

  10. Macucci M, Gattobigio M, Bonci L, Iannaccone G, Prins FE, Single C, Wetekam G, Kern DP: A QCA cell in silicon-on-insulator technology: theory and experiment.

    Superlattices Microstruct 2003, 34:205. Publisher Full Text OpenURL

  11. Lent CS, Tougaw PD: A device architecture for computing with quantum dots.

    Proc IEEE 1997, 85:541. Publisher Full Text OpenURL

  12. Nassiopoulou AG, Olzierski A, Tsoi E, Berbezier I, Karmous A: Ge quantum dot memory structure with laterally ordered highly dense arrays of Ge dots.

    J Nanosci Nanotechnol 2007, 7:316. PubMed Abstract OpenURL

  13. Pothier H, Lafarge P, Urbina C, Esteve D, Devoret MH: Single-electron pump based on charging effects.

    Europhys Lett 1992, 17:249. Publisher Full Text OpenURL

  14. Shin M, Lee S, Park KW: The study of a single-electron memory cell using coupled multiple tunnel-junction arrays.

    Nanotechnol 2001, 12:178. Publisher Full Text OpenURL

  15. Hirvi KP, Paalanen MA, Pekola JP: Numerical investigation of one‒dimensional tunnel junction arrays at temperatures above the Coulomb blockade regime.

    J Appl Phys 1996, 80:256. Publisher Full Text OpenURL

  16. Igarashi M, Tsukamoto R, Huang CH, Yamashita I, Samukawa S: Direct fabrication of uniform and high density sub-10-nm etching mask using ferritin molecules on Si and GaAs surface for actual quantum-dot superlattice.

    Appl Phys Express 2011, 4:015202. Publisher Full Text OpenURL

  17. Huang CH, Igarashi M, Woné M, Uraoka Y, Fuyuki T, Takeguchi M, Yamashita I, Samukawa S: Two-dimensional Si-nanodisk array fabricated using bio-nano-process and neutral beam etching for realistic quantum effect devices.

    Jpn J Appl Phys 2009, 48:04C187. Publisher Full Text OpenURL

  18. Huang CH, Igarashi M, Horita S, Takeguchi M, Uraoka Y, Fuyuki T, Yamashita I, Samukawa S: Novel Si nanodisk fabricated by biotemplate and defect-free neutral beam etching for solar cell application.

    Jpn J Appl Phys 2010, 49:04DL16. Publisher Full Text OpenURL

  19. Huang CH, Wang XY, Igarashi M, Murayama A, Okada Y, Yamashita I, Samukawa S: Optical absorption characteristic of highly ordered and dense two-dimensional array of silicon nanodiscs.

    Nanotechnol 2011, 22:105301. Publisher Full Text OpenURL

  20. Hirano R, Miyamoto S, Yonemoto M, Samukawa S, Sawano K, Shiraki Y, Itoh KM: Room-temperature observation of size effects in photoluminescence of Si0.8Ge0.2/Si nanocolumns prepared by neutral beam etching.

    Appl Phys Express 2012, 5:082004. Publisher Full Text OpenURL

  21. Budiman MF, Hu W, Igarashi M, Tsukamoto R, Isoda T, Itoh KM, Yamashita I, Murayama A, Okada Y, Samukawa S: Control of optical bandgap energy and optical absorption coefficient by geometric parameters in sub-10 nm silicon-nanodisc array structure.

    Nanotechnol 2012, 23:065302. Publisher Full Text OpenURL

  22. Igarashi M, Budiman MF, Pan W, Hu W, Tamura Y, Syazwan ME, Usami N, Samukawa S: Effects of formation of mini-bands in two-dimensional array of silicon nanodisks with SiC interlayer for quantum dot solar cells.

    Nanotechnol 2013, 24:015301. Publisher Full Text OpenURL

  23. Kuo DMT, Guo GY, Chang YC: Tunneling current through a quantum dot array.

    Appl Phys Lett 2001, 79:3851. Publisher Full Text OpenURL