古典成核理論
![本頁使用了標題或全文手工轉換](http://images.weserv.nl/?url=//upload.wikimedia.org/wikipedia/commons/thumb/c/cd/Zh_conversion_icon_m.svg/35px-Zh_conversion_icon_m.svg.png)
古典成核理論(Classical nucleation theory,CNT)是用於研究成核現象的最常見的理論模型。[1][2]
概述
基本的古典成核理論可被用於粗略地估計新相態在一系列相同成核點的成核速率R。我們可以通過下述例子來理解成核速率R:例如,在0.1m3的潮濕空氣中每1秒形成100滴水,則成核速率R=1000/s。根據現代的古典成核理論,成核速率R可以用下式預測:
其中
- 是吉布斯能。是平均熱能。表示絕對溫度,kB為波茲曼常數;
- 為成核點的數量;
- 表示的是吸附到核上的分子參與到成核過程中的速率;
- 也被稱作澤爾多維奇因子Z,表示勢壘頂部的核繼續生成新相態而不分解的機率。
這一表達式可以被分為兩部分:第一部分[3]。在某一成核點形成核的機率正比於。因此,若較大且為正,成核的機率就較低且成核速度慢。
表達式中剩下的是動力學部分。這裡的表示的是物質輸入的速率,而表示臨界大小(處於能量勢壘極大值)的核會繼續成長的機率。處於勢壘頂部的核既可以繼續成長得更大而成為一種新相態,也可能開始失去分子而縮小。
以下是古典成核理論在實際中應用的一個例子。Sanz與其同事[4]使用電腦模擬液態水中冰的成核,估計了上式中所有的物理量。他們使用的水的物理模型是TIP4P/2005。當溫度降至模型中水的凝固點以下19.5 °C時,他們預估冰的成核的自由能勢壘為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)成核,即冰是與容器表面的接觸而成核。
同質成核
![](http://images.weserv.nl/?url=//upload.wikimedia.org/wikipedia/commons/thumb/e/ee/Plot_of_classical_nucleation_theory_prediction_for_the_free_energy_of_a_nucleus_as_a_function_of_radius.png/250px-Plot_of_classical_nucleation_theory_prediction_for_the_free_energy_of_a_nucleus_as_a_function_of_radius.png)
相比於異質成核,同質成核比較少見[5],但後者較為簡單且容易理解。下面介紹如何計算同質成核中的勢壘。
根據古典理論的假設[6],即使是對於新相態的微觀核,自由能也可以被寫作體項(bulk term)和表面項(surface term)相加的總和。其中體項正比於核的體積,而表面項正比於其表面積:
第一項即為體積項。我們假設核為球形,因此其體積相當於半徑為的球體。表示的是每一單位體積中,已成核的相態與正在成核的相態之間的自由能差距。例如,水在過飽和空氣中成核的過程中,計算的是每一單位體積中過飽和空氣的自由能與同一氣壓下的液態水的自由能之差。由於空氣只有在過飽和的情況下才會出現成核現象,所以為負數。第二項來源於核表面形成的界面,為界面上的表面張力,是一個正數。
對於較小的,表達式的值主要取決於第二項且;當較大時,項的影響越來越大,最終自由能變為負數。因此當取中間值時,自由能有一個極大值,即當時,核的半徑若繼續增大,則其自由能開始減小。讓自由能取得最大值的半徑被稱作臨界半徑。自由能的最大值可以由下式給出:
由於和溫度相關,這一理論也說明通過調節溫度可以控制成核的臨界半徑。上述的變量可以表示為:
其中為熔點,是生成物質的焓。因此臨界半徑可進一步被表示為:
因此在附近升高溫度時,臨界半徑也會隨之增加。
異質成核
不同於同質成核,異質成核發生於表面上。之所以異質成核較同質成核常見,是因為有各種因素的影響[7];異質成核在動力學上更容易出現,是因為異質成核的成核勢壘ΔG*在表面上較同質成核的來得低。對於同質成核,核的形狀是用球體近似估計的;而在異質成核中則需要考慮更多因素。
![](http://images.weserv.nl/?url=//upload.wikimedia.org/wikipedia/commons/thumb/2/2c/Surface_Interactions_of_a_Droplet_%28Heterogeneous_Nucleation%29.svg/220px-Surface_Interactions_of_a_Droplet_%28Heterogeneous_Nucleation%29.svg.png)
![](http://images.weserv.nl/?url=//upload.wikimedia.org/wikipedia/commons/thumb/0/00/Surface_tension.svg/450px-Surface_tension.svg.png)
如左圖所示,影響異質成核的因素有液滴的大小、接觸角的大小,以及三種不同界面之間的交互作用(液體表面與固體表面之間、固體表面和空氣之間和液體表面和空氣之間)。在右圖中,可以看到液滴不是完全的球形,所以表面積小於。[8]
右圖中液滴的接觸角從左到右依次減小,其表面積也隨之減小。表面積的減小降低了勢壘,提高的成核速率。對於簡單的幾何形狀已有相關的理論推導[9]。
古典成核理論、電腦模擬和實驗之間的對比
古典成核理論中有不少的假設,例如理論中將微觀的核當作宏觀的、具有良好定義的表面的液滴,其自由能可以用一熱力學平衡下的性質——表面張力σ來估計。然而對於很小的、僅有數十分子的核,很難說是否能夠繼續用體積項加上表面項來處理。另外成核過程是一種隱含的熱力學平衡現象,因此能否用熱力學平衡下的性質來近似這一過程也是個未知數。
對於簡單的模型,現代電腦的計算能力已經足以精確計算出成核速率。一個例子是硬球體模型中晶相的成核。古典成核理論估計的成核速率和電腦計算的結果相符,證明了古典成核理論是一個十分有用的近似理論[10]。對於簡單的模型,古典成核理論能給出比較好的結果,但是不一定能同樣用於複雜系統。Jones等人通過數值計算古典水模型研究了水分子簇的成核。他們發現古典成核理論對於8到50個分子的水分子簇可以給出比較好的近似,然而對於更小的水分子簇卻失效了[11]。雖然古典成核理論可以通過一些更高精度的方法(例如量子化學計算)進行修正[12],但是古典成核理論在估計例如氬氣蒸汽成核實驗的時候還是會與實際結果相差好幾個數量級[13]。
參考資料
- ^ 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.