Home | About | Journals | Submit | Contact Us | Français |

**|**Nanoscale Res Lett**|**v.7(1); 2012**|**PMC3337235

Formats

Article sections

- Abstract
- 1. Introduction
- 2. The SWS and optical model
- 3. Results and discussion
- 4. Conclusions
- Abbreviations
- Competing interests
- Authors' contributions
- References

Authors

Related links

Nanoscale Res Lett. 2012; 7(1): 196.

Published online 2012 March 23. doi: 10.1186/1556-276X-7-196

PMCID: PMC3337235

Yiming Li: ymli/at/faculty.nctu.edu.tw; Ming-Yi Lee: mylee/at/mail.ymlab.org; Hui-Wen Cheng: hwcheng/at/mail.ymlab.org; Zheng-Liang Lu: zllu/at/mail.ymlab.org

Received 2011 November 17; Accepted 2012 March 23.

Copyright ©2012 Li et al; licensee Springer.

This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

In this study, we investigate the reflectance property of the cylinder, right circular cone, and square pyramid shapes of silicon nitride (Si_{3}N_{4}) subwavelength structure (SWS) with respect to different designing parameters. In terms of three critical factors, the reflectance for physical characteristics of wavelength dependence, the reflected power density for real power reflection applied on solar cell, and the normalized reflectance (reflected power density/incident power density) for real reflectance applied on solar cell, a full three-dimensional finite element simulation is performed and discussed for the aforementioned three morphologies. The result of this study shows that the pyramid shape of SWS possesses the best reflectance property in the optical region from 400 to 1000 nm which is useful for silicon solar cell applications.

Silicon solar cell is one of the promising renewable energy technologies in order to relieve the impact of the climate change. In semiconductor-based solar cells, electron-hole pairs are generated through absorption of impinging photons. Due to high refraction index of semiconductor materials, especially silicon, the incident sunlight power is largely reflected back, resulting in the reduction of light absorption and poor energy conversion efficiency. Antireflection coating (ARC) is mounted over absorption layers, resulting in three effects: (a) reduction in surface reflection, (b) increase in light absorption due to an increase in optical path length by diffraction, and (c) enhancement of internal reflection that reduces the amount of escaping light. Based on the theory of impedance matching, single layer (SLR) and multilayer of ARC are proposed for reduced reflectance property; however, the resulting reflectance spectra meet the demand only within a narrow spectral domain. Subwavelength structure's (SWS) dimensions are much smaller than the wavelengths of light; therefore, using ARC on the surface of silicon solar cells can substantially reduce the reflectivity and improve the capability of light trapping. It thus will achieve the enhanced efficiency according to our recent both numerical and experimental studies [1-3]. Compared with silicon solar cell with a SLAR, the efficiency of silicon solar cell with Si_{3}N_{4 }SWS is promising among various ARC layers in our recent work [4]. A rigorous coupled-wave analysis (RCWA) [1,5-7] has been reported to estimate the reflectance of Si_{3}N_{4 }SWS by approximating structural shapes with partitioned uniform homogeneous layers. RCWA is an exact solution of Maxwell's equations for the electromagnetic diffraction by grating structures which is generally applicable for 2D plane with 1D periodicity; however, RCWA may suffer numerical difficulties in presence of evanescent orders and it requires a large amount of calculation for retaining several diffraction orders. These factors limit flexible application of RCWA; in particular, for 3D problems with non-azimuthally symmetric structural shapes. Numerical simulation of 3D morphological effect on reflectance property has not been studied yet. Therefore, a full 3D finite-element (FE) analysis of Si_{3}N_{4 }SWS will be an interesting examination for quantitative understanding of the reflectance property.

In this study, 3D FE simulation for the reflectance of Si_{3}N_{4 }SWS with three types of structural shapes, the cylinder, the right circular cone, and the square pyramid shapes, is conducted with respect to different geometry parameters and lighting angles for quantitative understanding of reflectance property. First, proper selection on the boundary conditions can alleviate the computational load from simulating a holistic ARC. The reflectance of Si_{3}N_{4 }SWS on the silicon substrate is thus simulated using the 3D finite element method (FEM); consequently, in terms of three critical factors, the reflectance for physical characteristics of wavelength dependence, the reflected power density for real power reflection applied on solar cell, and the normalized reflectance (reflected power density/incident power density) for real reflectance applied on solar cell are calculated and discussed for the aforementioned three morphologies. The analysis of reflectance spectrum with wide-angle incidences of electromagnetic wave and the average reflectance with various heights are presented. Besides, according to our recent study, which presented the optimal design parameters of Si_{3}N_{4 }SWS based on RCWA [4], numerical verification and comparison is accomplished following the discussion. The engineering findings of this study show that the pyramid shape of SWS possesses the best reflectance property in the optical region from 400 to 1000 nm which is useful for silicon solar cell applications.

This rest of the article is organized as follows. In Section 2, we show the computational structure and model. In Section 3, we report the results and discussion. Finally, we draw conclusions and suggest future work.

Based upon our experimental characterization, Figure Figure1a1a illustrates a periodical structure of Si_{3}N_{4 }SWS which is used in our 3D FE simulation without loss of generality. We study Si_{3}N_{4 }SWS with the cylinder, the right circular cone, and the square pyramid shapes, as shown in Figure 1b-d, respectively. With a constant volume, the diameter of cylinder- and right circular cone-shaped Si_{3}N_{4 }SWS and the edge length of square pyramid are 130 nm, the heights (h) of the etched part of Si_{3}N_{4 }SWS are 200, 600, and 471.3 nm, the height (s) of the non-etched part is 70 nm, and the base (W) of a unit cell is 200 nm [4]. Note that the thickness of Si substrate is given 600 nm. Note that all structural parameters are adopted from our experimental studies [2-4,8]. Throughout the article, we consider time-harmonic fields assuming a time-dependence in *e ^{-jωt}*. The diffraction problem is governed by the well-known Maxwell equations

(1)

(2)

(3)

and

(4)

where *E *and *D *are electric field intensity and flux density, *H *and *B *are magnetic field intensity and flux density, *λ *is the corresponding frequency to the wavelength *λ, J*, and ρ are current density and charge density, *ε *is electric permittivity, *μ *is magnetic permeability. A repeated pattern is applicable to use periodic boundary conditions, thus the Floquet theorem is adopted to simulate the boundary condition of periodic structure. Floquet theorem asserts that the analysis region can be reduced significantly in one periodicity cell to characterize the propagation property. The electric fields in periodic structure are related as follows:

(5)

where *r *is position vector, *L *is the distance between the periodic boundaries, and *θ *is a phase factor determined by wave vector *k *and *L*:

(6)

The polarization of transverse electric (TE) mode, in which the electric field is normal to the direction of wave propagation, is excited as the normal incident light source with wavelengths sweeping from 400 to 1000 nm. The bottom region of Si substrate is assigned as perfect matched layer in avoidance of reflected wave. The refraction index of Si_{3}N_{4 }is 2.05, and the refraction index of Si is frequency dependent with the relation [1]:

(7)

where *λ *is the incident wavelength, *A *= 0.939816, *B *= 8.10461 × 10^{-3}, λ_{1 }= 1.1071 μm, and ε = 11.6858. The calculation settings of reflectance were reported and can be found in our recent studies [1,4].

In order to examine the effect of Floquet boundary condition in 3D FE analysis, as shown in Figure Figure2,2, we compare the difference between the simulated unit cells of 1 × 1 and 2 × 2 array of Si_{3}N_{4 }SWS. We find at the wavelengths above 600 nm, the reflectance of 1 × 1 array of Si_{3}N_{4 }SWS as unit cell is almost consistent with unit cell of 2 × 2 array, meanwhile insignificant discrepancy occurs at wavelengths shorter than 600 nm. Based on this consequence, it is enable us to do simulation with more computational efficient using 1 × 1 array of Si_{3}N_{4 }SWS as a simulated unit cell with engineering acceptable accuracy. According to our recent RCWA work [4], the reflectance spectra are first plotted in Figure Figure33 using the optimal design parameters [1,4]. Also, the spectra calculated by a full 3D FE analysis with the same design parameters are indicated by dashed lines. For the cylinder-shaped Si_{3}N_{4 }SWS, the reflectance spectra for RCWA and FE analysis are similar, but not agreed for the cone-shaped Si_{3}N_{4 }SWS due to existing evanescent orders along the top of structures. This comparison confirms the importance of 3D FEM simulation which is beyond the RCWA approach [4].

Figure 4a-c shows the reflectance spectra with incident angles of 0°, 15°, 30°, 45°, and 60° for the cylinder-, right circular cone-, and square pyramid-shaped Si_{3}N_{4 }SWS, respectively. For the normal incidence case, the lowest average reflectance among three structural shapes is 3.47% of square pyramid-shaped structure. The others are 6.86 and 4.42% for the cylinder- and the right circular cone-shaped Si_{3}N_{4 }SWS, respectively. Meanwhile, as shown in Figure 4a-c, one can observe that the reflectance increases significantly with larger incident angles, resulting in average reflectance beyond 50%. Table Table11 summarizes the average reflectance for various incident angles. Height effect on average reflectance of Si_{3}N_{4 }SWS at normal incident angle with *d *= 130 nm and *s *= 70 nm is also calculated, as shown in Figure Figure5.5. The resulting average reflectance of pyramid-shaped Si_{3}N_{4 }SWS nearly keeps lowest in comparison with the cylinder- and the right circular cone-shaped Si_{3}N_{4 }SWS as the structural height is ranging from 50 to 500 nm. Figure Figure66 shows the reflectance dependence on the structural height and wavelength. The pyramid-shaped Si_{3}N_{4 }SWS has lower reflectance and less sensitivity on structure height in comparison with the cylinder-shaped Si_{3}N_{4 }SWS. Hence, the impact of process variation of structure height on solar cell performance is smaller for pyramid-shaped Si_{3}N_{4 }SWS. Based on solar spectrum at the sea level revealed in American Society for Testing and Materials (ASTM) Standard Tables for Reference Solar Spectral Irradiances: Direct Normal and Hemispherical 37 Tilted Surface [9], we further estimate the reflected power density (W/m^{2}/nm) defined by reflectance times incident power density, as shown in Figure Figure7.7. The higher reflected power density of cylinder-shaped Si_{3}N_{4 }SWS (red line) indicates the less efficiency in the solar cell application. Therefore, normalized reflectance defined as

(8)

reveals the real power efficiency applied in the solar cell application. Figure Figure88 shows the normalized reflectance for the cylinder-, right circular cone-, and square pyramid-shaped Si_{3}N_{4 }SWS, respectively. The square pyramid-shaped Si_{3}N_{4 }SWS again shows the lowest normalized reflectance 3.13% while the cylinder- and the right circular cone-shaped Si_{3}N_{4 }SWSs have 6.66 and 4.12%, respectively.

In this study, the reflective property of unit cell with a validated Floquet boundary condition has been calculated using a full 3D FE simulation. Considering various incidence angles and height effect on three experimentally observed structural shapes of Si_{3}N_{4 }SWS, we have concluded that the pyramid-shaped Si_{3}N_{4 }SWS has best reflective property in the analysis of morphological effect. Compared with the results of RCWA, the reflective property calculated by the full 3D FEM is significantly deviated from the results from RCWA, giving the hint that a detailed and comprehensive methodology is dispensable for the design of Si_{3}N_{4 }SWS. The results of computed reflectance, reflected power density, and normalized reflectance have shown that the pyramid shape of SWS may have the best reflectance property in the optical region from 400 to 1000 nm which is useful for silicon solar cell applications. The optimized pyramid-shaped Si_{3}N_{4 }SWS is currently under plan for implementation with silicon solar cells.

3D: three-dimensional; ARC: antireflection coating; FEM: finite element method; RCWA: rigorous coupled-wave analysis; Si_{3}N_{4}: silicon nitride; SLAR: single layer; SWS: subwavelength structure.

The authors declare that they have no competing interests.

M-YL, H-WC, and Z-LL performed the numerical simulation and data analysis, YL conducted whole study including manuscript preparation. All the authors read and approved the final manuscript.

This study was supported in part by the Taiwan National Science Council under contract Nos. NSC-99-2221-E-009-175 and NSC-100-2221-E-009-018.

- Sahoo KC, Li Y, Chang EY. Numerical calculation of reflectance of sub-wavelength structures on silicon nitride for solar cell application. Comput Phys Commun. 2009;180:1721–1729. doi: 10.1016/j.cpc.2009.04.013. [Cross Ref]
- Sahoo KC, Lin MK, Chang EY, Tinh TB, Li Y, Huang JH. Silicon nitride nanopillars and nanocones formed by nickel nanoclusters and inductively coupled plasma etching for solar cell application. Jpn J Appl Phys. 2009;48:126508. doi: 10.1143/JJAP.48.126508. [Cross Ref]
- Sahoo KC, Chang EY, Lin MK, Li Y, Huang JH. Fabrication and configuration development of silicon nitride sub-wavelength structures for solar cell application. J Nanosci Nanotechnol. 2010;10:5692–5699. doi: 10.1166/jnn.2010.2553. [PubMed] [Cross Ref]
- Sahoo KC, Li Y, Chang EY. Shape effect of silicon nitride sub-wavelength structure on reflectance for solar cell application. IEEE Trans Electron Dev. 2010;57:2427–2433.
- Moharam MG, Gaylord TK. Rigorous coupled-wave analysis of planar-grating diffraction. J Opt Soc Am. 1981;71:811–818. doi: 10.1364/JOSA.71.000811. [Cross Ref]
- Moharam MG, Gaylord TK. Rigorous coupled-wave analysis of metallic surface-relief gratings. J Opt Soc Am A. 1986;3:1780–1787. doi: 10.1364/JOSAA.3.001780. [Cross Ref]
- Moharam MG, Pommet DA. Formulation for stable and efficient implementation of the rigorous coupled-wave analysis of binary gratings. J Opt Soc Am A. 1995;12:1077–1086. doi: 10.1364/JOSAA.12.001077. [Cross Ref]
- Sahoo KC, Lin MK, Chang EY, Lu YY, Chen CC, Huang JH, Change CW. Fabrication of antireflective sub-wavelength structures on silicon nitride using nano cluster mask for solar cell application. Nanoscale Res Lett. 2009;4:680–683. doi: 10.1007/s11671-009-9297-7. [PMC free article] [PubMed] [Cross Ref]
- ASTM Standard G173-03. Standard Tables for Reference Solar Spectral Irradiances: Direct Normal and Hemispherical on 37° Tilted Surface. 2008. http://www.astm.org/Standards/G173.htm

Articles from Nanoscale Research Letters are provided here courtesy of **Springer**

PubMed Central Canada is a service of the Canadian Institutes of Health Research (CIHR) working in partnership with the National Research Council's Canada Institute for Scientific and Technical Information in cooperation with the National Center for Biotechnology Information at the U.S. National Library of Medicine(NCBI/NLM). It includes content provided to the PubMed Central International archive by participating publishers. |