第1章: なぜOERがボトルネックなのか

水電解、4電子のアノード、そしてグリーン水素の価格を決める過電圧

📖 Reading Time: 20-25 minutes 📊 Difficulty: Beginner 💻 Code Examples: 0 📝 Exercises: 0

ビデオ講義

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

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

マテリアルズ・インフォマティクス道場 > OERの計算化学 > 第1章

グリーン水素のレシピは単純です。水に電気を入れ、水素を取り出す。この化学は2世紀にわたって理解されてきましたし、工学としても商業的に実用化されており、入力である再生可能エネルギー由来の電力は年々安くなっています。それでもなお、水素1キログラムあたりの電気代は熱力学が要求する値より頑固に高いままであり、その超過分の大半はたった1つの電極で費やされています。

その電極が アノード であり、そこで起きているのが 酸素発生反応（OER） です。本章では、副生成物 をつくるほうの反応がなぜ装置全体を律してしまうのかを説明し、本シリーズの残りが土台とする2つの量——平衡電位と 過電圧——を定義し、実際に走らせられるコードによって、1ミリボルトの上乗せがどれだけのエネルギーを要求するのかを正確に示します。

本シリーズは 触媒材料MI応用 シリーズに焦点を絞った姉妹編です。あちらのシリーズは機械学習とベイズ最適化が触媒の組成空間をどう探索するかを扱いますが、こちらはもう1段下の階層、すなわちそれらのモデルが予測しようとしている 物理量 そのものに踏み込みます。すでにあちらを読まれた方は、本シリーズをその記述子の下にある化学だと考えてください。

1.1 水を分解するには電極が2つ要る

水電解は、2つの別々の電極で起きる2つの半反応によって駆動される、1つの全体反応です。

\[ 2\text{H}_2\text{O} \;\longrightarrow\; 2\text{H}_2 + \text{O}_2 \]

電子は水の中を通っていくわけではありません。一方の電極を離れ、外部回路を通って電源を経由し、もう一方の電極に到着します。回路を閉じるのは、電解質の中をイオンが運ぶ電荷です。半反応を別々に、しかも 酸性条件の表記（可動イオンをプロトンとする）で書き下すと、その非対称性がただちに目に見えるようになります。

カソード では、水素発生反応（HER） が起きます。

\[ 2\text{H}^+ + 2e^- \;\longrightarrow\; \text{H}_2 \]

アノード では、酸素発生反応（OER） が起きます。

\[ 2\text{H}_2\text{O} \;\longrightarrow\; \text{O}_2 + 4\text{H}^+ + 4e^- \]

この2つを2:1の比で足し合わせるとプロトンが相殺され、先ほどの全体反応が残ります。この帳簿づけで失われるものも捏造されるものもありません——しかし、2つの半反応は難しさの点で比較になりません。

📚 2電子と4電子

HERは電子を2個動かし、結合を1本つくります。 2個のプロトンが2個の電子を受け取り、対になってH₂になります。良い表面の上では手順は短いものです。水素原子を1つ吸着させ、次にそれを別の吸着水素と結合させるか、あるいは別のプロトンで攻撃するかです。中間体は本質的に1つ——吸着したH原子——であり、したがって適正にすべき結合エネルギーも本質的に1つです。

OERは電子を4個動かし、O=O二重結合をつくります。 2分子の水から4個のプロトンと4個の電子を引き抜き、もともと別々の分子にあった2個の酸素原子を結合させなければなりません。これを表面上で1回の協奏的な事象としてやってのけることは物理的にありそうにありません。反応は 一連の段階 を経て進み、各段階でプロトン1個と電子1個を移動させながら、O₂が最終的に脱離するまで、表面に結合した中間体の連鎖をたどっていきます。

その帰結は偶発的なものではなく、構造的なものです。結合した移動が1つ増えるごとに、エネルギーを失う場所が1つ増えます。 どの中間体も表面に結合していなければなりません——弱く吸着しすぎれば表面はそれを生成できず、強く吸着しすぎれば手放してくれません。HERに向き合う触媒設計者が調整すべき結合エネルギーは、おおむね1つです。OERに向き合う設計者は、3つの中間体の結合エネルギーを、同じ表面の同じサイトで、同時に すべて適正にしなければなりません。

第3章では、これが聞こえよりもさらに厄介であることを示します。それらの結合エネルギーは互いに独立ではないからです。いまのところは、この数え上げの議論を押さえておいてください。結合した移動が増えれば中間体が増え、中間体が増えれば、たった1つの表面が一度に満たさなければならない制約が増えるのです。

1.2 電位は実際に何をしているのか

電気化学ではなくマテリアルズ・インフォマティクスから来た読者は、電位にデータセットの中の数値として出会うことがしばしばあります。その数値が物理的に何をしているのかに1段落を費やす価値があります。第2章の計算手法の全体が、そこに立脚しているからです。

電極電位とは、電極内の電子のエネルギーを操作するつまみ です。電極をより正にすればその電子のエネルギーは下がり、より負にすれば上がります。したがって、吸着種から電極への電子移動は、そのつまみをどこに設定したかによって上り坂にも下り坂にもなります。

定量的には、電子1個を電位差 \(U\) を横切って動かすと、その自由エネルギーは \(-eU\) だけ変化します。\(n\) 個の電子を移動させる反応では、自由エネルギー変化とセル電圧は次のように結びついています。

\[ \Delta G = -nFE \]

ここで \(F\) はファラデー定数、すなわち電子1モルの電荷です。マイナス符号は、自発的な反応（\(\Delta G\) が負）が正のセル電圧を生む、という慣習を符号化したものです。

📚 1.23 Vはどこから来るのか

この関係を逆向きに読んでください。水分解は自発的では ありません——押してやらなければなりません——ので、その \(\Delta G\) は正であり、必要な電位はそれをちょうど打ち消すものです。

液体の水1モルを標準条件で水素と酸素に分解するには正の標準ギブズ自由エネルギーがかかり、この反応はH₂分子1個あたり2個の電子を動かします。そのエネルギーを \(nF\) で割ると、水分解の標準平衡電位 が得られます。

\[ E^{0} = 1.23\ \text{V} \]

本シリーズは1.23 Vを 標準状態の定義 として扱います——あらゆるOERの議論がそれを基準に測られる、慣習上の参照値です。1.5節のコードは、その定義から \(\Delta G = nFE\) を順方向に走らせ、1.23 Vがキロジュールでいくらに相当するのかを示します。そして結果は、水の生成標準ギブズエネルギーが示すとおりの場所に着地します。これは新しい測定ではなく、整合性のチェックです。

同じ数値の2通りの読み方が、後の章で効いてきます。

📚 過電圧: 余分な押し込み

ちょうど1.23 Vでは、熱力学は水分解を許可します——そして、役に立つ速度では何も起きません。平衡電位が述べるのは反応が 禁止でなくなる 場所であって、それがどれだけ速く進むかではありません。実在の電極に実在の電流を流すには、もっと印加しなければなりません。

\[ \eta \;=\; U_{\text{applied}} - E^{0} \]

量 \(\eta\) が 過電圧 です。実在の電極が要求する、熱力学的な代価を超えた余分な電圧のことです。これは本シリーズ全体の中心的な性能指標です。計算OER研究におけるあらゆる触媒比較、あらゆる火山型プロット、あらゆるスクリーニングの試みは、突き詰めれば \(\eta\) をより小さくする競争にほかなりません。

実用的なセルは、そうした罰則の和を供給しなければなりません。アノードの過電圧、カソードの過電圧、そして電解質とハードウェアにおけるオーム損です。本シリーズがアノードに焦点を絞ることを正当化する事実が2つあります。第1に、OERの寄与は電気化学的な項の中で通常最大であり、その理由はまさに1.1節の4電子の話です。第2に、それは 触媒設計 によって最も直接的に手が届く項です。水の導電率を再設計することはできませんが、表面を再設計することはできるのです。

1.3 なぜ過電圧がお金になるのか

その経済性は、市場予測からではなく、保存則から導かれます。

水素1キログラムあたりに必要な電子の数は 化学量論によって固定されています。ファラデーの法則によれば、H₂ 1モルには電子がちょうど2モル必要であり、H₂ 1キログラムは決まったモル数です。したがって1キログラムあたりに押し込まなければならない 電荷 は、いかなる触媒によっても、いかなる技術者によっても、いかなる補助金によっても交渉の余地がありません。

一方、エネルギーは電荷と電圧の積です。電荷が固定されている以上、1キログラムあたりのエネルギーはセル電圧に正比例 します——そしてセル電圧は \(1.23 + \eta\) ボルトです。ここから、この分野全体を動機づける一文が出てきます。

過電圧の1ボルトは、まるごと1ボルト分の純然たる無駄であり、それはすべての電子1個ごとに、プラントの稼働寿命のあいだずっと支払われ続けます。

これはまた、収穫逓減を予想する人を驚かせるスケーリングの振る舞いも説明します。エネルギーは \(\eta\) について線形なので、過電圧を100 mV下げると、400 mVから300 mVへ下げる場合でも200 mVから100 mVへ下げる場合でも、キログラムあたり 同じ だけのキロワット時が節約されます。曲線上のどこにも、さらなる改良が割に合わなくなる点は存在しません。以下のコードがこれを厳密なものにします。

電力価格は一切述べません。価格は地域、契約、年によって変わり、勝手に1つでっち上げるのは本シリーズが避ける類の偽りの精密さです。1キログラムあたりのエネルギーは誠実に計算できる物理量です。あとは、ご自身の状況が抱える価格を掛けてください。

1.4 計算に何ができるのか

本シリーズの残りを実験ではなく量子化学に費やすことの論拠を、ここに示します。

実験のループは遅く、探索空間は膨大です。 OER触媒の候補を試すとは、それを合成し、セルに組み込み、測定するということです——そして、ありうる酸化物表面、ドーパント、終端、ファセットの空間は組合せ的に巨大です。よく運営されたハイスループット実験室でさえ、そのうちの消え入るほどわずかな割合しかサンプルできません。

計算は、合成の前に決定的な問いを立てることができます。 これを扱いやすくしている洞察は、OER活性が主として 表面が反応中間体をどれだけ強く吸着するか によって支配されている、というものです。それらの結合エネルギーは、まだ作られたことのない表面についても第一原理から計算できます。そこでスクリーニングの問いはこうなります。ある仮想的な表面が与えられたとき、それは4つの段階すべてを容易にするようなパターンで *OH、*O、*OOH を吸着するのか? 第2章は、計算水素電極（CHE）モデルを用いて、この問いを計算可能な数値——理論過電圧——へと変えます。

MIとの接続はここにあります。 記述子がひとたび計算可能になれば、それは予測可能になり、予測可能になれば、機械学習モデルは一度も計算されたことのない数百万の候補をスクリーニングできます。そのパイプライン——訓練データのためのDFT、網羅的な走査のためのモデル、生き残ったものに対する実験——を構築するのが 触媒材料MI応用 シリーズです。本シリーズは、そうしたパイプラインが必要とする、物理的に意味のある目的変数を供給します。誰にも解釈できない記述子で訓練されたモデルは、誰にも信頼できないモデルです。

📚 誠実な留保を、早めに述べておく

以上のすべてには条件がついており、それを後から発見するより、いま述べておくほうがましです。

これから構築するモデルが扱うのは理想化された表面です。平坦で完全で静的なファセットであり、溶媒は近似的に扱われ、電子移動は速度論的にではなく熱力学的に制限されていると仮定されます。実在の触媒にはステップ、キンク、欠陥があり、その表面は 稼働しながら再構成 します。私たちが計算する理論過電圧は、測定された分極曲線の予測ではなく、難しさの熱力学的な下限 です。

それでもなお、これは途方もなく有用です。候補を順位づけし、なぜ材料の一族まるごとが性能不足なのかを説明し、どちらの方向へ動くべきかを教えてくれます。測定の代わりにはなりません。第5章では、モデルが誤導すると分かっている場合も含めて、これらの限界をきちんと検討します。本シリーズは、写真ではなく良い地図を作るものとして読んでください。

1.5 ハンズオン: 過電圧の値段

帳簿を計算してみましょう。以下のコードは1.23 Vという定義とファラデーの法則だけを入力とし、それ以外のすべて——平衡電位の背後にあるエネルギー、水素1キログラムあたりに必要なエネルギー、そしてそれが過電圧とともにどうスケールするか——を導出します。

数値について一言。 以下で走査する過電圧の値は、ありうる範囲をまたぐように選ばれた 例示のための教材用の値 です。測定値ではなく、いかなる材料にも帰属されていません。定数——ファラデー定数、H₂のモル質量、1.23 Vという慣習——は定義であり、それ以外はすべてそこから導出されています。

import numpy as np

# ---------------------------------------------------------------
# 水を分解する熱力学的な代価と、過電圧がその請求書に
# 何を上乗せするのか。
#
# 以下の定数はSIで定義された値かCODATA値であり、加えて
# 慣習が1つだけ入っている。水電解の標準平衡セル電圧
# E0 = 1.23 V である。それ以外はすべて導出値。
# ---------------------------------------------------------------
FARADAY = 96485.33212        # C/mol、SIの定義により厳密値
M_H2 = 2.01588e-3            # kg/mol、H2のモル質量
E0_WATER = 1.23              # V、標準平衡セル電圧（慣習）
N_ELECTRONS_PER_H2 = 2       # 2 H+ + 2 e- -> H2

# --- 1. 1.23 Vはどこから来るのか: dG = n F E --------------------
dG_per_mol_H2 = N_ELECTRONS_PER_H2 * FARADAY * E0_WATER   # H2 1 molあたりのJ
print("Step 1: the 1.23 V convention, read as an energy")
print(f"  dG = n F E0 = {N_ELECTRONS_PER_H2} x {FARADAY:.2f} C/mol x {E0_WATER} V")
print(f"     = {dG_per_mol_H2/1000:.1f} kJ per mol H2 split from liquid water")
print(f"  Per electron transferred: {E0_WATER:.2f} eV")
print(f"  A full OER turnover moves 4 electrons: {4*E0_WATER:.2f} eV")
print()

# --- 2. 水素1キログラムあたりの電気エネルギー -------------------
def energy_per_kg_H2(cell_voltage):
    """与えられたセル電圧でH2を1 kgつくるのに要する電気エネルギー（kWh）。

    使うのはファラデーの法則だけ。H2 1 molには電子2 molが必要なので、
    H2 1 kgあたりの電気量は決まっており、エネルギーはセル電圧に対して
    厳密に線形になる。
    """
    charge_per_kg = N_ELECTRONS_PER_H2 * FARADAY / M_H2   # H2 1 kgあたりのC
    joules = charge_per_kg * cell_voltage                 # H2 1 kgあたりのJ
    return joules / 3.6e6                                 # -> 1 kgあたりのkWh

# --- 3. 例示用の過電圧（教材用の値） ----------------------------
# これらのeta値は、ありうる教材用の範囲をまたぐように選んだもの。
# 測定値ではなく、実在するいかなる材料にも結びついていない。
eta_values = np.array([0.00, 0.10, 0.20, 0.30, 0.40, 0.50])

cell_voltage = E0_WATER + eta_values
voltage_efficiency = E0_WATER / cell_voltage
energy = np.array([energy_per_kg_H2(u) for u in cell_voltage])
excess = energy - energy[0]

print("Step 2: cost of an overpotential (ILLUSTRATIVE eta values)")
print(f"{'eta (V)':>8} {'U_cell (V)':>11} {'efficiency':>11} "
      f"{'kWh per kg H2':>15} {'excess kWh/kg':>15}")
print("-" * 64)
for e, u, f, w, x in zip(eta_values, cell_voltage, voltage_efficiency, energy, excess):
    print(f"{e:8.2f} {u:11.2f} {100*f:10.1f}% {w:15.2f} {x:15.2f}")
print()

# --- 4. スケーリングは厳密に線形である --------------------------
slope = np.diff(energy) / np.diff(eta_values)
print("Step 3: is the cost linear in eta?")
print(f"  d(kWh/kg) / d(eta) for each interval: {np.round(slope, 4)}")
print(f"  spread across intervals: {slope.max() - slope.min():.2e} kWh/kg per V")
print(f"  -> every extra 0.10 V of overpotential costs "
      f"{0.10*slope.mean():.2f} kWh per kg H2, at any starting point")
print()

# --- 5. 過電圧を半分にすると、効率の観点ではどうなるか ----------
for e in [0.50, 0.25]:
    u = E0_WATER + e
    print(f"  eta = {e:.2f} V -> U_cell = {u:.2f} V, "
          f"voltage efficiency {100*E0_WATER/u:.1f}%, "
          f"{energy_per_kg_H2(u):.2f} kWh/kg")

出力:

Step 1: the 1.23 V convention, read as an energy
  dG = n F E0 = 2 x 96485.33 C/mol x 1.23 V
     = 237.4 kJ per mol H2 split from liquid water
  Per electron transferred: 1.23 eV
  A full OER turnover moves 4 electrons: 4.92 eV

Step 2: cost of an overpotential (ILLUSTRATIVE eta values)
 eta (V)  U_cell (V)  efficiency   kWh per kg H2   excess kWh/kg
----------------------------------------------------------------
    0.00        1.23      100.0%           32.71            0.00
    0.10        1.33       92.5%           35.37            2.66
    0.20        1.43       86.0%           38.02            5.32
    0.30        1.53       80.4%           40.68            7.98
    0.40        1.63       75.5%           43.34           10.64
    0.50        1.73       71.1%           46.00           13.30

Step 3: is the cost linear in eta?
  d(kWh/kg) / d(eta) for each interval: [26.5904 26.5904 26.5904 26.5904 26.5904]
  spread across intervals: 1.56e-13 kWh/kg per V
  -> every extra 0.10 V of overpotential costs 2.66 kWh per kg H2, at any starting point

  eta = 0.50 V -> U_cell = 1.73 V, voltage efficiency 71.1%, 46.00 kWh/kg
  eta = 0.25 V -> U_cell = 1.48 V, voltage efficiency 83.1%, 39.35 kWh/kg

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

効率の列について1つ注意があります。1.23 Vという参照は 自由エネルギー の参照なので、ここでの「効率100%」は「熱力学的な最小値に対して電気的な仕事の無駄がゼロ」という意味です——可逆的に運転される実在のセルは、周囲から熱を吸収する必要もあります。この比較は、私たちがやっていること（過電圧の順位づけ）に対しては厳密ですが、電解槽の完全な熱収支ではありません。

例示用の走査を、より細かいグリッド——np.arange(0, 0.61, 0.05)——に置き換えて、傾きが一定という結果が変わらないことを確かめてみてください。その不変性こそが要点です。結論は、どの教材用の値を選んだかに依存しないのです。

🎯 演習問題

  1. 移動を数える。 HERとOERを酸性条件の表記で書き下し、それらを2:1の比で組み合わせるとプロトンと電子がすべて相殺され、2H₂O → 2H₂ + O₂ が残ることを手計算で確かめてください。生成されるO₂分子1個あたり、電子が何個回路を横切るかを述べてください。

  2. つまみを読む。 ある電極が0.20 Vだけより正にされました。吸着種 から その電極 へ 電子を1個移動させる自由エネルギーは、eV単位でどれだけ変化しますか。符号を述べ、それを言葉で説明してください。

  3. 固定された電荷。 ファラデーの法則とH₂のモル質量だけを使って、水素1 kgを生成するのに必要な電荷をクーロン単位で計算してください。なぜいかなる触媒もこの数値を変えられないのか、そしてエネルギーの式に現れる因子のうち触媒が変え られる のはどれかを説明してください。

  4. 節約はどこにあるか。 コードを使って、過電圧を0.45 Vから0.35 Vに下げたときと、0.15 Vから0.05 Vに下げたときの、kWh/kgの節約量を計算してください。両者は同じですか。\(U = 1.23 + \eta\) の代数の観点から、なぜ答えがそうならざるをえなかったのかを説明してください。

  5. 主張を監査する。 あるプレスリリースが、新しい触媒が「効率95%を達成した」と述べています。その数値が意味を持つために答えを得ておかなければならない問いを4つ挙げてください——効率がどの参照電位に対して測られているのか、そしてそれが片方の電極を指すのか完全なセルを指すのか、を含めてください。

まとめ

水電解は2つの半反応に分かれますが、その難しさは等しくありません。カソードの HER は本質的に1つの中間体を経て電子を2個動かします。アノードの OER は電子4個とプロトン4個を、表面に結合した中間体の連鎖を通して動かし、しかも別々の水分子に由来する酸素原子からO=O結合をつくらなければなりません。結合した移動が増えれば中間体が増え、中間体が増えれば、たった1つの表面が一度に適正にしなければならない結合エネルギーが増えます——副生成物をつくる側の電極がボトルネックである理由がこれです。

電極電位とは電子のエネルギーを操作するつまみであり、\(\Delta G = -nFE\) によって反応の自由エネルギーと結びついています。水分解の標準平衡電位1.23 V はこの反応の熱力学的な代価であり、本シリーズではそれを標準状態の定義として採用します。私たちのコードはそこから \(\Delta G = nFE\) を順方向に走らせ、H₂ 1モルあたり237.4 kJを、水の熱化学が示すとおりの場所で回収しました。電子1個あたりで読めば、同じ数値が 1.23 eV となり、4電子のOERターンオーバーは 4.92 eV を説明しきらなければなりません。第2章が明示的に課すことになる制約です。

過電圧 \(\eta = U - E^0\) は実在の電極が要求する余分な押し込みであり、本シリーズの残りにとっての性能指標です。ファラデーの法則が水素1キログラムあたりの電荷を固定するので、エネルギーコストは \(\eta\) について厳密に 線形 です。私たちのコードは、傾きが \(10^{-13}\) の範囲で一定であり、0.10 Vごとにキログラムあたり2.66 kWhであることを見出しました。したがって過電圧は、プラントの寿命のあいだ電子1個ごとに支払われる無駄であり、その改良が割に合わなくなることはありません。計算がその居場所を得るのは、仮想的な表面がOER中間体をどう吸着するかを、誰かがそれを合成する 前に 予測できるからです——ただし、私たちがモデル化する理想化された表面が与えるのは難しさについての熱力学的な限界であって、測定された分極曲線ではない、という誠実な留保つきで。

次章では、そうした予測を可能にする道具を組み立てます。4つのプロトン共役電子移動の段階を明示的に書き下し、計算水素電極 ——溶媒和したプロトンを一度もシミュレートすることなくプロトン・電子対に値段をつけられるようにする仕掛け——を導入し、結合エネルギーの組を1つの数値、すなわち理論過電圧へと変換する自由エネルギー図を構築します。

← シリーズトップ 第2章: 計算水素電極 →

免責事項

シリーズトップへ戻る 次の章へ →