電極表面の第一原理計算
山本雅博
京都大学大学院工学研究科物質エネルギー化学専攻 〒 615-8510

京都市西京区京都大学桂Modified on October 25, 2003 1 : 33 pm
First-principles calculation of electrode surface
Masahiro Yamamoto
Department of Energy and Hydrocarbon Chemistry, Graduate School of Engineering, Kyoto University
Nishikyo-ku, Kyoto 615-8510
Electrode interface phenomena have multiscale nature in space and time, then the multiscale modelfrom quantum to continuum level should be considered. In this article we describe the structure of thePt(111)/dipolar liquid interface and the adsorption of sulfur on electrified Au(111) surface. The structureof the Pt(111)/dipolar liquid interface has been investigated by fully self-consistent combination of the firstprinciples calculation based on quantum mechanics for the metal and the reference hypernetted-chain (RHNC)
theory for the liquid. The electronic density profile for the metal, density and orientational structure of liquid molecules, and electrostatic potential across the interface are discussed in detail. A dense layer of liquid
molecules, which is orientationally ordered, is formed near the metal surface, but this surface-induced structure
extends about only three molecular diameters from the surface. This result is in agreement with the recent
experimental observations. For S adsorption on electrified Au(111) surface the screening of the electron at
metal surface and the change of the adsorption energy of sulfur are discussed briefly.

1

緒言を用いてモンテカルロ法や分子動力学法を用いて分子論的に取り扱う必要があること、３）電気２重層の遮蔽長
量子力学の基礎方程式から出発する第一原理計算は、 さが電解質濃度が小さい場合には µm にまで到達するの経験的パラメータのない計算であり、表面・界面の構造・ で Poisson-Boltzmann 方程式等を用いた連続体モデル
物性を予測することが可能である。第一原理計算による で近似しなくてはならないことがある。それにはこれら金属・半導体表面の構造・吸着子との相互作用等についの理論をシームレスに接続するいわゆるマルチスケールて表面科学の分野で多くの報告がなされてきた1 。しか理論が必要であろう。Fig.1 にマルチスケール理論の概しながら、現代の最先端の理論である第一原理計算を用念図を示す。マルチスケール理論は、現実の材料の物性

いても電極表面での現象を理論的に取り扱うことは未は電子の階層、原子・分子の階層、連続体レベルでの階だ完全になされてないと言って良い。既に、幾つかのグ層すべてを考慮して始めて説明可能であるという観点にループが金属｜水溶液界面の系を、第一原理計算（full

基づいた理論である。それぞれの階層の理論を矛盾無くquantum）で理論的に解析した報告はあるが3–6) 、水溶

繋ぐのは容易ではないが、電気化学系の場合はそれぞれ液中に電解質のない帯電していない電極｜純水界面につの階層で計算された静電ポテンシャルで自己無撞着的にいての報告である。これらの解析では計算量が膨大にな繋ぐのが最も単純な手法である。ただし、溶液側で水分るため、十分に統計的平均が得られているのかも不明で子の電極への吸着やハライドの特異吸着等は電子レベルある。

から解かないと正確には表せない。各階層での理論解析

電極の理論的な取り扱いが困難な理由は、１）電極材の例として、Fig.1 には、電子の階層では Au(111) 面上料である金属・半導体界面の電子、電極と吸着分子の結の硫黄の吸着と金上のメルカプトプロピオン酸 (MPA)

合を量子力学的に取り扱う必要があること、２）室温での吸着の第一原理計算による電子分布図、原子・分子のの溶媒を扱うために、電極近傍の溶媒分子、電解質イオ

階層では末端カルボキシル基が解離した MPA 自己組織

ンを量子力学的にあるいは経験的な分子間ポテンシャル

化単分子膜上のモンテカルロシミュレーションによる電1 ただし、表面系では最近指摘されたように Pt(111) 面上の一酸化炭素吸着サイトが実験結果と第一原理計算では一致しない問題がある1, 2) 。密度汎関数理論で電子間の交換・相関相互作用を近似しているところに根本の原因があるようだが、結合次数が変わるような系での計算には注意が必要であろう。

気２重層の解析の模式図、連続体レベルでの階層では、

Poisson-Boltzmann 方程式による末端カルボキシル基が解離した MPA 自己組織化単分子膜上の電気２重層の電1

似 (LDA) である Hedin-Lundqvist 型を用いた。18) 原子位分布の断面図を示した。

に局在化した d 電子をもつ白金を取り扱うために、15電極表面での酸化還元反応についてはバルク中でのイ
7)

Ry までの運動エネルギーをもつ平面波に加えて原子にオンの溶媒和の問題（酸化還元電位の問題 ）、イオンと電極間の電子移動の問題（電子移動速度の問題 ）等

局在する波動関数の Bloch 和の混合基底系を用いた19) 。

さらに複雑な問題を含んでいるため第一原理計算による表面の計算では、逆格子空間内の irreducible Brollouin

取り扱いはほとんどなされていない。

ゾーン内の１８個の点を用いた。

8)

電極界面の構造をマルチスケール理論で扱った例としては、ジェリウムモデルと積分方程式の結合理論9–12) 、水

3分子に分子擬ポテンシャル2 を用いた Halley-Price らの報
告,8)Pt(111) 面の第一原理計算と RHNC(Reference

hypernetted-chain) 積分方程式論を自己無撞着に結合させた我々の報告,14, 15) ３次元 RISM(Reference interac-

金属-液体界面の構造：積分方程式と第一原理計算の結合の自己無撞
着解
近年、金属｜液体界面の構造は、種々の分光法によりtion site model ) 理論と Cu(100) 面の第一原理計算を組

その微視的な構造が明らかになっている20–26) 。これらみ合わせた Kovalenko-平田の報告16, 17) がある。本文での結果を理解するために、第一原理計算のような信頼性

は、我々の研究結果と帯電した Au(111) 面上と硫黄原子のある理論モデルを構築することは重要である。ここでの吸着相互作用について報告する。

は、Pt(111) 面の第一原理計算と双極子液体の積分方程式論の結合の計算結果を述べる。
界面から無限に広がる溶液系を考慮できることに利点があるため、我々は溶液側に積分方程式を用いた。通常のcomputer simulation と異なり、積分方程式を用いた方法
は、系の有限サイズの効果や不十分なモンテカルロ・分子
動力学ステップ等の問題を考慮する必要がない9, 11, 12) 。
さらには、有限濃度での溶媒（例えば水）中の電解質イオンを取り入れることも容易に可能である10) 。Bérard ら
は、積分方程式論のひとつである reference hypernetted-

chain(RHNC) 理論を溶液側に適用し、金属表面と接している水溶液電解質と極性溶液に対して自己無撞着な解

Figure 1: Schematic diagram of multiuscale theory for elec- を与えた9–12) 。溶液に対しては分子モデルを採用した精
trode surface.

密な統計力学理論が用いられているのに対し，金属に対しては金属イオンを一様な正電荷バックグラウンドとみなす jellium モデルを用いるのにとどまっていることに

2

計算方法

問題があった。jellium モデルは，rs (Wigner-Seitz 半径)
値が約 2.5 以下で破綻する27) ことが知られている。しか固体表面の計算は、大別して LCAO(Linear Combina-

し，rs 値が小さい金属には，白金（rs ' 1.4 ）などの実

tion of Atomic Orbital) 基底系での量子化学クラスター

用的に重要なものが多い。本研究は，金属側を経験的パ計算と、平面波基底系を用いたスラブ周期格子系に対す

ラメータのない第一原理計算で行い，金属側・溶液側を静る第一原理バンド計算がある。クラスターモデルは、全電ポテンシャルで接続し，完全に ”self-consistent” な数

エネルギー・原子間力のクラスターサイズ依存性の問題値解を得た。白金 (111) 面−極性溶液系の電子密度、溶や電極表面が電場を遮蔽する効果を正確に計算できない

液構造、固液界面をよぎる静電ポテンシャルを議論する。

のでここでは、後者のスラブ周期境界条件を用いたバン

我々の方法は Bérard らによる方法と基本的に同じでド計算の手法を用いた。計算は、ノルム保存擬ポテンシャ

あるが、金属表面を第一原理計算で扱っていることに違ル法を用い、電子の交換相関相互作用には、局所密度近いがある。白金 (111) 面を 11 層持つ繰り返しスラブを金

2 通常の第一原理計算で用いる原子に対するノルム保存擬ポテン

属側の計算に用いた。金属８層分に相当する厚みをもつ

シャル（またはウルトラソフト擬ポテンシャル）とは異なり、水分子に対して擬ポテンシャルを構成する13) 。静電相互作用や分極、コア反発の効果は含まれるが、電極表面で化学吸着する表面近傍の水分子に対してこのポテンシャルを使う近似はおそらく正しくないだろう。

非金属相（真空あるいは溶液）がそのスラブ間にある。
溶液側の積分方程式と結合して行ったポテンシャルの自2

己無撞着計算より、ここで用いたスラブとその非金属領

locで定義する。vion
は擬ポテンシャルの局所成分32) であ域の厚みは十分な厚みを持つ。溶液側の計算では RHNCる。溶液の電極の電子への相互作用は１次元であるとす理論が使われているが、双極子液体が半無限のスラブとる。表面平行方向で平均化された１次元局所有効ポテン

相互作用するとして計算した。液体分子は直径 ds = 0.28

シャルおよび電子密度は、逆格子空間で自己無撞着に解nm もつ剛体球で, 剛体球の中心に双極子 µ をもつとしかれた Kohn-Sham 方程式の３次元ポテンシャルよりた。（以後液体分子を溶媒分子と呼ぶ。）溶媒分子に対

loc
veff
(z)

して金属は剛体壁として働くと仮定した。その位置は表

=

X

loc
veff
(0, 0, Gz ) exp(iGz z)(5)ρ(0, 0, Gz ) exp(iGz z)

(6)

Gz

面の白金原子から溶媒分子の中心まで 0.20nm とした。

ρ(z)これは Pt(111)/water(0.1M KOH solution) in-situ Ｘ線

=

X
Gz

散乱の実験に基づくものである22) 。その実験では、いかここで、Gz は逆格子ベクトルの表面垂直方向の成分で、なる電位でも表面再構成しないこと水分子の層は密度上式において平行方向の成分は Gx = Gy = 0 である。

0.8/Pt atom で表面白金原子から約 0.2nm のところにあ

電極上の溶媒（古典的双極子）と電極表面は静電ポテン

ると報告している。

シャルで相互作用するので、静電ポテンシャルを金属は、表面近傍の液体を分極させるのに十分な強い電場を与える。液体は逆に金属をポテンシャル（溶媒分

loc
velst (z) = veff(z) − vxc [ρ(z)](7)

子の配置と配向性に依存する）を通して分極させる。静電場を通した２相間の相互作用をのみを考え、相間の電locで求める。Fig.2 に 溶媒分子がない時の veff(z), vxc (z),

荷・物質の移動はないと仮定した。

velst (z) と ρ(z) を示す。例えば Al(111) 面のような単純金属の場合33) 、１次元電子密度 ρ(z) は金属内でほぼ一

まず、白金表面（真空）の第一原理計算を行った。解くべき Kohn-Sham 方程式は
h h̄2
X ps
v̂ion (|r − R − rj |) + vH (r)
−
∇2 +
2m
R,rj
i
+vxc [ρ(r)] + vdip (z) ψnk (r) = εnk ψnk (r)定となり jellium モデルが良い近似となるが、Pt(111) 面の場合は原子周りに局在したｄ軌道のため、jellium モ
デル良い近似にはなり得ない。表面から長距離離れた場
所では、有効ポテンシャルは漸近的に鏡像ポテンシャル(1)−1/z 近づくはずであるが、LDA の近似によりより早ここで、nk は第一ブリルアンゾーン内の k 点の n 番目

くゼロに減衰する。電極近傍にある溶媒の感じる静電ポのバンド指標である。vH (r) はZ
ρ(r0 )vH (r) = dr0|r − r0 |テンシャルは同様に早く減衰するが、z = ds においても kB T に比べ十分大きい。電極上の 0.2-0.3 nm の所で、

(2)

電場は 1010 V/m にもなることをここでは強調しておきたい。

で、電子密度は nk での重み因子 wnk を乗じた波動関数の絶対値の２乗より

ρ(r) =

X2

wnk |ψnk (r)|

(3)

nk
ps
である。vxc は交換相関ポテンシャルである。 v̂ion
は擬

ポテンシャルで、R は実格子ベクトル、rj は原子の位置ベクトルである。また、vdip は、双極子溶媒分子から電子に及ぼすポテンシャルである。
我々の遷移金属単結晶面及びその吸着構造の第一原
理計算は実験結果とよく一致することを既に報告している28–31) 。真空中での白金 (111) 清浄表面に関する計算の結果、表面エネルギー (11 eV/nm2 、実験値 11 eV/nm2 ),仕事関数 (6.0 eV、実験値 6.1 eV)、表面緩和 (表面層と第２層の面間はバルクの+1.4%、実験値 1.1%) は報告さ

Figure 2: xy-averaged effective local (solid line) , exchange-

れている実験結果と良く一致した。

correlation (dotted line), electrostatic potential (thick solid line),and pseudo-valance-charge (dashed line) of Pt(111) 11 layer slab
without solvent molecules. The atomic positions of platinum ion
are shown in the closed circle.

loc
今、有効局所ポテンシャル veff
を
loc
loc
veff
= vion+ vH + vxc(4)3電極と溶媒の相互作用は

½
ums =

∞,
(z)
cos θ,− µe dvelstdz|z| < 0|z| > 0(8)

ここで θ は溶媒の双極子ベクトルと表面法線のなす角である。溶媒間の相互作用は双極子-双極子相互作用であり、短距離の所での反発を直径 ds を持つ剛体球で近似した。積分方程式論は成書34) に詳しいので省略するが、溶媒の多体の相互作用は積分方程式論の一つである RHNC法を用いた15) 。溶媒のバルク密度は、ρ∗s = ρs d3s = 0.7√に, 双極子モーメントは µ∗ = µ/(d3s kB T )1/2 = 3 に,
比誘電率は ² ∼ 50 とした3 。
第一原理計算で得られた静電ポテンシャル velst (z) を用いて表面近傍での液体の密度、配向を求め、溶媒の配Figure 4: Self-consistent reduced density profile g000 and projection g 011 near the metal wall are shown as the function of the
distance from the surface. The distance from the suface is sclaedby the dimeter of the hard sphere in the figure. g 000 near theinert wall is also shown. g 011 near the inert wall is zero at all
separations. The values of g 000 (Inert Wall), g 000 (Metal Wall),and g 011 (Metal Wall) at z/ds = 0 are 4.28, 29.94, and −71.49,
respectively.

置、配向から双極子溶液（双極子ポテンシャル）と電子との相互作用 vdip (z) の大きさを決定した9) 。この双極子ポテンシャルを取り入れることにより金属側の電子が感じる電場は更新され、さらにその場を使って溶液側の計算を行うことをくり返して、自己無撞着場を決定した15) 。表面平行方向に平均化した電子の電荷分布の自己溶媒 (s) 間の対分布関数 gms を以下のように展開する。
X
gms =(−1)n g 0nn (z)Pn (cos θ)(9)n

こ こ で 、θ は 双 極 子 ベ ク ト ル と 表 面 垂 線 の な す 角 、Pn (cos θ) は Legendre 多項式である。g 000 は換算動径分布関数、g 011 は双極子の配向に関した分布関数に対応
し15) 、溶媒配向の平均と以下の関係がある。

hcos θi = −g 011 (z)3g 000 (z)(10)

Fig.4 に実線で示した固液界面での溶媒分子の動径分布
g 000 は金属界面がない（電場のない）単なる剛体壁（破Figure 3: Electron density profile for Pt(111)/liquid interface. 線）の密度分布と異なり、より界面近傍に溶媒分子が近The self-consistent density profile ρscf is shown in the dottedline and the density change ∆ρ(≡ ρscf − ρvac , where ρvac is thedensity in the vacuum case) caused by the presence of the liquid
is shown in the solid line. The atomic positions of platinum areshown in solid circles.

づいたことを示している。（１）吸着層に対応する顕著な３つのピークの数及びその位置、（２）第一ピーク層の面密度（バルクの密度に比較して約２倍) は、Pt(111)｜水溶液界面に対して得られている実験結果22) あるいはより詳しい解析がなされた Ag(111)｜水溶液界面系 (電極無撞着な計算結果を Fig.3 に示す。電荷密度は金属イオ

が負に帯電した時バルクの 1.4 倍、正に帯電した時 2.3ンを電子原子レベルで扱ったために金属内で強く振動し倍)20, 21) を良く再現した。4ている。液体の存在により金属表面から液側に向かって

4 Toney ら20, 21) の Ag(111) 電極の場合を除いて電極上の溶媒の構造の実験結果はほとんど得られてない。Tidswell らの Pt(111) 電極の場合溶媒構造は十分に決定されていない22) 。ただし、Toney ら20, 21)の結果において、特に電極が正に帯電した時、第一ピークに相当す
る表面水密度は過大評価しているのではないかというのが定説になりつつあるようだ。実験の分解能が良くないために起こった artifact
で、実際にはさらに二本のピークがあるとの計算結果が Halley8) 及び Izvekov ら6) によって報告されている。UHV-STM の測定においても ice-like bilayer から融解時表面の密度は若干増えるだけであると報告されている24) 。

電子がさらに浸みだしている様子が Fig.3 の ∆ρ より明確である。溶液側の構造を以下に議論する。電極 (m)3 ここでの理論解析において、溶媒の密度は ρ = 3.2 × 1028 分子/
m3 , 双極子モーメントは 1.65 D であり、溶媒として水を想定している。水の 300 K での密度は ρ = 3.3 × 10 28 分子 / m3 , 水分子の双極子モーメントは 1.86 D、誘電率は 78 である。また、酸素の半径は0.14 nm、水分子中での OH 間の距離は 0.0957 nm である。

4

Fig.5 に示すように界面では、溶媒分子がもつ双極子（孤立水分子の場合酸素から両水素の重心に向かって 1.84D = 6.13 ×10−30 C ｍの双極子がある。）はその向きを界面から外側に向けて配向した。この配向秩序に
より溶媒分子は強い引力を金属から感じ、電場のない場
合の剛体壁に比べより密な溶媒層を形成した。界面からさらに離れるに従い急速にバルクでのランダム配向、密度に近づく。Ataka らは赤外分光により Au(111) 面上の水分子の配向を明らかにした23) 。ゼロ電荷電位附近で、電極表面第一層の水分子は水素結合ネットワークを保持しようとして水分子はその双極子を表面平行方向にむけていると報告している。我々のモデルでは、溶媒分子間の水素結合的な相互作用は考慮されてないために直接のFigure 6: The electrostatic potential of the Pt(111)/liquid sur-

比較はできないが、最近の Izvekov ら6) による Ag(111)face. Note that the energy axes are scaled reverse in order tocompare the conventional electric potential. The total electrostatic potential(solid line) is the sum of the contributions fromthe liquid part[dashed line(1)] and the metal part[long-dashedline(2)]. The electrostatic potential for the Pt(111)/vacuum system is shown in the dotted line(3).

｜水溶液系の ab-initio 分子動力学計算では、酸素を電極側にして双極子は表面垂直に対して６０度の角度をもつ
水分子が一番多く分布していると報告している。

0.94eV 減少している。in-situ 電気化学 STM によるトンネル電流の tip 位置依存性から、見かけのポテンシャ
ル障壁は、真空の時に比べて固液界面では (数 meV から２ e Ｖ）減少することが報告されている。もし、電極の電子ポテンシャルが溶媒の有無にあまり依存しないな
ら、障壁の高さが減少することはポテンシャルドロップを引き起こすことに相当する。分子動力学計算では、このポテンシャルドロップは 0.2-0.3 Ｖであると報告されている35) 。我々の求めた値が過大評価しているのは、溶
媒に水素結合の影響が入っておらず、より強く電極の電子のしみ出しが作るポテンシャルに影響を受けたためであると考えられる。

Figure 5: Self-consistent orientational probability densitiesp(z, θ) for dipole vectors at z/ds = 0, 1, and 2 (the first, second, and third layers, respectively).

金属の第一原理計算と RHNC 法との結合の自己無撞
着な計算はここで述べた場合にだけ限られたものではな
固体真空表面、固液界面を横切る静電ポテンシャルをいということを指摘したい。(1) より洗練された水のモ

Fig.6 に示す。ポテンシャルは金属系を電子原子レベル

デル、(2) 電解質イオンの導入、(3) より現実的な金属−で扱っているため金属内で強く振動している。全ポテン

水、(4) 金属−イオン相互作用の考慮へも容易に拡張す

シャル (1)+(2) は金属 (2) と液体 (1) の部分の和である。 る事が可能である。さらには、(5) 表面の原子レベルで液体が存在すると金属のポテンシャルドロップに対するの凹凸の影響（例えば Step）、(6) 電極に電位を与えた

寄与 (2)-(3) は金属がない場合と比べて 0.37eV 大きく時の影響、(7) 金属表面構造の変化、(8) 電荷容量、(9)なった。これは双極子ポテンシャルが電子をさらに界面

水分子の吸着の量子力学的効果 (10) その電位依存性に

外側に引き出すことと一致する。溶液の存在による界面よる水の解離, 水素発生、酸化等の化学反応等多くのこをよぎる総計のポテンシャルドロップ (1)+(2)-(3) は、 とを今後解決してする必要がある。

54帯電した Au(111) 面上の硫黄の吸着これまでは、帯電した電極表面については取り扱わなかった。以下は、溶媒分子は存在しないが帯電した金属表面およびその吸着エネルギー変化について報告する36) 。

Andreasen らは37–41) 、Au(111) 上の硫黄原子の吸着組成・構造の電位依存性を報告している。被覆度 2/3 の S8構造から電位を pzc から約 0.4 - 0.5 V 負な電位にする√√と、被覆率 1/3 の 3 × 3R30˚に可逆的に相変化する。
この相変化が (1) 電極の帯電によるものか (2) 溶媒分子を含めた電気化学的な影響によるかを明らかにするため帯電の影響による変化をここでは報告する。
帯電の方法は既に報告されている方法を用いた42–44) 。
即ち、仮想的な電荷シート (ガウス分布させると逆格子

Figure 8: Screened potential profile ∆veff (z) for Au(111)
surface.
The difference ∆veff (z) is defined by ∆veff (z) =veff (z, σM ) − veff (z, σM = 0). The effective electron potentials
for σM = 0.0(broken line) and σM = +0.264 C / m2 (dottedline) are also shown.

空間で扱いやすい。）を金属スラブの中間に置く。Fig.7に周期セルと仮想電荷シートの配置図を示す。仮想的なは表面平行方向は必ずしも均一ではない状態で分布している。この分布の違いが、吸着子との相互作用に影響を与える可能性がある。金属表面を負に強く帯電させると、金属表面から仮想的な電荷シートへの電子の移動が(Fig.9 で σM = −0.132 C/m2 より負電荷密度の場合) あ
Figure 7: Schematic diagram of the electrified surface used in り、注意を要する。
the slab model calculation.

電荷が作る電場は金属表面の電子によって完全に遮蔽される。Au(111) 面上での遮蔽されたポテンシャル断面図
は Fig.8 のように、また遮蔽された電子の再配置は Fig.9のようになる。Fig.8 では、仮想的な電荷シートの電荷がゼロの時と負の電荷 (σ= - 0.264 C/m2 = −σM ）を置いた時の有効ポテンシャルの差をとった図である。ここで、σ, σM は電荷シートおよび金属の表面電荷密度である。電子のポテンシャルは仮想的な電荷シートから金属表面に向かってほぼ直線的に減少し、5 層の金属スラブ
内では、ポテンシャルの変化はほぼゼロである事がわかる。直線の部分から求めた電場は電磁気学から得られた

Figure 9: Screening charge distribution for Au(111) surface.

結果

σ
Ez =
2²0

2D map for σM = −0.044 C / m2 (the upper figure) and surfaceaveraged screening charge profile.

(11)に一致する。ここで、²0 は真空の誘電率である。電磁気学の教科書では、遮蔽電荷は表面にデルタ関数的に分√√電荷シートがないときの ( 3 × 3)R30˚構造の吸着エ

布してると近似しているが、微視的に見れば有効ポテン

ネルギー、清浄表面からの仕事関数変化、表面緩和、硫シャルの小さい表面 Au 原子の外側（矢印で screening

黄原子への電荷移行量を Table1 に示す。FCC サイトがcharge と示したところ）に多く存在し、また遮蔽電荷

最も安定であることは実験結果37–41) と一致した。S8 構

6√√
Table 1: Calculated physical properties for Au(111)-( 3 × 3)R30˚-S structure
Adsorption
Site
FCC
HCP
ATOP

Adsorption
Energy(eV / S atom)-4.950-4.797-3.440Workfunction
Change
(eV)
+0.22
+0.24
+0.85

Interlayer Distance (Å)
S-Au1 Au1-Au2 Au2-Au3

δQS
(e)

1.60
1.65
2.23

-0.15

2.38
2.39
2.41

2.38
2.36
2.39

-0.02

Table 2: Calculated physical properties for Au(111)-(1 × 1)-S structure
Adsorption
Site
FCC
HCP
ATOP

Adsorption
Energy(eV / S atom)-3.780-3.788-3.968Workfunction
Change
(eV)
+0.11
+0.14
+0.19

Interlayer Distance (Å)
S-Au1 Au1-Au2 Au2-Au3

δQS
(e)

2.12
2.12
2.40

-0.13
-0.12
-0.04

2.36
2.35
2.35

2.37
2.37
2.37

造は下地の金属と commensurate の関係にないので, 仮想的な (1 × 1) 構造を仮定して計算を行った。帯電していない時の (1 × 1)-S 構造の結果を Table2 に示す。
仮想的な電荷シートは人工的に挿入したので、系の全エネルギーは物理的な意味を持たない。ただし、２つの系で同じ電荷シートをもつ場合、その全エネルギーの差は、電荷シートの自己エネルギー等の付加的な部分
がキャンセルするため、帯電した２つの系の表面エネル
ギーの差として定義できる。今、Au(111) 面上の硫黄の吸着エネルギの帯電による変化を以下の式で定義する。

∆Ead (σM )=£E(σM ; S adsorption system)−E(σM ; Au(111) system)¤−2E(S atom) /2√√

Figure 10: Electrocapillary curves of Au(111)-( 3× 3)R30˚-S(12)structure and Au(111)-(1 × 1)-S structure.

ここで、σM は金属表面の表面電荷密度、上式右辺第一項は電極表面が σM に帯電した時の金-硫黄系スラブの全

実際は水溶液に脱離していく硫黄を考慮しなくてはなら

エネルギー、、第二項は Au(111) 面が帯電した時のスラないので、帯電状態だけで実験結果を説明することはでブの全エネルギーである。E(S atom) は硫黄（擬）原子

きないが、傾向は矛盾しないことが明らかとなった。

の全エネルギーである。表面が２つあるので、2 の係数√ √がついている。Au(111)-( 3 × 3)R30˚-S 構造と (1 × 1)5構造について ∆Ead (σM ) を計算し、帯電していない時のまとめ∆Ead (0) をゼロにした ∆∆Ead を Fig.10 に示す。Fig.10以上の白金上の溶液構造と金表面の帯電による吸着

は、電極を帯電させると、電荷の符号によって非対称でエネルギーの２例を電極の第一原理計算として紹介したはあるが吸着エネルギーは安定化することを示していが、現実の電気化学系との間には大きな隔たりがあるこる。実験結果は、電極の電位を負にすることにより高被√√覆率の吸着相から低被覆率の ( 3 × 3)R30˚相に変化√√している。Fig.10 では低被覆率の ( 3 × 3)R30˚相がとは確かである。よりミクロスコピックな観点にたった(1 × 1) 相にくらべて、負の帯電に対してエネルギーがよる固体表面の電気化学を新たな視点から見直すことがで分光実験を含む実験結果の積み重ねと理論の改良の突き合わせによって、既に理解しつくされたと考えられていり安定化している。このことは実験結果と矛盾しない。 きるかもしれない。

7

6

謝辞

25) H. Siegenthaler: “STM in Electrochemistry”, vol. II of Scanning Tunnelling Microscop (Springer, Berlin, 1992) ch.2,p.7.

この研究を行うにあたって、共同研究者の京都大学エ26) A. Peremans and A. Tadjeddine: J. Chem. Phys. 103, 7197
(1995).

ネルギー理工学研究所の木下正弘先生、京都大学大学院
工学研究科垣内隆先生、香港科学技術大学物理学科の C.

27) N. D. Lang and W. Kohn: Phys. Rev. B 1, 4555 (1970).

T. Chan 先生、Ames 研究所-アイオワ州立大物理学科の28) M. Yamamoto, C. T. Chan, and K. M. Ho: Phys. Rev. B
50, 7932 (1994).

K. M. Ho 先生、Max-Plank 金属研究所の C. Elsässer 29) M. Yamamoto, C. T. Chan, K. M. Ho, and S. Naito: Phys.
Rev. B 54, 14111 (1996).

先生に感謝いたします。

30) M. Yamamoto, C. T. Chan, K. M. Ho, M. Kurahashi, and
S. Naito: Phys. Rev. B 53, 13772 (1996).

文献31) M. Yamamoto, M. Kurahashi, C. T. Chan, K. M. Ho, and
S. Naito: Surf. Sci. 387, 300 (1997).

1) P. J. Feibelman, B. Hammer, J. K. Nørskov, F. Wagner,
M. Scheffler, R. Stumpf, R. Watwe, and J. Dumesic: J.
Phys. Chem. B 105, 4018 (2001).

32) M. L. Cohen and V. Heine: “The fitting of pseudopotentials to experimental data and their subsequent application”, vol. 24 of Solid State Physics (Academic, New York,1970) p.37.

2) I. Grinberg, Y. Yourdshahyan, and A. M. Rappe: J. Chem.
Phys. 117, 2264 (2002).

33) J. R. Chelikowsky, M. Schlüter, S. G. Louie, and M. L.
Cohen: Solid State Commun. 17, 1103 (1975).

3) C. P. Ursenbach, A. Calhoun, and G. A. Voth: J. Chem.
Phys. 106, 2811 (1997).

34) J.-P. Hansen and M. I. R.: “Theory of Simple Liquids 2nded.” (Academic, London, 1990).

4) A. Klesing, D. Labrenz, and R. A. van Santen: J. Chem.
Soc.-Faraday Trans. 94, 3229 (1998).

35) H. Spohr: “Computer Simulation of Electrochemical Interfaces”, vol. 6 of Advances in Electrochemical Science andEngineering (Wiley-VCH, Weinheim, 1999) p.1.

5) S. Izvekov, A. Mazzolo, K. VanOpdorp, and G. A. Voth: J.
Chem. Phys. 114, 3248 (2001).

36) M. Yamamoto and C. T. Chan: To be submitted.

6) S. Izvekov and G. A. Voth: J. Chem. Phys. 115, 7196
(2001).

37) G. Andreasen, C. Vericat, M. E. Vela, and R. C. Salvarezza:
J. Chem. Phys. 111, 9457 (1999).

7) 山本雅博: To be published in 電気化学会誌 Q&A (2003).

38) C. Vericat, G. Andreasen, M. E. Vela, and R. C. Salvarezza:
J. Phys. Chem. B 104, 302 (2000).

8) J. W. Halley, S. Walbran, and D. L. Price: “Chemical
Physics of the electrode-electrolyte interface”, vol. 116 ofAdv. Chem. Phys. (John Wiley & Sons, New York, 2001).

39) H. Martin, C. Vericat, G. Andreasen, A. H. Creus, M. E.
Vela, and R. C. Salvarezza: Langmuir 17, 2334 (2001).

9) D. R. Bérard, M. Kinoshita, X. Ye, and G. N. Patey: J.
Chem. Phys. 101, 6271 (1994).

40) C. Vericat, M. E. Vela, G. Andreasen, R. C. Salvarezza,
L. Vazquez, and J. A. Martin-Gago: Langmuir 17, 4919
(2001).

10) D. R. Bérard, M. Kinoshita, X. Ye, and G. N. Patey: J.
Chem. Phys. 102, 1024 (1995).
11) M. Kinoshita and D. R. Bérard: J. Comput. Phys. 124, 230(1996).

41) C. Vericat, M. E. Vela, G. A. Andreasen, R. C. Salvarezza,
F. Borgatti, R. Felici, T. L. Lee, F. Renner, J. Zegenhagen,and J. A. Martin-Gago: Phys. Rev. Lett. 90 (2003).

12) D. R. Bérard, M. Kinoshita, N. M. Cann, and G. N. Patey:
J. Chem. Phys. 107, 4719 (1997).

42) C. L. Fu and K. M. Ho: Phys. Rev. Lett. 63, 1617 (1989).
43) K. P. Bohnen and D. M. Kolb: Surf. Sci. 407, L629 (1998).

13) J. Schnitker and P. J. Rossky: J. Chem. Phys. 86, 3462.

44) J. G. Che and C. T. Chan: Phys. Rev. B 67 (2003).

14) M. Yamamoto and M. Kinoshita: Chem. Phys. Lett. 274,513 (1997).
15) M. Yamamoto, M. Kinoshita, and T. Kakiuchi:
trochim. Acta 46, 165 (2000).

Elec-16) A. Kovalenko and F. Hirata: J. Chem. Phys. 110, 10095(1999).
17) A. Kovalenko and F. Hirata: J. Mol. Liq. 90, 215 (2001).
18) L. Hedin and B. I. Lundqvist: J. Phys. C 4, 2064 (1971).
19) C. Elsässer, N. Takeuchi, K. M. Ho, C. T. Chan, P. Braun,and M. Fähnle: J. Phys.-Condes. Matter 2, 4371 (1990).
20) M. F. Toney, J. N. Howard, J. Richer, G. L. Borges, J. G.
Gordon, O. R. Melroy, D. G. Wiesler, D. Yee, and L. B.
Sorensen: Nature 368, 444 (1994).
21) M. F. Toney, J. N. Howard, J. Richer, G. L. Borges, J. G.
Gordon, O. R. Melroy, D. G. Wiesler, D. Yee, and L. B.
Sorensen: Surf. Sci. 335, 326 (1995).
22) I. M. Tidswell, N. M. Markovic, and P. N. Ross: J. Electroanal. Chem. 376, 119 (1994).
23) K. Ataka, T. Yotsuyanagi, and M. Osawa: J. Phys. Chem.
100, 10664 (1996).
24) M. Morgenstern, J. Müller, T. Michely, and G. Comsa: Z.
Phys. Chemie-Int. J. Res. Phys. Chem. Chem. Phys. 198,
43 (1997).

8