第2章: 計算水素電極

4回のプロトン共役電子移動、プロトンに値段をつけるトリック、そして過電圧を予測する自由エネルギー図

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

ビデオ講義

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

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

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

第1章は、1つの目標と1つの困難を残して幕を閉じました。目標は 過電圧 — 熱力学的な代価である1.23 Vを超えて、現実のアノードが要求する余分な電圧です。困難は、それを計算するには通電した界面をまるごとシミュレートしなければならないように見える、という点にあります。帯電した金属表面、液体の電解質、そこを泳ぎ回る溶媒和したプロトン、そして量子化学計算コードにとってそもそも自然な量ですらない電極電位。

本章は、それでもOERの計算スクリーニングを実用的なものにしたアイデアを紹介します。2000年代半ばにNørskov、Rossmeisl らによって導入された 計算水素電極（CHE）モデル は、たった1つのうまく選ばれた基準によって、この問題全体を迂回してしまいます。読み終えるころには、ある表面が中間体をどう吸着するかを表す4つの数値を受け取り、自由エネルギー図を描き、そこから理論過電圧をそのまま読み取れるようになっているはずです — しかも、どんな触媒も逃れられないただ1つの熱力学的な制約を、コードが強制してくれる形で。

2.1 4つの段階

表面上のOERは、1回の協奏的な事象として起こるわけではありません。プロトン共役電子移動（PCET） の連なりとして進行します。各段階はプロトン1個と電子1個を取り去り、活性サイトに結合した別の酸素含有種をあとに残していきます。

表面上の清浄な活性サイトを \(\ast\) と書き、そのサイトがヒドロキシ基、酸素原子、ヒドロペルオキソ基を担った状態をそれぞれ \(\ast\text{OH}\)、\(\ast\text{O}\)、\(\ast\text{OOH}\) と書きます。酸性条件の慣例における標準的な4段階機構は次のとおりです。

\[ \text{(1)}\quad \ast + \text{H}_2\text{O} \;\longrightarrow\; \ast\text{OH} + \text{H}^+ + e^- \]

\[ \text{(2)}\quad \ast\text{OH} \;\longrightarrow\; \ast\text{O} + \text{H}^+ + e^- \]

\[ \text{(3)}\quad \ast\text{O} + \text{H}_2\text{O} \;\longrightarrow\; \ast\text{OOH} + \text{H}^+ + e^- \]

\[ \text{(4)}\quad \ast\text{OOH} \;\longrightarrow\; \ast + \text{O}_2 + \text{H}^+ + e^- \]

4つを足し合わせてみましょう。3つの中間体は打ち消し合い — どれもある段階で生成され、次の段階で消費されます — サイト \(\ast\) は再生されます。残るのは

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

であり、これはまさに第1章のOER半反応そのものです。この機構は当該半反応の 分解 であって、そこに何かを付け加えたものではありません。そしてこの事実は、まもなく他のすべてを律する制約へと姿を変えます。

📚 機構を物理的に読む

計算を始める前に、3つの観察をしておく価値があります。

O–O結合が形成されるのは段階3です。 段階1と2は、サイトに吸着した水分子から水素を順に剥ぎ取っていきます。段階3では 2つ目の 水分子を引き込み、その酸素をすでにそこにある酸素と結び合わせて \(\ast\text{OOH}\) 種を作ります。これはOERが成し遂げねばならないもっとも化学的に困難なことであり、この反応がそもそも3つ目の中間体を必要とする理由でもあります。

どの段階も、ちょうどプロトン1個と電子1個を移動させます。 これは記法の偶然ではなく、CHEが利用するまさにその性質です。その数がどの段階でも同一であるため、印加電位は4つの段階すべてに同じように作用します。

サイトは出発点へ戻ります。 段階4のあと、表面は再び清浄になり、次のターンオーバーの準備が整います。触媒とは定義上、戻ってくるもののことです — そしてこのサイクルが閉じていることこそが、4つの \(\Delta G\) の値を1回のターンオーバーの完全な記述として扱うことを許してくれます。

2.2 CHEが解決する問題

たとえば段階2の自由エネルギーを計算するには、生成物の自由エネルギーから反応物のそれを引く必要があります。

\[ \Delta G_2 = \left[G(\ast\text{O}) + \mu(\text{H}^+) + \mu(e^-)\right] - G(\ast\text{OH}) \]

表面に関する項 \(G(\ast\text{O})\) と \(G(\ast\text{OH})\) はごく普通の電子状態計算です — 骨は折れますが、標準的な仕事です。厄介なのは残りの2項のほうです。

プロトン。 \(\mu(\text{H}^+)\) は溶液中のプロトンの化学ポテンシャルです。これを第一原理から得るということは、液体の水に溶媒和したプロトンを、その周りで水素結合ネットワークが再編成される様子ごとシミュレートし、多数の配置について平均を取ることを意味します — それ自体が1つの研究プロジェクトであり、その誤差幅は私たちが分解しようとしている差を飲み込んでしまうでしょう。電子。 \(\mu(e^-)\) は電極中の電子の化学ポテンシャルであり、印加電位に依存します — 標準的な周期系DFT計算には設定する手立てのない量です。

どちらの項も難しく、どちらも省略できません。そして注目してほしいのは、両者が 4つの段階のどれにも 、常に一緒に、常にちょうど1回ずつ現れるということです。

📚 トリック

この最後の観察がアイデアのすべてです。\(\mu(\text{H}^+)\) と \(\mu(e^-)\) が別々に必要になることは決してありません — 必要なのはその 和 だけです。そしてこの和には、それが自明に既知となるような参照系が存在します。

平衡状態にある水素電極を考えましょう。可逆水素電極（RHE）スケールの定義により、RHE基準で \(U = 0\) Vにおいて、反応

\[ \text{H}^+ + e^- \;\rightleftharpoons\; \tfrac{1}{2}\text{H}_2(g) \]

は平衡にあります。平衡にあるとは、両辺の化学ポテンシャルが等しいということです。

\[ \mu(\text{H}^+) + \mu(e^-) \;=\; \tfrac{1}{2}\,\mu(\text{H}_2) \]

そして \(\mu(\text{H}_2)\) は 気相中の水素分子 の自由エネルギー — 量子化学のあらゆる計算のなかでもっとも簡単な部類のものです。溶媒和もなく、電荷もなく、電極もありません。厄介な2項の組は、気相分子の半分に置き換えられたのです。

続いて電位です。 4つの段階はそれぞれ電子1個を電極へ放出します。電極電位を \(U\) まで上げるとその電子のエネルギーは \(eU\) だけ下がるので、どの段階の自由エネルギーも同じだけずれます。

\[ \Delta G_i(U) \;=\; \Delta G_i(U=0) - eU \]

これがCHEの後半であり、驚くほど安上がりです。各段階を \(U = 0\) で 1回だけ 計算すれば、電位依存性のすべては引き算で済みます。再計算も、帯電したスーパーセルも、電解質モデルも要りません。印加電位は算術として入ってくるのです。

📚 実際にΔGへ入るもの

以下で用いる \(\Delta G\) の値は自由エネルギーであって、生の電子エネルギーではありません。そしてこの違いは重要です。実際の研究では、各段階の \(\Delta G\) はいくつかの寄与から組み立てられます。

これらの補正は、各吸着種の振動解析と、気相の参照系についての標準的な熱化学から得られます。本章はそれらの数値を引用しません — 表面、吸着種、汎関数、そして振動の取り扱いに依存するものであり、代表値をでっち上げることは、本シリーズが避けている偽りの精度そのものだからです。持ち帰るべきなのは構造のほうです。以下のコードに与える4つの数値は 補正済みの自由エネルギー であり、生の電子エネルギーを自由エネルギーであるかのように報告する研究は、体裁の問題ではなく本物の誤りを犯していることになります。

2.3 自由エネルギー図

CHEを手にすれば、1回のターンオーバーは \(U = 0\) における4つの数値 \(\Delta G_1, \Delta G_2, \Delta G_3, \Delta G_4\) によって記述されます。その累積和をプロットしたものが 自由エネルギー図 — 4つのPCETを段とする階段です。

ここで、先に予告した制約です。2.1節で見たとおり、4つの段階を足し合わせるとOERの全半反応になります。その全体反応の自由エネルギーは熱力学によって固定されています — 第1章の水分解の自由エネルギーであり、電子1個あたりに直せば1.23 eVです。したがって、いかなる 触媒についても次が成り立ちます。

\[ \sum_{i=1}^{4} \Delta G_i(U=0) \;=\; 4 \times 1.23\ \text{eV} \;=\; 4.92\ \text{eV} \]

これは本章でもっとも重要な一文ですから、注意して読んでください。触媒が支配するのは 4.92 eVが4つの段階へどう分配されるか です。合計を支配することはできません。表面がOERを安上がりにすることはできず、できるのは4回の分割払いの額をより均等にすることだけです。

📚 図から過電圧へ

電位 \(U\) を印加すれば、どの段階も \(eU\) だけ下がります。反応が最後まで下り坂で進むためには — サイクルのどこにも熱力学的な障壁がないためには — すべての段階が \(\Delta G_i - eU \le 0\) を満たさねばならず、これは次を意味します。

\[ eU \;\ge\; \max_i \Delta G_i \]

つまり もっとも大きな1つの段階が、印加すべき電位を決めてしまう のです。これを \(\Delta G_{\max}\) と呼び、それが属する段階を 電位決定段階 と呼びます。サイクル全体が初めて下り坂になる電位は \(\Delta G_{\max}/e\) であり、理論過電圧 は、それが熱力学的な代価をどれだけ超えているかです。

\[ \eta \;=\; \frac{\Delta G_{\max}}{e} - 1.23\ \text{V} \]

合計が固定されていることから、2つの帰結が直ちに従います。

完璧な触媒とは、完璧に均等な触媒のことです。 4つの段階がすべて等しければ、各段階は 4.92/4 = 1.23 eV となり、\(\Delta G_{\max}\) は1.23 eV、\(\eta\) はゼロになります。均等さからの現実のずれはどれも、どこかの段階を1.23 eVより大きくし、\(\eta\) を押し上げます。

1つの段階を改善するのはただではありません。 合計が固定されているため、ある段階の \(\Delta G\) を下げれば、必然的に残りの段階の和が上がります。このモデルのもとでの触媒設計は「すべてを楽にする」ことではなく、再分配の問題 です。第3章では、この再分配が中間体どうしのスケーリング則によってさらに制約されることを示します。火山型プロット（volcano plot）が存在する理由も、\(\eta\) に頑固な下限がある理由も、そこにあります。

2.4 ハンズオン: 自由エネルギー図の計算機

以下のコードはこの枠組み全体を実装しています。\(U = 0\) における4つの段階の自由エネルギーを受け取り、熱力学的な総和則を検査し、任意の印加電位ですべての段階をずらし、累積プロファイルを組み立て、理論過電圧と電位決定段階を報告します。

数値についてひと言 — コードの前に読んでください。 触媒表面に対する本物のDFTはnumpyにできることではありません。したがって以下の各段階のエネルギーは、Catalyst A、B、C という架空の表面に割り当てた、例示のための教育用の値 です。合計が4.92 eVになるように、そして3つの異なる論点を示すように選ばれたものであり、測定値でもDFTの結果でもなく、いかなる実在の材料にも帰属できません。本章のどこにも、実在の触媒をランク付けしている箇所はありません。本物なのは仕組みのほうです — 本物の数値を4つ与えれば、本物の計算をしてくれます。

import numpy as np

# ---------------------------------------------------------------
# 計算水素電極（CHE）による自由エネルギー図。
#
# 活性サイト * の上で進むOERの4段階、酸性条件の慣例では:
#   1)  * + H2O    -> *OH  + H+ + e-
#   2)  *OH        -> *O   + H+ + e-
#   3)  *O + H2O   -> *OOH + H+ + e-
#   4)  *OOH       -> *  + O2 + H+ + e-
#
# どの段階もプロトン1個と電子1個を放出するので、印加電位Uは各段階を
# ちょうどeUだけ引き下げる。
#
# 以下の数値はすべて、架空の触媒に対する例示のための教育用の値で
# ある。測定値でもDFTの結果でもなく、それらはいかなる実在の材料
# にも帰属しない。
# ---------------------------------------------------------------
E0_WATER = 1.23          # V、標準平衡電位（慣例）
N_STEPS = 4
TOTAL_dG = N_STEPS * E0_WATER    # eV、熱力学によって固定される

STEP_LABELS = [
    "*  + H2O -> *OH  + H+ + e-",
    "*OH      -> *O   + H+ + e-",
    "*O + H2O -> *OOH + H+ + e-",
    "*OOH     -> *  + O2 + H+ + e-",
]

CATALYSTS = {
    "Catalyst A": np.array([1.00, 1.60, 1.42, 0.90]),
    "Catalyst B": np.array([0.85, 1.75, 1.35, 0.97]),
    "Catalyst C": np.array([1.20, 1.28, 1.25, 1.19]),
}


def check_sum(dG0, name):
    """4つの段階の和は、どの触媒でも 4 x 1.23 eV でなければならない。"""
    total = dG0.sum()
    ok = np.isclose(total, TOTAL_dG, atol=1e-9)
    print(f"  {name}: sum of steps = {total:.4f} eV "
          f"(required {TOTAL_dG:.2f} eV) -> {'OK' if ok else 'INCONSISTENT'}")
    return ok


def steps_at_U(dG0, U):
    """印加電位Uにおける各段階の自由エネルギー変化（各段階がeUだけ下がる）。"""
    return dG0 - U


def profile_at_U(dG0, U):
    """累積自由エネルギー曲線。清浄なサイトを0として出発する。"""
    return np.concatenate(([0.0], np.cumsum(steps_at_U(dG0, U))))


def theoretical_overpotential(dG0):
    """eta = 最大段階 / e - 1.23 V。最大の段階が電位を決定している。"""
    dG_max = dG0.max()
    limiting = int(dG0.argmax()) + 1
    return dG_max, dG_max - E0_WATER, limiting


print("Constraint check: the sum is fixed by thermodynamics, not by the catalyst")
for name, dG0 in CATALYSTS.items():
    check_sum(dG0, name)
print()

# --- Catalyst A を3つの電位で追いかける --------------------------
name = "Catalyst A"
dG0 = CATALYSTS[name]
dG_max, eta, limiting = theoretical_overpotential(dG0)
U_onset = dG_max          # どの段階も上り坂にならない電位

print(f"{name} (ILLUSTRATIVE teaching values) — step free energies")
for i, (lab, g) in enumerate(zip(STEP_LABELS, dG0), start=1):
    mark = "  <-- potential-limiting" if i == limiting else ""
    print(f"  step {i}  {lab:30s} dG(U=0) = {g:5.2f} eV{mark}")
print(f"  largest step dG_max = {dG_max:.2f} eV  (step {limiting})")
print(f"  theoretical overpotential eta = {dG_max:.2f} - {E0_WATER:.2f} = {eta:.2f} V")
print(f"  onset potential (all steps downhill) U = {U_onset:.2f} V")
print()

print(f"{'':>10} " + " ".join(f"{'step '+str(i):>9}" for i in range(1, N_STEPS + 1))
      + f" {'uphill?':>9}")
print("-" * 62)
for U in [0.00, E0_WATER, U_onset]:
    s = steps_at_U(dG0, U)
    uphill = "yes" if s.max() > 1e-12 else "no"
    print(f"U={U:5.2f} V " + " ".join(f"{v:9.2f}" for v in s) + f" {uphill:>9}")
print()

print("Cumulative free-energy profiles (eV), state 0 = clean site *")
print(f"{'':>10} " + " ".join(f"{'G'+str(i):>9}" for i in range(N_STEPS + 1)))
print("-" * 62)
for U in [0.00, E0_WATER, U_onset]:
    p = profile_at_U(dG0, U)
    print(f"U={U:5.2f} V " + " ".join(f"{v:9.2f}" for v in p))
print()

# --- 例示用の触媒をランク付けする --------------------------------
print("Theoretical overpotentials of the illustrative catalysts")
print(f"{'catalyst':>12} {'dG_max (eV)':>12} {'limiting step':>14} {'eta (V)':>9}")
print("-" * 51)
for cname, g in sorted(CATALYSTS.items(), key=lambda kv: kv[1].max()):
    m, e, lim = theoretical_overpotential(g)
    print(f"{cname:>12} {m:12.2f} {lim:14d} {e:9.2f}")
print()
print("An ideal catalyst would split 4.92 eV into four equal steps of "
      f"{TOTAL_dG/N_STEPS:.2f} eV, giving eta = 0.00 V.")

出力:

Constraint check: the sum is fixed by thermodynamics, not by the catalyst
  Catalyst A: sum of steps = 4.9200 eV (required 4.92 eV) -> OK
  Catalyst B: sum of steps = 4.9200 eV (required 4.92 eV) -> OK
  Catalyst C: sum of steps = 4.9200 eV (required 4.92 eV) -> OK

Catalyst A (ILLUSTRATIVE teaching values) — step free energies
  step 1  *  + H2O -> *OH  + H+ + e-     dG(U=0) =  1.00 eV
  step 2  *OH      -> *O   + H+ + e-     dG(U=0) =  1.60 eV  <-- potential-limiting
  step 3  *O + H2O -> *OOH + H+ + e-     dG(U=0) =  1.42 eV
  step 4  *OOH     -> *  + O2 + H+ + e-  dG(U=0) =  0.90 eV
  largest step dG_max = 1.60 eV  (step 2)
  theoretical overpotential eta = 1.60 - 1.23 = 0.37 V
  onset potential (all steps downhill) U = 1.60 V

              step 1    step 2    step 3    step 4   uphill?
--------------------------------------------------------------
U= 0.00 V      1.00      1.60      1.42      0.90       yes
U= 1.23 V     -0.23      0.37      0.19     -0.33       yes
U= 1.60 V     -0.60      0.00     -0.18     -0.70        no

Cumulative free-energy profiles (eV), state 0 = clean site *
                  G0        G1        G2        G3        G4
--------------------------------------------------------------
U= 0.00 V      0.00      1.00      2.60      4.02      4.92
U= 1.23 V      0.00     -0.23      0.14      0.33      0.00
U= 1.60 V      0.00     -0.60     -0.60     -0.78     -1.48

Theoretical overpotentials of the illustrative catalysts
    catalyst  dG_max (eV)  limiting step   eta (V)
---------------------------------------------------
  Catalyst C         1.28              2      0.05
  Catalyst A         1.60              2      0.37
  Catalyst B         1.75              2      0.52

An ideal catalyst would split 4.92 eV into four equal steps of 1.23 eV, giving eta = 0.00 V.

結果を読み解く。 表に並んだ3つの電位を順に追ってください。それぞれが違うことを語っています。

U = 0 Vでは、すべてが上り坂です。 累積プロファイルは4.92 eVまで単調に上っていきます — 1回のターンオーバーにかかる熱力学的コストの全額であり、まさに総和則が要求するとおりです。これは電源の助けをまったく借りない場合の図であり、語っていることは第1章と同じです。水はひとりでには分解しません。

U = 1.23 Vでは、両端は釣り合いますが途中は釣り合いません。 累積の行を見てください。0.00から始まり、0.00で 終わって います。これもまた総和則の現れです — 4つの段階がそれぞれ1.23 eVずつ引き下げられると、合計4.92 eVからちょうど4.92 eVが差し引かれるのです。熱力学的には、全体反応はいまや平衡にあります。触媒が理想的であれば、話はここで終わっていたはずでした。

そうはなりません。両端のあいだの 経路 が平坦ではないからです。段階2はいまだ+0.37 eVの上り坂であり、プロファイルは−0.23から+0.33まで上ってから下りてきます。この反応が1.23 Vで進めないのは、全体の熱力学がそれを禁じているからではなく、連なりのなかの1つの段階がいまだにそれを禁じているからです。 このたった1行こそが、行儀のよい触媒の上でさえ過電圧が存在する理由です。

U = 1.60 Vでは、最後の上り坂が平らになります。 段階2はちょうどゼロになり、他のすべての段階は負です — サイクル全体が下り坂です。この電位が \(\Delta G_{\max}/e\) であり、それと1.23 Vとの隔たりが理論過電圧、この例示用の表面では0.37 Vです。

ランキングの表は、再分配の問題を見せてくれます。 3つの架空の触媒はいずれも同じ4.92 eVの合計を抱えていますが、その分け方が異なり、過電圧は0.05 Vから0.52 Vまで広がります。Catalyst Cはほぼ理想に近い場合で — 4つの段階は1.20、1.28、1.25、1.19であり、完全に均等な各1.23 eVに近い値です — \(\eta\) もそれに応じて極めて小さくなります。Catalyst Bは1.75 eVの段階を1つ抱えており、そのために0.52 Vを支払っています。残り3つの段階はこの組のなかでもっとも容易であるにもかかわらず、です。このモデルにおいて、触媒の良し悪しは最悪の段階で決まります。他の場所での優秀さは何ももたらしません。

また、私たちの例示用の3つの表面では、たまたま3つとも段階2が電位決定段階になっていることにも注目してください。これは私たちが選んだ教育用の値の性質であって、自然法則ではありません — 吸着のパターンが違えば電位を決定する段階も変わりますし、ある表面でどの段階が電位を決定しているのかを特定することは、実際のCHE研究がもたらすもののなかでも行動に移しやすい成果の1つです。

先へ進む前に、正直な限界を2つ。第一に、これは 純粋な熱力学 です。どの段階も上り坂にならない電位を見つけるだけであり、活性化障壁については何も語りません。熱力学的に下り坂の段階であっても、障壁のせいで遅くなることはありうるのです。第二に、総和則はここでは 検査 として課されていますが、実際の計算では4つの \(\Delta G\) の値は独立したDFT計算から得られるので、その合計は4.92 eVからいくらかずれます。この残差は、手に入る診断のなかでもっとも有用なものの1つです — 残差が大きければ、参照系、補正、あるいは表面モデルのいずれかに誤りがあるということであり、第5章ではこれがどう使われるのかに立ち返ります。

CATALYSTS を編集して、あなた自身の表面を追加してみてください。合計が4.92 eVになる4つの段階を与えて、\(\eta\) をどこまで小さくできるか試してみましょう。次に、1つの段階を極端に小さくしてみて、その制約が残りの段階に何をするのかを眺めてみてください。

🎯 演習問題

  1. サイクルを閉じる。 4つのPCET段階を手で足し合わせ、\(\ast\)、\(\ast\text{OH}\)、\(\ast\text{O}\)、\(\ast\text{OOH}\) がすべて打ち消し合って、OERの半反応が残ることを確かめてください。この打ち消しこそが総和則を強制するのはなぜか、一文で説明してください。

  2. 参照系を言い直す。 RHE基準の0 Vにおいて \(\mu(\text{H}^+) + \mu(e^-) = \tfrac{1}{2}\mu(\text{H}_2)\) が成り立つ理由、そしてCHE法が必要とするのが常に2つの化学ポテンシャルの 和 だけである理由を、自分の言葉で説明してください。

  3. わざと規則を破る。 CATALYSTS を変更し、あるエントリーの合計が4.92 eVではなく5.20 eVになるようにしてください。check_sum がそれを検出することを確認したうえで、そのような合計を出す物理的に実在の計算があったとしたら、それが触媒についてではなく計算について何を教えているのかを説明してください。

  4. 理想的な場合を設計する。 4つの段階の合計が4.92 eVであるという制約を用いて、あらゆる可能な触媒について \(\Delta G_{\max} \ge 1.23\) eV であること、したがって常に \(\eta \ge 0\) であることを証明してください。等号が成り立つのはどのような条件のときでしょうか。

  5. 1つの悪い段階が招くコスト。 Catalyst Cを取り上げ、その段階3を0.30 eV上げ、合計が保たれるように段階1を同じだけ下げてください。コードを使って \(\eta\) を計算し直しましょう。いまやどの段階が電位を決定しているでしょうか。またその結果は、中間体を1つずつ最適化することについて何を語っているでしょうか。

  6. 第1章とつなげる。 Catalyst Aの理論過電圧は0.37 Vです。第1章のエネルギーコスト関数を用いて、このアノード過電圧だけで水素1キログラムあたり何kWh余計にかかることになるのかを計算し、それが現実のペナルティの下限にすぎない理由を明確に述べてください。

まとめ

表面上のOERは、中間体 \(\ast\text{OH}\)、\(\ast\text{O}\)、\(\ast\text{OOH}\) を経由する 4回のプロトン共役電子移動 として進行し、4つの段階はちょうどOERの半反応へと足し合わされます — 中間体は打ち消し合い、活性サイトは再生されるのです。各段階の自由エネルギーを計算するには、溶媒和したプロトンの化学ポテンシャルと、指定された電位における電子の化学ポテンシャルが必要に見えますが、そのどちらも標準的な計算では供給できません。

計算水素電極 は、ただ1つの参照系によってこの2つの障害を取り除きます。厄介な2項が常に一緒に現れること、そして \(\text{H}^+ + e^- \rightleftharpoons \tfrac{1}{2}\text{H}_2(g)\) がRHE基準の0 Vで平衡にあることから、その和は \(\tfrac{1}{2}\mu(\text{H}_2)\) — 気相の水素分子の自由エネルギー — に等しくなります。印加電位はそのあと、純粋な算術として入ってきます。各段階が \(-eU\) だけずれるので、\(U = 0\) での計算1回から、あらゆる電位における図が得られるのです。実際には、こうして組み立てられる \(\Delta G\) の値はDFTの電子エネルギーに零点エネルギー補正とエントロピー補正を組み合わせたものであり、これらは3つの中間体の間で異なるため、打ち消し合ってはくれません。

総和則 は、この手法全体を律する規律です。どんな触媒についても \(\sum \Delta G_i = 4 \times 1.23 = 4.92\) eV。表面が変えられるのは合計ではなく、4回の分割払いへの分け方だけであり、もっとも大きな段階が電位決定段階 となって、\(\eta = \Delta G_{\max}/e - 1.23\) V を与えます。私たちの計算機は、架空のCatalyst A（例示用の段階1.00、1.60、1.42、0.90 eV）を3つの電位で追いかけました。0 Vではどこもかしこも上り坂、1.23 Vでは両端は釣り合うものの段階2はなお+0.37 eVの上り坂であり — 1.23 Vでは足りない理由がまさにこれです — そして1.60 Vで完全に下り坂となり、\(\eta = 0.37\) V が得られました。3つの例示用の表面のあいだで、同じ固定された4.92 eVの合計のもとに \(\eta\) は0.05 Vから0.52 Vまで広がりました — 触媒の良し悪しは最悪の段階で決まり、設計とは削減の問題ではなく 再分配の問題 なのです。

次章では当然の続きを問い、そして居心地の悪い答えに行き着きます。理想的な触媒に必要なのが4つの等しい段階だけだというなら、なぜ誰もそれを作れていないのでしょうか。理由は、\(\ast\text{OH}\)、\(\ast\text{O}\)、\(\ast\text{OOH}\) の結合エネルギーが独立ではないことにあります — それらは スケーリング則 によって互いに縛り合わされており、4.92 eVのある種の分け方を物理的に到達不可能にしてしまうのです。この制約が 火山型プロット を生み、ある系統のなかでの組成のチューニングでは決して突き破れない過電圧の下限を課します。

← 第1章: なぜOERがボトルネックなのか 第3章: スケーリング則と火山型プロット →

免責事項

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