Accelerated 2D optical simulation of LEDs allows quicker validation of architectural improvements (MAGAZINE)

May 13, 2013
Finite-difference time-domain simulation in 3D is necessary to verify the light extraction efficiency of LED architectural enhancements, but a faster 2D method can allow for quicker iterations for testing new ideas during product development, write Chenglin Xu and Dan Herrmann.

This article was published in the April 2013 issue of LEDs Magazine.

View the Table of Contents and download the PDF file of the complete April 2013 issue, or view the E-zine version in your browser.


Fig. 1.LEDs, which began appearing in very limited applications more than 50 years ago, are now widely used for numerous applications, predominantly in the solid-state lighting (SSL) and display markets. Despite their widespread popularity, LEDs still exhibit some performance issues, the most notable of which is low extraction efficiency. Scientists and engineers are trying many approaches to increase the extraction efficiency with chip-level architectural enhancements, although traditionally it has been a slow iterative process. New 2D optical simulation techniques, however, can closely approximate real 3D analysis of point light sources and allow the engineers to test new ideas much more quickly with far less computer power and resources.

The limited extraction efficiency of LEDs is due to the fact that the light is generated inside media with a high refractive index and that total internal reflection occurs at the air interface. As a result, only a small portion of the light within the cone formed by the critical angle can be extracted into air, and the majority is trapped inside. The theoretical limit of the extraction efficiency is 1⁄4n2 (where n is the refractive index of the medium), which is only about 2% for GaAs-LEDs and 4% for GaN-LEDs.

To break through the theoretical limit, researchers are trying to manipulate the LED structure to extract more light. To accomplish this, many researchers rely on a trial-and-error approach to the design process. For example, S. David Roh from LG Innotek reported iterating as many as several hundred times in a paper at the Asia Communication and Photonics Conference last year. Although most designers may acknowledge that simulation software can reduce the design cycle and save time and money, not all designers use simulation tools in their practical design, perhaps because of the complexity of LED simulation problems or because they have not yet found an efficient, user-friendly simulation tool.

Fig. 2. It is well known that an LED is an incoherent and randomly polarized light source. It is temporally incoherent, spatially incoherent, as well as polarization incoherent. To deal with the complexity of the problem, we have to simulate the 3D structure by a finite-difference time-domain (FDTD) method at every position, every wavelength, and every polarization, and then sum the results of each simulation incoherently. To automate this complicated simulation process, RSoft (now part of the Synopsys Optical Solutions Group) developed the LED Utility that runs on top of the FullWave optical simulation tool.

Still, a full 3D FDTD simulation is very expensive in terms of computation time and computer memory. One simulation for a typical LED structure could use tens of gigabytes of RAM and take a few hours on a multi-core computer. For multiple scans over wavelength, position, and polarization, it is very common for a simulation to take a few days.

To increase efficiency and make the simulation more affordable, we recently proposed and briefly reported a simplified approach in a paper at the Asia Communication and Photonics Conference. Instead of simulating a full 3D structure, we assume that the LED device is circularly symmetric so that we can perform a 2D radial simulation. The simulation effort, including both computation time and computer memory, can be reduced by hundreds or even thousands of times.

General simulation procedure

Fig. 3. Let’s first consider the typical simulation procedures and then our simplified 2D methodology. Photons generated inside an LED by injected current are completely independent of each other, with random phase variations and polarizations. Hence, the LED light source is unpolarized as well as spatially and temporally incoherent. The cost of including incoherence via a Monte Carlo approach in an FDTD simulation, the most commonly used method for LED simulation, is very high. The time scale for a dipole’s phase variation is much larger than the optical period (T=2π/w) and would require prohibitively long simulations. A better approach is to analytically create the unpolarized and incoherent source by combining simulations of individual dipoles.

Unpolarized light is uniformly polarized in all directions. To create such a source, we can separately simulate the three orthogonal polarizations and sum the results. Shown in Fig. 1 are the far fields (also called radiation patterns) of a dipole polarized in three different directions, as well as the sum of all polarizations.

It is clearly shown that the radiation pattern of each dipole is polarization-dependent, and their sum is uniform in all directions and polarization-independent. Therefore, for each point source, three separate simulations for different orthogonal polarizations are needed. For certain types of LEDs, photons may not radiate evenly in all directions. For instance, in a GaN-based LED, there is no radiation in the crystal axis, so simulating two lateral polarizations is sufficient.

Fig. 4. The phase variation of each generated photon is independent of other photons and takes place on a time scale that is too long to model. Simulating multiple point sources together maintains a fixed phase difference between the point sources and results in a coherent simulation. Therefore, multiple point sources must be simulated separately to obtain spatially incoherent results.

Temporal coherence

The temporal coherence of light, which is measured by a coherent length, is inversely proportional to the spectral width. Due to the dominance of spontaneous emission, as shown in Fig. 2, an LED has a much wider emission spectrum than a laser. The extraction efficiency and radiation pattern of an LED are wavelength dependent, as shown for a simulated example in Fig. 2.

The final extraction efficiency R and far-field pattern ψ(θ,φ) can be obtained by averaging the wavelength-dependent results R(λi) and ψ(θ,φ,λi) at each sampled wavelength λi, weighted by the LED radiation spectrum S(λi), that is,

where N is the number of the sampled points. The combined results are shown in Fig. 3.

2D approximation

Fig. 5. As described in the previous section, a full 3D FDTD simulation of an LED structure involves scanning over wavelength, position, and polarization and can require a significant amount of computer memory. The computer resource issue is especially a problem for a patterned structure for which a much larger simulation domain is required. Patterned structures are being increasingly used in LEDs to improve light extraction. Fig. 4 compares flat and patterned structures. In a patterned structure, the trapped light could be extracted by the gratings while it propagates in the lateral directions. By comparison, a flat structure needs a much smaller computation domain, as long as it is large enough to cover the light cone formed by the critical angle.

For a typical 3D structure, a full simulation based on the preceding procedure could take a few days with tens of gigabytes of memory in order to obtain reasonable results. Such a simulation certainly limits the ability of design automation. Is there a more efficient way to simulate LEDs, even if it is an approximated method? The answer might be 2D approximation, which is appropriate for most photonic device simulations. For LEDs, however, 2D approximation oversimplifies the problem.

Fig. 6. As shown in Fig. 5, the 2D approximation changes the problem completely. It turns a 2D photonic crystal (PhC), a common, textured pattern for enhancing extraction efficiency, into 1D gratings. Most importantly, it turns a point source into a line source. As a result, the 2D simulation overestimates the extraction efficiency and projects a radiation pattern that is significantly different from the 3D results.

Radial approximation

We can improve the 2D results using a radial approximation. The commonly used 2D PhC patterns, either hexagonal or square lattices, are approximations to circular Bragg gratings, which create bandgaps in the lateral directions to force the trapped light to travel in a vertical direction. Therefore, 2D radial approximation is closer to the real problem geographically, as shown in Fig. 6.

To validate the 2D radial simulation, we first tested a flat structure, for which the 2D radial approximation is exactly the same as the original 3D structure, which is confirmed by the simulation results shown in Fig. 7. The 2D radial simulation produces almost the same results as the 3D simulation, and the discrepancy is purely numerical because of grid shapes. For the sake of comparison, the 2D results, with the line sources, are shown in the same chart, and they are very different from the 3D results.

Fig. 7. We also validated the 2D radial approximation for a patterned structure, and the results are shown in Fig. 8. As expected, there is some discrepancy between 2D radial and 3D results, and the discrepancy will be larger for a stronger patterned structure. Once again, the 2D and 3D simulation results are very different.

Still, a 2D radial simulation is a good approximation to real 3D LED problems, and it produces results that are similar to those of a 3D simulation. The computation effort, for both time and memory, is usually several orders of magnitude smaller. A 3D simulation for a typical LED structure could take a few days and use tens of gigabytes of memory, and a 2D radial simulation takes only minutes and uses very little memory.

Fig. 8. In conclusion, we developed a general simulation procedure for an LED structure that can handle an unpolarized and incoherent light source. To improve simulation efficiency, we demonstrated a 2D radial simulation of 3D LEDs that can reduce the complexity of the simulation task by hundreds or thousands of times. Using a tool such as FullWave and the LED Utility, engineers can quickly test ideas for new LED architectural structures that can improve extraction efficiency.

Although it is an approximation, the 2D method can still provide design insight and prove the design concept quickly. We do not expect that 3D simulation could be completely replaced by 2D radial simulation, as 3D simulation still serves as a verification tool for the final design. Also, the LED Utility can be used to perform 3D simulation as the LED designer closes on an architectural development project.