ビデオ講義
このビデオは以下のテキストと同じ内容をカバーしています。お好みの学習形式をお選びください。
🌐 JP | 🇬🇧 EN | Last sync: 2026-08-20
4章分の道具立てが、いま机の上に揃いました。第1章は、これらすべてが電子の受け渡しについての帳簿づけであること、そして同じ2つの電極が2つの向きに読めること——電気を与えてくれるセルと、電気を消費するセル——を確立しました。第2章は熱力学を与えました。電極電位とは電子1個あたりの自由エネルギーであり、\(\Delta G = -nFE\) が2つの言語のあいだを変換し、2つの電極電位の差が、何も起きていないときにセルが示す電圧である、と。第3章は、その「何も起きていない」ことこそが問題なのだと言い張り、何かを起こすために支払う代価として過電圧を導入しました。第4章は、その代価が実際に支払われる界面へ私たちを連れて行き、三電極測定がどうやってそれを切り分けるのかを示しました。
本章は、それらすべてを使い切ります。水電解、CO₂電解、そして電池という3つの技術を取り上げ、そのいずれもをすでに構築した語彙だけで読み解きます。ここに新しい理論はありません。そしてそれこそが要点です。もし前の4章がその務めを果たしたのなら、水素プラントも、炭素循環リアクターも、あなたのノートパソコンの中のセルも、ラベルが違うだけの同じ図であり、それらのあいだの違いは数値と工学の違いであって、原理の違いではないのです。
本章全体を貫く量は、封筒の裏に書ける1本の式です。
\[ E_{\text{cell}} = E_{\text{rev}} + \eta_{\text{anode}} + \left|\eta_{\text{cathode}}\right| + iR_{\text{cell}} \]
第1項は熱力学であり、これと交渉することはできません。残りは、この分野が実際に取り組んでいることのすべてです。5.5節では、この式をコードに変え、4つの項に現実の電圧収支のなかで場所を奪い合わせます。
5.1 水電解: 1.23 Vという床と、その上にあるすべて
水電解は、熱力学が許すことと速度論が請求してくることのあいだの隔たりを示す、電気化学で最も明快な例です。全反応は化学として最も単純な部類です。
\[ 2\text{H}_2\text{O} \;\longrightarrow\; 2\text{H}_2 + \text{O}_2 \]
これを2つの電極での半反応に分けます。酸性では、カソードが 水素発生反応(HER) を、アノードが 酸素発生反応(OER) を担います。
\[ \text{HER (cathode):}\quad 2\text{H}^+ + 2e^- \;\longrightarrow\; \text{H}_2 \]
\[ \text{OER (anode):}\quad 2\text{H}_2\text{O} \;\longrightarrow\; \text{O}_2 + 4\text{H}^+ + 4e^- \]
アルカリ条件では同じ2つの変換がプロトンではなく水酸化物イオンを使って書かれ——HERは \(2\text{H}_2\text{O} + 2e^- \rightarrow \text{H}_2 + 2\text{OH}^-\)、OERは \(4\text{OH}^- \rightarrow \text{O}_2 + 2\text{H}_2\text{O} + 4e^-\) になります——が、電子の数も全反応も同一です。変わるのは表記の約束であって、化学ではありません。
電子の収支を合わせれば、量論はただちに出てきます。OERはO₂ 1分子あたり4電子を放出し、HERはH₂ 1分子あたり2電子を消費するので、このセルは 酸素1モルに対して水素2モル を生成します。これは測定された比ではありません。電子の帳簿づけによって強制された比であり、第1章が電荷を質量に換算するのに用いたのと同じファラデーの法則の議論です。
📚 区別しなければならない2つの電圧
水分解については2つの数値が引用され、それらを取り違えることがこの分野で最もよくある誤りです。
- 可逆電圧、1.23 V。 これは25 °Cにおける水分解反応に \(\Delta G^\circ = -nFE^\circ\) を適用して得られる値です。反応が平衡にある電圧——第3章の言葉で言えば、順方向と逆方向の速度が釣り合い、正味 の電流がゼロになる点です。これより下では、電解は熱力学的に禁じられます。ちょうどこの電圧では、電解は無限にゆっくり進みます。つまり、まったく進まないということです。
- 熱中性電圧、およそ1.48 V。 これは自由エネルギーではなく エンタルピー から来る値です。水を分解する反応は吸熱的であり、生成物が持ち去るエネルギーの一部は、本来なら周囲から熱として供給されるはずのものでした。ちょうど熱中性電圧で動作するセルは、正味の熱を吸収も放出もしません。これより下では、稼働中のセルは外部からの熱を必要とします。これより上では、余剰分がスタック内部で廃熱として現れます。
どちらの数値も、狙って到達すべき目標ではありません。両方とも、自分がどこに立っているかを教えてくれる標識です。1.23 Vで運転するセルは何も生みませんし、1.48 Vを超えて運転するセルはその差分で自分自身を温めています。
現実のセル電圧は本章冒頭に書いた和であり、第1項より後のすべての項は損失です。
- \(\eta_{\text{OER}}\) — アノードにおける活性化過電圧。第3章が導いたとおり、その有効な範囲でターフェル挙動に従います。
- \(\eta_{\text{HER}}\) — カソードにおける同じもの。それ自身の交換電流密度と、それ自身のターフェル勾配を持ちます。
- \(iR_{\text{cell}}\) — 電解質、膜、接点、集電体を横切るオーム損。第4章が測定において補正することを教えたまさにその量ですが、ここでは補正して消し去ることはできません。稼働中のデバイスでは、そのエネルギーは本当に熱として失われているからです。
第3章の濃度過電圧もここに属します。反応物の供給や気泡の除去が電流に追いつかなくなる高速域で現れます。実際の電解装置では、これは大半が物質輸送とガス管理の工学的問題であり、5.5節ではそれらしい数値をでっち上げるのではなく、モデルから除外することにします。
📚 なぜOERがボトルネックなのか
2つの半反応を、式としてではなく プロセス として比べてみると、非対称性は際立っています。
HERは 2電子 を動かします。中間体は少なく——表面に吸着した水素があり、それが結合して放出されるだけ——良い触媒の上では経路は短く、寛容です。電気化学で知られているなかで最も容易な反応の1つです。
OERは 4電子と4プロトン を動かし、しかもどちらの水分子にも存在しなかった 酸素–酸素結合 を形成しなければなりません。標準的な機構像では、表面は結合した中間体の系列——ヒドロキシ種、オキソ種、ヒドロペルオキソ種——を経由し、各段階で電子1個とプロトン1個が引き抜かれます。4回の逐次的な電荷移動、3つの中間体、そして新しい結合の形成。それらすべてが1つのサイトで、すべて同じ電極電位に駆動されて起こります。
この最後の一節こそが、最も深い困難の源です。ただ1つのつまみ——電位——が、4つすべて の段階の駆動力を同時に決めてしまいます。段階3を、段階1を動かさずに調整することはできません。しかも中間体は同じ酸素原子を介して表面に結合するため、触媒を変えるとそれらの結合エネルギーは一緒に動く傾向があります。1つの結合を強めれば、他も強まってしまうのです。この相関こそが、数十年にわたる触媒開発をもってしてもOERの過電圧をゼロへ追い込めなかった理由であり、より良いアノードの探索が、材料をもっと試すという話ではなく、本当に難しい最適化である理由です。
この物語の計算科学の側面——それらの中間体の結合エネルギーがどう計算されるのか、なぜそれらは相関するのか、そしてその相関がある材料クラス全体の過電圧にどのように床を敷いてしまうのか——は、まさに OERの計算化学 の主題です。あのシリーズは、本節が止まるところから始まります。5.5節の収支表を見て、アノードの項が なぜ これほど頑固なのかを知りたくなったなら、次に行くべきはそこです。
実用上の帰結は、携えておく価値のある経験則です。よく作られた水電解装置では、アノード過電圧が電圧収支のなかで制御可能な最大の項であるのが普通です。 5.5節は、これを説明用のモデルで具体化し、各項を改善することにどれだけの価値があるのかを示します。
5.2 CO₂電解: 同じセル、より難しい問い
さて、反応物を変えてみましょう。カソードでプロトンを水素に還元する代わりに、セルに二酸化炭素を供給して それ を還元します。この反応は CO₂還元反応(CO₂RR) と呼ばれ、炭素循環社会への電気化学的な経路です。放出されるはずだったCO₂を取り込み、再生可能電力で駆動し、燃料や化学原料として使える分子として取り戻すのです。
アノードは変わりません。相変わらずOERを担い、相変わらず4電子を必要とし、相変わらずアノード側の収支を支配します。興味深いことはすべてカソードで起こります。
📚 電子数で並ぶ生成物の梯子
CO₂還元は1つの反応ではなく、電子とプロトンを何個供給するかで並ぶ一族です。
- 2電子 で 一酸化炭素 または ギ酸塩 — 最も浅い段であり、最も制御しやすい段でもあります。COは既存の工業化学の原料として直接役立ちます。
- 6〜8電子 で メタノール または メタン — 完全に還元された単一炭素分子であり、燃料として使えます。
- 12電子以上、そして 炭素–炭素結合の形成 を伴えば、エチレン、エタノール、その他の多炭素生成物 — 最も価値の高い標的であり、圧倒的に最も難しい標的でもあります。ここでは触媒は、1つのサイトへ十数回の電荷移動を届けるだけでなく、それをやりながら2つの炭素原子を結合させなければならないからです。
このリストを5.1節と照らして読めば、難しさはすでに見えています。OERは4電子と1本の新しい結合で難しかったのに、多炭素生成物へのCO₂RRは、その3倍の電子と炭素–炭素結合を、1つのサイトで、1つの電位に駆動されて実現せよと要求しているのです。
📚 選択性の問題を、率直に述べる
多電子であることの難しさは本物ですが、決定的な問題ではありません。決定的な問題は 競合 です。
CO₂還元は水の中で起こります。水にはプロトンが含まれています。そしてプロトンは水素に還元されます——5.1節のHER、電気化学で最も容易な反応の1つです——しかも CO₂還元経路の電位に近い 電位で起こるのです。さまざまなCO₂RR生成物の平衡電位は、可逆水素電極スケールで0 V付近に集まります。これはつまり、定義上HER自身の平衡電位の近くだということです。文献によって精確な値は異なりますが、本章はそれらを意図的に表にしません。重要なのは定性的な言明のほうであり、そちらには争いがないからです。すなわち、熱力学はこれらの反応を分離しない のです。
したがって、CO₂還元を駆動できるほど負の電位を印加すれば、同時に水素発生を駆動するには十分すぎるほど負の電位を印加していることになります。2つの反応が、同じ表面で、同じ電子を奪い合います。もしHERが勝てば、あなたの再生可能電力は水から水素を作っただけで、CO₂はそのまま素通りしていったことになります。
これは設計問題全体の枠組みを変えます。水電解では、触媒に 速い ことを求めます。CO₂電解では、触媒に 速く、かつ選り分けられる ことを求めます。そして選り分けのほうが難しい要求です。それは熱力学からは得られないからです。平衡電位が生成物を分離しないのなら、分離できるのは 速度論 だけです。CO₂RRにおける選択性は純粋に速度論的な現象です。表面がどの経路を安定化するのか、どの中間体を結合するのか、どの遷移状態を下げるのか。これは第3章のメッセージの最も強い形です。熱力学は何が許されるかを語る。ここではそれは何もかもを許してしまい、実際に何が起こるかは触媒が決めるのです。
だからこそ、この分野が拠って立つ指標は ファラデー効率 ——通した電荷のうち、狙った生成物になった割合——になります。これは純粋な電子の帳簿づけの量であり、第1章が電荷から生成物への変換を定義したのとまったく同じやり方で定義され、競合を監査可能にします。電流密度が高くてファラデー効率が低いセルとは、高価な水素発生装置です。そして効率と速度に加えて第3の軸があります。安定性 ——数百時間が経ったあとも触媒と電解質が同じ挙動を保つかどうかです。速度、選択性、安定性はトレードオフの面をなしており、工業規模でこの3軸すべてをクリアした材料はまだありません。
専用シリーズを準備中です。 CO₂電解は1つの節ではとても足りません——ガス拡散電極の構成、カーボネート問題、局所pHの効果、生成物の分離、そしてそれらすべてに意味があるかどうかを決める技術経済性。CO₂電解と炭素循環についての続編シリーズを計画しており、そこでは本章の5段落分は内容ではなく出発点の前提として扱われます。ここから持ち帰るべきは問題の形です。同じセル、同じOERアノード、同じ過電圧の語彙。そこに1つの新しい、そして支配的な要求——熱力学が消し去ることを拒む競合相手に対する選択性——が加わる。
5.3 電池を、第2章と第3章から読み直す
電池は別の主題のように感じられます。そうではありません。電池は同じ2電極セルを両方向に動かしたものであり、電解装置について知っていることのすべてが、符号を反転させただけで移植できます。
放電はガルバニ的な過程です。 セル反応は自発的で \(\Delta G < 0\) であり、セルは外部負荷に電流を押し流します。化学が回路に仕事をします。
充電は電解的な過程です。 水を分解させるときとまったく同じように、外部電源で反応を逆向きに押し戻します。回路が化学に仕事をします。
これは第1章のガルバニ/電解の区別であり、電池において特にアノード/カソードの名前が人を混乱させる理由でもあります。定義は 酸化と還元 に結びついていて物理的な端子に結びついてはいないので、放電時にアノードだった電極は充電時にはカソードになるのです。物理的な電極は動いていません。動いたのは反応の向きです。実務家はたいてい、電池の電極を放電時の役割で呼び、その約束を最後まで通すことでこの混乱を回避します——実際的な選択であり、書き下すときには明示的に述べておく価値のある選択です。
📚 セル電圧はどこから来るのか
第2章は、セルの平衡電圧は2つの電極電位の差であると述べました。この言明が、電池の設計原理のすべてです。電位スケール上で遠く離れた2つの酸化還元対を選ぶ のです。第2章は、−0.76 Vの亜鉛対と+0.34 Vの銅対から、1.10 Vを与えるダニエル電池をこのやり方で組み立てました。現代の電池はこれらの対を、単位質量あたりどれだけの電荷を蓄えられるか、どれだけ可逆にサイクルできるか、電解質と接触して生き延びられるか、といった多くの理由から同時に選ばれた材料で置き換えますが、電圧は依然として同じ引き算から来るのであり、他の何からでもありません。
そして第2章のネルンストの式は、携帯電話を使う人なら誰もが目にしていることを説明します。電池が空になるにつれて電圧が垂れる ということです。放電が進むと両電極の組成が変化し、ネルンストの関係は、電極電位がそれらの活量とともにずれると告げます。部分的に放電したセルの開回路電圧は、満充電のセルのそれとは本当に違うのです。これは劣化ではなく熱力学です。
リチウムイオン電池では、化学が特徴的な形をしており、定性的に名を与えておく価値があります。両方の電極が ホスト です。層状あるいは骨格構造をもつ材料で、そこへリチウムイオンが挿入され、そこから引き抜かれ、電極の酸化数がそれを補償するように変化します。リチウムイオンは電解質を通って2つのホストのあいだを行き来し、補償する電子は外部回路を旅します——この設計がしばしば ロッキングチェア型 セルと呼ばれる理由です。本章は、特定の化学系について電極電位も容量もセル電圧も引用しません。それらの数値は厳密な材料、充電状態、測定の約束に依存し、それらを微妙に取り違えることこそ、入門書が読者をミスリードする典型的な仕方だからです。一般化するのは説明の 構造 のほうであり、その構造とは、異なる電位にある2つのホスト、内側のイオン、外側の電子、というものです。
📚 グラフの上に見える過電圧
ここが第3章の投資が最も鮮やかに報われる場所です。あるセルを充電し、次に同じレートで放電して、容量に対して電圧をプロットしてみてください。2本の曲線は重なりません。充電曲線は放電曲線の上に位置し、両者の垂直方向の隔たりが 電圧ギャップ、すなわちヒステリシスです。
このギャップは謎ではありません。同じ3つの損失が、符号を反転して2度現れているだけです。
- 活性化過電圧 — 両電極で、電荷移動の速度が有限であることから生じます。
- 濃度過電圧 — 電解質中のイオン輸送と、ホスト粒子内部の拡散から生じます。高レートでは通常これが支配項であり、急速充電が難しい理由でもあります。
- オーム損 \(iR\) — 電解質、セパレータ、集電体、接点を横切る分です。
放電中、これらの損失は熱力学的な電圧から 引かれます。熱力学が約束したより少ししか取り出せません。充電中、これらは 加わります。より多く入れなければなりません。このギャップに対応するエネルギーは、どこか有用な場所へ行くわけではありません。熱になるのです。電池が激しく使うと温かくなる理由であり、熱管理が後付けではなく一級の工学的問題である理由です。
ここから2つの帰結が直接に導かれ、それらは電池のデータシートが伝えようとしていることの大半を説明します。
- 往復効率は速度論の数値であって、熱力学の数値ではありません。 出ていくエネルギーと入ってくるエネルギーの比は、あのギャップの大きさで決まり、そのギャップは過電圧に \(iR\) を足したものです。速度論を改善すれば効率は改善します。セル化学の電圧にはまったく手を触れずに。
- ギャップは電流とともに大きくなります。 その中のすべての項がそうです。活性化過電圧はターフェル領域で電流の対数とともに増え、オーム損は線形に増え、濃度過電圧は輸送が限界に近づくにつれてさらに速く増えます。これが、セルを速く放電したときに、ゆっくり放電したときより取り出せる有用エネルギーが少なくなる理由であり、「容量」がレートとともに引用されたときにのみ意味をもつ理由です。
まったく同じ文が電解装置を記述することに注意してください。印加しなければならない電圧は、熱力学的な最小値を過電圧と \(iR\) の分だけ超え、その超過分は熱になる。電解装置とは、充電しかしない電池です。 5.5節の収支表は水分解について書かれていますが、その会計は、往復損失がどこへ行ったのかを説明するときに電池技術者が行う会計とまったく同じものです。
5.4 この分野はどこへ向かうのか
3つの方向に名を与えておく価値があります。いずれも本シリーズの何かを直接に延長したものであり、いずれも新参者が実際に貢献できる場所だからです。
電極触媒の設計はデータの問題になりつつあります。 従来の方法——材料を合成し、ターフェルプロットを測定し、反復する——は、1つの研究室が作れる材料の数に縛られています。現代のアプローチは計算された記述子を加えます。ある表面が反応中間体をどれだけ強く結合するかを計算し、それを使って活性を予測し、合成できる数をはるかに超える候補をスクリーニングするのです。それが OERの計算化学 という取り組みを1つの特定の反応に適用したものであり、その上に載るより広いデータ駆動の層——不均一な触媒データセットからどうモデルを構築するのか、そしてそのモデルがどこであなたをミスリードするのか——は 触媒設計へのMI応用 の主題です。どちらも、肝心な限界について誠実です。計算された記述子は予測であり、予測は依然として第4章の三電極測定で検証されなければならないのです。
オペランド測定が、モデルと表面のあいだの隔たりを埋めつつあります。 第4章のあらゆる測定は電気的なもの——電流、電圧、時間——でした。そして電気的な測定は速度を教えてくれますが、表面に何があるのかという同定は教えてくれません。増えつつある オペランド 手法——セルが電位下で、電流を流しながら 稼働している最中に 行われる分光や回折——は、反応中に表面が実際にどう見えるのかを問います。その答えは繰り返し意外なものでした。触媒は稼働条件下で再構成し、溶解して再析出し、酸化数を変えるからです。運転の前と後だけで評価された触媒とは、それが稼働中には一度も取らなかった2つの状態について記述しただけの触媒なのです。
スケールアップはそれ自体が1つの科学であり、工学の後始末ではありません。 実験室の電極は小さく、平坦で、よく撹拌され、電位が一様です。工業用スタックはそのいずれでもありません。大きな電極では電流が不均一に分布し、気泡が活性面積を塞いで抵抗を加え、数百時間を生き延びた材料は年単位で生き延びなければなりません。5.5節の収支表は、なぜこれが些事ではないのかを示します。オーム損が0.1 V余分に増えることは、そのプラントが今後作るであろう生成物1キログラムごとに永久に課される税なのです。
5.5 ハンズオン: セル電圧の収支表を作る
5.1節のすべては算術に変えられます。水電解装置の電圧収支を項ごとに組み立て、どの項が支配的かを見て、電流を増やすにつれてそれらすべてが悪化する様子を眺め、3つの異なる改善策を互いに値付けしてみましょう。
数値を読む前に、この警告を読んでください。 コード中で確立された物理量はちょうど1つ、\(E_{\text{rev}} = 1.23\) V であり、これに定数 \(F\) と \(R\)、および比較に用いるおよそ1.48 Vの熱中性電圧が加わります。以下の速度論パラメータはすべて——交換電流密度も、ターフェル勾配も、セル抵抗も——説明用(ILLUSTRATIVE)です。 それらは桁のオーダーとして妥当であるように、そして収支の構造が見えるように選ばれています。特定の触媒、膜、セルの測定値ではなく、出力中のどの数値も実在するデバイスの性能として引用してはいけません。実在するのは結果の 形 です。どの項が大きいのか、それらが電流とともにどうスケールするのか、そしてどの改善にどれだけの価値があるのか、ということです。
import numpy as np
# ---------------------------------------------------------------
# A cell-voltage budget for water electrolysis.
#
# E_cell(i) = E_rev + eta_OER(i) + eta_HER(i) + i * R_cell
#
# E_rev = 1.23 V is the established thermodynamic value for water
# splitting at 25 C. EVERYTHING ELSE below is an ILLUSTRATIVE model
# parameter, chosen to be plausible in order of magnitude. They are
# NOT measured values of any specific catalyst or any specific cell.
# ---------------------------------------------------------------
F = 96485.0 # C/mol, Faraday constant
R_GAS = 8.314 # J/(mol K)
T = 298.15 # K
E_REV = 1.23 # V, reversible water-splitting voltage at 25 C
E_TN = 1.48 # V, thermoneutral voltage (approx)
# --- 0. The Tafel slope anchor from Chapter 3 -------------------
# b = 2.303 R T / (alpha F). With alpha = 0.5 this is the textbook
# 118 mV/decade. This one IS derived, not assumed.
alpha_ref = 0.5
b_ref = 2.303 * R_GAS * T / (alpha_ref * F)
print("Step 0: Tafel slope anchor (derived, not assumed)")
print(f" b = 2.303 R T / (alpha F), alpha = {alpha_ref}")
print(f" b = {b_ref * 1000:.1f} mV/decade at T = {T:.2f} K")
print()
# --- 1. ILLUSTRATIVE kinetic parameters -------------------------
# Units: exchange current densities in A/cm^2, Tafel slopes in V/dec,
# cell resistance in ohm cm^2.
I0_OER = 1.0e-7 # ILLUSTRATIVE: sluggish 4-electron anode
B_OER = 0.060 # ILLUSTRATIVE: 60 mV/dec
I0_HER = 1.0e-3 # ILLUSTRATIVE: fast 2-electron cathode
B_HER = 0.040 # ILLUSTRATIVE: 40 mV/dec
R_CELL = 0.20 # ILLUSTRATIVE: 0.20 ohm cm^2 (membrane + contacts)
def tafel_overpotential(i, i0, b):
"""Magnitude of the activation overpotential, Tafel form."""
return b * np.log10(i / i0)
# --- 2. The budget at one operating point -----------------------
I_OP = 1.0 # A/cm^2 (= 1000 mA/cm^2), a demanding duty
eta_o = tafel_overpotential(I_OP, I0_OER, B_OER)
eta_h = tafel_overpotential(I_OP, I0_HER, B_HER)
v_ir = I_OP * R_CELL
e_cell = E_REV + eta_o + eta_h + v_ir
print(f"Step 1: voltage budget at i = {I_OP:.2f} A/cm^2 (ILLUSTRATIVE parameters)")
print(f"{'term':<28}{'volts':>10}{'share':>10}")
print("-" * 48)
for name, val in [
("E_rev (thermodynamics)", E_REV),
("eta_OER (anode kinetics)", eta_o),
("eta_HER (cathode kinetics)", eta_h),
("i R_cell (ohmic)", v_ir),
]:
print(f"{name:<28}{val:>10.3f}{val / e_cell * 100:>9.1f}%")
print("-" * 48)
print(f"{'E_cell (total)':<28}{e_cell:>10.3f}{100.0:>9.1f}%")
print()
# --- 3. Efficiency ----------------------------------------------
eff_rev = E_REV / e_cell
eff_tn = E_TN / e_cell
print("Step 2: voltage efficiency")
print(f" vs E_rev = 1.23 V : {eff_rev * 100:.1f}%")
print(f" vs E_tn = 1.48 V : {eff_tn * 100:.1f}%")
print(f" overhead above thermodynamics: {e_cell - E_REV:.3f} V "
f"({(e_cell - E_REV) / E_REV * 100:.1f}% of 1.23 V)")
print()
# --- 4. Sweep the current density -------------------------------
currents = np.array([0.01, 0.05, 0.1, 0.5, 1.0, 2.0])
print("Step 3: how the budget changes with current density")
print(f"{'i (A/cm^2)':>11}{'eta_OER':>10}{'eta_HER':>10}{'iR':>8}"
f"{'E_cell':>9}{'eff(1.23)':>11}")
print("-" * 59)
for i in currents:
eo = tafel_overpotential(i, I0_OER, B_OER)
eh = tafel_overpotential(i, I0_HER, B_HER)
ir = i * R_CELL
ec = E_REV + eo + eh + ir
print(f"{i:>11.2f}{eo:>10.3f}{eh:>10.3f}{ir:>8.3f}{ec:>9.3f}"
f"{E_REV / ec * 100:>10.1f}%")
print()
# --- 5. Which knob is worth turning? ----------------------------
# Three separate one-at-a-time interventions at i = 1 A/cm^2.
print(f"Step 4: value of three interventions at i = {I_OP:.2f} A/cm^2")
base = e_cell
def budget(i0_oer=I0_OER, b_oer=B_OER, i0_her=I0_HER, b_her=B_HER,
r=R_CELL, i=I_OP):
return (E_REV + tafel_overpotential(i, i0_oer, b_oer)
+ tafel_overpotential(i, i0_her, b_her) + i * r)
scenarios = [
("baseline", budget()),
("OER i0 x100 (better anode catalyst)", budget(i0_oer=I0_OER * 100)),
("HER i0 x100 (better cathode catalyst)", budget(i0_her=I0_HER * 100)),
("R_cell halved (thinner membrane)", budget(r=R_CELL / 2)),
]
for name, val in scenarios:
print(f" {name:<40}{val:>7.3f} V "
f"(saves {base - val:+.3f} V, eff {E_REV / val * 100:.1f}%)")
print()
# --- 6. What the extra volts cost, per kilogram of hydrogen -----
# Two electrons per H2 molecule. Energy per mole = 2 F E_cell joules.
M_H2 = 2.016e-3 # kg/mol, molar mass of H2
JOULES_PER_KWH = 3.6e6
def kwh_per_kg(e):
return (2 * F * e / JOULES_PER_KWH) / M_H2
print("Step 5: electrical energy per kilogram of H2 (2 F E_cell per mole)")
for label, e in [("at E_rev = 1.23 V (thermodynamic floor)", E_REV),
("at E_tn = 1.48 V (thermoneutral)", E_TN),
(f"at E_cell = {base:.3f} V (this budget)", base)]:
print(f" {label:<42}{kwh_per_kg(e):>7.1f} kWh/kg")
print(f" overhead of this budget over the floor: "
f"{kwh_per_kg(base) - kwh_per_kg(E_REV):.1f} kWh/kg")
出力:
Step 0: Tafel slope anchor (derived, not assumed)
b = 2.303 R T / (alpha F), alpha = 0.5
b = 118.3 mV/decade at T = 298.15 K
Step 1: voltage budget at i = 1.00 A/cm^2 (ILLUSTRATIVE parameters)
term volts share
------------------------------------------------
E_rev (thermodynamics) 1.230 62.4%
eta_OER (anode kinetics) 0.420 21.3%
eta_HER (cathode kinetics) 0.120 6.1%
i R_cell (ohmic) 0.200 10.2%
------------------------------------------------
E_cell (total) 1.970 100.0%
Step 2: voltage efficiency
vs E_rev = 1.23 V : 62.4%
vs E_tn = 1.48 V : 75.1%
overhead above thermodynamics: 0.740 V (60.2% of 1.23 V)
Step 3: how the budget changes with current density
i (A/cm^2) eta_OER eta_HER iR E_cell eff(1.23)
-----------------------------------------------------------
0.01 0.300 0.040 0.002 1.572 78.2%
0.05 0.342 0.068 0.010 1.650 74.6%
0.10 0.360 0.080 0.020 1.690 72.8%
0.50 0.402 0.108 0.100 1.840 66.9%
1.00 0.420 0.120 0.200 1.970 62.4%
2.00 0.438 0.132 0.400 2.200 55.9%
Step 4: value of three interventions at i = 1.00 A/cm^2
baseline 1.970 V (saves +0.000 V, eff 62.4%)
OER i0 x100 (better anode catalyst) 1.850 V (saves +0.120 V, eff 66.5%)
HER i0 x100 (better cathode catalyst) 1.890 V (saves +0.080 V, eff 65.1%)
R_cell halved (thinner membrane) 1.870 V (saves +0.100 V, eff 65.8%)
Step 5: electrical energy per kilogram of H2 (2 F E_cell per mole)
at E_rev = 1.23 V (thermodynamic floor) 32.7 kWh/kg
at E_tn = 1.48 V (thermoneutral) 39.4 kWh/kg
at E_cell = 1.970 V (this budget) 52.4 kWh/kg
overhead of this budget over the floor: 19.7 kWh/kg
結果の読み方。 重要度の低いものから順に5つ挙げます。
-
ここで導出された速度論の数値は、ターフェルの基準値だけです。 Step 0は第3章の結果を再計算します。25 °Cで移動係数0.5のとき、\(b = 2.303RT/(\alpha F) = 118.3\) mV/decadeです。Step 0より後はすべて 説明用の 勾配を使っています。実際の勾配は機構と材料によって変わるからです。この基準値は、勾配とは何を意味するのか、そしておおよそどれくらいの大きさなのかを思い出させるために置かれています。
-
アノードが最大の損失です。 この収支ではOERの過電圧が0.420 Vを占め、HERの0.120 Vに対して3倍以上です。同じ電流、同じセルの中でのことです。この比は選んだパラメータの偶然ではありません。交換電流密度が4桁違うことから導かれるのであり、それは「4電子と新しいO–O結合 対 2電子と水素分子1個」というモデル上の言い換えです。OERがボトルネックであるという5.1節の主張が、帳簿の上に直接見えています。
-
熱力学が電気料金の大半を占める状態は長くは続きません。 1 A/cm²では可逆電圧が印加電圧の62.4%を占め、残りの37.6%が損失です。2 A/cm²まで押し上げると合計は2.200 Vに達し、電圧効率は55.9%まで落ちます。あらゆる工業用電解装置はこのトレードオフの上で生きています。電流密度を高くすれば設備単位あたりの水素は増えますが、その1キログラムごとの効率は悪くなります。両方が最適になる設定は存在せず、だからこの選択は科学的なものではなく経済的なものになるのです。
-
2種類の損失はスケールの仕方が異なり、それがどちらを重視すべきかを変えます。 活性化過電圧は電流の 対数 とともに増えます。0.01から2 A/cm²まででOERの項は0.300 Vから0.438 Vに上がるだけです。オーム項は 線形 に増え、0.002 Vから0.400 Vへ、実に200倍になります。低電流では触媒が事実上問題のすべてであり、高電流ではオーム損がそれに追いつき、なお登り続けています。これが、デバイスが工業的な電流密度へ向けてスケールするにつれて、この分野の関心が触媒作用からセルとスタックの工学へ移っていく定量的な理由です。
-
改善は見かけほどの価値がなく、それはターフェル勾配のしわざです。 アノード触媒を 100倍 活性にしても0.120 Vしか節約できません。セル抵抗を半分にする——ずっと控えめに聞こえる変更です——と0.100 V節約できます。理由は対数にあります。交換電流密度が2桁改善すれば、過電圧はちょうどターフェル勾配2つ分だけ下がるのであり、60 mV/decの勾配ならそれは0.120 Vです。抄録のなかで「100倍」がどれほど立派に響こうと関係ありません。触媒作用は対数で支払い、工学は線形で支払う。 研究プログラムの予算を組む人は、自分がその2つのどちらを買っているのかを知っておくべきです。
最後のステップは、この全体を物理単位で値付けします。水素分子1個あたり2電子と \(F = 96485\) C/mol から1モルあたりのエネルギーが直接求まり、モル質量がそれをキログラム基準に換算します。熱力学的な床は水素1キログラムあたり32.7 kWh、この説明用のセルは52.4 kWhを必要とします——19.7 kWh/kgの上乗せが、1キログラムごとに、永久に発生するのです。それが、1ボルトの過電圧に工業的な生産速度を掛けたときの代価であり、100ミリボルトが丸め誤差ではなく真面目な成果である理由です。
やってみよう。 B_OER = 0.120 と設定して——説明用の60ではなく、おおよそ118 mV/decの基準値です——再実行してみてください。アノードの項も収支全体も悪化しますが、もっと面白いことに、Step 4における100倍の触媒改善の 価値 が2倍になります。触媒の改善は、そのターフェル勾配が2倍急なとき、2倍の価値を持つのです。次に R_CELL = 0.05 と設定して、2 A/cm²の行がどの項に支配されるかを見てください。2箇所の編集で、電解装置の技術者が触媒化学者と交わす2つの主要な議論を再現したことになります。
🎯 演習問題
-
2つの電圧。 ちょうど1.23 Vで運転している水電解装置がなぜ水素を生成しないのか、そして1.23 Vと熱中性の1.48 Vのあいだで物理的に何が変わるのかを、振り返らずに自分の言葉で説明してください。そのうえで、実際のスタックが1.48 Vより下で運転されているとき、運転者は何を結論すべきかを述べてください。
-
電子を数える。 5.1節の2つの半反応だけを使い、全反応の式を参照せずに、H₂とO₂の2:1のモル比を導いてください。そのうえで、生成するモル数が違うにもかかわらず、なぜ 同じ 電流が両方の電極を流れるのかを説明してください。
-
なぜ選択性は熱力学からは得られないのか。 5.2節は、CO₂還元生成物の平衡電位が水素発生のそれの近くに集まっていると述べています。選択性が熱力学の問題になるためには、代わりに何が成り立っていなければならないかを説明し、したがって触媒が制御しなければならないのは第3章のどの量なのかを精確に述べてください。
-
電池のギャップ。 同じセルを同じレートで動かしたときの充電曲線と放電曲線を描き、ギャップにラベルを付け、それを5.3節の3つの損失項に帰属させてください。そのうえで、同じセルを10倍のレートでサイクルした場合にこのスケッチがどう変わるかを、理由とともに予測し、3つの項のうちどれが最も速く増えると考えるかを述べてください。
-
研究プログラムの予算を組む。 5.5節のコードを使って、1 A/cm²で0.300 Vを節約するために必要な
I0_OERの改善量と、同じ0.300 Vを節約するために必要なR_CELLの低減量を、それぞれ求めてください。そのそれぞれが妥当かどうかを述べ、その比較を使って、固定された研究予算を触媒探索とセル工学のあいだにどう配分するかを論じてください。(オーム側の答えが物理的に到達可能ですらあるのかどうか、そしてそれが何を意味するのかにも触れてください。)
まとめ
水電解は カソードでのHER(H₂ 1分子あたり2電子)と アノードでのOER(O₂ 1分子あたり4電子)に分かれ、電子の帳簿づけだけから2:1の生成比が強制されます。その熱力学的な床は 1.23 V ——何も起きない電圧——であり、およそ1.48 Vの熱中性電圧 は、セルが外部からの熱を必要としなくなり、熱を出し始める境目を示します。稼働中のセルは \(E_{\text{rev}} + \eta_{\text{OER}} + \eta_{\text{HER}} + iR_{\text{cell}}\) で運転され、第1項より後のすべての項は損失です。OERがボトルネックである のは、それが4電子と4プロトンを3つの表面中間体を経由して動かし、O–O結合を形成しなければならず、しかもそのすべてが、あらゆる段階の駆動力を一度に決めてしまう単一の電位に駆動されているからです。
CO₂電解 は同じアノードを保ったまま、カソードでの問いを変えます。その生成物は電子数で並ぶ梯子をなします——2電子でCOとギ酸塩、6〜8電子でメタノールとメタン、12電子以上に炭素–炭素結合を加えて多炭素生成物。決定的な困難は電子数ではなく 競合 です。水素発生がCO₂経路の電位の近くに集まった電位で進むため、熱力学はそれらを分離せず、選択性は触媒の純粋に速度論的な性質になり、ファラデー効率によって監査されます。
電池は同じセルを両方向に読んだもの です——放電はガルバニ的、充電は電解的——だからこそ電極のアノード/カソードのラベルは両者のあいだで入れ替わります。セル電圧は電極電位の差であり、放電時の電圧の垂れは組成の変化に作用するネルンストの式であり、充放電の電圧ギャップは、過電圧と \(iR\) が符号を反転して2度現れたもの です。したがって往復効率は速度論の数値であり、あらゆる損失項が電流とともに増えるため、レートを上げると悪化します。
ハンズオンの収支表 は、明確に説明用とした速度論パラメータでこの構造を定量的にしました。1 A/cm²でアノード過電圧は0.420 V、カソードの0.120 Vに対して大きく、印加電圧は1.970 V、1.23 V基準の電圧効率は62.4%でした。活性化損失は電流とともに対数的に増える一方、オーム損は線形に増えるので、デバイスがスケールするにつれて支配的な問題は触媒作用からセル工学へ移ります。100倍良いアノード触媒が買えたのは0.120 V、セル抵抗を半分にして買えたのは0.100 Vでした。触媒作用は対数で支払い、工学は線形で支払う ——そしてエネルギーで言えば、この説明用のセルは、32.7 kWhという熱力学的な床に対して水素1キログラムあたり52.4 kWhを必要としました。
シリーズの結び: 1つの議論、5つの章
5つの章を振り返れば、それらは1つの形をもつ1つの議論をなしています。
第1章 は会計を確立しました。電気化学とは、2つの半分を物理的に分離した酸化還元であり、そのおかげで電子は自分が所有する導線を通らなければならなくなります——そして通ってしまえば、それを数えることができ、数えることこそがファラデーの法則です。第1章はまた地図も描きました。ガルバニ電池は電気を与え、電解セルは電気を消費し、同じハードウェアが両方をこなすのです。
第2章 は何が 可能か を確立しました。電極電位は電子1個あたりの自由エネルギーであり、\(\Delta G = -nFE\) は化学の通貨と電気の通貨のあいだの為替レートであり、2つの電位の差がセル電圧であり、ネルンストの式は濃度が動いたときにその電圧がどう動くかを語ります。水分解の1.23 Vは、熱力学的な事実としてここで登場しました。
第3章 は何が 速いか を確立しました。平衡は速度ではなく、熱力学が許す反応も知覚できないほど遅く進むことがあります。過電圧はそれを進ませるために支払う代価であり、交換電流密度はある表面がどれだけ割のいい取引を提示しているかを測り、バトラー–フォルマー式は有用な領域でターフェル式になり、そして触媒は熱力学にはいっさい触れずに速度論を変えます。この最後の一文が、本シリーズで最も有用な単一の考えです。
第4章 は何が 測定できるか を確立しました。上記のすべては、電気二重層をもち、電位降下がナノメートル単位で起こり、単一の電極の電位を単独で測定する術がない界面で起きています——ゆえに三電極セル、参照電極、ボルタモグラム、そして界面が何をしたのかと電解質がいくら請求したのかを切り分ける \(iR\) 補正が必要になるのです。
第5章 はそれらを使い切りました。本章のあらゆる技術は、ラベルが違うだけのその4つの考えです。下回れない熱力学的な床、両電極で支払う速度論的な過電圧、野心とともに増えるオームの税、そしてどれがどれかを教えてくれる測定の作法。水電解装置とCO₂リアクターと電池が1つの章で教えられるのは、この記述の水準では それらが同じデバイスだから です。
ここから先へ、3つの道が続いています。OERの過電圧がなぜ改善に抵抗するのか——そして結合エネルギーの計算がそれをどう説明し、どう限界づけるのか——を知りたければ、OERの計算化学 へ進んでください。触媒探索の上に載るデータ駆動の層を、その手法と誠実な失敗様式とともに知りたければ、触媒設計へのMI応用 へ進んでください。そして5.2節を読んで、二酸化炭素を再び有用なものへ変える物語の全体が読みたくなったなら、そのシリーズは準備中であり、本シリーズはその前提となるように書かれています。
免責事項
- 本コンテンツは教育・研究・情報提供のみを目的としており、専門的な助言(法律・会計・技術的保証など)を提供するものではありません。
- 本コンテンツおよび付随するCode examplesは「現状有姿(AS IS)」で提供され、明示または黙示を問わず、商品性、特定目的適合性、権利非侵害、正確性・完全性、動作・安全性等いかなる保証もしません。
- 外部リンク、第三者が提供するデータ・ツール・ライブラリ等の内容・可用性・安全性について、作成者および東北大学は一切の責任を負いません。
- 本コンテンツの利用・実行・解釈により直接的・間接的・付随的・特別・結果的・懲罰的損害が生じた場合でも、適用法で許容される最大限の範囲で、作成者および東北大学は責任を負いません。
- 本コンテンツの内容は、予告なく変更・更新・提供停止されることがあります。
- 本コンテンツの著作権・ライセンスは明記された条件(例: CC BY 4.0)に従います。当該ライセンスは通常、無保証条項を含みます。