Multi-dimensional optimization of In0.53Ga0.47As thermophotovoltaic cell using real coded genetic algorithm

The optimization of thermophotovoltaic (TPV) cell efficiency is essential since it leads to a significant increase in the output power. Typically, the optimization of In0.53Ga0.47As TPV cell has been limited to single variable such as the emitter thickness, while the effects of the variation in other design variables are assumed to be negligible. The reported efficiencies of In0.53Ga0.47As TPV cell mostly remain < 15%. Therefore, this work develops a multi-variable or multi-dimensional optimization of In0.53Ga0.47As TPV cell using the real coded genetic algorithm (RCGA) at various radiation temperatures. RCGA was developed using Visual Basic and it was hybridized with Silvaco TCAD for the electrical characteristics simulation. Under radiation temperatures from 800 to 2000 K, the optimized In0.53Ga0.47As TPV cell efficiency increases by an average percentage of 11.86% (from 8.5 to 20.35%) as compared to the non-optimized structure. It was found that the incorporation of a thicker base layer with the back-barrier layers enhances the separation of charge carriers and increases the collection of photo-generated carriers near the band-edge, producing an optimum output power of 0.55 W/cm2 (cell efficiency of 22.06%, without antireflection coating) at 1400 K radiation spectrum. The results of this work demonstrate the great potential to generate electricity sustainably from industrial waste heat and the multi-dimensional optimization methodology can be adopted to optimize semiconductor devices, such as solar cell, TPV cell and photodetectors.

www.nature.com/scientificreports/ substrate, corresponds to E g and λ c of 0.74 eV and 1.68 µm, respectively. Moreover, In 0.53 Ga 0. 47 As is a promising TPV cell due to its high crystal quality and the cost-effectiveness of InP substrate, making it suitable for large scale production as compared to other TPV materials. It is worth mentioning that the existing epitaxy growth technology of metal-organic vapor-phase epitaxy (MOVPE) has the ability to produce In 0.53 Ga 0. 47 As/InP heterojunction with high crystal quality and low defect density 26,27 . The main structure of In 0.53 Ga 0. 47 As configuration includes the emitter, base, front surface field (FSF), back surface field (BSF), cap and buffer layers. In previous literature, the base thickness was reported between 1 and 5 µm, and emitter thickness was between 0.05 and 0.44 µm [27][28][29] . Several structures reported the use of highly doped In 0.53 Ga 0.47 As cap layer ~ 1 × 10 19 cm −3 and highly doped InP BSF/buffer layer ≥ 1 × 10 18 cm −3 . Table 1 reviews the design structure and output performances of In 0.53 Ga 0. 47 As photovoltaic (PV) cell under different testing conditions. Typical In 0.53 Ga 0. 47 As cells have open-circuit voltage (V oc ), short circuit current density (J sc ), fill factor (FF) and efficiency (η) ranging between 0.26 and 0.45 V, 18.8 and 64.5 mA/cm 2 , 59 and 74.2%, and 4.2 and 14.37%, respectively under air-mass 0 (AM0) and air-mass 1.5 (AM1.5) illuminations 26,28,30,31 . In 2019, Omair et al. 32 reported a TPV cell efficiency of 29.1% under 1480 K radiation temperature, by recycling sub bandgap photons to the radiator. It is possible to achieve > 50% cell efficiency by improving the series resistance, material quality and reflectivity using chamber with high mirror reflectivity. Recently, Fan et al. 33 presented an air-bridge structure that enables photons recycling, increasing the reflection of sub-bandgap photons to 99% and enhancing the efficiency up to 30%. As summarized in Table 1, the TPV testing conditions resulted in higher output performance as compared to solar spectrums. Moreover, commercially available TPV cell has the advantage of producing high output power density, ~ 26 times higher than the output power of solar PV cell 34,35 . All the reported work in the optimization of In 0.53 Ga 0. 47 As cell are based on the alteration of single design variable 36,37 . Recently, a multi-variable optimization of solar cell was used to optimize the physical properties of electron transport materials, hole transport materials, and metal contact and layers thickness. It was accomplished using the MATLAB optimization toolboxes incorporated with one-dimensional (1D) Solar Cell Capacitance Simulator (SCAPS) for device simulation 40 . However, the study did not take into account the impact of doping concentration on the cell performance.
The optimization of TPV cell structure is critical in getting the highest achievable η. A slight increment in η will significantly increase the output power and total energy. The simplest method to optimize a structure is the single-layer/variable optimization, which optimizes only a single parameter at a time while other parameters are kept constant. Several attempts were made to optimize single variable, especially on the emitter and base layers of In 0.53 Ga 0. 47 As cell 36,37 . However, device performance depends collectively on all the design variables 45,46 , and a more heuristic optimization that considers the effect of all important variables for the In 0.53 Ga 0. 47 As TPV cell is necessary to achieve the optimum cell efficiency. Therefore, this study investigates the effect of each variable through single variable optimization and performs multi-dimensional (simultaneous multi-variables) optimization using real coded genetic algorithm (RCGA) to obtain the optimum configuration of In 0.53 Ga 0.47 As TPV cell.  38 n/a n/a n/a n/a n/a n/a ~ 0.225 0. 5 27 . The input parameters are structure design and testing conditions. The structure design includes front and back gold contacts, thickness and doping concentration of the emitter, base, BSF, FSF, buffer and cap layers, as shown in Fig. 1a. No anti-reflective coating (ARC) was considered in the structure. The validation testing is performed under a standardized solar spectrum AM1.5 at room temperature (300 K). The output performance parameters that are of primary interest include JV-curve, J sc , V oc , FF and η. The simulation took into consideration the Auger, radiative and Shockley-Read-Hall (SRH) recombination as well as the carrier's lifetime and mobility concentration-dependent models. At a doping concentration of 1 × 10 17 cm −3 , lifetime of 16 ns for electrons and 40 ns for holes were used based on the model. While the In 0.53 Ga 0.47 As model and testing conditions remain constant over the validation process, the materials parameters of the In 0.53 Ga 0.47 As and InP were varied within the reported range in previous literature. Table 2 summarizes the electrical properties of In 1-x Ga x As as a function of x composition and concentration-dependent models. www.nature.com/scientificreports/ Using the identified material parameters in Table 2, the current density-voltage (JV) characteristics of the In 0.53 Ga 0.47 As cell are obtained. The generated JV characteristics of the simulation model can be seen in Fig. 1b. It can be observed that a close agreement was obtained between the performance parameters of the simulation model and the reported experimental data. A percentage error of less than 1% was achieved for each parameter. For example, a percentage error of 0.61% between the experimental and simulation results was calculated for η. In addition, it was reported by other literature that the performance parameters of η, FF, V oc and J sc , under similar testing condition (AM1.5), were in the range of 9.3-12.9%, 68-71%, 0.31-0.39 V and 21.5-42.8 mA/cm 2 , respectively 30,41,42 . The simulation results for the validation work achieved in this study are within the reported range; hence, validate the In 0.53 Ga 0.47 As cell model in this work.
Single layer/variable optimization. A single variable optimization gives an indication of the significance of a variable to the TPV cell performance, reveals the trend of performance variation for each variable and more importantly it identifies the range of the design parameters accurately, and that provides the RCGA with a much faster convergence speed as well as a higher solution accuracy. Based on the reported values for the thickness and doping concentration of each layer in Table 1, their maximum and minimum values were first estimated. However, since the majority of the In 0.53 Ga 0.47 As structures are used for PV application, several simulations were conducted to modify the upper and lower boundary conditions to suit TPV testing conditions. Table 3 shows the range of the design parameters for the In 0.53 Ga 0.47 As TPV cell. Individually, each variable was manipulated Table 2. Material parameters of In 1-x Ga x As and In 0.53 Ga 0.47 As at 300 K. * z represents mobility fitting parameter, *γ represents lifetime fitting parameter, * R SRH is SRH recombination, * C rad is electron and hole radiative coefficient, * C n,p are electron and hole radiative coefficient, * C is auger coefficient, * n is hole densities, * p represents hole densities, N D,A doner and accepter doping concentration and * n i represents intrinsic region. Hole density of states 26,47     www.nature.com/scientificreports/ while the rest of the design parameters in Fig. 1a remained constant. In the BSF analysis, base thickness was fixed at 10 µm to reduce the absorption in BSF, thereby allowing the investigation of the layer functionality as field generator. The optimization process was conducted under radiation temperatures from 800 to 2000 K at 50% beam illumination intensity. The 50% beam intensity was preferred as selective radiator usually emitted ~ 50% less power than an ideal blackbody 59 . In addition, a low pass optical filter at 2 µm was employed in the simulation following the reported TPV system 60-62 . Fourspring et al. 63 illustrated the use of the edge short-pass filter and emphasized that the optical interference material has high spectral efficiency which reduces the amount of energy that reaches the TPV cell for wavelengths higher than 2 µm. The cell temperature was maintained at 300 K, assuming that it is controlled using an effective thermal management system 64 . The consideration of heat transfer between the TPV cell and the environment can be a more crucial subject of investigation for TPV devices with nano-gap between the device and the heat emitter 65 .

Material properties at 300 K (x) is the composition and (N) is the doping level Range Numerical values determined in this work
Multi-dimensional optimization using real coded genetic algorithm method. The multi-dimensional optimization performs complete iterations for all possible combination of different variables to obtain the optimum values for all variables that achieve the highest efficiency. A flow chart of the multi-variable optimization of In 0.53 Ga 0. 47 As structure under different radiation temperatures is presented in Fig. 2. The optimization process consists of both the device simulation module and the numerical optimization. Firstly, the device simulation module, where ATLAS was used to simulate the 2D In 0.53 Ga 0.47 As model. An initial population size of 50 input vector X was implemented in the final optimization version. It was demonstrated that a higher population size improves the accuracy and reduces the number of generation (Gen) required to allow the iterations to converge to the optimum In 0.53 Ga 0.47 As configuration. X is given in Eq. (1) as the six-layer/twelve design variables to be optimized. www.nature.com/scientificreports/ where i = 1-50 is the initial population size, n = 1-12 is the design variables based on some lower and upper physical constraints Table 4. DBinternal tool is used to interfaced 50 sets of X to Deckbuild ATLAS, then solves the optically excited charge generation, 2D Poisson's Equation, transport Equation, and continuity Equations and calculates the efficiency (η(X)) of the cell based on the inputs. Moreover, the numerical optimization of RCGA method was performed. Real coded is a direct representation of the variables, where no coding and encoding is required 66 . The objective function is to maximize the η of the In 0.53 Ga 0.47 As TPV cell, as presented in Eq. (2).
where η is 20,67 : Q emit is the thermal emission of the blackbody at a fixed temperature (T BB ), and the Q filter is the thermal reflected emission of the low pass optical filter. The photon flux of an emitting blackbody, Φ, as a function of the emitted λ and T BB , is calculated via Planck's Law in Eq. (4).
where h is the Planck's constant, k B is the Boltzmann's constant and c o is the speed of light.
The P out of the cell is defined as: where the photocurrent, J m , as a function of voltage across the cell, V m , is the difference between the short circuit current and recombination loss, given by: where R rad , R SRH and R Aug , are the radiative, SRH and Auger recombination rates, respectively. Based on the principle of detailed balance, the J sc is calculated from the external quantum efficiency EQE(λ) and the Φ(λ): The fitness ratio is defined as the efficiency ratio (η r (X)), and it is shown in Eq. (8). The chromosomes are arranged based on their fitness from higher to lower (i` = 1-50), and then some evolution mechanisms like survivor selection, crossover and mutation were used to build the next generation (Gen) using Excel and Visual basic. A 20% (10-best fitted chromosomes: X i = X 1 to X 10 ) of the best chromosomes survived and directly passed to the next Gen as X i ′ while 80% (40-best fitted chromosomes) of the best-selected chromosomes go to the next step of crossover and mutation producing a new set of population 68 . This process is presented in Eq. (9).
(1) www.nature.com/scientificreports/ After some Gen's simulations, child chromosomes were created from the best performing parent chromosomes producing the optimum In 0.53 Ga 0.47 As configuration. A stopping criterion was decided after observing no significant change in the efficiency and efficiency ratio within the Gen (η r (X) = 2 ± 0.3). In this way, the optimization process was repeated for illumination source temperatures from 800 to 2000 K, with an interval of 200 K.

Results and discussion
Single layer/variable optimization. The effect of varying the thickness and doping concentration of cap, FSF, emitter, base, BSF and buffer layers on the performance parameters (J sc , V oc , FF, η) are tabulated in Table 5. The variables are classified into three categories: insignificant where the variation in η is ≤ 0.4%, significant where the change in η is between 0.4 and 3% and highly significant where the change in η is ≥ 3%.
As can be seen from Table 5, the V oc is not significantly affected by the doping concentration of the cap, FSF, BSF and buffer layers. Theoretically, V oc is influenced by dark current densities (J 01 and J 02 ), where J 01 is contributed to the dark current due to surface and bulk recombination losses, and J 02 is related to recombination due to traps in the space charge region (SCR) 69 . However, the implementation of InP in the front and rear side of the junction reduces the front surface recombination and back surface recombination, eliminating the dark current across the surface. The variation of cap, FSF, BSF and buffer doping concentration produces no effect on the V oc of the cell since no absorption and recombination have occurred. This is because the Type-II band alignment (staggered band) between InP/In 0.53 Ga 0.47 As and In 0.53 Ga 0.47 As /InP led to spatial separation of electrons and holes and passivated the surfaces 70 . However, the thickness increment of those layers will affect the V oc because of the losses due to the light aborption 71 .
The J sc is mainly related to the absorption and diffusion length of the photo-generated carriers, as shown in Eq. (10). 72 .
where G is the generation rate, q is the charge, and L n and L p are the electron and hole diffusion length, respectively.
Since In 0.53 Ga 0.47 As material has a long electron and hole diffusion lengths, the increment of emitter, base and BSF thicknesses will significantly improve the absorption, carrier's generation and J sc . For instance, the generated electrons (holes) at absorber layers has a diffusion length of 18.89 (5.88) µm at 1 × 10 17 cm −3 doping concentration. The increment of base layer thickness from 1 to 18.89 µm will increase the J sc since majority of the generated minority carriers (electrons) are able to reach the SCR before they are recombined. The doping concentrations of the absorber (emitter, base and BSF layers) affect the probability of carrier recombination as the mobility, lifetime and diffusion length of carriers decrease with higher doping concentration. However, since the emitter layer and BSF layer were remained constant at 0.05 µm and 0.025 µm, the variation of doping concentration gives a minor effect on the cell J sc .
The FF is defined as the ratio between the V oc I sc to the actual operating condition of the cell P mp = I mp V mp after considering the series resistance and shunt resistance in the structure 73 . The P mp , I mp and V mp are denoted as the cell maximum power, maximum current and maximum voltage, respectively. It can be seen that the FF and η performance parameters were affected by almost all of the design variables. This is due to the domination of resistance losses when manipulating the variable individually. Nevertheless, the manipulation of FSF and BSF (9) (10) I sc = qG(L n + L p ) Table 5. The summary of the result for single variable optimization for blackbody temperature from 800 to 2000 K. *√ Changing the variable affects the output parameter (percentage improvement by ± ≥ 1%) and * × Changing the variable does not affect the output parameter (percentage improvement by ± < 1%).

Designing variables
Optimum range www.nature.com/scientificreports/ doping concentration had no significant impact on cell performance since the thickness of those layers are kept at very thin (~ 0.02 µm), resulting in lower resistance losses and minor absorption. Based on the single variable optimization, it was found that the manipulation of the thickness and doping concentration variables for the base layer significantly affects all performance parameters. In particular, the highest simulated η is achieved at an optimum base layer thicknesses of 5 and 16 µm with a blackbody temperature of 2000 K and 800 K, respectively. The variation of the optimal base thicknesses is related to the radiation spectrum and its peak wavelength (λ p ). Based on Wien's displacement law, λ p of the spectrum shifts toward shorter IRs for higher blackbody temperature hence thinner base layer is required. Low doping in the base layer led to a high minority carrier lifetime and long diffusion length, which could increase the probability of carriers reaching the contact before recombining. However, it could also reduce the conductivity of the layer, electric field, and the built-in potential at the junction that decreases the V oc performance, offsetting the improvement in J sc . On the other hand, a thicker base layer will increase the absorption of IRs, which results in higher J sc . However, it could also increase the recombination and shunt resistance 74,75 . Increasing the thickness is detrimental to cell performance since it causes higher recombination, especially if the structure has high SRH rate 76 . Moreover, the minority carriers generated in the lower region of the cell with diffusion length shorter than the thickness has a higher probability of recombining before reaching the SCR 74 . Therefore, it is worth exploring the tradeoff relationship of the base layer thicknesses and doping concentrations to acquire the optimum configuration of the base layer that produces the highest cell efficiency. The correlation of cell efficiency with the base layer thickness and base doping concentration is depicted in Fig. 3. The base thickness (doping concentration) was varied between 1 (6 × 10 16 ) and 28 µm (1 × 10 19 cm −3 ), while the rest of structure design variables were kept at their baseline values.
Based on Fig. 3, it was found that an optimum η of 17.61% is obtained at an optimum base thickness (doping concentration) of 11 µm (1 × 10 17 cm −3 ). The increment in base thickness improves the absorption of IR, especially the band-edge photons 76 . Based on the Einstein relationship, Caughey-Thomas model and SRH model, the diffusion length of electron and hole are calculated as a function of doping concentration. At base doping of 1 × 10 17 cm −3 , the diffusion length of electron (hole) is equal to 18.89 (5.88) µm. The electron diffusion length of 18.89 µm is longer than the 11 µm base thickness. Since the diffusion length is longer than the base thickness, the probability that photo-generated carriers can reach the SCR is quite high 74 . Furthermore, the high diffusion coefficient of In 0.53 Ga 0.47 As carriers permits the formation of the active junction at low doping concentration. In 0.53 Ga 0.47 As cell ability to produce active junction at low doping concentration allows it to obtained high cell efficiency with minimum recombination rate.
Multi-dimensional optimization using real coded genetic algorithm method. The results of multi-dimensional optimization on In 0.53 Ga 0.47 As TPV cell model under 800-2000 K blackbody temperatures are summarized in Table 6. It was found that the optimized cell presented in this work has thicker base layer as compared to those TPV cells reported in Table 1. At 1800 K blackbody temperature, the optimized cell has an efficiency of 20.48%, compared to 15% reported for TPV cell tested under similar testing conditions 40 . Compared to the reported In 0.53 Ga 0.47 As TPV cell efficiencies which can be as high as 30% 32,33,77 , the optimized structure in this paper does not consider photons recycling and ARC. Instead of utilizing photon recycling, this work focuses on the TPV cell design structure. A thicker base layer between 16 and 18 µm is preferred to enhance the current density as it increases the absorption of near band-edge photons and free-carrier absorption (FCA), which highly impact TPV cell efficiency 78,79 . At base layer thickness of 16-18 µm, the maximum cell efficiency of 23.18% was reported at 1400 K radiation temperature, while the minimum efficiency of 18.41% was achieved at www.nature.com/scientificreports/ 800 K. The electrical losses due to the high photogenerated carriers increased at radiation temperatures > 1400 K, which reduced the η. Although maximum efficiencies were achieved with base layer thickness of 16-18 µm, considering the technical difficulties and cost in growing high-quality In 0.53 Ga 0.47 As layer at > 10 µm, the base layer thickness was reduced to 8 µm. With this reduction in base layer thickness, the efficiency of the TPV cell is reduced only by an average of 1% (Maximum reduction of 1.9% at 1000 K and minimum reduction of 0.1% at 2000 K). RCGA practical represents the RCGA results after taking into account the practicality in the growth of semiconductor wafer by reducing the base layer thickness to 8 µm. A thinner cap and buffer layers are designed with a higher doping concentration to form a better front and rear ohmic contacts, resulting in higher FF and lower series resistance. Besides, a thin FSF and BSF layers are designed with a higher doping concentration, resulting in a better front and rear junction passivation with minimum optical absorption. A proper cell design has to consider both optical and electronic losses. The optical losses were the dominant factor in determining the change in efficiencies of the non-optimized cell. Optimized cell efficiency, on the other hand, were dominated by the electrical losses, especially at high radiation temperature with high energy spectrum density. The In 0.53 Ga 0.47 As TPV cell has to be optically thick (i.e. to absorb all or most of the incident illumination) and electronically thin (i.e. to collect the photoexcited electron-hole pairs with little or no losses). These two requirements lead to an optimal configuration that maximizes efficiency. At 1400 K blackbody temperature, J sc and FF are increased respectively from 755.01 to 1719.35 mA/cm 2 and 62.65 to 68.63% after optimizing the entire cell configuration. Most of the optimization works focused on optimizing the electrical losses of cell to improve the V oc performance by reducing the thickness of the absorber layer 27,80 . However, the absorption of near band-edge photons is neglected. Alharbi et al. 81 highlighted that the leading cause of the reduced efficiency of a solar cell below the theoretical limit is the drop in the estimated V oc while usually, the obtained J sc is around the theoretically maximum values. It is important to mention that the solar spectrum is mainly concentrated around the visible region, and these photons do not require thicker base absorber. On the other hand, TPV illumination flux is usually concentrated at infrared wavelengths, and thicker absorber is needed to improve the absorption of IRs and significantly increase J sc . This statement was supported by the optimization study reported by Baudrit and Algora 46 , where increasing the absorber thickness of the bottom cell in GaInP/GaAs dual-junction increases the value of photocurrent density. Furthermore, the J sc and η were increased respectively from 13.85 to 15.62 A/ cm 2 and 32.6 to 36.4% under 1000 suns concentration.
For the validation of optimization, the In 0.53 Ga 0.47 As TPV structure reported from the RCGA optimization method was validated by two validation methods (VM), which are single-variable method (VM 1) and iteration method (VM 2). VM 1 utilized the results from the single-variable optimization where the optimum values of the 12 design parameters obtained through single-variable optimization were used to determine the efficiency of the TPV cell. VM 2 utilized the iteration method to simulate all possibilities for the 12 design parameters. However, it was estimated that the total number of iterations are 12 10 simulation runs. To reduce the computation time to a manageable level, we used a bigger step and only varied the important design variables while the insignificant variables were maintained at their baseline values. Due to these reasons, it is expected that the cell efficiencies obtained through VM 1 and VM 2 are slightly lower than those obtained through RCGA method.
The results obtained through RCGA, VM 1 and VM 2 are tabulated in Table 7. For different radiation temperatures, the average percentage differences between RCGA and VM 1, and between RCGA and VM 2, were 5.17% and 6.95%, respectively. On the other hand, RCGA method is able to converge to the optimum structure in a faster manner as compared to the long simulation time needed by the iteration method. Figure 4 illustrates the EQE of non-optimized and RCGA practical -optimized In 0.53 Ga 0.47 As TPV cells as a function photon λ. It can be seen from Fig. 4 that the EQE(λ) of the multi-variable optimized cell is higher than that of the non-optimized TPV cell. This is due to the enhancement of photocurrent generation at a wavelength higher than 0.85 µm. A possible explanation for this is that the optimized In 0.53 Ga 0.47 As TPV configuration tends to have higher absorption and collection of IRs photo-generated carriers. Despite the reduction in V oc by 6.12% after the optimization process, J sc is significantly increased by 127.73%, resulted in the improvement in the TPV cell efficiency by 138.23%.
Based on the optimized results, a significant increase in η is attained when the base layer increased from 1 to 8 µm. This finding can be supported by an analysis of the absorption coefficient and the absorption length of where k(λ) is the extinction coefficient of In 0.53 Ga 0.47 As 82,83 . On the other hand, the absorption length (α −1 ) is given as the inverse of the α and describes the penetration length of majority photons in semiconductor before being absorbed. Due to absorption in the material, the illumination intensity weakened with increasing penetration length and can be described by means of a decaying exponential function, as shown in Eq. (12). 74 .
where η Abs is absorption efficiency, the effective cell thickness (d) can be determined from Eq. (13) by maximizing the absorption in the cell η Abs = ~ 99% and R = 0 74 . Figure 5 shows the result of the α, α −1 and d of In 0.53 Ga 0.47 As as a function of λ. It can be seen that thicker absorber layer is needed to absorb IRs. For instance, photons with λ of 1.65 μm have an absorption length of 4.25 μm. To effectively absorb ~ 99% of the photons, the cell thickness will therefore be approximately 19.59 μm.
The In 0.53 Ga 0.47 As has a spectral response with wavelength at around 1.75 μm. It can be seen that infrared light requires a thick absorber, so that majority of the photons are absorbed. The illumination intensity of the blackbody is mainly concentrated at λ > 1.0 μm. For instance, 1400 K blackbody temperature has 69.25% of power density for wavelengths between 1 and 1.8 μm as compared to 3.43% of power density for wavelengths (11)  www.nature.com/scientificreports/ between 0.2 and 1 µm. An optical BSR is used to improve the light absorption, as the optical path distance can be doubled due to the back reflection 20 . Therefore, although the effective thickness to absorb photons up to 1.75 μm is approximately 32-36 μm (based on Fig. 5), the effective device thickness can be reduced to approximately 16-18 μm. However, based on additional simulation work upon completing the multi-dimensional optimization, an effective thickness higher than 8 μm (up to 18 μm) will only increase the efficiency by an average of about 1% for various blackbody temperatures due to the sharp decline of the absorption coefficient for λ ≥ 1.63 µm. Other than that, the light trapping method such as Lambertian rear reflector and textured surface can be used to improve the light absorption in the cell 74,84 . A further explanation for the high EQE is due to the effective separation and collection of generated carriers. In 0.53 Ga 0.47 As structure is frequently constructed with FSF and BSF layers to reduce the surface recombination and enhance the V oc 41,85 . High to low doping concentration between the FSF or BSF and the active junction is vital to generate a SCR similar to the SCR between n-p junction. For instance, holes diffused out of this highly doped BSF layer into the lower-doped base layer, leaving site-fixed negatively charged acceptor atoms behind 74 . The generated electrical field that acts like an electric mirror that returns the electrons generated through absorption in the direction of the SCR 74 . The probability of undesired recombination at the rear of the cell is thus significantly reduced. Belghachi et al. 86 described a mathematical model on the importance of high-low junction in the front and rear sides of GaAs cell, which play a crucial role in enhancing the light-generated free carriers' collection. The thickness of FSF and BSF layers should be as thin as possible while the doping concentration should be as high as 1 × 10 19 cm −386 . The band diagram is an alternative way to view the effect of FSF and BSF layers. The band diagram of the In 0.53 Ga 0.47 As/InP cell was extracted from TCAD model, with the FSF and BSF layers. Figure 6 presents the band diagram of In 0.53 Ga 0.47 As heterojunction. The band offset between In 0.53 Ga 0.47 As and InP forms type-II band (staggered) leads to spatial separation of electrons and holes. A discontinuity in conduction for base/(BSF/buffer) interface prevents the electrons from moving further up to the back contact while allowing the flow of holes. The discontinuity in valance for the FSF/emitter interface prevents the holes from moving further up to the front contact while allowing the flow of electrons 74 . The band alignment of In 0.53 Ga 0.47 As structure was investigated in several studies 87,88 . The n-InP will form a barrier to holes, but a sink for electrons at the FSF and visa verse at the BSF/buffer. Besides that, cell total voltage is now divided into the potential level at n-p junction and additional level at n + n and pp + interfaces 74 . Hence, FSF layer improves the EQE of shorter wavelengths photo-generated carriers, and the BSF/buffer layer increases the EQE of longer  www.nature.com/scientificreports/ wavelengths photo-generated carriers 86 . It should be stressed that BSF/buffer layer is very significant to enhance the J sc and V oc of the TPV cell since the long-wavelength photons tend to absorb at the deeper region of the structure. The combination of both thicker absorber layer and BSF/buffer significantly improves the collection of photocurrent collection of λ near to the band-edge of In 0.53 Ga 0.47 As TPV cell.

Conclusion
In summary, the simulation model is validated with reported experimental data, generating a low η percentage error of 0.61% between experimental and simulation. Research gap was identified based on comprehensive comparison of previous structure designs and performance of In 0.53 Ga 0.47 As cell. Single variable and multidimensional optimization (RCGA) of heterojunction In 0.53 Ga 0.47 As cell are developed and applied at 800-2000 K TPV radiation temperatures. The single variable optimization is used to investigate the effect of thickness and doping concentration of layers, which demonstrated the significant impact of the base layer to achieve high performance In 0.53 Ga 0.47 As TPV cell. Meanwhile, under radiation temperatures ranging from 800 to 2000 K, optimized In 0.53 Ga 0.47 As heterojunction TPV cell using RCGA increases the η by an average of 11.86% as compared to the reference structure. It was found that the incorporation of a thicker absorber with effective barrier layer BSF/Buffer layer improves the absorption and collection of photo-generated carriers near the band-edge, which produced higher output performance. The increment of the In 0.53 Ga 0.47 As cell efficiency leads to a significant increase in the generated output power, demonstrating great potential of TPV for industrial waste heat harvesting. Finally, the method of hybridizing the Silvaco TCAD software with RCGA for multi-dimensional optimization can be readily adopted in optimization work for semiconductor devices, such as solar cell, TPV cell and photodetectors. www.nature.com/scientificreports/