Expand this Topic clickable element to expand a topic
Skip to content
Optica Publishing Group

Optical center of a luminescent solar concentrator

Open Access Open Access

Abstract

This Letter introduces a novel, to the best of our knowledge, approach to estimate the power conversion efficiency (PCE) of a square luminescent solar concentrator (LSC) by point excitations on the “optical centers.” Predicted by theoretical calculations, photoluminescence emissions from these optical centers experience almost the same average optical path as those from the whole device under uniform illumination. This is experimentally verified by a 20 × 20 cm2 silicon quantum dot-based LSC, with a negligible error between the predicted PCE and the measured one. This method provides a convenient way to estimate the photovoltaic performance of large-area LSC devices with basic laboratory instruments.

© 2022 Optica Publishing Group under the terms of the Optica Open Access Publishing Agreement

A planar waveguide with uniformly embedded light sources is a generic concept used in laser phosphor displays [1], solar-pumped lasers [2], and luminescent solar concentrators (LSCs). The latter were first proposed in the early 1970s [3], with the main purpose to concentrate sunlight by means of photoluminescence (PL). The fluorophores embedded in the host matrices absorb a part of sunlight and re-emit light. A significant fraction of PL emission (∼75% for a slab with a refractive index of 1.5) is trapped in the waveguide by total internal reflection until it reaches the edges for conversion to electricity. In recent years, LSCs have gained much attention as semitransparent building-integrated photovoltaic (BIPV) units balancing PV performance and aesthetic quality, commonly known as “solar windows” [4].

The active light-converting materials are typically organic dyes [5], semiconductor quantum dots [4], rare-earth ions [6], or metal nanoclusters [7]. For efficient operation, a high photoluminescence quantum yield (PLQY) and a large Stokes shift are required to reduce reabsorption loss of the emitted light. The typical matrices are glass or polymers such as poly-(methyl-methacrylate) (PMMA) and off-stoichiometry thiol-ene (OSTE) [8]. A triplex laminate structure is also used [8], with a thin polymer interlayer sandwiched in between two glass pieces. The fluorophores should also be well dispersible inside, translating to low scattering and haze for better energy conversion efficiency and glazing function.

Over the years, efforts have been devoted to improving the power conversion efficiency (PCE) of LSCs [9]. In addition, extending the device size to large areas (≥100 cm2) becomes significant for practical applications. For the LSC PV performance characterization, following the conventional method requires a correspondingly large solar simulator. This approach, however, is not feasible when it comes to truly large LSCs (tens of cm) for in-the-lab device measurements. Recently, regional measurements were introduced to analyze the PV performance and the spectral response of large-area LSCs with a shifting mask and a piece of solar cell placed at different edge regions [10,11]. However, scanning numerous excitation points all over the large-area device may not always be practical either.

In this Letter, a new approach for estimating the PV performance of large-area LSCs from point excitations is introduced. These points are named as “optical centers” of an LSC in analogy with the center of mass in mechanics. It is defined as points from which the average optical path to the device edges is equal to that of the whole device under uniform excitation. Here, the response to point excitations for a planar waveguide is analytically derived, and the position of such optical centers is found to be at ∼82% of the center-edge distance for the square geometry. Using selective excitation to this spot, the short-circuit current ${I_{sc}}$ of the whole device can be obtained by linearly scaling with the excitation area. This hypothesis was then successfully verified by a set of measurements on a 20 × 20 cm2 LSC with silicon quantum dots (Si QDs) as fluorophores. It is shown that using the average ${I_{sc}}$ measured at the optical center, the PCE of this device can be very well predicted within a measurement error, confirming that this approach is practically useful for the experimental evaluation of the PV performance of large-area LSCs.

In this work, we consider square geometry, which is very common for LSCs. Previously, we found that uniform excitation square and circular geometries were very similar in terms of optical path distribution from an isotropic emitter [12]. Using this established similarity here, we can derive simple expressions for the average optical path for a point excitation.

Consider a square with side length 2a and a circle of the same area with radius R. These two quantities are related by $R = 2a/\sqrt \pi \approx 1.128a$. First, we derive an analytical formula for the average optical path inside a circle for an emitter placed at a distance $\rho$ from its center in 2D [Fig. 1(a), inset]. An optical ray emitted at an angle $\theta$ from the radial direction will experience path r until reaching the perimeter, where it is absorbed (all edges covered with solar cells). These two quantities are related as

$$\theta = \textrm{arccos}(\frac{{{R^2} - {\rho ^2} - {r^2}}}{{2\rho r}}). $$

The emitter is isotropic, so the angular distribution of its intensity $dI/d\theta$ is a constant. Then one can express emitted intensity distribution over the optical paths as

$$p(r) = \frac{{dI}}{{dr}} = \frac{{dI}}{{d\theta }}\frac{{d\theta }}{{dr}} = \frac{{{R^2} - {\rho ^2} + {r^2}}}{{\pi r\sqrt { - {\rho ^4} + 2({R^2} + {r^2}){\rho ^2} - {{({R^2} - {r^2})}^2}} }}, $$
where the constant was chosen to satisfy the normalization condition $\int_{R - \rho }^{R + \rho } {p(r)dr} = 1$, so that $p(r)$ has the meaning of a probability density function.

 figure: Fig. 1.

Fig. 1. (a) Average optical path as a function of the point excitation position. Square dots represent numerical simulation values for a square, while the line is for a circle of the same area from Eq. (3). Inset is a schematic illustration. (b) Waveguiding efficiencies for point (green, red, blue, orange solid lines) and uniform (black solid line) excitations of an LSC as a function of the dimensionless variable αmR. The vertical gray dashed line indicates the LSC device used in this work for the experimental validation.

Download Full Size | PDF

Then the average optical path is a first moment of this distribution. Integration yields

$${r_{avg}}(\rho \ast ) = \frac{{(1 + \rho \ast ) \cdot E(\gamma ) + (1 - \rho \ast ) \cdot K(\gamma )}}{\pi }R, $$
where $\rho \ast $ is the dimensionless distance from the center ($\rho \ast{=} \rho /R$), $\gamma = 2\sqrt {\rho \ast } /(1 + \rho \ast )$, and $K(\gamma )$, $E(\gamma )$ are complete elliptic integrals of the first and second kind, respectively. It is shown as a function of $\rho /a$ in Fig. 1(a) as a red line. It is seen that the average optical path shortens when the emitter approaches the absorbing perimeter, as would be expected.

To validate the applicability of this formula for a square geometry, we plotted numerically calculated average optical paths for an isotropic emitter placed at several points inside the same area square (details in the Supplementary material). Those are shown as square dots in the same units in Fig. 1(a), confirming similarity of these two geometries in this sense. So the first conclusion of this analysis is that the position of a point excitation in 2D geometry (x,y) for a square waveguide can be effectively reduced to the function of a single variable $\rho$.

Next, we postulate that there is a spot on the LSC (at a distance from the center ${\rho _0}$), where the average optical path for an isotropic emitter equals to that of the whole device under uniform illumination. The latter was derived in [12] for a circle as 8R/3π, therefore,

$${r_{avg}}({\rho _0}\ast ) = \frac{8}{{3\pi }}R. $$

The numerical solution of this equation can be found using Eq. (3). In the units of length, it is ${\rho _0} = 0.731R$. For a square with the side length $2a$, this is equivalent to

$${\rho _0} = 0.824a. $$

We assign these points to the optical centers of a square planar waveguide.

In general, however, the average optical path may not necessarily represent a valid parameter for light propagation assessment. For example, because the distribution of Eq. (2) is bimodal, the presence of attenuation may result in substantially different waveguiding efficiencies than those predicted using ${r_{avg}}$. Another necessary realistic aspect to consider is the 3D geometry of such planar waveguides.

Due to the total internal reflection, all photons emitted outside the escape cone (42° for a glass or plastic slab) will be waveguided to the device edge. The distribution of optical paths for such photons is therefore quite narrow and an average value for the out-of-plane emission can be used [13]. It was previously shown that all 2D optical paths can be simply extended by a factor $k \approx 1.144$ [12]. Changing variables ($l = kr$) from $p(r)$, the normalized probability density function $q(l)$ for the 3D case becomes

$$q(l) = \frac{{p(l/k)}}{k}. $$

Using this formula, we can compare the waveguiding efficiency for the single spot excitation with the total area illumination case in the presence of attenuation. Let the attenuation coefficient be $\alpha$ [cm−1], which in practice can represent a matrix absorption coefficient ($\alpha \approx {\alpha _m}$) when only considering the matrix absorption loss. Then the waveguiding efficiency for a spot excitation at position $\rho$ is

$${f_{sp}} = \int_{l\min }^{l\max } {\exp ( - {\alpha _m}L)q(l)dl}, $$
where ${l_{\min }} = k(R - \rho )$ and ${l_{\max }} = k(R + \rho )$. This function is shown for four different values of ρ (0, 0.8ρ0, ρ0, 1.2ρ0) in Fig. 1(b). Dimensionless variable ${\alpha _m}R$ is limited by the upper LSC critical size (${\alpha _m}R < 1.5$), as introduced in [12]. However, the expression for the waveguiding efficiency under uniform illumination was previously derived [12]:
$${f_{LSC}} ={-} \frac{{{M_1}(2{\alpha _m}Rk)}}{{{\alpha _m}Rk}}, $$
where ${M_1}$ is a modified Struve function of the second kind. The ${f_{LSC}}$ function is shown in Fig. 1(b) as a black line. The waveguiding efficiencies are shown to be nearly identical for point excitation to the optical center and for uniform excitation in the considered range (red and black lines). Thus, the concept of the LSC optical center is theoretically validated under realistic conditions of attenuation and 3D geometry. Note that for fluorophores with non-negligible reabsorption coefficients ${\alpha _{re}}$, this concept should still be applicable as long as the matrix absorption is the main loss mechanism (${\alpha _m} > {\alpha _{re}}$). In addition, for a certain device (given ${\alpha _m}R$), the further the excitation point is from the center will mean a higher waveguiding efficiency. Moreover, a larger device (larger ${\alpha _m}R$) will result in a larger waveguiding efficiency gap between the geometrical center and the optical center (or whole device). In practice, it implies that the photocurrent difference between exciting geometrical and optical centers becomes quite substantial for large-area devices. Consequently, it is reasonable to deduce the ${I_{sc}}$ value of the whole device from that of the optical center measurement.

To verify our hypothesis, a 20 × 20 cm2 (a = 10 cm, R ≈ 11.3 cm) low-haze LSC device based on luminescent Si QDs was used for the experimental validation. The device fabrication procedure can be found elsewhere [14]. In the triplex laminate structure, an off-stoichiometry thiol-ene (OSTE) polymer interlayer is of 3-mm thickness, and glass sheets are of 2-mm thickness each. Altogether ∼60 mg of near-infrared-emitting (PL peak position ∼870 nm, PLQY ∼50%) Si QDs are uniformly dispersed in the polymer interlayer. The transmittance, absorption, reflectance, and haze spectra of this device are shown in Fig. S1 (Supplement 1). Note that this device is highly transparent (visible light transparency of ∼90%) and thus not optimized for high power generation. Under the standard AM1.5G, ∼8.3% of solar power is absorbed for the ensuing conversion to electricity (Fig. S2, Supplement 1). As reported previously [12], the matrix absorption coefficient ${\alpha _m}$ is 0.04 cm−1 for OSTE and of similar value for the low-iron glass applied here. Therefore, the dimensionless variable for this device is ${\alpha _m}R \approx 0.45$ [marked as a dashed line in Fig. 1(b)]. Details of the attached solar cells and the solar simulator used for the test can be found in Supplement 1l.

For excitation position-dependent measurements, two directions (center-left and center-right) of the LSC were probed under a series of point excitations by shifting the 1 × 1 cm2 opening positions on the mask (details in Supplement 1). As shown in Fig. 2(a), a clear rising trend of ${I_{sc}}$ from the center to the edge along both directions was observed. An analytical curve [gray curve in Fig. 2(a)] can be drawn by considering the matrix absorption loss for the average optical path ${r_{avg}}$ under different point excitation positions [Fig. 1(a)] and the linearity of the short-circuit current with irradiance as ${I_{sc}}(\rho ) = {I_0}{e^{ - {\alpha _m}k{r_{avg}}(\rho )}}$. As can be seen, the experimental data match well with the predicted trend with proportionality coefficient ${I_0}$ being the only fitting parameter.

 figure: Fig. 2.

Fig. 2. (a) Short-circuit current under point excitations along two directions (center-left in blue and center-right in red) of a 20 × 20 cm2 LSC under 0.9 AM1.5G. The gray solid line indicates the predicted trend. Note that due to the fabrication procedure of layered polymerization, only the directions non-crossing interfacial optical defects were considered. (b) External quantum efficiency under the optical center illumination (red) and uniform illumination of the whole device (blue). Inset is a photograph of the 20 × 20 cm2 device.

Download Full Size | PDF

To independently verify that the observed trend indeed originates from the luminescence waveguiding and not from random stray light or direct light on solar cells, we performed external quantum efficiency (EQE) measurements under the optical center (8 cm from the geometrical center) and whole device illumination. For the LSC device, EQE(λ) is the ratio of generated electrons by the edge solar cells to incident photons on the device top area. To block the direct illumination of solar cells, four “roofs” (each 5 mm wide) were applied at the top of the device edges in the case of the whole device illumination (in the case of point excitation, those are redundant thanks to the mask). As shown in Fig. 2(b), the EQE curves are consistent with the shape of the absorption spectrum of Si QDs (Fig. S1, Supplement 1). More importantly, both EQE curves are essentially identical, indicating almost the same short-circuit current density (mA/cm2) after being integrated with the solar spectrum. Moreover, the EQEs are close to zero above 650 nm, implying almost no contribution from the direct illumination or stray light, validating previous results.

Verification of the signal origin being from QD PL is essential for correct evaluation and application of the optical center method. Another aspect to consider is the size of the solar cell modules, as discussed in detail in Supplement 1. Briefly, to correctly resolve optical path variations under the point excitation, small size solar cells units, connected in parallel, are necessary. Otherwise, the ${I_{sc}}$ of a solar cell will be limited by the lowest ${I_{sc}}$ among the in-series connected p-n junctions (longest optical path). An example is shown in Figs. S3 and S4 in Supplement 1, where four-fold longer solar cells units were attached to the device edges. In that case, the trend, which is clearly visible in Fig. 2(a), becomes disrupted when approaching the edges.

After the optical center concept has been experimentally verified, we are in the position to discuss practical implementations. The most straightforward application is the short-circuit current evaluation of the whole device from the optical center excitation. To minimize the impact from the solar simulator non-uniformity, each ${I_{sc\_oc}}$ value under the optical center excitation was measured at each orthogonal direction under the solar simulator by spinning both the device and the mask, yielding an average ${I_{sc\_oc}}$ value of ∼79 µA (±3 µA) under standard AM1.5G. Multiplying by the opening area ratio (400:1), one can predict ${I_{sc}}$ under uniform illumination to be ∼31.6 mA (±1.2 mA), which is very close to the measured value of ∼33 mA (±0.5 mA) by the conventional method, matching within measurement error (I-V curve shown in Fig. S5, left, Supplement 1). We attribute this small error to the measurements instead of the mathematical assumptions, most likely due to the slight non-uniformity of the irradiance from the solar simulator and the LSC device. This confirms the plausibility of single excitation spot measurements for the whole device evaluation in terms of photocurrents.

Next, we can extend the method of predictive ability to the power conversion efficiency (PCE), which is an important figure of merit for photovoltaic devices. To accomplish that, the open circuit voltage and the fill factor (FF) need to be estimated in addition to the short-circuit current. However, under optical center excitations, the ${V_{oc}}$ is limited by the lowest ${V_{oc}}$ among the solar cells connected in parallel (farthest solar cells from the open window), and we instead suggest two methods for the estimation of ${V_{oc}}$ and $FF$.

On average, for each piece of solar cell, ∼1 mA (32 pieces in parallel) of current can be coarsely estimated. This estimation does not introduce a large error since the open circuit voltage grows logarithmically to the incident irradiance (or to the short-circuit current, which is linear to irradiance) [15].

By tuning the irradiance under the solar simulator (either percentage of AM1.5G or monochromatic light), the I-V curve (Fig. S5, right, Supplement 1) of a single piece of solar cell was recorded to simulate device operation conditions delivering ${I_{sc\_s}}$ of ∼1 mA. Then the open circuit voltage was ∼2.01 V and the fill factor was ∼0.63, which are indeed almost the same as measured ones for the operating LSC (Fig. S5, left, Supplement 1). Alternatively, considering the known logarithmic relationship between the ${V_{oc}}$ of the solar cell and the incident irradiance similar values can be obtained, as detailed in the Supporting Information.

At last, the procedure of whole device PCE estimates from point excitation measurements can be summarized. (i) Measure the average ${I_{sc\_oc}}$ with only the window at the optical center open. Scaling with the whole device area gives the overall ${I_{sc}}$ under uniform illumination. (ii) Consider the solar cell connections and estimate the ${I_{sc\_s}}$ value for a single piece of solar cell. Place a single piece of solar cell under the solar simulator, and tune the irradiance until the solar cell delivers the same ${I_{sc\_s}}$. (iii) Extract the open circuit voltage ${V_{oc}}$ and fill factor $FF$ from the I-V curve under this irradiance. (iv) Apply these values to the whole device efficiency calculations: $PCE = {I_{sc}} \cdot {V_{oc}} \cdot FF/{P_{incident}}$.

To conclude, the concept of an optical center on a square LSC was proposed based on analytical derivations, where the average optical path of emitted photons from the optical center and the whole device are stipulated to be the same. It is found that the optical center is located at ∼82% of the half-length distance from the geometrical center of the planar waveguide. In practice, from the short-circuit current measured at the optical center and the ratio of illumination area, one can predict the short-circuit current of the whole device under uniform illumination. This hypothesis is successfully verified by a 20 × 20 cm2 Si QDs-based LSC both for ${I_{sc}}$ and PCE with a negligible error. This methodology can help access the PV performance of a large-area LSC device by point excitations with basic laboratory instruments.

Funding

Energimyndigheten (46360-2).

Acknowledgment

J.Z. acknowledges funding support from the China Scholarship Council (CSC).

Disclosures

The authors declare no conflicts of interest.

Data availability

Data underlying the results presented in this paper are not publicly available at this time but may be obtained from the authors upon reasonable request.

Supplemental document

See Supplement 1 for supporting content.

REFERENCES

1. I. Fujieda, S. Itaya, M. Ohta, Y. Hirai, and T. Kohmoto, J. Photonics Energy 7, 028001 (2017). [CrossRef]  

2. T. Masuda, M. Iyoda, Y. Yasumatsu, S. Dottermusch, I. A. Howard, B. S. Richards, J.-F. Bisson, and M. Endo, Commun. Phys. 3, 60 (2020). [CrossRef]  

3. A. Goetzberger and W. Greube, Appl. Phys. 14, 123 (1977). [CrossRef]  

4. F. Meinardi, S. Ehrenberg, L. Dhamo, F. Carulli, M. Mauri, F. Bruni, R. Simonutti, U. Kortshagen, and S. Brovelli, Nat. Photonics 11, 177 (2017). [CrossRef]  

5. A. Sanguineti, M. Sassi, R. Turrisi, R. Ruffo, G. Vaccaro, F. Meinardi, and L. Beverina, Chem. Commun. 49, 1618 (2013). [CrossRef]  

6. T. Wang, J. Zhang, W. Ma, Y. Luo, L. Wang, Z. Hu, W. Wu, X. Wang, G. Zou, and Q. Zhang, Sol. Energy 85, 2571 (2011). [CrossRef]  

7. K.-B. Cai, H.-Y. Huang, P.-W. Chen, X.-M. Wen, K. Li, J.-L. Shen, K.-P. Chiu, and C.-T. Yuan, Nanoscale 12, 10781 (2020). [CrossRef]  

8. J. Huang, J. Zhou, T. Haraldsson, A. Clemments, M. Fujii, H. Sugimoto, B. Xu, and I. Sychugov, Sol. RRL 4, 2000195 (2020). [CrossRef]  

9. S.-J. Ha, J.-H. Kang, D. H. Choi, S. K. Nam, E. Reichmanis, and J. H. Moon, ACS Photonics 5, 3621 (2018). [CrossRef]  

10. Y. Li, Y. Sun, and Y. Zhang, Renewable Energy 160, 127 (2020). [CrossRef]  

11. Y. Li, Y. Zhang, Y. Sun, and T. Ren, Appl. Opt. 59, 8964 (2020). [CrossRef]  

12. I. Sychugov, Appl. Opt. 59, 5715 (2020). [CrossRef]  

13. I. Sychugov, Optica 6, 1046 (2019). [CrossRef]  

14. J. Huang, J. Zhou, E. Jungstedt, A. Samanta, J. Linnros, L. A. Berglund, and I. Sychugov, ACS Photonics 9, 2499 (2022). [CrossRef]  

15. V. Lo Brano, A. Orioli, G. Ciulla, and A. Di Gangi, Sol. Energy Mater. Sol. Cells 94, 1358 (2010). [CrossRef]  

Supplementary Material (1)

NameDescription
Supplement 1       Supporting contents

Data availability

Data underlying the results presented in this paper are not publicly available at this time but may be obtained from the authors upon reasonable request.

Cited By

Optica participates in Crossref's Cited-By Linking service. Citing articles from Optica Publishing Group journals and other participating publishers are listed here.

Alert me when this article is cited.


Figures (2)

Fig. 1.
Fig. 1. (a) Average optical path as a function of the point excitation position. Square dots represent numerical simulation values for a square, while the line is for a circle of the same area from Eq. (3). Inset is a schematic illustration. (b) Waveguiding efficiencies for point (green, red, blue, orange solid lines) and uniform (black solid line) excitations of an LSC as a function of the dimensionless variable αmR. The vertical gray dashed line indicates the LSC device used in this work for the experimental validation.
Fig. 2.
Fig. 2. (a) Short-circuit current under point excitations along two directions (center-left in blue and center-right in red) of a 20 × 20 cm2 LSC under 0.9 AM1.5G. The gray solid line indicates the predicted trend. Note that due to the fabrication procedure of layered polymerization, only the directions non-crossing interfacial optical defects were considered. (b) External quantum efficiency under the optical center illumination (red) and uniform illumination of the whole device (blue). Inset is a photograph of the 20 × 20 cm2 device.

Equations (8)

Equations on this page are rendered with MathJax. Learn more.

θ=arccos(R2ρ2r22ρr).
p(r)=dIdr=dIdθdθdr=R2ρ2+r2πrρ4+2(R2+r2)ρ2(R2r2)2,
ravg(ρ)=(1+ρ)E(γ)+(1ρ)K(γ)πR,
ravg(ρ0)=83πR.
ρ0=0.824a.
q(l)=p(l/k)k.
fsp=lminlmaxexp(αmL)q(l)dl,
fLSC=M1(2αmRk)αmRk,
Select as filters


Select Topics Cancel
© Copyright 2024 | Optica Publishing Group. All rights reserved, including rights for text and data mining and training of artificial technologies or similar technologies.