表面と吸着:Langmuir等温式を占有率から導く

表面を「数の限られた席」と見なし、空席への吸着と占有席からの脱離のつり合いからLangmuir等温式を導く。低圧の比例則・高圧の飽和・半占有圧力を計算し、モデルの仮定とデータ解釈の限界を区別する。

この単元の目標
  • 吸着と吸収を区別し、占有率と空席率をサイト数から定義できる。
  • 吸着・脱離速度を占有率で表し、平衡条件からLangmuir等温式を導ける。
  • 吸着量・親和性パラメータ・半占有圧力を単位とともに計算できる。
  • 等価なサイト・単分子層・相互作用なしという仮定と、その破綻例を説明できる。
📑 目次
  1. 1.なぜ表面の量は、圧力に比例し続けないのか
  2. 2.「何個ついたか」を「何割埋まったか」へ変える
  3. 3.式の前に、席についての約束を置く
  4. 4.空席に入る速さと、占有席から出る速さ
  5. 5.比例・半占有・飽和を同じ式で読む
  6. 6.容量と親和性は違う情報を持つ
  7. 7.どこで単純な席モデルが破れるか

1.なぜ表面の量は、圧力に比例し続けないのか

空の駐輪場なら自転車を増やすほど入庫できますが、満車に近づくと空き場所が足りなくなります。表面への吸着も、分子がとどまれる場所に限りがあるなら同じ制約を受けます。一方、停めた自転車が出ていくように、表面の分子が離れる「脱離」も起こり得ます。入る数と出る数を別々に考えることが、飽和する曲線を理解する出発点です。この比喩は席数の制限だけを表し、実際の分子に意思があるという意味ではありません。

吸着は物質が界面付近に集まる現象で、物質の内部へ取り込まれる吸収とは区別します。内部に細孔がある固体では、外側だけでなく、分子が到達できる孔の壁も表面です。本記事は一定温度で、一種類の気体分子Aが固体表面に解離せず吸着する理想モデルを扱います。実験操作ではなく、架空の数値を用いた理論とデータ読解のレッスンです。

2.「何個ついたか」を「何割埋まったか」へ変える

駐輪場の大きさが違っても混雑を比べられるように、吸着できる場所(サイト)の総数を基準にします。総サイト数を Nₛ、占有サイト数を Nₐ とします。一分子が一サイトを使い、一サイトには一分子までなら、占有率 θ と空席率は次のように決まります。

θ=NaNs,1−θ=Ns−NaNs,0≤θ≤1\theta=\frac{N_a}{N_s},\qquad 1-\theta=\frac{N_s-N_a}{N_s},\qquad 0\leq\theta\leq1

例えば100個のサイトのうち30個が占有されていれば θ = 0.30、空席率は0.70です。θ は個数の比なので無次元です。ここでの「被覆率」は、この一分子一サイトのモデルでの占有率を意味します。異なる大きさの分子が複数のサイトを覆う状況では、面積の被覆率と単純なサイト占有率を同一視できません。

3.式の前に、席についての約束を置く

仮定式に効く理由
サイトは等価で、総数が変わらないどの席にも同じ吸着・脱離の係数を使える。
一分子一サイト、単分子層に限る占有済みの席の上には重ねて吸着せず、空席率は1 − θになる。
吸着分子どうしの横方向の相互作用を無視する周囲の混み具合で結合の強さや係数が変わらない。
一成分が解離せず可逆的に吸着する複数席を同時に必要とする解離や、他成分による席の競争を含めない。
温度一定で平衡を考え、気相は理想気体と近似する同じ温度の係数を使い、気体からの到来の効果を分圧に比例させる。
「表面に何かがつく」だけでは、この仮定一式を満たしません。またLangmuir式は、物理吸着か化学吸着かを名前だけで判定する式ではありません。解離、複数層、サイトの不均一性などを無視できる範囲かどうかが重要です。

4.空席に入る速さと、占有席から出る速さ

同じ温度で分圧を上げると、表面に到来するA分子は増えます。しかし到来しても空席がなければこのモデルでは入れません。したがって吸着による占有率の増加速度は「Aの分圧 p」と「空席率」の積に比例します。脱離は、すでにいる分子だけに起こるので占有率に比例します。p は平衡時のAの分圧であり、一般には全圧や投入前の圧力ではありません。

rad=kadp(1−θ),rde=kdeθ,dθdt=rad−rder_{\mathrm{ad}}=k_{\mathrm{ad}}p(1-\theta),\qquad r_{\mathrm{de}}=k_{\mathrm{de}}\theta,\qquad \frac{\mathrm{d}\theta}{\mathrm{d}t}=r_{\mathrm{ad}}-r_{\mathrm{de}}

r は物質量そのものの速度でなく、単位時間あたりの占有率の変化です。時間をs、pをkPaで表すなら、r の単位は s⁻¹、kₐd は kPa⁻¹ s⁻¹、kde は s⁻¹です。係数は正とし、同じ温度では一定とします。表面までの輸送の遅れや、表面での別反応はこの式には含めません。

🔍 平衡の占有率を一段ずつ解く

  1. 平衡では増加と減少がつり合う
    kadp(1−θ)=kdeθk_{\mathrm{ad}}p(1-\theta)=k_{\mathrm{de}}\theta

    左は空席への流入、右は占有席からの流出です。両方がゼロになる必要はありません。

  2. 脱離係数で割り、比をまとめる
    b=kadkde,bp(1−θ)=θb=\frac{k_{\mathrm{ad}}}{k_{\mathrm{de}}},\qquad bp(1-\theta)=\theta

    b はこの温度での吸着の親和性を表す圧力基準のパラメータです。

  3. θを含む項を一方へ集める
    bp=θ+bpθ=θ(1+bp)bp=\theta+bp\theta=\theta(1+bp)

    分母の1 + bpは、有限の空席を考えた結果として現れます。

  4. Langmuir(ラングミュア)等温式
    θ=bp1+bp\boxed{\theta=\frac{bp}{1+bp}}

    温度を一定にして、平衡分圧に対する占有率を表すので「等温式」と呼びます。時間発展そのものの式ではありません。

平衡は「吸着も脱離も止まる」ことではなく、平均的な占有率が変わらない動的なつり合いです。また、平衡曲線だけから平衡に達する時間は決まりません。kₐd と kde を両方同じ倍率で大きくしても、その比 b と平衡曲線は変わらないからです。

式の1と足し合わせる bp は無次元でなければなりません。したがって b の単位は圧力の逆数です。例えば b = 0.20 kPa⁻¹ は 2.0×10⁻⁴ Pa⁻¹と同じ値を表します。b 自体を無次元の熱力学的平衡定数と混同しないようにしましょう。標準圧力で無次元化する議論とは区別し、本記事では一貫して圧力基準の b を用います。

5.比例・半占有・飽和を同じ式で読む

まず bp を計算して、空席が多いか少ないかを見通します。bp が小さいときは分母がほぼ1で、占有率は圧力にほぼ比例します。bp = 1 なら占有と空席が半々です。bp が大きいと、分子が来ても入れる席が少なくなるため、占有率は1へ近づきます。

bp≪1: θ≈bp,p1/2=1b,bp≫1: 1−θ=11+bp≈1bpbp\ll1:\ \theta\approx bp,\qquad p_{1/2}=\frac{1}{b},\qquad bp\gg1:\ 1-\theta=\frac{1}{1+bp}\approx\frac{1}{bp}

「高圧で飽和」はこのモデル内での極限です。有限の正の bp では θ < 1 のままで、ちょうど1になるわけではありません。また低圧という言葉も圧力の数字だけでは決まりません。b が大きい表面なら、小さな分圧でも bp は大きくなります。現実の高圧域では気体の非理想性や多層吸着などにより、極限へ行く前に仮定が破れる場合があります。

席の埋まり具合だけでなく、どれだけの量を持てるかも必要です。固体の単位質量あたりの吸着量を q、全サイトが一層だけ埋まったときの容量を qₘ とすれば、一分子一サイトの仮定から q = qₘθ です。q と qₘ の単位を同じにすれば、例えばどちらも mmol/g として扱えます。

q=qmbp1+bp,p=θb(1−θ)(0≤θ<1)q=q_m\frac{bp}{1+bp},\qquad p=\frac{\theta}{b(1-\theta)}\quad(0\leq\theta<1)
例題1:分圧を2倍にしても、吸着量は2倍にならない

一定温度の架空モデルで b = 0.20 kPa⁻¹、qₘ = 3.0 mmol/g とする。p = 5.0 kPa と10 kPaの占有率と吸着量を求めよ。

解答・解説

5.0 kPaでは bp = 1.0、θ = 1/(1 + 1) = 0.50、q = 3.0×0.50 = 1.5 mmol/g。10 kPaでは bp = 2.0、θ = 2/3 ≈ 0.667、q = 2.0 mmol/g。吸着量の倍率は2.0/1.5 = 4/3 ≈ 1.33倍で、2倍ではありません。最初から半分の席が埋まっており、低圧近似 bp ≪ 1 は使えないためです。

例題2:半占有と90%占有に必要な分圧を比べる

同じ温度の架空表面について b = 0.50 kPa⁻¹とする。θ = 0.50 と θ = 0.90 を与える分圧を求め、その倍率を説明せよ。

解答・解説

半占有では p = 1/b = 2.0 kPa。90%では p = 0.90/[0.50×(1 − 0.90)] = 18 kPaで、半占有圧力の9倍です。90%占有なら占有席と空席の比は9:1なので bp = 9が必要になります。検算:0.50×18/(1 + 0.50×18) = 9/10 = 0.90。満席に近づくほど、残りの空席を埋めるための分圧の増加が大きくなります。

6.容量と親和性は違う情報を持つ

大きな駐輪場であることと、同じ混雑条件で埋まりやすいことは別です。qₘ は席の総量に対応し、b は半占有圧力の逆数に対応します。低圧だけでは q ≈ qₘbp となり、データの傾きからわかるのは積 qₘb です。傾きが大きいだけでは「席が多い」のか「埋まりやすい」のかを切り分けられません。

pq=1qmb+pqm(p>0, q>0)\frac{p}{q}=\frac{1}{q_m b}+\frac{p}{q_m}\quad(p>0,\ q>0)

理想モデルなら p/q を縦軸、p を横軸にすると直線で、傾きは1/qₘ、切片は1/(qₘb)です。p = 0では qも0で、p/qを直接計算できません。直線の切片は正の圧力側からの極限です。以下では丸め誤差のない架空データを使い、パラメータの意味を逆向きに確かめます。

例題3:二点から容量と親和性を逆算する

同じ温度の理想Langmuirモデルの二点 (p, q) が (2.0 kPa, 1.0 mmol/g)、(6.0 kPa, 1.8 mmol/g) であった。qₘ と b を求め、元の式で検算せよ。

解答・解説

p/q は順に2、10/3 kPa·g/mmol。傾きは[(10/3) − 2]/(6 − 2) = 1/3 g/mmolなので qₘ = 3.0 mmol/g。切片は2 − (1/3)×2 = 4/3 kPa·g/mmol。b = 1/[3×(4/3)] = 0.25 kPa⁻¹です。元の式では p = 2で q = 3×0.5/1.5 = 1.0、p = 6で q = 3×1.5/2.5 = 1.8 mmol/gとなり、両方再現します。

二点はいつでも一本の直線を決めるため、二点が合うだけではモデルの妥当性を証明できません。実際のデータでは、複数点の残差や誤差を検討する必要があります。p/qへの変換は誤差の重みを変え、特に小さなqの誤差を大きく見せることがあります。直線に近い見た目だけで機構を断定せず、元のqと予測値のずれも読むことが大切です。

7.どこで単純な席モデルが破れるか

  • 席の不均一性:結晶面・欠陥・細孔などで結合の強さが違うと、一つのbで全サイトを代表できない場合がある。
  • 分子間相互作用:隣に吸着分子があると付きやすくなる、または付きにくくなるなら、係数がθに依存しないという仮定が崩れる。
  • 多層吸着・細孔内の凝縮:占有済みの分子の上にも積み重なる場合、qₘは総吸着量の普遍的上限ではない。
  • 競争・解離:他成分が同じ席を使う場合や、一分子が解離して複数サイトを使う場合は、空席の数え方や速度式を変更する必要がある。
  • 非平衡・非理想性:まだ平衡に達していない値を等温式に入れても、その式は到達過程を説明しない。理想気体近似が不適切なら分圧に代わる熱力学的な量の検討も必要になる。
計算前に①同じ温度か、②pは吸着成分の平衡分圧か、③bpの単位は消えるか、④求めたθは0〜1の範囲か、⑤仮定を破る事情がないか、と順に確認しましょう。値が式に合うことと、その式の前提が現実に正しいことは別の判断です。
練習:空席・単位・逆算を確かめる
  1. 同じ温度で b = 0.10 kPa⁻¹、qₘ = 5.0 mmol/g、p = 2.0 kPaとする。θ、空席率、qを求めよ。

    ヒント

    まず無次元のbpを作る。

    解答・解説

    bp = 0.20、θ = 0.20/1.20 = 1/6 ≈ 0.167。空席率は5/6 ≈ 0.833。q = 5.0/6 ≈ 0.833 mmol/gです。吸着量が全容量より小さく、占有率と空席率の和が1になることも検算です。

  2. 架空の同じ温度のモデルで、p = 4.0 kPaのとき θ = 0.80だった。bと半占有圧力を求めよ。

    解答・解説

    bp = θ/(1 − θ) = 0.80/0.20 = 4。したがって b = 4/4.0 = 1.0 kPa⁻¹、半占有圧力は1/b = 1.0 kPaです。p = 4.0 kPaならθ = 4/(1 + 4) = 0.80に戻ります。

  3. bp = 0.020のとき、低圧近似 θ ≈ bp は正確なLangmuir値を過大評価するか。相対誤差を「(近似値 − 正確値)/正確値」で求めよ。

    解答・解説

    正確値は0.020/1.020 ≈ 0.01961、近似値は0.020なので過大評価です。相対誤差は[0.020 − 0.020/1.020]/[0.020/1.020] = 0.020 = 2.0%。一般にこの定義の相対誤差はbpなので、近似の許容範囲を自分で決められます。

  4. kₐd = 0.040 kPa⁻¹ s⁻¹、kde = 0.10 s⁻¹、p = 5.0 kPaのモデルで、平衡占有率と、そのときの吸着・脱離速度を求めよ。

    解答・解説

    b = 0.040/0.10 = 0.40 kPa⁻¹、bp = 2.0だからθ = 2/3。rₐd = 0.040×5.0×(1/3) = 1/15 ≈ 0.0667 s⁻¹、rde = 0.10×(2/3) = 1/15 s⁻¹。両速度は非ゼロで等しく、差だけがゼロになる動的平衡を確認できます。

確認クイズ:式の理由と適用範囲

問1. 吸着速度に1 − θを掛ける理由として正しいものはどれか。

問2. 温度とqₘを同じにしてbだけを2倍にしたとき、必ず正しいものはどれか。

問3. pをkPaで表すとき、式θ = bp/(1 + bp)のbの単位はどれか。

問4. 純粋な一成分の平衡データがLangmuir曲線に近かった。そこからの結論として最も適切なのはどれか。

理論の参照先:Atkins・de Paula・Keeler『Atkins’ Physical Chemistry』第11版(2018年)の固体表面・吸着等温式の議論。本文では単一成分・非解離・等価サイトの速度つり合いから式を導出しています。数値はすべて理解のための架空の教材値であり、特定の物質の吸着性能や安全性を示す測定値ではありません。温度を変えたbの予測や、実験操作は扱いません。
📝 メモこの端末のみに保存
メモ欄はページ読み込み後に使えるようになります。