Crystalline silicon (c-Si) solar cells have achieved enormous dominance in the photovoltaic market due to their high efficiencies, excellent material quality, and continual reduction in manufacturing costs. The current record efficiency of 26.7%1 is already 89% of the theoretical limit value for the power conversion efficiency for single-junction solar cells.2 Further improvement in this technology will become increasingly difficult. However, to further boost the distribution of photovoltaics in the future, higher efficiencies at low costs are required.3
We aim to surpass the Shockley–Queisser limit4—the most important technological limit for single-junction solar cells—with tandem solar cells, which combine two solar cells with different bandgaps. The incident sunlight first hits the top cell, which has a higher bandgap and harvests the high-energy photons at a higher voltage, whereas the low-energy photons are transferred to the bottom cell, which has a lower bandgap and corresponding lower voltage. In this way, high-energy photons are able to contribute more voltage to the device instead of losing their excess energy by thermalization.
Lead halide perovskite materials have excellent optical properties for tandem applications due to their steep absorption edge5 and a tunable bandgap.6,7 As a result, tandem cells with a perovskite top cell and a silicon bottom cell have the potential to reach efficiencies beyond 30%,89.–10 if optical losses are adequately adressed.11
In monolithic (two-terminal) tandem cells, the perovskite top cell and the c-Si bottom cell are electrically connected in series. Hence, for high efficiencies the photocurrent density of the top and bottom cell must be matched. Recently, we numerically maximized the achievable photocurrent density for different planar monolithic perovskite-silicon tandem device architectures with antireflective coatings by optimizing the layer thicknesses. The architecture with the electron-selective contact and -doped layers on the front of the perovskite and silicon cells, respectively,12 enables efficiencies exceeding 30%, when the perovskite bandgap is optimized as well.13 However, for fully planar devices, reflection losses cause a significant limitation for the matched photocurrent density .
Textured interfaces can reduce reflection losses via (1) enhancing coupling of light into the structure and (2) scattering light such that the average path length in the absorber is increased, leading to increased absorption especially in weakly absorbing regimes. In recent years, much effort has been put into developing two-dimensional (2-D) structures for light trapping.14220.127.116.11.19.20.–21 For nanostructures, which texture the electrically active layers of the solar cell, it is mandatory that they have no detrimental effect on the electric solar cell performance. Therefore, we decided to focus on hexagonal sinusoidal nanotextures,22,23 which allow for the combination of a strong antireflective effect with good electrical performance, as we demonstrated for liquid-phase crystallized silicon thin-film solar cells.24,25
In this work, we investigate how sinusoidal nanotextures affect the optical performance of perovskite–silicon tandem solar cells. To efficiently optimize the perovskite layer thickness for current matching, we use Newton’s method. As motivation, we demonstrate the experimental feasibility of spin-coating perovskite layers onto sinusoidally nanotextured substrates.
So far, little experimental evidence of nanostructured perovskite layers has been published: in some studies, different nanostructures were imprinted into the perovskite top surface.2627.–28 Although these top surface patterning techniques demand postprocessing of the perovskite film after fabrication, other studies report perovskite grown onto different kinds of nanophotonically patterned charge selective layers.29,30
Here, we demonstrate how sinusoidally nanotextured films can be fabricated with spin coating. This deposition method is chosen since it has shown to yield the best performing perovskite absorber layers for photovoltaic application so far. As model systems, we use two kinds of glass substrates, which are covered with sinusoidal nanotextures. The periods are and ; the respective valley-to-peak heights are and , which correspond to aspect ratios of and , respectively. The nanotextures were manufactured with nanoimprint lithography using a hybrid polymer with glass-like properties (OrmoComp by Micro Resist Technology GmbH). Details on our nanoimprint lithography process can be found in Ref. 31.
For better imaging capabilities, we thermally evaporate 30 nm of silver onto the nanotextured substrates. In addition, we prepare a flat reference substrate by spin-coating a thin layer of polytriarylamine (PTAA) onto an indium tin oxide (ITO) substrate. PTAA is a common hole transporting material (HTM) in p-i-n perovskite solar cells,3233.–34 which supports the uniform growth of perovskite.
For the perovskite layer fabrication in our experiment, we use the mixed halide perovskite composition using a recipe by Saliba et al.35
We spin-coated the as-prepared 0.8-M solution according to the antisolvent method to promote the crystallization into perovskite. This method enables crystallization starting from the top interface of the perovskite36 film, which allows for the growth of films with a high optoelectronic quality, which is mostly independent of the (nanotextured) substrate.37
We investigate the morphology of the perovskite layers with scanning electron microscopy (SEM). Figures 1(a) and 1(b) show cross-sectional SEM images of the perovskite on the two sinusoidal textures with and . The perovskite films show thicknesses of 420 to 500 nm and 420 to 600 nm on the and nanotextured substrates, respectively. The layers are thicker than those required for good current matching, as we will discuss in Sec. 4 and Table 2. However, they can be easily adjusted by reducing the molarity of the perovskite precursor solution.
The perovskite fills the sinusoidal structures for both periodicities and can compensate a height difference of up to 200 nm on the substrate. On the top surface, the underlying sinusoidal nanotexture is not visible, which is in accordance with crystallization starting from the top surface.
We further investigate the perovskite surface with atomic force microscopy (AFM) measurements. Figure 2(a) shows a top view AFM micrograph () of the bare nanotextured substrate with . The hexagonal arrangement of the sinusoidal structure is clearly visible. In contrast, the micrographs in Figs. 2(b) and 2(c) show the disordered surfaces of perovskite layers on this structured substrate and on (flat) PTAA, respectively. The perovskite films grown on these different substrates have a root mean square (RMS) roughness of 16 and 22 nm for the structured and flat substrates, respectively. The RMS roughness values of the perovskite layers in this study are larger than reported elsewhere,38,39 which can be attributed to differences in compositions, processing, and substrates.
Figures 2(d)–2(f) show the Fourier transforms (FTs) corresponding to the AFM height profiles from Figs. 2(a)–2(c). The FT of the bare substrate features six hexagonally arranged spots of high intensity around the center and low intensity spots of higher order at larger wavenumbers. In contrast, the FTs of both perovskite layers in Figs. 2(e) and 2(f) feature a multitude of different wavevectors distributed over all angles. The highest amplitudes are at low frequencies at low wavenumbers, reflecting some unordered base signal, which is modified by low amplitude, high-frequency components to describe the small grains, which were also found in Figs. 2(b) and 2(c). The FT of the perovskite layer on structured substrate [Fig. 2(e)] resemblances that of the flat perovskite layer on PTAA [Fig. 2(f)] much better that of the sinusoidal texture [Fig. 2(a)]. Hence, the hexagonal texture of the substrate is not transferred to the perovskite top surface.
As we now have established that spin coating allows for the fabrication of perovskite films of a high morphological quality on sinusoidally nanotextured substrates, we can proceed with the optical simulations of nanotextured perovskite–silicon tandem solar cells.
Simulated Tandem Cell Architectures
For this study, we numerically investigated the three architectures shown in Fig. 3: Figure 3(a) shows the architecture with the “double-side” textured perovskite top cell. In Fig. 3(b), the “front side” architecture is shown, where the front interface of the perovskite absorber and all layers above are textured. In the “rear-side” textured architecture, Fig. 3(c), only the intermediate layers between the Si absorber and perovskite absorber are textured, whereas the interfaces on top of the perovskite layers remain flat. In all structures we assume conformal growth, meaning that all nanotextured interfaces have the same texture. For all architectures, we use a period of , which is very well suited for solar cells.22,23,40
As we have seen in Sec. 2, spin-coating allows for the growth of perovskite layers on nanotextured substrates. These layers are almost flat on top and hence resemble the architecture with the rear-side texture shown in Fig. 3(c). As demonstrated by other groups, a front-side textured architecture can be realized with nanoimprint processes.2627.–28 Paetzold and coworkers demonstrated that pattering the front surface does not affect the open circuit voltage and fill factors of solar cells.26 One strategy to realize double-sided architectures is using vacuum deposition, which can lead to conformal growth of the perovskite layer.41
In detail, the perovskite top cells shown in Fig. 3 consist of a lithium fluoride (LiF) antireflective coating, an ITO front contact, a tin oxide () buffer layer, an electron-selective phenyl-C61-butyric acid methyl ester (PCBM) contact layer, the perovskite absorber, a hole-selective nickel oxide (NiO) layer, and the ITO interconnecting layer.12,13 The silicon bottom cell consists of -doped and intrinsic (i) hydrogenated amorphous silicon (a-Si:H) layers and the c-Si wafer absorber, which is assumed to be infinitely thick in the finite-element method (FEM) simulations, as explained in Sec. 3.3. All other layers take their thicknesses from previous work on optimizing a planar tandem device42 and are shown in Table 1.
The thicknesses of the layers in our simulated tandem cell. These are the set of optimal layer thicknesses numerically calculated in Ref. 42. The perovskite thickness has been adapted to take into account the different (n,k) data used in this work. The air and c-Si layer are treated as infinite half spaces in the simulation. All layer thicknesses except the one of the perovskite are kept constant during the simulation. All values are in nm.
|LiF||ITO||SnO2||PCBM||Perovskite||NiO||ITO||n+ a-Si:H||i a-Si:H|
Mathematical Description of Hexagonal Lattices
The hexagonal sinusoidal nanotextures are mathematically described with22, where the texture described in Eq. (1) is referred to as “negative cosine.” For the solar cells discussed in this work, we also tested the “positive cosine” textures and just as in Ref. 22, they are outperformed by the negative cosine textures. Therefore, we only report about solar cells with negative cosine textures in this article.
We can scale the nanotexture to the desired period , which is the side length of the rhombus-shaped unit cell, with the following substitutions:
The optical simulations are performed with the finite element method (FEM) solver JCMsuite,43 which provides a rigorous solution to Maxwell’s equations for a given structure. As shown in Fig. 3, the three-dimensional (3-D) structures are meshed using tetrahedral and prismoidal elements.
The complex refractive index spectra used for the simulations were determined as follows: perovskite data were retrieved using ellipsometry and transmittance/reflectance spectrophotometry.44 For 45 and the sputtered ITO layers, ellipsometry and the program RIGVM were used.46 Data for PCBM were extracted from reflectance/transmittance measurements with the method described in Refs. 47 and 48. The layers were deposited using plasma-enhanced atomic layer deposition and characterized with ellipsometry.49 Data for the RF-PECVD a-Si:H layers50 were extracted using SCOUT.51 For LiF52 and spiro-OMeTAD,8 we used data from the literature.
We set the side-length constraint of the 3-D elements (prisms and tetrahedrals) for the FEM simulations to , where is the incident wavelength and is the refractive index of the material. In this way, one grid was generated for every wavelength interval of 100-nm width. Employing higher polynomial degrees (between 3 and 6) rather than finer meshes in an adaptive -FEM implementation allows to achieve accurate results at optimal computational costs.53 A convergence study confirmed the accuracy of the pre-set accuracy value. To accurately reconstruct the sinusoidal interfaces, we set their surface side-length constraint to between 55 and 17 nm, depending on the aspect ratio. To obtain the angular dependence of the hexagonal unit cell, the azimuth of the incident wave is set to 0 deg and 30 deg and an average is taken.22
We model the c-Si absorber and the air layer as infinite half spaces, which is numerically realized using perfectly matched layers on the top and bottom boundaries. For the other faces of the computational domain, periodic boundary conditions are applied. Light is incident from the air half space.
To cross check the accuracy of the simulations, we compared simulations for the planar tandem device with JCMsuite (3-D FEM) and GenPro4, which can combine wave optics for coherent thin layers and ray optics for thick incoherent layers.54 As can be seen in Fig. 4, the results obtained with the two methods are in excellent agreement.
However, the Si absorber of our tandem cell has a realistic thickness of . To account for this, we approximate the absorption in the Si absorber by assuming that the rear side of the Si behaves as a Lambertian scatterer. In practice this means that the wavelength-dependent absorbtion coefficient of the infinite Si layer —which is obtained using the light transmitted into the Si in the FEM simulation—is converted to the absorption of a finite layer with Lambertian scattering using the the Tiedje–Yablonovitch limit55
As a figure of merit, we use the photocurrent density , which can be calculated for the ’th layer from the absorption profile using56 and is the angle between the incident light and the solar cell normal. The photocurrent densities calculated for the perovskite and silicon layers are the maximum achievable current densities because we assume that all the absorbed light leads to the generation of electron-hole pairs, which are then extracted from the solar cell. The current densities for the other layers are losses due to parasitic absorption.
To estimate the potential power conversion efficiency of the tandem device, the open circuit voltage () of each subcell and the fill factor (FF) of the overall device are required. Additionally, the logarithmic dependence of the on the photocurrent density must be taken into account. We assume the fill factor to be 81%.57 Furthermore, we assume that a single junction perovskite and Si solar cell under standard illumination conditions have a short circuit current density of 22 and with a related open circuit voltage of 1.130 and 0.730 V, respectively. These values will then be used to estimate the values for the subcells in the tandem device.
To determine the perovskite thickness for current matching, we use Newton’s method. The method finds the roots of a function with derivative beginning at an inital value and repeats the iteration58
As inital value for the architecture with the smallest aspect ratio, we choose a perovskite thickness of 333 nm, which is the optimal value for the planar architecture.42 For each new value of aspect ratio, we use the optimized value for the previous aspect ratio as initial value. By employing this method, we reduce the number of simulations because no extensive thickness parameter scans are required. An example is shown in Fig. 5, where we see the current densities in the perovskite and silicon subcells after each iteration for the back side textured structure with aspect ratio . To reach current matching, the perovskite thickness had to be increased from 412 to 478 nm.
Figure 6 shows the matched maximum achievable current density as a function of the aspect ratio for the three architectures shown in Fig. 3. For all three textures, the largest aspect ratio leads to the maximum current densities, which are 21.3, 20.9, and for the double-, front-, and back-sided textures, respectively. For the double- and back-side textures that was , while for the front-side texture that was . An aspect ratio of 1.0 was not achievable with the front-sided texture because the layer thickness required for current matching would be smaller than the texture height. As discussed in Sec. 2, the architecture with the back-side textured solar cell (Fig. 6, red lines) is closest to the experimentally realizable architecture. Note that in this case the difference in the current densities between , 0.4, and 0.5 is .
The effective perovskite thicknesses at current match conditions are shown in Table 2; is the thickness a flat perovskite layer would have with the same volume and hence the same amount of material. The values are obtained by dividing the volumes of the perovskite absorber by the surface of the base. Comparing the effective thicknesses with the associated current densities we can see that a high effective thickness does not necessarily relate to a high absorption rate. For the architecture where only the layers above the perovskite absorber are textured (front-side), we find the highest proportion of absorption rate in the perovskite layer per effective thickness for all aspect ratios.
The effective perovskite thicknesses of the three different architectures illustrated in Fig. 3 with various aspect ratios after thickness optimization for current matching. All values are in nm.
Starting from the photocurrent density and of single junction cells given in Sec. 3.3, the matched photocurrent density values for the planar tandem of 19.7 and result in values of 1.127 and 0.711 V for the perovskite and Si subcells, respectively. This leads to a power conversion efficiency of 29.3%. The matched photo-current density for a device with a double-sided texture with is , which results in values of 1.129 and 0.713 V for the perovskite and Si subcells, respectively, leading to a power conversion efficiency of 31.8%. This represents an increase in 2.5% efficiency absolute at the standard perovskite bandgap of 1.56 eV. For the realistic back-textured architecture, the matched current density is at , leading to values of 1.128 and 0.711 V for the perovskite and Si cells, respectively. That results in an overall power conversion efficiency of 30.2%.
Figure 7 shows the absorption profile for silicon perovskite solar cells with (a) a planar and (b) a fully textured perovskite top cell with aspect ratio . Additionally, the corresponding current densities are shown. Note that only the current densities from the perovskite and silicon layers can contribute to electricity generation, whereas the other values are parasitic absorption or reflection losses. For silicon, we assume Lambertian scattering at the back using Eq. (3). Compared with the planar architecture, an additional current density of about per subcell can be collected by implementing the double-sided texture with .
We see that the planar device suffers from large reflective losses in the wavelength range from 400 to 1100 nm. These losses oscillate in magnitude, which is typical for Fabry–Perot type reflections from a planar thin-film layer stack. In contrast, the fully textured device shows negligible reflection losses in the same wavelength region. Hence, the overall reflection losses are reduced by . The sinusoidal texture scatters the light, thereby stopping multiple reflections necessary for Fabry–Perot oscillations, which is evident from the reflection curve being largely independent of wavelength. At the same time, the texture reduces the total reflection, which directly increases the absorption in the absorber layers.
The angular-dependent performance is also an important factor for solar cells without active tracking. Figure 8 shows how the maximum matched current density varies with the angle of incidence for both the optimized planar device and the back textured device with aspect ratio of 0.5. As the incident power falls off as , where is the incident angle, we divide the maximum matched current density by this factor to highlight the angular-dependent difference in light absorption between the planar- and back-textured device, independent of the incident power. It can be seen that the advantage in absorption and therefore current density provided by the nanotexture is maintained over a wide-angle range up to 60 deg. This supports the application of sinusoidal nanotextures for solar applications.
In this work, we numerically studied nanophotonic light management concepts for perovskite–silicon tandem solar cells. As nanophotonic structures we used hexagonal sinusoidal nanotextures. To achieve current matching, we optimized the thickness of the perovskite absorber with Newton’s method, which typically converged after only a single step. To experimentally motivate this numerical study, we demonstrated the deposition of perovskite layers onto sinusoidally structured nanophotonic substrates with spin coating.
The highest maximum achievable photocurrent density in each subcell was obtained with the fully textured structure with aspect ratio of , which corresponds to a power conversion efficiency of 31.8% at the standard perovskite bandgap. Compared with the planar architecture, we have an increase in of for each absorber layer and a potential increase of power conversion efficiency of 2.5% absolute. The efficiency can be increased even further using perovskite materials with higher bandgaps, where the optimal bandgap is around 1.7 eV.13,57
Future experimental efforts should focus on the implementation of double-side textured perovskite layers onto high-aspect-ratio substrates as well as to the experimental realization of textured perovskite single-junctions. Special care has to be taken to investigate the implications of structured substrates for perovskite crystal quality and thus electrical performance parameters were neglected in this study. Ultimately, future efforts should be directed to fully textured and highly efficient perovskite–silicon tandem solar cells.
The PTAA layer for the perovskite reference sample is prepared by spin-coating of a 2-mg PTAA (Aldrich) per mL of Toluene (Aldrich) solution onto an ITO-coated substrate (Automatic Research). After spin-coating at 4000 rpm for 30 s, the sample is annealed at 100°C for 10 min.
For the mixed cation, mixed halide perovskite, two 1.5-M precursor solutions of (TCI) and (TCI) in dimethylformamide (DMF, Aldrich) and dimethyl sulfoxide (DMSO, Aldrich) ( ) are mixed with formamidinium iodide (FAI, Dyenamo) and methylammonium bromide (MABr, Dyenamo), respectively. The resulting and solutions with a Pb-to-cation ratio of are mixed in a ratio and supplemented by addition of 5 vol% of a 1.5 cesium iodide (CsI, abcr) in DMSO solution. This solution is then diluted down to a 0.8-M solution to obtain films in the desired thickness range. For this, DMF:DMSO is added to the as-prepared solution.
For fabricating the perovskite layers, the antisolvent method was used. For this, the perovskite precursor solution was spin-coated onto the different substrates at 1000 and 6000 rpm for 10 and 20 s, respectively. About 5 s prior to the end of process, of chlorobenzene (Aldrich) was dripped as antisolvent. After this, the sample was directly transferred to a hot plate and baked at 100°C for 1 h.
The numerical results were obtained at the Berlin Joint Lab for Optical Simulations for Energy Research (BerOSE) of Helmholtz-Zentrum Berlin für Materialien und Energie, Zuse Institute Berlin and Freie Universität Berlin. D.E., G.K., C.B., and K.J. acknowledge the German Federal Ministry of Education and Research (BMBF) for funding the research activities of the Nano-SIPPE group within the program NanoMatFutur (Grant No. 03X5520). P.M. is funded by the Helmholtz Innovation Lab HySPRINT, which is financially supported by the Helmholtz Association. M.H. and S.B. acknowledge support by Einstein Foundation Berlin through ECMath within subproject SE6. S.A. acknowledges the BMBF within the project “Materialforschung für die Energiewende” for funding of his Young Investigator Group (Grant No. 03SF0540); together with P.T. he acknowledges the German Federal Ministry for Economic Affairs and Energy (BMWi) for funding of the “PersiST” project (Grant No. 0324037C). We thank Johannes Sutter for conducting the AFM measurements and Florian Ruske for supporting us with the determination of optical properties. The authors have no relevant financial interests in the article and no other potential conflicts of interest to disclose.
Duote Chen studied physics at the Freie Universität Berlin from where she received a bachelor’s degree in 2016. Currently, she studies medical physics at the Beuth University of Applied Sciences Berlin. She conducted the work presented in this paper as a research assistant at Helmholtz–Zentrum Berlin.
Phillip Manley completed a master’s degree in physics at Durham University in 2012. He went on to complete a doctorate in physics at the Freie Universität Berlin in cooperation with the Helmholtz–Zentrum Berlin in 2016 on the topic of nanoplasmonics for solar energy. He now works at the Helmholtz–Zentrum Berlin continuing to work on light management for solar energy.
Philipp Tockhorn is currently a PhD student at Helmholtz–Zentrum Berlin, working on the optimization of perovskite solar cells for the application in monolithic perovskite–silicon tandem solar cells. Before that, he studied physics at the Freie Universität Berlin and at the University of Amsterdam with a focus on photovoltaics and other renewable energy technologies.
David Eisenhauer studied physics at Heidelberg University and Technische Universität Berlin. He is currently a doctoral researcher at Helmholtz–Zentrum Berlin, where his research focusses on the development of tailored nanostructures for thin-film silicon solar cells and perovskite–silicon tandem solar cells.
Grit Köppel studied natural sciences in the Information Society at the Technische Universität Berlin and material science and engineering at Kiel University, where she completed a master’s degree in 2013. In 2017, she received his PhD at Technische Universität Berlin for her work on “Liquid phase crystallized silicon on sinusoidal textured glass substrates,” which she conducted at Helmholtz–Zentrum Berlin.
Martin Hammerschmidt received his Diploma in applied mathematics from the Technical University of Munich in 2011 and his PhD in the same field from the Freie Universität Berlin in 2016. His research interests include advanced finite element techniques, model order reduction methods and their applications in rigorous nano-optical simulations of photonic devices with a focus on energy harvesting devices.
Sven Burger leads the Computational Nano Optics Group at Zuse Institute Berlin. He also holds a position at JCMwave GmbH. He graduated in physics at the Leibniz Universität Hannover, where he also received a PhD in experimental physics. He further was a postdoctoral researcher at the European Laboratory for Nonlinear Spectroscopy in Florence and a lecturer at Technische Universität Berlin.
Steve Albrecht leads a young investigator group for perovskite-based multijunction photovoltaics at Helmholtz–Zentrum Berlin. He received a PhD in physics from the University of Potsdam. For his work on organic solar cells he was awarded with the Carl–Ramsauer–Prize and the Young Researcher Prize of the Berlin Physical Society and the Leibniz–Kolleg Potsdam, respectively. In 2014, he joined at Helmholtz–Zentrum Berlin as a postdoc until his young investigator group was founded.
Christiane Becker studied physics at the University of Karlsruhe and Heidelberg University, from where she obtained her diploma (MSc) in physics with distinction. After graduation as a PhD from University of Karlsruhe in 2006, she joined Helmholtz–Zentrum Berlin. Since 2012 she is head of the Nanostructured Silicon Group. She received an appointment as a professor at the University of Applied Sciences HTW Berlin in 2014. Her research primarily focusses on silicon nanostructures for photovoltaic and photonic applications.
Klaus Jäger received his diploma (MSc) as a physicist from the Swiss Federal Institute of Technology Zurich in 2008. In 2012, he graduated as a PhD with distinction from Delft University of Technology. After working in industry and a research position at Delft University of Technology, currently, he is a senior scientist at the Helmholtz–Zentrum Berlin. His research activities lie primarily with optical simulations for solar energy.