经典成核理论

经典成核理论（Classical nucleation theory，CNT）是用于研究成核现象的最常见的理论模型。 基本的经典成核理论可被用于粗略地估计新相态在一系列相同成核点的成核速率R。我们可以通过下述例子来理解成核速率R：例如，在0.1m3的潮湿空气中每1秒形成100滴水，则成核速率R=1000/s。根据现代的经典成核理论，成核速率R可以用下式预测： R   =   N S Z j exp ⁡ ( − Δ G ∗ k B T ) {\displaystyle R\ =\ N_{S}Zj\exp \left({\frac {-\Delta G^{*}}{k_{B}T}}\right)} 其中 Δ G ∗ {\displaystyle \Delta G^{*}} 是吉布斯能。 k B T {\displaystyle k_{B}T} 是平均熱能。 T {\displaystyle T} 表示絕對温度，kB为玻尔兹曼常数； N S {\displaystyle N_{S}} 为成核点的数量； j {\displaystyle j} 表示的是吸附到核上的分子参与到成核过程中的速率； Z {\displaystyle Z} 也被称作泽尔多维奇因子Z，表示势垒顶部的核继续生成新相态而不分解的概率。 这一表达式可以被分为两部分：第一部分 N S exp ⁡ ( − Δ G ∗ / k B T ) {\displaystyle N_{S}\exp \left(-\Delta G^{*}/k_{B}T\right)} 。在某一成核点形成核的概率正比于 exp ⁡ [ − Δ G ∗ / k T ] {\displaystyle \exp[-\Delta G^{*}/kT]} 。因此，若 Δ G ∗ {\displaystyle \Delta G^{*}} 较大且为正，成核的概率就较低且成核速度慢。 表达式中剩下的是动力学部分 Z j {\displaystyle Zj} 。这里的 j {\displaystyle j} 表示的是物质输入的速率，而 Z {\displaystyle Z} 表示临界大小（处于能量势垒极大值）的核会继续成长的概率。处于势垒顶部的核既可以继续成长得更大而成为一种新相态，也可能开始失去分子而缩小。 以下是经典成核理论在实际中应用的一个例子。Sanz与其同事使用电脑模拟液态水中冰的成核，估计了上式中所有的物理量。他们使用的水的物理模型是TIP4P/2005。当温度降至模型中水的凝固点以下19.5 °C时，他们预估冰的成核的自由能势垒为 Δ G ∗ = 275 k B T {\displaystyle \Delta G^{*}=275k_{B}T} j = 1011/s，且泽尔多维奇因子Z = 10−3。1 m3水大概包含1028的水分子。将所有物理量代入公式可得，成核速率大概是10−83/s。这意味着在-20 °C的条件下，我们平均需要等待1083秒（1076年）才能看到一个冰核在体积为1 m3的水中出现！ 需要注意的是，上述例子是水的一个理论模型中的同质（homogeneous）成核速率的预测，不是实际的水——在实验中，势垒ΔG*或动力学因子j都是无法直接从实际的水中测得的。然而，这一理论结果也表明-20 °C左右条件下的同质成核速率之慢；因此实际生活中的水在-20 °C或更高温度结冰主要是缘于异质（heterogeneous）成核，即冰是与容器表面的接触而成核。 相比于异质成核，同质成核比较少见，但后者较为简单且容易理解。下面介绍如何计算同质成核中的势垒 Δ G ∗ {\displaystyle \Delta G^{*}} 。 根据经典理论的假设，即使是对于新相态的微观核，自由能 Δ G {\displaystyle \Delta G} 也可以被写作体项（bulk term）和表面项（surface term）相加的总和。其中体项正比于核的体积，而表面项正比于其表面积： Δ G = 4 3 π r 3 Δ g + 4 π r 2 σ {\displaystyle \Delta G={\frac {4}{3}}\pi r^{3}\Delta g+4\pi r^{2}\sigma } 第一项即为体积项。我们假设核为球形，因此其体积相当于半径为 r {\displaystyle r} 的球体。 Δ g {\displaystyle \Delta g} 表示的是每一单位体积中，已成核的相态与正在成核的相态之间的自由能差距。例如，水在过饱和空气中成核的过程中， Δ g {\displaystyle \Delta g} 计算的是每一单位体积中过饱和空气的自由能与同一气压下的液态水的自由能之差。由于空气只有在过饱和的情况下才会出现成核现象，所以 Δ g {\displaystyle \Delta g} 为负数。第二项来源于核表面形成的界面， σ {\displaystyle \sigma } 为界面上的表面张力，是一个正数。 对于较小的 r {\displaystyle r} ，表达式的值主要取决于第二项且 Δ G ( r ) > 0 {\displaystyle \Delta G(r)>0} ；当 r {\displaystyle r} 较大时， r 3 {\displaystyle r^{3}} 项的影响越来越大，最终自由能变为负数。因此当 r {\displaystyle r} 取中间值时，自由能有一个极大值，即当 d G d r = 0 {\displaystyle {\frac {dG}{dr}}=0} 时，核的半径若继续增大，则其自由能开始减小。让自由能取得最大值 Δ G ∗ {\displaystyle \Delta G^{*}} 的半径被称作临界半径。自由能的最大值可以由下式给出： Δ G ∗ = 16 π σ 3 3 ( Δ g ) 2 {\displaystyle \Delta G^{*}={\frac {16\pi \sigma ^{3}}{3(\Delta g)^{2}}}} 由于 Δ G {\displaystyle \Delta G} 和温度相关，这一理论也说明通过调节温度可以控制成核的临界半径。上述的变量 Δ g {\displaystyle \Delta g} 可以表示为： H f ( T m − T ) T m {\displaystyle {\frac {H_{f}(T_{m}-T)}{T_{m}}}} 其中 T m {\displaystyle T_{m}} 为熔点， H f {\displaystyle H_{f}} 是生成物质的焓。因此临界半径可进一步被表示为： r c   =   − 2 σ T m Δ H f 1 T m − T {\displaystyle r_{c}\ =\ {\frac {-2\sigma T_{m}}{\Delta H_{f}}}{\frac {1}{T_{m}-T}}} 因此在 T m {\displaystyle T_{m}} 附近升高温度时，临界半径也会随之增加。 不同于同质成核，异质成核发生于表面上。之所以异质成核较同质成核常见，是因为有各种因素的影响；异质成核在动力学上更容易出现，是因为异质成核的成核势垒ΔG*在表面上较同质成核的来得低。对于同质成核，核的形状是用球体近似估计的；而在异质成核中则需要考虑更多因素。 如左图所示，影响异质成核的因素有液滴的大小、接触角的大小，以及三种不同界面之间的相互作用（液体表面与固体表面之间、固体表面和空气之间和液体表面和空气之间）。在右图中，可以看到液滴不是完全的球形，所以表面积小于 4 π r 2 {\displaystyle 4\pi r^{2}} 。 右图中液滴的接触角从左到右依次减小，其表面积也随之减小。表面积的减小降低了势垒，提高的成核速率。对于简单的几何形状已有相关的理论推导。 经典成核理论中有不少的假设，例如理论中将微观的核当作宏观的、具有良好定义的表面的液滴，其自由能可以用一热力学平衡下的性质——表面张力σ来估计。然而对于很小的、仅有数十分子的核，很难说是否能够继续用体积项加上表面项来处理。另外成核过程是一种隐含的热力学平衡现象，因此能否用热力学平衡下的性质来近似这一过程也是个未知数。 对于简单的模型，现代电脑的计算能力已经足以精确计算出成核速率。一个例子是硬球体模型（英语：hard spheres）中晶相的成核。经典成核理论估计的成核速率和电脑计算的结果相符，证明了经典成核理论是一个十分有用的近似理论。对于简单的模型，经典成核理论能给出比较好的结果，但是不一定能同样用于复杂系统。Jones等人通过数值计算经典水模型（英语：water model）研究了水分子簇的成核。他们发现经典成核理论对于8到50个分子的水分子簇可以给出比较好的近似，然而对于更小的水分子簇却失效了。虽然经典成核理论可以通过一些更高精度的方法（例如量子化学计算）进行修正，但是经典成核理论在估计例如氩气蒸汽成核实验的时候还是会与实际结果相差好几个数量级。 Microphysics of Clouds and Precipitation 1. Springer Science + Business Media B.V. ISBN 978-0-7923-4211-3.  Metastable liquids : concepts and principles. Princeton University Press. 1996. ISBN 9780691085951.  Frenkell, Daan; Smit, Berent. Understanding Molecular Simulation, Second Edition: From Algorithms to Applications. 2001: Academic Press. ISBN 0122673514.  Sanz, Eduardo; Vega, Carlos; Espinosa, J. R.; Cabellero-Bernal, R.; Abascal, J. L. F.; Valeriani, Chantal. Homogeneous Ice Nucleation at Moderate Supercooling from Molecular Simulation. Journal of the American Chemical Society. 2013, 135 (40): 15008–15017. PMID 24010583. arXiv:1312.0822 . doi:10.1021/ja4028814.  Sear, Richard P. Quantitative Studies of Crystal Nucleation at Constant Supersaturation: Experimental Data and Models (PDF). CrystEngComm. 2014, 16 (29): 6506–6522 [2018-12-12]. doi:10.1039/C4CE00344F. （原始内容 (PDF)存档于2014-10-25）.  Homogeneous nucleation theory; the pretransition theory of vapor condensation. Academic Press. ISBN 978-0-12-038361-0.  Liu, X. Y. Heterogeneous nucleation or homogeneous nucleation?. The Journal of Chemical Physics. 31 May 2000, 112 (22): 9949–9955 [16 November 2017]. Bibcode:2000JChPh.112.9949L. ISSN 0021-9606. doi:10.1063/1.481644.  Sear, R. P. Nucleation: theory and applications to protein solutions and colloidal suspensions. J. Phys.: Condens. Matter. 2007, 19 (3): 033101. Bibcode:2007JPCM...19c3101S. doi:10.1088/0953-8984/19/3/033101.  Sholl, C. A.; N. H. Fletcher. Decoration criteria for surface steps. Acta Metall. 1970, 18 (10): 1083–1086 [2018-12-12]. doi:10.1016/0001-6160(70)90006-4. （原始内容存档于2015-09-24）.  Auer, S.; D. Frenkel. Numerical prediction of absolute crystallization rates in hard-sphere colloids. J. Chem. Phys. 2004, 120 (6): 3015 [2018-12-12]. Bibcode:2004JChPh.120.3015A. doi:10.1063/1.1638740. （原始内容存档于2014-03-27）.  Merikanto, Joonas; Zapadinsky, Evgeni; Lauri, Antti; Vehkamäki, Hanna. Origin of the Failure of Classical Nucleation Theory: Incorrect Description of the Smallest Clusters. Physical Review Letters. 4 April 2007, 98 (14). Bibcode:2007PhRvL..98n5702M. doi:10.1103/PhysRevLett.98.145702.  Temelso, Berhane; Morrell, Thomas E.; Shields, Robert M.; Allodi, Marco A.; Wood, Elena K.; Kirschner, Karl N.; Castonguay, Thomas C.; Archer, Kaye A.; Shields, George C. Quantum Mechanical Study of Sulfuric Acid Hydration: Atmospheric Implications. The Journal of Physical Chemistry A. 22 February 2012, 116 (9): 2209–2224. Bibcode:2012JPCA..116.2209T. doi:10.1021/jp2119026.  Fladerer, Alexander; Strey, Reinhard. Homogeneous nucleation and droplet growth in supersaturated argon vapor: The cryogenic nucleation pulse chamber. The Journal of Chemical Physics. 2006-04-28, 124 (16): 164710. doi:10.1063/1.2186327.