(57)【要約】    【課題】有機電界発光素子の発光特性を精度よくシミュレーションする。    【解決手段】有機電界発光素子の取り出し光の発光特性を求めるシミュレーション方法であって、発光層内に存在する発光分子が励起状態から無輻射で基底状態に遷移する単位時間あたりの確率をｋｔ、角度（θ，ψ）方向を含む微小立体角ｄΩ内へ波長λからλ＋ｄλのS偏光の光子を放出して基底状態に遷移する単位時間あたりの確率をｋｒｓ（θ，ψ，λ）ｄΩｄλ、角度（θ，ψ）方向を含む微小立体角ｄΩ内へ波長λからλ＋ｄλのP偏光の光子を放出して基底状態に遷移する単位時間あたりの確率をｋｒｐ（θ，ψ，λ）ｄΩｄλとして、励起状態にある発光分子から（θ，ψ）方向を含む微小立体角ｄΩ内へ波長λからλ＋ｄλの光子が放出される確率φ（θ，ψ，λ）ｄΩｄλを数式（１）より求める。【選択図】  なし【特許請求の範囲】
【請求項１】
有機電界発光素子の取り出し光の発光特性を求めるシミュレーション方法であって、発光層内に存在する発光分子が励起状態から無輻射で基底状態に遷移する単位時間あたりの確率をｋｔ、角度（θ，ψ）方向を含む微小立体角ｄΩ内へ波長λからλ＋ｄλのS偏光の光子を放出して基底状態に遷移する単位時間あたりの確率をｋｒｓ（θ，ψ，λ）ｄΩｄλ、角度（θ，ψ）方向を含む微小立体角ｄΩ内へ波長λからλ＋ｄλのP偏光の光子を放出して基底状態に遷移する単位時間あたりの確率をｋｒｐ（θ，ψ，λ）ｄΩｄλとして、励起状態にある発光分子から（θ，ψ）方向を含む微小立体角ｄΩ内へ波長λからλ＋ｄλの光子が放出される確率φ（θ，ψ，λ）ｄΩｄλを数式（１）より求めることを特徴とする有機電界発光素子のシミュレーション方法。
【数１】
2007095500.tif 000056
（λｍｉｎは発光分子から放出されることが可能な最短波長を、λｍａｘは発光分子から放出されることが可能な最長波長を表す。）
【請求項２】
発光波長オーダーの距離にミラーの存在しない、あるいはミラーによる反射が微小な位置に置かれた（θＴＭ，ψＴＭ）方向を向いた発光分子のｋｒｓ（θ，ψ，λ）、およびｋｒｐ（θ，ψ，λ）の値から光学薄膜理論を用いて有機電界発光素子の発光層内に置かれた（θＴＭ，ψＴＭ）方向を向いた発光分子のｋｒｓ（θ，ψ，λ）、およびｋｒｐ（θ，ψ，λ）の値を算出して用いることを特徴とする請求項１記載の有機電界発光素子のシミュレーション方法。
【請求項３】
有機電界発光素子の取り出し光の発光特性を求めるシミュレーション装置であって、発光層内に存在する発光分子が励起状態から無輻射で基底状態に遷移する単位時間あたりの確率をｋｔ、角度（θ，ψ）方向を含む微小立体角ｄΩ内へ波長λからλ＋ｄλのS偏光の光子を放出して基底状態に遷移する単位時間あたりの確率をｋｒｓ（θ，ψ，λ）ｄΩｄλ、角度（θ，ψ）方向を含む微小立体角ｄΩ内へ波長λからλ＋ｄλのP偏光の光子を放出して基底状態に遷移する単位時間あたりの確率をｋｒｐ（θ，ψ，λ）ｄΩｄλとして、励起状態にある発光分子から（θ，ψ）方向を含む微小立体角ｄΩ内へ波長λからλ＋ｄλの光子が放出される確率φ（θ，ψ，λ）ｄΩｄλを数式（１）より求める手段を含むことを特徴とする有機電界発光素子のシミュレーション装置。
【請求項４】
発光波長オーダーの距離にミラーの存在しない、あるいはミラーによる反射が微小な位置に置かれた（θＴＭ，ψＴＭ）方向を向いた発光分子のｋｒｓ（θ，ψ，λ）、およびｋｒｐ（θ，ψ，λ）の値から光学薄膜理論を用いて有機電界発光素子の発光層内に置かれた（θＴＭ，ψＴＭ）方向を向いた発光分子のｋｒｓ（θ，ψ，λ）、およびｋｒｐ（θ，ψ，λ）の値を算出して用いる手段を含むことを特徴とする請求項３記載の有機電界発光素子のシミュレーション装置。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、表示素子、バックライト、照明、インテリア、標識、看板、電子写真機、光信号発生器などの分野に利用可能な有機電界発光素子の発光特性を予測するシミュレーション方法およびその装置に関する。
【背景技術】
【０００２】
基板上に形成された第一電極（一般には陽極）から注入された正孔と、第二電極（一般には陰極）から注入された電子とが、両極に挟まれた有機蛍光体や有機リン光体内で再結合する際に発光する原理を利用する有機電界発光素子の研究が近年活発に行われている。この素子は、薄型化が可能であり、低駆動電圧下での高輝度発光が可能であり、蛍光材料を選ぶことにより多色発光が可能であるという特徴を有している。
【０００３】
有機電界発光素子が高輝度に発光することは、コダック社のＣ．Ｗ．Ｔａｎｇらによって初めて示された（例えば、非特許文献１参照）。コダック社の提示した有機電界発光素子の代表的な構成は、ＩＴＯガラス基板上に正孔輸送性のジアミン化合物、発光層であり、かつ電子輸送性も併せ持ったトリス（８−キノリノラト）アルミニウム、および陰極としてのＭｇ：Ａｇを順次設けたものであり、１０Ｖ程度の駆動電圧で１０００ｃｄ／ｍ２の緑色発光が可能であった。現在の有機電界発光装置は、基本的にはコダック社の上記構成を踏襲しており、基板上に第一電極と、発光層を含む薄膜層と、第二電極とが順次積層された構造を有している。そして、薄膜層は、発光層のみの単層構造であってもよいが、多くの場合には正孔輸送層や電子輸送層を設けた、複数層からなる積層構造である。
【０００４】
ところで、有機電界発光素子においては、素子内部で光が干渉するため、同一材料を用いて素子を作製しても、素子構造によって発光スペクトル、色度、輝度などの発光特性が変化することが知られている。そして、この光学干渉による発光特性の変化をシミュレーションする方法が提案されている（特許文献１参照）。
【０００５】
特許文献１の方法は、発光面の法線方向に対してある角度θで発光面から観察側に放射される直接光と、前記発光面から放射されて光反射層で一度以上反射された後、前記直接光と同じ角度θで観察側に向う反射光との振幅比もしくは位相差を、複素屈折率を用いてフレネルの式から導かれる有機電界発光素子を構成する各層界面での反射、透過時の振幅、位相の変化と、光吸収と、光路差とに加えて、発光面から放射された光が光反射層によって反射される際に正反射以外の角度に散乱されて生じる散乱光を考慮して求めるシミュレーション方法であり、発光分子から放射される光の振幅に異方性はなく等方的であるとし、光反射層によって反射される際に正反射以外の角度に散乱されて生じる散乱光の存在を考慮して任意の角度方向の発光特性を予測しようというものである。
【特許文献１】特開２００５−３８６５９号公報
【非特許文献１】アプライド・フィジックス・レターズ（Ａｐｐｌ．Ｐｈｙｓ．Ｌｅｔｔ．）５１（１２）２１，９１３（１９８７）
【発明の開示】
【発明が解決しようとする課題】
【０００６】
しかしながら、上記のように散乱光をも考慮する特許文献１のシミュレーション方法では、発光量子収率の大きな発光材料を用いた時、シミュレーションで得た基板法線方向の輝度の有機電界発光素子を構成する各層膜厚依存性と実際の発光素子で得られる基板法線方向の各層膜厚依存性との間には大きな違いがあり、シミュレーションとしての精度は必ずしも高いものではなかった。
【０００７】
本発明は、励起状態から基底状態へ光子を放出して遷移する単位時間あたりの確率を求め、さらに励起状態から基底状態へ無輻射で遷移する単位時間あたりの確率を考慮することで、実測値と合致する精度のよいシミュレーション結果を得ることができるシミュレーション方法およびその装置を提供することを目的としている。
【課題を解決するための手段】
【０００８】
すなわち本発明は、構成が既知の有機電界発光素子の取り出し光の発光特性を求めるに当たって、発光層内に存在する発光分子が励起状態から無輻射で基底状態に遷移する単位時間あたりの確率をｋｔ、角度（θ，ψ）方向を含む微小立体角ｄΩ内へ波長λからλ＋ｄλのS偏光の光子を放出して基底状態に遷移する単位時間あたりの確率をｋｒｓ（θ，ψ，λ）ｄΩｄλ、角度（θ，ψ）方向を含む微小立体角ｄΩ内へ波長λからλ＋ｄλのP偏光の光子を放出して基底状態に遷移する単位時間あたりの確率をｋｒｐ（θ，ψ，λ）ｄΩｄλとして、励起状態にある発光分子から（θ，ψ）方向を含む微小立体角ｄΩ内へ波長λからλ＋ｄλの光子が放出される確率φ（θ，ψ，λ）ｄΩｄλを数式（１）より求めることを特徴とする有機電界発光素子のシミュレーション方法である。
【０００９】
【数１】
2007095500.tif 000002
【００１０】
λｍｉｎは発光分子から放出されることが可能な最短波長を、λｍａｘは発光分子から放出されることが可能な最長波長を表す。
【００１１】
また、発光層内に存在する発光分子が励起状態から無輻射で基底状態に遷移する単位時間あたりの確率をｋｔ、角度（θ，ψ）方向を含む微小立体角ｄΩ内へ波長λからλ＋ｄλのS偏光の光子を放出して基底状態に遷移する単位時間あたりの確率をｋｒｓ（θ，ψ，λ）ｄΩｄλ、角度（θ，ψ）方向を含む微小立体角ｄΩ内へ波長λからλ＋ｄλのP偏光の光子を放出して基底状態に遷移する単位時間あたりの確率をｋｒｐ（θ，ψ，λ）ｄΩｄλとして、励起状態にある発光分子から（θ，ψ）方向を含む微小立体角ｄΩ内へ波長λからλ＋ｄλの光子が放出される確率φ（θ，ψ，λ）ｄΩｄλを数式（１）より求める手段を含むことを特徴とする有機電界発光素子のシミュレーション装置である。
【発明の効果】
【００１２】
本発明によれば、実際に得られる有機電界発光素子の発光特性と高精度に合致するシミュレーション結果を得ることができる。
【発明を実施するための最良の形態】
【００１３】
以下、添付図面を参照して、本発明の有機電界発光素子のシミュレーション方法およびその装置の実施の形態を詳細に説明する。
【００１４】
本発明でいう取り出し光の発光特性とは、ある任意の膜厚構成の素子から発光面の法線方向に対してある角度で外部に出射される光の発光スペクトル、色度、輝度の少なくとも１つを含む概念として使用される。換言すれば、取り出し光の発光特性は、外部に出射される光の発光スペクトル、色度、輝度の少なくとも１つの膜厚依存性である。
【００１５】
以下では、まず発光量子収率の小さな発光材料を用いて素子を作製した時に利用できる双極子モデルについて述べ、続いて本発明である発光量子収率の大きな発光材料を用いて素子を作製した時に利用できる輻射遷移の失活速度定数算出モデルについて述べる。
【００１６】
発光面（図１中の２を参照）の法線方向に対してある角度θで発光面から観察側（図１中の５を参照）に放射される直接光（図１中の４を参照）と、発光面から放射されて光反射層（図１中の１を参照）で一度以上反射された後、直接光と同じ角度θで観察側に向う反射光（図１中の３を参照）とは、一般に振幅と位相が異なるが、それは下記三要素の影響である。
【００１７】
本発明における第一の要素とは、複素屈折率を用いてフレネルの式から導かれる、有機電界発光素子を構成する各層界面での反射、透過時の振幅、位相の変化であり、以下の数式（２）から数式（６）により表される。
【００１８】
複素屈折率Ｎ０の媒質０（図２、図３中の６を参照）から複素屈折率Ｎ１の媒質１（図２、図３中の７を参照）へ角度θ０で光が入射して角度θ１で透過する場合（図２、図３参照、なお、矢印はベクトルを表し、丸印中に黒丸の記号は、紙面に対して垂直に紙面手前側へ向かうベクトルを表し、丸印中に×印の記号は、紙面に対して垂直に紙面向こう側へ向かうベクトルを表す。これらの表記規約は、以下、他の図においても同様に適用される。）（また、図２、図３中、８は入射光を、９は反射光を、１０は透過光を、それぞれ表す）、添え字ｓ、ｐにより光の電場の入射面に垂直な成分（Ｓ偏光）と平行な成分（Ｐ偏光）とをそれぞれ表すと、振幅反射率及び振幅透過率は下記数式（２）から数式（５）のようになる。なお、ｒｓはＳ偏光の光の振幅反射率、ｔｓはＳ偏光の光の振幅透過率、ｒｐはＰ偏光の光の振幅反射率、ｔｐはＰ偏光の光の振幅透過率である。また、図においてＥは電場、Ｈは磁場をそれぞれ表し、図２に示すようにＳ偏光における電磁場の正方向を定義し、図３に示すようにＰ偏光における電磁場の正方向を定義する。数式（６）はスネルの法則である。
【００１９】
【数２】
2007095500.tif 000003
【００２０】
複素屈折率はＮｊ＝ｎｊ−ｉ・ｋｊで表される。ここでｎｊは屈折率、ｋｊは消衰係数であり、ｉは虚数単位である。そして、屈折率と消衰係数をまとめて光学定数と呼ぶ。光反射層はもちろん、有機電界発光素子を構成する各層の消衰係数は一般には０ではないので、一般に複素屈折率Ｎｊは複素数である。そのため振幅反射率、振幅透過率も複素数となり、反射、透過時の位相のずれの情報を含む。入射光の振幅に対する反射光の振幅の比、透過光の振幅の比はそれぞれ振幅反射率、振幅透過率の絶対値である。また、振幅反射率、振幅透過率の偏角を求めれば反射、透過時の位相のずれの値を求めることができる。有機電界発光素子の各層界面での反射、透過時の振幅、位相の変化を求めるには有機電界発光素子を構成する各層の光学定数の値が必要であるが、これは例えば偏光解析法（エリプソメトリ）等により求めることができる。
【００２１】
第二の要素とは光吸収である。有機電界発光素子を構成する各層の消衰係数は一般には０ではないので、光波の進行に伴ない光が吸収される。波長λの光波が消衰係数ｋの媒質中を距離ｄだけ進行すると振幅はｅｘｐ（−２πｋｄ／λ）倍となる。また光反射層での反射時には一部の光は反射されずに吸収される。反射時のエネルギー反射率は振幅反射率の絶対値の二乗であり、エネルギー吸収率は１からエネルギー反射率を差し引けば求まる。これら光吸収の結果、光の振幅が減少する。
【００２２】
第三の要素とは光路差である。光路差の分だけ位相が異なることになる。例えば図４に示すように屈折率ｎの媒質中にある発光面（図４中の１２を参照）の法線方向に対して角度θで発光面から観察側（図４中の１５を参照）に放射される直接光（図４中の１４を参照）と、発光面から放射されて光反射層（図４中の１１を参照）で反射された後、角度θ
で観察側に向う反射光（図４中の１３を参照）とは光路差２ｎｄ・ｃｏｓθを有する（ｄは発光面から光反射層までの距離である）。したがって直接光と反射光とは位相差
４πｎｄ・ｃｏｓθ／λを有する（ここでλは光の波長である）。
【００２３】
本発明においては、上記三要素に加えて、第四の要素として、発光分子から放射される光の振幅の異方性を考慮する。これまでの検討によると、発光分子から放射される光の振幅は等方的ではなく、異方的であることが判明した。そこで本発明は、従って実際と類似の、あるいは合致したシミュレーション結果を得るには、発光分子から放射される光の振幅の異方性を考慮することとした。この考え方に基づく、下記に説明する各数式を用いるシミュレーション方法では、発光面の法線方向に対してある角度θで発光面から観察側に放射される直接光と、前記発光面から放射されて光反射層に向かう光（光反射層で反射された後は前記直接光と同じ角度θで観察側に向かう光）との振幅比は一般に１ではなくなり、また、出射角度に応じて放射される光の振幅の値が異なることになる。従来の散乱光を考慮したシミュレーション技術では、発光面の法線方向に対してある角度θで発光面から観察側に放射される直接光と、前記発光面から放射されて光反射層に向かう光との振幅比は常に１であり、また、出射角度に応じて放射される光の振幅の値も常に一定としていたが、これは発光分子から光が放射されるという現象を物理的にとらえるならば正しい結果をもたらすことはできなかった。
【００２４】
有機電界発光素子を構成する各層の膜厚、光学定数、発光層内の膜厚方向の発光強度分布、発光層から発せられる光スペクトル、発光分子数の方向分布を用い（これらは、設計データとして与えられているので、実際には、図示しないメモリに保持させておき、必要に応じてメモリから読み出してシミュレーションのための処理部に供給する）、上記四要素を考慮することで、素子内部の全界面での反射、透過を考慮した内部発光干渉モデルを構成することができ、取り出し光の発光特性を求めることができる。計算の際には、Ｓ偏光とＰ偏光はそれぞれ独立に計算することができ、最後に足し合わせればよい。発光分子から放射される光の電場振幅は一般に異方的であり、本発明では電気双極子遷移による発光の場合を説明する。この場合、発光分子の遷移双極子モーメントが発光面の法線方向に対して垂直である場合を除き、発光面（図５中の１８を参照）の法線方向に対してある角度θで発光面から観察側（図５中の２１を参照）に放射される直接光（図５中の２０を参照）と、前記発光面から放射されて光反射層（図５中の１７を参照）に向かう光との振幅比は１ではなく、遷移双極子モーメント（図５中の２２を参照）の方向に応じて変化を受けることになる（図５）。図５には発光分子からの放射パターンの概念図を示したが、ここではＳ偏光の光とＰ偏光の光を足し合わせた全体の光エネルギーの放射パターンとして示した。後述の数式（２０）から（２３）に示すように、本発明に記載の座標系での計算式においてはＳ偏光については発光面から観察側に放射される直接光と前記発光面から放射されて光反射層に向かう光との振幅比は常に１となるが、Ｐ偏光についてはこの振幅比が遷移双極子モーメントの方向および光の出射方向に応じて変化する。
【００２５】
発光スペクトルを求めるには、各波長の発光エネルギーを計算するが、その際、まず発光層の屈折率（一般に波長依存性を有する）と発光界面の湧き出し光の電場振幅の二乗との積が、発光分子の遷移双極子モーメントの向きに対して湧き出し光の向きがなす角度に応じて定まった大きさになるようにして各波長の発光エネルギーを求める。次にこの値に発光層から発せられるスペクトル値（着目する波長のエネルギー）を掛けて発光スペクトルを求める。発光層から発せられるスペクトルとしては、発光層と同一の材料からなる膜の光励起発光スペクトルを使用すればよい。なお、発光分子は一般に様々な方向を向いているので、以下のようにする。
【００２６】
まず、遷移双極子モーメントの２つの偏角（方向）をそれぞれいくつかに分割し、遷移双極子モーメントが各角度方向を向いている場合の発光スペクトルを前記方法でそれぞれ独立に求める。次に各発光スペクトルにその角度方向を向いている発光分子数を乗じた後、この値を角度方向に亘って積分して、全方向からの寄与を足し合わせた全体での発光スペクトルを求める。
【００２７】
また、発光強度分布は以下のように取り扱えばよい。まず発光層を何層かに分割し、各層界面のみが発光している場合の発光スペクトルを前記方法でそれぞれ独立に求める。次に各発光スペクトルにその界面での発光強度を乗じた後、この値を厚み方向に亘って積分して、全界面からの寄与を足し合わせた全体での発光スペクトルを求める。このスペクトルから色度・輝度を求めることができる。
【００２８】
以上のようにすれば双極子放射による内部発光干渉モデル（双極子モデル）を構成することができる。Ａｌｑ３を発光層に用いた素子を作製し、素子の様々な視野角方向の発光効率を測定すると、視野角変化に対する発光効率の変化の仕方は双極子モデルによる計算と実験でよく一致し、実際の発光素子とほぼ同じ数値を得ることができた。これはＡｌｑ３の発光量子収率が小さいためである。しかし、発光量子収率の高いクマリン５４５Ｔをゲスト材料としてＡｌｑ３をホスト材料として発光層に用い、下記構造の化合物（化１）を電子輸送材料に用い、電子輸送層の膜厚を様々に変化させた素子を作製し、素子の基板法線方向の発光効率を測定すると、膜厚変化に対する発光効率の変化の仕方は計算と実験とでは大きな差が見られ、散乱を考慮する従来技術（散乱モデル）と同様の結果であった。
【００２９】
【化１】
2007095500.tif 000004
【００３０】
双極子モデルにおいては、発光分子がどのような向きにあっても、発光分子から基板正面方向へ放出される光と電極側に放出される光の振幅は同じになる。（後述の数式（２０）から数式（２３）参照）このため、基板正面方向の発光特性の膜厚依存性に関する限り、双極子モデルを用いても前記第一〜第三の要素のみからなる内部発光干渉モデルと同じ結果になる。散乱モデルとの相違は散乱光の存在の有無であるが、散乱光の比率が非散乱光と比べて小さい場合（特許文献１の実施例に記載のＡｌｑ３素子の場合など）干渉により発光効率が低くなった部分での発光特性に違いが出るに留まるのである。しかし、双極子モデルは発光量子収率が小さな発光材料を用いて作製した素子においては、その基板法線方向に対して任意の視野角方向の発光特性を正確に予測することができ、有用である。
【００３１】
本発明においては、発光量子収率の大きな発光材料を用いて作製した素子の発光特性を予測するために、励起状態にある発光分子が基底状態へ遷移する単位時間当たりの確率（失活速度定数）の素子構成に応じた変化を考慮する。励起状態から基底状態への遷移には光子放出を伴うものと伴わないものとがある。本発明では前者の輻射遷移の失活速度定数が素子構成に応じて変化することを特に考慮する。後者の無輻射遷移の失活速度定数が変化することも考慮してもよいが、本発明においては、素子内部で電界発光する場合の発光分子の無輻射遷移の失活速度定数は、波長オーダーの距離にミラーの存在しない、あるいはミラーによる反射が微小な位置に置かれた励起状態の発光分子の無輻射遷移の失活速度定数と同じであるとして話をすすめる。
【００３２】
励起状態にある発光分子から光子が放出される時、どの方向に放出されるかを決定することは原理的にできないが、方向を定めれば、その方向へ光子が放出される確率は求めることができる。発光分子が励起状態から角度（θ，ψ）方向を含む微小立体角ｄΩ内へ波長λからλ＋ｄλのS偏光の光子を放出して基底状態に遷移する単位時間あたりの確率をｋｒｓ（θ，ψ，λ）ｄΩｄλ、角度（θ，ψ）方向を含む微小立体角ｄΩ内へ波長λからλ＋ｄλのP偏光の光子を放出して基底状態に遷移する単位時間あたりの確率をｋｒｐ（θ，ψ，λ）ｄΩｄλとする。ｋｒｓ（θ，ψ，λ）ｄΩｄλとｋｒｐ（θ，ψ，λ）ｄΩｄλを加えて、全方向全波長で積分すると、発光分子が励起状態から光子を放出して基底状態に遷移する単位時間あたりの確率ｋｒｔｏｔａｌ（輻射遷移の失活速度定数）が求まる。ここで（θ，ψ）は球面座標系における二つの偏角である。
【００３３】
【数３】
2007095500.tif 000005
【００３４】
無輻射遷移の失活速度定数をｋｔとすると全失活速度定数はｋｔ＋ｋｒｔｏｔａｌとなる。従って、励起状態にある発光分子から（θ，ψ）方向を含む微小立体角ｄΩ内へ波長λからλ＋ｄλの光子が放出される確率φ（θ，ψ，λ）ｄΩｄλは（θ，ψ）方向を含む微小立体角ｄΩ内へ波長λからλ＋ｄλの光子が放出される輻射遷移の失活速度定数（ｋｒｓ（θ，ψ，λ）ｄΩｄλ＋ｋｒｐ（θ，ψ，λ）ｄΩｄλ）を全輻射遷移の失活速度定数（ｋｔ＋ｋｒｔｏｔａｌ）で割ることにより求まり、前記数式（１）となる。
【００３５】
以下、まず双極子モデルについて詳細に述べる。その後、本発明の、励起状態から基底状態へ光子を放出して遷移する単位時間あたりの確率を求め、さらに励起状態から基底状態へ無輻射で遷移する単位時間あたりの確率を考慮するモデル（輻射遷移の失活速度定数算出モデル）について詳細に述べる。
【００３６】
以下双極子モデルの説明であるが、まず発光分子の遷移双極子モーメントがある特定の方向を向いていて厚さ方向に発光分布がなく、発光がある界面に集中している場合、すなわちある方向を向いた発光分子が一つだけある場合について述べ、次にこれを用いて様々な方向を向いた発光分子が多数存在し発光分布が厚さ方向に広がっている場合について拡張する。
【００３７】
基板上に、ｍ層の薄膜層と陰極とを有する有機電界発光素子において、図6、図7のように層番号をつける。ｊ＝１〜ｍについて、第ｊ層の膜厚をｄｊとおく。同様にｊ＝０〜ｍ＋１について、波長λにおける屈折率をｎｊ（λ）、消衰係数をｋｊ（λ）とおき、複素屈折率をＮｊ（λ）＝ｎｊ（λ）−ｉ・ｋｊ（λ）とおく。（ここでｉは虚数単位である。）素子外部の波長λにおける屈折率をｎ（λ）とおく。（通常は空気であり、減衰率（ｋ（λ））＝０と近似し、又、真空の屈折率と同じであると近似して、Ｎ（λ）＝ｎ（λ）＝１である。）。
【００３８】
まず、発光が第（ｈ−１）層／第ｈ層界面に集中している場合を考え、発光する面を第（ｈ−１）層／第ｈ層界面とする。ここでは図８、図９に示すように第（ｈ−２）層、第（ｈ−１）層、第h層、第（ｈ＋１）層の計４層を発光層とし、Ｎ（ｈ−１）（λ）およびＮｈ（λ）を実数にとる。すなわち、Ｎ（ｈ−２）（λ）＝Ｎ（ｈ＋１）（λ）を発光層本来の複素屈折率とし、Ｎ（ｈ−１）（λ）＝Ｎｈ（λ）＝Ｒｅ［Ｎ（ｈ−２）（λ）］＝Ｒｅ［Ｎ（ｈ＋１）（λ）］（ｋ（ｈ−１）（λ）＝ｋｈ（λ）＝０）とする。（ここでＲｅは実数部分を示す。）そのかわり、複素屈折率の値を変化させた（吸収を０にした）ことによる影響を少なくするために、第（ｈ−１）層および第ｈ層の膜厚を０とする。すなわち、ｄ（ｈ−１）（λ）＝ｄｈ（λ）＝０とする。発光層の消衰係数は一般に０ではないが、発光界面の消衰係数が値をもつ場合、θhを実数にとると、基板に出射される方向θ0や外部に出射される方向θが複素数となってしまい、所望の実数の角度方向を計算することができないからである。θ方向の取り出し光に寄与する波長λの光の各界面での電場をＥｓｊα（λ）、Ｅｐｊｊα（λ）で表す｛ここでｓ、ｐはそれぞれｓ偏光、ｐ偏光を表す（ｓ偏光とは電場がｙｚ面に垂直な光であり、ｐ偏光とは磁場がｙｚ面に垂直な光である）。なお、下付き添え字ｊは層番号を表す。また、上付き添え字αは０か１の値を持ち、０はｚ軸の負の方向に向う光を、１はｚ軸の正の方向に向う光を示す｝。θ方向の波長λの取り出し光に寄与する光の各層での方向θｊ（θ，λ）に対して、ｓｉｎθｊ（θ，λ）、ｃｏｓθｊ（θ，λ）を数式（８）、数式（９）から（１１）を用いて求める。
【００３９】
【数４】
2007095500.tif 000006
【００４０】
次に発光分子から発せられる光の振幅を求める。発光分子の遷移双極子モーメントが図１０に示すようにＭ方向（θＴＭ，ψＴＭ）を向いている場合を考える。θＴＭおよびψＴＭはＭの偏角である。（図中において太字はベクトルであることを示す。以下も同様である。）着目するのはｋｈ０方向、ｋｈ１方向に向かって発せられる光である。ここでｋｈ０方向に向かうＳ偏光の光の電場方向の単位ベクトルをｅｓ０、P偏光の光の電場方向の単位ベクトルをｅｐ０としている。また、ｋｈ１方向に向かうＳ偏光の光の電場方向の単位ベクトルをｅｓ１、Ｐ偏光の光の電場方向の単位ベクトルをｅｐ１としている。ｋｈ０、ｋｈ１は波数ベクトルである。それぞれのベクトルのｘｙｚ成分は数式（１２）から数式（１８）のようになる。
【００４１】
【数５】
2007095500.tif 000007
【００４２】
なお、遷移双極子モーメントＭとは量子力学で光放射を伴う遷移の遷移確率を計算する際に現れる量で、双極子モーメントの演算子の始状態と終状態との間の行列要素として定義されるものであり、数式（１９）で与えられる。
【００４３】
【数６】
2007095500.tif 000008
【００４４】
ここでｅｋ、ｘｋは、ｋ番目の粒子の電荷および位置ベクトルを表す。またψｆ、ψｉは、それぞれ系の終状態および始状態の波動関数を表す。
【００４５】
本発明は、第四の要素として発光分子から放射される光の振幅の異方性を考慮するが、この要素は以下の数式（数式（２０）から数式（２３））で表され、シミュレーションの一要素として含まれる。
【００４６】
励起分子の、ある角度方向への光子放出の失活速度は、放出される光子の電場方向の単位ベクトルをｅとした場合｜ｅ・Ｍ｜２に比例するので、発せられる光の振幅は｜ｅ・Ｍ｜に比例すると考えることができる。すなわち、ｍ０＝１とすると、ｋｈ０方向に向かうＳ偏光の光の電場振幅をｖｓ０、Ｐ偏光の光の電場振幅をｖｐ０、ｋｈ１方向に向かうＳ偏光の光の電場振幅をｖｓ１、Ｐ偏光の光の電場振幅をｖｐ１として数式（２０）から数式（２３）のようになる。
【００４７】
【数７】
2007095500.tif 000009
【００４８】
次に実際のシミュレーションについてその手順を含め、以下に説明する。下記（１）〜（４）の場合における方程式（数式（２４）から数式（３９））を満たすａｓ（θ，λ）、ａｐ（θ，λ）を求め、その時のＥｓ００（θ，λ）、Ｅｐ００（θ，λ）を求める。これらは多層膜内から発光して基板に到達した光の電場である。
【００４９】
【数８】
2007095500.tif 000010
【００５０】
【数９】
2007095500.tif 000011
【００５１】
次に数式（４０）から数式（４３）を用いて角度θ０（θ，λ）で基板内をｚ軸の負の方向に向う波長λのｓ偏光の光の、基板から外部への透過率Ｔｓ（λ）、反射率Ｒｓ（λ）、ｐ偏光の光の基板から外部への透過率Ｔｐ（λ）、反射率Ｒｐ（λ）を求める。
【００５２】
【数１０】
2007095500.tif 000012
【００５３】
次に外部／基板界面に到達して反射された光が再度多層膜で反射される際の多層膜反射率を求める。図１１、図１２に、波長λの光が角度θ０（θ，λ）で基板から多層膜に入射した場合を示す。各界面での電場をＥ（１）ｓｊα（θ，λ）、Ｅ（１）ｐｊα（θ，λ）で表す。
【００５４】
次に下記（５）（６）の場合における方程式（数式（４４）から数式（５２））を満たすＥ（１）ｓ００（θ，λ）、Ｅ（１）ｓ０１（θ，λ）、Ｅ（１）ｐ００（θ，λ）、Ｅ（１）ｐ０１（θ，λ）を求める。Ｅ（１）ｓ０１（θ，λ）、Ｅ（１）ｐ０１（θ，λ）は多層膜への入射光の電場を表し、Ｅ（１）ｓ００（θ，λ）、Ｅ（１）ｐ００（θ，λ）は多層膜で反射された反射光の電場を表す。
【００５５】
【数１１】
2007095500.tif 000013
【００５６】
【数１２】
2007095500.tif 000014
【００５７】
角度θ０（θ，λ）で基板内をｚ軸の正の方向に向う波長λのｓ偏光の光の多層膜での反射率ｒｓ（θ，λ）、ｐ偏光の光の多層膜での反射率ｒｐ（θ，λ）は数式（５３）から数式（５４）のようになる。
【００５８】
【数１３】
2007095500.tif 000015
【００５９】
基板の厚さは通常光の波長よりも遥かに大きいため外部／基板界面で一度以上反射された光は他の光との干渉性がなくなると考えてよい。外部／基板界面と多層膜とで順次反射された光はインコヒーレントに加えることができ、全経路を足し合わせるとθ方向の発光スペクトルＩＥＬ’（λ）は数式（５５）で表すことができる。ここでＩＦ（λ）は発光層から発する光のスペクトルである。
【００６０】
【数１４】
2007095500.tif 000016
【００６１】
ｎ０（λ）ｃｏｓθ０（θ，λ）／ｎｈ（λ）ｃｏｓθｈ（θ，λ）の項が入る理由について以下に説明する。発光界面においてθｈ（θ，λ）方向へ湧き出すＳ偏光の光の電場の振幅はｖｓ０であるので、進行方向に垂直な面を通過する単位面積あたりのエネルギーの時間平均はポインティングベクトルの時間平均であり数式（５６）に比例する。
【００６２】
【数１５】
2007095500.tif 000017
【００６３】
一方、基板に取り出されたＳ偏光の光の、進行方向（θ０（θ，λ）方向）に垂直な面を通過する単位面積あたりのエネルギーの時間平均はポインティングベクトルの時間平均であり数式（５７）で表される。
【００６４】
【数１６】
2007095500.tif 000018
【００６５】
また、図１３に示すように、発光分子から放出された光の光線としての太さを考慮すると、発光界面においてひとつの発光分子からθｈ（θ，λ）方向へ放出される光の光線の太さはｃｏｓθｈ（θ，λ）に比例し、基板内をｃｏｓθ０（θ，λ）方向へ進行する光の光線の太さはｃｏｓθ０（θ，λ）に比例するので、ひとつの発光分子からθｈ（θ，λ）方向へ放出される光のエネルギー（時間平均）は数式（５８）に比例する。
【００６６】
【数１７】
2007095500.tif 000019
【００６７】
ひとつの発光分子から放出され基板内に到達する光エネルギー（時間平均）は数式（５９）に比例する。従ってひとつの発光分子から放出された光エネルギーのうち基板内に到達する光エネルギーの割合は、数式（６０）に比例する。
【００６８】
【数１８】
2007095500.tif 000020
【００６９】
ここで発光分子から放出される光エネルギーの波長分散を考え、発光波長オーダーの距離にミラーの存在しない、あるいはミラーによる反射が微小な位置に置かれた発光材料から発する光スペクトルをＩＦ（λ）とすると、上記有機電界発光素子内のひとつの発光分子からθｈ（θ，λ）方向（ｋｈ０方向）へ放出される光のエネルギーはその遷移双極子モーメントの方向まで考慮するとＩＦ（λ）｜ｖｓ０｜２に比例する。従って、基板に取り出されるＳ偏光の光のスペクトルは数式（６１）に比例する。
【００７０】
【数１９】
2007095500.tif 000021
【００７１】
これらはＰ偏光においても同様である。以上を考慮するとひとつの発光分子から放出され外部に到達する波長λの光のエネルギーの時間平均（発光スペクトル）は、数式（５５）で表すことができる。
【００７２】
ところで、発光界面に多数の発光分子が並んでいる場合、図１４に示すように、同じスポット径が見込む発光分子数は角度θｈ（θ，λ）に応じて変化し、１／ｃｏｓθｈ（θ，λ）に比例する。従って、全体としての発光スペクトルＩＥＬ（λ）は数式（６２）で表すことができる。
【００７３】
【数２０】
2007095500.tif 000022
【００７４】
発光分子の向きが様々であり、遷移双極子モーメントが様々な方向を向いている場合は以下のようにする。まず、任意の遷移双極子モーメントの方向M（θＴＭ，ψＴＭ）を向いた発光分子数をｇ（θＴＭ，ψＴＭ）とおく。遷移双極子モーメントの向きをM（θＴＭ，ψＴＭ）とした場合について数式（６２）を用いてθ方向の発光スペクトルを求めてＩＥＬ’（θ，θＴＭ，ψＴＭ，λ）とおく。さらに数式（６３）を用いて角度方向で積分し、θ方向の発光スペクトルＩＥＬ（θ，λ）を求める。
【００７５】
【数２１】
2007095500.tif 000023
【００７６】
発光分子が配向しておらず、向きがランダムな場合はｇ（θＴＭ，ψＴＭ）を定数とすればよい。
【００７７】
発光領域が厚さ方向（Ｚ軸方向）に広がりを有している場合は以下のようにする。発光領域が第（ｈ−３）層／第（ｈ−２）層界面から第ｈ層／第（ｈ＋１）層界面の間に広がっているとし、この領域内の任意の面ｄｈ−２｛第（ｈ−１）層／第ｈ層界面｝における発光強度をｆ（ｄｈ−２）とおく。｛ｄｈ−２は第（ｈ−３）層／第（ｈ−２）層界面から第（ｈ−１）層／第ｈ層界面までの距離である。ここでは第（ｈ−１）層／第ｈ層界面を可動としている。｝発光面｛前記発光領域が第（ｈ−１）層／第ｈ層界面に限られる場合の第（ｈ−１）層／第ｈ層界面｝を数式（６３）の領域内の面ｄｈ−２とした場合について、数式（６３）を用いてθ方向の発光スペクトルを求めてＩＥＬ’(θ，ｄｈ−２，λ)とおく。さらに数式（６４）を用いて発光領域内で積分し、θ方向の発光スペクトルＩＥＬ（θ，λ）を求める。ここでｄｅｍｌは発光層のトータル膜厚であり、ｄｅｍｌ＝ｄｈ−２＋ｄｈ−１＋ｄｈ＋ｄｈ＋１である。
【００７８】
【数２２】
2007095500.tif 000024
【００７９】
遷移双極子モーメントの向きおよび発光領域に関する積分をまとめると以下のようである。すなわち、数式（６２）を用いて面ｄｈ−２における遷移双極子モーメントM（θＴＭ，ψＴＭ）の発光分子からの発光スペクトルを求めてＩＥＬ’（θ，θＴＭ，ψＴＭ，ｄｈ−２，λ）とおく。さらに数式（６５）を用いて角度方向および発光領域内で積分し、θ方向の発光スペクトルＩＥＬ（θ，λ）を求める。
【００８０】
【数２３】
2007095500.tif 000025
【００８１】
発光スペクトルＩＥＬ（θ，λ）から数式（６６）から数式（６９）を用いてＣＩＥ１９３１表色系における３刺激値Ｘ（θ）、Ｙ（θ）、Ｚ（θ）を求め、θ方向の色度座標（ｘ（θ），ｙ（θ））、相対輝度Ｌ（θ）を求めることができる。なおバー付きのｘ（λ），ｙ（λ），ｚ（λ）はＣＩＥ１９３１表色系における等色関数である。
【００８２】
【数２４】
2007095500.tif 000026
【００８３】
次に本発明である、励起状態から基底状態へ光子を放出して遷移する単位時間あたりの確率（輻射遷移の失活速度定数）を求め、さらに励起状態から基底状態へ無輻射で遷移する単位時間あたりの確率（無輻射遷移の失活速度定数）を考慮するモデル（輻射遷移の失活速度定数算出モデル）について述べる。
【００８４】
発光分子の遷移双極子モーメントが図１０のようにM（θＴＭ，ψＴＭ）方向を向いている場合を考える。発光波長オーダーの距離にミラーの存在しない、あるいはミラーによる反射が微小な位置に発光分子が置かれた場合、発光分子が励起状態にある時に、ｋｈ０方向を含む微小立体角ｄΩ内にS偏光の光子を放出して基底状態に遷移する単位時間あたりの確率は、数式（２０）から数式（２３）に記載のｖｓ０を用いて表すと、｜ｖｓ０｜２ｄΩ／２に比例する。同様にｋｈ１方向を含む微小立体角ｄΩ内にS偏光の光子を放出して基底状態に遷移する単位時間あたりの確率は、｜ｖｓ１｜２ｄΩ／２に比例する。また、ｋｈ０方向を含む微小立体角ｄΩ内にＰ偏光の光子を放出して基底状態に遷移する単位時間あたりの確率は、｜ｖｐ０｜２ｄΩ／２に比例する。同様にｋｈ１方向を含む微小立体角ｄΩ内にＰ偏光の光子を放出して基底状態に遷移する単位時間あたりの確率は、｜ｖｐ１｜２ｄΩ／２に比例する。
【００８５】
以上はミラーが存在しない場合であるが、発光波長オーダーの距離にミラーが存在する場合、一般にこれら輻射遷移の単位時間当たりの確率は変化する。前述の双極子モデルにおいても実はこれらの変化が盛り込まれている。双極子モデルは、下記に述べる輻射遷移の失活速度定数算出モデルにおいて輻射遷移の発光失活速度定数よりも無輻射遷移の発光失活速度定数が非常に大きい場合、すなわち発光量子収率が非常に小さい場合の極限にあたる。
【００８６】
まず、発光分子から外部に放出される経路に光子が放出される場合を考える。すなわちｋｈ０方向に放出された光子は素子内部で吸収されなければ、素子外部に取り出されることが可能な経路であるとする。この時、双極子モデルで行ったように、まず、発光分子から数式（２０）から数式（２３）のような電場振幅をもった光波が放出されるとして計算を行う。計算に必要なのはｋｈ０およびｋｈ１方向に放出される光波である。これらはミラーが存在すると反射して互いに同方向になり干渉するからである。
【００８７】
発光分子からｋｈ０方向、ｋｈ１方向（θｈ方向）に放出された光による、波長λの光の各界面での電場をＥｓｊα（θｈ，λ）、Ｅｐｊα（θｈ，λ）で表す（図６、図７参照）｛ここでｓ、ｐはそれぞれｓ偏光、ｐ偏光を表す（ｓ偏光とは電場がｙｚ面に垂直な光であり、ｐ偏光とは磁場がｙｚ面に垂直な光である）。なお、下付き添え字ｊは層番号を表す。また、上付き添え字αは０か１の値を持ち、０はｚ軸の負の方向に向う光を、１はｚ軸の正の方向に向う光を示す｝。θ方向の波長λの取り出し光に寄与する光の各層での方向θｊ（θｈ，λ）に対して、ｓｉｎθｊ（θｈ，λ）、ｃｏｓθｊ（θｈ，λ）を数式（７０）、数式（７１）を用いて求める。
【００８８】
【数２５】
2007095500.tif 000027
【００８９】
次に発光分子から発せられる光の振幅を求める。ｋｈ０方向に湧き出すＳ偏光の光の電場振幅をｖｓ０、Ｐ偏光の光の電場振幅をｖｐ０、ｋｈ１方向に湧き出すＳ偏光の光の電場振幅をｖｓ１、Ｐ偏光の光の電場振幅をｖｐ１として数式（２０）から数式（２３）のようになる。
【００９０】
次に下記（１）〜（４）の場合における方程式（数式（７４）から数式（８９））を満たすａｓ（θｈ，λ）、ａｐ（θｈ，λ）を求める。
【００９１】
【数２６】
2007095500.tif 000028
【００９２】
【数２７】
2007095500.tif 000029
【００９３】
数式（７４）から数式（８９）を用いるとＥｓ（ｈ−１）０（θｈ，λ）、Ｅｓ（ｈ−１）１（θｈ，λ）、Ｅｐ（ｈ−１）０（θｈ，λ）、Ｅｐ（ｈ−１）１（θｈ，λ）が求まる。発光分子からｋｈ０方向、ｋｈ１方向への光の湧き出しエネルギーは、発光界面の取り出し側の面およびその反対面でポインティングベクトルの時間平均を差し引けば求まる。ｋｈ０方向へ湧き出す波長λのS偏光の光のエネルギーをＳｓ０（θｈ，θＴＭ，ψＴＭ，λ）、ｋｈ１方向へ湧き出す波長λのS偏光の光のエネルギーをＳｓ１（θｈ，θＴＭ，ψＴＭ，λ）、ｋｈ０方向へ湧き出す波長λのP偏光の光のエネルギーをＳｐ０（θｈ，θＴＭ，ψＴＭ，λ）、ｋｈ１方向へ湧き出す波長λのP偏光の光のエネルギーをＳｐ１（θｈ，θＴＭ，ψＴＭ，λ）とすると、
【００９４】
【数２８】
2007095500.tif 000030
【００９５】
と書くことができる。なお、ここではエネルギーの波長依存性をなくすために、ポインティングベクトルであれば本来掛けておかなければならない発光層の屈折率ｎｈを右辺に含めていない。（ｎｈは一般に波長に応じて変化する。）また、数式（２０）から数式（２３）のｖｓ０等をθｈ，θＴＭ，ψＴＭの関数をしてｖｓ０（θｈ，θＴＭ，ψＴＭ）等と記した。以降もこの表記を用いる。
【００９６】
ミラーが存在しない時に発光材料から発せられる発光スペクトルをＩＦ（λ）とすると、有機電界発光素子内の発光分子から実際に湧き出す光のエネルギーは数式（９０）から数式（９３）にＩＦ（λ）を掛けたものに比例すると考えることができる。また、さらにｈｃ／λで割るとエネルギーを光子数に変換することができ、これが励起状態にある発光分子がｋｈ０方向、およびｋｈ１方向に波長λの光子を放出して基底状態に遷移する単位時間あたりの確率に比例するものとなる。
【００９７】
以上をまとめると以下のようになる。すなわち、遷移双極子モーメントがM（θＴＭ，ψＴＭ）を向いた励起状態にある発光分子から、角度ｋｈ０方向を含む微小立体角ｄΩ内へ波長λからλ＋ｄλのS偏光の光子を放出して基底状態に遷移する単位時間あたりの確率をｋｒｓ０（θｈ，θＴＭ，ψＴＭ，λ）ｄΩｄλ、角度ｋｈ１方向を含む微小立体角ｄΩ内へ波長λからλ＋ｄλのS偏光の光子を放出して基底状態に遷移する単位時間あたりの確率をｋｒｓ１（θｈ，θＴＭ，ψＴＭ，λ）ｄΩｄλ、角度ｋｈ０方向を含む微小立体角ｄΩ内へ波長λからλ＋ｄλのＰ偏光の光子を放出して基底状態に遷移する単位時間あたりの確率をｋｒｐ０（θｈ，θＴＭ，ψＴＭ，λ）ｄΩｄλ、角度ｋｈ１方向を含む微小立体角ｄΩ内へ波長λからλ＋ｄλのＰ偏光の光子を放出して基底状態に遷移する単位時間あたりの確率をｋｒｐ１（θｈ，θＴＭ，ψＴＭ，λ）ｄΩｄλと記載すると下記数式（９４）から数式（９７）のようになる。
【００９８】
以上は外部に放出される経路に光子が放出される場合（外部放出モード）であったが、基板内には到達するが基板／外部界面で全反射を繰り返し基板内を横方向に導波する経路の場合（基板導波モード）についても同様に扱うことができる。さらに本発明においては、発光層を含む薄膜内部に閉じ込められ、基板／第1電極および第2電極で反射を繰り返し、薄膜層内を横方向に導波する経路の場合（薄膜導波モード）についても同様に扱う。
【００９９】
【数２９】
2007095500.tif 000031
【０１００】
遷移双極子モーメントがM（θＴＭ，ψＴＭ）を向いた励起状態にある発光分子がその方向を問わず光子を放出して基底状態に遷移する単位時間あたりの確率ｋｒ（θＴＭ）（輻射遷移の失活速度定数）は数式（９８）となる。
【０１０１】
ｋｒ（θＴＭ）を求めるには発光分子の遷移双極子モーメントの向きをある方向に固定してｋｈ０およびｋｈ１の方向を様々に変化させて、各方向へ輻射遷移の失活速度定数を求めて足し合わせればよいが、ここでは以下のように計算を行った。すなわち、ｋｈ０およびｋｈ１に関してはｚｘ平面内に固定しθhの値を０〜９０°まで変化させ、ｚ軸を軸に回転はさせないが、代わりに遷移双極子モーメントMをｚ軸を軸に３６０°回転させる。また、輻射遷移の失活速度定数はθＴＭのみに依存しψＴＭには依存しないのでｋｒ（θＴＭ）と記載した。
【０１０２】
発光分子が励起状態から無輻射で基底状態に遷移する単位時間あたりの確率（無輻射遷移の失活速度定数）をｋｔとすると、遷移双極子モーメントがM（θＴＭ，ψＴＭ）を向いた発光分子が励起状態になった後に光子を放出する確率φａｌｌ(θＴＭ)（発光量子収率）は数式（９９）となる。
【０１０３】
【数３０】
2007095500.tif 000032
【０１０４】
有機電界発光素子内の発光分子の発光量子収率φａｌｌ（θＴＭ）は発光分子が発光波長オーダーの距離にミラーの存在しない、あるいはミラーによる反射が微小な位置に置かれた場合の発光量子収率φｏｕｔと一般に値は異なっている。なお、本発明においては無輻射遷移の失活速度定数ｋｔの値は、発光分子が光励起された場合と、有機電界発光素子内で正孔と電子が再結合することにより励起された場合とで値は変化しないとして話を進める。
【０１０５】
ｋｔの値を定めるには例えば以下のようにすればよい。すなわち、発光分子が発光波長オーダーの距離にミラーの存在しない、あるいはミラーによる反射が微小な位置に置かれた場合の輻射遷移の失活速度定数ｋｒｏｕｔを数式（１００）により求める。
【０１０６】
【数３１】
2007095500.tif 000033
【０１０７】
さらに発光分子が発光波長オーダーの距離にミラーの存在しない、あるいはミラーによる反射が微小な位置に置かれた場合の発光量子収率φｏｕｔから数式（１０１）によりｋｔを求める。
【０１０８】
【数３２】
2007095500.tif 000034
【０１０９】
遷移双極子モーメントがM（θＴＭ，ψＴＭ）を向いた励起状態にある発光分子が角度ｋｈ０方向を含む微小立体角ｄΩ内へ波長λからλ＋ｄλの光子を放出する確率φ０（θｈ，θＴＭ，ψＴＭ，λ）ｄΩｄλは、
【０１１０】
【数３３】
2007095500.tif 000035
【０１１１】
同様に、遷移双極子モーメントがM（θＴＭ，ψＴＭ）を向いた励起状態にある発光分子が角度ｋｈ１方向を含む微小立体角ｄΩ内へ波長λからλ＋ｄλの光子を放出する確率φ１（θｈ，θＴＭ，ψＴＭ，λ）ｄΩｄλは、以下のようになる。
【０１１２】
【数３４】
2007095500.tif 000036
【０１１３】
ｋｈ０方向が取り出し側の外部放出モードである場合、遷移双極子モーメントがM（θＴＭ，ψＴＭ）を向いた励起状態にある発光分子が角度ｋｈ０方向を含む微小立体角ｄΩ内へ波長λからλ＋ｄλの光子を放出しその光子が基板に到達する確率φｓｕｂ（θｈ，θＴＭ，ψＴＭ，λ）ｄΩｄλは数式（１０４）のようになる。
【０１１４】
【数３５】
2007095500.tif 000037
【０１１５】
ここでｎ０（λ）ｃｏｓθ０（θｈ，λ）／ｎｈ（λ）ｃｏｓθｈの項が含まれることについて説明をする。ｃｏｓθ０（θｈ，λ）／ｃｏｓθｈの項は図１３に示すように屈折により光線の太さが変化することにより挿入される項である。そしてポインティングベクトルの時間平均が例えばS偏光の基板取り出し光、光源の湧き出し光について、それぞれ数式（１０５）、数式（１０６）
【０１１６】
【数３６】
2007095500.tif 000038
【０１１７】
で表されるため、光源でｋｈ０方向に湧き出した波長λのＳ偏光の光エネルギーのうちどれだけの割合が基板に到達するかを計算するには数式（１０７）を求めればよい。Ｐ偏光の光についても同様である（数式（１０８））。
【０１１８】
【数３７】
2007095500.tif 000039
【０１１９】
以上の考察から数式（１０４）にｎ０（λ）ｃｏｓθ０（θｈ，λ）／ｎｈ（λ）ｃｏｓθｈの項が含まれることが説明された。
【０１２０】
基板に到達した光は全てが外部に到達するわけではなく、一部は基板／外部界面で反射されて素子内部に戻り、戻った光のうち一部がまた再度多層膜で反射されて基板に戻る。これらが繰り返されることになる。
【０１２１】
そこで、双極子モデルで述べたのと同様に、数式（１０９）から数式（１１２）を用いて角度θ０（θｈ，λ）で基板内をｚ軸の負の方向に向かう波長λのｓ偏光の光の、基板から外部への透過率Ｔｓ（θｈ，λ）、反射率Ｒｓ（θｈ，λ）、ｐ偏光の光の、基板から外部への透過率Ｔｐ（θｈ，λ）、反射率Ｒｐ（θｈ，λ）を求める。
【０１２２】
【数３８】
2007095500.tif 000040
【０１２３】
次に、外部／基板界面に到達して反射された光が再度多層膜で反射される際の多層膜反射率を求める。波長λの光が角度θ０（θｈ，λ）で基板から多層膜に入射した場合を示す。各界面での電場をＥ（１）ｓｊα（θｈ，λ）、Ｅ（１）ｐｊα（θｈ，λ）で表す（図１１、図１２）。
【０１２４】
次に下記（５）（６）の場合における方程式（数式（１１３）から数式（１２１））を満たすＥ（１）ｓ００（θｈ，λ）、Ｅ（１）ｓ０１（θｈ，λ）、Ｅ（１）ｐ００（θｈ，λ）、Ｅ（１）ｐ０１（θｈ，λ）を求める。Ｅ（１）ｓ０１（θｈ，λ）、Ｅ（１）ｐ０１（θｈ，λ）は多層膜への入射光の電場を表し、Ｅ（１）ｓ００（θｈ，λ）、Ｅ（１）ｐ００（θｈ，λ）は多層膜で反射された反射光の電場を表す。
【０１２５】
【数３９】
2007095500.tif 000041
【０１２６】
【数４０】
2007095500.tif 000042
【０１２７】
角度θ０（θｈ，λ）で基板内をｚ軸の正の方向に向う波長λのｓ偏光の光の多層膜での反射率ｒｓ（θｈ，λ）、ｐ偏光の光の多層膜での反射率ｒｐ（θｈ，λ）は数式（１２２）から数式（１２３）のようになる。
【０１２８】
【数４１】
2007095500.tif 000043
【０１２９】
これらを考慮すると、遷移双極子モーメントがM（θＴＭ，ψＴＭ）を向いた励起状態にある発光分子が角度ｋｈ０方向（θｈ）を含む微小立体角ｄΩ内へ波長λからλ＋ｄλの光子を放出しその光子が外部に到達する確率φＥＬ（θｈ，θＴＭ，ψＴＭ，λ）ｄΩｄλは数式（１２４）のようになる。ここで外部の屈折率ｎ＝１、消衰係数ｋ＝０とした。
【０１３０】
【数４２】
2007095500.tif 000044
【０１３１】
φＥＬ（θｈ，θＴＭ，ψＴＭ，λ）は光子数表示のＥＬスペクトルに対応する。これを通常の表示であるエネルギー表示のＥＬスペクトルＩＥＬ（θｈ，θＴＭ，ψＴＭ，λ）に変換するには、以下のようにすればよい。
【０１３２】
【数４３】
2007095500.tif 000045
【０１３３】
発光界面の光の進行方向がθｈの場合、外部出射光の進行方向はθ（θｈ，λ）となり、波長によって出射角度θ（θｈ，λ）が異なるが、外部出射光の角度をθに固定するなら、発光界面の光の進行方向はθｈ（θ，λ）と書くことができる。改めて書き直すと、遷移モーメントがM（θＴＭ，ψＴＭ）を向いた発光分子から角度θ方向に光が取り出される時、そのＥＬスペクトルは、数式（１２６）のようになる。
【０１３４】
【数４４】
2007095500.tif 000046
【０１３５】
発光界面に多数の発光分子が並んでいる場合、図１４に示すように、同じスポット径が見込む発光分子数は角度θｈ（θ，λ）に応じて変化し、１／ｃｏｓθｈ（θ，λ）に比例する。従って、数式（１２６）のひとつの発光分子による発光スペクトルをＩＥＬ’（θ，θＴＭ，ψＴＭ，λ）と書けば、全体としての発光スペクトルＩＥＬ（θ，θＴＭ，ψＴＭ，λ）は数式（１２７）で表すことができる。
【０１３６】
【数４５】
2007095500.tif 000047
【０１３７】
発光分子の向きが様々であり、遷移双極子モーメントが様々な方向を向いている場合は以下のようにする。
【０１３８】
まず、任意の遷移双極子モーメントの方向M（θＴＭ，ψＴＭ）を向いた発光分子数をｇ（θＴＭ，ψＴＭ）とおく。遷移双極子モーメントの向きをM（θＴＭ，ψＴＭ）とした場合について数式（１２７）を用いてθ方向の発光スペクトルを求めてＩＥＬ’（θ，θＴＭ，ψＴＭ，λ）とおく。さらに数式（１２８）を用いて角度方向で積分し、θ方向の発光スペクトルＩＥＬ（θ，λ）を求める。
【０１３９】
【数４６】
2007095500.tif 000048
【０１４０】
発光分子が配向しておらず、向きがランダムな場合はｇ（θＴＭ，ψＴＭ）を定数とすればよい。
【０１４１】
発光領域が厚さ方向（Ｚ軸方向）に広がりを有している場合は以下のようにする。発光領域が第（ｈ−３）層／第（ｈ−２）層界面から第ｈ層／第（ｈ＋１）層界面の間に広がっているとし、この領域内の任意の面ｄｈ−２｛第（ｈ−１）層／第ｈ層界面｝における発光強度をｆ（ｄｈ−２）とおく。｛ｄｈ−２は第（ｈ−３）層／第（ｈ−２）層界面から第（ｈ−１）層／第ｈ層界面までの距離である。ここでは第（ｈ−１）層／第ｈ層界面を可動としている。｝発光面｛前記発光領域が第（ｈ−１）層／第ｈ層界面に限られる場合の第（ｈ−１）層／第ｈ層界面界面｝を発光領域内の面ｄｈ−２とした場合について、数式（１２５）を用いてθ方向の発光スペクトルを求めてＩＥＬ’(θ，ｄｈ−２，λ)とおく。さらに数式（１２６）を用いて発光領域内で積分し、θ方向の発光スペクトルＩＥＬ（θ，λ）を求める。ここでｄｅｍｌは発光層のトータル膜厚であり、ｄｅｍｌ＝ｄｈ−２＋ｄｈ−１＋ｄｈ＋ｄｈ＋１である。
【０１４２】
【数４７】
2007095500.tif 000049
【０１４３】
遷移双極子モーメントの向きおよび発光領域に関する積分をまとめると以下のようである。すなわち、数式（１２７）を用いて面ｄｈ−２における遷移双極子モーメントM（θＴＭ，ψＴＭ）の発光分子からの発光スペクトルを求めてＩＥＬ’（θ，θＴＭ，ψＴＭ，ｄｈ−２，λ）とおく。さらに数式（１３０）を用いて角度方向および発光領域内で積分し、θ方向の発光スペクトルＩＥＬ（θ，λ）を求める。
【０１４４】
【数４８】
2007095500.tif 000050
【０１４５】
発光スペクトルＩＥＬ（θ，λ）から数式（１３１）から数式（１３４）（数５９参照）を用いてＣＩＥ１９３１表色系における３刺激値Ｘ（θ）、Ｙ（θ）、Ｚ（θ）を求め、θ方向の色度座標（ｘ（θ），ｙ（θ））、相対輝度Ｌ（θ）を求めることができる。なおバー付きのｘ（λ），ｙ（λ），ｚ（λ）はＣＩＥ１９３１表色系における等色関数である。
【０１４６】
【数４９】
2007095500.tif 000051
【０１４７】
本発明は透明基板上に透明な第一電極、有機化合物からなる薄膜層、光反射層が積層された構造の素子において効果的であるが、基板上に光反射層、有機化合物からなる薄膜層、透明な第二電極が積層されたトップエミッション構造の有機電界発光素子に対しても効果的である。
【０１４８】
以上のモデルを用いて実際にシミュレーションするには、コンピュータを用いて本発明のシミュレーション方法を反映したプログラムを実行すればよい。具体的には、例えば以下のようにして実施する。まず、有機電解発光素子を構成する各層の膜厚、光学定数、発光層内の膜厚方向の発光強度分布、発光層から発せられる光スペクトル、発光分子数の方向分布、着目する視野角をキーボード等の入力手段により数値データとしてコンピュータのハードディスク内に記憶させる。
【０１４９】
次にあらかじめハードディスク内に記憶させている本発明のシミュレーションを反映したプログラムを起動し、ハードディスク内に記憶されたデータを使用し、計算処理を実行させる。計算処理は、コンピュータ内のハードディスク、メモリ、ＣＰＵを用いて行われるが、その内容はおおよそ次の通りである。まず、発光層内のある位置におけるある方向を向いた発光分子による着目する視野角方向の発光スペクトルを求め、そこにその発光分子の存在確立を乗じる。これを種々の位置、種々の方向の発光分子について順次求めて行き、最後に、求めたものをすべて足し合わせることにより、素子の着目する視野角方向の発光スペクトルが得られる。また、その発光スペクトルから輝度および色度が得られる。得られた結果については、モニタ等の出力手段により表示させることで確認できる。また、本発明は発光量子収率の高い発光材料を用いた素子の発光特性のシミュレーションに有効である。
【０１５０】
なお、本発明による数値計算を行う際には、以下の２点の近似を行うこととした。１点目は、発光分子が励起状態からある方向へ光子を放出して基底状態へ遷移する単位時間あたりの確率は、発光分子が波長オーダーの距離にミラーの存在しない、あるいはミラーによる反射が微小な場合の遷移確率の２５倍以下とする近似である。図１５に示すように符号２５の位置にある発光分子が反射率１００％の完全なミラーに挟まれた状態にあり、かつ点線で示す右側に放射された光波がミラーに反射されて左側に進行する際、実線で示す左側に放射された光波と位相が一致し強めあう場合、光の吸収がないとすると図の符号２５の位置において初めの光波の電場振幅はＥであったが、図の符号２６の位置において５Ｅになる。図には記載していないが、反射を繰り返すにつれ電場振幅はより大きくなって行く。もしミラーが存在しなければ、発光分子から図の左の角度方向へ光子が放出される単位時間あたりの遷移確率（失活速度）は｜Ｅ｜２に比例するが、ミラーが存在する場合、符号２６の位置で観測を行うと、左の角度方向へ光子が放出される失活速度は２５｜Ｅ｜２となる。観測位置が図の上の方になる程、失活速度は大きくなって行く。本発明においては、実際の有機電界発光素子においては大きく強めあうような理想的な干渉が起こらないであろうと推測し、ミラーが存在しない場合の失活速度の２５倍より大きくはならないと近似した。なお、これはＳ偏光の光について、Ｐ偏光の光についてそれぞれ独立に近似した。具体的には下記数式（１３５）の場合数式（１３６）、数式（１３７）の場合数式（１３８）、数式（１３９）の場合数式（１４０）、数式（１４１）の場合数式（１４２）とするよう近似した。
【０１５１】
【数５０】
2007095500.tif 000052
【０１５２】
２点目は、発光分子が励起状態からθｈの大きな方向すなわち基板法線方向に対して垂直に近い方向へ光子を放出して基底状態へ遷移する単位時間あたりの確率は、発光分子が波長オーダーの距離にミラーの存在しない、あるいはミラーによる反射が微小な場合の遷移確率の１倍以下とする近似である。図１６に示すようにθｈが大きい場合発光分子から放出された光波はミラーで反射されるまでに長い距離を進む必要があり、ミラーで反射されるまでにホスト分子に吸収される確率が高くなる。ミラーで反射されなければ干渉も起きず、その方向への輻射遷移の失活速度が変化することもない。従って、θｈが大きい場合その方向への輻射遷移の失活速度はミラーがない場合と同じになる可能性が高く、ここではθｈ≧８５°の場合に、ｋｈ０方向またはｋｈ１方向のどちらかへの輻射遷移の失活速度がミラーがない場合の速度を上回る場合に、ｋｈ０方向およびｋｈ１方向両方についてその輻射遷移の失活速度をミラーがない場合の速度と同じ値になるよう修正することとした。
【０１５３】
具体的には下記数式（１４３）または数式（１４４）の場合、数式（１４５）および数式（１４６）とするように、また、下記数式（１４７）または数式（１４８）の場合、数式（１４９）および数式（１５０）とするように近似した。
【０１５４】
【数５１】
2007095500.tif 000053
【実施例】
【０１５５】
以下、実施例および比較例をあげて本発明を説明するが、本発明はこれらの例によって限定されるものではない。
【０１５６】
スパッタリング法によりＩＴＯ透明導電膜を１０９ｎｍ積層したガラス基板を３８×４６ｍｍに切断した後、ＩＴＯの不要部分をエッチング除去した。得られた基板をアルカリ洗浄液で１０分間超音波洗浄してから、超純水で洗浄した。この基板を、素子を作製する直前に１時間ＵＶ／オゾン処理し、真空装置内に設置して、装置内の真空度が５×１０−４Ｐａ以下になるまで排気した。抵抗加熱法によって、まず銅フタロシアニン（ＣｕＰｃ）を６ｎｍ蒸着し、続いて正孔輸送材料としてＮ，Ｎ’−ジ−（ナフタレン−１−イル）−Ｎ，Ｎ’−ジフェニル−ベンジジン（ＮＰＤ）を６２ｎｍ蒸着し、続いて発光層として、ホスト材料トリス（８-キノリノラト）アルミニウム（Ａｌｑ３）に、ゲスト材料１０−（２−ベンゾチアゾリル）−１，１，７，７−テトラメチル−２，３，６，７−テトラヒドロ−１Ｈ，５Ｈ，１１Ｈ−〔ｌ〕ベンゾピラノ〔６，７，８−ｉｊ〕キノリジン−１１−オン（クマリン５４５T）を０．５ｗｔ％ドーピングした膜を２５ｎｍ蒸着し、続いて電子輸送材料として上記に示した化１で表される化合物を４．７５ｎｍ蒸着した。陰極用マスクを装着し、フッ化リチウムの蒸気にさらしてドーピングした後、アルミニウムを１５０ｎｍ蒸着して陰極とした。このようにして作製した有機電界発光素子を素子１と呼ぶ。電子輸送材料の蒸着膜厚を下記表１のように４．７５ｎｍ〜２０９ｎｍの間の値に変えた以外は同様にして素子２〜素子２０を作製した。
【０１５７】
【表１】
2007095500.tif 000054
【０１５８】
このようにして作製した有機電界発光素子（素子１〜素子２０）を１０ｍＡ／ｃｍ２の電流密度で発光させ、基板表面の法線方向の輝度を測定したところ、図１７の（１）のようであった。輝度は電子輸送層膜厚０〜１００ｎｍにおいて最も高い輝度を示した素子の輝度を１として規格化している。
【０１５９】
実施例１
上記構成の有機電界発光素子（素子１〜素子２０）について本発明のモデルによるシミュレーションを行い、ミラーが存在しない時の発光層の発光量子効率が０．９９９であるとして、基板表面の法線方向の輝度を求めたところ図１７の（２）のようになり、実際の発光素子に近い数値が得られた。電子輸送層膜厚を厚くしていった時に、輝度は一度高くなり（第１のピーク）、より厚くすれば次第に低くなってゆく。さらに厚くすると、再び高くなって行き第２のピークが現れる。実測データでは第２ピークが第１ピークより高くなっているが、実施例においても同様に第２ピークが第１ピークより高くなっている。なお、輝度は電子輸送層膜厚０〜１００ｎｍにおいて最も高い輝度を示した素子の輝度を１として規格化して示した。
【０１６０】
ここで発光分子の向きはランダムとして数式（１２８）におけるｇ（θＴＭ，ψＴＭ）＝１とし、θＴＭについてはθＴＭ＝０、１５、３０、４５、６０、７５、９０°の場合を、ψＴＭについてはψＴＭ＝０、１０、２０、３０、４０、５０、６０、７０、８０、９０、１００、１１０、１２０、１３０、１４０、１５０、１６０、１７０、１８０°の場合をそれぞれ独立に計算して足し合わせた。空間的対称性によりψＴＭの範囲は０〜１８０°でも０〜３６０°と同じ結果が得られるので０〜１８０°の範囲とした。また、ここで発光層内での発光強度分布ｆ（ｄｈ−２）はＮＰＤ／発光層界面で強く発光層／電子輸送層界面に近づくに従って減少し、膜厚に関する指数関数分布として数式（１５１）のように表せると仮定した。
【０１６１】
【数５２】
2007095500.tif 000055
【０１６２】
ここでｆ０は定数であり、ｄｈ−２はＮＰＤ／発光層界面を０として発光層／電子輸送層界面に向って測った距離のｎｍ単位の数字である。なお、発光スペクトルは最終的には規格化されるからｆ０は如何なる値であってもよいが、ここでは１とした。また指数関数の分母の１０という値は、発光領域の広がりが１０ｎｍ程度であると仮定することにより定めた。
【０１６３】
すなわちＮＰＤ／発光層界面から発光層／電子輸送層界面に向って１０ｎｍ進んだ地点で発光強度は１／ｅに減少するとした。計算では発光層を等間隔に４つに分割し、計５面が発光界面である場合をそれぞれ独立に計算して足し合わせた。
【０１６４】
なお、波長については４００〜８００ｎｍについて５ｎｍ毎に計算を行った。また発光面から放射される光スペクトルとして、石英板上にＡｌｑ３にクマリン５４５Ｔを０．５ｗｔ％ドーピングした膜を２５ｎｍ積層し、光励起スペクトルを測定して得られたスペクトル（図１８）を用いた。また、有機電界発光素子を構成する各層の光学定数として、ガラス基板上に各層をそれぞれ１００ｎｍ積層し偏光解析により光学定数を測定して得られた値を用いた。
【０１６５】
比較例１
この構成の有機電界発光素子について、特許文献１の段落００６６〜０１０３に記載の方法でシミュレーションを行い、素子１〜素子２０の基板表面の法線方向の輝度を求めたところ図１７の（３）のようになった。なお、光反射層界面全体の面積に対して正反射が起こる面積の比ｒ２は特許文献１に記載の０．８７とした。輝度は電子輸送層膜厚０〜１００ｎｍにおいて最も高い輝度を示した素子の輝度を１として規格化して示した。本実施例と比べると第１ピークと第２ピークの輝度の値の大小関係において実測との差が大きい。
【図面の簡単な説明】
【０１６６】
【図１】角度θの直接光と反射光を示した概略図。
【図２】Ｓ偏光の反射と透過を示した概略図。
【図３】Ｐ偏光の反射と透過を示した概略図。
【図４】角度θの直接光と反射光の光路差を示した概略図。
【図５】双極子遷移による角度θの直接光と反射光を示した概略図。
【図６】ある方向を向いた単独の発光分子からのS偏光の内部発光干渉モデル図。
【図７】ある方向を向いた単独の発光分子からのP偏光の内部発光干渉モデル図。
【図８】Ｓ偏光の内部発光干渉モデル図のうち発光層部分を示したモデル図。
【図９】Ｐ偏光の内部発光干渉モデル図のうち発光層部分を示したモデル図。
【図１０】遷移双極子モーメント、波数ベクトルおよび各電場の単位ベクトルの方向を示した概略図。
【図１１】基板から多層膜へＳ偏光の光が入射した時の電場を示した概略図。
【図１２】基板から多層膜へＰ偏光の光が入射した時の電場を示した概略図。
【図１３】基板および発光界面における光線の太さを示した概略図。
【図１４】同じスポット径が見込む発光分子数を示した概略図。
【図１５】発光分子が完全なミラーで挟まれた場合の干渉による強めあいを示した概略図。
【図１６】発光分子から基板法線方向に対して垂直に近い方向に放たれた光波が吸収される様子を示した概略図。
【図１７】Ａｌｑ３にクマリン５４５Ｔをドープした素子の基板法線方向の輝度の電子輸送層膜厚依存性を示す概略図。
【図１８】Ａｌｑ３にクマリン５４５Ｔをドープした薄膜の光励起スペクトルを示す概略図。
【符号の説明】
【０１６７】
１、１１、１７  光反射層
２、１２、１８  発光面
３、１９  反射光
４、２０  直接光
５、１５、２１  観察側
６  媒質０:複素屈折率Ｎ０
７  媒質１:複素屈折率Ｎ１
８  入射
９  反射
１０  透過
１３  反射光Ｂ
１４  直接光Ａ
１６  屈折率ｎ
２２  遷移双極子モーメント
２３  ｃｏｓθ０
２４  ｃｏｓθｈ
２５  発光分子の位置
２６  光波が進行した位置
２７  電場振幅Ｅ
２８  電場振幅５Ｅ
２９  光吸収
【図１】
2007095500.tif 000057
【図２】
2007095500.tif 000058
【図３】
2007095500.tif 000059
【図４】
2007095500.tif 000060
【図５】
2007095500.tif 000061
【図６】
2007095500.tif 000062
【図７】
2007095500.tif 000063
【図８】
2007095500.tif 000064
【図９】
2007095500.tif 000065
【図１０】
2007095500.tif 000066
【図１１】
2007095500.tif 000067
【図１２】
2007095500.tif 000068
【図１３】
2007095500.tif 000069
【図１４】
2007095500.tif 000070
【図１５】
2007095500.tif 000071
【図１６】
2007095500.tif 000072
【図１７】
2007095500.tif 000073
【図１８】
2007095500.tif 000074