Temperature and color management of silicon solar cells for building integrated photovoltaic

. Color management of integrated photovoltaics must meet two criteria of performance: provide maximum conversion ef ﬁ ciency 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 in ﬂ uence of the heat transfer and temperature. In this article, we focus on the color space and brightness achieved by varying the antire ﬂ ective properties of ﬂ at 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 ef ﬁ ciency, 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 ﬂ exibility 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 ef ﬁ ciency at the same time.


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 anti-reflective 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 electro-optical 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 anti-reflective 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][7][8][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 thermo-electro-optical properties.

Modeling of solar cells performances and colors 2.1 Solar cell simulation
For cell simulation, we use our in-house 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 ), non-radiative 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] mm, 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 ).

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 on 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 short-middle-long wavelength of "standard" human cone cells, this color space relies on 3 color matching functions xðlÞ, yðlÞ, and zðlÞ (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(l) on the visible range from 380 to 780 nm [12].For this article, the optical signal S(l) will be the reflected light (R(l)) coming from a normalized power density spectrum AM x (l), leading to S(l) = R(l) ⋅ AM x (l): with: 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: 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.

Simulation inputs parameters 3.1 Description of the reference solar cell
A non-encapsulated flat photovoltaic cell is considered with a generic structure composed of a standard p-n junction with a negatively doped upper zone (n-type, emitter), a positively doped lower zone (p-type, 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 u.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 multi-stack 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 .

Optical and thermophysical properties
In this article, we use two types of data: experimental measurements and databases (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 l = 15 mm.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 mm.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 LF-PECVD with different stoichiometries [21,22].One of the SiN x :H layers (n = 2.08 at l = 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, free-carriers 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].

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 (u = 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 u 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 u = 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 u increases.This particular trend is verified whatever the thickness of the ARC.Consequently, drawing the chromaticity for u = 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 dark-blue 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 l = 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.

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 u = 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.

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 sub-stoichiometric 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 ffiffiffiffiffiffiffiffiffiffiffiffiffiffiffi ffi n Si ⋅n air p (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.

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 u = 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 t Si 3 N 4 ARC > 650 nm while for t SiO 2 ARC > 650 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.

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 non-radiative 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   The anti-reflecting layers are usually designed in order to maximize the short-circuit 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 h 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 h; 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 h 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 ).

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 electrothermo-optical 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 u = 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 multi-color 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).

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.

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 ffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffi n Si ⋅n EV A p . 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 green-cyan and green-yellow 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.

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]: Table 2. Optimal thickness as a function of the heat coefficient to maximize the conversion efficiency.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: 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 non-encapsulated cell.The only case which is not treated correctly through the non-encapsulated 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.

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  sub-stoichiometric 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.

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

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

Fig. 8 .
Fig. 8. CIE1931 color space chromaticity for a Si 3 N 4 ARC with an AM 1.5G irradiance.For u = 0, arrows represent the direction for which t ARC is increasing and round marks correspond to ARC thicknesses multiples of 25 nm.

Fig. 7 .
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. 10 .
Fig. 10.Influence of the illumination on the perceived color (Si 3 N 4 ARC) at normal incidence for different irradiances.

Fig. 14 .
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. 12 .
Fig. 12. Perceived color and brightness converted into sRGB for a SiO 2 /Si 3 N 4 double layer ARC at normal incidence.

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

Fig. 17 .
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. 15 .
Fig. 15.Heat and sink sources as a function of the ARC thickness (Solid lines: Si 3 N 4 , dashed lines: SiO 2 ).

Fig. 20 .
Fig. 20.CIE1931 color space chromaticity for different SiN x :H ARC layers with an AM 1.5G irradiance (optical index at l = 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.

Table 1 .
Parameters and their default values used in simulation.
Absolute and relative conversion efficiencies as a function of the ARC thickness for both SiO 2 and Si 3 N 4 .