JP2013164331A - 電極の評価装置および評価方法 - Google Patents

Abstract

【課題】正確なインピーダンス評価を行うことができる電極の評価装置等を提供する。
【解決手段】インピーダンス取得手段は、特定の直流動作条件の下で前記電極のインピーダンス特性を取得する。電流値取得手段は、前記特定の直流動作条件の下における直流電流値の経時的変化を取得する。規格化インピーダンス算出手段は、前記電流値取得手段で取得された前記直流電流値の経時的変化を適用して、前記インピーダンス取得手段により取得されたインピーダンスに前記直流電流値を乗じた規格化インピーダンスを算出する。
【選択図】図１

Description

本発明は、電気化学特性に基づいて電極の特性を評価する電極の評価装置等に関し、とくに正確なインピーダンスの評価が可能な評価装置等に関する。

燃料電池等における電極反応は、一般的に、（１）反応物質の電極表面への物質輸送（拡散）、（２）電極表面上での触媒による電荷移動反応、の過程を経て進行する。電極反応においてどの過程がどの程度支配的であるのかを知ることは電極の特性や性能劣化の要因などを評価する上で非常に重要である。そのための手法として、直流（ＤＣ）特性およびインピーダンスによる評価がある。

特開２０１２−３９２２号公報 特開２００９−４８８１４号公報 特開２００９−４８８１３号公報

燃料電池の一般的な触媒であるＰｔ（白金）は、酸性電解液中において安定で酸素還元活性の高い純物質とされているが、電位履歴により表面酸化物の量や質が変化することや、表面の酸化物生成や有機物などの表面吸着によって触媒活性が大きな影響を受けることが一般的に知られている。そのため、電位履歴をもった状態から測定を開始すると、Ｐｔの表面状態は履歴によってバラバラであるため、直流特性もインピーダンス特性も再現性のないものとなってしまう。また、電位によっては酸化物生成やアニオン・有機物などの吸着が経時的に進行するため、触媒活性（直流電流値）は刻々と変化してしまう。そのため、触媒活性の評価結果は測定時の条件等によって大きく異なり、その難しさ故に国際的な評価基準は未だ明確に定まっていない。このような理由から、インピーダンスの評価の際に、電位履歴をキャンセルすることや、測定中の経時的な触媒活性変化まで加味することが望ましいと考えられる。

本発明の目的は、正確なインピーダンス評価を行うことができる電極の評価装置等を提供することにある。

本発明の電極の評価装置は、電気化学特性に基づいて電極の特性を評価する電極の評価装置において、電極に与える電位制御を行う電位制御手段と、前記電位制御手段により与えられる特定の直流動作条件の下で前記電極のインピーダンス特性を取得するインピーダンス取得手段と、前記特定の直流動作条件の下における直流電流値の経時的変化を取得する電流値取得手段と、前記電流値取得手段で取得された前記直流電流値の経時的変化を適用して、前記インピーダンス取得手段により取得されたインピーダンスに前記直流電流値を乗じた規格化インピーダンスを算出する規格化インピーダンス算出手段と、を備えることを特徴とする。
この電極の評価装置によれば、電流値取得手段で取得された直流電流値の経時的変化を適用して、インピーダンス取得手段により取得されたインピーダンスに直流電流値を乗じた規格化インピーダンスを算出するので、正確な規格化インピーダンスを求めることができる。

前記電流値取得手段は、前記特定の直流動作条件の下におけるクロノ・アンペロメトリー、または前記インピーダンス取得手段により前記インピーダンス特性が取得される間における、前記電極に接続されたポテンションスタットからの制御された直流電流に基づいて前記直流電流値の経時的変化を取得してもよい。

前記電位制御手段は、前記特定の直流動作条件への移行に先行して、高速掃引を繰り返す工程、前記特定の直流動作条件まで電位制御する工程、およびバックグラウンド電流を開放するために電位を維持して待機する工程を順次実行してもよい。

高速掃引直後におけるリニア・スイープ・ボルタンメトリーの測定を複数の掃引速度について行うことにより、前記特定の直流動作条件に至るまでの掃引速度を決定する掃引速度決定手段を備え、前記電位制御手段は、前記掃引速度決定手段により決定された前記掃引速度で前記特定の直流動作条件まで電位制御してもよい。

前記掃引速度決定手段は、前記複数の掃引速度のうち定常状態が得られる最も遅い掃引速度を前記特定の直流動作条件に至るまでの掃引速度として決定してもよい。

前記特定の直流動作条件の前後にサイクリック・ボルタンメトリーを実施し前後の波形を比較することで、前記特定の直流動作条件における表面吸着についての評価を行う手段を備えてもよい。

本発明の電極の評価方法は、電気化学特性に基づいて電極の特性を評価する電極の評価方法において、電極に与える電位制御を行う電位制御ステップと、前記電位制御ステップにより与えられる特定の直流動作条件の下で前記電極のインピーダンス特性を取得するインピーダンス取得ステップと、前記特定の直流動作条件の下における直流電流値の経時的変化を取得する電流値取得ステップと、前記電流値取得ステップで取得された前記直流電流値の経時的変化を適用して、前記インピーダンス取得手段により取得されたインピーダンスに前記直流電流値を乗じた規格化インピーダンスを算出する規格化インピーダンス算出ステップと、を備えることを特徴とする。
この電極の評価方法によれば、電流値取得ステップで取得された直流電流値の経時的変化を適用して、インピーダンス取得ステップにより取得されたインピーダンスに直流電流値を乗じた規格化インピーダンスを算出するので、正確な規格化インピーダンスを求めることができる。

本発明の電極の評価装置によれば、電流値取得手段で取得された直流電流値の経時的変化を適用して、インピーダンス取得手段により取得されたインピーダンスに直流電流値を乗じた規格化インピーダンスを算出するので、正確な規格化インピーダンスを求めることができる。

また、本発明の電極の評価方法によれば、電流値取得ステップで取得された直流電流値の経時的変化を適用して、インピーダンス取得ステップにより取得されたインピーダンスに直流電流値を乗じた規格化インピーダンスを算出するので、正確な規格化インピーダンスを求めることができる。

電極の評価方法を機能的に示すブロック図。 制御演算装置の動作を示すフローチャート。 電気化学的前処理（ステップＳ１）およびＬＳＶ特性測定時（ステップＳ２）における電位制御を示す図。 酸化還元反応におけるＬＳＶの掃引速度依存性を示す図であり、（ａ）は縦軸を電流値、（ｂ）は縦軸を電流値の対数とした図。 定常状態になった際のＬＳＶ特性から得られた「電位−ターフェル勾配」特性を示す図。 インピーダンス測定プロトコルを示す図。 インピーダンス特性を示す図であり、（ａ）は、各測定電位（Ｖ１，Ｖ２，・・・Ｖｘ）ごとのインピーダンス特性（コール−コールプロット）を示す図、（ｂ）は（ａ）に示すインピーダンスに対して、周波数ごとにＤＣ電流値ｉｆを乗じることで得られる、触媒反応に由来する規格化されたインピーダンス（ｉｆ×Ｚ）特性を示す図。 工程「Ｅ」における測定電位Ｖ１，Ｖ２，Ｖ３におけるＤＣ電流値の経時的変化を示す図。 燃料電池セルの測定システムを示すブロック図。 反応抵抗の抽出およびインピーダンスの規格化処理を示すフローチャート。 ｂ１、ｂ２領域における測定に基づき取得されたインピーダンスを示す図。 直流補正を施した後のインピーダンスを示す図。 図４に示すインピーダンス（直流抵抗補正済）の|Z|-f特性、θ-f特性を示す図。 モデルフィッティング処理の方法を示す図。 分布定数域分を考慮し原点補正したインピーダンスを示す図。 抽出された反応抵抗に対して規格化処理して得たインピーダンス特性を示す図。 測定システムの動作を示すフローチャート。 拡散領域まで測定したターフェルプロットを示す図。 ターフェルプロットを示す図であり、（ａ）は電流−ターフェル勾配特性を、（ｂ）は電流の対数−ターフェル勾配特性を、それぞれ示す図。 ｂ２の領域から拡散の影響がでている場合の電流−ターフェル特性を示す図。 ｂ１の領域から拡散の影響がでている場合の電流−ターフェル特性を示す図。 拡散補正を施した電流−ターフェル勾配特性を示す図。 各動作点で測定した規格化インピーダンスを示す図。 算出された拡散性能評価パラメータを例示する図。 拡散の補正を施した規格化インピーダンスを示す図。 図６に示すインピーダンス測定プロトコルにＣＶを実施する工程を挿入した例を示す図。

以下、本発明による電極の評価方法の実施形態について説明する。

図１は、本実施形態の電極の評価方法を機能的に示すブロック図である。

図１に示すように、本実施形態の電極の評価方法は、電気化学測定の対象となる３電極系電解セル１００、３電極系電解セル１００に接続されるポテンションスタット１、ポテンションスタット１による測定データを記録するレコーダ２、および評価に必要な制御や演算を実行する制御演算装置３を用いて行われる。

３電極系電解セル１００は、作用極１０１、対極１０２および参照極１０３を備え、作用極１０１上ではバトラー・ボルマー式に基づく反応が起きているものとする。参照極１０３はルギン管１０４を介して作用極１０１の近傍に接続されている。作用極１０１、対極１０２および参照極１０３は、それぞれ、ポテンションスタット１に接続される。

なお、図１では３電極系電解セルを測定対象電極とする例を示しているが、本発明における測定対象は電解セルに限定されず、例えば電池アッセンブリなどでもよい。電気化学測定は原則として３電極系で行うが、対極のインピーダンスが小さく、安定な平衡電位を持つ場合には参照極を対極と兼用し、２電極系での測定としてもよい。

図１に示すように、制御演算装置３は、ポテンションスタット１による電位規制を指定の条件に従って制御する負荷制御部３１と、サイクリック・ボルタンメトリー（ＣＶ：Cyclic Voltammetry）、リニア・スイープ・ボルタンメトリー（ＬＳＶ：Linear Sweep Voltammetry）、クロノ・アンペロメトリー（ＣＡ：Chrono Amperometry）などの直流（ＤＣ）特性を取得する直流特性取得部３２と、直流特性取得部３２により得られたＤＣ特性を解析する直流特性解析部３３と、インピーダンス特性（ＥＩＳ：Electrochemical Impedance Spectroscopy）を取得するインピーダンス取得部３４と、インピーダンス取得部３４により得られた結果に各種補正を含む一連の規格化処理を施す規格化インピーダンス算出部３５と、直流特性解析部３３による解析結果および規格化インピーダンス算出部３５により得られる周波数特性値をモニタ等に提示する提示部３６と、を備える。なお、レコーダ２において、３電極系電解セル１００から得られるデータの他、直流特性解析部３３による解析結果および規格化インピーダンス算出部３５により得られる周波数特性値を保存できるようにしてもよい。

次に、制御演算装置３の動作を説明する。

図２は、制御演算装置３の動作を示すフローチャートである。

図２のステップＳ１では、電極（作用極１０１）の電気化学的前処理を実施し、続いてステップＳ２ではＬＳＶ特性を取得する。これらのステップＳ１およびステップＳ２は、ステップＳ２におけるＬＳＶの掃引速度を変更しながら複数回繰り返される。

図３は、電気化学的前処理（ステップＳ１）およびＬＳＶ特性測定時（ステップＳ２）における電位制御を示す図である。横軸が時間、縦軸が電位（対参照極）を示している。図３では、電気化学的前処理（ステップＳ１）およびＬＳＶ特性測定時（ステップＳ２）を、１回目、２回目、３回目・・・と繰り返す様子が表現されている。また、図３の例では、１回目、２回目、３回目・・・と測定が繰り返される度に、掃引速度が徐々に遅くなっている。

電気化学的前処理では、電極表面の不純物除去や電極の電位履歴のキャンセルを目的として定常状態となるまで、すなわち電流の応答状態が一定となるまで高速電位掃引をサイクリックに実施し、続いてＬＳＶ特性を測定する。図３に示された上限電位（ＶＨ）や下限電位（ＶＬ）および掃引速度・方向・回数は電極劣化を防ぐなどの条件に応じて変更することができる。掃引速度を変えて再びＬＳＶ測定を行う際には、電極の表面状態の再現性を得るために毎回、サイクリックに高速電位掃引を行う。

図４は、酸化還元反応におけるＬＳＶの掃引速度依存性を示す図であり、図４（ａ）は縦軸を電流値、図４（ｂ）は縦軸を電流値の対数としたものである。横軸はいずれも電位（対参照極）である。

図４に示すように、触媒活性は掃引速度が大きくなるにつれて高くなり、やがて定常値に収束する。これは掃引速度が大きくなるにつれ、電極反応の妨げとなる不純物等の表面吸着が追いつかなくなり、その影響が小さくなっているからだと考えられる。すなわち、不純物等の表面吸着の影響を受けずに目的の測定電位まで到達させるには、上記の定常状態となる掃引速度以上で電位掃引する必要がある。一方、掃引速度が大き過ぎると電気二重層の充放電電流などのバックグラウンド電流の影響が大きくなり、これが誤差要因となる。そのため、目的のインピーダンス測定電位まで到達させるために最適な掃引速度の決定が重要となる。

なお、より厳密な分析のため、電気二重層の充放電電流など不活性雰囲気下において測定されたバックグラウンド電流を用いた補正や、インピーダンスから得られる直流抵抗値を用いたＩＲ補正（電流値と直流抵抗値の積による電圧補正）をあらかじめ施して評価してもよい。

また、ここで定常状態になった際のＬＳＶ特性を用いて触媒活性評価（ＤＣにおける触媒活性評価）を行ってもよい。

図５は上記の定常状態になった際のＬＳＶ特性から得られた「電位−ターフェル勾配」特性を示す図である。この特性から、電荷移動律速領域と拡散律速移動領域をおおまかに分離でき、目的に応じた測定電位領域を判断することができる。

次に、ステップＳ３では、インピーダンス測定プロトコルを決定する。ここでは、ステップＳ１〜ステップＳ２の繰り返しで得られた掃引速度の異なるＬＳＶ特性に基づいて、例えば、上記の定常状態が得られたものの中で最も遅い掃引速度を「目的の測定電位まで到達させるための掃引速度」として決定する。また、図５の「電位−ターフェル勾配」特性に基づいて、目的に応じてインピーダンスを測定する「電位レンジ」と「測定電位」（電位レンジ内における測定電位Ｖ１，Ｖ２，Ｖ３等）を決定する。

次に、ステップＳ４では、電極の電気化学的前処理を実施し、続いてステップＳ５ではＣＡ特性を取得する。さらに、ステップＳ６では、再び電極の電気化学的前処理を実施し、続いてステップＳ７ではインピーダンス特性を取得する。これらのステップＳ４〜ステップＳ７は、決定されたインピーダンス測定プロトコル（ステップＳ３）に従って、測定電位を切り替えつつすべての測定電位Ｖ１，Ｖ２，・・・Ｖｘについて繰り返される。

図６は、ステップＳ４〜ステップＳ７を繰り返す手順に対応するインピーダンス測定プロトコルを示す図である。図６では、インピーダンス測定プロトコルの各工程をＡ，Ｂ，Ｃ，Ｄ，Ｅ，Ｆ，Ｇ，Ｈに区分して示している。

まず、工程「Ａ」では、電極表面に付着している不純物等を取り除くため、できるだけ広い電位レンジ（上限電位Ｖ_H−下限電位Ｖ_L間）で定常状態となるまで高速電位掃引を繰り返し、電極表面を清浄にする。

次に、工程「Ｂ」では、開始電位Ｖｂ−下限電位Ｖ_L間で高速電位掃引を数回繰り返し、電極の電位履歴をキャンセルする。

次に、工程「Ｃ」では、上記のＬＳＶ特性に基づきステップＳ３において決定された「目的の測定電位まで到達させるための掃引速度」において、開始電位Ｖｂから最初の測定電位Ｖ１まで掃引する。

次に、工程「Ｄ」では、測定電位Ｖ１において数秒待機し、工程「Ｃ」の掃引に起因する電気二重層の充放電電流、表面酸化物の生成電流や還元電流などのバックグラウンド電流を開放する。

次に、工程「Ｅ」では、インピーダンスの規格化処理に用いるためのＤＣ電流値の経時変化をＣＡにより測定、記録する。

次に、工程「Ｆ」では、インピーダンス測定終了後直後のＤＣ電流値を確認するためＣＡにより電流値を測定、記録する。

次に、工程「Ｇ」では、電位を制御しながらゆっくりと開始電位Ｖｂ（図６）まで戻す。

次に、再び工程「Ｂ」に戻り、工程「Ｅ」を工程「Ｈ」に置換した手順（工程「Ｂ」→「Ｃ」→「Ｄ」→「Ｈ」→「Ｆ」→「Ｇ」）を繰り返す。工程「Ｈ」では、工程「Ｅ」のＣＡに替えてインピーダンス測定（ＥＩＳ）を実行する。ここでは、各測定電位での測定点において交流電圧または交流電流を重畳し、その応答電流または応答電圧に基づいてインピーダンスを算出する。

以上で、測定電位Ｖ１における一連の手順は完了し、測定電位Ｖ２，Ｖ３，・・・Ｖｘについて、ＤＣ電流値の計測のための手順（工程「Ｂ」→「Ｃ」→「Ｄ」→「Ｅ」→「Ｆ」→「Ｇ」）と、インピーダンス測定（ＥＩＳ）のための手順（工程「Ｂ」→「Ｃ」→「Ｄ」→「Ｈ」→「Ｆ」→「Ｇ」）とを繰り返す。

ここで、図６の上限電位Ｖ_H、下限電位Ｖ_L、開始電位Ｖｂ、電位掃引速度、電位掃引回数、電位掃引方向は、電極劣化を防ぐなどの目的に応じて任意に設定できる。

次に、ステップＳ８では、反応円弧の規格化処理を行う。

図７（ａ）は、インピーダンス測定プロトコル（ステップＳ４〜ステップＳ７）において得られる各測定電位（Ｖ１，Ｖ２，・・・Ｖｘ）ごとのインピーダンス特性（コール−コールプロット）を示している。また、図７（ｂ）は図７（ａ）に示すインピーダンスに対して、周波数ごとにＤＣ電流値ｉｆを乗じることで得られる、触媒反応に由来する規格化されたインピーダンス（ｉｆ×Ｚ）特性を示している。

ステップＳ８では、インピーダンス測定プロトコル（ステップＳ４〜ステップＳ７）において得られたインピーダンス（図７（ａ））に対して、周波数ごとにＤＣ電流値ｉｆを乗じることにより触媒反応に由来する規格化されたインピーダンス（ｉｆ×Ｚ）特性（図７（ｂ））を算出する。ここで、本実施形態では、ステップＳ５において取得された経時的なＤＣ電流値の変化を、各周波数のインピーダンス測定時における時間と対応付けることで、周波数ごとにインピーダンス測定時（ステップＳ７）のＤＣ電流値ｉｆを取得する。すなわち、図６に示す工程「Ｅ」と工程「Ｈ」において、ＤＣ電流値が経時的に同様に変化すると見做すことにより、工程「Ｈ」におけるＤＣ電流値ｉｆを推定している。

図８は、工程「Ｅ」における測定電位Ｖ１，Ｖ２，Ｖ３におけるＤＣ電流値の経時的変化を示している。図８に示すように、各測定電位においてＤＣ電流（絶対値）は時間経過とともに減少する。ステップＳ８では、工程「Ｅ」で取得されたＤＣ電流値の経時変化を工程「Ｈ」のインピーダンス取得時に当てはめることで、ＤＣ電流値を抽出している。

なお、インピーダンス取得時にリアルタイムでＤＣ電流値を出力可能なポテンションスタットを使用する場合には、ステップＳ４〜ステップＳ５の手順を省略することもできる。

ここで、インピーダンスとインピーダンスを測定したときの測定電流値との関係を説明する。

特に拡散の影響が現れない領域（濃度非依存領域）では電流Ｉと電圧Ｅの間にターフェルの式（（１）式）の関係がある。

Ｅ＝ａ＋ｂ・logＩ ・・・（１）式

（１）式中、ｂはターフェル勾配である。一方、インピーダンス測定では入力電流のＡＣ変調に対する出力電圧の応答に線形関係を仮定して（２）式にてインピーダンスを算出する。

Ｚ＝ｄＥ／ｄＩ ・・・（２）式

（１）式より、濃度非依存領域においては（３）式の関係がある。

ｄＥ／ｄＩ＝ｂ／Ｉ ・・・（３）式

したがって、濃度非依存領域では、インピーダンスＺに測定電流値Ｉを乗じた積Ｉ・Ｚは一定の値となる。

この関係は例えば文献「多孔質電極内の電流分布とインピーダンスの関係の解析 光島重徳、神谷信行、太田健一郎 Electrochemistry P.810-814、2006」に言及されている。この文献中では、周波数を変化させたときのＩ・Ｚ円弧のつぶれ具合を、イオン抵抗／電子抵抗の比を考察することに利用することが記載されている。

このように、インピーダンスＺと測定電流値Ｉの積Ｉ・Ｚのプロットによって電流密度の影響が除外され、視覚的に濃度非依存領域又は濃度依存領域、あるいは特定のターフェル勾配を示す領域の判別ができる。したがって、例えば、燃料電池セルの濃度非依存領域や特定のターフェル勾配を示す領域が、燃料電池セルの動作条件におけるどの範囲に相当するのかを把握することができる。

また、積Ｉ・Ｚのプロットはターフェル勾配などの定常状態におけるパラメータに留まらず、周波数特性という動的なパラメータとして提示される。このため、積Ｉ・Ｚの情報は、例えば、単に電極性能の劣化を判断するパラメータとしてのターフェル勾配を示すだけでなく、例えば、燃料電池の性能を多面的に示すパラメータとして広く利用できる。

次に、ステップＳ９では、ステップＳ８において得られた規格化されたインピーダンス（ｉｆ×Ｚ）特性（図７（ｂ））を用いて、モデルフィッティングや電流分布の解析、拡散評価パラメータの抽出など種々の特性評価を実行することができる。

以下、本発明の電極の評価方法を燃料電池セルにおける規格化インピーダンスの取得および拡散評価パラメータの抽出に適用した例について説明する。

（燃料電池セルにおける規格化インピーダンスの取得）
図９は、本発明の電極の評価装置としての燃料電池セルの測定システムを示すブロック図である。

図９に示す燃料電池セル２００は単セルであり、Ｐｔ／Ｃ触媒が表面に修飾された電解質膜２０１を、ガス拡散層をそれぞれ有するアノード極２０２およびカソード２０３で挟み込んだ構造をとる。またアノード２０２の側には、アノード２０２から分割された分割電極２０４が設けられている。アノード２０２および分割電極２０４は互いに絶縁された状態にある。

測定時には、アノード２０２に水素、カソード２０３に空気もしくは酸素を供給し発電状態とする。電気化学測定は３電極系で行い、カソード特性を評価する際は作用極をカソード２０３、対極をアノード２０２とする。また、アノード極２０２およびカソード極２０３とは絶縁されたアノード側の別の分割電極２０４を参照極(RHE)とする。図９に示すように、アノード２０２、カソード２０３および分割電極２０４は、それぞれポテンショスタット１に接続され、ポテンショスタット１による計測結果はレコーダ２に転送される。

図９に示すように、ポテンショスタット１およびレコーダ２には、制御演算装置３Ａが接続される。制御演算装置３Ａはポテンショスタット１の制御、レコーダ２からの上記計測結果の取得、および、取得された上記計測結果に基づく演算を実行する。

燃料電池セル２００のアノード極２０２には水素が、カソード極２０３には空気もしくは酸素が供給される。これらの供給ガスは図９に示すガス供給装置４を介して与えられる。ガス供給装置４は制御演算装置３Ａに接続され、燃料電池セル２００に供給されるガスのガス濃度、流量、温度等は、ガス供給装置４を介して制御演算装置３Ａにより制御される。

図１０は、燃料電池セルにおける反応抵抗の抽出およびインピーダンスの規格化処理を示すフローチャートである。このフローチャートに示された手順に従って、測定対象である燃料電池セルの反応抵抗を抽出し、規格化インピーダンスを求めることができる。

まず、ステップＳ１０１において、ターフェル勾配が異なるｂ１、ｂ２領域のそれぞれについて、上記のインピーダンス測定プロトコル（ステップＳ４〜ステップＳ７）を用いて、インピーダンスを測定、記録する。ターフェル勾配は領域ｂ１、ｂ２のそれぞれで異なる。領域ｂ１では、白金表面のほとんどが吸着酸素種(酸化皮膜)に覆われ、Temkin-Typeと呼ばれる反応に起因するターフェル勾配が観測される。一方、ｂ２領域では白金表面の酸化物はすべて還元され、Langmuir-Typeと呼ばれるベアな白金表面の反応に起因するターフェル勾配が観測される。

図１１は、ｂ１、ｂ２領域における測定に基づきインピーダンス取得部３３により取得されたインピーダンスを示す図である。図１１において、Ｚｉ１は電流ｉ１におけるインピーダンスを、Ｚｉ２は電流ｉ２におけるインピーダンスを、Ｚｉ３は電流ｉ３におけるインピーダンスを、Ｚｉ４は電流ｉ４におけるインピーダンスを、それぞれ示す。また、電流ｉ１および電流ｉ２は領域ｂ１に、電流ｉ３および電流ｉ４は領域ｂ２に、それぞれ対応し、ｉ１＜ｉ２＜ｉ３＜ｉ４である。

ここでインピーダンスの実軸と交わる点をゼロクロスと定義する。この図１１に示す例では、高周波側ゼロクロス値が０とはなっておらず、直流抵抗成分（Ｒ_DC）が存在していることが分かる。この直流抵抗成分は反応抵抗とは直接関係しないので、ステップＳ１０２において、高周波でゼロクロスした値を直流抵抗値とし、その分を差し引いて原点補正を実施する。

図１２は、ステップＳ１０２において直流補正を施した後のインピーダンスを示す図である。図１２において、高周波ゼロクロス値が原点となっていることが確認できる。

また、図１３（ａ）は、図１２に示すインピーダンス（直流抵抗補正済）の|Z|-f特性、θ-f特性を示す図である。

次に、ステップＳ１０３では、分布定数域の判別を行う。図１２および図１３から判るように、高周波側ではインピーダンスに直線領域が現れていることが分かる。一般的に多孔質電極ではこのような45°の直線が高周波側で現れることが知られている（分布定数域）。分布定数域に関して成立する等式eq.1、eq.2、eq.3を図１３（ｂ）に示す。

これらの等式から、分布定数域は反応抵抗Ｒ_ctには依存していないことが分かる。したがって、分布定数域も除かなければ、後述する規格化によって精度良く反応を分離することができない。

図１３（ａ）には、ステップＳ１０３において、分布定数域と判断する領域の例が示されている。図１３（ａ）の例では、触媒中のプロトン移動抵抗（Ｒ_ion）や電気2重層容量（Ｃ_dl）が反応過程で異なっている可能性も許容し、ゼロクロス周波数から「θ−f特性が45°±5°の領域（グレー帯の部分）」かつ「同じ反応過程で動作点が異なるインピーダンス|Z|の比較において、低周波ゼロクロス値の差の5%以内」に該当する周波数までを分布定数域と判断している。ただし、触媒中のプロトン移動抵抗は異なる反応過程でも変化しないという思想や、直線は必ずしも45°とは限らないなどの思想のもと、分布定数域の具体的な判断基準は適宜変更が可能である。

次に、ステップＳ１０４において、分布定数域を除いたモデルフィッティング処理を行う。

図１４は、ステップＳ１０４におけるモデルフィッティング処理の方法を示す図である。ここでは、フィッティングソフトを用いて、分布定数域以外の周波数帯、すなわち反応抵抗の影響が表れている低周波の周波数帯でのみモデルフィッティングを行う。図１４では、Ｚｉ３についてのみ例示しており、Ｚｉ３についての分布定数域の上限周波数をｆ_ixとしている。ここで、分布特性（伝送線モデル）に起因している直流抵抗成分を総括してＲ_poreと定義する。等価回路モデルとして、Ｒ_poreと反応インピーダンスモデル(Ｚ_f)とを直列に接続した回路（図１４中に記載）を使用する。反応インピーダンスモデル(Ｚ_f)として、実際の反応過程を表現していると考えられるモデルを用いることが望ましく、モデルは適宜選択できる。ここでは、簡易的にR//Cの２次モデルをＺ_fとして用いてフィッティングを行った例を示す。このフィッティングによりＲ_ctとＲ_poreの値が求まり、反応抵抗を反応円弧として簡易的に抽出することができる。

図１５は、分布定数域分を考慮し原点補正したインピーダンスを示す図である。図１５では、図１２に示すインピーダンスに対し、上記手順にて求めたＲ_pore分を差し引き、原点を取り直している。

ここで、ステップＳ１０２およびステップＳ１０４の処理は、いずれも直流抵抗分を取り除く処理であるため、ステップＳ１０４においてまとめて（Ｒ_pore＋Ｒ_DC）として算出してもよい。

次に、ステップＳ１０５では、直流抵抗分を取り除き反応抵抗のみを抽出したインピーダンスに対して規格化処理を施す。ここでは、原点補正したインピーダンス（図１５）にインピーダンス測定動作点の電流値を乗じることでインピーダンスを規格化する。ここでは、上記のステップＳ８の手順を用いて、電流値を各測定点および測定周波数ごとに抽出し、インピーダンスに乗じる。

図１６は、抽出された反応抵抗に対して規格化処理して得たインピーダンス特性（周波数特性）を示す図である。

図１６に示すように、規格化インピーダンスが精度良く反応過程（領域ｂ１、ｂ２）ごとに分離できていることが確認できる。すなわち、Ｚｉ１とＺｉ２に対応するプロット、およびＺｉ３とＺｉ４に対応するプロットが、それぞれ、ほぼ同一曲線上に重なり合うことが判る。したがって、インピーダンス特性に基づいて領域ｂ１、ｂ２の反応過程にあるのかが判定できるとともに、異なる電流値におけるインピーダンスを推定可能となる。また、この規格化インピーダンスの形・大きさなどから反応過程それぞれの特性について論じることが可能となる。また、燃料電池の性能を多面的に示すパラメータとして広く利用できる。

（燃料電池セルにおける拡散評価パラメータの抽出）
燃料電池セルにおける拡散評価パラメータの抽出は、図９に示す測定システムを用いて行われる。また、燃料電池セル２００の構成は図９に示す通りである。

図１７は、測定システムの動作を示すフローチャート、図１８〜図２５は、制御演算装置３Ａによる表示例を示す図である。

図１７のステップＳ２０１では、電圧掃引もしくは電流掃引によりＩＶ特性を測定し、電流の対数と電圧との関係を表現するターフェルプロットを取得する。

図１８は、拡散領域まで測定したターフェルプロットを示す図である。

図１８に示すように、酸素還元反応には一般的に２つのターフェル勾配(ｂ１、ｂ２)が存在することが知られている。拡散領域では直線領域は存在せず、電流が大きくなるにつれ、勾配も変化する。

ここで、反応物の拡散を考慮した電極反応式について述べる。分極が十分大きく、逆反応の電流が無視できる場合、電流ｉおよび過電圧ηは、図１８に示す等式eq.11およびeq.12のように表すことができる。

さらに、等式eq.2を次のように変形し、図１８に示す等式eq.13を得る。

等式eq.13の左辺はインピーダンス(dη/di)に電流値ｉを乗じたものとなっており、ターフェル勾配に相当する値である。インピーダンスに電流値を乗じる操作をここでは規格化処理と呼ぶ。

限界拡散電流ｉ_Lが測定電流ｉよりも十分大きい時、等式eq.13のように近似することができる。この式から拡散の影響がでている領域では反応のターフェル勾配(2.3RT/αnF)に拡散起因の電流に比例する項が加わっていることが分かる。

ここで、上記の（３）式の関係から、濃度非依存領域では、インピーダンスＺに測定電流値Ｉを乗じた積Ｉ・Ｚは一定の値となる。

次に、ステップＳ２０２では、拡散領域の抽出を行う。

ここでは、まず、図１８に示すターフェルプロットを各測定点で微分した値（ターフェルプロットの傾き）を算出する。

図１９（ａ）は、図１８に示すターフェルプロットを各測定点で微分し、その結果の絶対値をターフェル勾配と定義して(電流−ターフェル勾配)特性として表したグラフを示している。ターフェルプロット上の直線領域は、ここではプラトー領域として表され、ｂ１、ｂ２の領域を示す２つのプラトーが確認できる。また、等式eq.13に示されるように、拡散の影響がでている領域（ｂ２を越える領域）では電流に比例する項がターフェル勾配に加算され、比例領域を形成する。ただし、等式eq.13は限界電流ｉ_Lが測定電流ｉよりも十分大きい時に成立する近似式であるので、反応物質(酸素)の濃度が高く、利用率が小さいほどこの比例領域は広くなる。

図１９（ｂ）は、図１９（ａ）のグラフと同一の特性を（電流の対数−ターフェル勾配）特性として表現したものを示す。図１９では測定電流値がリニア軸表示のため、電流が小さい領域で現れるｂ１領域が、グラフ上で見えない場合がある。この場合、図１９（ｂ）のように測定電流値を対数軸表記することでｂ１領域を判断し易くなり、例えば的確にインピーダンスの測定動作点を定めることが可能となる。

図２０は、ｂ２の領域から拡散の影響がでている場合の（電流−ターフェル特性）を示す図である。図２０に示す比例領域を探すことで拡散項を求めることができる。ここで、拡散の影響が出始めている電流をｉ_C,0と定義する。図２０において、比例領域における拡散項に相当する部分を網掛けで示している。後述の図２１においても同様である。

図２０に示すように、ターフェル勾配は電流ｉ_C,0を境界として、ｂ２の領域に相当するターフェル勾配に、電流に比例する拡散項が付加された値をとる。

図２１は、ｂ１の領域から拡散の影響がでている場合の（電流−ターフェル特性）を示す図である。この場合にも、網掛けで示される比例領域を探すことで拡散項を求めることができる。図２１に示すように、ターフェル勾配は電流i_c,0を境界として、ｂ１の領域に相当するターフェル勾配およびｂ２の領域に相当するターフェル勾配に、それぞれ電流に比例する拡散項が付加された値をとる。拡散項により生ずる直線の傾きは、すべての領域で一定である。

図２２は、拡散補正を施した（電流−ターフェル勾配）特性を示す図である。拡散補正により、図２０および図２１の拡散項（網掛けで示す部分）が除去され、ｂ１およびｂ２の領域を示すプラトーのみが現れる。なお、拡散補正は後述するステップＳ２０５における処理に相当する。

次に、ステップＳ２０３では、インピーダンス測定および測定されたインピーダンスに対する規格化処理を行う。ここでは、上記のインピーダンス測定プロトコル（ステップＳ４〜ステップＳ７）を用いて、インピーダンスを測定、記録するとともに、記録されたインピーダンスに対して規格化処理を施す。

図２３は、ｉ１〜ｉ７の各動作点で測定した規格化インピーダンスＺｉ１〜Ｚｉ７を示す図である。図１９〜図２２には、個々の規格化インピーダンスＺｉ１〜Ｚｉ７を測定した動作点が矢印（↓）により表されている。図１９および図２０では、規格化インピーダンスＺｉ１〜Ｚｉ２がｂ１領域に、規格化インピーダンスＺｉ３〜Ｚｉ４がｂ２領域に、規格化インピーダンスＺｉ５〜Ｚｉ７が領域ｂ２および拡散領域に、それぞれ対応している。また、図２１では、規格化インピーダンスＺｉ１がｂ１の領域に、規格化インピーダンスＺｉ２がｂ１領域に、規格化インピーダンスＺｉ３〜Ｚｉ３がｂ１領域および拡散領域に、それぞれ対応している。

なお、図２３は、拡散がｂ２領域の途中から現れる場合（図１９〜図２０に相当）について示している。

図２３に示すように、ｂ１の領域に属するインピーダンスおよびｂ２の領域に属するインピーダンスは、規格化によってそれぞれ共通の曲線群にまとまり、拡散の影響が加わると動作点電流が大きくなるにつれ規格化インピーダンスも大きくなっていることがわかる。すなわち、拡散の影響により、ｂ２の領域の曲線群から離れた位置に曲線が移動する。

次に、ステップＳ２０４では、拡散性能評価パラメータ（Ｃ_(f))を算出する。

ここで、拡散性能評価パラメータ（Ｃ_(f))の算出方法について述べる。図２３に、等式eq.14として、規格化インピーダンスにおける評価パラメータの算出式を示す。

この式は、拡散領域にある２点の規格化インピーダンス（ｉＺ_(f)）の差を、そのそれぞれの動作点電流値の差で除している。したがって図２０あるいは図２１に示す比例領域における直線の傾きを求めていることに相当する。すなわち電流に依存しない評価パラメータ値として抽出することができる。

図２４は、算出された拡散性能評価パラメータを例示する図である。このように拡散のみに起因する周波数情報を持った評価パラメータを求めることで、拡散性能についてより広い視点から議論することが可能となる。また、拡散性能評価パラメータを用いることで、実際に測定していない動作点におけるインピーダンスを推定することもできる。このような処理は推定手段の機能に相当する。

次に、ステップＳ２０５では、拡散性能評価パラメータを用いて拡散項の補正を行う。図２４に示す等式eq.15は、拡散性能評価パラメータを用いた拡散項の補正方法を示す補正式である。

この式は評価パラメータ（Ｃ_(f)）に測定動作点の電流値（ｉ_X）と、拡散影響が始まる電流値（ｉ_C,0）の差を乗じたものを、元の規格化インピーダンス（（ｉＺ_(f)）_C,X）から差し引く補正を行うものである。

すなわち拡散領域の任意の動作点の拡散項を求めて差し引くことで、元の規格化インピーダンスを補正している。

図２５は、拡散の補正を施した規格化インピーダンスを示す図である。図２３と比較すれば明らかなように、上記手法によって拡散項の影響が排除され、拡散領域での測定結果を示す曲線群がｂ２の領域の規格化インピーダンスを示す位置に移動することがわかる。これにより、拡散性能が悪くｂ２領域が狭い場合やｂ２領域が見えない場合においても、的確にｂ２領域の規格化インピーダンスを算出し評価を行うことができる。

また、等式eq.14のＣ_(f)と等式eq.13の右辺第2項の−ｉの係数とは等価（直線領域の傾きに相当）であるため、図２５に示す等式eq.16と表され、等式eq.16を変形して図２５に示す等式eq.17を得ることで、限界拡散電流（ｉ_L）についても周波数情報を持った特性値（ｉ_L(f)）として推測することができる。
以上のように、ターフェルプロットの微分値（ターフェル勾配）を、（電流−ターフェル勾配）特性として表現し、プラトーおよび電流比例領域を探すことで、ｂ１領域、ｂ２領域、および拡散領域を簡便に分離することができるほか、拡散影響開始電流についても求めることができる。

また、電流に依存しない周波数情報を持った拡散性能評価パラメータを算出することで、従来評価が困難であった拡散性能の広い視点からの議論が可能となる。また、限界拡散電流を、周波数情報をもった特性値として推測できることから、最適運転制御時の新たな情報として活用できる可能性がある。

さらに、拡散性能評価パラメータを用いて規格化インピーダンスの拡散補正ができることから、拡散性能が悪くｂ２領域が小さい、もしくはｂ２領域が見えない電極反応においても、反応過程（領域）ごとに規格化インピーダンスを分離することができる。

さらにまた、ターフェル式に従う反応と拡散の影響とが混合する電極反応において、拡散性能評価パラメータの算出および規格化インピーダンスの拡散補正ができるため、燃料電池の評価のみならず溶液系などの酸素還元反応以外の電極反応においても広く応用できる。

以上説明したように、上記実施形態の評価方法によれば、高速電位掃引に続いて、ＬＳＶの掃引速度依存性から決定した掃引速度で電位制御するため、触媒表面が清浄で、かつ電位履歴がキャンセルされた状態にあるので、触媒活性サイトが最大限に活かされ、かつバックグラウンドの影響による誤差が少なく、再現性のよい触媒活性評価（ＤＣ特性およびＡＣ特性についての評価）が可能となる。

また、インピーダンス測定中の酸化物やアニオン・有機物等の表面吸着による経時的な触媒活性サイトの変化に由来する触媒活性変化を加味し、測定中のＤＣ電流値を用いて測定周波数ごとに規格化処理を施すことで、触媒反応に由来するインピーダンスのみをより精度よく抽出することが可能となる。

本発明の電極の評価方法は、バトラー・ボルマー式に基づいた反応であれば、酸素還元反応以外の電極反応系についても広く適用することができる。

図２６は、図６に示したインピーダンス測定プロトコルにＣＶを実施する工程を挿入した例を示す図である。

図２６に示すように、図６に示したインピーダンス測定プロトコルにおける工程「Ｆ」に代えて、工程「Ｊ」を挿入することにより、工程「Ｅ」におけるＣＡおよび工程「Ｈ」におけるＥＩＳの前後に、工程「Ｂ」におけるＣＶと、工程「Ｊ」におけるＣＶとを配置している。このため、それぞれのＣＡまたはＥＩＳ前後のＣＶ波形を比較することにより、酸化物由来や有機物由来のピーク強度変化やピーク電位変化を捉えることができ、その測定電位においてどのような表面吸着種がどの程度付着しているかを判断することができる。

以上説明したように、本発明の電極の評価装置等によれば、電流値取得手段で取得された直流電流値の経時的変化を適用して、インピーダンス取得手段により取得されたインピーダンスに直流電流値を乗じた規格化インピーダンスを算出するので、正確な規格化インピーダンスを求めることができる。

本発明の適用範囲は上記実施形態に限定されることはない。本発明は、電気化学特性に基づいて電極の特性を評価する電極の評価装置等に対し、広く適用することができる。

３１ 負荷制御部（電位制御手段）
３２ 直流特性取得部（直流値取得手段）
３４ インピーダンス取得部（インピーダンス取得手段）
３５ 規格化インピーダンス算出部（規格化インピーダンス算出手段）

Claims (7)

1. 電気化学特性に基づいて電極の特性を評価する電極の評価装置において、
   電極に与える電位制御を行う電位制御手段と、
   前記電位制御手段により与えられる特定の直流動作条件の下で前記電極のインピーダンス特性を取得するインピーダンス取得手段と、
   前記特定の直流動作条件の下における直流電流値の経時的変化を取得する電流値取得手段と、
   前記電流値取得手段で取得された前記直流電流値の経時的変化を適用して、前記インピーダンス取得手段により取得されたインピーダンスに前記直流電流値を乗じた規格化インピーダンスを算出する規格化インピーダンス算出手段と、
   を備えることを特徴とする電極の評価装置。
2. 前記電流値取得手段は、前記特定の直流動作条件の下におけるクロノ・アンペロメトリー、または前記インピーダンス取得手段により前記インピーダンス特性が取得される間における、前記電極に接続されたポテンションスタットからの制御された直流電流に基づいて前記直流電流値の経時的変化を取得することを特徴とする請求項１に記載の電極の評価装置。
3. 前記電位制御手段は、前記特定の直流動作条件への移行に先行して、高速掃引を繰り返す工程、前記特定の直流動作条件まで電位制御する工程、およびバックグラウンド電流を開放するために電位を維持して待機する工程を順次実行することを特徴とする請求項１または２に記載の電極の評価装置。
4. 高速掃引直後におけるリニア・スイープ・ボルタンメトリーの測定を複数の掃引速度について行うことにより、前記特定の直流動作条件に至るまでの掃引速度を決定する掃引速度決定手段を備え、
   前記電位制御手段は、前記掃引速度決定手段により決定された前記掃引速度で前記特定の直流動作条件まで電位制御することを特徴とする請求項１〜３のいずれか１項に記載の電極の評価装置。
5. 前記掃引速度決定手段は、前記複数の掃引速度のうち定常状態が得られる最も遅い掃引速度を前記特定の直流動作条件に至るまでの掃引速度として決定することを特徴とする請求項４に記載の電極の評価装置。
6. 前記特定の直流動作条件の前後にサイクリック・ボルタンメトリーを実施し前後の波形を比較することで、前記特定の直流動作条件における表面吸着についての評価を行う手段を備えることを特徴とする請求項１〜５のいずれか１項に記載の電極の評価装置。
7. 電気化学特性に基づいて電極の特性を評価する電極の評価方法において、
   電極に与える電位制御を行う電位制御ステップと、
   前記電位制御ステップにより与えられる特定の直流動作条件の下で前記電極のインピーダンス特性を取得するインピーダンス取得ステップと、
   前記特定の直流動作条件の下における直流電流値の経時的変化を取得する電流値取得ステップと、
   前記電流値取得ステップで取得された前記直流電流値の経時的変化を適用して、前記インピーダンス取得手段により取得されたインピーダンスに前記直流電流値を乗じた規格化インピーダンスを算出する規格化インピーダンス算出ステップと、
   を備えることを特徴とする電極の評価方法。