Combined optical simulation of the ray-tracing technique and the finite difference time domain method was used to investigate the effect of the emitter orientation and the photonic crystal layer on the outcoupling efficiency (OCE) of bottom-emission type organic light emitting diodes (OLEDs). The OLED with a horizontal emitter exhibited an opposite interference effect to that of one with a vertical emitter, which suggested that the OCE would be very sensitive to the emitter orientation at a fixed emitter-cathode distance. The OLED with a horizontal emitter exhibited much larger OCE than that with a vertical emitter did, which was due to the substantial difference in the radiation pattern along with the different coupling with the surface plasmon excitation. The OCE with a horizontal emitter was increased by approximately 1.3 times by inserting a photonic crystal layer between the indium tin oxide layer and the glass substrate. The present study suggested that appropriate control of the emitter orientation and its combination to other outcoupling structures could be used to enhance the OCE of OLEDs substantially. Organic light-emitting diodes (OLEDs) are widely used as display devices and recently began to enter the market of general lighting [1 , 2] . Since the light is emitted from the organic emitting layer having a high refractive index, a substantial portion of the generated light cannot escape from the organic layers and the glass substrate due to total internal reflection (TIR) [3] . The trapped light in OLED can be categorized into the glass mode and the organic (or waveguide) mode depending on the location where it is guided by TIR. In addition, coupling of the light with the metal cathode and the resulting plasmon excitation or absorption is another loss mechanism, in particular, at a short cathodeemitter distance [4] . Therefore, enhancing light extraction efficiency has been one of the most important strategies for increasing wall-plug efficiency of OLEDs. There have been many theoretical and experimental efforts to improve the outcoupling efficiency (OCE) of OLEDs. For example, the glass mode can be extracted from the substrate by forming either a microlens array or a diffusing layer on the glass substrate [5 - 12] . The internal bulk or interface structures were modified to extract the light trapped in the organic layers by, for example, reflector modification, insertion of low-index grids in the organic layers, interface or structure modification of the organic and ITO layers, and patterning the ITO electrodes along with an insertion of a high-conductivity polymer layer [13 - 20] . Hybrid structures may be used to extract both the trapped light in the waveguide mode and the glass mode [21 , 22] . In addition, transparent electrodes can be used to realize metal-free OLEDs, where the surface plasmon and absorption losses can be reduced significantly [23] . Recent studies showed that it was possible to control the emitting molecular orientation, thus the emitter orientation, with respect to the normal of the OLED layers to some degree [24 - 29] . The emission process in OLEDs can be regarded as a dipole transition of the emissive molecules. Achieving horizontal orientation of emitting dipoles is a very favorable condition for obtaining higher OCE than vertical or isotropic orientations. The efficiency values of these oriented OLEDs were found to be much larger compared to those of conventional OLEDs. The purpose of this simulation study is to investigate the effect of the emitter orientation in the emitting layer on the OCE of bottom-emission OLED with various thicknesses of ETL. In particular, the combined effect of the emitter orientation and photonic crystals on the OCE has been investigated in detail by using the finite difference time domain (FDTD) method [30 - 33] and the ray-tracing technique. The detailed dependence of the OCE of OLEDs with a photonic crystal layer on the ETL thickness may be a new aspect of this study and can serve as a new, useful data set which can be considered in the optimization of the OLED layer structure. Figure 1 shows a schematic structure of the OLED studied by the FDTD method. Commercially-available software (FDTD Solutions, Lumerical Solutions, Vancouver, Canada) was used to carry out the electromagnetic wave simulation. The OLED in Fig. 1 is bottom-emitting type consisting of a metal cathode (aluminum), an electron-transporting layer (ETL), a hole-transporting layer (HTL), an indium tin oxide (ITO) layer and a glass substrate. The simulation area was 10×10 µm ² . The refractive indices of these layers ( n ) are shown in the same figure. The thickness of the cathode, HTL, and ITO were fixed to 100, 50, and 150 nm, respectively. A point-dipole radiation source at a wavelength of 550 nm was placed at the center of the interface between the HTL and ETL. The oscillating direction of the dipole was set to be either horizontal (a) or vertical (b) with respect to the interface. The far-field angular radiant intensity distribution, detected in the glass substrate over the ITO layer at a distance of 95 nm, was imported into the ray-tracing software program (ASAP, Breault Research Organization, Tucson, USA), which was treated as the emitting ray distribution of an imaginary flat light source. This imaginary light source, with an area of 0.5×0.5 mm ² , was located at the center of the glass substrate. The area and the thickness of the glass substrate is 2×2 mm ² and 0.7 mm, respectively. According to this procedure, we could obtain the OCE escaped from organic/ITO layers into the glass substrate (denoted as “ITO Mode”) and that escaped from the substrate into the air (denoted as “Air Mode”). In addition to the investigation of the basic OLED structure, a photonic crystal layer was formed on the ITO layer. The SiN cylinder ( n =1.9) was put on a square lattice with a period of 500 nm. This period was found to be the optimum value for the wavelength of 550 nm [33] . Between the SiN cylinders was the glass material ( n =1.48). It is well known that the interference effect between the emitter radiation with its reflected field from the cathode is critically dependent on the distance between the emitter and the cathode [2] . This two-beam interference makes the OCE oscillate with the thickness of ETL [30 , 32 , 33] . In addition, the surface plasmon excitation and absorption in the metal cathode is another pathway of the generated light in OLEDs into heat dissipation, in particular when the emitter is close to the cathode. The plasmon excitation loss can easily be checked if we replace the aluminum cathode by an ideal perfect electrical conductor (PEC), the reflectance of the latter being 100% without any loss. Figure 2 shows the comparison of the light intensity (proportional to the squared electric field) propagated from the vertically oscillating electric dipole. It is clearly seen that the propagation of the electromagnetic waves is localized along the cathode-ETL interface for the aluminum cathode indicating substantial excitation of the surface plasmon, while some portion of the electromagnetic waves penetrates into the ITO and the glass substrate when the ideal reflector is used as a cathode. The radiation profiles reported in Ref. 17 also revealed that the radiation emission from vertical emitters was mainly limited to large angles below the critical angle resulting in a low extraction efficiency [17] . Figure 3(a) and (b) show the polar plots of the far-field intensity distributions proportional to | E | ² , where E is the electric field, detected in the glass substrate as a function of ETL thickness for both horizontal and vertical emitters. Figure 3(a) shows that the intensity along the normal direction increases and then decreases with increasing ETL thickness. On the other hand, the far-field intensity distribution of the vertical emitter shows the opposite effect, i.e., the intensity along the normal direction decreases with the increase in the ETL thickness from 20 to 120 nm. Figures 4(a) and (b) show the dependence of OCE escaped from organic/ITO layers into the glass substrate (ITO Mode) and that escaped from the substrate into the air (Air Mode) on ETL thickness for the horizontal and the vertical emitter, respectively. It becomes clear from the comparison between Fig. 4(a) and (b) that the two orthogonal emitters exhibit opposite interference conditions. For example, the first maximum of OCE (Air Mode) is located at ~ 80 nm for the horizontal emitter, which is close to the minimum of OCE for the case of the vertical emitter. This different interference effect can be explained by considering the phase relationship between the forward (towards the substrate) and backward (towards the cathode) electric fields in addition to the phase shift at the cathode. In the case of vertical emitter, the forward and the backward electric fields are initially out of phase with a phase difference of π , while they are in phase for the horizontal emitters [25] . These results clearly demonstrate that the emitter orientation is the most important factor for the OCE at a fixed emitter-cathode distance. The OCE in the “ITO Mode” and the “Air Mode” for the case of the horizontal emitter orientation are 87% and 43%, respectively, which are much larger than those for the case of the vertical emitter orientation. Figures 5(a) and (b) show the polar plots of the far-field intensity distributions detected in the glass substrate as a function of the ETL thickness for both horizontal and vertical emitters when the photonic crystal layer was inserted between the ITO layer and the glass substrate. (The dimensions of the SiN cylinders were optimized values as described below.) The overall change in the intensity as a function of ETL thickness is similar to the results shown in Fig. 3 . However, the intensity distribution exhibits a unique pattern having a four-fold symmetry, which corresponds to the symmetry of the photonic crystal placed on a square lattice. This directional property was also observed from the intensity distribution in the “Air Mode”. The selective enhancement of the light along particular directions is one of the disadvantages of the photonic crystal structure applied to OLEDs. Addition of scattering particles in the glass substrate may reduce this directional enhancement to some degree, making the intensity distribution closer to the Lambertian distribution and thus more suitable for illumination applications [34] . The OCE of OLED was investigated as functions of the radius and of the height of the SiN cylinders in the photonic crystal. (The period was fixed to 500 nm according to previous study [33] .) The height of the cylinder did not have substantial effect on OCE, while the radius has some effect on it. Figures 6(a) and (b) show the dependence of OCE in the “ITO Mode” on the radius of the SiN cylinder for the horizontal and the vertical emitters. According to this simulation result, optimized dimensions of the SiN cylinders could be obtained for both emitter orientations. The radius and the height of the cylinder were fixed to 130 nm and 100 nm for the horizontal orientation and 200 nm and 120 nm for the vertical orientation. Figures 7(a) and (b) display the dependence of OCE on ETL thickness for the horizontal and vertical emitters, respectively, obtained from the simulation for OLEDs having the optimized photonic crystal layer. The oscillating behavior of OCE reflects the two-beam interference effect, and the two cases exhibit almost an opposite behavior to each other. The OCE in the “Air Mode” is 54% at the ETL thickness of 100 nm in the case of horizontal orientation, while it is less than 13% over the whole ETL thickness in the case of vertical orientation. The application of the photonic crystal to the bottom-emission type OLED increases the OCE in the “Air Mode” from 43% to 54%, which indicates that the nanoscale photonic crystal layer is effective in extracting the light trapped in the organic/ITO layer. If the emitter orientation is not optimized, the OCE for the OLED with the optimized photonic crystal layer is reduced from 54% to 34.2% [33] . This study suggests that the control of the emitter orientation combined with other outcoupling structures such as the photonic crystal may be an effective way to increase the OCE of OLEDs. The effect of the emitter orientation and the photonic crystal layer on the OCE of bottom-emission type OLED was investigated by combined simulation of the ray-tracing technique and the FDTD method. The OLED with a horizontal emitter exhibited much larger OCE than that with a vertical emitter did, which was due to the substantial difference in the radiation pattern along with the different coupling with the surface plasmon excitation. The application of the photonic crystal between the organic/ITO layer and the glass substrate further increased the OCE by approximately 1.3 times, however, the intensity distribution was characterized by four-fold symmetry due to the square-lattice structure of the photonic crystal. The present study suggests that OCE of OLED may substantially be enhanced by controlling the emitter orientation and using other outcoupling structures such as the photonic crystal.