二維熔化的研究進(jìn)展

席航波+張海波

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

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

中圖分類號(hào)：TF841 文獻(xiàn)標(biāo)識(shí)碼：A 文章編號(hào)：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

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

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

1 實(shí)驗(yàn)研究

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

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

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

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

2 計(jì)算機(jī)模擬

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

2.1 蒙特卡羅

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

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

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

2.2 分子動(dòng)力學(xué)

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

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

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

3 結(jié)束語(yǔ)

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

參考文獻(xiàn)

[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.

作者簡(jiǎn)介

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