小蜗熊的蜂蜜罐
塔板理论简介及色谱流出曲线方程的推导
发布于: 2020-04-29 更新于: 2020-05-25 分类于: 杂谈 阅读次数: 

塔板理论是色谱分析中的一种理论,可以用于解释色谱曲线的形状及其影响因素。塔板理论的思路是将色谱柱类比为蒸馏塔,将色谱分离过程类比为蒸馏过程。具体过程是把色谱柱分成许多小段(即理论塔板),每个理论塔板内都存在固定相与流动相。经过多次分配平衡,不同分配系数的组分会先后分离出来。根据塔板理论,待分离组分流出色谱柱时的浓度与时间之间符合二项分布,当塔板数很高时,二项分布趋向于正态分布。这一曲线方程称作色谱流出曲线方程。

塔板理论的基本假设

  • 流动相进入色谱柱不是连续的而是脉动式的,每次进入量未一个塔板体积
  • 初始时所有组分存在与0号塔板上
  • 组分在两相内瞬间达到平衡
  • 塔板间无纵向扩散(柱向扩散)
  • 所有塔板上的分配系数为常数

色谱流出曲线方程的推导

流动相在通过色谱柱的过程可以视作一个塔板体积的单向多次运动,每次运动中某一分子可能处于流动相或是固定相。

在塔板上某一分子出现在流动相中的概率等于流动相中分子数与全部分子数之比,即:

$$ P_m=\frac{M_m}{M_m+M_s}=\frac{1}{1+\frac{M_s}{M_m}}=\frac{1}{1+\kappa} \tag{1} $$
其中$\kappa$是平衡时组分在固定相和流动相中的物质的量之比,因此分子出现在固定相中的概率为:
$$ P_s=1-P_m=\frac{\kappa}{1+\kappa} $$
设流动相运动r次,若要使某分子出现在第n个塔板上,需要在这r次运动中有n次位于流动相中,r-n次位于固定相中。概率为:
$$ P=P_{m}^{n}\cdot(1-P_m)^{r-n}=\frac{\kappa^{r-n}}{(1+\kappa)^r} $$
对于n和r存在以下关系:
$$ n=\frac{L}{H} \tag{2} $$ $$ r=\frac{V}{Hqw} \tag{3} $$
其中L为色谱柱长度,H为塔板间隔长度,V为通过色谱柱的流动相体积,q为柱内横截面积,w为流动相横截面积所占比例。

在r次运动中n次位于流动相相当于一个排列组合问题,所以某分子出现在第n个塔板上的概率为

$$ P(r,n)=P \cdot C^n_r=\frac{\kappa^{r-n} \cdot r!}{(1+\kappa)^r \cdot n! \cdot (r-n)!} $$
将上式两侧同时取对数可得
$$ \ln P(r,n)=n\ln P_{m}+(r-n)\ln (1-P_m)+\ln (r!)-\ln (n!)-\ln (r-n!) $$
使用斯特林公式替换掉阶乘后可得:
$$ \ln P(r,n)=n\ln P_{m}+(r-n)\ln (1-P_m)+r\ln r-n\ln n-(r-n)\ln (r-n) $$
对r求偏导数可得:
$$ \frac{\partial \ln P(r,n)}{\partial r}=\ln r-\ln (r-n)+\ln (1-P_m) $$
令偏导数为0,可得在$r=n/P_m$处曲线取得极值 继续对$\ln P(r,n)$取二阶偏导数可得:
$$ \frac{\partial^2 \ln P(r,n)}{\partial r^2}=\frac{1}{r}-\frac{1}{r-n} $$
在极值点处二阶导数取值为:
$$ \frac{\partial^2 \ln P(r,n)}{\partial r^2} \Bigg|_{r=n/P_m} = \frac{P_m^2}{n(1-P_m)} $$
将$\ln P(r,n)$在极值点$r=n/P_m$处进行二阶泰勒展开可得:
$$ \begin{split} \ln P(r,n) &= \ln P(\frac{n}{P_m} ,n) + \frac{r-\frac{n}{P_m}}{1!} \cdot \frac{\partial \ln P(r,n)}{\partial r} + \frac{(r-\frac{n}{P_m})^2}{2!} \cdot \frac{\partial^2 \ln P(r,n)}{\partial r^2} + o \left[\left(r-\frac{n}{P_m}\right)^3 \right] \\\\ &= \ln P(\frac{n}{P_m},n) + 0 - \frac{(r-\frac{n}{P_m})^2}{2!} \cdot \frac{P_m^2}{n(1-P_m)} + o \left[\left(r-\frac{n}{P_m}\right)^3 \right] \\\\ \end{split} $$
两侧同时取e的对数可得:
$$ P(r,n) = P(\frac{n}{P_m} ,n) \cdot exp\left[ - \frac{(rP_m-n)^2}{2n(1-P_m)} \right] $$
将$(1)(2)(3)$式代入上式,同时令$V_R=nHqw(1+\kappa)$,得:
$$ P(r,n) = P(\frac{n}{P_m} ,n) \cdot exp\left[ - \frac{(V-V_R)^2}{2HL\kappa(1+\kappa)q^2w^2} \right] \tag{*} $$
其中$V_R$为保留体积,指柱后浓度达到最大时流过的流动相体积 由于体积无穷大时一定会出现:
$$ \int_{0}^{+\infty}P(r,n)dV=1 $$
带入$(*)$式可得:
$$ \int_{0}^{+\infty}P(\frac{n}{P_m},n)\cdot exp\left[ - \frac{(V-V_R)^2}{2HL\kappa(1+\kappa)q^2w^2} \right]dV=P(\frac{n}{P_m},n)\cdot\sqrt{2\pi HLq^2w^2\kappa(1+\kappa)}=1 $$
整理后再与$(*)$式联立,并用浓度替换单个分子的概率可得:
$$ c=\frac{m}{\sqrt{2\pi HLq^2w^2\kappa(1+\kappa)}}exp\left[ - \frac{(V-V_R)^2}{2HL\kappa(1+\kappa)q^2w^2} \right] \tag{**} $$
或者写成下面的形式:
$$ c=\frac{m\sqrt{n}}{V_R\sqrt{2\pi}}exp\left[ - \frac{n(V-V_R)^2}{2{V_R}^2} \right] \tag{***} $$
从$(***)$式的形式中可以看出,色谱流出曲线是以$V_R$为对称轴的正态分布曲线。
--- 本文结束 The End ---