第4章: 電気化学界面

電解液の最初の1ナノメートルで何が起きているのか、それを見るのになぜ電極が3本要るのか、そしてサイクリックボルタモグラムをどう読むか

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

ビデオ講義

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

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

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

第2章と第3章は、無造作に使われながら一度も吟味されなかった1つの言い回しの上に築かれていました。電極の電位 です。熱力学は \(E_{\text{eq}}\) を与え、速度論は \(\eta = E_{\text{applied}} - E_{\text{eq}}\) を与えました。どちらも、電極には明確に定義された電位が存在し、その値を私たちは知っており、それを設定できると仮定していました。これらの仮定のどれもが、聞こえるよりずっと繊細です。

「その電位」は正確にはどこにあるのでしょうか。等電位体である金属の中ではなく、同じくほぼ等電位であるバルク溶液の中でもありません。それは両者のあいだの遷移領域に住んでおり、その遷移領域は驚くほど薄い——事実上すべての電位差がその上で降下する、数ナノメートルの領域です。その領域が 電気二重層 であり、第3章の物理のすべてが実際に起きている場所です。

単一の界面における電位差をどうやって測るのか。測れません。少なくとも2本の配線では測れません。どんな回路にも界面は2つ含まれており、電圧計はその和を報告するからです。これを回避するには第3の電極と、ある特定の仕掛けが必要であり、その仕掛けを理解しているかどうかが、電気化学データを解釈できる人と、ただ集めることしかできない人とを分けます。

そして単一の界面を制御できるようになったとき、それで何をすべきでしょうか。圧倒的に、誰もが最初に走らせる実験は サイクリックボルタンメトリー です。電位を上下に掃引し、電流を記録し、その形を眺める。得られるグラフはこの分野の万能診断です。同時にそれは日常的に過剰解釈されてもおり、本章の最後の部分は、それを誠実に読むことについてのものです——拡散とNernst式だけから、あの古典的なピークをもつ形を生み出すNumPyシミュレーションも含めて。

4.1 最初の1ナノメートル: 電気二重層

金属電極をとり、わずかに負に帯電させます。電解液は何をするでしょうか。

バルク溶液は電気的に中性です——どのカチオンも統計的にアニオンによって釣り合っています。しかし帯電した表面の近くでは、その中性が破れます。負の電極はカチオンを引きつけアニオンを退けるので、金属のすぐ隣の溶液の層は、金属上の負電荷を鏡映する正味の正電荷を帯びます。互いに逆符号の2枚の電荷のシートが、ごく小さな距離を隔てて並ぶ。ゆえに 二重層 です。

この配置の2つの特徴がきわめて重要であり、どちらも幾何から従います。

隔たりは分子スケールです。 対イオンは自らを金属に平らに押しつけることができません。それらは溶媒分子の殻をまとって到来しますし、金属そのものも吸着した溶媒に覆われています。したがってそれらのイオンの中心が最も近づける位置は、表面からおよそ溶媒和殻1つ分だけ外側になります——慣習的に Helmholtz面 と呼ばれる面です。関わる距離は1ナノメートルの何分の一かのオーダーです。これは「小さい」ことの比喩ではありません。実在する、物理的で、分子的な長さです。

したがって電場は途方もないものになります。 1ナノメートルの何分の一かにわたって1ボルトでも降下させれば、電場強度は \(10^9\) V/m のオーダーになります。これは通常の化学のいかなる基準から見ても極端な環境であり、多くのことを説明します。

二重層はコンデンサであり、それには避けようのない帰結があります。 電極電位を変えるたびに、そのコンデンサに電荷を出し入れしなければなりません。その電荷の移動は電流です。それは化学反応が起きているかどうかにかかわらず流れます。4.5節でこの点に戻ります。化学を行う電流からそれを分離することが、この分野の繰り返し現れる実務的問題の1つだからです。

📚 これはどこまで通用するのか? 率直な限界

いま述べた描像——分子直径ほど離れた2枚の電荷のシート——は、最も単純で有用なモデルであり、本シリーズが留まる水準です。実際の二重層は、少なくとも2つの点でより複雑であり、その名を挙げておく価値があります。

この領域については、順次より詳細になっていくモデルの長い歴史がありますが、私たちはそれを飛ばします。追加の構造は本章のどの結論も変えないからです。持ち帰るべきは、分子スケールのコンデンサ、途方もない電場、そして電位が変わるたびに流れる充電電流 です。

4.2 なぜ電極2本では足りないのか

ここに、電気化学を実験的に厄介な科学にしている問題があります。

あなたは1つの界面の電位を知りたい——たとえば、関心のある化学が起きている白金作用電極の電位を。そこで電圧計をつなぎます。一方のリードは白金へ。もう一方のリードもどこかへ行かなければならず、行ける先は溶液の中だけで、それには別の金属片を通す必要があります。

さて、回路の中の界面を数えてみましょう。2つあります。白金–溶液と、第2の金属–溶液です。電圧計は両方にわたる電位差の 和 を報告します。その読みをどちらか一方に帰属させる手立てはなく、どれほど注意深く配線してもこれは変わりません。回路は閉じなければならないからです。

これは会計上の問題よりも悪いのです。反応を駆動するために電流を流せば、その電流は2つの界面を直列に流れ、それぞれが自らのButler–Volmer関係に従って自らの過電圧を発生させます。第2の界面は受動的な参照ではありません。あなたが測ろうとしているあいだにそれは動く のであり、しかもその量は、たまたま引き出している電流に依存します。そして電流は電解液も横切らなければならず、電解液には抵抗があるので、電圧計が報告する何かに \(iR\) 降下——第3章のオーム項——が加わります。

つまり2電極測定は、欲しい界面と、欲しくないうえに電流とともに漂う界面と、オーム項とを与えます。未知数は3つ、数値は1つ。これは精度の問題ではありません。構造的な不可能性です。

三電極セルは、第2の電極が下手にこなしていた2つの仕事を分割することでこれを解決します。

電極 仕事 電流を流すか? 設計上の優先事項
作用極（WE） 研究対象の界面 はい それが試料であり、他のすべてはそれに仕えるために存在する
対極（CE） 電流が流れるよう回路を閉じる はい——そのすべてを 大面積、不活性、そしてWEを汚染しないこと
参照極（RE） 測定の基準となる固定電位を与える いいえ——事実上ゼロ 安定で再現性があり、WEの近くに配置されること

要点は3行目にあります。ポテンショスタットは作用極と参照極のあいだの電位を、入力インピーダンスが非常に高い状態で測るので、事実上 参照電極には電流が流れません。電流が流れないので、参照極は正味の反応を進めず、過電圧を発生させず、自らのオーム降下も寄与しません。それは平衡電位に、動かずに座り続け、基準となる目盛を提供します。

一方、電流はすべて作用極と対極のあいだを流れます。対極は確かに大きく未知の過電圧を発生させます——誰も気にしません。誰もそれを測っていないからです。ポテンショスタットはただ、作用極–参照極間の電圧をあなたが指定した値に保つのに必要なだけ、作用極–対極間の電圧を連続的に調整します。そのフィードバックループこそがポテンショスタット そのもの であり、だからこそ、どれだけ電流が流れているかとは独立に、単一の界面の電位を絶対的な尺度の上で指定できるのです。第2章と第3章のすべてが測定可能なのは、この仕組みのおかげです。

4.3 実務における参照電極

第2章は 標準水素電極（SHE） を使って電位の尺度を定義しました。標準圧の水素ガスが、水素イオン活量が1の酸中の白金黒めっき白金の上を泡立ち、慣習によって \(0\) V と宣言されるものです。定義としては申し分ありませんが、装置としては徹底的に非実用的です。水素ガスの供給を要し、被毒しやすい白金黒表面を注意深く維持する必要があり、その活量条件は何かを秤で量って調製できるものではなく、そして日常的な実験の隣に水素ボンベを置きたい人などいません。

そこで尺度はSHEで定義され、測定は別のものに対して行われます。実務では、

参照電極 何であるか SHEに対するおおよその位置 典型的な用途
Ag/AgCl AgClで被覆した銀線を塩化物溶液に浸したもの 塩化物濃度に依存するが \(+0.2\) V 前後 日常の働き手。小型、頑丈、安価
SCE（飽和カロメル） 水銀／塩化水銀(I) を飽和KCl中に置いたもの やはり \(+0.2\) V 前後 歴史が長く非常に安定。水銀の廃棄が難点
RHE（可逆水素） 同じ水素反応を、あなたの 電解液の その pHで行うもの 構成上、どのpHでも \(0\) V 水分解およびpHを変える研究の標準

中央の列の値が意図的に「前後」と述べていることに注意してください。実際のAg/AgCl参照極の電位は、その内部充填液の塩化物濃度に依存し、よく使われる充填液は何種類かあります。報告値は温度によっても、参照極の充填液とあなたの電解液とのあいだの液絡によっても、わずかにずれます。充填液を述べずに精密な数値を引用するのは偽りの精度であり、本シリーズはそれをしません。信頼できるのは実務上の規則です。Ag/AgClまたはSCEの参照極はSHEより数百ミリボルト正側に座り、報告するどんな電位も、どの参照極に対して測定したかを述べなければならない。 参照極が述べられていない電位はデータではありません。

📚 RHEはなぜ存在し、なぜ水分解の論文はそれを使うのか

可逆水素電極は、きちんと理解しておく価値のあるものです。その目的は微妙であり、実在する問題を解いているからです。

水素反応 \(2\text{H}^+ + 2e^- \rightleftharpoons \text{H}_2\) はプロトンを含むので、Nernst式によりその平衡電位はpHに依存します。25 °C ではそのずれはおなじみの 1 pH単位あたり59 mV ——本シリーズのすべての章にいまや登場したあの \(2.303RT/F\) です——であり、pHが上がるにつれて負側へ動きます。

厄介な帰結はこうです。水の酸化と水の還元は どちらも プロトンを含むので、その平衡電位は 両方とも 同じ1 pH単位あたり59 mVだけずれます。水分解触媒をAg/AgClのような固定参照極に対して測ると、電解液のpHを変えたときに、測定される電位 と 比較対象となる目標とが、同じ方向に同じだけ動きます。本当に関心のある量である過電圧はまったく変わらないのに、生の数値はさまよい、pH間の比較は帳簿づけになってしまいます。

RHEは、尺度のゼロ点を化学とともに動かすことで、その帳簿づけをなくします。RHEは あなた自身の電解液の中の 水素反応なので、関心のある反応と同じ1 pH単位あたり59 mVだけずれます。したがってRHE基準ではpH依存性が相殺されます。水素発生はどのpHでも \(0\) V に、水の酸化は \(1.23\) V に座ります。RHE基準で引用された電位は、良い近似において、一定のオフセットを伴う過電圧です ——異なる電解液で測られた触媒を比較するときに、まさに欲しいものです。

留保が2つ。RHEがpH依存性を相殺するのは、水素発生と同じプロトン／電子比をもつ反応についてだけであり、他の量論をもつ反応はRHE尺度上でもpH依存性を示します。そして論文中の「RHE」は、作用電解液中の本物の水素電極である場合もあれば、測定または仮定されたpHを使ってAg/AgClの読みに事後的に適用された換算である場合もあります。この2つは等しく信頼できるわけではなく、注意深い実験手法の節はどちらを行ったかを述べます。

4.4 サイクリックボルタモグラムを読む

三電極セルがあれば、電位を指令できます。その能力でできる最も情報量の多いこと、そして圧倒的に最もよく行われることは、それを掃引することです。

サイクリックボルタンメトリー は、名前がそのまま述べているとおりのものです。何も起こらない電位から出発し、折り返し電位まで直線的に傾斜させ、直線的に戻し、そのあいだじゅう電流を記録する。電流を電位に対してプロットします——時間に対してではありません。これが新参者を最初に混乱させる点で、時間はループを回る隠れた変数だからです。結果は閉じた曲線になり、還元されてから再酸化される単純な溶存化学種については、電気化学者がアヒルと呼ぶほど特徴的な形をもちます。特徴を1つずつ読んでいきましょう。

出発点の平らな部分。 あなたは、その化学種が現在の酸化状態のまま安定でいられる電位にいます。反応はありません。唯一の電流は、二重層を充電する小さな容量性電流で、4.5節で扱います。

立ち上がり。 電位がその酸化還元対の形式電位に近づくにつれ、Nernst式は表面での酸化体と還元体の比を別の値にせよと要求し始めます。表面は反応することでそれに応じ、電流が流れます。要求が電位に対して指数的なので、立ち上がりは急峻です。

ピーク——人を驚かせる特徴。 電流は最大値を通り過ぎ、そして 下がります。あなたはますます強く押しているというのにです。理由は、あなたが反応物を消費しているからです。形式電位の近くでその表面濃度はゼロへと追い込まれるので、電流は新しい物質がどれだけ速く到達するかによって決まり、そしてそれは、実験が進むにつれて着実に厚くなる枯渇層を横切る拡散によって到達します。層が厚くなれば勾配は緩やかになり、供給は遅くなります。ピークを過ぎると、電位はもはや何も制限していません——輸送が制限しています ——そして電流は枯渇層が広がるにつれて減衰します。これが第3章の濃度過電圧が、形として見えているものです。

逆掃引と、第2のピーク。 電位を折り返します。いま作ったばかりの生成物はまだ電極の近くにあり、拡散して去ってはいません。電位が戻って形式電位を通過すると、それは元に戻され、逆符号のピークを生じます。したがって逆ピークは 生成物が生き延びたかどうか についての直接の報告になります。分解したり、反応したり、逃げ出したりしていれば、ピークは小さいか、まったくありません。化学者はこれを絶えず利用します——サイクリックボルタンメトリーは、ある分子の酸化体または還元体が実験の時間スケールで安定かどうかを知る、最も手早い方法の1つなのです。

📚 ボルタモグラムから読み取られる3つの数値

特徴 何であるか 何に使われるか 率直な留保
ピーク位置 \(E_{pc}\), \(E_{pa}\) カソード側とアノード側の極大の電位 その中点が酸化還元対の形式電位を近似する 拡散係数が等しい、行儀の良い可逆な対についてのみ
ピーク分離 \(\Delta E_p\) \(E_{pa} - E_{pc}\) 電子移動がどれだけ速いかの診断 遅い速度論と同じくらいもっともらしく、未補償抵抗によっても膨らむ——4.6節を参照
ピーク電流 \(i_p\) ピークの高さ 濃度と、掃引速度の平方根とにスケールする 有用なのは掃引速度の指数の部分。高さを拡散係数に換算するには、4.7節のシミュレーションが明示する諸仮定が必要

ピーク分離には特別なコメントが必要です。この分野で最も過剰に読まれる数値だからです。速い（可逆な）1電子対については、\(\Delta E_p\) は特徴的な小さな値をとり、私たちのシミュレーションはそれを第一原理から計算します。遅い（準可逆な）電子移動はそれを広げ、その広がりは掃引速度とともに増大します。速い掃引ほど界面が追いつく時間が少なくなるからです。これは本物で有用な診断です。困るのは、未補償抵抗が \(\Delta E_p\) をきわめてよく似た仕方で広げる ことです——こちらも掃引速度とともに増大します。速い掃引は大きな電流を意味し、大きな電流は大きな \(iR\) を意味するからです。抵抗に手を打たないまま、広い \(\Delta E_p\) を根拠に遅い速度論を報告することは、文献において最もよくある誤りの1つです。

4.5 容量性電流とファラデー電流

あらゆる電気化学測定において、物理的に異なる2つの電流が流れており、それらを見分けることは恒久的な実務上の関心事です。

ファラデー電流 は界面を横切ります。電子が化学種へ、あるいは化学種から移動し、反応が起こり、通った電荷は、変換された物質量とファラデーの法則によって結びつきます——第1章が用いた関係です。これがあなたの望むことを行う電流です。

容量性（非ファラデー）電流 は界面を横切りません。それは4.1節の二重層コンデンサの充放電です。金属の上に電荷が、溶液中に対電荷が蓄積しますが、何も反応せず、何も生成されません。変換ではなく、再配置です。

両者の掃引速度に対するふるまいはまったく異なり、その違いこそが、あなたが両者に対してもつ手がかりです。直線的な電位傾斜で充電されるコンデンサでは、電流は静電容量と電圧の変化率の積です。

\[ i_{\text{cap}} \;=\; C_{\text{dl}} A \, v \]

すなわち 掃引速度に線形 です。一方、拡散律速のファラデーピークでは、ピーク電流は次のようになります。

\[ i_p \;\propto\; \sqrt{v} \]

4.7節のシミュレーションは、これを仮定することなく小数点以下3桁まで確認します。物理的な理由は、速い掃引ほど枯渇層が成長する時間が短くなるからです。層が薄ければ勾配は急になり、その勾配は \(1/\sqrt{Dt}\) のように振る舞います。

その帰結は不愉快で、避けようがありません。欲しい信号と欲しくない背景との比は \(1/\sqrt{v}\) として下がります。短寿命の中間体を捉えるためにまさに必要となる速い掃引は、同時に、容量性電流がファラデーピークを最もよく埋もれさせる条件でもあるのです。4.7節はこのトレードオフに数値を与えます。

4.6 iR補償: あなたが決して印加しなかった電圧

吟味すべき最後の仮定は、第3章が静かに依拠していたものです。あなたが指令した電位が、界面が感じた電位である、という仮定です。

厳密にはそうではありません。作用電極表面と参照電極の先端とのあいだにはある量の電解液が座っており、その電解液には抵抗があります。そこを電流が流れれば、オームの法則が取り分を持っていきます。界面が実際に経験する電位は次のとおりです。

\[ E_{\text{true}} \;=\; E_{\text{applied}} \;-\; i R_u \]

ここで \(R_u\) は 未補償抵抗 ——溶液抵抗のうち、作用電極と参照極先端のあいだに横たわり、したがって測定ループの内側にある部分——です。

この誤差の3つの性質が、それを厄介なものにしています。

電流とともに増大します。 マイクロアンペアでは目に見えません。同じセルで数百ミリアンペアになれば、数百ミリボルトになりうる——あなたが測ろうとしている過電圧の全体に匹敵する大きさです。4.7節がこれを表にしており、その数値は身の引き締まるものです。

歪みの向きは常に同じです。 界面は常に、あなたが指令したよりも小さい駆動力しか感じないので、与えられた印加電位における測定電流は低く出て、Tafel勾配は大きすぎに出て、ボルタンメトリーのピークは本来より大きく分離します。あらゆるアーティファクトが「実際より悪い速度論」の方向を指すのです。

それは化学と同じくらい幾何に依存します。 \(R_u\) は電解液の導電率、作用電極から参照極先端までの距離、そして電流経路の形によって定まります。したがって第一の防衛線は物理的なものです。参照極の先端を作用電極の近くに置くこと ——古典的には、表面近くに配置した細く引き伸ばしたキャピラリで——そして導電性の高い支持電解質を使うことです。ループから排除した1ミリメートルの溶液は、補正しなくて済む抵抗です。

物理的に取り除けないものは、2通りの方法で扱われます。正帰還 は、ポテンショスタットが推定した \(iR_u\) を実時間で出力に足し戻す方法で、速いのですが、かけすぎると制御ループを発振させかねません。測定後補正 は、\(R_u\) を独立に測定し——高周波インピーダンスが通常の手段です——各データ点を事後的に補正する方法で、より安全で、完全に監査可能です。

どちらを使うにせよ、専門家としての義務は同じです。\(R_u\) を述べ、それをどう決定したかを述べ、そのうちどれだけを補償したかを述べること。 その情報なしに報告されたTafel勾配やピーク分離は、化学の測定であると同じくらい、セルの測定でもあるのです。

4.7 ハンズオン: 拡散だけからボルタモグラムを組み立てる

以下のシミュレーションは、サイクリックボルタモグラムを生み出しうる最小限の仮定の集合から、それを組み立てます。1電子還元を受ける溶存化学種が1つ、そして電子移動は 速い と仮定されるので、表面はあらゆる瞬間にNernst式に従います——つまり、コードの中には速度論のモデルがまったく存在しません。残る唯一の物理は半無限拡散であり、1次元格子上の陽的差分法で解かれます。

その禁欲さこそが要点です。もし特徴的なピークをもつ形も、ピーク分離も、掃引速度に対するスケーリングも、すべてが 拡散とNernst式、それだけ から現れるのなら、それらの特徴はそれ以上に風変わりな何かの証拠ではありません。それらは、実際のボルタモグラムが比較されるべき基準線なのです。

コードを読む前に、2つの約束事を。ここでは還元電流を正としてプロットしています。これは使われている2つの流儀のうちの1つです——出会ったボルタモグラムの軸は必ず確認してください。そして輸送のパラメータ（拡散係数 \(10^{-5}\) cm²/s、1 mM溶液、1 cm² の電極）は、水中の小さなイオンについての慣習的な桁のオーダーです。それらは電流の スケール を定めますが、形 についての結論はどれも定めません。形の結論はすべて無次元だからです。

import numpy as np

# ---------------------------------------------------------------
# A minimal cyclic voltammogram, from diffusion alone.
#
# One dissolved species O, one electron, O + e- <-> R.
# Electron transfer is assumed FAST, so the electrode surface obeys
# the Nernst equation at every instant. The only physics left is
# semi-infinite diffusion, solved by explicit finite differences.
#
# Nothing here is fitted to an experiment. The duck shape, the peak
# separation and the sqrt(scan-rate) scaling all fall out of
# diffusion + Nernst and nothing else.
# ---------------------------------------------------------------
F = 96485.0        # C/mol
R_GAS = 8.314      # J/(mol K)
T = 298.15         # K
f = F / (R_GAS * T)

D = 1.0e-5         # cm^2/s   conventional order of magnitude for a small ion in water
C_BULK = 1.0e-6    # mol/cm^3 = 1 mM
AREA = 1.0         # cm^2
E0 = 0.0           # V, formal potential placed at zero for convenience

E_START = 0.4      # V, well positive of E0: nothing happens yet
E_SWITCH = -0.4    # V, well negative of E0: reduction has run its course


def simulate_cv(scan_rate, n_time=40000, lam=0.40):
    """One full cycle. Returns (E, i, cO at the switching potential, dx).

    lam = D dt/dx^2 is the explicit-scheme stability number and must stay <= 0.5.
    """
    span = abs(E_START - E_SWITCH)
    t_total = 2.0 * span / scan_rate
    dt = t_total / n_time
    dx = np.sqrt(D * dt / lam)

    # The depletion layer never grows past a few sqrt(D t); add margin.
    n_x = int(6.0 * np.sqrt(D * t_total) / dx) + 3
    cO = np.full(n_x, C_BULK)
    cR = np.zeros(n_x)

    E = np.empty(n_time)
    i = np.empty(n_time)
    cO_switch = None

    for k in range(n_time):
        t = (k + 1) * dt
        # triangular sweep: down to E_SWITCH, then straight back up
        if t <= t_total / 2.0:
            e = E_START - scan_rate * t
        else:
            e = E_SWITCH + scan_rate * (t - t_total / 2.0)

        # interior: one explicit diffusion step for each species
        cO[1:-1] += lam * (cO[2:] - 2.0 * cO[1:-1] + cO[:-2])
        cR[1:-1] += lam * (cR[2:] - 2.0 * cR[1:-1] + cR[:-2])

        # surface: equal D means O and R trade one-for-one, so cO+cR is
        # continuous across the boundary; Nernst fixes their ratio.
        theta = np.exp(f * (e - E0))          # = cO(0)/cR(0)
        total = cO[1] + cR[1]
        cR[0] = total / (1.0 + theta)
        cO[0] = total - cR[0]

        # surface gradient, three-point (second-order) form
        grad = (-3.0 * cO[0] + 4.0 * cO[1] - cO[2]) / (2.0 * dx)
        i[k] = F * AREA * D * grad            # A, reduction taken as positive
        E[k] = e
        if k == n_time // 2 - 1:              # snapshot at the switching potential
            cO_switch = cO.copy()

    return E, i, cO_switch, dx


# --- 1. one voltammogram, and its landmarks --------------------
v = 0.100                                     # V/s
E, i, cO_switch, dx = simulate_cv(v)
half = len(E) // 2

kc = int(np.argmax(i[:half]))                 # cathodic peak, forward sweep
ka = int(np.argmin(i[half:])) + half          # anodic peak, reverse sweep
ipc, Epc = i[kc], E[kc]
ipa, Epa = i[ka], E[ka]

print("Step 1: landmarks of the simulated voltammogram")
print(f"  scan rate                 : {v * 1000:.0f} mV/s")
print(f"  cathodic peak ipc         : {ipc * 1e6:.2f} uA at Epc = {Epc * 1000:+.1f} mV")
print(f"  anodic  peak ipa          : {ipa * 1e6:.2f} uA at Epa = {Epa * 1000:+.1f} mV")
print(f"  peak separation dEp       : {(Epa - Epc) * 1000:.1f} mV")
print(f"  midpoint (Epa+Epc)/2      : {(Epa + Epc) / 2 * 1000:+.1f} mV  (E0 was set to {E0 * 1000:+.0f} mV)")
print(f"  RT/F                      : {R_GAS * T / F * 1000:.1f} mV")
print(f"  (Epc - E0) / (RT/F)       : {Epc / (R_GAS * T / F):.2f}")
print(f"  raw |ipa / ipc|           : {abs(ipa / ipc):.3f}  (measured from zero, NOT baseline-corrected)")
print()

# --- 2. why the current falls after the peak -------------------
print("Step 2: after the cathodic peak the current decays")
for offset_mV in [0, -25, -50, -100, -200]:
    j = int(np.argmin(np.abs(E[:half] - (Epc + offset_mV / 1000.0))))
    print(f"  E = Epc {offset_mV:+5d} mV : i = {i[j] * 1e6:7.2f} uA "
          f"({100.0 * i[j] / ipc:5.1f} % of peak)")
t_switch = 0.8 / v
depleted = int(np.argmax(cO_switch > 0.99 * C_BULK)) * dx
print(f"  at the switching potential, O is depleted out to ~{depleted * 1e4:.0f} um")
print(f"  for comparison sqrt(D t) at that moment = {np.sqrt(D * t_switch) * 1e4:.0f} um")
print()

# --- 3. does the peak scale with the square root of scan rate? -
print("Step 3: peak current vs scan rate")
print(f"{'v (mV/s)':>10} {'ipc (uA)':>12} {'ipc/v':>12} {'ipc/sqrt(v)':>14} {'dEp (mV)':>10}")
print("-" * 62)
rates = [0.025, 0.050, 0.100, 0.200, 0.400]
peaks = []
for vr in rates:
    Er, ir, _, _ = simulate_cv(vr)
    h = len(Er) // 2
    c = int(np.argmax(ir[:h]))
    a = int(np.argmin(ir[h:])) + h
    peaks.append(ir[c])
    print(f"{vr * 1000:10.0f} {ir[c] * 1e6:12.3f} {ir[c] / vr * 1e6:12.1f} "
          f"{ir[c] / np.sqrt(vr) * 1e6:14.3f} {(Er[a] - Er[c]) * 1000:10.1f}")
peaks = np.array(rates), np.array(peaks)
const = peaks[1] / np.sqrt(peaks[0])
print(f"  spread of ipc/sqrt(v) over a 16x range of scan rate: "
      f"{100.0 * (const.max() - const.min()) / const.mean():.2f} %")
slope = np.polyfit(np.log10(peaks[0]), np.log10(peaks[1]), 1)[0]
print(f"  fitted exponent of ipc vs v (log-log): {slope:.3f}")
print()

# --- 4. capacitive current does not behave that way ------------
# A double layer of capacitance Cdl charged at scan rate v draws
# i_cap = Cdl * A * v -- linear in v, not sqrt(v).
CDL = 20.0e-6      # F/cm^2, textbook order of magnitude for a metal in aqueous electrolyte
print("Step 4: faradaic peak vs capacitive background")
print(f"{'v (mV/s)':>10} {'faradaic (uA)':>15} {'capacitive (uA)':>17} {'ratio':>9}")
print("-" * 55)
for vr, ipk in zip(*peaks):
    icap = CDL * AREA * vr
    print(f"{vr * 1000:10.0f} {ipk * 1e6:15.3f} {icap * 1e6:17.3f} {ipk / icap:9.1f}")
for vr in [10.0, 100.0, 1000.0]:
    ipk = const.mean() * np.sqrt(vr)
    icap = CDL * AREA * vr
    print(f"{vr * 1000:10.0f} {ipk * 1e6:15.3f} {icap * 1e6:17.3f} {ipk / icap:9.1f}")
print()
print("  faradaic ~ sqrt(v), capacitive ~ v, so their RATIO falls as 1/sqrt(v):")
print("  fast scans progressively bury the peak in charging current.")
print()

# --- 5. how much voltage does the solution eat? ----------------
# Ohmic drop is a property of the cell, not of the reaction.
print("Step 5: uncompensated resistance turns current into lost potential")
print(f"{'Ru (ohm)':>10} {'10 uA':>12} {'1 mA':>12} {'100 mA':>12}")
print("-" * 48)
for Ru in [1.0, 10.0, 100.0]:
    a1, a2, a3 = (Ru * cur * 1000.0 for cur in (10e-6, 1e-3, 100e-3))
    print(f"{Ru:10.0f} {a1:9.3f} mV {a2:9.3f} mV {a3:9.3f} mV")

出力:

Step 1: landmarks of the simulated voltammogram
  scan rate                 : 100 mV/s
  cathodic peak ipc         : 268.65 uA at Epc = -28.5 mV
  anodic  peak ipa          : -208.78 uA at Epa = +29.0 mV
  peak separation dEp       : 57.5 mV
  midpoint (Epa+Epc)/2      : +0.2 mV  (E0 was set to +0 mV)
  RT/F                      : 25.7 mV
  (Epc - E0) / (RT/F)       : -1.11
  raw |ipa / ipc|           : 0.777  (measured from zero, NOT baseline-corrected)

Step 2: after the cathodic peak the current decays
  E = Epc    +0 mV : i =  268.65 uA (100.0 % of peak)
  E = Epc   -25 mV : i =  249.03 uA ( 92.7 % of peak)
  E = Epc   -50 mV : i =  214.49 uA ( 79.8 % of peak)
  E = Epc  -100 mV : i =  161.81 uA ( 60.2 % of peak)
  E = Epc  -200 mV : i =  115.99 uA ( 43.2 % of peak)
  at the switching potential, O is depleted out to ~232 um
  for comparison sqrt(D t) at that moment = 89 um

Step 3: peak current vs scan rate
  v (mV/s)     ipc (uA)        ipc/v    ipc/sqrt(v)   dEp (mV)
--------------------------------------------------------------
        25      134.324       5373.0        849.540       57.5
        50      189.963       3799.3        849.540       57.5
       100      268.648       2686.5        849.540       57.5
       200      379.926       1899.6        849.540       57.5
       400      537.296       1343.2        849.540       57.5
  spread of ipc/sqrt(v) over a 16x range of scan rate: 0.00 %
  fitted exponent of ipc vs v (log-log): 0.500

Step 4: faradaic peak vs capacitive background
  v (mV/s)   faradaic (uA)   capacitive (uA)     ratio
-------------------------------------------------------
        25         134.324             0.500     268.6
        50         189.963             1.000     190.0
       100         268.648             2.000     134.3
       200         379.926             4.000      95.0
       400         537.296             8.000      67.2
     10000        2686.481           200.000      13.4
    100000        8495.398          2000.000       4.2
   1000000       26864.808         20000.000       1.3

  faradaic ~ sqrt(v), capacitive ~ v, so their RATIO falls as 1/sqrt(v):
  fast scans progressively bury the peak in charging current.

Step 5: uncompensated resistance turns current into lost potential
  Ru (ohm)        10 uA         1 mA       100 mA
------------------------------------------------
         1     0.010 mV     1.000 mV   100.000 mV
        10     0.100 mV    10.000 mV  1000.000 mV
       100     1.000 mV   100.000 mV 10000.000 mV

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

指摘しておく価値のある、率直な引っかかりが1つあります。本物の罠だからです。印刷された \(|i_{pa}/i_{pc}|\) は \(1.000\) ではなく \(0.777\) です。このシミュレートされた酸化還元対は完全に可逆で完全に安定であるにもかかわらず、です——生成物は分解しえません。モデルにはそれを許す化学が含まれていないからです。この比が \(1\) でないのは、逆ピークが 電流ゼロから 測られている一方、実際にはそれが、まだ減衰しつつある順方向反応の尾の上に乗っているからです。正しく測るには、逆ピークの下で減衰していく順方向電流を外挿し、それをベースラインとして使う必要があります。したがって、素朴にゼロから測って \(0.78\) 前後という比は、安定な 生成物の見え方なのであり、生のプロットから \(1.0\) を期待する人は、生成物が分解していないのに分解していると結論してしまいます。ベースラインは測定の一部 であり、これこそ、隠れた化学をもたないシミュレーションが得意とする類のことです。ずれがあればそれは測定の仕方のアーティファクトに違いありません。他に何でもありえないからです。

このコードで走らせてみる価値のある実験を2つ。E_SWITCH = -0.1 と設定して、還元が進みきる前に折り返してみてください。順方向ピークは \(268.65\) μAのまま手つかずである一方、生のピーク比は \(0.777\) から \(0.572\) へ落ちます——折り返し電位が帰結を伴う実験的選択であることの証であり、生のピーク比に測定条件を添える必要があるもう1つの理由です。次に lam = 0.6 と設定し、陽的スキームの安定条件を破ってみてください。濃度配列は発散します。オーバーフロー警告が出て、それから至るところが nan になります。2つ目はそれ自体として有用な教訓です——ある閾値までは問題なく、それを超えると破滅的に誤る数値手法は、事前に警告してはくれません。そのチェックはあなた自身が行わなければならないのです。

🎯 演習問題

  1. 界面を数える。 2電極セルを描き、電圧計の端子間に電位差が存在するあらゆる場所に印をつけてください。次に三電極版を描き、あなたの図に即して、それらの電位差のうちどれをポテンショスタットが測定から取り除いたのか、そしてどんな仕組みでそうしたのかを正確に説明してください。

  2. どの参照極か? 水の酸化触媒をpH 1、pH 7、pH 13で比較するよう依頼されました。Ag/AgCl基準で報告するかRHE基準で報告するかを決め、1 pH単位あたり59 mVのずれを使ってその選択を定量的に正当化してください。そのうえで、逆の選択のほうが良い状況を1つ述べてください。

  3. 掃引速度シリーズを読む。 同僚が10、40、160 mV/sでピーク電流を測定し、12、24、48 μAを得ました。この化学種は自由に拡散しているのでしょうか、それとも表面に固定されているのでしょうか。4.5節の2つのべき乗則を使って推論を示し、結論を裏づけるにはどんな追加測定が必要かを述べてください。

  4. 広いピーク分離を診断する。 あるボルタモグラムが100 mV/sで \(\Delta E_p = 150\) mV を示し、1 V/sでは320 mVに増大しました。この観察と整合する説明を2つ挙げ、次にそれらを区別する単一の実験を設計してください。（化学を変えずにセルについて何を変えられるかを考えてみてください。）

  5. 高速掃引の壁。 Step 4の出力を使って、この系についてファラデーピークが容量性背景の10分の1にまで落ちる掃引速度を見積もってください。次に、背景が再現性をもつならファラデーピークは容量性背景から回収できるという事実を踏まえ、背景差し引きが測定を救うためには二重層のどんな性質が成り立たなければならないかを論じ、そしてなぜその性質を保証するのが難しいのかを述べてください。

まとめ

本章は、前の2つの章が乗っていた仮定を吟味しました。「電極の電位」が、私たちが知り、制御している、明確に定義されたものである、という仮定です。

電位差はほぼ全体が 電気二重層 ——金属上の電荷を鏡映するように対イオンが蓄積する、分子直径数個分の厚さの領域——に住んでいます。溶媒和したイオンの最接近距離が Helmholtz面 を定義し、1ボルト程度が1ナノメートルの何分の一かにわたって降下するので、そこでの電場は \(10^9\) V/m のオーダーになります。その極端な電場こそが、電極電位を活性化エネルギーのつまみにしているものであり、同じ幾何が界面を効率的な コンデンサ にもしています——だから、何かが反応するかどうかにかかわらず、電位が変わればいつでも電流が流れるのです。

2電極測定は1つの界面を切り出すことができません。回路には界面が2つ含まれ、さらにオーム経路があり、電圧計はその和を報告するからです。三電極セル は仕事を分割することでこれを解決します。対極 は電流のすべてを運び、それ自身の測られない過電圧は無関係です。一方 参照極 は事実上まったく電流を運ばず、したがって平衡電位に固定されたままでいます。ポテンショスタットは、作用極–参照極間の電圧をあなたが求めた値に保ち、そのために必要なだけ作用極–対極間の電圧を調整します。

尺度はSHEによって定義され、測定は実用的な何かに対して行われます。Ag/AgCl と SCE はSHEに対しておよそ \(+0.2\) V に座り、正確な値は充填液と温度に依存します——参照極を添えずに引用された電位はデータではありません。RHE は重要な特別な場合です。水素反応が水分解とまったく同じく 1 pH単位あたり59 mV だけずれるので、RHE基準での報告はpH依存性を相殺し、どのpHでも水素発生を \(0\) V に、水の酸化を \(1.23\) V に置きます。

サイクリックボルタンメトリー は電位を掃引し、電流を記録します。立ち上がりは新しい表面組成を求めるNernst的な要求から来ます。ピーク は駆動力のいかなる変化からでもなく、反応物の枯渇 から来ます。拡散層が厚くなり、電極に供給する勾配が平坦になるからです。逆ピーク は、生成物が元に戻されるまで生き延びたかどうかを報告します。

私たちのシミュレーションは、これらすべてを 拡散とNernst式、それだけ から作り上げました。速度論のモデルもなく、フィットしたパラメータもありません。それはカソードピークを \(-28.5\) mV、形式電位からちょうど \(1.11 \times RT/F\) のところに置き、形式電位をピークの中点として \(0.2\) mV の精度で復元し、\(57.5\) mV のピーク分離を生み、16倍の掃引速度範囲にわたって \(i_p/\sqrt{v}\) が \(0.00\)% の精度で一定であること、そしてフィットした指数が \(0.500\) であることを示し、ピークを200 mVも過ぎたところで電流がなおピークの \(43.2\)% にあり、反応物が232 μmまで枯渇していることを示し、そしてファラデー／容量性比が25 mV/sの269から1000 V/sの \(1.3\) まで落ちること——高速掃引研究の厳しい限界——を実証しました。さらにそれは、完全に安定な生成物について生のピーク電流比 \(0.777\) を生み出しました。ベースラインは測定の一部 であるという警告です。

最後に、未補償抵抗 は、あなたが印加するあらゆる電位から静かに \(iR_u\) を差し引きます。同じ10 Ω が10 μAでは \(0.1\) mV、100 mAでは1000 mVの代価となり、それが引き起こす歪みはすべて同じ方向を指します。実際より悪く見える速度論、です。参照極を近づけ、導電性の高い電解液を使い、残りを補償し、そして自分が行ったことを報告してください。

第5章は、これらすべてを使い切ります。水電解装置の完全な電圧収支——第2章の \(1.23\) V に、第3章の2つの活性化過電圧を加え、本章の \(iR\) を加えたもの——を組み立て、どの項が支配的か、そしてなぜ酸素側が高くつくのかを見ます。多電子経路のせいで速度ではなく選択性が中心問題になる \(\text{CO}_2\) 電解を眺め、電池の充放電を、速度論的な税を伴う熱力学量として読み直します。その途上で OERの計算化学 シリーズにつなぎます。そこでは、誰かが酸素発生触媒を合成する前に、その過電圧を第一原理から計算します。

← 第3章: 速度論 — 過電圧とTafel解析 第5章: 応用 — 電解から電池へ →

免責事項

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