聚电解质微凝胶的泊松-玻尔兹曼-弗洛里理论 _psi

参考文献：J. Chem. Phys. 141,(2014)
聚电解质微凝胶由\(N\)条链组成，每条链链长为\(m\) ，链带电分率为\(\alpha\) ，电离出反离子数目为\(Z=fNm\) 。微凝胶半径为\(a\) 。溶液中平均每个凝胶占据的体积为\(\frac{4\pi}{3}R^3\)
体系自由能（以\(k_BT\)约化）为：
\begin{} \begin{split} F=& F_{ela}+\int f(\vec{r})\ d\vec{r}\\ =& F_{ela}+\int [f_{ele}(\vec{r})+f_{tr}(\vec{r})+f_{mix}(\vec{r})]\ d\vec{r} \end{split} \label{} \end{}
\eqref{}式中\(F_{ela}\)为链的熵弹性能：
\begin{} F_{ela}=\frac{3N}{2}\left [\left (\frac{\phi_0}{\phi}\right )^{2/3}-\ln\left (\frac{\phi_0}{\phi}\right ) \right ] \label{Fela} \end{}
其中\(\phi_0\)、\(\phi\)分别为凝胶在参考态和当前态的链节体积分数。
\eqref{}式中\(f_{ele}(\vec{r})\)为静电自由能密度：
\begin{} f_{ele}(\vec{r})=-\frac{1}{8\pi l_B}\mid \nabla \psi(\vec{r}) \mid^2+[n_+(\vec{r})-n_-(\vec{r})-\alpha\phi/v_0]\psi(\vec{r}) \label{fele} \end{}
其中\(\psi(\vec{r})\)为电势（以\(k_BT/e\)约化）， \(n_{\pm}\)为小离子的数密度。
\eqref{}式中\(f_{tr}(\vec{r})\)为小离子的平动熵：
\begin{} f_{tr}(\vec{r})=n_+(\vec{r})[\ln (n_+\^3)-1]+n_-(\vec{r})[\ln (n_-\^3)-1]-\mu_+n_+(\vec{r})-\mu_-n_-(\vec{r}) \label{ftr} \end{}
其中\(\\)为德布罗意热波长， \(\mu_{\pm}\)为正负小离子的化学势。
\eqref{}式中\(f_{mix}(\vec{r})\)为溶剂和凝胶的混合自由能：

文章插图
\begin{} f_{tr}(\vec{r})=\frac{1}{v_0}[(1-\phi)\ln (1-\phi)+\chi \phi(1-\phi)] \label{fmix} \end{}
其中\(v_0\)为溶剂分子和凝胶链节的体积。
\eqref{}式对\(\psi(\vec{r})\)变分，得
\begin{} \nabla^2 \psi(\vec{r})=-4\pi l_B[n_+(\vec{r})-n_-(\vec{r})-\alpha \phi/v_0] \label{PB} \end{}
【聚电解质微凝胶的泊松-玻尔兹曼-弗洛里理论】\eqref{}式对\(n_{\pm}(\vec{r})\)变分，并整理，得
\begin{} n_{\pm}(\vec{r})=\mu_{\pm}e^{\mp\psi(\vec{r})} \label{nden} \end{}
在本体溶液\(\psi=0\) ， \(n_{\pm}=c_s\) ，这里\(c_s\)为本体溶液盐浓度。所以上式可化为
\begin{} n_{\pm}(\vec{r})=c_se^{\mp\psi(\vec{r})} \label{ndis} \end{}
\eqref{Fela}式对\(V\)求导，得熵弹性力
\begin{} \Pi_{ela}=-\frac{\ dF_{ela}}{\ dV}=-\frac{N}{V_0}\left [\left ( \frac{\phi}{\phi_0}\right )^{1/3} -\frac{2}{3}\frac{\phi}{\phi_0} \right ] \label{Piela} \end{}
凝胶内小分子渗透压为
\begin{} \begin{split} \Pi_{g}=&n_+\frac{\ f}{\ n_+}+n_-\frac{\ f}{\ n_-}+\mid \nabla \psi\mid\frac{\ f}{\ \nabla \psi}+\phi\frac{\ f}{\ \phi}-f\\ =&-\frac{1}{8\pi l_B}\mid \nabla \psi(\vec{r}) \mid^2+n_+(\vec{r})+n_-(\vec{r}) \\ & -\frac{1}{v_0}\left [\ln(1-\phi)+\chi\phi^2+\phi\right ] \end{split} \label{Pig} \end{}
本体溶液小分子渗透压为
\begin{} \Pi_{b}=2c_s \label{Pib} \end{}
体系渗透压差为

文章插图
\begin{}
\begin{split}
\Pi_{os}=\Pi_{g}-\Pi_{b}=&-\frac{1}{8\pi l_B}\mid \nabla \psi(\vec{r}) \mid^2+n_+(\vec{r})+n_-(\vec{r})-2c_s \
& -\frac{1}{v_0}\left [\ln(1-\phi)+\chi\phi^2+\phi\right ]
\end{split}
\label{Pios}
\end{}
渗透压差与熵弹性力二者平衡：
\begin{} \begin{split} \Pi_{ela}+\Pi_{os}=&-\frac{N}{V_0}\left [\left ( \frac{\phi}{\phi_0}\right )^{1/3} -\frac{2}{3}\frac{\phi}{\phi_0} \right ]\\ &-\frac{1}{8\pi l_B}\mid \nabla \psi(\vec{r}) \mid^2+n_+(\vec{r})+n_-(\vec{r})-2c_s \\ & -\frac{1}{v_0}\left [\ln(1-\phi)+\chi\phi^2+\phi\right ]\\ =&0 \end{split} \label{Pieq} \end{}
此式中\(\vec{r}\)可取凝胶内部任意一点，显然取\(r=0\)处最为方便。
体系状态由\eqref{PB}、\eqref{ndis}和\eqref{Pieq}组成的方程组给出。体系具有球对称性，以上方程在球坐标系求解。要求解的方程即为：
\begin{*} \begin{cases} \frac{\ d^2}{\ dr^2}\psi(r)+\frac{2}{r}\frac{\ d}{\ dr}\psi(r)=-4\pi l_B[n_+(r)-n_-(r)-\alpha\phi/v_0]\\ \\ n_{\pm}(r)=c_se^{\mp\psi(r)}\\ \\ -\frac{N}{V_0}\left [\left ( \frac{\phi}{\phi_0}\right )^{1/3} -\frac{2}{3}\frac{\phi}{\phi_0} \right ]+n_+(0)+n_-(0)-2c_s\\ -\frac{1}{v_0}\left [\ln(1-\phi)+\chi\phi^2+\phi\right ]=0 \\ \phi=\frac{mNv_0}{\frac{4\pi a_3}{3}} \Theta(a-r) \end{cases} \end{*}
其中\(\Theta(x)\)为阶跃函数。
边界条件：
\begin{*} \begin{cases} \frac{\ d}{\ dr}\psi(r)\Big|_{r=0}=0 \\ \psi(r=a^-)=\psi(r=a^+)\\ \frac{\ d}{\ dr}\psi(r)\Big|_{r=R}=0 \end{cases} \end{*}