Issue 
EPJ Photovolt.
Volume 9, 2018
Special Issue: French PV Days 2016



Article Number  1  
Number of page(s)  11  
Section  Modelling  
DOI  https://doi.org/10.1051/epjpv/2017008  
Published online  29 January 2018 
https://doi.org/10.1051/epjpv/2017008
Regular Article
Temperature and color management of silicon solar cells for building integrated photovoltaic
^{1}
Univ Lyon, CNRS, INSALyon, Université Claude Bernard Lyon 1, CETHIL UMR5008, 69621 Villeurbanne, France
^{2} Université de Lyon, Institut des Nanotechnologies de Lyon INL – UMR5270 CNRS INSA Lyon, 69621 Villeurbanne, France
^{*} email: fabien.mandorlo@insalyon.Fr
Received:
27
March
2017
Received in final form:
12
July
2017
Accepted:
7
September
2017
Published online: 29 January 2018
Color management of integrated photovoltaics must meet two criteria of performance: provide maximum conversion efficiency and allow getting the chosen colors with an appropriate brightness, more particularly when using side by side solar cells of different colors. As the cooling conditions are not necessarily optimal, we need to take into account the influence of the heat transfer and temperature. In this article, we focus on the color space and brightness achieved by varying the antireflective properties of flat silicon solar cells. We demonstrate that taking into account the thermal effects allows freely choosing the color and adapting the brightness with a small impact on the conversion efficiency, except for dark blue solar cells. This behavior is especially true when heat exchange by convection is low. Our optical simulations show that the perceived color, for single layer ARC, is not varying with the position of the observer, whatever the chosen color. The use of a double layer ARC adds flexibility to tune the wanted color since the color space is greatly increased in the green and yellow directions. Last, choosing the accurate material allows both bright colors and high conversion efficiency at the same time.
Key words: CIE 1931 color space / BIPV / thermal management / operating temperature / ARC / conversion efficiency / silicon solar cell
© M. Amara et al., published by EDP Sciences, 2018
This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.
1 Introduction
Building integrated photovoltaics (BIPV) [1] is one of the solutions to support and develop renewable and nonpolluting energy. As a rule, performances of solar panels depend on the material quality and solar cells architecture, but also on cell temperature condition since high temperature degrades the conversion efficiency [2]. Nowadays, under the present conditions, the BIPV are not favorable to the cooling of solar cells. In addition, the esthetics of the PV panels must also be taken into account, including their color related to the urban environment.
In order to determine the cell or module operating temperature, it is possible to use experimental correlation such as those listed and discussed in the review of Skoplaki and Palyvos [3] and Kurnik et al. [4]. Even if these correlations calculate the module temperature as function of environmental conditions, they do not make it possible to link the thermal equilibrium as a function of the cells characteristics, in particular the antireflective layers. Moreover, the optimization of solar cells is done at standard test conditions (STC, 25 ^{°}C) which ignores the actual temperature of the cell. It is necessary to use a thermal model in addition to the electrooptical ones to take into account heat exchange in order to calculate the operating temperature.
To adjust the perceived color of a solar cell observed outside a building, we can control the reflection on its top surface, and so adjust the thickness of the antireflective coating (ARC) layers [5]. Such a solution is both easy and flexible since we only need to control the duration of a deposited layer. However, we have to balance between the brightness/contrast of the reflected light, and the number of photons producing the electrons and holes. Indeed, maximizing the short circuit current J_{sc}, which corresponds to the number of carriers that are photogenerated and collected in electrodes, is a simple way to optimize the conversion efficiency at STC [6–9] for which the solar cell is cooled at the ambient temperature 25 ^{°}C, with no heat transfer. In real configurations, thermal conditions and solar irradiation are major varying factors that strongly affect the conversion efficiency of photovoltaic cells. As a consequence, thermal effects are likely to be required for the optimization analysis of ARC layers of PV cells. It is clear from this observation that the obvious criterion is the maximization of the electrical conversion efficiency. Ideally, the aim is to satisfy two conflicted requirements: obtain both the wanted colors and a high conversion efficiency.
In this article, we propose to determine the influence of the ARC thickness on the color and conversion efficiency taking into account the thermoelectrooptical properties.
2 Modeling of solar cells performances and colors
2.1 Solar cell simulation
For cell simulation, we use our inhouse computer simulation tool [10] based on the resolution of three interdependent physical aspects. The first one is the heat transfer, which includes the following source terms: local heat due to radiation absorption and thermal emission (Q_{R}), nonradiative recombination (Q_{NRR}), excess power released by thermalization (Q_{T}), Joule effect (Q_{J}) which corresponds to the interaction between the charge carriers current and the electrical field inside the cell, and Thomson effect (Q_{TM}) due to the heat carried and extracted from the cell by the charge carriers. The second set is a spectral radiative transfer, which is solved through a gray per band approach and separated treatments of the direct and diffuse components of radiation fluxes. The third set of equations is related to generation and diffusion of electrons and holes in the device.
All the equations are solved with spectral and thermal variation of all used parameters. So, in range [0.2–15] μm, the incident spectrum is decomposed in several spectral bands. For the boundary conditions, in addition to the front and back surface recombination velocity respectively, we take into account current flow at metal–semiconductor interfaces thanks to the Peltier effect (Q_{P}) and the heat exchange by convection (Q_{Conv}).
2.2 Optical calculation: conversion of an optical spectrum to a color in the CIE 1931 color space
In this part we focus on the color prediction of a solar cell, by using its reflected spectrum. The CIE 1931 XYZ color space has been created by the International Commission Illumination (CIE) [11] in 1931 in order to link the visible spectrum with the perceived colors of the human eye. To match the sensitivity of the shortmiddlelong wavelength of “standard” human cone cells, this color space relies on 3 color matching functions , , and (Fig. 1).
These functions are used to calculate the tristimulus values X, Y and Z by projection of the optical signal to be converted S(λ) on the visible range from 380 to 780 nm [12]. For this article, the optical signal S(λ) will be the reflected light (R(λ)) coming from a normalized power density spectrum AM_{x}(λ), leading to S(λ) = R(λ) ⋅ AM_{x}(λ): (1) with: (2)
Among coordinates (X, Y, Z), Y represents the brightness of the signal, since the human eyes perceive light mainly through the middle range cone cells, which correspond to the green color detection. This parameter necessarily varies between 0 (no visible light) and 1 (all visible light is reflected). The chromaticity of a color, that quantifies the balance between each elementary wavelengths of the signal, is described by two derived parameters x and y defined such as: (3)
Then, x, y and Y scalars correspond to the CIE 1931 xyY color space. Its chromaticity diagram, also called gamut (Fig. 2) represents the domain covered by x and y for any visible signal: it is delimited by the single wavelength optical signals and any linear combination of the 380 and 780 nm ones (also called the “purple line”).
We can also use the tristimulus values to convert them into the sRGB system using [13] to calculate the Red, Green and Blue (RGB) coordinates for computers and displays. As this color space is smaller than the CIE 1931 xyY one (Fig. 2), colors out of the sRGB domain are approximated: the perceived color on a display (and so on this article) is slightly different from the real color of the simulated structure.
Fig. 1 CIE 1931 standard observer color matching functions. 
Fig. 2 CIE1931 color space chromaticity diagram including the sRGB one (used for displays). 
3 Simulation inputs parameters
3.1 Description of the reference solar cell
A nonencapsulated flat photovoltaic cell is considered with a generic structure composed of a standard p–n junction with a negatively doped upper zone (ntype, emitter), a positively doped lower zone (ptype, base) and a depletion zone or charge space (Fig. 3).
In addition, the front surface is covered by a stack of dielectric thin films (ARC) in order to minimize the reflection of the incident light. The cell’s upper face is subject to collimated and diffuse isotropic solar irradiations. The incidence of collimated beam is characterized by the incidence angle θ. Additional diffuse isotropic radiation fluxes are considered at both boundaries and refer to surroundings’ infrared radiation emission. The surroundings are considered as black bodies at T_{+∞,top} and T_{+∞,bot} and the ambient temperature is noted T_{amb}. To get the optical response of the ARC stack on the solar cell, we use the transfer matrix formalism [14] that has been validated by comparison with [15] to calculate the reflected spectrum of a multistack system and convert it into the CIE 1931 xyY color space. The studied cell has the parameters given in Table 1.
The used incident spectrum is a standard AM_{1.5G} [16]. The influence of air mass is investigated by also using AM_{1G} and AM_{3G} spectra (Fig. 4). As for the coefficient of convection on each side of the module h_{conv}, its value depends on the configuration of the module integration. Thus, there are many empirical correlation depending on the wind velocity (forced convection) or ambient temperature (free convection). For example, Jones and Underwood analysis [17] indicates a large variation of the heat coefficient in between 1.2 and 9.6 W⋅m^{−2}⋅K^{−1}. The last value corresponds to forced convection condition which is greatly favorable to the cooling of the modules. Inorder to evaluate the impact of ARC on the cell temperature, we choose two values of h_{conv}: 2 and 5 W⋅m^{−2}⋅K^{−1}, which corresponds respectively to disadvantageous and favorable cooling conditions. In case of standard conditions, h_{conv} is equal to + ∞ leading to an operating temperature of 25 °C equal to the temperatures T_{+∞,top}, T_{+∞,bot} and T_{amb}.
Fig. 3 Description of the simulated solar cell and its boundary conditions. 
Parameters and their default values used in simulation.
Fig. 4 Spectral irradiances for different air masses. 
3.2 Optical and thermophysical properties
In this article, we use two types of data: experimental measurements and database (Fig. 5). To calculate the infrared radiation emission and exchange with the atmosphere, we need to know the dispersion curve of all the materials until λ = 15 μm. Only few materials are known from visible to long wavelength infrared. Due to the lack of experimental data for other materials, we decide to use Philipp [18] and Kischkat et al. [19] for SiO_{2}. As for Si_{3}N_{4}, we use the data given by Gao et al. [20] and Kischkat et al. [19]. It should be noted that the concatenation procedure of these data creates a discontinuity for Si_{3}N_{4} at 1.5 μm. However, around this value, there is no impact on optical absorption and thermal behaviors.
Nevertheless, industrial solar cells are based on hydrogenated SiN_{x} and SiO_{x}. In order to predict colors with these materials, we will use experimental data in the visible range. For this purpose, both SiO_{x}:H and SiN_{x}:H layers are deposited by LFPECVD with different stoichiometries [21,22]. One of the SiN_{x}:H layers (n = 2.08 at λ = 600 nm) is optimized for both passivation and optical efficiency for standard flat silicon solar cells with no glass and no EVA. Then, a Horiba Jobin Yvon UVISEL ellipsometer is used to extract their optical indexes in the visible range (Fig. 6).
To take into account the temperature variation of the semiconductor models, we implement Fermi–Dirac free carriers statistics, freecarriers mobility [23], intrinsic carrier density [24], bandgap and bandgap narrowing [25,26] with temperature dependent models. Last, the transport of photogenerated carriers is governed by a standard diode equation. Details related to the determination of the band gap, optical and radiation properties, minority carrier diffusion coefficients, recombination rates and intrinsic concentrations can be found in our previous paper [24].
Fig. 5 Refractive index of SiO_{2} and Si_{3}N_{4} to cover both the visible and long wavelength infrared domains (data extracted from [18–20]. 
Fig. 6 Experimental optical indexes of different silicon oxides and nitrides (refractive index n at λ = 600 nm in brackets). 
4 Results and discussion
4.1 Achievable color space depending on the position of the observer
In order to design colored solar cells, we propose to study both the thickness t_{ARC} of a single layer ARC and the influence of the observation angle of an illuminated solar cell by an AM_{1.5G} spectrum.
The theoretical result for a Si_{3}N_{4} ARC layer is shown in Figure 7a. The chromaticity is compatible with Figure 7 in [27] or Figure 2b in [9] for the normal reflection (θ = 0). We observe that for a thin ARC (single mode regime), the color (Fig. 7a) and the brightness (Fig. 7b) are almost independent from the inclination of the solar cell. When the optical path of photons in the ARC layer is very long (i.e. for very high values of θ or t_{ARC}), the higher number of interferences lowers the color selectivity and the emitted light is close to a dark gray (Y close to 25%).
From Figure 7b, the brightness (and so the contrast), gets its minimum and maximum values for θ = 0. We can also observe that colored cells are quite dark since the domain where the perceived colors are different from gray (Fig. 7a) correspond to the minimum domains of the brightness Y (Fig. 7a). Figure 8 represents the chromaticity of the calculated colors. The dotted black line corresponds to the evolution of the color for a 75 nm thick ARC, from an inclination from 0^{°} to 85^{°} . With such a thickness, the solar cell is blue when observed with an AM_{1.5G} irradiance and its color varies in direction of the gray point x ≈ y ≈ 0.33 when θ increases. This particular trend is verified whatever the thickness of the ARC. Consequently, drawing the chromaticity for θ = 0 can be used to get all the achievable colors for a single layer ARC.
Figure 8 also shows that a small part of the visible colors can be obtained with a silicon nitride ARC, with ability to reproduce darkblue colors. Indeed, a standard 75 nm thick ARC (point B_{1} in Fig. 8) reflects blue light with a low brightness (Y = 1.8% in Fig. 7), that corresponds to the perceived light when a standard solar cell is observed. On the contrary, red and green colors cannot be easily obtained, since the extension of the gamut is very poor in direction to the 520 and 680 nm wavelengths.
Figure 9 shows the reflected power at normal incidence in the visible range for points R, G, B_{1} and B_{2} in the gamut diagram. These curves correspond to t_{ARC} = 60, 390, 75 and 360 nm respectively.
We can clearly observe that curves for R, G and B_{1} are well accurate to maximize X, Y and Z respectively and obtain red, green and blue solar cells. Indeed, in order to get a green color, the reflected power must be maximum around λ = 550 nm and minimum around 450 nm and 600 nm: for this reason, we need a thick multimode ARC, leading to t_{ARC} = 390 nm. It is important to note that the chromaticity diagram is not sufficient to infer the perceived color. Indeed, for the R and B points, the brightness is so low (Y < 3%) that side of a solar cell with a thickness corresponding to G (Y = 21%), the human eye will perceive them as almost black. This is confirmed by Figure 7a, even with the sRGB approximation.
Last, the comparison of spectrum for the points B_{1} and B_{2} shows that both are accurate to get a high value of X, but in case of B_{2} the Y and Z variables are higher (see Eq. (1)), leading to x and y values closer to the (0.33,0.33) point in the chromaticity diagram (Eq. (3)). The reflected power in case of B_{2} is multiplied by a factor 6.8, even if this color is less “pure” in terms of blue content.
Fig. 7 Perceived color converted into sRGB (a) and brightness (b) depending on the observation angle and ARC thickness for Si_{3}N_{4}. 
Fig. 8 CIE1931 color space chromaticity for a Si_{3}N_{4} ARC with an AM_{1.5G} irradiance. For θ = 0, arrows represent the direction for which t_{ARC} is increasing and round marks correspond to ARC thicknesses multiples of 25 nm. 
Fig. 9 Reflected power for different thicknesses at normal incidence. 
4.2 Influence of the irradiance on the perceived color
In this part, we propose to study whether the irradiance influences or not the color of the reflected light. As the atmosphere absorbs more light in the short wavelengths, the balance between short (blue color) and long wavelengths (red color) varies in disfavor of the blue content when the air mass increases. Consequently, the AM_{1D,G} irradiance spectra relatively contain more power in the short wavelengths than AM_{3D,G}. Here, we propose to calculate the response of the same Si_{3}N_{4} layer as in the previous part for AM_{1}, AM_{1.5} and AM_{3} spectrum (direct and global).
As previously explained, only the simulation for θ = 0 is necessary to extract the reachable chromaticity domain for a given single ARC layer. Consequently, after the simulation, we applied a Jarvis based method [28] to extract the convexhull of the data at normal incidence. Results are presented in Figure 10.
The shapes and delimited surfaces by the convex hull of the achievable colors are close for each spectrum. However, each domain is shifted and rotated in order to reduce the coverage in the blue side and increase it in the green and red ones when the air mass increases (i.e. from AM_{1D,G} to AM_{3D,G}). This is a consequence of a lower amount of blue photons reaching the surface of the solar cell (Fig. 4). The comparison of the direct and global spectra shows that direct spectra have a lower blue content than the global ones for a given air mass. Consequently, the gamut of a direct spectrum is always shifted in the red direction compared to that one of a global spectrum.
In all cases, the color of a given optical system is not affected a lot by the incident spectrum.
Fig. 10 Influence of the illumination on the perceived color (Si_{3}N_{4} ARC) at normal incidence for different irradiances. 
4.3 Gamuts using silicon nitrides and oxides
Figure 11 shows the color behavior for both silicon nitrides and oxides taking into account the experimental dispersion curves (Fig. 6).
As we can observe in Figure 11a, the difference between Si_{3}N_{4} and SiN_{x}:H (2.08) is quite limited: the stoichiometric nitride has a better extension in the green direction while the standard passivating ARC allows more orange hues. We can also observe that increasing the Si content reduces the gamut since the optical index gets closer to the Si one, leading to a less selective reflection at the front surface.
On the contrary, for substoichiometric oxides (Fig. 11b), increasing the amount of Si increases the surface of the chromaticity diagram in every direction (Fig. 11b). By comparing the behavior of all the SiO_{x} layers, we can deduce that lowering the optical index reduces the size of the convex hull in every direction. For Si rich oxides we get similar gamuts as for Si_{3}N_{4} and SiN_{x}:H since their optical indexes are close. However, in case of silicon oxides, increasing the Si concentration in the ARC helps designing the wanted color, but also increases the parasitic absorption due to a higher extinction coefficient in Figure 6.
From both figures, we can deduce that the widest gamut corresponds to an optical index close to n = 2. Consequently, SiN_{x}:H is a very interesting layer since it combines a wide range of colors with low absorption and good electrical passivation. This property is compatible with the study of the reflected spectra in the previous part, since the optical index of this standard SiN_{x}:H is close to (optical index of silicon and air respectively), and so, well adapted to maximize the transmission into the active material. As a consequence, it is also accurate to avoid partial reflection of unwanted wavelengths, leading to a better color selectivity.
Fig. 11 CIE 1931 color space chromaticity for different ARC layers with an AM_{1.5G} irradiance (optical index at λ = 600 nm in brackets). The green point in the center of the gamut corresponds the simulation with no ARC. 
4.4 Optical benefits of using a double layer ARC
In order to extend the reachable gamut, we propose to study a double layer ARC made of both SiO_{2} and Si_{3}N_{4}. This particular combination is not optimized to reduce the reflectivity [29].
The obtained colors with θ = 0, approximated in the sRGB domain, are presented in Figure 12a. The corresponding CIE 1931 xyY color space is shown in Figures 12b and 13. On this second figure, we also add the convex hull of the gamut for both single layer SiO_{2} and Si_{3}N_{4} ARCs. It must be mentioned that these two layers have the same behavior when t_{ARC} increases with the same sequence of colors to define the convex hull, but this latter is smaller for SiO_{2}.
The behavior of the double layer ARC is quite similar as the single one studied previously: when the optical path increases (either in SiO_{2} or Si_{3}N_{4}), the system becomes multimode and the color sensitivity decreases, converging toward a gray color. However, in the case of this study, the convergence is quicker with the increase of the silicon nitride thickness than for the silica one, due to its higher optical index. For instance, we observe limited variations of the brightness Y for nm while for nm, the color and the brightness still significantly vary.
The main advantage of using a double layer ARC is the wider extension of the chromaticity (Fig. 13) compared to the single one (Fig. 11), more particularly in the cyan, green, yellow and orange directions. On the contrary, in the red, magenta and blue directions, the chromaticity seems to be limited by the Si_{3}N_{4} layer.
Fig. 12 Perceived color ARC layer and brightness converted into sRGB for a SiO_{2}/Si_{3}N_{4} double at normal incidence. 
Fig. 13 CIE xyY 1931 color space for a doublestack ARC layer (SiO_{2}/Si_{3}N_{4}) with an AM_{1.5G} irradiance. 
4.5 Thermal and electrical behaviors
Figure 14 depicts plots of the cell’s temperature as a function of the ARC thickness for Si_{3}N_{4} (solid lines) and SiO_{2} (dashed lines). In the range 0–200 nm, the cell temperature is strongly sensitive to the ARC thickness regardless of the composition of the ARC layer and the heat coefficient h_{conv}. The peak temperature is achieved with d = 70 nm with a Si_{3}N_{4} coating and about 100 nm for SiO_{2}. Besides, it is obvious that increasing the heat coefficient induces more cooling of the cell temperature. Beyond d = 200 nm, the temperature shows less sensitivity to the ARC thickness even if it presents a second local maximum around 350 nm for Si_{3}N_{4} and 400 nm for SiO_{2}.
To explain the variation of the temperature with the ARC thickness, we display in Figure 15 the variation of the heat sources and sinks in case of h_{conv} = 5 W⋅m^{−2}⋅K^{−1}. As we can see, all sources exhibit the same trend: an increase followed by a decrease, and again, an increase and decrease. However, the second peak is less pronounced. In addition, the maximum of the curves does not occur at the same thickness of the ARC. It can also be observed that the heating is dominated by the thermalization, the Joule and the Peltier sources with less contribution for nonradiative recombination. As for cooling, it is mainly due to convection and radiation. The contribution of the Thomson effect is negligible.
It should be pointed out that the thermal sources are higher with Si_{3}N_{4} compare to SiO_{2}, thanks to a lower reflection at the front surface (Fig. 16). Actually, as explained in Section 4.3, the optical index of Si_{3}N_{4} is close to . The same effect may explain the trend of the Joule and the Thomson sources, since both phenomena are correlated to the photocurrent.
The antireflecting layers are usually designed in order to maximize the shortcircuit current density J_{sc} at STC where the cell temperature is maintained at 25 ^{°}C. This is illustrated in Figure 17 in the case of SiO_{2} and Si_{3}N_{4} layers. Our own finding indicates that whatever the ARC layer, the optimum thickness does not depend on the heat coefficient h_{conv}. This is readily explained by the fact that the bandgap is weakly influenced by temperature, and so the absorption range. Thus, we find an optimum thickness for SiO_{2} around 105 nm and 75 nm for Si_{3}N_{4}. It should be emphasized that the maximum of J_{sc} occurs very close to maximum of the temperature.
Now, considering the electrical conversion efficiency η with real environmental conditions as criterion (rather than J_{SC} at STC) leads to different conclusions as illustrated in Figure 18:
the decrease of the heat coefficient h_{conv} lowers the conversion efficiency η;
the optimum thickness depends on the thermal environment of the solar cell (Tab. 2) since the thermal heat and sources vary with the ARC thickness (see Fig. 15);

at high temperature (small value of h_{conv} in Figure 14), the variation of the conversion efficiency η is less pronounced.
Table 2 summarizes the optimal characteristic of the ARC regarding the thermal coefficient. It should be emphasized that calculated values with h = + ∞ is equal to those obtained with STC optimization criteria.
We can clearly observe that whatever the ARC material, the optimal thickness must be increased when the cooling conditions deteriorates, i.e. when h_{conv} decreases, even for realistic conditions with wind (h_{conv} = 5 W⋅m^{−2}⋅K^{−1}).
Fig. 14 Solar cell temperature as a function of the ARC thickness under AM_{1.5G} illumination. For h_{conv} = + ∞, temperature is 25 ^{°}C whatever t_{ARC}. Solid lines: Si_{3}N_{4}, dashed lines: SiO_{2}. 
Fig. 15 Heat and sink sources as a function of the ARC thickness (Solid lines: Si_{3}N_{4}, dashed lines: SiO_{2}). 
Fig. 16 Reflectivity of the solar cell at the maximum of the thermalization sources. 
Fig. 17 Short current circuit current density as a function of the ARC thickness (solid lines: Si_{3}N_{4}, dashed lines: SiO_{2}). 
Fig. 18 Absolute and relative conversion efficiencies as a function of the ARC thickness for both SiO_{2} and Si_{3}N_{4}. 
Optimal thickness as a function of the heat coefficient to maximize the conversion efficiency.
4.6 Color selection and its impact on the thermal and electrical performances
From previous analysis, we observed that the thickness of an ARC layer can be easily used to tune the color of a solar cell. However, taking into account the thermal effects, we have shown that beyond t_{ARC} = 125 nm, the electrothermooptical balance of the solar cell leads to similar conversion efficiencies.
This behavior is mainly the consequence of a lower thermalization due to absorption of blue photons, and less heat production, that lowers the conversion efficiency.
In Figure 19, we draw the limit conversion efficiencies with h = + ∞ (infinite convection) and h = 2 W⋅m^{−2}⋅K^{−1} (no wind), and superimposed the sRGB colors calculated for AM_{1.5G} and θ = 0. We can observe that the contrast of the colors achieved by the SiO_{2} ARC is lower than for Si_{3}N_{4}, that is compatible with the smaller chromaticity diagram for SiO_{2} (Fig. 12). In fact, this behavior is a direct consequence of the higher curvature (Fig. 16) of the reflectivity spectrum that helps in color selection. We can also notice that higher conversion efficiencies are achieved for the standard Si_{3}N_{4} layer, making it the ideal candidate to satisfy both performance constraint and color flexibility. In both cases, the peak of the maximum conversion efficiency only leads to dark orange, red, magenta and blue colors (sequence obtained thanks to the gamut in Fig. 8), and the maximum conversion efficiency does not correspond to the darkest blue. Indeed, these dark colors correspond to a high thermalization level that reduces the conversion efficiency due to a high operating temperature. Then, increasing the ARC thickness first leads to a gray color and, then, different colors are available for even thicker ARCs. As a consequence, except for few colors, the maximum conversion efficiency cannot be obtained and it may be difficult to balance the relative brightness between solar cells with different colors. However, for higher ARC thicknesses (t_{ARC} > 125 nm), the conversion efficiency is lower but do not vary a lot with the ARC thickness since the ARC also favors the absorption of infrared photons. Then, designers are almost free to tune t_{ARC} and get the color they need since cyan green and yellow colors can be also obtained. This can be done at the cost of a relative 5–10% decrease of the conversion efficiency. If the blue color is necessary for a multicolor BIPV project, one may consider the range 200–225 nm to balance its brightness with the other colors for Si_{3}N_{4}.
Last, the chromaticity diagram (Fig. 8) shows different intersections: for instance, the yellow color at point (x = 0.41, y = 0.45) can be obtained using a 155 nm or a 290 nm thick Si_{3}N_{4} layer. The respective brightnesses are Y = 0.32 and Y = 0.30, meaning that the thinner ARC will lead to little bit brighter solar cells, with a slightly lower conversion efficiency (Fig. 19a).
Fig. 19 Conversion efficiency and perceived color for different single layer ARC thicknesses (with θ = 0). 
5 Impact of the EVA and glass for Si flat solar cells
A realistic solar cell is actually protected by a front glass with an EVA layer to prevent its degradation. This structure may affect both optical and thermal behavior.
5.1 Optical impact
From Section 4.3, the best candidate for an extended gamut is expected to correspond to the material for which the transmission in silicon is very efficient. Taking into account the optical index of EVA [30], the optimal index is expected to be around . This is what we observe in Figure 20 since the widest gamut is obtained for n between 2.34 and 2.48 for SiN_{x}:H layers [21]. Then, the gamut extension in the red and blue directions are very close to the previous study without considering the EVA and glass materials (Fig. 11a), but the extension in the greencyan and greenyellow directions is better.
The main issue in using Si rich SiN_{x}:H layers is the increase of the parasitic absorption in the ARC layer. A standard ARC layer optimizing the conversion efficiency at STC leads to an optical index close to 2.22 [30]. Then, the gamut is mainly reduced in the red and green directions.
Fig. 20 CIE1931 color space chromaticity for different SiN_{x}:H ARC layers with an AM_{1.5G} irradiance (optical index at λ = 600 nm in brackets) when taking the EVA and glass layers into account. The green point in the center of the gamut corresponds the simulation with no ARC. 
5.2 Thermal impact of the encapsulation
The encapsulation of the solar cells affects the thermal balance by adding thermal resistances and modifies the thermal emission of the module. Regarding silicon/EVA/glass stack structure, the thermal energy balance can be expressed through an effective thermal resistance, R_{eff}, between the cell and the environment [31]: (4)where H is the total heat source, T_{c} the temperature of the cell and T_{amb} the ambient temperature. It is possible to calculate the value of R_{eff} through the total of its components for one side of the module: (5)where h_{rad} is the heat coefficient induced by radiative cooling, h_{conv} represents the convection effect, and R_{EVA} and R_{glass} are the thermal resistances of glass and EVA layers. With typical EVA and glass characteristics [32], R_{EVA} ≈ 2 × 10^{−3} W⋅m^{−2}⋅K^{−1}, R_{glass} ≈ 2 × 10^{−3} W⋅m^{−2}⋅K^{−1} and h_{rad} is around 7 W⋅m^{−2}⋅K^{−1}. When h_{conv} > 2 W⋅m^{−2}⋅K^{−1}, the contribution of R_{EVA} + R_{glass} is less than 3% in the final value of R_{eff}. This result is consistent with Weiss et al. [33] and shows that the thermal behavior of a solar module is dominated by heat convection and heat radiation of the solar cell.
Since glass and EVA absorb slightly solar radiation in the UV part which is mostly thermalized in a cell without encapsulation, the heat source, H, does not change much. Therefore, only the emission by thermal radiation is modified.
The variation of h_{rad} caused by the increase of the emissivity of the module from 0.8 to 0.9 [34] induces a variation of R_{eff} below 10%. This variation has the same impact on the thermal energy balance that a similar variation of h_{conv} which has been analyzed in Figures 14 and 19 on the nonencapsulated cell. The only case which is not treated correctly through the nonencapsulated cell is h_{conv} = 2 W⋅m^{−2}⋅K^{−1}. Indeed, the increase of the emissivity can slightly decrease the operating temperature. However, the trends described in Figures 14 and 19 are sufficient to understand the implications of such higher emissivity because the impact is limited to only few ^{°}C. Hence, this analysis can only be improved marginally with the consideration of the glass and the EVA layers whereas texturisation of the surfaces is an interesting challenge to take up.
6 Conclusion
In our study, we focused on a solar cell without any encapsulating materials (glass and EVA) on the front surface. However, these materials have a little impact on our conclusions as for colors, the standard ARC suitable for this kind of encapsulation should be accurate for the same reasons as the SiN_{x}:H we used. From a thermal point of view, the incorporation of these elements leads to use an effective value of the convection coefficient h slightly higher to that used in our demonstration.
The balance between the conversion efficiency and the possibility to tune the perceived color of a solar cell has been studied. We have demonstrated that increasing the ARC thickness, has a limited impact on the energy production, which implies that the color can be almost freely chosen. This behavior is the consequence of an almost constant thermalization level, and small variations of the thermal sources. When the convection cooling conditions deteriorate, the impact of choosing the solar cell’s color decreases further.
Concerning the color properties, we demonstrated that the thinner the ARC, the wider is the inclination for which the perceived color is not degraded. This color is almost independent from the irradiance spectrum, and so the latitude on Earth. For a single layer ARC, the standard Si_{3}N_{4} leads to more contrasted colors compared to substoichiometric SiO_{x} materials. As its optical index is close to that one of the standard SiN_{x}:H used for both surface passivation and optical transmission into silicon, we can deduce that this material is also a great candidate to cover a wide range of colors. Last, using a double layer ARC mainly extends the available colors in the cyan, green, yellow and orange directions. These conclusions are obtained for stoichiometric Si_{3}N_{4} and SiO_{2} which are not the standard materials in PV industry. For more realistic results, it is mandatory to extend the characterization of the optical properties to long wavelength infrared.
Considering the glass and the EVA layers on front of the solar cell slightly extents the gamut around the green color. As for thermal effect, the heat source is not affected since the absorbed light in the active material is marginally impacted. As for cooling, the thermal resistance of glass and EVA is not increased enough to increase significantly the operating temperature.
As a conclusion, taking into account realistic thermal conditions helps in reducing the difference in terms of conversion efficiency between two solar cells with different colors. This may be very handy for BIPV since making drawings or colored patterns with different color solar cells may not influence significantly the overall conversion efficiency.
Author contribution statement
All authors have contributed to the writing of the manuscript. More precisely, Romain COUDERC and Mohamed AMARA have refined the homemade code of the 1D electrothermooptical model. Mustapha LEMITI and Fabien MANDORLO worked on the optical behavior for any incident signal and color properties of the simulated solar cells. Felix GERENTON has fabricated and characterized all the SiN_{x} and SiO_{x} layers to get realistic estimation of the gamut when simulating the optical stack.
Acknowledgments
The authors are thankful to the NanoLyon platform for processing and characterizing different PECVD deposited layers.
References
 B.P. Jelle, C. Breivik, Energy Procedia 20, 78 (2012) [Google Scholar]
 A.J. McEvoy, L. Castaner, T. Markvart, Solar cells: materials, manufacture and operation (Academic Press, 2012) [Google Scholar]
 E. Skoplaki, J. Palyvos, Sol. Energy 34, 23 (2009) [Google Scholar]
 J. Kurnik, M. Jankovec, K. Brecl, M. Topic, Sol. Energy Mater. Sol. Cells 95, 373 (2011) [Google Scholar]
 I. Tobias, A.E. Moussaoni, A. Luque, IEEE Trans. Electron Devices 46, 1858 (1999) [Google Scholar]
 P. Nubile, Thin Solid Films 342, 257 (1999) [Google Scholar]
 H. Nagel, A.G. Aberle, R. Hezel, Prog. Photovolt.: Res. Appl. 7, 245 (1999) [CrossRef] [Google Scholar]
 B.S. Richards, Prog. Photovolt.: Res. Appl. 12, 253 (2004) [CrossRef] [Google Scholar]
 M. Li, L. Zeng, Y. Chen, L. Zhuang, X. Wang, H. Shen, Int. J. Photoenergy 2013, 352473 (2013) [Google Scholar]
 R. Couderc, M. Amara, M. Lemiti, IEEE J. Photovolt. 6, 1123 (2016) [CrossRef] [Google Scholar]
 International Commission on Illumination. [Online]. Available at: www.cie.co.at [Google Scholar]
 J. Schanda, Colorimetry: understanding the CIE system (John Wiley & Sons, 2007) [Google Scholar]
 International Electrotechnical Commission et al., Multimedia systems and equipment – colour measurement and management – part 21: Colour management – default RGB colour space – sRGB, IEC 6196621, Tech. Rep., 1999 [Google Scholar]
 F. Bernoux, J.P. Piel, B. Castellon, C. Defranoux, J.H. Lecat, P. Boher, J.L. Stehle, Ellipsométrie. Théorie, in Techniques de l’ingénieur. Mesures et contrôle (2003), p. R64901 [Google Scholar]
 D.L. Windt, Comput. Phys. 12, 360 (1998) [CrossRef] [Google Scholar]
 C.A. Gueymard, Sol. Energy 71, 325 (2001) [Google Scholar]
 A.D. Jones, C.P. Underwood, Sol. Energy 70, 349 (2001) [Google Scholar]
 H.R. Philipp, J. Electrochem. Soc. 120, 295 (1973) [Google Scholar]
 J. Kischkat, S. Peters, B. Gruska, M. Semtsiv, M. Chashnikova, M. Klinkmüller, O. Fedosenko, S. Machulik, A. Aleksandrova, G. Monastyrskyi, et al., Appl. Opt. 51, 6789 (2012) [CrossRef] [PubMed] [Google Scholar]
 L. Gao, F. Lemarchand, M. Lequime, J. Eur. Opt. Soc. Rapid Publ. 8, 13010 (2013) [Google Scholar]
 J.F. Lelièvre, E. Fourmond, A. Kaminski, O. Palais, D. Ballutaud, M. Lemiti, Sol. Energy Mater. Sol. Cells 93, 1281 (2009) [Google Scholar]
 F. Gérenton, F. Mandorlo, E. Fourmond, M. Le Coz, D. BlancPélissier, M. Lemiti, J. Vac. Sci. Technol. A 34, 051201 (2016) [Google Scholar]
 D.B.M. Klaassen, SolidState Electron. 35, 953 (1992) [CrossRef] [Google Scholar]
 R. Couderc, M. Amara, M. Lemiti, J. Appl. Phys. 115, 093705 (2014) [Google Scholar]
 R. Pässler, Phys. Rev. B 66, 085201 (2002) [Google Scholar]
 A. Schenk, J. Appl. Phys. 84, 3684 (1998) [Google Scholar]
 J. Henrie, S. Kellis, S. Schultz, A. Hawkins, Opt. Express 12, 1464 (2004) [CrossRef] [PubMed] [Google Scholar]
 R.A. Jarvis, Inf. Process. Lett. 2, 18 (1973) [Google Scholar]
 D.N. Wright, E.S. Marstein, A. Holt, Double layer antireflective coatings for silicon solar cells, in Photovoltaic Specialists Conference, 2005. Conference Record of the Thirtyfirst IEEE (IEEE, 2005), pp. 1237–1240 [CrossRef] [Google Scholar]
 H. Nagel, A.G. Aberle, R. Hezel, Prog. Photovolt.: Res. Appl. 7, 245 (1999) [CrossRef] [Google Scholar]
 A. Royne, C.J. Dey, D.R. Mills, Sol. Energy Mater. Sol. Cells 86, 451 (2005) [Google Scholar]
 G. Notton, C. Cristofari, M. Mattei, P. Poggi, Appl. Therm. Eng. 25, 2854 (2005) [Google Scholar]
 L. Weiss, M. Amara, C. Ménézo, Prog. Photovolt.: Res. Appl 24, 12 (2015) [CrossRef] [Google Scholar]
 E. Kaplani, S. Kaplanis, Sol. Energy 107, 443 (2014) [Google Scholar]
Cite this article as: M. Amara, F. Mandorlo, R. Couderc, F. Gérenton, M. Lemiti, Temperature and color management of silicon solar cells for building integrated photovoltaic, EPJ Photovoltaics 9, 1 (2018)
All Tables
Optimal thickness as a function of the heat coefficient to maximize the conversion efficiency.
All Figures
Fig. 1 CIE 1931 standard observer color matching functions. 

In the text 
Fig. 2 CIE1931 color space chromaticity diagram including the sRGB one (used for displays). 

In the text 
Fig. 3 Description of the simulated solar cell and its boundary conditions. 

In the text 
Fig. 4 Spectral irradiances for different air masses. 

In the text 
Fig. 5 Refractive index of SiO_{2} and Si_{3}N_{4} to cover both the visible and long wavelength infrared domains (data extracted from [18–20]. 

In the text 
Fig. 6 Experimental optical indexes of different silicon oxides and nitrides (refractive index n at λ = 600 nm in brackets). 

In the text 
Fig. 7 Perceived color converted into sRGB (a) and brightness (b) depending on the observation angle and ARC thickness for Si_{3}N_{4}. 

In the text 
Fig. 8 CIE1931 color space chromaticity for a Si_{3}N_{4} ARC with an AM_{1.5G} irradiance. For θ = 0, arrows represent the direction for which t_{ARC} is increasing and round marks correspond to ARC thicknesses multiples of 25 nm. 

In the text 
Fig. 9 Reflected power for different thicknesses at normal incidence. 

In the text 
Fig. 10 Influence of the illumination on the perceived color (Si_{3}N_{4} ARC) at normal incidence for different irradiances. 

In the text 
Fig. 11 CIE 1931 color space chromaticity for different ARC layers with an AM_{1.5G} irradiance (optical index at λ = 600 nm in brackets). The green point in the center of the gamut corresponds the simulation with no ARC. 

In the text 
Fig. 12 Perceived color ARC layer and brightness converted into sRGB for a SiO_{2}/Si_{3}N_{4} double at normal incidence. 

In the text 
Fig. 13 CIE xyY 1931 color space for a doublestack ARC layer (SiO_{2}/Si_{3}N_{4}) with an AM_{1.5G} irradiance. 

In the text 
Fig. 14 Solar cell temperature as a function of the ARC thickness under AM_{1.5G} illumination. For h_{conv} = + ∞, temperature is 25 ^{°}C whatever t_{ARC}. Solid lines: Si_{3}N_{4}, dashed lines: SiO_{2}. 

In the text 
Fig. 15 Heat and sink sources as a function of the ARC thickness (Solid lines: Si_{3}N_{4}, dashed lines: SiO_{2}). 

In the text 
Fig. 16 Reflectivity of the solar cell at the maximum of the thermalization sources. 

In the text 
Fig. 17 Short current circuit current density as a function of the ARC thickness (solid lines: Si_{3}N_{4}, dashed lines: SiO_{2}). 

In the text 
Fig. 18 Absolute and relative conversion efficiencies as a function of the ARC thickness for both SiO_{2} and Si_{3}N_{4}. 

In the text 
Fig. 19 Conversion efficiency and perceived color for different single layer ARC thicknesses (with θ = 0). 

In the text 
Fig. 20 CIE1931 color space chromaticity for different SiN_{x}:H ARC layers with an AM_{1.5G} irradiance (optical index at λ = 600 nm in brackets) when taking the EVA and glass layers into account. The green point in the center of the gamut corresponds the simulation with no ARC. 

In the text 
Current usage metrics show cumulative count of Article Views (fulltext article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.
Data correspond to usage on the plateform after 2015. The current usage metrics is available 4896 hours after online publication and is updated daily on week days.
Initial download of the metrics may take a while.