Takuya Kawanishi

6. 接触角現象と濡れ

6.1 ヤング方程式

接触角

  • 完全な濡れ(perfect wetting, complete wetting, total wetting) \(\theta = 0\)

  • 3 相接触線(濡れ線)の形成, 部分的な濡れ(partial wetting)

部分的な濡れ

  • \(\theta < 90{}^\circ\): hydrophilic

  • \(\theta > 90{}^\circ\): hydrophobic

  • 「濡れが好まれる」, 不自然な訳, Wetting is of the surface is favorable. むしろ, 「濡れは親和的」, といった訳が

6.1.1 接触角

  • 液・気・固 の 3 相が存在する場合

\[\begin{split}\begin{cases} \gamma_\mathrm L & \textsf{気液界面} \\ \gamma_\mathrm S & \textsf{固気界面} \\ \gamma_\mathrm {SL} & \textsf{固液界面} \end{cases}\end{split}\]
  • 気相 G の添字は省略されるのが普通. このほうが見やすい.

  • ヤング方程式

\[\gamma_\mathrm L \cos \theta = \gamma_\mathrm S - \gamma_\mathrm{SL} \tag{6.1}\]

6.1.2 導出

  • 力のつりあいから,

  • 図 6.2

\[\gamma_\mathrm{S} \, \mathrm d l - \gamma_\mathrm{SL} \, \mathrm d l - \gamma_\mathrm L \, \mathrm d l\, \mathrm \cos \theta = 0 \tag{6.10}\]
  • これからヤングの式が導出できる.

\[\gamma_\mathrm S - \gamma_\mathrm{SL} - \gamma_\mathrm L \, \mathrm \cos \theta = 0\]

液滴形状から

  • 平坦面状の液滴

  • 液滴が無限小だけ広がるときのギブズエネルギー変化

  • 重力の影響が無視できる場合, 液滴の形は, 球形帽子(球を平面で切り取った形)となる(なぜか考えよ).

  • 球の半径を \(R\) (変数)とする.

  • 球形帽子の表面積は

\[A_\mathrm L = \pi (a^2 + h^2) \tag{6.2}\]
  • 液滴の接触面(円形になる)の半径\(a\) から \(a + \mathrm d a\), 液滴の高さが \(h\) から \(h + \mathrm d h\) に変化したと仮定する.

  • \(\mathrm d a > 0\) ならば \(\mathrm d h < 0\) となる.)

  • 固液接触面積の変化は,

\[\boxed{ \mathrm d A_\mathrm{SL} = 2 \pi a \mathrm d a } \tag{dASL}\]
  • 液体の表面(気液接触面)の変化は,

\[\mathrm d A_\mathrm L = \frac{\partial A_\mathrm L}{\partial a} \, \mathrm d a + \frac{\partial A_\mathrm L}{\partial h} \, \mathrm d h = 2 \pi a \, \mathrm d a + 2 \pi h \, \mathrm d h \tag{6.3}\]
  • 液滴の体積一定の条件から(教科書 (6.4)–(6.6)), \(a\)\(h\) につぎの関係が成立

    \[\frac{\mathrm d h}{\mathrm d a} = - \frac{a}{R} \tag{6.7}\]
  • この関係を使って, (6.3) を変形して

\[\boxed{ \mathrm d A_\mathrm L = 2 \pi a \cos \theta \, \mathrm d a } \tag{6.8}\]
  • 表面のギブズエネルギーは, 界面張力 :math:`boldsymbol times` 面積で表されるので,

  • 液表面積(気液界面積)の変化にともなうギブズエネルギー変化は

\[\mathrm d G_\mathrm L = \gamma_\mathrm L \mathrm d A_\mathrm L \tag{dGA}\]
  • 固液接触面積の変化に伴うギブズエネルギー変化は

\[\mathrm d G_\mathrm{SL} = (\gamma_\mathrm{SL} - \gamma_\mathrm S) \, \mathrm d A_\mathrm {SL}\tag{dGSL}\]
  • (dGA), (dGSL) に (6.8), (dASL) を代入して, 平衡点(トータルのギブズエネルギー変化 = 0)を求める.

\[0 = \mathrm d G = 2 \pi a (\gamma_\mathrm{SL} - \gamma_\mathrm S) \, \mathrm d a + 2 \pi a \gamma_\mathrm L \cos \theta \, \mathrm d a\]
  • 整理して

\[\boxed{ \gamma_\mathrm{SL} - \gamma_\mathrm S + \gamma_\mathrm L \cos \theta = 0 }\]

導出の詳細 1

\[V = \frac{\pi}{6}( 3 a^2 h + h^3) \tag{6.4}\]
  • \(V\)\(a\)\(h\) との関数だから)

\begin{align*} \mathrm d V &= \frac{\partial V}{\partial a} \, \mathrm d a + \frac{\partial V}{\partial h} \, \mathrm d h \\ & = \frac\pi a h \, \mathrm d a + \frac{\pi}{2}(a^2 + h^2) \, \mathrm d h \, \tag{6.5} \end{align*}

\[- \pi a h \, \mathrm d a = \frac{\pi}{2}(a^2 + h^2), \quad \frac{\mathrm d h}{\mathrm d a} = - \frac{2ah}{a^2 + h^2} \tag{6.6}\]

導出の詳細 2

\[\frac{\mathrm d h}{\mathrm d a} = - \frac{2ah}{2Rh - h^2 + h^2} = - \frac{a}{R} \tag{6.7}\]

\begin{align*} \mathrm d A_\mathrm L & = 2 \pi a \, \mathrm d a - 2 \pi h \frac{a}{R} \, \mathrm d A = 2\pi a \left( 1 - \frac{h}{R}\right) \, \mathrm d a \\ & = 2 \pi a \frac{R- h}{R} \, \mathrm d a = 2 \pi a \, \cos \theta \, \mathrm d a \tag{6.8} \end{align*}

\begin{align*} \mathrm d G & = (\gamma_\mathrm{SL} - \gamma_\mathrm L) \, \mathrm d A_\mathrm {SL} \\ & = 2\pi a (\gamma_\mathrm{SL} - \gamma_\mathrm L ) \, \mathrm d a + 2 \pi a \gamma_\mathrm L \, \cos \theta \, \mathrm d a \tag{6.9} \end{align*}

6.1.3 線張力

  • 面積が広がるとき, 3 相接触線の長さも増加する.

  • 新しい接触線を形成する場合にエネルギーを必要とする.

  • 1 \(\mu\mathrm{m}\) 以下の液滴では, 線張力も考慮する必要がある.

  • \(\kappa_l\): 線張力

  • \(a\): 3 相接触線の曲率半径

\[\gamma_\mathrm L \, \cos \theta = \gamma_\mathrm S - \gamma_\mathrm{SL} - \frac{\kappa_l}{a} \tag{6.11}\]
\[\kappa_l = \frac{\Delta_\mathrm{vap} U}{3 N_\mathrm A a_\mathrm M}\]

具体例 6.2

  • シクロヘキサン液体の線張力を気化熱 \(\Delta_\mathrm{vap} U = 30.5 \, \mathrm{kJ \, mol^{-1}}\) から見積もれ.

  • 密度 \(\rho = 773 \, \mathrm {kg \, m^{-3}}\), 分子量 \(M = 84.2 \, \mathrm{g \, mol^{-1}}\), \(a_\mathrm M = 0.565 \, \mathrm{nm}\) とする.


解答例 6.2

\begin{align*} \kappa_l & = \frac{\mathrm{\Delta_\mathrm{vap} U}}{3 N_\mathrm A a_\mathrm M} \\ & = \frac{(30500)}{(3)(6.02 \times 10^{23})(0.565\times 10^{-9})} \\ & = 2.9890532 \times 10^{-11} = 2.99 \times 10^{-11} \ \mathrm {J/m} \end{align*}

6.1.4 完全なぬれと濡れ転移

  • 広がり係数(spreading coefficient)

\[S = \gamma_\mathrm S - \gamma_\mathrm{SL} - \gamma_\mathrm L \tag{6.12}\]
  • 平衡状態で完全濡れの場合は, \(S = 0\).

  • 非平衡状態では, \(S > 0\) になりうる.

  • \(S < 0\) では, 有限の接触角が形成される.


具体例 6.3

  • 実験室で湿度 100%でなくても, 水滴近くは湿度 100% とみなせる.

  • 空気中の水蒸気の拡散係数が十分遅いため, 表面付近は湿度 100 % に保たれる.


温度による表面張力・接触角の変化

\[\cos \theta = \frac{\gamma_\mathrm S - \gamma_\mathrm{SL}}{\gamma_\mathrm L}\]
  • 通常, 温度上昇により, \(\gamma_\mathrm L\)\(\gamma_\mathrm S - \gamma_\mathrm{SL}\) よりかなり速く減少し, 接触角 \(\theta\) は小さくなる.

  • 液体が沸点に達しない場合, \(S=0\) \((\theta=0)\) となる温度, 濡れ温度(wetting temperature)が存在する.

  • 濡れ転移

    • 濡れ温度 \(T_\mathrm w\) 以下では, 接触角が観測され,

    • 濡れ温度 \(T_\mathrm w\) 以上では, 液体膜が形成される.

6.1.5 接触角の理論的側面

  • 多くの場合, \(\gamma_\mathrm L\) の評価は可能, \(\gamma_\mathrm S\) の評価は難しいが可能.

  • \(\gamma_\mathrm{SL}\) は評価できない.

  • \(\gamma_\mathrm{SL}\)\(\gamma_\mathrm L\), \(\gamma_\mathrm S\) で表現する.

\[\gamma_{12} = \gamma_1 + \gamma_2 - w \tag{6.13}\]
\[w = \frac{A_{12}}{12 \pi D_0^2} \tag{6.14}\]
\[A_{12} \approx \sqrt{A_{11} A_{22}} \tag{6.15}\]
\[\gamma_1 = \frac{A_{11}}{24 \pi D_0^2}, \quad \gamma_2 = \frac{A_{22}}{24 \pi D_0^2} \tag{6.16}\]
\[\gamma_{12} = \gamma_1 + \gamma_2 - 2 \sqrt{\gamma_1 \gamma_2} \tag{6.17}\]
\[\gamma_\mathrm{SL} \approx \gamma_\mathrm S + \gamma_\mathrm L - 2 \sqrt{\gamma_\mathrm S^\mathrm d \gamma_\mathrm L^\mathrm d} - 2 \sqrt{\gamma_\mathrm S^\mathrm p \gamma_\mathrm L^\mathrm p} \tag{6.18}\]
  • \(\gamma_\mathrm S^\mathrm d\), \(\gamma_\mathrm S^\mathrm p\) の上付き文字は, べきの指数でなくて分散, 極性を意味する添字.

6.3 重要な濡れの幾何配置

6.2.1 毛管上昇

\[h = \frac{2\gamma_\mathrm L \cos \theta}{r_\mathrm c g \rho} \tag{6.19}\]

具体例 6.5

  • 半径 \(5 \mu \mathrm m\) の完全濡れ性の毛管での水の毛管上昇高を求めよ.


解答例 6.5

\begin{align*} h & = \frac{(2)(0.072)}{(5 \times 10^{-6})(9.81)(997)} \\ & = 2.9446137 = 2.94 \ \mathrm{m} \tag{6.20} \end{align*}

式 (6.19) の導出

  • 毛管上昇により, 固気界面が固液界面に置き換わる. これによるギブズエネルギー変化は

\[\mathrm d G_{\Gamma} = - 2 \pi r_\mathrm c \, \mathrm d h \, (\gamma_\mathrm - \gamma_\mathrm{SL})\]
  • 毛管上昇により, (管断面積×高さ)分の水が重力に抗して上昇する. 微小高 \(\mathrm d h\) の変化に対応するギブズエネルギー変化は

\[\mathrm d G_\mathrm H = \pi r^2_\mathrm c \rho g h \, \mathrm d h\]
  • 併せて

\[\mathrm d G = \mathrm d G_\Gamma + \mathrm d G_\mathrm H = - 2 \pi r_\mathrm c (\gamma_\mathrm S - \gamma_\mathrm{SL}) + \pi r_\mathrm c^2 \rho g h \, \mathrm d h \tag{6.21}\]
  • 釣り合いから

\[\frac{\mathrm d G}{\mathrm d h } = - 2\pi r_\mathrm c (\gamma_\mathrm S - \gamma_\mathrm{SL}) + \pi r_\mathrm c \rho gh = 0\]
\[2(\gamma_\mathrm S - \gamma_\mathrm{SL}) = r_\mathrm c \rho g h \tag{6.22}\]
  • ヤングの式より, \(\gamma_\mathrm S - \gamma_\mathrm{SL} = \gamma_\mathrm L \cos \theta\) だから

\[2 \gamma_\mathrm L \cos \theta = r_\mathrm c \rho g h\]

毛管上昇の速度

  • Washburn の式(重力の影響を無視した場合)

\[l(t) = \sqrt{ \frac{r_\mathrm c \gamma_\mathrm Lt \cos \theta}{2 \eta}} \tag{6.23}\]

メニスカスの高さ

\[h = \sqrt{\frac{2\gamma_\mathrm L}{\rho g}(1 - \sin \theta)} \tag{6.24}\]

6.2.2 界面の粒子

\begin{align*} \pi( r^2 + h^2) & = \pi R^2 \left( \sin^2 \theta + 1 - 2 \cos \theta + \cos^2 \theta \right) \\ & = 2\pi R^2 ( 1 - \cos \theta) \tag{6.25} \end{align*}

\[\Delta G = 2 \pi R^2 (1 - \cos \theta) (\gamma_\mathrm{SL} - \gamma_\mathrm S) + \pi R^2 \gamma_\mathrm L \sin^2 \theta \tag{6.26}\]
  • 変形

\[\Delta G = \pi R^2 \gamma_\mathrm L (\cos \theta - 1)^2 \tag{6.27}\]
\[F = 2 \pi R \gamma_\mathrm L \sin^2 \frac{\theta}{2} \tag{6.28}\]

具体例 6.6

  • 水表面(\(\gamma_\mathrm L = 0.072 \ J \, \mathrm{J \, m^{-2}}\)), \(5 \, \mathrm{\mu m}\) の疎水性微小球(\(\theta = 90{}^\circ\))を取り除くのに必要な仕事と粘着力


解答例 6.6

  • (6.26), (6.27) の逆向きの操作であるが, 計算は同じ(符号のみ逆だが, この場合絶対値が求まればよい).

\begin{align*} \Delta G & = \pi R^2 \gamma_\mathrm L ( \cos \theta - 1) \\ & = \pi (5\times 10^{-6})^2 (0.072) (\cos \frac{\pi}{2}- 1)^2 \\ & = 5.654866 \times 10^{-12} = 5.7 \times 10^{-12} \ \mathrm{J} \end{align*}

6.3 接触角の計測

液滴の顕微鏡観察

\[\tan \left( \frac{\theta}{2}\right)= \frac{h}{a} \tag{6.29}\]

毛管上昇法

\[\Delta P_\mathrm{ref} = \frac{2 \gamma_\mathrm L^\mathrm{ref}}{\gamma_\mathrm c^\mathrm{eff}} \tag{6.30}\]
\[\Delta P = \frac{2 \gamma_\mathrm L \cos \theta}{\gamma_\mathrm c^\mathrm{eff}} \tag{6.31}\]

接触角測定での履歴

  • 前進接触角(advancing contact angle): 濡れ線が前進する直前の接触角

  • 後退接触角(receding contact angle): 濡れ線が後退する直前の接触角

  • 通常 \(\theta_\mathrm{adv}\)\(\theta_\mathrm{rec}\) よりも大きい.

  • \(\theta_\mathrm{adv} - \theta_\mathrm{rec}: \ 5 \sim 20 ^\circ\), さらに大きい場合もある.

接触角履歴の物理的起源

  • 表面粗さ

  • 表面の化学的・構造的不均一性

  • 溶解物質の 3 相接触線への吸着

  • 3 相接触線での表面張力の作用

\[\gamma = \frac{F}{l} = \frac{P l \delta}{l} = P \delta \tag{6.32}\]

表面の粗さと不純物

  • 表面粗さの影響

\[\cos \theta_\mathrm{app} = R_\mathrm{rough} \cos \theta \tag{6.33}\]
  • 不均一な平坦面での接触角

\[\cos \theta_\mathrm{app} f_1 \cos \theta_1 + f_2 \cos \theta_2 \tag{6.34}\]

具体例 6.8

  • パターニングされた表面での水滴の接触角.

  • 金の表面にヘキサデカンチオール \(\mathrm{HS(CH_2)_{15}CH_3}\) とジウンデカンジスルフィドカルボン酸 \(\mathrm{[S(CH_2)_{11}COOH]_2}\) を各々 2.5 と 3.0 \(\mu \mathrm{m}\) 幅で平行ストライプ状にパターニング.

  • 純粋なヘキサデカンチオール単分子膜(9.3.1 項)は疎水性で, pH 7.0 での前進接触角は \(107.8 {}^\circ\).

  • ジスルフィドはカルボン酸基が表面に出るため, より親水性である(前進接触角は \(50.1 {}^\circ\)

  • ストライプに平行な方向の有効接触角を求めよ.


解答例 6.8

\[\cos \theta_\mathrm{app} = \frac{2.5}{2.5 + 3.0} \cos (107.8^\circ) + \frac{3.0}{2.5 + 3.0} \cos (50.1^\circ)\]
  • 表面の粗さと不純物つづき

  • 表面構造の長さスケールと形状

  • 巨視的には接触角履歴, 微視的にはヤング方程式

  • 巨視的スケールでの接触角履歴の観測は必ずしもヤング方程式の破れを意味しない.

  • マイクロピラー表面での研究

超撥水性表面

  • 超撥水性表面(superrepellent surface)

  • 以下の 2 つの条件を満たす.

    • 水のみかけの接触角が 150° 以上である.

    • 基板状に水滴を置き, 基板を傾けた場合, 10° 以下で水滴がころがり落ちる.

[20]:
#### 6.8
import numpy as np

coste = (2.5) / (2.5 + 3.0) * np.cos (np.pi * 107.8 / 180) + (3.0) / (2.5 + 3.0) * np.cos(np.pi * 50.1 / 180)
theta = np.arccos(coste)
theta_deg = theta * 180
print('{:.7e}'.format(coste))
print('{:.7e}'.format(theta))
print('{:.7e}'.format(theta_deg))
2.1092921e-01
1.3582709e+00
2.4448876e+02
[21]:
#### 6.6
import numpy as np


r = 5e-6
gmm_l = 0.072
theta = np.pi/2
dg = np.pi * r**2 * gmm_l * (np.cos(theta) - 1)**2
f = 2 * np.pi * r * gmm_l * (np.sin(theta/2))**2
print('{:.7e}'.format(dg))
print('{:.7e}'.format(f))
5.6548668e-12
1.1309734e-06
[2]:
#### 6.5

h = 2 * 0.072 / 5e-6 / 9.81 /997
print('{:.7e}'.format(h))
2.9446137e+00
[1]:
#### 6.2

import numpy as np

kl = 30500 / 3 / 6.02e23 / 0.565e-9
print('{:.7e}'.format(kl))
2.9890532e-11