Precise optical modeling for silicate-based white LEDs - OSA Publishing

Precise optical modeling for silicate-based white LEDs Ching-Cherng Sun, Ching-Yi Chen, Hsin-Yin He, Cheng-Chien Chen, Wei-Ting Chien, Tsung-Xian Lee and Tsung-Hsin Yang Institute of Lighting and Display Science, National Central University, Chung-Li 320, TAIWAN Department of Optics and Photonics, National Central University, Chung-Li 320, TAIWAN [email protected]

Abstract: In this paper, as to our best knowledge, we propose and demonstrate the first precise phosphor modeling scheme to simulate the chromatic performance of white LEDs with silicate phosphors. The phosphor model is useful to accurately simulate the power ratio of the blue and yellow lights emitted by the white LEDs and is important in white LED package. ©2008 Optical Society of America OCIS codes: (230.3670) light-emitting diodes; (220.4298) Nonimaging Optics.

__________________________________________________________________________ References and links 1. 2. 3. 4. 5. 6. 7. 8. 9. 10. 11. 12. 13. 14. 15. 16. 17. 18. 19. 20.

A. Zukauskas, M. S. Shur, and R. Caska, Introduction to Solid-state Lighting (John Wiley & Sons, New York, 2002). D. A. Steigerwald, J. C. Bhat, D. Collins, R. M. Fletcher, M. O. Holcomb, M. J. Ludowise, P. S. Martin, and S. L. Rudaz, "Illumination with Solid State Lighting Technology," IEEE J. Selected Topics in Quantum Electron 8, 310-320 (2002). Optoelectronics Industry Development Association (OIDA), Light emitting diodes (LEDs) for general illumination: An OIDA technology roadmap update 2002. Optoelectronics Industry Development Assn., Washington DC (2002). S. J. Duclos, J. Jansma, J. C. Bortscheller, and R. J. Wojnarowski, “Phosphor Coating with Self-adjusting Distance from LED Chip,” United States Patent,US 6635363 B1 (2003).. M. Kerker, H. Chew, P. J. McNulty, J. P. Kratohvil, D. D. Cooke, M. Sculley, and M. P. Lee, “Light scattering and fluorescence by small particles having internal structure,” J. Histochem. Cytochem. 27, 250-263 (1979). C. C. Chang, R. Chern, C. C. Chang, C. Chu, J. Y. Chi, J. Su, I-Min Chan, and J. T. Wang, “Monte Carlo Simulation of Optical Properties of Phosphor-Screened Ultraviolet Light in a White Light-Emitting Device,” Jpn. J. Appl. Phys. 44, 6056-6061 (2005). S. J. Lee, “Analysis of light-emitting diodes by Monte-Carlo photon simulation,” Appl. Opt. 40, 14271437 (2001). M. S. Kaminski, K. J. Garcia, M. A. Stevenson, M. Frate, and R. J. Koshel, “Advanced Topics in Source Modeling,” Proc. SPIE 4775, 46 (2002). Z.-Y. Ting and C. McGill, “Monte Carlo simulation of light-emitting diode light-extraction characteristics,” Opt. Eng. 34, 3545-3553 (1995). Á. Borbély and S. G. Johnson, “Performance of phosphor-coated light-emitting diode optics in ray-trace simulations,” Opt. Eng. 44, 111308 (2005). A. Doicu and T. Wriedt, “Equivalent refractive index of a sphere with multiple spherical inclusions,” Appl. Opt. 3, 204-209 (2001). D. Toublanc, “Henyey-Greenstein and Mie phase functions in Monte Carlo radiative transfer computations,” Appl. Opt. 35, 3270-3274 (1996). C. F. Boren and D. R. Huffmarn, Absorption and scattering of Light by Small Particles (Wiley, 1983). A. A. Stelur, A. M. Srivastava, H. A. Comanzo, and D. D. Doxsee, “Phosphor Blends for Generating White Light from Near-UV/Blue Light-Emitting Devices,” United States Patent, US 6685852 B2 (2004). Intematix SY450-B, http://www.intematix.com/images/Catalog-2008_V2.pdf. I. Moreno, “Spatial distribution of LED radiation,” Proc. SPIE 6342, 634216 (2006). C. C. Sun, T. -X. Lee, S. -H. Ma, Y. -L. Lee, and S. -M. Huang, “Precise optical modeling for LED lighting verified by cross correlation in the midfield region,” Opt. Lett. 31, 2193-2195 (2006). W. T. Chien, C. C. Sun, and I. Moreno, "Precise optical model of multi-chip white LEDs," Opt. Express 15, 7572-7577 (2007). D. Kang, E. Wu, and D. Wang, “Modeling white light-emitting diodes with phosphor layers,” Appl. Phys. Lett. 89, 231102 (2006). N. Narendran, Y. Gu, J. P. Freyssinier-Nova, and Y. Zhu, “Extracting Phosphor-Scattered Photons to Improve White LED Efficiency,” Phys. Status Solidi A, 202, R60–R62 (2005).

#102345 - $15.00 USD

(C) 2008 OSA

Received 3 Oct 2008; revised 13 Nov 2008; accepted 18 Nov 2008; published 20 Nov 2008

24 November 2008 / Vol. 16, No. 24 / OPTICS EXPRESS 20060

21. K. Yamada, Y. Imai, and K. Ishii, “Optical Simulation of Light Source Devices Composed of Blue LEDs and YAG Phosphor,” J. Light Vis. Env. 27, 70-74 (2003). ______________________________________________________________________________________________

1. Introduction LEDs have been regarded as the most important and promising solid-state light source for the next generation lighting due to the advantages in efficiency, life time, chromatic performance, reliability and environmental protection [1-2]. Recently, white LEDs which are made by GaN covered with yellow phosphors have performed high energy efficiency caused by the great improvement in GaN epitaxy and package. The available commercialized high-power white LEDs achieve 100 lm/W in luminous efficiency [3], which is larger than most light sources in general lighting. However, the chromatic performance of LEDs is still not well controlled due to the variation in GaN chips and package structure [4]. In the package side, the geometry as well as phosphors are the both main issues impacting on the optical and chromatic performance. If we focus on the phosphors, the key parameters include the particle size, concentration, absorption coefficient [5-6], quantum efficiency and thickness. So far, there are several phosphors proposed to apply to GaN-based LED. The most attentions paid in the past years include YAG and silicate. Since different phosphors may perform different characteristics, the package parameters may completely different. Therefore, to the practical application, there is always a demand to white LEDs that should provide a precise and stable power and chromatic performance. Therefore, a precise optical model to simulate the optical effects caused by the phosphors and to describe the power efficiency and color presentation is necessary in LED package. In this paper, as to our best knowledge, this is the first time to propose and demonstrate a precise phosphor model which accurately predicts the optical and chromatic performance of silicate phosphors applied to GaN-based white LEDs. In the following, we will start from the measurement of optical parameters of the phosphors, followed by the simulation with MonteCarlo ray tracing [7-10] incorporated with Mie scattering [11-13]. Then we make the corresponding white LEDs with the parameters that we have figured out. Finally, a comparison between the experimental measurement and the simulation from the design will be made. 2. Optical modeling

：

The phosphor in the study is a silicate-based chemistry of (SrBaMg)2SiO4 Eu2+, which absorbs blue light and emits the lights covering the spectrum in visible range but with greenyellow dominant and thus is useful for white light generation with blue GaN dies [14]. The optical modeling process in this study is illustrated in Fig. 1, where we did the measurement of the scattering lights, absorption coefficient, and conversion efficiency of the phosphor plate in different concentrations and thickness. Then we build up an optical simulation model with Mie scattering in a simulation program with Monte Carlo ray tracing. The model is aimed to simulate the outside power ratio of the blue light by the GaN die and the yellow light by the phosphor pumped by a blue light source. Finally, the successful model is applied to predict the color behavior of LEDs. In the following, we will describe the details. The characterization of the phosphor starts from the measurement of the particle size, which is a key parameter in determining the fraction obscuration and the cross-sectional area of the phosphor mixed with silicone. The measured particle size is in a Gaussian distribution from 2 um to 30 um with the peak at around 17 um. Then we calculated the particle number and the concentration of the phosphor with the density of Dp=4.83g/cm3 [15]. Thus the parameters for Mie scattering are determined.

#102345 - $15.00 USD

(C) 2008 OSA



Scattering modeling

phosphor plate phosphor plate

Absorption coefficient phosphor plate

Conversion Efficiency Bad

Successful phosphor model

Verification with real LEDs

Simulation Yellow

Blue

Good

1.0 0.8

Exp Sim

10%

0.6 0.4 0.2 0.0 -200 -150 -100 -50

0

50 100 150 200

Radiation angle (deg.)

Normalized intensity (a.u.)

Normalized intensity (a.u.)

Fig. 1. The diagram of the process for phosphor modeling.

1.2 1.0

Exp Sim

25%

0.8 0.6 0.4 0.2 0.0 -200 -150 -100 -50

0

50 100 150 200

Radiation angle (deg.)

Fig. 2. The measured scattering distribution and the corresponding simulation of the phosphor plates for the phosphor concentrations of 10% and 15%, respectively.

The second step is to evaluate the validation of the scattering model in the simulation. Accordingly, measurement of the scattering light distribution of phosphor plates in various concentrations and thickness should be done. In the experiment, we made various phosphor plates with different concentrations (from 10% to 30%) and thickness (from 0.2 mm to 1.5 mm). In considering wavelength-independent property of the scattering and no photon conversion in the phosphor, we used a He-Ne laser as the light source, and the in-plane scattered light was measured. Figure 2 shows two of the simulated results for the concentration of 10% and 25% based on the measurement. Then the measured scattering characteristic is applied to fine-tune the phosphor parameters so the scattering model can describe the scattered light distribution for the same phosphor but at different concentration and thickness. The result as in Fig. 2 shows that the scattering model can be used to predict various phosphor plates with the same phosphor but with different concentration and thickness. The third step is for the measurement of the absorption coefficient and conversion efficiency of the phosphors. In the experiment, we focused on the concentration of 10% and 15% respectively with the thickness from 0.2 mm to 1.5 mm. The measurements were done with phosphor plates attached on a GaN die attached on a metal slug (here called bare LED), as shown in Fig. 3. In order to effectively figure out these two parameters, we had to build up a precise optical model of the bare LED. The optical model of the bare LED was built with the optical modeling process described in Ref. 16-18. In order to make the LED be easily modeled, we darkened some reflection surfaces to reduce the number of the parameters in the LED. Then we obtained the power of the blue and yellow lights with the measurement by an integrating sphere, as shown in Fig. 4 [19]. The blue light power with and without phosphor plate were used to determine the absorption coefficient of the phosphor plate in the phosphor model with a fitting process. The fitting absorption coefficient for 10% and 15% #102345 - $15.00 USD

(C) 2008 OSA



concentration obtained in the experiment was obtained 0.51 and 0.88 mm-1, respectively. With the measurement result of the yellow light, we further determined the conversion coefficient, which was 0.58 for among our silicate phosphor samples. Thus all the optical parameters required in the phosphor model were determined.

phosphor plate

Fig. 3. The phosphor plate covered a GaN die attached on a metal slug and its computer model in the simulation.

0.0024 10% 15%

0.0022

Output power (W)

Output power (W)

0.005 0.004 0.003 0.002 0.001

0.0020 0.0018 0.0016 0.0014 0.0012 0.0010

10% 15%

0.000 0.0

0.2

0.4

0.6

0.8

1.0

1.2

1.4

0.0008 0.0

1.6

0.2

0.4

0.6

0.8

1.0

1.2

1.4

1.6

Thickness (mm)

Thickness (mm)

(a)

(b)

0.7

Blue-T Blue-R Yellow-T Yellow-R

0.6 0.5 0.4 0.3 0.2 0.1

0.0 0.0 0.2 0.4 0.6 0.8 1.0 1.2 1.4 1.6

Thickness (mm)

(a)

Relative intensity (a.u.)


Fig. 4. The measured blue (a) and yellow (b) lights with the structures shown in Fig. 3.

0.7

Blue-T Blue-R Yellow-T Yellow-R

0.6 0.5 0.4 0.3 0.2 0.1

0.0 0.0 0.2 0.4 0.6 0.8 1.0 1.2 1.4 1.6

Thickness (mm)

(b)

Fig. 5. The forward/backward blue/yellow lights of a phosphor plate of the concentration of 15% for (a) experimental measurement and (b) the corresponding simulation.

#102345 - $15.00 USD

(C) 2008 OSA



The final step in verifying the parameters is to compare forward and backward blue and yellow lights of the phosphor plates in the experiment [20-21] and in the simulation. In the simulation, the blue lights were absorbed by the phosphors and partially turned to yellow lights. The optical paths of the non-absorbed lights, including blue and yellow lights, were described by Mie scattering model. The forward lights were measured by locating the phosphor plates on the entrance face of an integrating sphere, where a laser at 448 nm illuminated the phosphor plates to obtain enough blue and yellow lights in the measurement. In contrast, the phosphor plates were located in center of the integrating sphere to measure the backward lights, where an absorption box attached on the back of the phosphor plates to remove all the forward lights. Figure 5 shows an example of the experimental result and the corresponding simulation with the phosphor model, where the concentration is 15% and the thickness was from 0.2 mm to 1.5 mm. Though there are some differences between the simulation and experimental measurement, we think that the model is good enough to describe the power distribution of the blue and yellow lights. 3. Verification with white LEDs The phosphor model is aimed to precisely simulate the optical and color behavior of a practical LED with silicate phosphor. Therefore, we made a real package with GaN die covered with silicate phosphor, where the concentration was 10% and 15% respectively and the thickness of the phosphor was well controlled. In contrast to the model of the bare LEDs, here we had to put all the parameters of the package and phosphor in the optical model. In the simulation, the blue lights were emitted from the active layer and partially escaped from the die and went into the phosphor mixed with silicone. With the phosphor model based on Monte Carlo ray tracing and Mie scattering, we calculated the power of the blue and yellow lights. Fig. 6 shows both the LEDs with and without an encapsulating lens and the corresponding models in computer simulation. The emitting spectrum of the LEDs in the real samples and the corresponding simulation is shown in Fig. 7. Figure 8 and Table 1 show the comparison of the color coordinates of both cases in the experimental measurement and the simulation. We can find that the color differences (ΔxΔy) between the simulation and experimental measurement are almost within 0.004. Such small differences show that the developed phosphor model is accurate and successful.

(a)

(b)

Fig. 6. The picture and optical model of the white LEDs with (a) without (b) with encapsulation with a glass lens.

#102345 - $15.00 USD

(C) 2008 OSA



5. Summary In this paper, we present a phosphor modeling process to precisely simulate the power and chromatic performance with silicate phosphor covered on blue GaN dies. The process starts from the measurement of the particle size of the phosphor to determine some key parameters in Mie scattering model. Then scattering model of the phosphor plates in different concentrations and thickness was evaluated with the corresponding experimental measurement. The absorption coefficients and conversion efficient of the phosphor sample were determined with the measurement of the blue and yellow lights in a special experiment. The determined parameters incorporated with Monte Carlo ray tracing and Mie scattering model form the final phosphor model. In order to evaluate the developed model, we made comparisons between the simulated and measured the power ratio of the yellow and blue lights of the phosphor plates pumped by a blue laser. Finally, we applied the phosphor model to simulate the power and chromatic performance of packaged white LEDs with silicate phosphor covered on blue GaN LEDs. We can find that the color differences between the simulation and experimental measurement are almost within 0.004. The simulation by the phosphor model well predicts the color coordinates of the LEDs with and without lens encapsulation. Such a model will be helpful in determining the phosphor parameters in LED package and can also be useful to design a white LED in high optical and chromatic performance. w/o Lens

0.6 0.4 0.2 0.0 300

400

500

600

700

0.8 0.6 0.4 0.2 0.0 300

800

(b)

400

Wavelength (nm)

w/o Lens

800

0.6 0.4 0.2 0.0 600

700

W avelength (nm )

800

Experim ent Sim ualtion

1.0



0.8

(c)

700

with Lens Experim ent Sim ulation

500

600

(b)

1.0

400

500

Wavelength (nm)

(a)

300

Experiment Simulation

1.0



(a)

0.8

with Lens Experiment Simulation

1.0

0.8 0.6 0.4 0.2 0.0 300

400

500

600

700

800

W avelength (nm)

(d)

Fig. 7. The measured and simulated spectrum of the white LEDs with the phosphor concentrations of (a) 10% without encapsulation, (b) 10% with encapsulation with a glass lens,, (c) 15% without encapsulation, (d) 15% with encapsulation with a glass lens.

#102345 - $15.00 USD

(C) 2008 OSA



0.30 0.28 0.26

10%-Exp 10%-Sim 15%-Exp 15%-Sim

w/o Lens

0.24

with Lens

0.22 0.20 0.18

w/o Lens with Lens

0.16 0.20 0.22 0.24 0.26 0.28 0.30 0.32 0.34

Fig. 8. The measured and simulated chromatic performance of the white LEDs for the concentrations of 10% and 15%. Table 1. The comparison of the chromaticity coordinates for the LEDs with and without encapsulation with lens between simulation and experimental measurement. Concentration

10 %

15 %

Package Type

Simulation

x=0.241 y=0.200

x=0.230 y=0.182

x=0.281 y=0.283

x=0.248 y=0.223

Measurement

x=0.240 y=0.202

x=0.223 y=0.184

x=0.284 y=0.279

x=0.244 y=0.223

Deviation

Δ x= +0.001 Δ y= -0.002

Δ x= +0.007 Δ y= -0.002

Δ x= -0.003 Δ y= +0.004

Δ x= +0.004 Δ y= 0.000

Acknowledgment This study was sponsored by the Ministry of Economic Affairs of the Republic of China with the contract no. 94-EC-17-A-07-S1-043 and the National Science Council with the contract no. 97-2221-E-008-025-MY3.

#102345 - $15.00 USD

(C) 2008 OSA

