查看“经典成核理论”的源代码

{{NoteTA
|G1=Physics
|G2=Chemistry}}

'''经典成核理论'''（Classical nucleation theory，'''CNT'''）是用于研究[[成核]]现象的最常见的理论模型。<ref>{{cite book |title=Microphysics of Clouds and Precipitation |publisher=Springer Science + Business Media B.V |isbn=978-0-7923-4211-3 |edition=1}}</ref><ref name="debenedetti_book">{{cite book |title=Metastable liquids : concepts and principles |publisher=Princeton University Press |isbn=9780691085951}}</ref>

== 概述 ==
基本的经典成核理论可被用于粗略地估计新相态在一系列相同成核点的成核速率R。我们可以通过下述例子来理解成核速率R：例如，在0.1m<sup>3</sup>的潮湿空气中每1秒形成100滴水，则成核速率R=1000/s。根据现代的经典成核理论，成核速率R可以用下式预测：

:: <math>R\ =\ N_S Zj\exp \left( \frac{-\Delta G^*}{k_BT} \right)</math>

其中

:* <math>\Delta G^*</math><span>T表示温度，</span>''k''<sub>B</sub>为[[波茲曼常數|玻尔兹曼常数]]；
:* <math>N_S</math>为成核点的数量；
:* <math>j</math>表示的是吸附到核上的分子参与到成核过程中的速率；
:* <math>Z</math>也被称作泽尔多维奇因子Z，表示势垒顶部的核继续生成新相态而不分解的概率。

这一表达式可以被分为两部分：第一部分<math> N_S \exp \left( -\Delta G^*/k_BT \right)</math><ref name="frenkel_book">{{Cite book|title=Understanding Molecular Simulation, Second Edition: From Algorithms to Applications|last=Frenkell|first=Daan|last2=Smit|first2=Berent|date=2001|isbn=0122673514|page=Academic Press}}</ref>。在某一成核点形成核的概率正比于<math> \exp[-\Delta G^*/kT]</math>。因此，若<math>\Delta G^*</math>较大且为正，成核的概率就较低且成核速度慢。

表达式中剩下的是动力学部分<math>Zj</math>。这里的<math> j</math>表示的是物质输入的速率，而<math>Z</math>表示临界大小（处于能量势垒极大值）的核会继续成长的概率。处于势垒顶部的核既可以继续成长得更大而成为一种新相态，也可能开始失去分子而缩小。

以下是经典成核理论在实际中应用的一个例子。Sanz与其同事<ref>{{Cite journal|title=Homogeneous Ice Nucleation at Moderate Supercooling from Molecular Simulation|last=Sanz|first=Eduardo|last2=Vega|first2=Carlos|journal=Journal of the American Chemical Society|issue=40|doi=10.1021/ja4028814|year=2013|volume=135|pages=15008–15017|arxiv=1312.0822|pmid=24010583|last3=Espinosa|first3=J. R.|last4=Cabellero-Bernal|first4=R.|last5=Abascal|first5=J. L. F.|last6=Valeriani|first6=Chantal}}</ref>使用电脑模拟液态水中冰的成核，估计了上式中所有的物理量。他们使用的水的物理模型是TIP4P/2005。当温度降至模型中水的凝固点以下19.5&#x202F;°C时，他们预估冰的成核的自由能势垒为<math> \Delta G^*=275k_BT</math>''j'' = 10<sup>11</sup>/s，且泽尔多维奇因子''Z'' = 10<sup>−3</sup>。1 m<sup>3</sup>水大概包含10<sup>28</sup>的水分子。将所有物理量代入公式可得，成核速率大概是10<sup>−83</sup>/s。这意味着在-20 °C的条件下，我们平均需要等待10<sup>83</sup>秒（10<sup>76</sup>年）才能看到一个冰核在体积为1 m<sup>3</sup>的水中出现！

需要注意的是，上述例子是水的一个理论模型中的同质（homogeneous）成核速率的预测，不是实际的水——在实验中，势垒''&#x394;G*''或动力学因子j都是无法直接从实际的水中测得的。然而，这一理论结果也表明-20&#x202F;°C左右条件下的同质成核速率之慢；因此实际生活中的水在-20&#x202F;°C或更高温度结冰主要是缘于异质（heterogeneous）成核，即冰是与容器表面的接触而成核。

=== 同质成核 ===
[[File:Plot_of_classical_nucleation_theory_prediction_for_the_free_energy_of_a_nucleus_as_a_function_of_radius.png|缩略图|250像素|左|绿色曲线表示总体自由能随半径的变化关系。自由能的势垒被标记为<math>\Delta G^*</math>。总自由能是体自由能（红色曲线）和界面自由能（黑色曲线）之和。]]
相比于异质成核，同质成核比较少见<ref name="sear_review14">{{Cite journal|title=Quantitative Studies of Crystal Nucleation at Constant Supersaturation: Experimental Data and Models|url=http://personal.ph.surrey.ac.uk/~phs1rs/crys_eng_review.pdf|last=Sear|first=Richard  P.|journal=CrystEngComm|issue=29|doi=10.1039/C4CE00344F|year=2014|volume=16|pages=6506–6522|format=PDF|author=|access-date=2018-12-12|archive-url=https://web.archive.org/web/20141025202306/http://personal.ph.surrey.ac.uk/~phs1rs/crys_eng_review.pdf|archive-date=2014-10-25|dead-url=yes}}</ref>，但后者较为简单且容易理解。下面介绍如何计算同质成核中的势垒<math>\Delta G^*</math>。

根据经典理论的假设<ref>{{cite book |title=Homogeneous nucleation theory; the pretransition theory of vapor condensation |publisher=Academic Press |isbn=978-0-12-038361-0}}</ref>，即使是对于新相态的微观核，自由能<math> \Delta G</math>也可以被写作体项（bulk term）和表面项（surface term）相加的总和。其中体项正比于核的体积，而表面项正比于其表面积：

: <math>\Delta G = \frac{4}{3} \pi r^3 \Delta g + 4 \pi r^2 \sigma</math>

第一项即为体积项。我们假设核为球形，因此其体积相当于半径为<math>r</math>的球体。<math>\Delta g</math>表示的是每一单位体积中，已成核的相态与正在成核的相态之间的自由能差距。例如，水在过饱和空气中成核的过程中，<math>\Delta g</math>计算的是每一单位体积中过饱和空气的自由能与同一气压下的液态水的自由能之差。由于空气只有在过饱和的情况下才会出现成核现象，所以<math>\Delta g</math>为负数。第二项来源于核表面形成的界面，<math>\sigma</math>为界面上的表面张力，是一个正数。

对于较小的<math>r</math>，表达式的值主要取决于第二项且<math>\Delta G(r)>0</math>；当<math>r</math>较大时，<math>r^3</math>项的影响越来越大，最终自由能变为负数。因此当<math>r</math>取中间值时，自由能有一个极大值，即当<math> \frac{dG}{dr} = 0 </math>时，核的半径若继续增大，则其自由能开始减小。让自由能取得最大值<math>\Delta G^*</math>的半径被称作临界半径。自由能的最大值可以由下式给出：

:: <math>\Delta G^* = \frac{16 \pi \sigma ^3}{3(\Delta g)^2}</math>

由于<math> \Delta G</math>和温度相关，这一理论也说明通过调节温度可以控制成核的临界半径。上述的变量<math> \Delta g</math>可以表示为：

: <math> \frac{H_f(T_m) - T}{T_m}</math>

其中<math> T_m </math>为熔点，<math>H_f</math>是生成物质的焓。因此临界半径可进一步被表示为：

: <math> r_c\ =\ \frac{-2\sigma T_m}{\Delta H_f}\frac{1}{T_m -T}</math>

因此在<math> T_m </math>附近升高温度时，临界半径也会随之增加。

=== 异质成核 ===
不同于同质成核，异质成核发生于表面上。之所以异质成核较同质成核常见，是因为有各种因素的影响<ref>{{Cite journal|title=Heterogeneous nucleation or homogeneous nucleation?|url=https://doi.org/10.1063/1.481644|last=Liu|first=X. Y.|date=31 May 2000|journal=The Journal of Chemical Physics|accessdate=16 November 2017|issue=22|doi=10.1063/1.481644|volume=112|pages=9949–9955|bibcode=2000JChPh.112.9949L|issn=0021-9606}}</ref>；异质成核在动力学上更容易出现，是因为异质成核的成核势垒&#x394;G*在表面上较同质成核的来得低。对于同质成核，核的形状是用球体近似估计的；而在异质成核中则需要考虑更多因素。
[[File:Surface_Interactions_of_a_Droplet_(Heterogeneous_Nucleation).svg|左|缩略图|图中标记了影响异质成核的各种因素。]]
[[File:Surface_tension.svg|右|450像素|缩略图|固体表面上的三滴液滴，其接触角从左到右依次减小。]]
如左图所示，影响异质成核的因素有液滴的大小、[[接觸角|接触角]]的大小，以及三种不同界面之间的相互作用（液体表面与固体表面之间、固体表面和空气之间和液体表面和空气之间）。在右图中，可以看到液滴不是完全的球形，所以表面积小于<math> 4\pi r^2 </math>。<ref name="sear_review">{{Cite journal|title=Nucleation: theory and applications to protein solutions and colloidal suspensions|url=http://iopscience.iop.org/0953-8984/19/3/033101|last=Sear|first=R. P.|journal=J. Phys.: Condens. Matter|issue=3|doi=10.1088/0953-8984/19/3/033101|year=2007|volume=19|page=033101|bibcode=2007JPCM...19c3101S}}</ref>

右图中液滴的接触角从左到右依次减小，其表面积也随之减小。表面积的减小降低了势垒，提高的成核速率。对于简单的几何形状已有相关的理论推导<ref>{{Cite journal|title=Decoration criteria for surface steps|url=http://www.sciencedirect.com/science/article/pii/0001616070900064|last=Sholl|first=C. A.|last2=N. H. Fletcher|journal=Acta Metall.|issue=10|doi=10.1016/0001-6160(70)90006-4|year=1970|volume=18|pages=1083–1086}}</ref>。

[[File:Hethomnucdifference.JPG|右|缩略图|不同类型的成核表现出不同的能量势垒。]]

== 经典成核理论、电脑模拟和实验之间的对比 ==
经典成核理论中有不少的假设，例如理论中将微观的核当作宏观的、具有良好定义的表面的液滴，其自由能可以用一[[热力学平衡]]下的性质——表面张力σ来估计。然而对于很小的、仅有数十分子的核，很难说是否能够继续用体积项加上表面项来处理。另外成核过程是一种隐含的热力学平衡现象，因此能否用热力学平衡下的性质来近似这一过程也是个未知数。

对于简单的模型，现代电脑的计算能力已经足以精确计算出成核速率。一个例子是{{le|硬球体模型|hard spheres}}中晶相的成核。经典成核理论估计的成核速率和电脑计算的结果相符，证明了经典成核理论是一个十分有用的近似理论<ref name="auer04">{{Cite journal|title=Numerical prediction of absolute crystallization rates in hard-sphere colloids|url=http://scitation.aip.org/content/aip/journal/jcp/120/6/10.1063/1.1638740|last=Auer|first=S.|last2=D. Frenkel|journal=J. Chem. Phys.|issue=6|doi=10.1063/1.1638740|year=2004|volume=120|page=3015|bibcode=2004JChPh.120.3015A}}</ref>。对于简单的模型，经典成核理论能给出比较好的结果，但是不一定能同样用于复杂系统。Jones等人通过数值计算经典{{le|水模型|water model}}研究了[[水分子簇]]的成核。他们发现经典成核理论对于8到50个分子的水分子簇可以给出比较好的近似，然而对于更小的水分子簇却失效了<ref>{{Cite journal|title=Origin of the Failure of Classical Nucleation Theory: Incorrect Description of the Smallest Clusters|last=Merikanto|first=Joonas|last2=Zapadinsky|first2=Evgeni|date=4 April 2007|journal=Physical Review Letters|issue=14|doi=10.1103/PhysRevLett.98.145702|volume=98|bibcode=2007PhRvL..98n5702M|last3=Lauri|first3=Antti|last4=Vehkamäki|first4=Hanna}}</ref>。虽然经典成核理论可以通过一些更高精度的方法（例如[[量子化学]]计算）进行修正<ref>{{Cite journal|title=Quantum Mechanical Study of Sulfuric Acid Hydration: Atmospheric Implications|last=Temelso|first=Berhane|last2=Morrell|first2=Thomas E.|date=22 February 2012|journal=The Journal of Physical Chemistry A|issue=9|doi=10.1021/jp2119026|volume=116|pages=2209–2224|bibcode=2012JPCA..116.2209T|last3=Shields|first3=Robert M.|last4=Allodi|first4=Marco A.|last5=Wood|first5=Elena K.|last6=Kirschner|first6=Karl N.|last7=Castonguay|first7=Thomas C.|last8=Archer|first8=Kaye A.|last9=Shields|first9=George C.}}</ref>，但是经典成核理论在估计例如[[氩气]]蒸汽成核实验的时候还是会与实际结果相差好几个数量级<ref>{{cite journal |last1=Fladerer |first1=Alexander |last2=Strey |first2=Reinhard |title=Homogeneous nucleation and droplet growth in supersaturated argon vapor: The cryogenic nucleation pulse chamber |journal=The Journal of Chemical Physics |date=2006-04-28 |volume=124 |issue=16 |pages=164710 |doi=10.1063/1.2186327}}</ref>。

== 参考资料 ==
{{reflist|2}}

[[Category:粒子探测器]]
[[Category:自组织]]
[[Category:物理化学]]