第3章: 速度論 — 過電圧とTafel解析

熱力学が許した反応がそれでも起こらないのはなぜか、そして急がせるための代価をどう読むか

📖 読了時間: 20-25分 📊 難易度: 初級 💻 コード例: 0個 📝 演習問題: 0問

ビデオ講義

このビデオは以下のテキストと同じ内容をカバーしています。お好みの学習形式をお選びください。

🌐 JP | 🇬🇧 EN | Last sync: 2026-08-20

材料科学道場 > 電気化学入門 > 第3章

第2章は、心地よい数値で終わりました。水は \(1.23\) V で分解する。ダニエル電池は \(1.10\) V を供給する。Nernst式は、濃度が変わったときにそれらの数値がどれだけ動くかを厳密に教えてくれる。熱力学は この反応は許されるのか、そして電圧をどれだけ要求し、あるいは支払ってくれるのか という問いに答えました——しかも完全に、表に載った量だけから、電極が何でできているかにはいっさい触れずに。

その完全さこそが罠です。水電解装置を組み立て、\(1.23\) V をかけてみてください。事実上、何も起こりません。「少しは起こる」のではありません——気づくようなことは何も起こらないのです。有用な電流を引き出すには、熱力学的な要求を超えて0.5 V以上を押し込まなければならないこともあり、しかもどれだけ余分に必要かは、電極が何でできているかに完全に依存します。第2章にはどこにも現れなかった変数に、です。

本章はその隔たりについての章です。この隔たりには名前があります——過電圧——そしてそれは第2章の理論の欠陥ではありません。熱力学は速度について何かを言うと主張したことは一度もないのです。最終状態にいくらかかるかを述べたのであって、そこへどれだけ速く到達するかを述べたのではなく、この2つは本当に独立した問いです。丘の頂上にある巨岩は転がり落ちます。熱力学がそれを保証します。それが今日の午後に転がるのか1万年後なのかは、丘の高さではなく、経路の形についての問いです。

電気化学において、経路の形こそが工学のすべてが住む場所です。あらゆる触媒、あらゆる電極材料、白金やイリジウムに費やされるあらゆる1ドルは、速度論に費やされています。ですから、速度論が何を買えて何を買えないのかについては、精確であるだけの価値があります。

3.1 平衡は速度の約束ではない

まず、電極において平衡が実際に何を意味するのかから始めましょう。この語の日常的な読み方は誤解を招くからです。

亜鉛電極を硫酸亜鉛溶液に入れて放置します。第2章によれば、そこにはある電位が確立されます——Nernst式が与える平衡電位 \(E_{\text{eq}}\) です。これを静的な状況として思い描くのは容易です。何も起きておらず、系は落ち着き、計測器は1つの数値を示している、と。

その像は誤りであり、正しい像こそが本章のあらゆるものの土台です。

平衡においては、反応の両方向がフルスピードで走っています。亜鉛原子は金属を離れて \(\text{Zn}^{2+}\) として溶液に入っていき、\(\text{Zn}^{2+}\) イオンは表面に到達して金属として析出しています。どちらも止まってはいません。平衡において真であるのは、その2つの速度が 厳密に等しい ということだけであり、だからこそ正味の変化はゼロになり、正味しか見ることのできない外部の計測器は電流ゼロを示すのです。

\[ \text{Zn} \;\rightleftharpoons\; \text{Zn}^{2+} + 2e^- \]

平衡とは つり合い であって、停止 ではありません。そしてこれは直ちに、熱力学には尋ねようのない問いを引き起こします。そのつり合った往来は、どれくらい速いのか?

2つの電極が厳密に同じ平衡電位に座ることはありえます——熱力学的に区別不能で、Nernst式も同一、\(\Delta G\) も同一です——それでいて、界面を両方向に横切る往来の量はまるで違いうるのです。一方はイオンが絶えず飛び交っているかもしれず、もう一方は界面があまりに鈍く、ほんのわずかしか横切らないかもしれません。どちらも同じ電圧を示します。どちらも平衡にあります。それらは同じ電極ではなく、その違いが、あなたがそれらで行おうとするあらゆることを支配することになります。

その違いこそ、続く数節の主題です。しかしその前に、どれだけ強く押しているかを測る量を導入しましょう。

3.2 過電圧: どこかへ行くための代価

過電圧 は、ほとんど気恥ずかしいほど単純に定義されます。電極に実際に印加した電位と、そこで起きている反応の平衡電位との差です。

\[ \eta \;=\; E_{\text{applied}} \;-\; E_{\text{eq}} \]

\(\eta = 0\) なら、あなたはちょうど平衡に座っており、正味電流はゼロです。\(\eta\) を正に押せば酸化方向を駆動し、負に押せば還元方向を駆動します。符号の約束は第1章の還元電位の約束に従います。正の \(\eta\) は アノード的（酸化、電子が溶液から電極へ出ていく）、負の \(\eta\) は カソード的（還元）です。

\(\eta\) について重要なのは定義ではなく、その会計です。過電圧は 無駄になったエネルギー です。その1ボルト1ボルトが、通した電荷に掛け算され、生成物ではなく熱として出ていきます。熱力学的要求の \(1.23\) V に過電圧 \(0.6\) V を加えて運転している電解装置は、電気入力のおよそ3分の1を部屋を暖めることに変換しています。過電圧を下げることは細かな改良ではありません——それは、存在するあらゆる電気化学デバイスの効率に対する主要なてこなのです。

そして過電圧は1つのものではありません。それは少なくとも3つの物理的に異なる罰金の和であり、それぞれ原因も、処方箋も、特徴的な兆候も異なります。

📚 過電圧への3つの寄与

寄与 物理的原因 電流への依存の仕方 減らし方
活性化 \(\eta_{\text{act}}\) 電荷移動段階そのものが界面にエネルギー障壁をもつ 対数的——これが本章のTafel挙動 より良い触媒、より大きな実表面積、より高い温度
濃度 \(\eta_{\text{conc}}\) 反応物が、輸送で補給できるより速く表面で消費される 限界電流に近づくまでは小さく、そこで発散する 撹拌する、濃度を上げる、拡散層を薄くする、フローセルにする
オーム \(iR\) 電解液、膜、接点、配線のすべてが抵抗をもつ 電流に厳密に線形 より導電性の高い電解液、より短い経路、より太い接点、電子的補償

先へ進む前に、この表から取り出しておく価値のあることが3つあります。

それらは足し合わさります。 実際のセルが必要とする総駆動電圧は、熱力学的な最小値に、両方の電極における3つの罰金すべてを加えたものです。

\[ E_{\text{cell}} \;=\; E_{\text{thermo}} \;+\; |\eta_{\text{anode}}| \;+\; |\eta_{\text{cathode}}| \;+\; iR \]

第5章では、水電解装置についてまさにこの積み上げを構築し、どの項が支配的かを示します。

それらは電流依存性が異なり、それこそが見分ける手立てです。 これは聞こえるより有用です。電流を2倍にして余分な罰金がちょうど2倍になるなら、あなたが見ているのはオーム損失です。どこから始めたかによらず、電流を2倍にする代価が一定のミリボルト数なら、あなたが見ているのは活性化損失です——対数だけがその性質をもち、他のどれももちません。そして押し込むにつれて罰金が突然無限へ走り去るなら、輸送の限界に当たったのです。よく設計された測定は、これらの兆候を利用して3つを分離します。

そのうち化学であるものは1つだけです。 オーム損失はセルの配管の性質です。濃度損失は流体力学の性質です。活性化過電圧だけが その特定の表面上の反応 の性質であり、したがって触媒が扱えるのはこれだけです。誰かが触媒は「過電圧を下げる」と言うとき、それは活性化項のことを指しており、本章の残りはまさにその項についてのものです。

3.3 交換電流密度: 平衡はどれくらい賑わっているか

3.1節が提起した問いに戻りましょう。平衡において、界面をどれだけの往来が横切っているのか?

その往来には名前があります。交換電流密度 \(i_0\) は、電極が平衡に座っているときに 各 方向へ流れている一方向電流の大きさです。アノード分極電流とカソード分極電流はそれぞれ \(i_0\) に等しく、打ち消し合い、正味はゼロになります。

\[ i_{\text{anodic}} = i_0, \qquad i_{\text{cathodic}} = -i_0, \qquad i_{\text{net}} = 0 \]

単位は単位面積あたりの電流——A/cm² が通常の選択です——というのも、重要なのは界面の単位あたりの往来であって、総量ではないからです。総量は、あなたの電極がどれだけ大きいかを教えてくれるにすぎません。

覚えておくべき直観はこうです。\(i_0\) は平衡がどれくらい賑わっているかを測ります。 大きな \(i_0\) は、反応が絶えず両方向に走っている界面——往来がちょうど釣り合った、混雑した対面通行の道路——を表します。小さな \(i_0\) は、横断が稀な出来事である界面を表します。同じつり合いですが、どちらの方向にもほとんど何も起きないことによって達成されたつり合いです。

さてその帰結こそが、この量が重要である理由のすべてです。賑わった平衡は、崩しやすいのです。 巨大な双方向の流れがすでに走っているなら、電位をわずかに小突けばそれが少し傾き、直ちに大きな正味電流が生じます。平衡がほとんど死んでいるなら、正味の流れが測定可能になるまでに非常に強く押さなければなりません。したがって、

これが電気化学における触媒作用の操作的定義です。触媒とは、あなたが関心をもつ反応について大きな \(i_0\) をもつ表面のことなのです。

📚 数値についての率直な注記

教科書や論文は、特定の金属上の特定の反応について交換電流密度を表にしており、それを引用したくなります。本シリーズはそうしません。その理由は、ただ受け入れるよりも理解しておく価値があります。

交換電流密度は 表面にきわめて敏感 です。同じ金属の2つの試料が、どう研磨したか、どう洗浄したか、あるいはどれだけ長く空気中に置かれていたかだけで異なっていても、桁が違う \(i_0\) 値を与えうるのです。電解液中の微量不純物は、まさに化学を担う席に吸着します。結晶面も効きます。文献に報告される値はそれに応じて散らばり、表面調製の全体を添えずに引用された特定の数値は、ほとんど無意味に近いのです。真剣な電極触媒の研究が、電極の公称組成を信じるのではなく、電気化学測定を表面敏感なキャラクタリゼーション——とりわけXPSで、分光分析入門 シリーズで扱われています——と組み合わせるのも、これが理由です。

頑健であるもの——そして持ち帰るべきもの——は その幅 です。同じ反応について、優れた触媒と貧弱な触媒のあいだで、\(i_0\) は 何桁も 違いうるのです。それが分野全体を駆動している定性的事実であり、3.8節ではそのような幅がボルト単位でいくらの代価になるかを厳密に計算します。個々の桁の数字は、誰かの特定の表面についての誰かの丁寧な測定であって、転用は効きません。

3.4 Butler–Volmer式

これで2つの材料が揃いました。平衡がどれくらい賑わっているかの尺度（\(i_0\)）と、どれだけ強く押しているかの尺度（\(\eta\)）です。Butler–Volmer式 は、それらを正味電流に結びつけます。

\[ i \;=\; i_0 \left[ \exp!\left( \frac{(1-\alpha)\, n F \eta}{RT} \right) \;-\; \exp!\left( \frac{-\alpha\, n F \eta}{RT} \right) \right] \]

いかめしく見えます。そうではありません。これは実のところ、2つの言明を並べて書いただけのものだからです。

2つの項として、方向ごとに1つずつ読んでください。 第1の指数関数はアノード（酸化）分極電流です。第2はカソード（還元）分極電流です。正味電流はその差です。両者は電荷を逆向きに運ぶからです。それ以外はすべて、各方向が電位にどれだけ強く応答するかについての帳簿づけです。

なぜ指数関数か? 電荷移動段階が活性化障壁をもち、障壁を越える速度がその高さに指数的に依存するからです——化学におけるあらゆる熱活性化過程を支配するのと同じArrheniusの論理です。電気化学的な 障壁が特別なのは、装置の前面にあるつまみでその高さを変えられる という点です。電極電位を \(\eta\) だけずらすと電子のエネルギーがずれ、それが電子を動かす反応の障壁をずらします。化学の反応で活性化エネルギーを連続的にダイヤルできるものはごくわずかであり、これこそ電気化学をあれほど精密な実験科学にしている性質です。

\(\alpha\) とは何か? 移動係数、または対称因子と呼ばれる0と1のあいだの数で、印加した過電圧が2つの方向にどう分配されるかを述べます。\(\eta\) を印加すると、一方の方向の障壁は \(nF\eta\) のうち \(\alpha\) の割合だけ下がり、他方の障壁はその \((1-\alpha)\) だけ上がります。物理的には、反応経路に沿って遷移状態がどこに座るかを符号化しています。\(0.5\) 付近の値は遷移状態がおよそ中間にあることを意味し、したがって電位は、逆方向を妨げるのとほぼ同じだけ順方向を助けます。\(\alpha = 0.5\) は慣習的な出発点の仮定であり、それは モデル化の選択 であって、測定された自然定数ではありません——実際の系はずれますし、多段階機構は、単一段階の描像が許す範囲をはるかに外れた実効値を生みうるのです。

\(\eta = 0\) でこの式が何を言うかに注目してください。両方の指数関数が1に等しくなり、角括弧は消え、\(i = 0\) となります。平衡がそのまま落ちてきます。しかし \(i_0\) は依然として前に座り続けており、正味電流には見えない往来を記述しているのです。

📚 Butler–Volmer式が仮定していること

この式は働き者であり、たいていの働き者と同様、尋ねられない限りほとんど何も打ち明けません。仮定しているのは、

これらの留保のどれも、この式が正しい第一のモデルであることを妨げません。ただしそれらは、Butler–Volmerフィットが挙動の記述であって、機構の証明ではないことを意味します。

3.5 2つの極限: 直線と、対数の直線

Butler–Volmer式は、書かれたままの形では扱いやすいとは言えません。その価値は、重要な2つの領域において、グラフから読み取れる何かへと潰れることにあります。

小さな過電圧: 線形領域。 \(|\eta|\) が熱的スケール \(RT/F\) ——室温で数十ミリボルト——に比べて小さいとき、両方の指数関数を \(e^x \approx 1 + x\) と展開できます。2つの項のあいだで1どうしが打ち消し合い、\(\alpha\) と \((1-\alpha)\) は足して1になり、残るのは直線です。

\[ i \;\approx\; i_0 \, \frac{nF}{RT} \, \eta \]

過電圧に比例する電流とは、変装したオームの法則であり、その比例定数はコンダクタンスの単位をもちます。その逆数は 電荷移動抵抗 と呼ばれます。

\[ R_{\text{ct}} \;=\; \frac{RT}{nF\, i_0} \]

これは本当に有用な結果です。平衡の近傍では、電気化学界面は抵抗器のようにふるまい、その抵抗値は \(i_0\) に反比例します。小振幅の手法——とりわけインピーダンス分光——が、系を平衡から遠ざけることも、その過程で表面を変えることもなく \(i_0\) を取り出せるのは、これが理由です。3.8節では、この近似が破れるまでにどこまで引き伸ばせるかを測ります。

大きな過電圧: Tafel領域。 \(\eta\) を一方向に大きく押すと、一方の指数関数が成長し、他方は潰れます。大きな正の \(\eta\) では、カソード項が無視できるようになり、アノード項だけが生き残ります。

\[ i \;\approx\; i_0 \exp!\left( \frac{(1-\alpha) n F \eta}{RT} \right) \]

両辺の対数をとれば、指数関数は直線になります。

\[ \log_{10}|i| \;=\; \log_{10} i_0 \;+\; \frac{(1-\alpha) n F}{2.303\, RT}\, \eta \]

これが Tafel式 であり、\(\log_{10}|i|\) を \(\eta\) に対してプロットしたものが Tafelプロット です。これは電気化学速度論で最もよく使われる単一のグラフであり、その理由はまったく実務的です。2つの未知数をもつ曲線を、傾き が \(\alpha\) を与え 切片 が \(i_0\) を与える直線に変えてしまうのです。両方のパラメータが、1回の測定から、線を1本引くことで得られます。

3.6 Tafel勾配と、118 mVはどこから来るのか

電気化学者は通常、Tafelの関係を逆向きに、すなわちボルトあたりの電流の桁数ではなく、電流1桁あたりのボルトとして引用します。そうして得られる数値が直接的に意味をもつ量だからです。電流を10倍にするには、どれだけ余分な過電圧が要るのか?

\[ \eta \;=\; a \;+\; b \log_{10}|i|, \qquad b \;=\; \frac{2.303\, RT}{\alpha n F} \]

量 \(b\) が Tafel勾配 で、単位はミリボルト毎デケードです。\(R = 8.314\) J/(mol·K)、\(F = 96485\) C/mol、\(T = 298.15\) K、そして \(\alpha = 0.5\) と \(n = 1\) を入れれば、本章のコードが確認することになる数値が得られます。およそ118 mV/decade です。

その数値の部品がどこから来るのかを見ておく価値があります。同じ部品が第2章のNernst勾配を作ったからです。

つまり有名な118 mVは、どんな材料についての実験事実でもありません。それは \(59 / 0.5\) であり、印加電位が順方向と逆方向の障壁に均等に分配されると仮定したことの直接の帰結です。仮定を変えれば、この数値も完全に予測可能な仕方でそれに従って変わります——3.8節でその対応を表にします。

📚 現実のTafel勾配を読む

勾配は \(\alpha\) と、律速段階に先立つ電子の数に依存するため、測定されたTafel勾配はしばしば 機構の指紋 として使われます。118、59、39、30 mV/decade はいずれも、見分けのつく組み合わせに対応します。これは正当で広く使われる推論ですが、同時に過剰に読み込みやすくもあります。3つの警告を。

Tafel勾配の誠実な使い方は比較的です。同じセル、同じ電解液、同じフィッティング窓で、電極材料だけを変える。そのときはじめて、差が意味をもちます。

3.7 触媒が実際に変えるもの、そして決して触れられないもの

これで電気化学における最も重要な構造的事実を述べることができます。それは一文で済みます。

触媒は速度論を変える。熱力学を変えることはできない。

Butler–Volmer式を振り返ってください。触媒が影響を与えるものはすべて——\(i_0\)、\(\alpha\)、したがってTafel勾配——電流の式の内側に座っています。\(\eta\) 軸の原点である平衡電位 \(E_{\text{eq}}\) は、調整可能な量としてはそもそも現れません。それは全反応の \(\Delta G\) によって固定されており、\(\Delta G\) は始状態と終状態の差です。触媒はそのどちらでもありません。

したがって、具体的には、

この分業こそ、この分野が今のような姿をしている理由です。熱力学はどの反応が試みる価値をもつかを教え、目標を設定します。速度論は目標にどこまで近づけるかを教え、そこに研究努力の事実上すべてが注がれます。OERの計算化学 シリーズはまさにこの枠組みを取り、さらに押し進めます。候補となる表面上の吸着中間体の自由エネルギーを計算し、そこから理論過電圧を導くことで、「この触媒はどれだけ良いか」を、材料が合成される前に走らせられる計算に変えるのです。

もう1つ帰結があり、それは静かに残酷なものです。Tafelの関係は対数的なので、触媒を改良しても、電流については収穫逓減、電圧については収穫一定 になります。\(i_0\) を1桁稼げば、過電圧はちょうどTafel勾配1つ分——対称な場合でおよそ118 mV——だけ買えます。良い触媒から出発しようと、ひどい触媒から出発しようと、変わりません。逆に言えば、10倍というのは実現するには非常に大きな化学的改良であり、それで買えるのは0.1 Vです。この算術が、工業電解における過電圧が緩やかに、少しずつしか下がってこなかった理由と、それが決してゼロにならなかった理由を説明します。

3.8 ハンズオン: Butler–Volmer式からTafel勾配を取り戻す

以下のコードは、1つの誠実な仕事をします。Butler–Volmer式から電流–過電圧データを生成し、高過電圧の枝以外をすべて捨て、実験家がそうするように対数空間で直線をフィットし、返ってきた勾配が代数の予言したものかどうかを確かめます。その途上で、線形近似がどこまで信頼できるかを測り、\(i_0\) の1桁がボルト換算でいくらの価値をもつかを計算します。

入力のうち2つは測定されたものではなく選ばれたものであり、そのことはコード中に明記されています。\(\alpha = 0.5\) は慣習的な対称性の仮定であり、交換電流密度は代役です。代役でも害がないのは、\(i_0\) が 勾配からは完全に消える からです。それは切片を定めるのであって、傾きを定めるのではありません。Tafel勾配が交換電流密度より持ち運びのきく量であるのは、まさにそのためです。

import numpy as np

# ---------------------------------------------------------------
# Butler-Volmer -> Tafel slope, recovered by fitting
#
# Inputs are only universal constants plus two CHOSEN model
# parameters (alpha and i0). Nothing here is an experimental
# measurement: i0 is a stand-in whose value cancels out of the
# slope entirely, which is exactly the point of the exercise.
# ---------------------------------------------------------------
F = 96485.0        # C/mol      Faraday constant
R = 8.314          # J/(mol K)  gas constant
T = 298.15         # K          25 degrees Celsius
n = 1              # electrons per elementary step

ALPHA = 0.5        # symmetry factor: chosen, not measured
I0 = 1.0e-3        # A/cm^2, a stand-in exchange current density

f = n * F / (R * T)
print("Step 1: the thermal voltage scale")
print(f"  RT/F            = {R * T / F * 1000:.2f} mV")
print(f"  f = nF/(RT)     = {f:.3f} 1/V")
print(f"  2.303 RT/F      = {2.303 * R * T / F * 1000:.2f} mV/decade   (Nernst slope)")
print()

# --- 2. The Butler-Volmer current-overpotential curve -----------
eta = np.linspace(-0.4, 0.4, 4001)           # V
i_anodic = I0 * np.exp((1.0 - ALPHA) * f * eta)
i_cathodic = -I0 * np.exp(-ALPHA * f * eta)
i_total = i_anodic + i_cathodic

print("Step 2: both directions, and their sum")
print(f"{'eta (mV)':>10} {'i_anodic':>14} {'i_cathodic':>14} {'i_total':>14} {'|back/fwd|':>12}")
print("-" * 68)
for e_mV in [0, 10, 25, 50, 100, 200, 300]:
    k = int(np.argmin(np.abs(eta - e_mV / 1000.0)))
    ratio = abs(i_cathodic[k] / i_anodic[k])
    print(
        f"{eta[k] * 1000:10.1f} {i_anodic[k]:14.4e} {i_cathodic[k]:14.4e} "
        f"{i_total[k]:14.4e} {ratio:12.3e}"
    )
print()

# --- 3. The linear (low-overpotential) regime -------------------
# For |eta| << RT/F, exp(x) ~ 1 + x and the two exponentials
# collapse to a straight line: i ~ i0 * f * eta. The slope of that
# line has units of a conductance; its reciprocal is the
# charge-transfer resistance Rct = RT/(nF i0).
print("Step 3: near equilibrium the curve is a straight line")
Rct = R * T / (n * F * I0)
print(f"  Rct = RT/(nF i0) = {Rct:.3f} ohm cm^2   (for the stand-in i0 above)")
for e_mV in [1, 2, 5, 10, 20, 50]:
    e = e_mV / 1000.0
    exact = I0 * (np.exp((1 - ALPHA) * f * e) - np.exp(-ALPHA * f * e))
    linear = e / Rct
    err = 100.0 * (linear - exact) / exact
    print(f"  eta = {e_mV:3d} mV: exact {exact:.4e}, linear {linear:.4e}, error {err:+6.2f} %")
print()

# --- 4. The Tafel (high-overpotential) regime -------------------
# Once one exponential dominates, log10|i| is linear in eta.
# Fit that branch and see what slope comes back.
mask = eta >= 0.15                            # anodic branch, well past the linear region
x = eta[mask]
y = np.log10(np.abs(i_total[mask]))
slope, intercept = np.polyfit(x, y, 1)        # decades per volt, and log10(i) at eta = 0

fitted_tafel_mV = 1000.0 / slope              # mV per decade
predicted_tafel_mV = 1000.0 * 2.303 * R * T / ((1.0 - ALPHA) * n * F)
i0_recovered = 10.0**intercept

print("Step 4: fit log10|i| vs eta on the high-overpotential branch")
print(f"  fit window            : eta = {x[0] * 1000:.0f} to {x[-1] * 1000:.0f} mV")
print(f"  fitted slope          : {slope:.4f} decades/V")
print(f"  fitted Tafel slope    : {fitted_tafel_mV:.2f} mV/decade")
print(f"  predicted 2.303RT/(anF): {predicted_tafel_mV:.2f} mV/decade")
print(f"  intercept -> i0       : {i0_recovered:.4e} A/cm^2 (input was {I0:.4e})")
print()

# --- 5. What the symmetry factor does to the slope --------------
print("Step 5: Tafel slope vs transfer coefficient (anodic branch)")
print(f"{'alpha_a':>9} {'mV/decade':>12}")
print("-" * 22)
for a in [0.2, 0.3, 0.4, 0.5, 0.6, 1.0, 1.5, 2.0]:
    print(f"{a:9.2f} {1000.0 * 2.303 * R * T / (a * F):12.1f}")
print()

# --- 6. What a catalyst buys you --------------------------------
# Same equation, same equilibrium potential, different i0.
# How much overpotential does it take to reach 10 mA/cm^2?
TARGET = 10.0e-3                              # A/cm^2, a conventional benchmark current
b = 2.303 * R * T / ((1.0 - ALPHA) * n * F)   # V/decade
print(f"Step 6: overpotential needed to reach {TARGET * 1000:.0f} mA/cm^2")
print(f"{'i0 (A/cm^2)':>14} {'eta needed (mV)':>18}")
print("-" * 34)
for i0_try in [1e-3, 1e-5, 1e-7, 1e-9, 1e-11]:
    eta_needed = b * np.log10(TARGET / i0_try)
    print(f"{i0_try:14.0e} {eta_needed * 1000:18.0f}")
print()
print(f"  every factor of 10 lost in i0 costs exactly one Tafel slope: {b * 1000:.0f} mV")
print("  the equilibrium potential never moved.")

出力:

Step 1: the thermal voltage scale
  RT/F            = 25.69 mV
  f = nF/(RT)     = 38.924 1/V
  2.303 RT/F      = 59.17 mV/decade   (Nernst slope)

Step 2: both directions, and their sum
  eta (mV)       i_anodic     i_cathodic        i_total   |back/fwd|
--------------------------------------------------------------------
       0.0     1.0000e-03    -1.0000e-03     0.0000e+00    1.000e+00
      10.0     1.2148e-03    -8.2315e-04     3.9170e-04    6.776e-01
      25.0     1.6267e-03    -6.1475e-04     1.0119e-03    3.779e-01
      50.0     2.6461e-03    -3.7791e-04     2.2682e-03    1.428e-01
     100.0     7.0020e-03    -1.4282e-04     6.8591e-03    2.040e-02
     200.0     4.9027e-02    -2.0397e-05     4.9007e-02    4.160e-04
     300.0     3.4329e-01    -2.9130e-06     3.4328e-01    8.486e-06

Step 3: near equilibrium the curve is a straight line
  Rct = RT/(nF i0) = 25.691 ohm cm^2   (for the stand-in i0 above)
  eta =   1 mV: exact 3.8926e-05, linear 3.8924e-05, error  -0.01 %
  eta =   2 mV: exact 7.7867e-05, linear 7.7848e-05, error  -0.03 %
  eta =   5 mV: exact 1.9493e-04, linear 1.9462e-04, error  -0.16 %
  eta =  10 mV: exact 3.9170e-04, linear 3.8924e-04, error  -0.63 %
  eta =  20 mV: exact 7.9828e-04, linear 7.7848e-04, error  -2.48 %
  eta =  50 mV: exact 2.2682e-03, linear 1.9462e-03, error -14.20 %

Step 4: fit log10|i| vs eta on the high-overpotential branch
  fit window            : eta = 150 to 400 mV
  fitted slope          : 8.4547 decades/V
  fitted Tafel slope    : 118.28 mV/decade
  predicted 2.303RT/(anF): 118.33 mV/decade
  intercept -> i0       : 9.9813e-04 A/cm^2 (input was 1.0000e-03)

Step 5: Tafel slope vs transfer coefficient (anodic branch)
  alpha_a    mV/decade
----------------------
     0.20        295.8
     0.30        197.2
     0.40        147.9
     0.50        118.3
     0.60         98.6
     1.00         59.2
     1.50         39.4
     2.00         29.6

Step 6: overpotential needed to reach 10 mA/cm^2
   i0 (A/cm^2)    eta needed (mV)
----------------------------------
         1e-03                118
         1e-05                355
         1e-07                592
         1e-09                828
         1e-11               1065

  every factor of 10 lost in i0 costs exactly one Tafel slope: 118 mV
  the equilibrium potential never moved.

結果の読み方。 重要度の低いものから順に5つ挙げます。

ALPHA を \(0.3\) に変えて再実行してみてください。Step 4のフィット勾配は、予測値 \(84.52\) に対して \(84.49\) mV/decade で返ってきます。これは \(2.303RT/(0.3F)\) ではなく \(2.303RT/(0.7F)\) に一致します。アノード側 の枝は \(\alpha\) ではなく \((1-\alpha)\) に支配されているからです——何世代もの学生をつまずかせてきた符号の約束であり、罠でなくなるように一度は手で導いておく価値があります。

🎯 演習問題

  1. 2通りの平衡。 2つの電極が同じ平衡電位をもちながら、交換電流密度は \(10^4\) 倍だけ異なっているとします。両方の電流–過電圧曲線を同じ軸の上に、次にTafelプロットの上にスケッチしてください。それぞれのプロットでどの特徴が一致し、どれが一致しないかを述べ、2つの触媒を比較するのにどちらの測定を選ぶかについて、それが何を意味するかを説明してください。

  2. 過電圧を診断する。 あるセルを3つの電流で測定し、熱力学値を超える余分な電圧が、電流を2回2倍にするにつれて120、240、480 mVであったとします。どの寄与が支配的で、なぜそう分かりますか。3つの値が120、155、190 mVである場合について、同じ推論を繰り返してください。

  3. 線形領域を、手で。 \(R_{\text{ct}} = RT/(nF i_0)\) を使って、25 °C、\(n = 1\) で \(i_0 = 10^{-6}\) A/cm² の電荷移動抵抗を計算してください。次に、Step 3の誤差表を使って、この方法で測った抵抗を1%の精度で信頼できる最大の過電圧を述べてください。

  4. 大きすぎる勾配。 ある学生がTafel勾配210 mV/decadeを測定し、\(\alpha \approx 0.28\) と結論しました。それを受け入れる前に、見かけのTafel勾配を膨らませる実験的アーティファクトを3つ挙げ、アーティファクトと真の速度論とを区別できる対照実験を1つ記述してください。（3.6節と第4章の両方に関連する材料があります。）

  5. 1デケードの価値。 Step 6の関係を使い、118 mV/decadeの勾配を仮定して、過電圧を400 mVから250 mVに下げるために触媒開発者が \(i_0\) で何桁稼ぐ必要があるかを計算してください。次に、代わりにTafel勾配を下げることで同じ150 mVを得るほうが容易かどうかを論じ、勾配が動くためには物理的に何が変わらなければならないかを述べてください。

まとめ

熱力学は、反応が起こりうるか、そしていくらかかるかを述べます。どれだけ速いかについては、いっさい何も述べません。本章はその隔たりを埋めました。

平衡において電極は静止していません。反応の両方向がフルスピードで走り、厳密に打ち消し合っている のであり、だからこそ往来が界面を絶えず横切っているあいだも計測器はゼロを示すのです。その往来の大きさが 交換電流密度 \(i_0\) ——平衡がどれくらい賑わっているか——であり、良い電極触媒と悪い電極触媒を分かつ唯一の数値です。賑わった平衡は崩しやすく、小さな押し込みで大きな正味電流を生みます。ほとんど死んだ平衡は強い一押しを要求します。同じ反応についての良い触媒と悪い触媒のあいだで、\(i_0\) は何桁にもわたって広がり、表にされた特定の値は表面にあまりに敏感なので、研究室間で転用は効きません。

押し込みそのものが 過電圧 \(\eta = E_{\text{applied}} - E_{\text{eq}}\) であり、それはその全体が無駄になったエネルギーです。それは物理的に異なる3つの罰金へと、3つの異なる電流依存性とともに分解されます。活性化（電流に対して対数的、触媒が触れられる唯一のもの）、濃度（輸送が尽きるまでは無視でき、そこから発散する）、そして オーム \(iR\)（厳密に線形で、化学ではなくセルの性質）です。それらの異なる依存性こそが、よく設計された実験がそれらを見分ける手立てです。

Butler–Volmer式 は \(i_0\) と \(\eta\) を正味電流に結びつけ、方向ごとに1つずつ、2つの指数関数として読めます。その差こそが計測器の見るものです。指数関数が現れるのは、電荷移動が活性化障壁を越えるからであり、その高さは電極電位——化学の他のどの分野もほとんど提供しないつまみ——によって定まります。2つの極限では式は単純化します。平衡近傍では実効的な抵抗器 \(R_{\text{ct}} = RT/(nF i_0)\) へと線形化し、平衡から遠くでは一方の指数関数が死に、生き残ったほうが Tafelプロット、すなわち \(\eta\) に対する \(\log_{10}|i|\) の直線を与えます。その傾きは \(\alpha\) を、切片は \(i_0\) を与えます。

私たちのコードはButler–Volmerのデータを生成し、150 mVより上の枝だけを残し、実験家がそうするとおりに対数空間でフィットしました。予測される \(2.303RT/(\alpha nF) = 118.33\) mV/decade に対して 118.28 mV/decade を復元し、入力した交換電流密度を \(0.2\)% 以内で取り戻しました。逆向きの分極電流が100 mVで順方向の \(2 \times 10^{-2}\)、300 mVで \(8 \times 10^{-6}\) まで落ちること——Tafel領域が自ら生まれる様子——を示し、線形近似が10 mVでは \(0.63\)% の精度で成り立つのに50 mVでは \(14\)% ずれることを示しました。

最後に、分野全体を組織する構造的事実です。触媒は \(i_0\) とTafel勾配を変え、\(E_{\text{eq}}\) を変えることは決してできません。これまでに作られたどんな電極の上でも、水は1.23 V未満では分解しません。私たちのベンチマーク計算はその差を厳密に値付けしました。10 mA/cm² に到達するのに、\(i_0 = 10^{-3}\) A/cm² なら118 mV、\(i_0 = 10^{-11}\) A/cm² なら1065 mVかかり、\(i_0\) の10倍はちょうどTafel勾配1つ分の価値 をもちます。この対数こそが、触媒の進歩が実在し、漸進的で、そして決して無料ではない理由です。

本章のすべては、界面が実際に感じている電位を私たちが知っていると仮定していました。第4章はその仮定を吟味し、それが擁護を必要とすることを見出します。電解液の最初の1ナノメートルに形成される 電気二重層 を眺め、2端子測定が単一の界面を制御できない理由と、したがって 三電極セル が必須である理由を説明し、実際に使われる参照電極と、pHを変える研究がなぜ RHE 基準で報告されるのかを概観します。そのうえで サイクリックボルタモグラム を正しく読みます——軸が何を意味するのか、なぜピークが現れるのか、そして容量性電流とファラデー電流をどう分離するのか。私たちはNumPyで拡散だけからボルタモグラムを組み立て、Fickの法則とNernst式だけから、あの古典的なアヒルの形が何もないところに現れるのを見届けます。

← 第2章: 電極電位と熱力学 第4章: 電気化学界面 →

免責事項

シリーズ目次に戻る 次の章 →