席航波+張海波

摘 要 晶體的熔化一直是凝聚態(tài)物理中廣泛關(guān)注的研究方向，盡管有上百年的研究歷史，但是仍然沒有一個完整的基礎(chǔ)層面上的理論。大量實驗和計算模擬研究工作表明，熔化符合KTHNY理論。晶體在熔化過程中的結(jié)構(gòu)和動力學(xué)也已成為大家關(guān)注的熱點。

關(guān)鍵詞 二維熔化；KTHNY；膠體；計算機模擬

中圖分類號：TF841 文獻標(biāo)識碼：A 文章編號：1671-7597（2014）07-0067-02

Research Progress on Two-Dimensional Melting

Xi Hangbo Zhang Haibo

（Center for Soft Condensed Matter Physics and Interdisciplinary Research， Soochow University， Suzhou 215006， China）

Abstract： Crystal melting has been widely studied in condensed matter physics for more than 100 years， but a complete theory about their fundamental mechanisms still remains to be answered. Over the years， a number of experiments and computer simulations have been consistent with KTHNY theory. Crystal structureand dynamicsin the meltingprocesshasbecome afocus of attention.

Key words： Two-dimensional melting KTHNY Colloid Computer simulation

熔化是自然界中一種常見的現(xiàn)象，通?？醋鞴獭旱囊患壪嘧?，這個相變過程，在凝聚態(tài)物理、材料物理以及納米技術(shù)中具有十分重要的研究意義。盡管有上百年的研究歷史，但是仍然缺少一個完整的基礎(chǔ)層面上的理論。在早期的熔化理論中，空間維度決定熔化行為，長波漲落對低維度晶體熔化的影響非常大，任何一個很小的熱擾動就會破壞一維和二維晶格的長程平移有序性，從而破壞一維和二維晶體。然而隨后的研究發(fā)現(xiàn)，雖然在二維晶體中，不存在長程的平移有序性，卻可以存在長程的取向有序性。二維晶體由于處在臨界維度，其熔化機制相比于三維時的情況又有所不同。目前比較普遍被人們接受的一個理論是KTHNY理論，它是最早由Kosterlitz和Thouless提出，經(jīng)過Halperin，Nelson和Young三人的完善最終形成的一個理論[1]。KTHNY理論指出二維熔化不是通常認為的固—液一級相變，而是由固相—中間相—液相這樣的一個連續(xù)的二級相變。這里的中間相也叫六角相（hexatic相），具有短程的平移有序性和準(zhǔn)長程的取向有序性。缺陷分離在KTHNY熔化理論中起到非常重要的作用，通常認為位錯（dislocation）的分離破壞平移序，固相熔化進入六角相，而向錯（disclination）的分離破壞取向序，六角相熔化進入液相。此外，二維熔化的主要理論還有一級相變理論，以及晶界熔化理論等。

本文主要從近年來國內(nèi)外關(guān)于二維晶體熔化的研究方法以及研究熱點這兩個方面，介紹一些二維晶體熔化的實驗和計算模擬工作，總結(jié)二維晶體熔化的研究進展。

1 實驗研究

熔化過程涉及大量粒子，原子、分子晶體在長度尺度以及時間尺度上，都不適合實驗觀察，然而，膠體體系在這方面卻有得天獨厚的優(yōu)勢。膠體（colloid）是一種不連續(xù)介質(zhì)分散在另一種連續(xù)介質(zhì)的均勻混合物。膠體粒子的大小在10 nm到10 μm之間，在溶液中有較強的布朗運動。大小均勻的微米級膠體小球，懸浮在溶液中，通過改變粒子的體積分數(shù)，可以形成液相，玻璃態(tài)相，固相，可以用來模擬原子、分子的相變行為。通過光學(xué)顯微鏡可直接觀察到膠體粒子的熱運動，再經(jīng)過圖像處理得到粒子的實時運動軌跡，從而進行各種定量測量。膠體作為原子分子的模型系統(tǒng)已經(jīng)被廣泛認可。

Han等人利用微米大小的熱敏膠體NIPA（N-isopropyl acrylamide）球懸浮液，封入兩層玻璃蓋玻片之間，形成單層的樣品，研究二維膠體晶體的熔化[2]。這種微凝膠球是由NIPA高分子鏈通過交聯(lián)形成，這種網(wǎng)狀的鏤空的結(jié)構(gòu)使其內(nèi)部含有大量的水。當(dāng)溫度升高，NIPA球收縮將水?dāng)D出，直徑會變小。在一個特定的溫度段，NIPA球的直徑與溫度呈線性關(guān)系，從而可以調(diào)節(jié)溫度來精確地改變體積分數(shù)，驅(qū)動體系相變。研究發(fā)現(xiàn)這種NIPA球的相互作用勢是短程軟排斥，其相圖與硬球的相圖十分相近。通過定量計算平移關(guān)聯(lián)和取向關(guān)聯(lián)函數(shù)，發(fā)現(xiàn)六角相這一中間相的存在，與KTHNY理論預(yù)測的一致。

Zheng和Grieve則通過另一個體系同樣發(fā)現(xiàn)六角相的存

在[3]。他們采用毫米級的鋼球放在金屬基板和透明的導(dǎo)電板之間，加上電壓后，讓鋼球之間產(chǎn)生庫倫相互作用，整個基板掛在數(shù)根彈簧上，這樣基板能夠自由地垂直運動，基板上附加一個電動馬達，驅(qū)動偏心輪來給整個系統(tǒng)補充動能。這一系統(tǒng)用來模擬真空中的二維原子晶體以及帶電粒子。通過調(diào)節(jié)庫侖力，可以改變整個體系的狀態(tài)，他們發(fā)現(xiàn)六角相中的缺陷密度與剪切阻率是連續(xù)變化，整個體系的熔化行為符合KTHNY理論預(yù)測。

2 計算機模擬

實驗方法存在一定的局限性，在有些情況下，以實驗手段難以完成，這時就采用計算機模擬來有效地解決這些問題。本文主要介紹兩種常用的計算機模擬方法：蒙特卡羅（Monte Carlo）和分子動力學(xué)（Molecular Dynamics）。endprint

2.1 蒙特卡羅

MC是以概率統(tǒng)計理論為指導(dǎo)的一種數(shù)值模擬方法，通常將粒子移動一個隨機分布的距離，再通過比較移動前后的能量變化來判斷是否接受這個移動，在判斷過程中需要比較隨機數(shù)與系統(tǒng)的玻爾茲曼因子，直到最后獲得能量最低的構(gòu)型。理論上MC更加靈活，可以通過修改相應(yīng)的統(tǒng)計權(quán)重來實現(xiàn)不同的系宗。MC能夠更真實的描述粒子的隨機運動，原則上比MD具有更短的平衡時間，但缺點是無法獲得動力學(xué)信息。

早期的MC模擬，通常研究二維硬球（hard disks）體系的熔化行為，其結(jié)果更微弱偏向一級相變。Lee和Strandburg采用恒壓MC模擬研究二維硬球晶體的熔化，發(fā)現(xiàn)系統(tǒng)的體積漲落分布出現(xiàn)雙峰結(jié)構(gòu)，在轉(zhuǎn)變點附近，自由能勢壘隨系統(tǒng)大小增加，表明系統(tǒng)是一級相變[4]。

Janke和Kleinhert采用MC模擬研究缺陷模型以及粒子間的彈性相互作用，發(fā)現(xiàn)作用勢的“軟硬”以及長度尺度對熔化機制有很大的影響[5]。通常長程的軟作用勢體系熔化是連續(xù)相變，遵循KTHNY理論，而短程的硬作用勢體系熔化是一級相變。二維硬球熔化是一級相變還是連續(xù)相變?nèi)源嬖跔幾h，這還需要更多的大體系模擬來進行驗證。

2.2 分子動力學(xué)

MD主要是求解牛頓運動方程，以粒子的初始位置和速度，求解新的粒子位置和速度，不用的算法需要不同的初始條件。在MD中，模型的設(shè)定十分重要，首先確定合理的勢函數(shù)來描述物理過程中粒子間的相互作用，然后需要一個能量較低的初始構(gòu)型，賦予各個粒子初始速度，這個速度要符合玻爾茲曼統(tǒng)計，接著調(diào)整各個粒子的速度，要保證各個方向上的動量總和為零，最后給系統(tǒng)加上合適的邊界條件。MD是從系統(tǒng)相空間中抽取樣本進行統(tǒng)計分子動力學(xué)計算，抽樣間隔就是時間步長。時間步長的選取關(guān)系到運動方程的數(shù)值積分，太長會丟失精度，太短則會增加模擬時間。MD的優(yōu)勢在于可以研究體系真實的動力學(xué)。

Chen等人采用Lennard-Jones作用勢[6]。L-J作用勢同硬球作用勢一樣，也是一個經(jīng)典作用勢，常用來模擬兩個分子之間的相互作用，特別是在描述惰性氣體分子間的相互作用時尤為精確。L-J是一個二體作用勢模型，它包含短程的排斥和長程的吸引兩個部分，由于作用勢含有長程部分，其結(jié)果更偏向連續(xù)相變。Chen分析了L-J體系焓的變化，發(fā)現(xiàn)在熔化過程中，存在穩(wěn)定的中間相，結(jié)果符合KTHNY理論。

Shiba等人關(guān)注了L-J體系的結(jié)構(gòu)和動力學(xué)，分析了熔化過程中體系的無序度、結(jié)構(gòu)因子，均方位移和粒子構(gòu)型的變化，發(fā)現(xiàn)在二維熔化中，體系的結(jié)構(gòu)和動力學(xué)均呈現(xiàn)異質(zhì)性[7]。Zhang等人發(fā)現(xiàn)了熔化過程中的存在鏈狀的協(xié)同運動粒子團簇，與玻璃體系中的協(xié)同重排區(qū)域相似[8]。分析體系的均方位移和結(jié)構(gòu)弛豫，同樣發(fā)現(xiàn)動力學(xué)呈異質(zhì)性。此外，他們發(fā)現(xiàn)了間隙缺陷往往會誘發(fā)這種鏈狀的協(xié)同運動。

3 結(jié)束語

二維熔化的物理實質(zhì)目前尚未完全清楚，二維熔化理論仍然存在爭議，KTHNY理論和一級相變理論在二維熔化中的適用性一直是大家關(guān)注的焦點。膠體體系和顆粒體系是很好的模型體系，能夠清晰地反映晶體熔化的微觀過程，為尋求熔化過程的本質(zhì)機理提供參考。利用計算機模擬研究二維熔化一直是一種非常靈活有效的方法，大量的系統(tǒng)化模擬研究為后續(xù)的實驗工作鋪平道路，也是構(gòu)建理論基礎(chǔ)的一個主要方法。本文總結(jié)了近些年關(guān)于二維熔化的實驗和模擬工作：不同實驗體系關(guān)于熔化理論的驗證，不同的作用勢對于熔化機制的影響，以及熔化過程中結(jié)構(gòu)和動力學(xué)的異質(zhì)性。豐富和完善熔化理論，研究相變行為的微觀過程，不僅對凝聚態(tài)物理有基礎(chǔ)意義，也對材料物理，納米技術(shù)有一定的指導(dǎo)意義。

參考文獻

[1]K. J. Strandburg. Two-dimensional melting. Rev. Mod. Phys.， 1988， 60（1）： 161-207.

[2]Y. Han， N. Y. Ha， A. M. Alsayed and A. G. Yodh. Melting of two-dimensional tunable-diameter colloidal crystals. Phys. Rev. E， 2008， 77， 041406.

[3]X. H. Zheng and R. Grieve. Melting behavior of single two-dimensional crystals. Phys. Rev. B， 2006， 73， 064205.

[4]J. Lee and K. J. Strandburg. First-order melting transition of the hard-disk system. Phys. Rev. B， 1992， 46（17）： 11190-11193.

[5]W. Janke and H. Kleinhert. From First-Order to Two Continuous Melting Transition： Monte Carlo Study of a New 2D Lattice-Defect Model. Phys. Rev. Lett.， 1988， 61（20）： 2344-2347.

[6]K. Chen， T. Kaplan and M. Mostoller. Melting in Two-Dimensional Lennard-Jones Systems： Observation of a Metastable Hexatic Phase. Phys. Rev. Lett.， 1995， 74（20）： 4019-4022.

[7]H. Shiba， A. Onuki and T. Araki. Structural and dynamical heterogeneities in two-dimensional melting. EPL， 2009， 86， 66004.

[8]H. Zhang， M. Khalkhali， Q. Liu and J. F. Douglas. String-like cooperative motion in homogeneous melting. J. Chem. Phys.， 2013， 138， 12A538.

作者簡介

席航波（1989-），男，漢族，江蘇宿遷人，蘇州大學(xué)軟凝聚態(tài)物理及交叉研究中心，碩士研究生，研究方向：膠體晶體的結(jié)構(gòu)和動力學(xué)。endprint