A variety of light-trapping strategies capitalizing on structures engineered on the order of the wavelength of solar light have been examined both experimentally and theoretically to enhance the efficiencies of solar cells.1 These strategies include texturing the front surface of the solar cell,2,3 incorporating a periodically corrugated metallic back reflector,45.–6 coating the front surface with an “antireflection” (AR) layer,7,8 embedding nanoparticles inside the light-absorbing layer(s),910.–11 and using light concentrators.1213.14.–15
An attractive approach for boosting the solar-cell efficiency requires the use of photon-absorbing component materials whose electromagnetic properties are periodically nonhomogeneous along the thickness direction. Recently, it was demonstrated that the incorporation of a periodically nonhomogeneous intrinsic layer (i.e., layer), along with a periodically corrugated back reflector, in amorphous silicon junction solar cells can improve overall efficiency by up to 17%.16 This improvement is likely due to
i. the periodic corrugation of the metallic back reflector45.–6,17 facilitating the excitation of surface-plasmon-polariton waves1819.–20 and waveguide modes21 to intensify the electric field inside the semiconductor region, leading to an increase in the electron–hole pair (EHP) generation rate;
ii. the periodic nonhomogeneity of the -layer that may facilitate the excitation of multiple surface-plasmon-polariton waves22 and waveguide modes,17 thereby further boosting the EHP generation rate; and
The study described herein concerns an especially simple type of solar cell that has received scant attention from researchers to date: a Schottky-barrier solar cell.25 The Schottky barrier is provided by a metal–semiconductor junction.26,27 Recent theoretical studies,28,29 along with an earlier experimental study,30 have suggested that Schottky-barrier solar cells may be a particularly promising proposition if the semiconductor layer (i.e., the absorbing layer) was made from alloys of indium gallium nitride (), since the bandgap for these alloys can closely match the range of energies of photons across the entire solar spectrum (i.e., 0.70 to 3.42 eV) by varying the relative proportions of indium and gallium through the parameter .31 Specifically, indium nitride (i.e., ) has a bandgap of 0.7 eV32,33 and absorbs efficiently across the infrared regime in the solar spectrum, while gallium nitride (i.e., ) has a bandgap of 3.42 eV and absorbs efficiently across the near-ultraviolet portion of the solar spectrum.
With current technologies, there are significant challenges to overcome in the production of alloys for all values of , particularly for smaller-bandgap alloys. When the proportion of indium is large (i.e., ), poor crystal growth plagues the realization of solar-cell applications. Such poor growth results in decreased carrier transport, background doping due to Fermi pinning above the conduction band edge,34 and a bandgap that is greater than expected.35,36 By using a Schottky-barrier junction with -doped , as opposed to a junction, the difficulty of doping the material is avoided.
In the following sections, a two-dimensional (2-D) numerical simulation is described for a Schottky-barrier solar cell made with . It is essential that the coupled processes of optical absorption and electrical-current generation are simultaneously accommodated in the simulation. This is because a focus on the computation of only the optical (i.e., at full quantum efficiency) short-circuit current density, but not of the open-circuit voltage, overplays the importance of the EHP generation rate by not taking the EHP recombination rate into account.16 In our simulation, particular attention is paid to the role of periodic nonhomogeneity of the absorbing material. Thus, the effect of periodically varying the fractional composition parameter of in the direction perpendicular to the mean plane of the back reflector is explored, the back reflector being periodically corrugated in one direction. The model is briefly described in Secs. 2.1 and 2.2, with further details being available elsewhere.16,37 Numerical results are presented in Secs. 3.1 and 3.2. Closing remarks are presented in Sec. 4.
Summary of the Model
A schematic diagram of the simulated Schottky-barrier solar cell is shown in Fig. 1. The back reflector is corrugated along the -axis, and the axis is normal to the plane of the Schottky contact and the mean plane of the back reflector. In the remainder of this paper, the term “width” refers to extent in the -direction, while “thickness” refers to extent in the -direction. Insolation is provided via an excitation port at . A planar AR window made of an insulating material occupies the region . To align with our earlier work, the optical permittivity of this layer was taken to be identical to that of aluminum-doped zinc oxide.38
The region is occupied by -doped , forming a Schottky junction and two ohmic junctions with a metal in the region . For calculations, the metal is assumed to be silver.39 The Schottky contact of width is centered at . The two ohmic contacts, each of width , are centered about . Note that . The two regions between the metal contacts for are occupied with the same material as the AR window that occupies .
A silver back reflector, covered with an insulating window (made of the same material as the window that occupies ), occupies the region . This back reflector is periodically corrugated in the -direction with period . The region is filled with the window material, while the region is filled with silver. The corrugation in the unit cell is specified by the function
Absorption of the normally incident solar flux with AM1.5G spectrum40 is calculated by solving the frequency-domain Maxwell postulates. The semiconductor charge-carrier drift-diffusion equations model the electron and hole density spatial distributions.41,42 Because of the nonhomogeneity of the semiconductor (i.e., ), the effective dc electric field acting on (a) electrons includes a contribution from gradients in the electron affinity and (b) holes includes contributions from gradients in both the electron affinity and the bandgap. Direct, mid-gap Shockley–Read–Hall, and Auger recombination are all included in our simulation. The current density , which is averaged over the Schottky contact (or, identically, both of the ohmic contacts), is calculated for a range of external biasing voltages .
By varying the proportion of indium relative to that of gallium, the bandgap of can be engineered to take any value from (i.e., for InN when ) continuously through to (i.e., for GaN when ). Thus, allowing for nonhomogeneity in the -direction, the alloy has bandgap given by43 To estimate for a material with a specific bandgap , Eq. (2) may be solved as
The electron affinity for is modeled in an analogous manner to the bandgap. HenceTable 1, Vegard’s law of linear interpolation44 is assumed to apply, i.e., Table 1.
Electronic data used for GaN and InN. The composition of InξGa1−ξN was estimated using Eq. (3), with linear interpolation used to estimate data for the semiconductor-filled region 0<z<Lz with bandgaps not presented here in all cases, except for the electron affinity χ0 which uses Eq. (4).
|Density of states (conduction band)|
|Density of states (valence band)|
|Electron mobility 1||295||1030|
|Electron mobility 2||1460||14150|
|Caughey–Thomas doping power (electrons)||—||0.71||0.6959|
|Caughey–Thomas critical doping density (electrons)|
|Hole mobility 1||3||3|
|Hole mobility 2||170||340|
|Caughey–Thomas doping power (holes)||—||2||2|
|Caughey–Thomas critical doping density (holes)|
|Auger recombination factor (electrons)|
|Auger recombination factor (holes)|
|Direct recombination factor|
|Slotboom reference energy||eV|
|Slotboom reference concentration|
|Adachi refractive-index parameter||—||9.31||13.55|
|Adachi refractive-index parameter||—||3.03||2.05|
The narrowing of the bandgap associated with doping was incorporated through the Slotboom model.45 While this model was developed for silicon, similar narrowing behavior under heavy-doping conditions has been observed in GaN.46 When doped, the bandgap of the semiconductor narrows to , while the electron affinity reduces to , with the conduction-band fraction being a material-specific parameter. Here, the bandgap narrowing is estimated as
The real part of the optical refractive index of is29 where and are parameters given in Table 1. The photon energy is , where is the Planck constant, is the speed of light in free space, and is the free-space wavelength. The imaginary part of the refractive index is 29 47
An empirical low-field mobility model—called either the Caughey–Thomas 29 or the Arora 48 mobility model—describes the variations of the electron mobility and the hole mobility with temperature and doping. Thus29 in lieu of experimental data, we assumed that for all simulations reported here. With these assumptions, the carrier mobilities simplify to Table 1 were provided by Hamady et al.,29 who also demonstrated that a relatively large metal work-function for the Schottky-barrier contact improves efficiency. Accordingly, the relatively large value of was chosen here.
The bandgap profile of for the solar cells simulated in this study is given byFig. 2 for , when and .
A 2-D finite-element optoelectronic model was implemented in the COMSOL Multiphysics (V5.1) software package48 in two major steps, as described now. Let it be noted that terms in block capitals are COMSOL Multiphysics48 terms.
In the first major step, the Electromagnetic Waves, Frequency Domain module was used to calculate the 2-D EHP generation rate as a function of and . Normally incident monochromatic light sampled at 5-nm intervals on the -scale across the AM1.5G spectrum, 50% polarized and 50% polarized, was activated by the Periodic ports option at . Diffraction Order ports for diffraction orders were added. The inclusion of diffraction ports for even higher orders has little impact on the resulting EHP generation rate, predominantly due to strong absorption of shorter wavelength photons by , the majority of which are absorbed close to the surface of the device before they can be scattered. The boundaries parallel to the -axis have Floquet periodicity with the wavevector provided by the periodic ports. The region behind the back reflector was taken to be a perfectly matched layer. The semiconductor region was first meshed with a Mapped mesh with a 10-nm Distribution in both the - and -directions and then split into to a triangular distribution using the insert center points conversion. The back-reflector region was covered with an Extra fine, Delaunay, and Free Triangular mesh. Data for a spatial map of the computed EHP generation rate were stored in an external file.
The EHP generation rate can be used to compute the optical short-circuit current density , assuming that every absorbed photon creates an EHP in the layer and that no recombination takes place. Neglect of recombination implies that is necessarily larger than the short-circuit current density , which is the electronically simulated current density that flows when the solar cell is illuminated and no external bias is applied (i.e., when ).
In the second major step, the Semiconductor module was used to calculate the electron and hole densities in the semiconductor region, and thereby the current densities. Due to the symmetry of the unit cell, only its right half (i.e., ) needs to be electrically simulated. Fermi-Dirac carrier statistics, along with continuous quasi-Fermi levels at any internal boundary, were employed. Finite volume (constant shape function) discretization was employed as this inherently conserves current throughout the solar cell.48 COMSOL uses a Scharfetter–Gummel upwinding scheme for solving the charge carrier transport equations. The Free triangular, Delaunay mesh has a maximum element size of 15 nm.
A potential difference of was applied between the ohmic and Ideal Schottky contacts. Thermionic currents, with standard Richardson coefficients of , and were applied at the Schottky barrier.29,48 Insulator Interfaces were placed at the remaining external electrical boundaries. External Mathematica™ or MATLAB™ codes were used to calculate the User-Defined Generation from the output of the first major step. Recombination was incorporated via Auger, Direct and Trap-Assisted (Midgap Shockley-Read-Hall) pathways, with parameters as provided in Table 1. To facilitate convergence, the nonhomogeneity, EHP generation, and EHP recombination physics were slowly activated as the solver progressed by use of a continuation parameter. An analytic doping model was used to set the donor concentration to . The Free triangular, Delaunay mesh used has a maximum element size of 15 nm.
Numerical Simulations of Schottky-Barrier Solar Cells
Design of Periodically Corrugated Back Reflector
The numerical results presented here are for devices with parameters listed in Table 2. A preliminary study was carried out to ascertain reasonable values for the dimensions of the device. Figure 3 presents as a function of the back-reflector period for a Schottky-barrier solar cell with a homogeneous layer of thickness . For this set of simulations, , , and were determined as optimal for a solar cell with a homogeneous layer characterized by , , and . While this optimization is far from exhaustive, significant optical variation was seen to occur for thinner solar cells: a 6% variation in was calculated for a solar cell with , with a maximum of when . Importantly, for solar cells with layers thicker than about 600 nm, saturates as increases beyond a threshold value. Therefore, was chosen to be optimal. While the peak value of of a Schottky-barrier solar cell with a 600-nm-thick layer is obtained for larger periods than this (), falls more rapidly for devices with back reflectors of overly large periods than for those with periods that are too small.
Summary of parameters used for the simulations.
|Device width||600 nm|
|-doped thickness||600 nm|
|Back reflector insulating layer thickness||135 nm|
|Corrugation height||130 nm|
|Corrugation duty cycle||0.5|
|Metal thickness||100 nm|
|Contact region thickness||5 nm|
|Ohmic contact width||140 nm|
|Schottky contact width||440 nm|
|AR window thickness||100 nm|
|Air thickness||500 nm|
|Predoping bandgap||[0.7,3.42] eV|
|Bandgap-nonhomogeneity shaping parameter||(0,32)|
Effect of Periodic Nonhomogeneity of Layer
The essential effects of including periodic nonhomogeneity in the thickness direction in the layer may be inferred from Fig. 4. For a fixed bandgap-nonhomogeneity amplitude , the baseline bandgap is shown to dramatically affect the simulated efficiencyFig. 4, it may be inferred that the difference in the efficiencies with a homogeneous layer and a periodically nonhomogeneous layer is greatest when
Indeed, at , the efficiency is 13.26% for and 16.81% for , i.e., the relative increase in efficiency attributable to the periodic nonhomogeneity of the layer is 26.8%, which is this paper’s most significant result. For , is optimal, yielding an efficiency of 15.8%.
The importance of the periodicity in the bandgap variation when seeking maximal efficiency may be inferred from Figs. 5 and 6. In Fig. 5, the thickness of the layer is fixed at while the ratio varies from 0.5 to 4. The efficiency of the solar cell is seen to reach a maximum value of about 17% at and . Between these two values of , only a slight decrease in was found. Significantly, the efficiency for all values of is greater than it is when the layer is homogeneous.
Figure 6 illustrates the effect of varying the thickness of the layer when the nonhomogeneity period is fixed. Irrespective of the period, the solar cell with the nonhomogeneous layer performs better than the one with a homogeneous layer. Furthermore, the solar cell with the smaller period () is found to perform more efficiently than the one with larger period (). Distinct maximums can be seen at and , which correspond to and , respectively. The efficiencies at these points, respectively, are 17.28% and 16.88%, which correspond to relative increases in efficiency of 16.7% and 28.5% as compared to the analogous solar cell with a homogeneous layer.
By varying the bandgap-nonhomogeneity shaping parameter , the spatial profile of the nonhomogeneity can be varied. Thus, holds for a sinusoidal bandgap profile, while smaller and larger values of yield bandgap profiles with steeper gradients, as illustrated in Fig. 2. For all three values of in Fig. 7, the efficiency of the solar cell is substantially greater than it is for the corresponding solar cell with a homogeneous layer. The maximum increase in efficiency is found for a nearly sinusoidal bandgap profile; indeed, the efficiency increases from 13.2% for a Schottky-barrier solar cell with a homogeneous layer to 16.88% for its analog with a nonhomogeneous layer with , which is a relative increase in efficiency of 28.5%.
Let us note here that as becomes much larger or smaller than unity, the spatial gradients of the bandgap increase in magnitude and the peaks in Fig. 2 become narrower. Accordingly, at extreme values of , the semiclassical carrier-transport equations implemented in our model become less appropriate as quantum processes become significant. Then, it may be necessary to take into account the quantization of allowed energy states between the peaks of and tunneling through those peaks, especially for large amplitudes .
A 2-D finite-element model was devised to simulate the combined optical and electrical performances of Schottky-barrier solar cells. First, it was found that a periodically corrugated back reflector of a period of 600 nm is optimal for photon absorption in a solar cell containing a homogeneous layer of . Second, the effects of periodic nonhomogeneity of the layer were elucidated. The nonhomogeneity was directed perpendicular to the mean plane of the periodically corrugated back reflector. For the particular model investigated here, the efficiency of a solar cell with a 600-nm-thick layer of was found to increase by 26.8% when suitable periodic nonhomogeneity was incorporated. Thus, the incorporation of a periodically nonhomogeneous layer in a Schottky-barrier solar cell can substantially increase the efficiency as compared to the analogous solar cell with a homogeneous layer. A comprehensive optimization of study of material and design parameters may well yield even greater nonhomogeneity-induced increases in efficiency. However, such an optimization study—which is justified by the substantial efficiency boosts reported herein—represents a major undertaking that lies beyond the scope of this paper. Similarly, it would be of some value to delineate the photon-absorption mechanisms that underpin the boost in the light-to-electricity conversion efficiency that arises following the introduction of periodic nonhomogeneity, but this too is matter for future study.
The feasibility of producing alloys with the prescribed spatial variation in to achieve efficiency boosts in Schottky-barrier solar cells is a matter for our experimentalist colleagues to shed light on. On the basis of our numerical study and the relatively huge boosts in efficiency that may be attained, it would be worthwhile for major efforts to be directed toward the production of appropriate alloys.
Last, the principal finding of our study is that the introduction of periodic nonhomogeneity can, in principle, substantially boost the efficiency of a Schottky-barrier solar cell. To focus on efficiency, our attention was restricted to only normally incident solar radiation for all cell configurations considered. The influence of the angle of incidence is planned for a future study in which the mechanisms for efficiency boosting will be more fully explored.
T.H.A. thanks the Charles Godfrey Binder Endowment for partial financial support during a six-month stay at Pennsylvania State University. T.G.M. acknowledges the support of EPSRC Grant No. EP/M018075/1. A.L. thanks the National Science Foundation for partial financial support under Grant No. DMS-1619901, and he is grateful to the Charles Godfrey Binder Endowment for ongoing support of his research.
Tom H. Anderson received his BSc degree in mathematics and physics from the University of Edinburgh in 2011. He received his MSc degree in fusion energy from the University of York in 2012 and his PhD in applied mathematics from the University of Edinburgh in 2016. His research interests include the optical and electrical modeling of thin-film solar cells.
Tom G. Mackay is a reader at the School of Mathematics, the University of Edinburgh, and an adjunct professor in the Department of Engineering Science and Mechanics at Pennsylvania State University. He is a graduate of the Universities of Edinburgh, Glasgow, and Strathclyde and a fellow of the Institute of Physics (UK) and SPIE. His research interests include the electromagnetic theory of complex materials, including homogenized composite materials.
Akhlesh Lakhtakia received degrees from the Banaras Hindu University and the University of Utah. He is the Charles Godfrey Binder professor of engineering at the Pennsylvania State University. His research interests include surface multiplasmonics, bioreplication, forensic science, solar cells, sculptured thin films, and mimumes. He is a fellow of Optical Society of America, SPIE, Institute of Physics, American Association for the Advancement of Science, American Physical Society, Institute of Electrical and Electronics Engineers, and Royal Society of Chemistry. He received the 2010 SPIE Technical Achievement Award and the 2016 Walston Chubb Award for Innovation.