Exact Electrostatic Theory of Fermi-Energy Response at Metallic Interfaces

技術要約：金属界面におけるフェルミエネルギー応答の厳密な静電理論

問題提起
外部摂動に対するフェルミエネルギー（$E_F$）の応答は、仕事関数、電極電位、電気化学的オンセット電位を含む、金属界面における物理的観測量を決定する基本的な要因である。概念的密度汎関数理論（DFT）は、外場ポテンシャル $v_{ext}(\mathbf{r})$ に対する $E_F$ の汎関数微分として定義されるフクイ関数 $f(\mathbf{r})$ を通じて、この現象に対する局所的な線形応答フレームワークを提供する。これは、電子密度 $\rho(\mathbf{r})$ の電子数 $N$ に対する微分とも等価なマックスウェル関係式によって表現される。

しかし、周期的な凝縮相系においてフクイ関数を評価することは、大きな計算上の障壁となる。$N$ に対する微分を計算するには、シミュレーションセル内に有限の電荷（$\pm \Delta N$）を導入する必要がある。周期的なDFTにおいて、これは正味の電荷を生じさせ、静電エネルギーの発散を引き起こす。均質な背景電荷（HBC）を付加するといった標準的な補正は、非物理的な分極アーティファクトを導入し、計算されたフクイ関数を著しく歪ませる。局所状態密度（LDOS）に基づく代替手法も、$E_F$ における電子状態の緩和を適切に考慮できないことが多い。その結果、既存の手法では、固有の近似なしにこの問題に対する厳密な解を提供できていない。

手法
著者らは、電子導体の完全な静電遮蔽特性を利用した厳密な解決策を提案している。核心となる洞察は、十分に厚い金属スラブに対しては、電子の添加または除去によって誘起される電荷の変化は完全に表面に局在し、金属バルクの電荷は中性を保つという点にある。

1. 有限電荷と有限電場の等価性：
   著者らは、系を左側（$z<0$）と右側（$z>0$）の寄与に分割する。完全な遮蔽により、系に加えられた任意の正味の電荷 $\Delta N$ は、両方の表面に等分される（$Q_L = Q_R = \Delta N/2$）。これにより、フクイ関数は、表面電荷を生成するグローバルなメカニズムに依存しない、局所的な表面特性であることが確立される。

2. 有限電場の経路：
   この局所性を活用し、非導電方向への印加電場 $\varepsilon_z$ によって表面電荷の蓄積を独立に誘起できると仮定する。この際、スラブの反対側を無限の電子リザーバーとして扱う。これにより、連鎖律を用いてフクイ関数を以下のように書き換えることができる：
   $$f(\mathbf{r}) = \frac{1}{2} \frac{\partial \rho(\mathbf{r})}{\partial Q} = \lim_{\varepsilon_z \to 0} \frac{\Delta \rho_{+\varepsilon_z}(\mathbf{r}) - \Delta \rho_{-\varepsilon_z}(\mathbf{r})}{4 Q \varepsilon_z}$$
   決定的なのは、この有限電場アプローチがグローバルな電気的中性を維持するため、有限電荷法に伴う非物理的なアーティファクトを排除できる点である。印加された電場は金属を分極させ、電場をバルク内で完全に遮蔽する反対称な電荷分布を生成する。

主要な結果
本理論は、以下の計算テストを通じて検証された：

• 精度と誤差のスケーリング： 本手法は、3つの代表的なPt(111)表面（裸のスラブ、3原子アイランド欠陥を持つスラブ、および吸着CO分子を持つスラブ）に対してテストされた。計算の結果、7層の白金があれば完全な静電遮蔽に十分であることが示された。線形外場下でのフェルミ準位シフト（$\Delta E_F$）を予測する際、本手法は0.1 V/Åまでの電場に対してサブmeVの精度を達成した。誤差は電場の大きさに伴い、厳密に二次的なスケーリング（$\Delta \Delta E_F \sim O(\varepsilon_z^2)$）を示し、選択された電子構造理論（PBE）内での線形応答定式化の厳密性を裏付けた。
• 仕事関数の変化： 本フレームワークは、分子摂動によって誘起される仕事関数の変化（$\Delta W$）を再現することに成功した。具体的には、3 Åの位置にある水分子の様々な配向に対する、欠陥を持つPt(111)表面の応答をモデル化した。本理論は、期待される双極子ジャンプの周囲の $\sim 0.1$ eV の変動を正確に捉え、四重極摂動を含むケースにおいても、参照となるDFT計算と一致した。
• 電極電位と容量： 著者らは、古典的分子動力学軌道から得られた熱平均フクイ関数を用いて、完全なAg(111)/NaF(aq)界面に本手法を適用した。零電荷電位（PZC）に対する平均場電極電位のシフトを計算することで、微分容量－電圧（$C-V$）曲線を作成した。本モデルは、実験的に観察される特徴的な「キャメルシェイプ（ラクダ型）」の微分容量を再現することに成功した。特に、正の電圧における $F^-$ 吸着に関連するピークを再現した。著者らは、純粋な古典モデルは電子のスピルオーバー（溢れ出し）の欠如によりここで失敗することを指摘しており、フクイ関数アプローチは、古典的な電位変動をフィルタリングして正確な電極電位を得ることで、この量子効果を自然に取り入れている。

意義と主張
本論文は、金属界面における静電遮蔽とフェルミエネルギーの変化を結びつける、局所理論の厳密な基礎を確立している。主な意義は以下の通りである：

1. 厳密性： 有限電荷法やHBC法に伴う近似やアーティファクトのない、周期系におけるフクイ関数計算の厳密な解決策を提供する。
2. 理論的統一： 電子の感受率を通じた間接的な計算を必要とする従来の非局所的なDFTの線形応答定式化とは対照的に、導体における静電遮蔽と界面電位ステップの変化との間の直接的な物理的関連性を明らかにした。
3. 実用性： 本手法により、複雑なグランドカノニカル定式化や、特定の現象に対する明示的な溶媒構造を必要とすることなく、標準的な電子構造フレームワーク（半局所的DFTなど）を用いて、フェルミエネルギーの変化、仕事関数の変化、および電極電位を定量的に決定することが可能になる。

著者らは、「厳密」という用語は、金属的な遮蔽を正しく捉える選択された電子構造理論のレベルを超えた追加の近似が存在しないことを指すと明言している。また、完全な遮蔽を前提としているため、絶縁体や、遮蔽が完全に局在化していないナノスケールの金属構造への直接的な適用は制限されるという限界についても認めている。