Optical Power Distribution in Graphene OLEDs

Lee et al. (2016): optical power distributions for TiO2, graphene, and low-index hole-injection layers at 550 nm
AuthorLuke Cole

Synergetic electrode architecture for efficient graphene-based flexible organic light-emitting diodes

Authors: Jaeho Lee, Tae-Hee Han, Min-Ho Park, Dae Yool Jung, Jeongmin Seo, Hong-Kyu Seo, Hyunsu Cho, Eunhye Kim, Jin Chung, Sung-Yool Choi, Taek-Soo Kim, Tae-Woo Lee, and Seunghyup Yoo

Journal: Nature Communications 7, 11791 (2016) · Comparison target: Supplementary Table 2 and Table 5 · License: CC BY 4.0

Adding 55 nm of TiO₂ under the graphene anode lifts the light that reaches air from 31.5% to 43.4%. In the same stack, swapping the low-index hole-injection layer GraHIL for PEDOT:PSS pushes evanescent loss from 22.4% to 28.6%.

This case builds four models at 550 nm and reproduces both changes with the Emission Mode result, channel by channel against the paper's two supplementary tables.

Lee Supplementary Table 2 — air, substrate, waveguide, evanescent, absorption, and other-loss fractions for three anodes at 550 nm.Lee et al., Nature Communications 7, 11791 (2016), Supplementary Table 2CC BY 4.0

Every row closes to 100%. TiO₂/graphene/GraHIL has the highest air fraction, graphene/GraHIL the largest waveguide fraction, and ITO/GraHIL the largest absorption.

Background

In a bottom-emitting OLED the light generated by the emitter splits between air, the glass substrate, waveguide modes in the organic layers, and the evanescent field near the metal cathode. The paper places high-index TiO₂ under the graphene to push more light toward air with cavity resonance, and inserts low-index GraHIL between graphene and the organic stack to weaken evanescent loss at the cathode.

Structure

The paper's optical model merges every organic layer except the hole-injection layer into one equivalent uniform film, 90 nm thick for the TiO₂ device. This case follows that convention.

The primary stack runs from the glass outcoupling side to the Al cathode.

LayerThicknessOptical constants at 550 nmSource
Glass, incoherent1 mm1.52Thick substrate
TiO₂55 nm2.50Optical constants published in the paper
Four-layer graphene1.36 nm2.60 + 1.30iGraphene optical model cited by the paper; four sheets of 0.34 nm
GraHIL70 nm1.42 + 0.0013iOptical constants published in the paper
Equivalent organic emitting layer90 nm1.80The paper's uniform-layer convention
LiF1 nm1.39Transparent-layer approximation at one wavelength
Al150 nm0.92 + 6.40iOptical constants published in the paper; optically thick cathode

Entered as above, the structure page should look like this. Check the 55/1.36/70/90/1/150 nm layer order, the incoherent flag on Glass, and the Emis. toggle on the equivalent organic layer.

The other three models change only the anode or the hole-injection layer; every remaining layer stays as entered.

ModelChange from the primary stackQuestion it answers
Graphene / GraHILRemove the 55 nm TiO₂What TiO₂ contributes to air outcoupling
ITO / GraHILReplace TiO₂ and graphene with 185 nm ITO, 1.90 + 0.02iHow the graphene anode differs from a conventional ITO anode
TiO₂ / graphene / PEDOT:PSSReplace GraHIL with 50 nm PEDOT:PSS, 1.56 + 0.003iWhat a low-index hole-injection layer does to evanescent loss

The stack with TiO₂ removed:

TiO₂ and graphene together replaced by a single 185 nm ITO layer:

Only the 70 nm GraHIL replaced by 50 nm PEDOT:PSS, anode side untouched:

The paper does not publish the complex refractive index of its ITO at 550 nm. This case uses 1.90 + 0.02i, a common value for transparent conducting oxides, so absorption and cavity phase in the ITO row carry more input uncertainty than the other three.

Optical and Emitter Settings

All four models share the same emitter and detector settings; once the structure is built nothing else changes.

SettingValue
EmitterIr(ppy)₂acac
SpectrumUnit White (no spectral weighting in a single-wavelength run)
Dipole orientationCustom, vertical dipole fraction 0.23, i.e. about 77% horizontal
Depth distributionDelta, 70 nm from the organic/Al interface
Relative position in the layer0.2222 measured from the outcoupling side
Intrinsic radiative quantum efficiency0.97
Electron–hole balance1
DetectorMode, Single, 550 nm

Filled in, the emitter panel looks like this. The 0.2222 relative position places the dipole plane 70 nm from the Al:

Enable Mode only, and set the wavelength mode to Single:

A single-wavelength run reads no broadband emission spectrum, so the Ir(ppy)₂acac EL curve does not need to be digitized. Run each of the four models once.

If the only parameter notices at the bottom of the page concern Optimize, go ahead and Run; they do not affect a fixed-parameter Mode calculation.

Simulation Results and Comparison with Supplementary Tables 2 and 5

The Mode result reports seven channels, which map onto the paper's six as follows.

Paper channelDreapex TMM channelNote
Outcoupled / AirTop OutcouplingSame quantity
SubstrateSubstrateSame quantity
W/G (waveguide)WaveguideSame quantity
EvanescentEvanescentSame quantity
AbsorptionAbsorptionSame quantity
Other lossesNon-radiativeThe paper merges non-radiative decay with residual terms; the software lists it separately
Not listed separatelyBottom OutcouplingNear zero behind the thick Al cathode

Anode structure: Supplementary Table 2

The published table:

Lee Supplementary Table 2 — air, substrate, waveguide, evanescent, absorption, and other-loss fractions for three anodes at 550 nm.Lee et al., Nature Communications 7, 11791 (2016), Supplementary Table 2CC BY 4.0

The paper's values (grey bars) and this run (blue bars), drawn in the same channel order with shared axes across all four panels:

Direct comparison of the paper and Dreapex TMM optical power distributions at 550 nm; all four panels use the same axes and channel order.Lee et al. 2016 Supplementary Tables 2 and 5 and the completed simulations in this case; comparison graphic independently drawnIndependent comparison graphic: CC BY 4.0

The bars track channel by channel. The key numbers for the three anodes:

ModelAir: paper / DreapexWaveguide + evanescent: paper / Dreapex
TiO₂ / graphene / GraHIL43.4% / 45.96%29.10% / 26.74%
Graphene / GraHIL31.5% / 28.08%30.20% / 29.18%
ITO / GraHIL28.9% / 28.53%32.40% / 31.74%

Adding TiO₂ gives the highest air fraction, removing it gives the highest waveguide fraction, and switching to ITO gives the largest absorption — the same ordering as the paper.

The three raw Mode results from the app follow, in the order of the table above. TiO₂/graphene/GraHIL:

With TiO₂ removed, the Top Outcoupling band narrows visibly:

With the ITO anode, absorption is the largest of the three:

Hole-injection layer: Supplementary Table 5

A separate table compares the two hole-injection layers:

Lee Supplementary Table 5 — effect of GraHIL and PEDOT:PSS on the optical power distribution at 550 nm.Lee et al., Nature Communications 7, 11791 (2016), Supplementary Table 5CC BY 4.0

The Mode result after switching to PEDOT:PSS, with a visibly wider evanescent band:

Hole-injection layerAir: paper / DreapexEvanescent: paper / Dreapex
GraHIL, 70 nm43.4% / 45.96%22.4% / 11.99%
PEDOT:PSS, 50 nm40.2% / 38.52%28.6% / 29.23%

Switching to PEDOT:PSS raises evanescent loss and lowers air outcoupling, the same direction as the paper. The PEDOT:PSS row tracks channel by channel; the low evanescent value in the GraHIL row is explained below.

All four runs close to 100.0000% across seven channels, with Bottom Outcoupling below 0.000001% and Non-radiative between 1.93% and 2.24%, the same order as the paper's 2.1–2.3% other losses.

Deviation Notes

The paper publishes the optical constants of TiO₂, GraHIL, PEDOT:PSS, and Al only as curves, so reading single-wavelength values introduces error in cavity phase and absorption split. The ITO row additionally uses a proxy refractive index.

The clearest deviation is the boundary between the waveguide and evanescent channels. Both hold high in-plane-wavevector power; the software divides them by mode character and the paper by a different boundary. Summed, the two agree to within 0.66–2.36 percentage points; split apart, the individual gap grows to about 10 points.

The paper merges non-radiative decay and assorted residuals into one "other losses" entry, while the software lists Non-radiative and Bottom Outcoupling separately, so those two can only be compared by order of magnitude.

Back to Case Studies

Copyright © 2026 Dreapex