Simulation of High Efficiency Heterojunction Solar Cells with AFORS-HET

Similar documents
Materials, Electronics and Renewable Energy

Thin film silicon technology. Cosimo Gerardi 3SUN R&D Tech. Coordinator

Light-Induced Degradation of Thin Film Silicon Solar Cells

N-PERT BACK JUNCTION SOLAR CELLS: AN OPTION FOR THE NEXT INDUSTRIAL TECHNOLOGY GENERATION?

Available online at ScienceDirect. Energy Procedia 84 (2015 ) 17 24

TWO-DIMENSIONAL MODELING OF EWT MULTICRYSTALLINE SILICON SOLAR CELLS AND COMPARISON WITH THE IBC SOLAR CELL

Thin film CdS/CdTe solar cells: Research perspectives

Study of a-sige:h Films and n-i-p Devices used in High Efficiency Triple Junction Solar Cells.

PROSPECTS OF MOLYBDENUM DISULFIDE (MoS 2 ) AS AN ALTERNATIVE ABSORBER LAYER MATERIAL IN THIN FILM SOLAR CELLS FROM NUMERICAL MODELING

Preparation and Characterization of Micro-Crystalline Hydrogenated Silicon Carbide p-layers

Amorphous Silicon Solar Cells

Available online at Energy Procedia 00 (2009) Energy Procedia 2 (2010) E-MRS Spring meeting 2009, Symposium B

EFFECT OF HYDROGEN, CERIUM AND TUNGSTEN DOPING ON INDIUM OXIDE THIN FILMS FOR HETEROJUNCTION SOLAR CELLS

24th European Photovoltaic Solar Energy Conference and Exhibition, September 2009, Hamburg, Germany.

Solar cells conversion efficiency enhancement techniques

xsi PV Technologies: Status and Outlook of Passivated Contacts and Rear Contacted Solar Cells

Photovoltaics under concentrated sunlight

INDUSTRIALLY FEASIBLE >19% EFFICIENCY IBC CELLS FOR PILOT LINE PROCESSING

HIGH EFFICIENCY INDUSTRIAL SCREEN PRINTED N-TYPE SOLAR CELLS WITH FRONT BORON EMITTER

Schottky-barrier and MIS solar cells

Silver Diffusion Bonding and Layer Transfer of Lithium Niobate to Silicon

OPTICAL, ELECTRICAL AND STRUCTURAL PROPERTIES OF PECVD QUASI EPITAXIAL PHOSPHOROUS DOPED SILICON FILMS ON CRYSTALLINE SILICON SUBSTRATE

Organic Solar Cells. Green River Project

Available online at ScienceDirect. Energy Procedia 55 (2014 )

Transmission Mode Photocathodes Covering the Spectral Range

Schottky Tunnel Contacts for Efficient Coupling of Photovoltaics and Catalysts

"Plasma CVD passivation; Key to high efficiency silicon solar cells",

Basics of Solar Photovoltaics. Photovoltaics (PV) Lecture-21

Topics Relevant to CdTe Thin Film Solar Cells

Effect of Growth Process on Polycrystalline Silicon Solar Cells Efficiency.

PV research in Neuchâtel: from high efficiency crystalline cells to novel module concepts

A Study towards the Possibility of Ultra Thin CdS/CdTe High Efficiency Solar Cells from Numerical Analysis

Simulation of Performance of Cadmium Telluride Solar Cell Using AMPS-1D Program

Basic efficiency limits, recent experimental results and novel light-trapping schemes in a-si:h, lc-si:h and micromorph tandem solar cells

The Effects of the Adding V2O5 on the Oxide Semiconductor Layer of a Dye-sensitized Solar Cell

Hei Wong.

Device Simulation of Grain Boundaries in Lightly Doped Polysilicon Films and Analysis of Dependence on Defect Density

Band-gap grading in Cu(In,Ga)Se 2 solar cells

Solar 101 for the Duke Energy Academy

Surface Preparation Challenges in Crystalline Silicon Photovoltaic Manufacturing

Temporal stability of a-si:h and a-sin x :H on crystalline silicon wafers

Back Surface Recombination Effect on the Ultra- Thin CIGS Solar Cells by SCAPS

2. High Efficiency Crystalline Si Solar Cells

OPTIMIZATION OF A FIRING FURNACE

Impurity Photovoltaic Effect in Multijunction Solar Cells

HF last passivation for high efficiency a-si:h/c-si heterojunction solar cells

The Effect of Front and Rear Surface Recombination Velocities on the Photocurrent of Buried Emitter Silicon Solar Cell

Review Article High-Efficiency Crystalline Silicon Solar Cells

Available online at ScienceDirect. Energy Procedia 92 (2016 ) Mathieu Boccard*, Nathan Rodkey, Zachary C.

HOMEWORK 4 and 5. March 15, Homework is due on Monday March 30, 2009 in Class. Answer the following questions from the Course Textbook:

Passivating contacts for crystalline silicon solar cells

Effects of CdCl 2 treatment on ultra-thin MOCVD-CdTe solar cells

Study of ICP-CVD grown Amorphous and Microcrystalline Silicon thin films in HIT structure

Presented at the 29th European PV Solar Energy Conference and Exhibition, September 2014, Amsterdam (NL)

Point-contacting by Localised Dielectric Breakdown: A new approach for contacting solar cells

Correlation Between Energy Gap and Defect Formation of Al Doped Zinc Oxide on Carbon Doped Silicon Oxide

I. GaAs Material Properties

Citation for the original published paper (version of record):

Silicon Thin Films. Hydrogenated amorphous silicon: hydrogen related structure and metastability Arno Smets (Delft University of Technology)

Project III. 4: THIN FILM DEVICES FOR LARGE AREA ELECTRONICS

Thin film PV Technologies Thin film Silicon PV Technology

Griddler 2.5: Full Area Bifacial Cell Design and Simulation for the Bifi PV Community

METALLIZATION OF PASSIVATING AND CARRIER SELECTIVE CONTACTS: STATUS AND PERSPECTIVES AT FRAUNHOFER ISE

Introduction to Solar Cell Materials-I

Production of PV cells

The Pennsylvania State University. The Graduate School. Graduate Program in Materials Science and Engineering

Solar Spectrum. -Black body radiation. Light bulb 3000 K Red->Yellow->White Surface of Sun 6000 K

Supporting Online Material for

High Transmittance Ti doped ITO Transparent Conducting Layer Applying to UV-LED. Y. H. Lin and C. Y. Liu

SUPPLEMENTARY INFORMATION

Anodic Aluminium Oxide for Passivation in Silicon Solar Cells

Available online at ScienceDirect. Energy Procedia 77 (2015 )

Recap of a-si and a-si cell technology Types of a-si manufacturing systems a-si cell and module manufacturing at Xunlight. Xunlight Corporation

Figure 16.31: Two-dimensional representations of (a) a quartz crystal and (b) a quartz glass.

ARTICLE IN PRESS. Solar Energy Materials & Solar Cells

Epitaxial Silicon Solar Cells

Device Architecture and Lifetime Requirements for High Efficiency Multicrystalline Silicon Solar Cells

Hsin-Ying Lee, 1 Ting-Chun Wang, 1 and Chun-Yen Tseng Introduction

Nanoscience in (Solar) Energy Research

Light Trapping Enhancement in Thin Film Silicon Solar Cell with Different Back Reflector

Application of infrared thermography to the characterization of multicristalline silicon solar cells

Halbleiter Prof. Yong Lei Prof. Thomas Hannappel

Inline-high-rate thermal evaporation of aluminum for novel industrial solar cell metallization

Photon Enhanced Thermionic Emission for Solar Energy Harvesting Progress Report to the Global Climate and Energy Project

Towards scalable fabrication of high efficiency polymer solar cells

R&D ACTIVITIES AT ASSCP-BHEL,GURGAON IN SOLAR PV. DST-EPSRC Workshop on Solar Energy Research

Research and Development of High-Voltage CIS-Based Thin Film Solar. Cells for Industrial Technology

Simulation study on the active layer thickness and the interface of a-igzo-tft with double active layers

School of Photovoltaic and Renewable Energy Engineering

DEVELOPMENT OF HIGH EFFICIENCY FLEXIBLE CdTe SOLAR CELLS

RECORD EFFICIENCIES OF SOLAR CELLS BASED ON N-TYPE MULTICRYSTALLINE SILICON. Center, ISC-Konstanz, Rudolf-Diesel-Str. 15, D Konstanz, Germany

Polycrystalline CdS/CdTe solar cells

Radial pn junction silicon devices have been extensively

Hall coefficient, mobility and carrier concentration as a function of composition and thickness of Zn-Te thin films

Semiconductor Technology

Materials Characterization

Assignment Questions

Synthetic antiferromagnet with Heusler alloy Co 2 FeAl ferromagnetic layers

Suitable Semiconductor Properties for Photoelectrochemistry Evaluated using Nitrides

Solar Energy. Photovoltaics Theory and Application. ME 430 Queen s University

Transcription:

Journal of Physics: Conference Series Simulation of High Efficiency Heterojunction Solar Cells with AFORS-HET To cite this article: Wang Lisheng et al 2011 J. Phys.: Conf. Ser. 276 012177 View the article online for updates and enhancements. Related content - Improvement of parameters in a-si(p)/c- Si(n)/a-Si(n) solar cells Mohammed Moustafa Bouzaki, Michel Aillerie, Sidi Ould Saad Hamady et al. - Heterojunction silicon solar cells performance optimization and sensitivity reduction to the interface defect states Hadi Bashiri, Mohammad Azim Karami and Shahramm Mohammad Nejad - Influence of interface states, conduction band offset, and front contact on the performance of a-sic:h(n)/c-si(p) heterojunction solar cells Zhi Qiao, Jian-Li Ji, Yan-Li Zhang et al. Recent citations - Simulation approach for optimization of device structure and thickness of HIT solar cells to achieve 27% efficiency Neeraj Dwivedi et al This content was downloaded from IP address 148.251.232.83 on 19/06/2018 at 19:12

Simulation of High Efficiency Heterojunction Solar Cells with AFORS-HET Wang Lisheng 1, Chen Fengxiang and Ai Yu Department of physics science and technology, Wuhan University of Technology, Wuhan city, Hubei province, P.R.China. 430070 E-mail:wang_lesson@yahoo.com.cn Abstract: In this paper, the high efficiency TCO/a-Si:H (n)/a-si:h(i)/c-si(p)/uc-si(p + )/Al HIT (heterojunction with intrinsic thin-layer) solar cells was analyzed and designed by AFORS-HET software. The influences of emitter, intrinsic layer and back surface field (BSF) on the photovoltaic characteristics of solar cell were discussed. The simulation results show that the key role of the intrinsic layer inserted between the a-si:h and crystalline silicon substrate is to decrease the interface states density. If the interface states density is lower than 10 10 cm -2 V -1, thinner intrinsic layer is better than thicker one. The increase of the thickness of the emitter will decrease the short-current density and affect the conversion efficiency. Microcrystalline BSF can increase conversion efficiency more than 2 percentage points compared with HIT solar cell with no BSF. But this BSF requires the doping concentration must exceed 10 20 cm -3. Considered the band mismatch between crystalline silicon and microcrystalline silicon, the optimal band gap of microcrystalline silicon BSF is about 1.4-1.6eV. 1. Introduction Silicon HIT solar cell is now attracting much interest, because very high efficiencies (above 22%) have been demonstrated [1]. This technology is based on Sanyo s successful technologies for fabricating high-quality a-si film and solar cells with low plasma damage processes. It combines many advantages, such as (1) it simultaneously enables an excellent surface passivation and p-n junction, resulting in high efficiency, (2) its low-temperature processed (<200 ) can prevent any degradation of bulk quality that happen with high-temperature cycling processes in low-quality silicon materials such as solar grade Czochralski (CZ) Si, and (3) compared with conventional diffused cells, a much 1 This paper is supported by The Fundamental Research Funds for the Central Universities (Project number:2010-la-019) Published under licence by Ltd 1

better temperature coefficient can be obtained with high V oc cells [2]. Though SANYO s original design used an n-type substrate as the absorber for the HIT solar cell, current researches concentrate on developing the HIT solar cell on a p-type substrate, because of its popularity in the photovoltaic industry [3]. However, inferior performance was observed for devices fabricated on c-si(p) as compared with those on c-si(n) except SANYO, which suggests that a further investigation is necessary to fully understand the factors that affect the performance of the HIT solar cells. Owing to the large numbers of processing variables, such as the thickness of the emitter, the interface states density, the band gap of the back surface field. etc, it is a formidable task to scrutinize the effect of each variable on the performance of the solar cell experimentally. Fortunately, numerical simulation software AFORS-HET provides a convenient way to evaluate the role of the various parameters present in the fabrication processing of HIT solar cells, and its reliability is proved by many references [3-4]. In this paper, the HIT solar cell TCO/a-Si:H (n)/a-si:h(i)/c-si(p)/uc-si(p + )/Al is designed and optimized by AFORS-HET. 2. Simulation Figure 1 The gap state distribution in n-type a-si:h layer. AFORS-HET is a very general computer simulation code for analyzing and designing the HIT solar cells, device operation may be in dark or under light illumination. In our simulation process, the solar AM1.5 radiation was adopted as the illumination source with a power density of 100mW/cm 2. The front and back contacts are assumed as flat band ones to neglect the contact potential influence. The surface recombination velocities of both electrons and holes are set as 10 7 cm/s. For the density of localized states in the band gap of amorphous silicon it has been assumed that there are both acceptor-like states and donor-like states modeled by exponential band tails (Urbach tails) and Gaussian mid-gap states (associated to silicon dangling bonds) [5]. The detailed localized states distributions are shown in figure 1. Other parameters used in the simulation are given in Table 1. During the simulations, all the parameters were adopted as the above setting values except for the specific declared one. 2

Table 1 The parameters used for HIT solar cell simulation Parameter and units e,h for electrons and holes, respectively a-si(n) a-si(i) c-si(p) uc-si(p+) Thickness(nm) variable variable 30000 variable Dielectric constanct 11.9 11.9 11.9 11.9 Electron affinity (ev) 3.9 4 4.05 4 Bandgap(eV) 1.72 1.72 1.12 variable Effective conduction band density (cm -3 ) 10 20 10 20 2.8 10 19 10 20 Effective valence band density (cm -3 ) 10 20 10 20 1.04 10 19 10 20 Electron mobility (cm 2 V -1 s -1 ) 5 5 1041 10 Hole mobility (cm 2 V -1 s -1 ) 1 1 412 3 Acceptor concentration (cm -3 ) 0 0 1.5 10 16 variable Donor concentration (cm -3 ) 2.49 10 19 100 0 0 Thermal velocity of electrons(cm.s -1 ) 10 7 10 7 10 7 10 7 Thermal velocity of holes(cm.s -1 ) 10 7 10 7 10 7 10 7 Layer density (g.cm -3 ) 2.328 2.328 2.328 2.328 Auger recombination coefficient for electron(cm 6.s -1 ) 0 0 2.2 10-31 0 Auger recombination coefficient for hole(cm 6.s -1 ) 0 0 9.9 10-32 0 Direct band-to-band recombination coefficient (cm 3.s -1 ) 0 0 0 0 3. Results and Discussion Figure 2 The J-V curve of TCO/a-Si:H(n)/c-Si(p)/Al cell. Firstly, the HIT solar cell with structure TCO/a-Si:H(n)/c-Si(p)/Al is considered. The emitter thickness is 5nm and the related J-V curve is shown in figure 2. The detailed photovoltaic properties are V oc =646.1mV, J sc =35.32mA/cm 2, FF=82.68% and E ff =18.87%. 3.1. Optimization of the intrinsic layer While SANYO develops HIT solar cells with a very thin intrinsic a-si:h(i) layer inserted between different type a-si:h and c-si, it should be noted that there is a controversy on the need of such a layer 3

[1]. Some authors claim that it is beneficial, while others get good results without it and do not see significant improvements if they introduce it. One reasonable explanation for the benefit of this undoped buffer layer is that the density of states in undoped a-si:h is weaker than in doped a-si:h, so we can expect to have less interface defects when the heterointerface is formed with undoped rather than doped a-si:h. 3.1.1. The influence of the intrinsic layer thickness. Figure3 shows the photovoltaic characteristics of HIT solar cell varies with the thickness of the intrinsic layer. As can be seen from figure 3, the conversion efficiency of HIT is up to 19.27% with 1nm intrinsic layer is inserted, which is 0.4% higher compared with no intrinsic layer solar cell. But as I layer thickness increases, the conversion efficiency decreases. When the thickness of intrinsic layer reaches 5nm, the conversion efficiency of HIT solar cell with intrinsic layer is equivalent to solar cell with no intrinsic layer. Compromising the production processes and the conversion efficiency, the optimal intrinsic layer thickness should be set at 3nm. In addition, in figure 3, we can see that with the increase of intrinsic layer thickness, the open circuit voltage keeps almost unchanged, while the short circuit density decreases. This is because, as I layer thickness increases, the electric field strength of space charge region decreases and the short-spectrum absorption of amorphous silicon increases. The corresponding light-induced carriers cannot be effectively collected, resulting in the decrease of short circuit current density. Figure 3 The performances of TCO/a-Si:H(n)/a-Si:H(i)/c-Si(p)/Al solar cell vary with the intrinsic layer thickness. 4

3.1.2. The influence of the interface states density. In the actual solar cell production processes, the impact of interface states density can not be ignored. In this paper, the interface defect states is assumed to insert between the I-type a-si:h layer and the p-type crystalline silicon substrate layer. The interface states is assumed as continuous donor-like states and acceptor-like states with average distribution in the band gap, with the capture cross sections of electron and hole both are 10-14 cm 2. The interface state density of a-si:h/c-si interface varies from 10 9 cm -2.eV -1 to 10 13 cm -2.eV -1, and its effects on the photovoltaic performances of the solar cell are shown in figure 4. Figure 4 The impact of interface states density on the photovoltaic performances of HIT solar cell. As shown in figure 4, when the interface state density D it is lower than10 10 cm -2.V -1, the solar cell performance is almost unaffected by the interface states. Under this condition, the efficiency is still maintained at 18.97 %, which is very close to the one without interface states. But when the interface state density D it is higher than 10 10 cm -2.V -1, the open circuit voltage decreased rapidly, while short-circuit current remained almost unchanged, even increased slightly. When the interface state density varies from 10 10 cm -2.V -1 to 10 12 cm -2.V -1, the open circuit voltage decrease from 644.5mV down to 549.2mV and the short-circuit current density increased from 35.52mA/cm 2 to 35.56mA/cm 2. The open circuit voltage reduction results from the increase of p-n junction reverse saturation current. As the interface defect density increases, the carrier recombination probability at interface increases, leading to the increase in reverse saturation current and reducing the open circuit voltage and fill factor. In order to obtain high-efficiency HIT solar cells, surface passivation methods is essential, such as plasma-assisted H passivation, in order to control interface state defect density as down to 10 10 cm -2.V -1 as possible. 3.2. Optimization of the emitter 5

Because of the structural disorder and high doping concentration, the diffusion length of carriers in amorphous emitter is shorter. And there is only drift current and no diffusion current in emitter. Furthermore, since the amorphous silicon layer is a higher-doping layer, the space charge region of p-n junction has fallen mainly on crystalline silicon side, and even there is no electric field in amorphous silicon area. Therefore, emitter should be prepared as thin as possible with higher-doping. Figure 5 shows the influence of emitter thickness on the photovoltaic properties of solar cell. Figure 5 The impact of emitter thickness on the photovoltaic properties of HIT solar cell. As can be seen from figure 5, with the emitter thickness increases, the open circuit voltage changed little, while the short circuit current is dramatically reduced. This is because with the emitter thickness increases, the absorption of the photon in emitter has increased. Considering the large amounts of recombination centers and the feature of no electric field in emitter, the photo-induced carriers are impossible to reach the edge of space charge region and contribute to light current. On the contrary, they will be recombined in the region and disappeared, resulting in the reduced short-spectrum response and short-circuit current. Fill factor also decreases with the n region thickness increases, for as the n layer thickness increases, series resistance will increases, which will reduce the fill factor [6]. Taking into account the technical production and conversion efficiency, the emitter thickness should be chosen about 5nm. 3.3. Optimization of the back surface field The back surface field refers to a region which will be a barrier of minority carrier, thereby reducing the recombination rates of photo-induced carriers in the back surface. This barrier can not only enhance the photocurrent, but also increase photo-voltage to some extent [7]. Similar with the crystalline silicon solar cells, the back surface field can be achieved by a layer doped with the same type, but with higher doping concentration. Considering the microcrystalline silicon is a coexistence 6

mixed phase of nano-crystalline silicon, amorphous silicon and grain boundaries, and its band gap can be adjusted by crystallization phase, so microcrystalline silicon was used to simulate the back-field structure [8], and the corresponding parameters were listed in Table 1. 3.3.1. The effect of the BSF thickness. Table 2 lists the effect of the BSF thickness on the J-V characteristics of HIT solar cell. We can see that when the thickness increases, the open circuit voltage and short circuit current are nearly unchanged, and only a slight change in fill factor and cell efficiency. So the efficiency can be considered as independent with the thickness [7], which is very beneficial. In view of technology, the difficulty in fabricating microcrystalline silicon thin films is that its growth rate is very low. If the required thickness is very thin, production time can be dramatically saved. Therefore, the BSF thickness can be set to 5nm. Table 2 The J-V characteristics of solar cell with the variation of the BSF thickness thickness/nm V oc /mv J sc /ma.cm -2 F.F/% /% 2 646.1 35.53 82.65 18.97 5 646.1 35.53 82.67 18.98 10 646.1 35.54 82.69 18.98 15 646.1 35.54 82.71 18.98 20 646.1 35.55 82.73 19 3.3.2. The effect of the BSF doping concentration. Figure 6 shows the influence of BSF doping concentration on the photovoltaic characteristics of solar cell. In figure 6, we can see that the doping Figure 6 The influence of the doping concentration of BSF. concentration must reach a certain value, preferably more than 10 20 cm -3 before the conversion 7

efficiency increased by 0.5 percentage points. With the doping concentration increases, the open circuit voltage, short circuit current, fill factor and so on all are increasing. This can attributed to the BSF band structure. When the doping concentration is low, the reflection role of BSF is not clear, and the barrier on the carrier transport can be reduced by increasing the doping concentration, which explains the high doping concentration is guarantee of good BSF. Microcrystalline silicon, due to its special crystallization phase and possibility of high doping, is selected as the BSF here. 3.3.3. The effect of the BSF band gap. With the development of thin film deposition process, silicon films with different band gap can be prepared under different deposition temperature and crystallization phase. In general, the band gap of nano-crystalline silicon is larger, up to 1.8eV and above; followed by amorphous silicon, the band gap of crystalline silicon is smaller. Therefore, we supposed that the band gap of microcrystalline silicon varies between 1.12eV-2eV. Figure 7 shows the parameters of HIT solar cells vary with the band gap of microcrystalline silicon BSF. Figure 7 The effect of BSF band-gap. As can be seen from figure 7, with the band gap increases, the open circuit voltage increased firstly and then saturated; when the band gap is above 1.8eV, the fill factor began to decrease, leading to decreased conversion efficiency. The reason is that with the band gap broader, although the BSF provides barrier for the minority carriers, but also will impede the transport of majority carriers. So an optimum band gap should be selected in the 1.4-1.6 ev. 3.3.4. The comparison between HIT solar cells with and without BSF. Figure 8 shows the effciency comparison of HIT solar cells before and after the BSF. In figure 8, we can see that with the BSF is added, all parameters of solar cells are improved, indicating a reasonable design of the back field, especially the optimization of doping concentration and the band gap, will increase the conversion efficiency by 2 percentage points. The detailed information also shown in figure. 8

Type V oc J sc FF n/i/p 647.7 35.73 82.55 19.1 n/i/p/p + 678.9 37.35 83.97 21.29 Figure 8 The J-V curves for HIT solar cells before and after BSF. As shown in figure 8, when a 5nm thick microcrystalline silicon layer with doping 10 20 cm -3 and band gap 1.5eV is added as BSF, it can improve V oc 30mV and raise J sc 1.5mA/cm 2. 4. Conclusion In this paper, the TCO / a-si: H (n) / a-si: H (i) / c-si (p) / uc-si (p + ) / Al HIT solar cell is simulated with AFORS-HET software. It was found that with the increase of the n layer thickness, the short-spectrum response and the short-circuit current density of solar cell will decrease; the introduction of intrinsic amorphous silicon layer is mainly to reduce the interface defect state density. If the interface states density is lower, the better intrinsic thickness is no more than 5nm. The thickness of BSF has little effect on the conversion efficiency, while reasonable doping and band gap design can improve efficiency more than 2 percentage points. At that time, the photovoltaic parameters of HIT solar cells are V oc =678.9mV, J sc =37.35mA/cm 2,F.F=83.97%, =21.29%. References: [1] Kleider J P, Chouffot R, Gudovskikh A S, Cabarrocas P R and Labrune M 2009 Thin Solid Films 517 6386-6391 [2] Mikio T, Hitoshi S, Yukihiro Y, Eiji M, Akira T and Makoto T 2005 PVSEC 501 [3] Zhao L, Zhou C L, Li H L, Diao H W and Wang W J 2008 Solar Energy Materials & Solar cells. 92 673-681 [4] Ren B Y, Zhang Y, Guo B, Zhang B, Li H Y, Xu Y and Wang W J 2008 Acta Energy Solaris Sinica 29(9) 1112-1116 [5] Norberto H C and Arturo M A 2010 Solar Energy Materials & Solar cells 94 62-67 [6] Ren B Y, Wang M H, Liu X P and Li Y L 2008 Acta Energy Solaris Sinica 29(2) 126-129 [7] Zhao L, Zhou C L, Li H L, Diao H W and Wang W J 2008 Acta Physcia Sinica 57(5) 3212-3217 [8] Li L M, Zhou B Q, Chen X, Han B and Hao L Y 2009 Information Recording Materials 10(3) 18-21 9