Effect of the Incoherent Encapsulation Layer and Oblique Sunlight Incidence on the Optical and Current-Voltage Characteristics of Surface-Textured Cu(In,Ga)Se 2 Solar Cells Based on the Angle-Dependent Equispaced Thickness Averaging Method

Sunlight on the Optical and Current-Voltage Characteristics of Surface-Textured Cu(In,Ga)Se 2 Solar Cells Based on the Angle-Dependent Thickness Abstract: In general, the optical and electrical characteristics of Cu(In,Ga)Se 2 (CIGS) solar cells have been studied under the condition that sunlight is normally incident from the air to the CIGS solar cell having no thick front encapsulation layers. To obtain the calculation results in a realistic module application, we calculate the optical and current–voltage (J–V) characteristics of surface-textured CIGS solar cells by simultaneously considering the thick front encapsulation layers and oblique sunlight incidence. Using the proposed angle-dependent equispaced thickness averaging method (ADETAM), we incoherently model two successive front encapsulation layers of a cover glass layer and an ethylene vinyl acetate (EVA) layer, whose respective thicknesses are greater than the coherence length of sunlight (~0.6 µ m). The angular dependences of reﬂectance spectrum and J–V curves are calculated and compared in a surface-textured CIGS solar cell with and without the inclusion of the two front encapsulation layers. We show that the optical absorption improvement of the surface-textured CIGS solar cell over the planar CIGS solar cell can be over-predicted when the thick front encapsulation layers are not considered in the optical modeling.


Introduction
Thin-film solar cells have been greatly studied as one of the promising alternatives to traditional silicon-based solar cells due to their potential to low cost and high efficiency [1]. Thin-film solar cells have been fabricated in various absorption materials such as amorphous silicon, organic compound, and Cu(In, Ga)Se 2 (CIGS) [2][3][4]. Among them, the CIGS solar cell has reached a power conversion efficiency of beyond 20% owing to high optical absorption coefficient, wide absorption spectrum, and long material stability [5]. Furthermore, it has reached a power conversion efficiency of 23.35% with double buffer layers [6]. Generally, thin-film solar cells, including the CIGS solar cell, comprise thin multilayer on the order of sub-micrometer and a thick substrate on the order of one millimeter. Because of the optical interference effect, the absorption efficiency of the thin-film solar cell is greatly affected by the nano-meter-order thickness of the multilayer structure. Thus, optimizing the thickness of each thin layer based on accurate optical modeling is crucial to maximize the optical absorption in the active region and to improve the power conversion efficiency.
In the optical modeling of the thin-film solar cell, the coherence length of sunlight (~0.6 µm) [7], within which the phase information of sunlight is preserved, should be carefully considered. The sub-micrometer multilayer, covering an absorptive active region, a carrier transport layer, and a transparent electrode, can be modeled as optically coherent layers, where their reflectance and absorbance are affected by optical interference effect [3]. In contrast, the millimeter-order thick substrate has to be treated as an optically incoherent layer, within which the phase of light is randomly changed at each coherence length and optical interference effect disappears due to the averaging of optical responses over the random phase changes [8,9].
To calculate the reflectance or absorbance of this mixed coherent-incoherent thinfilm solar cells, various analytical methods have been used for planar structures. In the generalized transfer matrix method (GTMM), the propagation behavior of the electric field intensity in planar incoherent layers was separately described by matrix formalism similar to the transfer matrix method (TMM) used in coherent planar multilayers [10,11]. In the case of the random phase method (RPM), coherent calculation results obtained by the TMM were averaged over thousands of random phase shifts added to a selected incoherent layer [12]. To reduce the computation time, coherent simulation results with respect to a set of equidistant phase shifts added to an incoherent layer were averaged in the equidistant phase method (EPM) [13]. Using the GTMM, our group successfully calculated and analyzed the optical absorption characteristics of a planar thin-film organic solar cell having an incoherent thick glass substrate in both normal and oblique incidence of sunlight [14,15].
Light absorption efficiency of planar thin-film solar cells can be improved by use of surface-textured structures, of which absorption characteristics cannot be calculated based on the above-mentioned analytical methods applied to planar structures. In this case, the numerical methods of the rigorous coupled wave analysis (RCWA) [16], the finite-difference time-domain (FDTD) [17,18], and the finite element method (FEM) [19][20][21] should be used, where wave equations are directly solved in a small mesh or grid for optical modeling. Because these numerical methods provide only coherent simulation results, an additional averaging procedure such as the so-called "spectral averaging method" (SAM) [9], the first-principle calculation [22], and "one-pass coherent calculation" [23] is required to include the incoherent characteristics of light propagation. However, these methods have very complicated mathematical and computational procedures, which result in a long computation time and implementation complexity [21]. Recently, we proposed a so-called "equispaced thickness averaging method" (ETAM) as a simple and efficient optical modeling method to calculate the incoherent characteristic of surface-textured thinfilm solar cell [20]. In the numerical simulation of the ETAM, a millimeter-order incoherent layer could be converted into a micrometer-order coherent layer with a set of equispaced phase thicknesses and corresponding coherent simulation results were averaged. Using the ETAM, we successfully calculated the effect of the three successive incoherent layers on the absorption characteristics of a surface-textured CIGS solar cell in only normal incidence of sunlight [21].
The above-mentioned numerical modeling methods to consider the mixed coherentincoherent characteristics of surface-textured thin-film solar cells have been performed under the condition that sunlight is normally incident from the air, which did not reflect a realistic deployment situation of thin-film solar cells. In a general deployment situation such as a rooftop application, the incidence angle of sunlight gradually varies during the daytime such that an optimized device structure of the surface-textured solar cell has to be determined in the consideration of the variation in the incidence angle and polarization of sunlight [15]. Moreover, in a realistic module configuration, millimeter thick front encapsulant and cover glass layers, which can provide a higher mechanical durability and operational stability, are deposited on the top contact layer. These front optically-incoherent encapsulation layers will affect the light absorption and power conversion efficiencies of surface-textured thin-film solar cell, which also has to be considered to obtain the accurate optical modeling results. However, the effects of both incoherent encapsulation layers and oblique sunlight incidence have not yet been simultaneously considered on the optical modeling of surface-textured CIGS solar cells.
Appl. Sci. 2021, 11, 2121 3 of 18 In this paper, we calculate the optical and current-voltage (J-V) characteristics of surface-textured CIGS solar cells by simultaneously considering the front incoherent encapsulation layers and oblique sunlight incidence. We propose an angle-dependent equispaced thickness averaging method (ADETAM) that can calculate the incident-angledependent optical characteristics of the mixed coherent-incoherent solar cell by modifying the currently-used ETAM only applicable to the normal sunlight incidence. The layer structure of the CIGS solar cell includes two successive front encapsulation layers of a cover glass and an ethylene vinyl acetate (EVA), whose respective thicknesses are greater than the coherence length of sunlight (~0.6 µm). We verify the calculation accuracy of the ADETAM in a planar CIGS solar cell through a comparison with exact analytical results obtained by the GTMM. In a surface-textured CIGS solar cell, the reflectance spectra and J-V characteristics depending on the incidence angle and polarization of sunlight are calculated and compared when two successive encapsulation layers are included or not in the simulation. According to the calculation results, the output performance enhancement of the surface-textured CIGS solar cell over the planar CIGS solar cell can be over-predicted when the thick encapsulation layers are not considered in the optical modeling.

Theory of ADETAM
The theoretical formulation of the ADETAM, which is an averaging procedure for coherent simulation results to consider the incoherent characteristics in a mixed coherentincoherent multiplayer, can be applied to both planar and non-planar structures. For convenience, we derive the theoretical formulation of the ADETAM in a planar mixed coherent-incoherent thin-film multilayer structure shown in Figure 1. of surface-textured thin-film solar cell, which also has to be considered to obtain the accurate optical modeling results. However, the effects of both incoherent encapsulation layers and oblique sunlight incidence have not yet been simultaneously considered on the optical modeling of surface-textured CIGS solar cells.
In this paper, we calculate the optical and current-voltage (J-V) characteristics of surface-textured CIGS solar cells by simultaneously considering the front incoherent encapsulation layers and oblique sunlight incidence. We propose an angle-dependent equispaced thickness averaging method (ADETAM) that can calculate the incident-angledependent optical characteristics of the mixed coherent-incoherent solar cell by modifying the currently-used ETAM only applicable to the normal sunlight incidence. The layer structure of the CIGS solar cell includes two successive front encapsulation layers of a cover glass and an ethylene vinyl acetate (EVA), whose respective thicknesses are greater than the coherence length of sunlight (∼0.6 μm). We verify the calculation accuracy of the ADETAM in a planar CIGS solar cell through a comparison with exact analytical results obtained by the GTMM. In a surface-textured CIGS solar cell, the reflectance spectra and J-V characteristics depending on the incidence angle and polarization of sunlight are calculated and compared when two successive encapsulation layers are included or not in the simulation. According to the calculation results, the output performance enhancement of the surface-textured CIGS solar cell over the planar CIGS solar cell can be over-predicted when the thick encapsulation layers are not considered in the optical modeling.

Theory of ADETAM
The theoretical formulation of the ADETAM, which is an averaging procedure for coherent simulation results to consider the incoherent characteristics in a mixed coherentincoherent multiplayer, can be applied to both planar and non-planar structures. For convenience, we derive the theoretical formulation of the ADETAM in a planar mixed coherent-incoherent thin-film multilayer structure shown in Figure 1. The multiple incoherent layers (j = 1, ···, m) are assumed to be transparent and the remaining coherent layers (j = m+1, ···, n) can be either transparent or absorptive. To optically model random phase variations in the incoherent layers, each incoherent layer (v = 1, ···, m) is divided into an initial layer of the original thickness dv and an additional phase-shift layer of the equispaced phase thickness bv. A plane wave is obliquely incident from the semi-infinite transparent ambient on the left (j = 0) with the incidence angle θ0. The directions of the s-and p-polarized light are orthogonal to or parallel to the plane of incidence.
We assume that multiple incoherent layers, representing millimeter-order front encapsulation layers, lie in front of the coherent multilayer structure, which can range from Figure 1. Schematic diagram of a planar mixed coherent-incoherent thin-film multilayer structure. The multiple incoherent layers (j = 1, ···, m) are assumed to be transparent and the remaining coherent layers (j = m + 1, ···, n) can be either transparent or absorptive. To optically model random phase variations in the incoherent layers, each incoherent layer (v = 1, ···, m) is divided into an initial layer of the original thickness d v and an additional phase-shift layer of the equispaced phase thickness b v . A plane wave is obliquely incident from the semi-infinite transparent ambient on the left (j = 0) with the incidence angle θ 0 . The directions of the s-and p-polarized light are orthogonal to or parallel to the plane of incidence.
We assume that multiple incoherent layers, representing millimeter-order front encapsulation layers, lie in front of the coherent multilayer structure, which can range from a transparent conducting oxide (TCO) top contact to a metal back contact. The mixed coherent-incoherent multilayer is assumed to be sandwiched between a semi-infinite transparent ambient medium on the left (j = 0) and the right (j = n + 1). The multiple incoherent layers (j = 1, ···, m) are assumed to be transparent and the remaining coherent layers (j = m + 1, ···, n) can be either transparent and absorptive. All the layers are assumed to be Appl. Sci. 2021, 11, 2121 4 of 18 isotropic and homogeneous with planar and parallel interfaces. When the thickness of each layer is d j , the corresponding complex refractive index is denoted by n j = n j + iκ j , where n j indicates the refractive index and κ j represents the extinction coefficient, respectively. To optically model random phase variations in the incoherent layers, each incoherent layer (v = 1, ···, m) is divided into an initial layer of the original thickness d v and an additional phase-shift layer of the equispaced phase thickness b v [21].
We assume that a plane wave, which has the optical intensity of I 0 at the wavelength of λ, is obliquely incident from the semi-infinite transparent ambient on the left (j = 0) with the incidence angle θ 0 . According to the Snell's law, the refraction angle of each layer is determined by sin θ j = (n 0 / n j ) sin θ 0 and k 0 = 2π/λ is the magnitude of the wave number in the ambient (j = 0). When the light propagating from the left to the right denotes a positive direction, the light propagation in the positive and negative directions are expressed as + and − superscripts. The directions of the s-and p-polarized light are orthogonal to or parallel to the plane of incidence.
where β j = k 0 n j cos θ j is the propagation constant at the layer j. For the multiple incoherent layers (v = 1, · · · , m), light propagation in each additional phase-shift layer can be expressed as where represents a random phase shift that light experience in the incoherent layer v. If the number of equispaced phase shifts for each incoherent layer is set to be X, the equispaced phase thickness b v can be determined by Because of the assumption that the incoherent layers are transparent, the values of the refractive index n v and refraction angle θ v are real, which results in a real-valued equispaced phase thickness b v . In addition, the equispaced phase thickness b v are the same for both polarizations.
In Figure 1, the light propagation through the incoherent layers are written as where S s(p) inc (θ q1 1 , θ q2 2 , · · ·, θ q m m ) indicates the system transfer matrix for the whole incoherent multilayer, of which the matrix elements depend on a combination of the equispaced phase shifts of θ q 1 1 (q 1 = 1, · · · , X),···, and θ q m m (q m = 1, · · · , X). In the remaining coherent multiple layers (j = m + 1, ···, n), the electric field amplitude propagating from the righthand boundary of the rightmost incoherent layer (j = m) to the left-hand boundary of the layer j is described by Based on the matrix elements in Equation (7), the front transmission coefficient propagating from the layer 0 to the layer j is determined by .
Similarly, the back reflection coefficient propagating from the layer j to the layer 0 is given by When we focus on the absorption characteristics in the layer j within the coherent multilayer, the electric field intensity at the position z (0 ≤ z ≤ d j ) can expressed as [4] Here, r s(p) j(n+1) is the reflection coefficient propagating from the coherent layer j to the right-hand ambient layer n + 1, which can be determined similarly to Equation (9) except the absence of the equispaced phase shifts. Correspondingly, the spatial distribution of light absorption within the layer j is calculated as [4,15] Q s(p) j (z; θ where c is the speed of light in free space, ε 0 is the electric permittivity in free space, and α j = 4πκ j /λ is the absorption coefficient of the layer j. In the ADETAM, the incoherent characteristics of light propagation within the incoherent multilayer (v = 1, ···, m) can be considered by averaging coherent calculation results over a set of equispaced phase thicknesses. The ADETAM-applied electric field intensity is given by where X m is the number of total coherent calculations required to include the incoherent effect of m multiple incoherent layers. Hence, the ADETAM-applied spatial distribution of light absorption within the layer j is expressed as Then, the ADETAM-applied spectral absorbance at the layer j is where I s(p) 0 is the optical power of the s(p)-polarized incident light at the wavelength of λ. Finally, the ADETAM-applied total reflectance in Figure 1 can be calculated as where we assume that the incident sunlight is totally absorbed in the active region or metal electrode without transmitting to the ambient medium on the right (j = n + 1), which is valid if the metal back contact layer is sufficiently thick Figure 2a shows a schematic diagram of the planar CIGS solar cell used in this calculation [24]. The sunlight with the incidence angle of θ 0 is incident from the air to the CIGS solar cell, which can be considered as a mixed incoherent-coherent system. The spectral irradiance of sunlight is assumed to follow the air mass (AM) 1.5 sunlight with a total power of S 0 = 100 mW/cm 2 , which is shown in inset of Figure 2a. Because sunlight is unpolarized, the incidence optical powers of s-and p-polarized light are equally distributed at each incidence angle of θ 0 . The CIGS solar cell is composed of a 3 mm cover glass at the top, a 0.5 mm EVA layer as a surface flattening layer, 50 nm aluminum-doped zinc oxide (AZO) as a TCO, 10-nm zinc sulfide (ZnS) as an n-doped layer, 2 µm CIGS as a p-doped active layer, and 400 nm molybdenum (Mo) as a bottom contact layer [24]. The bottom glass substrate and EVA layer are not included in this calculation because they have no effect on the optical and electrical characteristics of the CIGS solar cell. Figure 2b shows the complex refractive index spectra of the materials used in the CIGS solar cell, which are taken from the literature (CIGS, ZnS, AZO, and Mo) [25][26][27]. It is well known that the complex refractive index of semiconductor materials depends on doping concentration. In particular, the variation of Al dopant concentration significantly changes the complex refractive index of AZO thin films [28,29]. In this study, the donor concentration of the ZnS and acceptor concentration of the CIGS are set to be 5 × 10 16 cm −3 and 3 × 10 16 cm −3 , respectively [24]. According to the complex refractive index spectra of the AZO shown in Figure 2b, the Al concentration of the AZO is speculated to be 2 wt% in reference to the measured complex refractive index spectra of the AZO in [28]. that the complex refractive index of semiconductor materials depends on doping concen-tration. In particular, the variation of Al dopant concentration significantly changes the complex refractive index of AZO thin films [28,29]. In this study, the donor concentration of the ZnS and acceptor concentration of the CIGS are set to be 5 × 10 16 cm -3 and 3 × 10 16 cm -3 , respectively [24]. According to the complex refractive index spectra of the AZO shown in Figure 2(b), the Al concentration of the AZO is speculated to be 2 wt% in reference to the measured complex refractive index spectra of the AZO in [28].

Planar CIGS Solar Cells
(a) Schematic diagram of the planar CIGS solar cell together with the spectral irradiance of the AM 1.5 sunlight ranging between 300 and 1100 nm. Both the 3 mm cover glass and 0.5 mm EVA layers, the thicknesses of which are greater than the coherence length of sunlight (~0.6 μm), are treated as incoherent. The remaining layers are considered as coherent. Because the glass cover and the EVA flattening layers are assumed to be transparent, both layers are modeled as the 1 μm thick initial layer plus the respective phase-shift layer whose thicknesses b1 and b2 are determined based on Eq. (4) when the ADETAM is applied to calculate the effect of the incoherent front glass and EVA layers in the simulation. (b) Complex refractive index spectra of the materials used in the CIGS solar cell.
In this calculation, the front encapsulation layers of the cover glass and EVA are assumed to be transparent at the wavelength range between 300 and 1100 nm. In general, the cover glass and EVA has high absorption at the wavelength range between 300 and 350 nm [30,31]. However, our assumption that the cover glass and EVA are transparent between 300 and 350 nm will not significantly change the calculation results of optical absorption at the ZnS/CIGS active layers, where the generated electron-hole pairs contribute to the short-circuit current density (JSC). The first reason is that the AZO, of which complex refractive index spectrum is shown in Figure 2, also has very high extinction coefficient at the wavelength range between 300 and 350 nm. Hence, most of incident sunlight within these wavelength ranges is absorbed at the AZO TCO, which results in very low external quantum efficiency at the wavelength range between 300 and 350 nm [24]. Secondly, the amount of the AM1.5 sunlight irradiance at the wavelength range between 300 and 350 nm, shown in the inset of Figure 2(a), is relatively lower than that of the AM1.5 sunlight irradiance at the other wavelengths, which indicates that the contribution of JSC generated at the wavelength range between 300 and 350 nm is relatively small.
Both the 3 mm cover glass and 0.5 mm EVA layers, the thicknesses of which are greater than the coherence length of sunlight (~0.6 μm), are treated as incoherent. Because Figure 2. (a) Schematic diagram of the planar CIGS solar cell together with the spectral irradiance of the AM 1.5 sunlight ranging between 300 and 1100 nm. Both the 3 mm cover glass and 0.5 mm EVA layers, the thicknesses of which are greater than the coherence length of sunlight (~0.6 µm), are treated as incoherent. The remaining layers are considered as coherent. Because the glass cover and the EVA flattening layers are assumed to be transparent, both layers are modeled as the 1 µm thick initial layer plus the respective phase-shift layer whose thicknesses b 1 and b 2 are determined based on Equation (4) when the ADETAM is applied to calculate the effect of the incoherent front glass and EVA layers in the simulation. (b) Complex refractive index spectra of the materials used in the CIGS solar cell.
In this calculation, the front encapsulation layers of the cover glass and EVA are assumed to be transparent at the wavelength range between 300 and 1100 nm. In general, the cover glass and EVA has high absorption at the wavelength range between 300 and 350 nm [30,31]. However, our assumption that the cover glass and EVA are transparent between 300 and 350 nm will not significantly change the calculation results of optical absorption at the ZnS/CIGS active layers, where the generated electron-hole pairs contribute to the short-circuit current density (J SC ). The first reason is that the AZO, of which complex refractive index spectrum is shown in Figure 2, also has very high extinction coefficient at the wavelength range between 300 and 350 nm. Hence, most of incident sunlight within these wavelength ranges is absorbed at the AZO TCO, which results in very low external quantum efficiency at the wavelength range between 300 and 350 nm [24]. Secondly, the amount of the AM1.5 sunlight irradiance at the wavelength range between 300 and 350 nm, shown in the inset of Figure 2a, is relatively lower than that of the AM1.5 sunlight irradiance at the other wavelengths, which indicates that the contribution of J SC generated at the wavelength range between 300 and 350 nm is relatively small.
Both the 3 mm cover glass and 0.5 mm EVA layers, the thicknesses of which are greater than the coherence length of sunlight (~0.6 µm), are treated as incoherent. Because the cover glass and the EVA flattening layers are assumed to be transparent, both layers are modeled as the 1 µm thick initial layer plus the respective phase-shift layer whose thicknesses b it can be treated as a coherent layer because no optical interference effect is observed in the light absorption profiles of the CIGS layer.The optical modeling of the planar CIGS solar cell is numerically performed based on the two-dimensional (2D) FEM, where the RF module of COMSOL Multiphysics is used as a simulation platform [32]. The wave equation is numerically solved at very tiny meshes of the simulation domains with appropriate boundary conditions. We use the port boundary condition at the uppermost boundary of the air domain, where a plane wave having the wavelength of λ and the incidence angle of θ 0 , is generated for either s and p polarization. The wavelength-dependent input power of the incident plane wave for each polarization is determined by (1/2)S(λ)dλ, where S(λ) indicates the spectral irradiance of the AM 1.5 sunlight shown in Figure 2a. Because the AM 1.5 sunlight is unpolarized, a factor of 1/2 is included to represent the incident optical power for each polarization. In this simulation, the spectral irradiance of the AM 1.5 sunlight between 300 and 1100 nm is divided with the wavelength spacing of dλ = 10 nm. The port boundary condition applied to the lowermost boundary suppresses the artificial reflection at the bottom boundary. In addition, the Floquet periodic boundary conditions are used at the leftmost and rightmost boundaries of the simulation domain.
To show how the light absorption takes place within the planar CIGS solar cell, two FEM simulation results, calculated at the single wavelength of 400 nm, are represented in Figure 3. Here, spectral 2D coherent calculation results of the electric field amplitude and optical absorption are obtained at the wavelength of 400 nm when the equispaced phase thicknesses in the phase-shift layers are zero (b 1 1 = b 1 2 = 0). The 400 nm input plane wave with the surface power density of 5.57 W/m 2 has the incidence angle of 45 • for both polarizations. The 2D spatial profiles of the electric field amplitude have a tilted wavy pattern, which is ascribed to the interference between the obliquely incident and reflected lights. On the other hand, the spatial profiles of optical absorption depend only on the longitudinal z-direction and have no spatial dependence in the transverse x-direction, which results from the fact that optical power only flows into the longitudinal z direction.
are modeled as the 1 μm thick initial layer plus the respective phase-shift layer whose thicknesses 1 1 q b and 2 2 q b are determined based on Eq. (4) when the ADETAM is applied to calculate the effect of the incoherent front glass and EVA layers in the simulation. The remaining layers from the AZO TCO to the Mo bottom contact are considered as coherent. As will be shown in Figure 3, most of light are absorbed near the ZnS/CIGS interface and a part of light reflected from the Mo bottom contact layer is negligible. Although the 2 μm thickness of the CIGS layer is greater than the coherence length of sunlight (~0.6 μm), it can be treated as a coherent layer because no optical interference effect is observed in the light absorption profiles of the CIGS layer.The optical modeling of the planar CIGS solar cell is numerically performed based on the two-dimensional (2D) FEM, where the RF module of COMSOL Multiphysics is used as a simulation platform [32]. The wave equation is numerically solved at very tiny meshes of the simulation domains with appropriate boundary conditions. We use the port boundary condition at the uppermost boundary of the air domain, where a plane wave having the wavelength of λ and the incidence angle of θ0, is generated for either s and p polarization. The wavelength-dependent input power of the incident plane wave for each polarization is determined by (1/2)S(λ)dλ, where S(λ) indicates the spectral irradiance of the AM 1.5 sunlight shown in Figure 2(a). Because the AM 1.5 sunlight is unpolarized, a factor of 1/2 is included to represent the incident optical power for each polarization. In this simulation, the spectral irradiance of the AM 1.5 sunlight between 300 and 1100 nm is divided with the wavelength spacing of dλ = 10 nm. The port boundary condition applied to the lowermost boundary suppresses the artificial reflection at the bottom boundary. In addition, the Floquet periodic boundary conditions are used at the leftmost and rightmost boundaries of the simulation domain.
To show how the light absorption takes place within the planar CIGS solar cell, two FEM simulation results, calculated at the single wavelength of 400 nm, are represented in Figure 3. Here, spectral 2D coherent calculation results of the electric field amplitude and optical absorption are obtained at the wavelength of 400 nm when the equispaced phase thicknesses in the phase-shift layers are zero ( = = ). The 400 nm input plane wave with the surface power density of 5.57 W/m 2 has the incidence angle of 45° for both polarizations. The 2D spatial profiles of the electric field amplitude have a tilted wavy pattern, which is ascribed to the interference between the obliquely incident and reflected lights.
On the other hand, the spatial profiles of optical absorption depend only on the longitudinal z-direction and have no spatial dependence in the transverse x-direction, which results from the fact that optical power only flows into the longitudinal z direction.  When the ADETAM is applied to the incoherent modeling of a mixed coherentincoherent multilayer, the accuracy of the calculation results converges if the number of the equispaced phase thicknesses is greater than five [20,21]. In this calculation, we set the number of equispaced phase thicknesses in the phase-shift layer to be ten. According to Equation (4), the equispaced phase thickness of the phase-shift layer at the wavelength of λ can be b 1 , · · ·, b 9 1 = 8λ 20n g cos θ 1 , b 10 1 = 9λ 20n g cos θ 1 , Here, n g and θ 1 are the refractive index and refraction angle of the cover glass while n EVA and θ 2 are those of the EVA layer. At each combination of equispaced phase thicknesses in the two incoherent layers, FEM-based numerical solutions of the optical absorption profiles such as those in Figure 4 are obtained at every mesh point. Then, the ADETAM-applied optical absorption in Equation (13) is obtained by averaging a total of 10 2 = 100 different coherently calculated electric field intensities at each mesh point.
To verify the validity of the calculation results based on the proposed ADETAM, the ADETAM-applied calculation results are compared with those obtained by the GTMM. Figure 4 shows the calculated spectral spatial distributions of light absorption obtained by the ADETAM and the GTMM when the input plane wave with the wavelength of λ = 400 nm and the surface power density of 5.57 W/m 2 has the incidence angle of θ 0 = 0 • , 30 • , 45 • , and 60 • for both polarizations. The calculated spectral spatial absorption profiles obtained by the ADETAM are well matched with the exact analytical results calculated by the GTMM, which includes the optically incoherent characteristics of the cover glass and EVA layers. Figure 5 shows a comparison of the calculated reflectance spectra of the planar CIGS solar cell between the ADETAM and the GTMM at various incidence angles for sand p-polarized light. The reflectance spectra calculated by the ADETAM are well matched with those obtained by GTMM for all the incidence angles and polarizations of sunlight.
Here, ng and θ1 are the refractive index and refraction angle of the cover glass while nEVA and θ2 are those of the EVA layer. At each combination of equispaced phase thicknesses in the two incoherent layers, FEM-based numerical solutions of the optical absorption profiles such as those in Figure 4 are obtained at every mesh point. Then, the ADETAM-applied optical absorption in Eq. (13) is obtained by averaging a total of 10 2 = 100 different coherently calculated electric field intensities at each mesh point.
To verify the validity of the calculation results based on the proposed ADETAM, the ADETAM-applied calculation results are compared with those obtained by the GTMM. Figure 4 shows the calculated spectral spatial distributions of light absorption obtained by the ADETAM and the GTMM when the input plane wave with the wavelength of λ = 400 nm and the surface power density of 5.57 W/m 2 has the incidence angle of θ0 = 0°, 30°, 45°, and 60° for both polarizations. The calculated spectral spatial absorption profiles obtained by the ADETAM are well matched with the exact analytical results calculated by the GTMM, which includes the optically incoherent characteristics of the cover glass and EVA layers. Figure 5 shows a comparison of the calculated reflectance spectra of the planar CIGS solar cell between the ADETAM and the GTMM at various incidence angles for s-and p-polarized light. The reflectance spectra calculated by the ADETAM are well matched with those obtained by GTMM for all the incidence angles and polarizations of sunlight.      Figure 6a shows a schematic diagram of the surface-textured CIGS solar cell including the front incoherent encapsulation of the cover glass/EVA layers. The interfaces of the EVA/AZO/ZnS/CIGS multilayer are assumed to be randomly textured with the rootmean-square (RMS) surface roughness of σ RMS , which enhances light absorption caused by the scattering effect. The thickness of each layer, along with its complex refractive index, is the same as that of the planar CIGS solar cell in Figure 2. Figure 6b shows incoherent multiple optical reflections inside the planar cover glass layer. A 2D FEM simulation is performed in the surface-textured CIGS solar cell with the same boundary conditions used for the planar CIGS solar cell. To investigate how the degree of the surface roughness affects the absorption characteristics of the surface-textured CIGIS solar cell, calculation results are obtained at four values of σ RMS = 25, 50, 75, and 100 nm.

Surface-Textured CIGS Solar Cells
the ADETAM and the GTMM at various incidence angles for (a) s-and (b) p-polarized lig calculated reflectance spectra obtained by the ADETAM are well matched with the exact a results calculated by the GTMM. Figure 6(a) shows a schematic diagram of the surface-textured CIGS solar cell ing the front incoherent encapsulation of the cover glass/EVA layers. The interfac EVA/AZO/ZnS/CIGS multilayer are assumed to be randomly textured with th mean-square (RMS) surface roughness of σRMS, which enhances light absorption by the scattering effect. The thickness of each layer, along with its complex refra dex, is the same as that of the planar CIGS solar cell in Figure 2.  transparent, they are modelled as the 1 μm initial layer plus the respective phase-shift layer whose thicknesses b1 and b2 are determine on Eq. (4) in the ADETAM-based calculation. The interfaces of the EVA/AZO/ZnS/CIGS m layer are assumed to be randomly textured with the RMS surface roughness, which varies = 25, 50, 75, and 100 nm. (b) Incoherent multiple optical reflections inside the planar incoh cover glass layer to determine the minimum reflectance of the surface-textured CIGS solar The reflectance from the cover glass to the air is Rg-a and the reflectance from the cover glas EVA is Rg-EVA. The refraction angle at the cover glass is θ1 when the incidence angle of sunl θ0 in the air. Figure 6. (a) Schematic diagram of the surface-textured CIGS solar cell. Because the 3 mm cover glass and the 0.5 mm EVA flattening layers are transparent, they are modelled as the 1 µm thick initial layer plus the respective phase-shift layer whose thicknesses b 1 and b 2 are determined based on Equation (4) in the ADETAM-based calculation. The interfaces of the EVA/AZO/ZnS/CIGS multilayer are assumed to be randomly textured with the RMS surface roughness, which varies as σ RMS = 25, 50, 75, and 100 nm. (b) Incoherent multiple optical reflections inside the planar incoherent cover glass layer to determine the minimum reflectance of the surface-textured CIGS solar cell. The reflectance from the cover glass to the air is R g-a and the reflectance from the cover glass to the EVA is R g-EVA . The refraction angle at the cover glass is θ 1 when the incidence angle of sunlight is θ 0 in the air.

Surface-textured CIGS solar cells
Because the 3 mm cover glass and the 0.5 mm EVA flattening layers are transparent, these incoherent layers are modelled as the 1 µm thick initial layer plus the respective phase-shift layer whose thicknesses b 1 and b 2 are determined based on Equation (4) in the ADETAM-based calculation. To show how the light absorption takes place within the surface-textured CIGS solar cell, two FEM simulation results, calculated at the single wavelength of 400 nm, are represented in Figure 7. Here, spectral 2D coherent calculation results of the optical absorption in the surface-textured CIGIS solar cell are obtained at the wavelength of 400 nm when the equispaced phase thicknesses in the phase-shift layers are zero. The RMS surface roughness is set to be σ RMS = 100 nm. The 400 nm input plane wave with the surface power density of 5.57 W/m 2 has the incidence angle of 45 • for both s and p polarizations. To clarify the spectral optical absorption profile near the ZnS/CIGS interface, all the calculation results are shown near the ZnS/CIGS interface. The absorption distributions of the surface-textured CIGS solar cell are not uniform due to the optical scattering, which is different from those of the planar CIGS solar cell in Figure 4. At each combination of equispaced phase thicknesses, FEM-based coherent numerical results of the optical absorption profiles such as those in Figure 7 are calculated at each combination of the wavelength, polarization, and RMS surface roughness. In the ADETAM-based calculation, a total of 10 2 = 100 optical absorption profiles are averaged at each mesh point. In addition, at each RMS surface roughness, ADETAM-based calculation results are averaged over three different arbitrarily-textured surfaces with the same RMS roughness to simulate CIGS solar cells with random surface-texture.
wavelength of 400 nm when the equispaced phase thicknesses in the phase-shift layers are zero. The RMS surface roughness is set to be σRMS = 100 nm. The 400 nm input plane wave with the surface power density of 5.57 W/m 2 has the incidence angle of 45° for both s and p polarizations. To clarify the spectral optical absorption profile near the ZnS/CIGS interface, all the calculation results are shown near the ZnS/CIGS interface. The absorption distributions of the surface-textured CIGS solar cell are not uniform due to the optica scattering, which is different from those of the planar CIGS solar cell in Figure 4. At each combination of equispaced phase thicknesses, FEM-based coherent numerical results o the optical absorption profiles such as those in Figure 7 are calculated at each combination of the wavelength, polarization, and RMS surface roughness. In the ADETAM-based cal culation, a total of 10 2 = 100 optical absorption profiles are averaged at each mesh point In addition, at each RMS surface roughness, ADETAM-based calculation results are aver aged over three different arbitrarily-textured surfaces with the same RMS roughness to simulate CIGS solar cells with random surface-texture.  Figure 8 shows the calculated ADETAM-applied reflectance spectra of the surface textured CIGS solar cell at various incidence angles of sunlight when the RMS surface roughness is σRMS = 25, 50, 75, and 100 nm. Because we assume that unpolarized sunligh has the equal contribution of s and p polarizations, the total reflectance spectrum at each incidence angle are obtained by the average over the reflectance spectrum of s-polarized light and that of p-polarized light. For comparison, the reflectance spectra of the planar CIGS solar cell are also plotted at each incidence angle of sunlight. In Figure 8, the spectra responses of the reflectance in the surface-textured CIGS solar cell are lower than those in the planar CIGS solar cell for all the incidence angles, which verifies the enhancement o light absorption in the surface-textured CIGS solar cell. As the RMS surface roughness increases, the total reflectance becomes lower, which agrees with the experimental results in the literature [33].  Figure 8 shows the calculated ADETAM-applied reflectance spectra of the surfacetextured CIGS solar cell at various incidence angles of sunlight when the RMS surface roughness is σ RMS = 25, 50, 75, and 100 nm. Because we assume that unpolarized sunlight has the equal contribution of s and p polarizations, the total reflectance spectrum at each incidence angle are obtained by the average over the reflectance spectrum of s-polarized light and that of p-polarized light. For comparison, the reflectance spectra of the planar CIGS solar cell are also plotted at each incidence angle of sunlight. In Figure 8, the spectral responses of the reflectance in the surface-textured CIGS solar cell are lower than those in the planar CIGS solar cell for all the incidence angles, which verifies the enhancement of light absorption in the surface-textured CIGS solar cell. As the RMS surface roughness increases, the total reflectance becomes lower, which agrees with the experimental results in the literature [33].
In Figure 8, as the RMS surface roughness of the surface-textured CIGS solar cell increases, the reflectance of the surface-textured CIGS solar cell decreases and converges to a certain value, which results from the incoherent multiple optical reflections inside the planar cover glass layer, as shown in Figure 6b. The scattering effect of the surface-textured structure cannot affect the incoherent multiple reflection caused by the front encapsulation layers, which will limit the absorption enhancement of the surface-textured CIGS solar cell. Because we assume two front encapsulation layers to be transparent, the reflectance caused by the incoherent multiple optical reflections can be calculated by [34] Here, R g-a is the reflectance from the cover glass to the air and R g-EVA is the reflectance from the cover glass to the EVA. The values of R g-a and R g-EVA are determined by the Fresnel equation, which depends on the refraction angle and light polarization [15]. The black dotted lines in Figure 8 show the calculated values of the incidence-angle-dependent reflectance of the incoherent multiple optical reflections, which are well matched with the minimum reflectance values of the surface-textured CIGS solar cell with σ RMS = 100 nm.
Appl. Sci. 2021, 11, x FOR PEER REVIEW In Figure 8, as the RMS surface roughness of the surface-textured CIGS solar creases, the reflectance of the surface-textured CIGS solar cell decreases and conve a certain value, which results from the incoherent multiple optical reflections ins planar cover glass layer, as shown in Figure 6(b). The scattering effect of the surfa tured structure cannot affect the incoherent multiple reflection caused by the front sulation layers, which will limit the absorption enhancement of the surface-texture solar cell. Because we assume two front encapsulation layers to be transparent, the tance caused by the incoherent multiple optical reflections can be calculated by [34 Here, Rg-a is the reflectance from the cover glass to the air and Rg-EVA is the refle from the cover glass to the EVA. The values of Rg-a and Rg-EVA are determined by the equation, which depends on the refraction angle and light polarization [15]. Th dotted lines in Figure 8 show the calculated values of the incidence-angle-depend flectance of the incoherent multiple optical reflections, which are well matched w minimum reflectance values of the surface-textured CIGS solar cell with σRMS = 100

Electrical Modeling Procedure
The J-V characteristics of the CIGS solar cell are calculated based on the 2D numerical solution of the drift-diffusion model, which is expressed as [24] J n ( where q is the basic electric charge. In the current equation of Equations (19) and (20), J n ( → r ) and J p ( → r ) are electron and hole current densities when µ n(p) and D n(p) represent electron (hole) mobility and diffusion coefficient, respectively. The terms n(  (22) and (23), G( → r ) represents the generation rate of the electron-hole pair and can be obtained by where h is the Plank constant, and c is the speed of light in a free space. The ADETAMapplied spatial distribution of light absorption for s(p)-polarized light is Q s(p) (λ, → r ), whose values are averaged to represent unpolarized light. Finally, R( → r ) represents the recombination rate of the electron-hole pair and is given by where R SRH ( → r ) and R rad ( → r ) are the non-radiative Shockley-Read-Hall (SRH) recombination and radiative recombination rates [33].
The 2D drift-diffusion model, shown in Equations (19)- (25), are numerically solved through the semiconductor module of the COMSOL Multiphysics [32]. In the 2D simulation domains of the planar CIGS solar cell in Figure 2a and surface-textured CIGS solar cell in Figure 6a, both the top anode terminal of the AZO/ZnS interface and the bottom cathode terminal of the Mo/CIGS interface are assumed to be ideal Schottky contacts. The J-V characteristics are calculated at only the regions containing the ZnS and CIGS layer. To minimize the error of a mismatch, the electrical modeling is performed under the same mesh configuration as the optical modeling. The periodic conditions are applied at the leftmost and rightmost boundaries of the simulation domain. The simulation parameters used in the electrical modeling are obtained in [33].

Comparison of the J-V Characteristics without and with the Inclusion of the Incoherent Front Encapsulation Layers
To investigate the effect of the front incoherent encapsulation layers on output performance of the surface-textured CIGS solar cell, the angular dependences of reflectance spectrum and J-V curves are calculated and compared between the planar and the surfacetextured CIGS solar cells without and with the inclusion of the two front encapsulation layers, as shown in Figure 9. In Figure 9a,c, the multilayer structures of the planar and surface-textured CIGS solar cells are the same as those in Figures 2a and 6a except the incoherent front encapsulation layers. In this case, sunlight is assumed to be directly incident from the air to the AZO/ZnS/CIGS/Mo layer, all of which is treated coherently. The device structures of the planar and surface-textured CIGS solar cell with two front encapsulation layers are the same as those in Figures 2a and 6a. The RMS surface roughness of the surface-textured structure is set to be σ RMS = 100 nm, which show the minimum reflectance in Figure 8. Here, we will focus on how the absorption and short-circuit current enhancement of the surface-textured CIGS solar cell over the planar CIGS solar cell can be changed whether the front incoherent encapsulation layers are included or not in the optical modeling. Figure 10 shows the calculated reflectance spectra of the surface-textured CIGS solar cell without and with the inclusion of the front incoherent layers when the incidence angle of the AM 1.5 sunlight is 0 • , 30 • , 45 • , and 60 • . For reference, the reflectance spectra of the planar CIGS solar cell are also calculated. According to the Fresnel equation, the reflectance of the planar interface between two media is proportional to the refractive-index difference of the two media [35]. In reference to the complex refractive index spectra of the materials in Figure 2b, the refractive index of the cover glass, which is very close to that of the EVA, is smaller than that of the AZO, ZnS, or CIGS. Because of the optical interference effect in the AZO/ZnS/CIGS/Mo multilayer, the reflectance of the planar CIGS solar cell without the front encapsulation layers can be lower than that of the CIGS solar cell with the front encapsulation layers in a relatively narrow spectral range around 500 nm. However, overall reflectance of the planar CIGS solar cell with the front encapsulation layers is lower than that of the CIGS solar cell without the front encapsulation layers. This is ascribed to the fact that the refractive-index difference between the air and the cover glass is less than the refractive-index difference between the air and the AZO TCO.
Appl. Sci. 2021, 11, x FOR PEER REVIEW 14 of 18 minimum reflectance in Figure 8. Here, we will focus on how the absorption and shortcircuit current enhancement of the surface-textured CIGS solar cell over the planar CIGS solar cell can be changed whether the front incoherent encapsulation layers are included or not in the optical modeling.  Figure 10 shows the calculated reflectance spectra of the surface-textured CIGS solar cell without and with the inclusion of the front incoherent layers when the incidence angle of the AM 1.5 sunlight is 0°, 30°, 45°, and 60°. For reference, the reflectance spectra of the planar CIGS solar cell are also calculated. According to the Fresnel equation, the reflectance of the planar interface between two media is proportional to the refractive-index difference of the two media [35]. In reference to the complex refractive index spectra of the materials in Figure 2(b), the refractive index of the cover glass, which is very close to that of the EVA, is smaller than that of the AZO, ZnS, or CIGS. Because of the optica interference effect in the AZO/ZnS/CIGS/Mo multilayer, the reflectance of the planar CIGS solar cell without the front encapsulation layers can be lower than that of the CIGS solar cell with the front encapsulation layers in a relatively narrow spectral range around 500 nm. However, overall reflectance of the planar CIGS solar cell with the front encapsulation layers is lower than that of the CIGS solar cell without the front encapsulation layers. This is ascribed to the fact that the refractive-index difference between the air and the cover glass is less than the refractive-index difference between the air and the AZO TCO.
When the AM 1.5 sunlight is normally incident (θ0 = 0°) in Figure 10(a), the reflectance of the surface-textured CIGS solar cell without the front encapsulation layers has a very small value of less than 5% in nearly all the wavelength ranges. This results from the significant reduction of reflectance caused by light scattering effect at the surface-textured interface between the air and the AZO TCO. On the other hand, the reflectance of the surface-textured CIGS solar cell with the front encapsulation layers, which is determined by the incoherent multiple optical reflections as shown in Figure 8, is greater than that of the surface-textured CIGS solar cell without the front encapsulation layers. As the incidence angle of sunlight increases from θ0 = 0°, the degree of increase in reflectance of the surface-textured CIGS solar cell without the front encapsulation layers is relatively higher than that of the surface-textured CIGS solar cell with the front encapsulation layers. Thus when the incidence angle of the AM 1.5 sunlight is θ0 = 60° in Figure 10(d), the reflectance   When the AM 1.5 sunlight is normally incident (θ 0 = 0 • ) in Figure 10a, the reflectance of the surface-textured CIGS solar cell without the front encapsulation layers has a very small value of less than 5% in nearly all the wavelength ranges. This results from the significant reduction of reflectance caused by light scattering effect at the surface-textured interface between the air and the AZO TCO. On the other hand, the reflectance of the surface-textured CIGS solar cell with the front encapsulation layers, which is determined by the incoherent multiple optical reflections as shown in Figure 8, is greater than that of the surface-textured CIGS solar cell without the front encapsulation layers. As the incidence angle of sunlight increases from θ 0 = 0 • , the degree of increase in reflectance of the surfacetextured CIGS solar cell without the front encapsulation layers is relatively higher than that of the surface-textured CIGS solar cell with the front encapsulation layers. Thus, when the incidence angle of the AM 1.5 sunlight is θ 0 = 60 • in Figure 10d, the reflectance of the surface-textured CIGS solar cell without the front encapsulation layers is either lower than or higher than that of the surface-textured CIGS solar cell with the front encapsulation layers in reference to the wavelength of 690 nm. Figure 11 shows the calculated J-V curves of four different device structures when the incidence angles of the AM 1.5 sunlight are 0 • , 30 • , 45 • , and 60 • . The corresponding values of J SC are summarized in Table 1. In the case of surface-textured CIGS solar cells, optical calculation results are averaged over three different arbitrarily-textured surfaces at each RMS roughness to simulate CIGS solar cells with random surface-texture. This averaging procedure for the surface-textured CIGS solar cells can cause the deviation of the calculation results of J SC in Table 1. The calculated standard deviations of J SC for the surface-textured CIGS solar cell without the front encapsulation layer are 0.18, 0.23, 0.22, and 0.34 at the incidence angles of 0 • , 30 • , 45 • , and 60 • , respectively. In the case of the surface-textured CIGS solar cell with the front encapsulation layer, the calculated standard deviations of J SC are 0.04, 0.08, 0.08, and 0.06 at the incidence angles of 0 • , 30 • , 45 • , and 60 • , respectively. Because of the larger refractive-index difference between the air and the TCO, the standard deviation of J SC is speculated to be greater for the surface-textured CIGS solar cell without front encapsulation layer. According to the calculated standard deviations for the surface-textured CIGS solar cells, all the calculation results of J SC are speculated to be within the error range. The surface-textured CIGS solar cells have a larger value of J SC than the planar CIGS solar cells due to the reflectance reduction and absorption enhancement caused by the scattering effect of the surface-textured structure. It is noticeable that the improvement percentage of J SC in the surface-textured CIGS solar cell over the planar CIGS solar cell is different whether the front incoherent encapsulation layers are included or not. According to the calculation results in Table 1, the improvement percentages of J SC in the surface-textured CIGS solar cell without the front encapsulation layers are 9.31%, 9.17%, 8.55%, and 6.38% at the AM 1.5 sunlight incidence angle of 0 • , 30 • , 45 • , and 60 • , respectively. On the other hand, the improvement percentages of J SC in the surface-textured CIGS solar cell with the front encapsulation layers are 5.58%, 5.74%, 5.76%, and 5.51% at the same incidence angle. Thus, the absorption enhancement of the surface-textured CIGS solar cell over the planar CIGS solar cell can be over-predicted when the thick encapsulation layers are not considered in the numerical calculation. Table 1. Calculated short-circuit current density (J SC ) of the planar and surface-textured CIGS solar cells without and with the inclusion of the front encapsulation layers, which are taken from Figure 11. The unit of J SC is mA/cm 2 . The improvement percentage of J SC in the surface-textured CIGS solar cell over the planar CIGS solar cell is denoted in the parenthesis whether the front incoherent encapsulation layers are included or not.

Conclusions
We numerically investigated the effect of the incoherent encapsulation lay oblique sunlight incidence on the reflectance spectra and J-V curves of surface-te CIGS solar cells. The optical modeling of the two incoherent encapsulation laye performed based on the proposed ADETAM, which provided the incident-angle-d ent reflectance spectra of the mixed coherent-incoherent planar and surface-texture solar cells. The calculation accuracy of the ADETAM was verified in a planar CIG cell through a very good agreement with exact analytical results obtained by the G In addition, The J-V curve of the CIGS solar cell was calculated based on the 2D

Conclusions
We numerically investigated the effect of the incoherent encapsulation layer and oblique sunlight incidence on the reflectance spectra and J-V curves of surface-textured CIGS solar cells. The optical modeling of the two incoherent encapsulation layers was performed based on the proposed ADETAM, which provided the incident-angle-dependent reflectance spectra of the mixed coherent-incoherent planar and surface-textured CIGS solar cells. The calculation accuracy of the ADETAM was verified in a planar CIGS solar cell through a very good agreement with exact analytical results obtained by the GTMM. In addition, The J-V curve of the CIGS solar cell was calculated based on the 2D FEM-based numerical solution of the drift-diffusion model.
To investigate the effect of the front incoherent encapsulation layers on output performance of the CIGS solar cell, the angular dependences of reflectance spectrum and J-V curves were calculated and compared in both planar and surface-textured CIGS solar cells with and without the inclusion of the two front encapsulation layers. The calculated reflectance of the planar CIGS solar cell with the front encapsulation layers was greater than that of the planar CIGS solar cell without the front encapsulation layers because the refractive-index difference between the air and the cover glass was less than the refractiveindex difference between the air and the AZO TCO. In contrast, the calculated reflectance of the surface-textured CIGS solar cell with the front encapsulation layers was lower than that of the surface-textured CIGS solar cell without the front encapsulation layers. This was ascribed to the fact that the reflectance of the surface-textured CIGS solar cell with the front encapsulation layers, which was determined by the incoherent multiple optical reflections, was greater than that of the surface-textured CIGS solar cell without the front encapsulation layers, where the significant reflectance reduction caused by light scattering effect occurred at the surface-textured interface between the air and the AZO TCO. Correspondingly, the improvement percentages of J SC in the surface-textured CIGS solar cell over the planar CIGS solar cell with the front encapsulation layers are lower than that in the surface-textured CIGS solar cell over the planar CIGS solar cell without the front encapsulation layers. Thus, the output performance of the surface-textured CIGS solar cell over the planar CIGS solar cell can be over-predicted when the thick encapsulation layers are not considered in the optoelectronic modeling. We expect that the proposed ADETAM, together with comprehensive numerical modeling and analysis, will give a deeper understanding on the simultaneous effect of the thick front encapsulation layers and oblique sunlight incidence on output performance of the CIGS solar cells and help to provide a more efficient optical design of thin-film CIGS solar cells in the consideration of realistic deployment situations.