(57)【要約】【課題】各種ダイオードの開発や性能向上のために、単一のI-V特性から様々な特性を有する太陽電池やダイオードに対応した全パラメータを安定して抽出できる新しい評価方法を提供する。【解決手段】ダイオードにおける電流Ｉと電圧Ｖとの関係を示したダイオード特性式のパラメータを評価する方法であって、仮定に基づいてダイオード特性式を近似した近似式を求めるステップと、ＩとＶの変化関係を測定した実験値に基づいて、前記近似式におけるパラメータを仮決定するステップと、仮決定したパラメータを前記ダイオード特性式に代入し、前記実験値と前記ダイオード特性式のずれを算出することにより、パラメータを修正する修正ステップと、を含むことを特徴とする評価方法。【選択図】図２
2009168559.tif 000002
【特許請求の範囲】
【請求項１】
ダイオードにおける電流Ｉと電圧Ｖとの関係を示したダイオード特性式のパラメータを評価する方法であって、
仮定に基づいてダイオード特性式を近似した近似式を求めるステップと、
ＩとＶの変化関係を測定した実験値に基づいて、前記近似式におけるパラメータを仮決定するステップと、
仮決定したパラメータを前記ダイオード特性式に代入し、前記実験値と前記ダイオード特性式のずれを算出することにより、パラメータを修正する修正ステップと、
を含むことを特徴とする評価方法。
【請求項２】
請求項１に記載の評価方法において、
前記仮決定するステップでは、候補パラメータを用いた近似式に対し前記実験値を適用して候補パラメータを修正する過程を繰り返すことにより、修正された候補パラメータを前記近似式におけるパラメータとして仮決定することを特徴とする評価方法。
【請求項３】
請求項１または２に記載の評価方法において、
前記仮決定するステップと、前記修正ステップとを交互に繰り返して、ダイオード特性式のパラメータの評価を行うことを特徴とする評価方法。
【請求項４】
請求項１に記載の評価方法において、
前記ダイオードの特性式は、単一ダイオードモデルの等価回路である評価方法。
【請求項５】
請求項４に記載の評価方法において、
前記パラメータには、並列抵抗Ｒｓｈが含まれ、
前記近似式を求めるステップでは、並列抵抗Ｒｓｈを無視することなく前記仮定を行う評価方法。
【請求項６】
請求項１に記載の評価方法において、
前記ダイオード特性式は、
【数１】
2009168559.tif 000020
（ここで、 Iph, I0, Rs, Rsh, q, ,n, kB, T は、それぞれ、光電流、ダイオードの飽和電流、直列抵抗、並列抵抗、電子素量、ダイオードの理想係数、ボルツマン定数及び温度）であり、
前記仮定をするステップにおいては、
【数２ａ】
2009168559.tif 000021
（ここでVoc, Iscは、それぞれ、解放電圧及び短絡電流である）の近似を行うことを特徴とする評価方法。
【請求項７】
請求項１に記載の評価方法において、
前記ダイオード特性式は、
【数１】
2009168559.tif 000022
（ここで、 Iph, I0, Rs, Rsh, q, ,n, kB, T は、それぞれ、光電流、ダイオードの飽和電流、直列抵抗、並列抵抗、電子素量、ダイオードの理想係数、ボルツマン定数及び温度）であり、
前記仮定をするステップにおいては、
【数２ｂ】
2009168559.tif 000023
（ここでVoc, Iscは、それぞれ、解放電圧及び短絡電流である）の仮定を行うことを特徴とする評価方法。
【請求項８】
請求項６または請求項７に記載の評価方法において、
前記近似式は、
【数９】
2009168559.tif 000024
（ここで、Ｒｓｏとn oは一時的なパラメータ）であり、
パラメータを仮決定をするステップにおいては、Ｒｓｏとn oの値がそれぞれＲｓとnの値に近いことを条件として、パラメータを仮決定することを特徴とする評価方法。
【請求項９】
ダイオードにおける電流Ｉと電圧Ｖとの関係を示したダイオード特性式のパラメータを評価する装置であって、
仮定に基づいてダイオード特性式を近似した近似式に対し、ＩとＶの変化関係を測定した実験値に基づいて、前記近似式におけるパラメータを仮決定する手段と、
仮決定したパラメータを前記ダイオード特性式に代入し、前記実験値と前記ダイオード特性式のずれを算出することにより、パラメータを修正する手段と、
を含むことを特徴とする評価装置。
【請求項１０】
コンピュータを、請求項９に記載の評価装置として機能させることを特徴とする評価プログラム。
【請求項１１】
請求項１に記載の評価方法を用いて、その他の太陽電池の解析を可能としたことを特徴とする太陽電池解析評価方法。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、各種ダイオード、太陽電池の特性を解析、評価するのに最適な解析評価方法およびその方法を搭載した解析評価装置に関するものである。
【背景技術】
【０００２】
太陽電池はクリーンな発電のために非常に有効なデバイスである為、その性能向上のための研究・開発が広くなされている。一般的に太陽電池の特性は、光電流Iph、ダイオードの飽和電流I0、直列抵抗Rs、並列抵抗Rsh、電子素量q、ダイオードの理想係数n、ボルツマン定数kB及び温度Tにより解析、評価される。したがって、その等価回路にも多くのパラメータを含んでいる。これまでのシリコン太陽電池の解析の場合、一般的には、並列抵抗Rsh は無視されてきた（非特許文献１，２）。
【０００３】
しかしながら、シリコン太陽電池でさえ並列抵抗Rshは無視できるほど大きくなく、それは、理想係数nや直列抵抗Rsの見積に大きな影響を与える（非特許文献３，４，５）。それゆえ、太陽電池の性能をさらに向上させるには、太陽電池の全パラメータを精確に見積もることが極めて重要である。 これまで、太陽電池の性能を解析するために多くの方法が提案されている（非特許文献２）。
【０００４】
しかしながら、それらの多くは、並列抵抗Rshを無視するものであったり、特別な仮定を必要としたりものであった。 例えば、直列抵抗Rs は2つの異なる照射量で測定されたI-V特性から近似的に見積もることができることが知られており（非特許文献１，６）、この手法は、“different illumination level method”と呼ばれている。しかしながら、照射量を変化させることによる温度変化が大きな誤差を生むことが指摘されている（非特許文献２）。したがって、単一のI-V特性から全パラメータを抽出することが重要である。
【０００５】
このような手法は、“constant illumination level method”と呼ばれている。constant illumination level methodの1つとして、(Isc-I)-1を横軸とした-(dV/dI)プロットのy切片及び傾きから直列抵抗Rsと理想定数n を見積もる手法が知られている。ここで、Iscは短絡電流である。この作図を用いる方法はとてもシンプルであるが、これを用いるためには、Rshが無視できるほど十分大きくなければならない。したがって、この条件が成立しない場合には、見積もられた理想係数nと直列抵抗Rsには大きな誤差を含むことになる。
【０００６】
その他の手法として、カーブフィッティングのような計算機を用いる方法も提案されている（非特許文献３，４，７，８，９）。しかしながら、カーブフィッティングなどに広く用いられている最小2乗法を太陽電池の解析に直接用いることはできない（非特許文献３，４，７，８，９）。また、その結果は、初期に強く依存することが知られている。初期値の依存性や不安定性はフィッティング法が確立していない大きな原因の1つである（非特許文献４，７，８）。
【０００７】
【非特許文献１】Martin Wolf and Hans Rauschenbach, Adv. Energy Convers. 3, 455 (1963).
【非特許文献２】P. Mialhe, A. Khoury, and J. P. Charles, phys. stat. sol. (a) 83, 403 (1984).
【非特許文献３】T. Easwarakhanthan, J. Bottin, I. Bouhouch and C. Boutrit, Int. J. Solar Energy 4, 1 (1986).
【非特許文献４】M. Chegaar, Z. Ouennouchi, and A. Hoffmann, Solid-State Electron. 45, 293 (2001).
【非特許文献５】V. N. Singh and R. P. Singh, J. Phys. D 16, 1823 (1983).
【非特許文献６】R. J. Handy, Solid-State Electron. 10, 765 (1967).
【非特許文献７】A. Kaminski, J.J. Marchand, A. Laugier, Sold-State Eltectron. 43, 741 (1999).
【非特許文献８】Masaki Murayama and Tatsuo Mori, Jpn. J. Appl. Phys. 45, 542 (2006).
【非特許文献９】Masaki Murayama and Tatsuo Mori, Thin Solid Films 509, 123 (2006).
【非特許文献１０】A. Yakimov and S. R. Forrest, Appl. Phys. Lett. 80, 1667 (2002).
【発明の開示】
【発明が解決しようとする課題】
【０００８】
有機太陽電池や色素増感太陽電池（DSSC）などの新しい太陽電池、または、有機ＥＬなどの各種ダイオードの開発や性能向上のためには、単一のI-V特性から様々な特性を有する太陽電池やダイオードに対応した全パラメータを安定して抽出できる新しい方法の開発が極めて重要である。そこで、本発明では、constant illumination level methodとフィッティング法を改良した適用可能範囲が広く、安定した方法を提供する。
【課題を解決するための手段】
【０００９】
一般的に太陽電池やフォトダイオードの特性は、図1に示す単一ダイオードモデルの等価回路を用いて説明され、電流I、電圧Vとの関係は、
【数１】
2009168559.tif 000003
によって記述される。ここで、 Iph, I0, Rs, Rsh, q, ,n, kB, Tは、それぞれ、光電流、ダイオードの飽和電流、直列抵抗、並列抵抗、電子素量、ダイオードの理想係数、ボルツマン定数及び温度である。また、ダイオードの特性は、逆方向電流を正として、式(1)において、Iph=0とすることにより表される。等価回路は多くのパラメータを含んでいることがわかる。
【００１０】
本発明による全パラメータを抽出するための方法のフローチャートを図2に示す。本方法での仮定は、次に示すたった1つの仮定、
【数２ａ】
2009168559.tif 000004
あるいは、
【数２ｂ】
2009168559.tif 000005
のみである。ここで、Voc, Iscは、それぞれ、解放電圧及び短絡電流である。仮定(2a)及び(2b)は通常の太陽電池やダイオードに対して一般的に成立する。
【００１１】
式(1)から、dV/dIは、
【数３】
2009168559.tif 000006
のように与えられる。
【００１２】
ここで仮定(2a)及び、式(1)から
【数４】
2009168559.tif 000007
が得られ、それを式(3)に代入して、
【数５】
2009168559.tif 000008
を得る。
【００１３】
また、ダイオード(Iph = 0)の場合においても、式(1)及び仮定(2b)から同様に式(5)が得られる。ただし、Isc = Isc = 0である。I = Isc (V = 0)の場合、式(5)は、
【数６】
2009168559.tif 000009
のように近似される。ここで用いた近似Rsh>>Rsは一般的に成立する。式(6)は、RshはI-V特性のI =Isc(V=0)の点における傾きから単独に求めることができることを意味している。
【００１４】
一方、理想係数nと直列抵抗Rsを決定するために、一時的なパラメータRs0、n0を式(5)に導入する。それにより、式(5)は、
【数７】
2009168559.tif 000010
のように変形される。
【００１５】
ここで|δ|<<1は、
【数８】
2009168559.tif 000011
と定義する。|δ|<<1のとき、式(7)は、
【数９】
2009168559.tif 000012
のように与えられる。
【００１６】
条件|δ|<<1は、Rs0とn0の値が、それぞれ、直列抵抗Rs及び理想係数nの値に近いとき成立する。このとき、式(9)右辺の分母から直列抵抗Rs及び理想係数n は消えており、直列抵抗Rs及び理想係数nは、(Isc-I-((V-Rs0(Isc-I)-n0kBT/q)/Rsh))-1を横軸とした-(dV/dI)プロットy切片及び傾きから、それぞれ求めることができる。最適な一時パラメータを見積もるために、図2に示すとおり、それぞれの一時パラメータの初期値を0として、式(9)を用いた繰り返し計算を行う。
【００１７】
すなわち、初回のループで見積もられた直列抵抗Rs及び理想係数nには、一時パラメータの値として0を用いているために大きな誤差を含んでいる。しかしながら、それらの値は初期値である0よりも実際の値に近づいている。そこで、一時パラメータRs0及びn0 の候補として、これらの値をそれぞれ用いる。これにより、次のループで見積もられる直列抵抗Rs及び理想係数nは初回のループよりも実際の値に近い値を示すことになる。
【００１８】
これを図2に示すようにRs=Rs0,n=n0、すなわち|δ|=0になるまで繰り返すことにより、直列抵抗Rs及び理想係数nを決定することができる。ここで、式(1)におけるV =Voc、I =0及びV=0、I =Iscの境界条件と仮定（2a）、(2b)から
【数１０】
2009168559.tif 000013
式(10)が得られる。
【００１９】
得られた式(10)と式(4b)より
【数１１】
2009168559.tif 000014
及び
【数１２】
2009168559.tif 000015
がそれぞれ得られる。
【００２０】
すでに求まっている直列抵抗Rs及び理想係数n、並列抵抗Rshの値を式(11)に代入することにより、ダイオードの飽和電流を求めることができる。さらに、これまでに求まっているパラメータと飽和電流I0の値を式(12)へ代入することにより、光電流Iphが決まる。
【００２１】
一方、ダイオード(Iph=0)の場合は、Isc=Iph =0であり、飽和電流I0は、式(1)から得られる次式
【数１３】
2009168559.tif 000016
により求まる。
【００２２】
この段階で、すべてのパラメータの候補が出揃い、I-V特性を式(1)を用いて計算することができる。最適なパラメータの値は、図2に示すように、実験値とその計算値との平均2乗誤差(MSE)が最小(MMSE)となるようにパラメータの見積を最初の並列抵抗Rshと短絡電流のIscの見積から繰り返すことにより得られる。最後に、仮定(2a)または(2b)を確認することにより、見積もられたパラメータの妥当性をチェックする。
【発明の効果】
【００２３】
太陽電池の特性を解析するためには、その太陽電池の全パラメータを精確に見積もることが重要である。そこで、我々は、それら全パラメータを安定して決定するための方法を提供した。本方法で必要とする仮定は、太陽電池であれば通常成立する条件1つのみである。本発明では、すでに報告されているシリコン太陽電池や有機太陽電池、色素増感太陽電池（DSSC）のI-V特性を本方法に適用した。その結果、本方法で決定されたパラメータは以前の報告よりもそれらのI-V特性を忠実に表現していた。つまり、我々の提供する方法は、様々な特性をもつ太陽電池に適用可能かつ安定な方法であり、新しい太陽電池の性能を向上させるために非常に有効な方法である。
【発明を実施するための最良の形態】
【００２４】
本方法を、様々な太陽電池のI-V 特性の解析に適用した。用いた太陽電池は、シリコン(Si)、3つのヘテロ接合を有する有機太陽電池、及び色素増感太陽電池（DSSC）の3種類である。
【００２５】
本方式では、離散的な実験データからdV/dIの値を見積もる必要がある。そこで、本発明では、
【数１４】
2009168559.tif 000017
を用いてdV/dIとした。
【００２６】
ここで、Vm及びImはI-V 特性のデータ中におけるm番目の電圧、電流の値である。また、Mは、データ点の総数である。電流Iの値が短絡電流Iscのそれに近いとき、たとえIの値の誤差が小さいとしても、(Isc-I)の値に対するIの誤差が大きくなるため、式(9)右辺の計算値には大きな誤差を含むことになる。
【００２７】
そこで、式(9)右辺の計算には、次に示す条件を満たすデータ点のみを用いる。
【数１５】
2009168559.tif 000018
実験値からその微分した値を見積もるためには、滑らかなデータが必要である。そこで、滑らかでないデータについては、条件(15)を満たすデータを多項式によって近似し、その近似多項式から式(9)左辺dV/dIの値を見積もった。本発明では、色素増感太陽電池（DSSC）の解析において9次の多項式を用いる。
【実施例１】
【００２８】
Si太陽電池の解析に本方法を適用した。図3 (a)に、実験によるI-V特性(○)と本方法によって計算されたI-V特性(実線)を示す。両者のI-V特性はよく一致していることがわかる。表１は、本方法及び以前の方法によって見積もられた太陽電池のパラメータの値とΔ1の値を示す。表1に本方法によって見積もられたパラメータと、すでに報告されている他の方法による見積及びΔ1の値を示す。本方法によって得られたパラメータの値は他の方法のそれとよく一致していることがわかる。また、Δ1は１に比べて十分小さく、仮定(2a)を十分満足していることがわかる。
【表１】
2009168559.tif 000019
【００２９】
これらの結果は、本方法の妥当性を支持している。図3(b)にパラメータが最適化されているとき(Rsh=32Ω)の(Isc-I-((V-Rs(Isc-I)-nkBT/q)/Rsh))-1を横軸としたときの-(dV/dI)及び(Isc-I-V/Rsh)/Iscのプロットを示す。本方法を用いた場合の-(dV/dI)プロットは、よい直線性を示しており、Rshを無視した場合には非線形な特性となっていることがわかる。
【００３０】
これは、Si太陽電池でさえRshを無視できないことを示している。式(15)の下、理想係数nと直列抵抗Rsを見積もるとそれぞれ1.4及び37Ωとなった。図4に(a)理想係数nと直列抵抗Rsの収束特性と(b)MSEの並列抵抗Rsh依存性を示す。図4(a)に示すとおり、理想係数nと直列抵抗Rsは初期値として0を用いているにも関わらず、数回の反復計算のみで急速に収束していくことがわかる。また、図4(b)に示すように、MSEはRshの変化によって急速にMMSEに近づいていくことがわかる。これらの結果は、本方法が高速な収束特性を有していることを示している。
【００３１】
Si太陽電池モジュールについても本方法を用いて解析を行った。本方法によって得られたパラメータの値は、表1に示すように、非特許文献3のものとよく一致し、以前に報告されているもう1つの方法（非特許文献４）では、本方法と比較して小さな並列抵抗Rshの値及び大きな理想係数nが見積もられている。以前の方法による見積の失敗は、初期値に強く依存する一般的なフィッティング法によるパラメータ抽出法の不安定性、困難さを表している。
【００３２】
図5にさらに別のSi太陽電池のI-V特性の実験データ(○)と本方法(実線)及び他の方法(破線)（非特許文献５）によって計算されたI-V特性を示す。破線は全くI-V特性の実験データを反映していないことがわかる。非特許文献５は、理想係数nは2.5かつ並列抵抗Rshは無視できるほど十分大きいとの仮定の下、I-V特性の実験データ中で最大電力となる1点のみを用いて見積り、直列抵抗Rsが3.3Ωと見積もっている。その結果、それらの値から計算されたI-V特性は実験結果から逸脱する結果となった。一方、本方法(実線)は、完全に実験データと一致している。これは、本方法によって見積もられたパラメータは実際の太陽電池の特性をよく反映しており、理想係数nを固定したり、並列抵抗Rshを無視したりすることが大きな誤差を生じることを示している。
【実施例２】
【００３３】
本方法は、有機太陽電池の解析にも適用可能である。図6に有機太陽電池のI-V特性の実験結果(実線)（非特許文献１０）及び本方法によって見積もられたパラメータから計算されたI-V特性(破線)を示す。計算結果と実験結果は、見分けのつかないほど完璧に一致していることがわかる。また、表1に示すように仮定(2a)も十分満足している。実験結果と計算結果の完全な一致は、本方法がシリコン太陽電池だけでなく様々な特性持つ他の種類の太陽電池についても適用可能であることを示している。
【実施例３】
【００３４】
図7(a)に非特許文献８で”C4”と呼ばれるDSSCのI-V特性の実験結果(○)及び以前の方法(破線)及び本方法(実線)によって得られたパラメータから計算されたI-V特性を示す。この場合、図7(b)に示すようにdV/dIを直接実験データから見積もるにはその滑らかさが不十分である。本発明では、条件(15)を満たす実験データを9次の多項式によって近似し、その近似多項式からdV/dIを見積もった。図7(a)の挿入図は、条件(15)を満たす実験データ(○)及び9次多項式による近似曲線(実線)の拡大図である。近似曲線は、実験データをよく反映していることがわかる。図7(b)に示すとおり、実験データを直接用いた場合、-(dV/dI)プロットは散らばっており、直線を引くことが不可能であるが、多項式近似を用いた場合のそれは、良い直線性を有しており、容易に理想係数n及び直列抵抗Rsを見積もることが可能であることがわかる。図7(a)における実線は、多項式近似を導入して本方法から得られた計算結果であるが、破線よりもよく一致していることがわかる。破線は、非特許文献８で理想係数nを2.5に固定してフィッティングした結果である。先述したとおり、不明なパラメータを固定しなければならない方法は決して精確なパラメータの値を抽出することはできない。
【００３５】
なお、参考として、本願における発明の別の態様について以下に記す。
【００３６】
１．ダイオードの特性を評価する方法において、並列抵抗Ｒｓｈを無視することなく、解析に必要な全パラメータを一意的に安定かつ精確に取得し、実データに極めて近似させることを可能としたダイオード特性解析評価方法。
【００３７】
２．太陽電池の特性を、単一ダイオードモデルの等価回路を用いて解析し、評価する方法において、並列抵抗Ｒｓｈを無視することなく、解析に必要な全パラメータを一意的に安定かつ精確に取得し、実データに極めて近似させることを可能とした太陽電池解析評価方法。
【００３８】
３．Δ（ある演算式）が極めて０に近いことのみを条件として、シリコン系太陽電池のみならず、その他の太陽電池の解析をも可能とした方法。
【００３９】
４．Δ（ある演算式）が極めて０に近いことを条件として、ＲｓｈをＩ＝Ｉｓｃ（Ｖ＝０）の点における傾きから単独に求める方法。
【００４０】
５．一時パラメータのＲｓｏと、noの値がそれぞれ直列抵抗Ｒｓ、及び理想係数nの値に近いことを条件として、直列抵抗Ｒｓ、及び理想係数nを、−ｄＶ/ｄＩプロットのｙ切片及び傾きからそれぞれ求める方法。
【００４１】
６．実データをｎ次多項式（ｎ＝１，２，・・・）により近似し、その近似式から、ｄＶ/ｄＩを得る手法を用いた方法。
【００４２】
７．直列抵抗Ｒｓ及び理想係数nが、Ｒｓ＝Ｒｓｏ、n＝no、δ（ある演算式）（絶対値）＝０になるまで繰り返し演算することにより、直列抵抗Ｒｓ及び理想係数nを決定する方法。
【００４３】
８．平均二乗誤差（ＭＳＥ）が最小（ＭＭＳＥ）となるようにパラメータの見積もりを、最初の並列抵抗Ｒｓｈと短絡電流Ｉｓｃの見積もりから繰り返し演算を行うことによって求める方法。
【００４４】
９．ダイオードの特性を解析し評価する解析評価装置において、請求項１から８記載のアルゴリズムをプログラム化して搭載し、ダイオードのパラメータの値を決定することを可能ならしめた装置。
【００４５】
１０．太陽電池の特性を解析し評価する解析評価装置において、請求項１から７記載のアルゴリズムをプログラム化して搭載し、太陽電池のパラメータの値を決定することを可能ならしめた装置。
【産業上の利用可能性】
【００４６】
本発明は、単一のI-V特性から太陽電池やダイオードの全パラメータを安定して決定するための方法を提供する。本方法を用いるためには、通常の太陽電池に対して成立する仮定を1つしか必要としない。本方法を用いることにより、以前の方法ではパラメータの値を抽出できないような特性を持つ様々な太陽電池でさえ、精確なパラメータの値を決定することができる。また、本方法は、高速な収束特性を有している。適用範囲が広く、安定して動作する本方法を搭載した太陽電池解析評価装置は、様々な新しい太陽電池、ダイオードの解析、評価に役立ち、研究・開発を促進するとともに、太陽電池、有機ELなどのダイオードの特性向上にもつながるものである。
【図面の簡単な説明】
【００４７】
【図１】単一ダイオードモデルによる太陽電池の等価回路。
【図２】本方法のフローチャート。
【図３ａ】Si太陽電池(R.T.C France)のI-V特性の実験結果(○)（非特許文献３）及び、本方法によって見積もられたパラメータを用いて計算されたI-V特性(実線)。
【図３ｂ】最適化された(Rso=Rs,n0= n, and Δ=0)後の[Isc-I-{V-Rs(Isc-I)-nkBT/q}/Rsh]-1を横軸としたときの-dV/dIと(Isc-I-V)/Isc プロット(Rsh = 32Ω)。
【図４ａ】理想係数nと直列抵抗Rsを見積もるループの収束特性とMSEの並列抵抗Rsh依存性。
【図４ｂ】理想係数nと直列抵抗Rsを見積もるループの収束特性とMSEの並列抵抗Rsh依存性。
【図５】Si太陽電池モジュールのI-V特性の実験結果(○)（非特許文献５）と、本方法(実線)及び以前の方法(破線)（非特許文献５）によって見積もられたパラメータから計算されたI-V特性。
【図６】有機太陽電池のI-V特性の実験結果(実線)（非特許文献１０）と本方法によって見積もられたパラメータから計算されたI-V特性(破線)。
【図７ａ】DSSC(C4)のI-V特性の実験結果(○)（非特許文献８）と本方法(実線)及び以前の方法(破線)（非特許文献８）によって得られたパラメータから計算されたI-V特性。挿入図は、I-V特性の実験結果(○)と多項式による近似曲線(実線)。
【図７ｂ】最適化された後の[Isc-I-{V-Rs(Isc-I)-nkBT/q}/Rsh]-1を横軸としたときの-dV/dI及び(Isc-I-V)/Iscプロット。
【図１】
2009168559.tif 000025
【図２】
2009168559.tif 000026
【図３ａ】
2009168559.tif 000027
【図３ｂ】
2009168559.tif 000028
【図４ａ】
2009168559.tif 000029
【図４ｂ】
2009168559.tif 000030
【図５】
2009168559.tif 000031
【図６】
2009168559.tif 000032
【図７ａ】
2009168559.tif 000033
【図７ｂ】
2009168559.tif 000034