算實(shí)踐)
2.2° 的魔法當(dāng)兩張?jiān)印皾O網(wǎng)”疊在一起材料科學(xué)被改寫(xiě)了如果你對(duì)材料科學(xué)的印象還停留在“摻雜”、“高溫高壓”和“元素周期表排列組合”那么這兩年最值得關(guān)注的一個(gè)技術(shù)方向可能會(huì)徹底顛覆你的認(rèn)知僅僅把一張?jiān)雍穸鹊牟牧闲D(zhuǎn)一個(gè)特定的角度疊在另一張同樣材料的上面就能讓絕緣體變成超導(dǎo)體讓普通金屬變成磁性材料。這不是科幻設(shè)定而是魔角莫爾材料Magic Moiré Materials研究領(lǐng)域正在發(fā)生的真實(shí)突破。很多人第一次聽(tīng)到“魔角石墨烯”時(shí)第一反應(yīng)是“這又是材料學(xué)家在炒概念”。但實(shí)際上這個(gè)方向解決的是一個(gè)長(zhǎng)期困擾凝聚態(tài)物理和材料科學(xué)的根本問(wèn)題我們能不能不用改變材料的化學(xué)成分僅僅通過(guò)幾何堆疊和轉(zhuǎn)角主動(dòng)“設(shè)計(jì)”出全新的量子物態(tài)傳統(tǒng)方法調(diào)整材料性質(zhì)靠的是“往鍋里加料”。而魔角莫爾材料提供了一條完全不同的路徑——“我什么都不加只是把兩張網(wǎng)錯(cuò)開(kāi)一點(diǎn)角度疊起來(lái)”。這個(gè)思路的優(yōu)雅之處在于它把材料設(shè)計(jì)從“化學(xué)配方”時(shí)代推進(jìn)到了“幾何工程”時(shí)代。這篇文章我會(huì)從一個(gè)開(kāi)發(fā)者比較容易理解的“干涉條紋”現(xiàn)象講起拆解莫爾超晶格的形成機(jī)制、魔角為什么“魔”以及為什么這個(gè)方向會(huì)成為當(dāng)前凝聚態(tài)物理和材料計(jì)算領(lǐng)域最熱門的賽道之一。同時(shí)我會(huì)給出一個(gè)基于 Python 的莫爾勢(shì)場(chǎng)模擬示例、TB 模型構(gòu)建思路和能帶計(jì)算流程幫助剛接觸這個(gè)領(lǐng)域的讀者快速建立從物理概念到計(jì)算實(shí)踐的完整鏈路。1. 這篇文章真正要解決的問(wèn)題如果你不是凝聚態(tài)物理專業(yè)出身接觸“魔角莫爾材料”時(shí)通常會(huì)有三個(gè)困惑第一個(gè)困惑莫爾條紋我見(jiàn)過(guò)但為什么它是材料科學(xué)的大事件你小時(shí)候可能見(jiàn)過(guò)兩張窗紗或蚊帳疊在一起產(chǎn)生的那種明暗相間的花紋。兩根絲線本來(lái)間距很小但重疊后出現(xiàn)了一組間距大得多的條紋。莫爾超晶格本質(zhì)上是同一個(gè)物理現(xiàn)象兩個(gè)周期接近但略有偏差的原子晶格疊加后會(huì)形成一個(gè)全新的、周期大得多的公共周期結(jié)構(gòu)。這個(gè)看似簡(jiǎn)單的幾何效應(yīng)在電子世界里會(huì)引發(fā)劇烈的物理性質(zhì)變化。第二個(gè)困惑為什么“轉(zhuǎn)一個(gè)角度”就能改變材料的導(dǎo)電性關(guān)鍵在于“平帶”flat band。當(dāng)兩層石墨烯的轉(zhuǎn)角接近某個(gè)特定值理論預(yù)言約為 1.1°時(shí)電子能帶會(huì)被急劇壓縮形成幾乎不隨動(dòng)量變化的“平帶”。平帶意味著電子的動(dòng)能幾乎被“凍結(jié)”了此時(shí)電子與電子之間的相互作用庫(kù)侖排斥在總能量中的占比極大提升。換句話說(shuō)關(guān)聯(lián)效應(yīng)被放大了。第三個(gè)困惑這個(gè)方向跟我一個(gè)寫(xiě)代碼的有什么關(guān)系如果你做的是材料計(jì)算、量子化學(xué)模擬、機(jī)器學(xué)習(xí)勢(shì)函數(shù)或者高性能計(jì)算魔角莫爾材料是一個(gè)極佳的研究與應(yīng)用場(chǎng)景。莫爾超晶格的原子數(shù)量從幾千到幾萬(wàn)甚至幾十萬(wàn)起步遠(yuǎn)超傳統(tǒng) DFT密度泛函理論能處理的規(guī)模。這意味著必須依賴 TB緊束縛模型、深度學(xué)習(xí)勢(shì)函數(shù)、連續(xù)介質(zhì)模型等新的計(jì)算范式。換句話說(shuō)這個(gè)領(lǐng)域天然就是計(jì)算驅(qū)動(dòng)和數(shù)據(jù)驅(qū)動(dòng)的計(jì)算工具鏈的開(kāi)發(fā)者、算法工程師在這里有大量發(fā)揮空間。本文的目標(biāo)是讓有 Linux 環(huán)境下 Python 使用經(jīng)驗(yàn)和一定線性代數(shù)基礎(chǔ)的讀者在閱讀完這篇文章后能夠理解魔角莫爾材料背后的核心概念并跑通一個(gè)最小化的莫爾超晶格模擬與能帶計(jì)算流程。2. 基礎(chǔ)概念與核心原理從莫爾條紋到魔角平帶2.1 莫爾條紋與莫爾超晶格莫爾條紋的物理本質(zhì)是空間頻率的“差頻”現(xiàn)象。想象兩把間距分別為 (a_1) 和 (a_2) 的梳子疊在一起當(dāng)它們的間距很接近時(shí)重疊區(qū)域會(huì)出現(xiàn)一個(gè)遠(yuǎn)大于二者間距的“拍頻”周期[ L \approx \frac{a_1 \cdot a_2}{|a_1 - a_2|} ]在二維材料中當(dāng)我們把兩層六角晶格例如石墨烯以一個(gè)小角度 (\theta) 疊加時(shí)會(huì)形成一個(gè)莫爾周期[ L_M \frac{a}{2 \sin(\theta/2)} ]其中 (a) 是晶格常數(shù)。當(dāng) (\theta \approx 1.1^\circ) 時(shí)(L_M) 可以達(dá)到晶格常數(shù)的約 55 倍即約 13.4 納米。這個(gè)納米級(jí)的周期結(jié)構(gòu)對(duì)于電子來(lái)說(shuō)如同一個(gè)“人工原子晶格”其上每個(gè)周期格點(diǎn)都像一個(gè)巨大的“人造原子”。這個(gè)概念之所以重要是因?yàn)樗锩缘馗淖兞恕罢{(diào)控材料”的方式以前要實(shí)現(xiàn)對(duì)電子的約束必須依靠物理刻蝕或摻雜形成勢(shì)壘現(xiàn)在只需“疊放 旋轉(zhuǎn)”就能讓電子感受到周期性勢(shì)場(chǎng)。2.2 魔角與平帶在莫爾超晶格中電子的能帶結(jié)構(gòu)與轉(zhuǎn)角 (\theta) 有著強(qiáng)依賴關(guān)系。當(dāng)轉(zhuǎn)角較大時(shí)例如 (5^\circ)兩層材料之間的層間耦合較弱電子行為還接近單層石墨烯的線性色散。但隨著轉(zhuǎn)角逐漸減小到某個(gè)“魔幻”的角度層間耦合會(huì)把能帶強(qiáng)烈地重整化導(dǎo)致費(fèi)米面附近的電子速度趨近于零——這就是平帶。需要特別強(qiáng)調(diào)的是魔角石墨烯中的“魔角”并不止一個(gè)。理論計(jì)算表明在約 (1.1^\circ) 附近有一系列魔角候選但 (1.1^\circ) 左右的第一個(gè)魔角效應(yīng)最顯著。在魔角位置平帶帶寬可以小于 10 meV這意味著電子的動(dòng)能尺度meV量級(jí)遠(yuǎn)低于相互作用能量尺度10-30 meV。此時(shí)平帶中的電子不能再被當(dāng)作“自由電子氣”處理而必須用強(qiáng)關(guān)聯(lián)物理的語(yǔ)言來(lái)描述。2.3 莫爾材料 vs 傳統(tǒng)摻雜半導(dǎo)體理解魔角莫爾材料為什么“顛覆認(rèn)知”最好的方式是與傳統(tǒng)半導(dǎo)體調(diào)控手段做對(duì)比。維度傳統(tǒng)摻雜半導(dǎo)體魔角莫爾材料調(diào)控方式化學(xué)摻雜、外加電場(chǎng)、應(yīng)力工程幾何堆疊、轉(zhuǎn)角控制、垂直電場(chǎng)低溫量子物態(tài)需要極低溫和復(fù)雜外場(chǎng)在轉(zhuǎn)角精確時(shí)可以自發(fā)涌現(xiàn)能帶尺度eV 量級(jí)meV 量級(jí)關(guān)聯(lián)效應(yīng)弱關(guān)聯(lián)為主單電子近似成立強(qiáng)關(guān)聯(lián)主導(dǎo)單電子近似失效規(guī)?；苽涑墒旃饪獭⑼庋庸に嚩S材料堆疊/干法轉(zhuǎn)移尚處實(shí)驗(yàn)室階段計(jì)算難度常規(guī) DFT 可很好處理需要超越 DFT 的模型和算法這個(gè)對(duì)比可以幫助理解魔角莫爾材料為什么讓凝聚態(tài)物理學(xué)家興奮它把原本只能在極低溫極高壓等極端條件下觀察到的量子現(xiàn)象變成了“桌面實(shí)驗(yàn)”中可以主動(dòng)調(diào)節(jié)的研究平臺(tái)。2.4 Simon Becker 團(tuán)隊(duì)的工作定位從公開(kāi)資料看Simon Becker 團(tuán)隊(duì)的工作主要集中在“魔角莫爾材料”的精確制備、表征和物性調(diào)控方面。這類研究的核心價(jià)值在于把“魔角”概念從石墨烯拓展到更廣泛的二維材料體系同時(shí)通過(guò)掃描隧道顯微鏡/輸運(yùn)測(cè)量等手段證實(shí)平帶中的強(qiáng)關(guān)聯(lián)電子行為。這類工作其實(shí)是近年來(lái)凝聚態(tài)物理諾貝爾獎(jiǎng)級(jí)發(fā)現(xiàn)的延續(xù)——從實(shí)驗(yàn)上證實(shí)了“轉(zhuǎn)角自由度”可以作為一種強(qiáng)大的材料設(shè)計(jì)維度。對(duì)于開(kāi)發(fā)者而言可以這樣理解過(guò)去的材料工程師是“配方師”現(xiàn)在這個(gè)領(lǐng)域的研究者是“建筑師”他們不改變材料的原子種類只改變?cè)訉拥呐帕蟹绞健?. 環(huán)境準(zhǔn)備與前置條件在進(jìn)行莫爾超晶格模擬前需要準(zhǔn)備好計(jì)算環(huán)境。本文沒(méi)有依賴特定的商業(yè)軟件思路是用 Python 生態(tài)完成“莫爾勢(shì)場(chǎng)可視化 緊束縛模型構(gòu)建 能帶計(jì)算”的最小閉環(huán)。3.1 操作系統(tǒng)與硬件建議推薦使用 Ubuntu 20.04 LTS 或更新版本因?yàn)楹罄m(xù)如果需要擴(kuò)展到 GPU 加速的深度學(xué)習(xí)勢(shì)函數(shù)計(jì)算Linux 環(huán)境兼容性最好。Windows 用戶可以通過(guò) WSL2 搭建同樣的環(huán)境。硬件方面僅跑本文的示例代碼4 核 CPU 8GB 內(nèi)存就足夠。但如果要構(gòu)建超過(guò) 1 萬(wàn)個(gè)原子的莫爾超晶格并進(jìn)行 TB 模型對(duì)角化建議至少 16GB 內(nèi)存并考慮使用稀疏矩陣求解器。如果跑機(jī)器學(xué)習(xí)勢(shì)函數(shù)分子動(dòng)力學(xué)則需要 NVIDIA GPU建議顯存 16GB 以上。3.2 Python 環(huán)境與依賴包推薦使用 Anaconda 或 Miniconda 管理環(huán)境。創(chuàng)建一個(gè)專用的虛擬環(huán)境conda create -n moire python3.10 -y conda activate moire核心依賴包如下pip install numpy scipy matplotlib pip install ase # Atomic Simulation Environment用于構(gòu)建原子結(jié)構(gòu)版本要求以實(shí)際安裝為準(zhǔn)本文代碼基于 Python 3.10 及以上版本編寫(xiě)NumPy 2.x 或 1.26 均可。3.3 輔助工具可視化結(jié)構(gòu)可以配合 VESTA 或 ASE GUI 查看生成的原子結(jié)構(gòu)。如果需要做第一性原理驗(yàn)證可以從 PWscfQuantum ESPRESSO或者 VASP 開(kāi)始但本文不涉及 DFT 實(shí)操。如果需要處理 TB 模型推薦 PythTB 包pip install pythtb或自行用 NumPy 實(shí)現(xiàn)。4. 核心流程拆解從物理問(wèn)題到計(jì)算實(shí)現(xiàn)一個(gè)完整的莫爾材料模擬流程包含以下幾個(gè)步驟。4.1 構(gòu)建雙層的莫爾超晶格模型這是最關(guān)鍵的步驟。構(gòu)建莫爾超晶格有兩種常見(jiàn)方案方案一直接構(gòu)建大超胞。根據(jù)層間轉(zhuǎn)角確定超胞大小一個(gè)能顯示清晰莫爾周期的超胞往往包含數(shù)千到數(shù)萬(wàn)原子適合后續(xù)用機(jī)器學(xué)習(xí)勢(shì)函數(shù)做結(jié)構(gòu)弛豫和分子動(dòng)力學(xué)。方案二使用連續(xù)介質(zhì)模型Continuum Model。不顯式構(gòu)建原子結(jié)構(gòu)而是用一個(gè)周期性勢(shì)場(chǎng)描述兩層之間的相互作用。這種方法計(jì)算量極小適合快速解析平帶和電子結(jié)構(gòu)。后面的示例代碼采用這個(gè)思路。4.2 確定轉(zhuǎn)角與失配轉(zhuǎn)角決定了莫爾波長(zhǎng)。對(duì)于石墨烯晶格常數(shù) 2.46 ?(1.1^\circ) 對(duì)應(yīng)的莫爾周期約為 13.4 nm。為了在有限尺寸的超胞內(nèi)精確實(shí)現(xiàn)目標(biāo)轉(zhuǎn)角通常使用公度角commensurate angle的概念即找到一個(gè)整數(shù) ((m, n)) 使得超胞能夠在兩層晶格上同時(shí)閉合。這套方法在文獻(xiàn)中已經(jīng)很成熟感興趣的讀者可以搜索“石墨烯莫爾超晶格公度角構(gòu)建”。4.3 構(gòu)建緊束縛TB模型TB 模型是連接原子尺度與電子結(jié)構(gòu)的橋梁。在莫爾體系中可以只考慮石墨烯的 (p_z) 軌道并引入層間跳躍interlayer hopping隨距離變化的近似層間跳躍在莫爾周期上會(huì)呈現(xiàn)出局域增強(qiáng)或減弱的圖案這恰好對(duì)應(yīng)了“AA 堆垛區(qū)域”和“AB/BA 堆垛區(qū)域”的周期性分布。4.4 計(jì)算能帶與態(tài)密度在 TB 模型下能帶計(jì)算轉(zhuǎn)化為對(duì)一個(gè)大稀疏矩陣的哈密頓量 (H(k)) 進(jìn)行對(duì)角化。因?yàn)槟獱柍艽髮?shí)空間哈密頓量維數(shù)很容易超過(guò)數(shù)萬(wàn)通常選擇使用稀疏矩陣特征值求解器scipy.sparse.linalg計(jì)算最低若干條能帶或者沿高對(duì)稱路徑掃描。4.5 驗(yàn)證與可視化能帶圖是否展現(xiàn)出“平帶”是判斷模型是否捕捉到魔角效應(yīng)的核心指標(biāo)。可視化部分要輸出能帶色散圖和某一能量窗口內(nèi)的局域態(tài)密度圖。5. 完整示例與代碼實(shí)現(xiàn)下面用三個(gè)代碼示例完整演示“莫爾超晶格勢(shì)場(chǎng)可視化 — TB 模型能帶計(jì)算 — 平帶態(tài)密度驗(yàn)證”的流程。5.1 示例一莫爾超晶格勢(shì)場(chǎng)可視化這個(gè)示例用 NumPy 生成一個(gè)二維勢(shì)場(chǎng)模擬兩層六角晶格在不同轉(zhuǎn)角下的莫爾干涉圖案。# 文件路徑moire_potential.py import numpy as np import matplotlib.pyplot as plt def hex_lattice_potential(x, y, wavelength1.0): 生成一個(gè)簡(jiǎn)化的六角晶格周期勢(shì)場(chǎng)。 六角晶格可以用兩組成 60 度夾角的格矢描述 我們疊加三個(gè)方向上的余弦波來(lái)近似。 k1 np.array([1.0, 0.0]) k2 np.array([0.5, np.sqrt(3) / 2.0]) k3 k1 - k2 phase1 2 * np.pi * (k1[0] * x k1[1] * y) / wavelength phase2 2 * np.pi * (k2[0] * x k2[1] * y) / wavelength phase3 2 * np.pi * (k3[0] * x k3[1] * y) / wavelength return np.cos(phase1) np.cos(phase2) np.cos(phase3) def rotate_coords(x, y, theta_deg): 將坐標(biāo)旋轉(zhuǎn) theta_deg 度角度制。 theta np.deg2rad(theta_deg) x_rot x * np.cos(theta) - y * np.sin(theta) y_rot x * np.sin(theta) y * np.cos(theta) return x_rot, y_rot # 構(gòu)建網(wǎng)格 nx, ny 1200, 1200 x np.linspace(-15, 15, nx) y np.linspace(-15, 15, ny) X, Y np.meshgrid(x, y) # 第一層晶格 V1 hex_lattice_potential(X, Y, wavelength1.0) # 第二層晶格旋轉(zhuǎn) theta_deg theta_deg 3.0 # 使用 3 度轉(zhuǎn)角莫爾條紋清晰可見(jiàn) X_rot, Y_rot rotate_coords(X, Y, theta_deg) V2 hex_lattice_potential(X_rot, Y_rot, wavelength1.0) # 莫爾總勢(shì)場(chǎng)干涉圖案 V_moire V1 V2 plt.figure(figsize(8, 8)) plt.pcolormesh(X, Y, V_moire, cmapviridis, shadingauto) plt.title(fMoiré Pattern at θ {theta_deg}°) plt.xlabel(x) plt.ylabel(y) plt.colorbar(labelPotential) plt.savefig(moire_pattern.png, dpi150) plt.show()運(yùn)行方式python moire_potential.py如果一切正常圖像中會(huì)看到明顯的菱形或多邊形莫爾超晶格條紋。調(diào)整theta_deg為更小的角度如 1.1°莫爾周期會(huì)變得更大。5.2 示例二基于 PythTB 的魔角石墨烯能帶計(jì)算PythTB 是一個(gè)輕量的緊束縛模型工具包適合新手快速驗(yàn)證物理模型。下面的示例構(gòu)建一個(gè)雙層石墨烯的簡(jiǎn)化莫爾模型只考慮一個(gè)魔角附近的格子。# 文件路徑tb_band_structure.py # 需要先安裝 pythtbpip install pythtb import pythtb as ptb import numpy as np import matplotlib.pyplot as plt # 定義晶格常數(shù) a 2.46 # 石墨烯原始晶格常數(shù)單位埃 # 構(gòu)建雙層石墨烯的簡(jiǎn)化模型每個(gè)原胞包含4個(gè)原子 # 這里僅演示沿高對(duì)稱路徑計(jì)算能帶的方法 lat [[a, 0.0], [0.0, a]] # 4個(gè)基矢坐標(biāo)下層A、下層B、上層A、上層B orb [[0.0, 0.0], [a/2, 0.0], [0.0, a/2], [a/2, a/2]] model ptb.tb_model(2, 2, lat, orb) # 最近鄰跳躍 model.set_hop(0, 0, 1, 2.8) # 層內(nèi)最近鄰跳躍 model.set_hop(1, 1, 0, 2.8) model.set_hop(0, 0, 3, 2.8) model.set_hop(1, 2, 3, 2.8) # 層間跳躍簡(jiǎn)化不考慮轉(zhuǎn)角的空間調(diào)制 model.set_hop(0, 0, 2, 0.3) # 主要層間跳躍 model.set_hop(0, 1, 3, 0.3) # 沿 K → Gamma → M 路徑計(jì)算 path [[0.0, 0.0], [1/3, 1/3], [0.0, 0.0], [0.5, 0.0]] (kpts, labels, paths) model.k_path(path, 201) evals model.solve_all(kpts) # 作圖 fig, ax plt.subplots(figsize(6, 5)) for band in evals: ax.plot(paths, band, colorblack, linewidth1.2) ax.set_ylabel(Energy (eV)) ax.set_xticks(labels) ax.set_xticklabels([r$\Gamma$, r$K$, r$\Gamma$, r$M$]) ax.set_title(TB band structure (schematic)) plt.savefig(tb_bands.png, dpi150) plt.show()注意這是一個(gè)簡(jiǎn)化的雙層石墨烯模型并沒(méi)有真正實(shí)現(xiàn) 1.1° 轉(zhuǎn)角下的莫爾超胞。如果要真實(shí)刻畫(huà)魔角需要建立數(shù)千原子的超胞并將層間跳躍設(shè)置為隨莫爾周期呈空間調(diào)制形式。這段代碼的作用是幫助讀者理解 PythTB 的基本使用流程。5.3 示例三魔角平帶態(tài)密度驗(yàn)證當(dāng)模型正確包含莫爾周期和層間耦合調(diào)制時(shí)能帶會(huì)在費(fèi)米能級(jí)附近出現(xiàn)電子速度接近零的平帶。此時(shí)態(tài)密度會(huì)出現(xiàn)一個(gè)尖銳的峰。下面的代碼演示如何基于 TB 模型輸出態(tài)密度。# 文件路徑dos_check.py import numpy as np import matplotlib.pyplot as plt def gaussian_broadening(energies, smearing0.005): 將一組能級(jí)用高斯展寬得到態(tài)密度。 energies: 所有 k 點(diǎn)所有能帶的能量數(shù)組 smearing: 高斯展寬寬度單位 eV e_min np.min(energies) - 5 * smearing e_max np.max(energies) 5 * smearing grid np.linspace(e_min, e_max, 800) dos np.zeros_like(grid) for e in energies.flatten(): dos np.exp(-0.5 * ((grid - e) / smearing) ** 2) dos * 1.0 / (np.sqrt(2 * np.pi) * smearing * len(energies.flatten())) return grid, dos # 模擬一個(gè)包含平帶的能量集合 # 假設(shè)在 -0.05 eV 附近存在大量平帶電子態(tài) np.random.seed(42) flat_band_energies np.random.normal(-0.05, 0.003, size3000) scattered_energies np.random.normal(0.15, 0.1, size5000) all_energies np.concatenate([flat_band_energies, scattered_energies]).reshape(-1, 1) grid, dos gaussian_broadening(all_energies, smearing0.005) plt.figure(figsize(7, 5)) plt.plot(grid, dos, colordarkred, linewidth2) plt.axvline(0, colorgray, linestyle--, linewidth1) plt.xlabel(Energy (eV)) plt.ylabel(DOS (arb. units)) plt.title(Density of States with Van Hove Singularity / Flat Band) plt.savefig(dos_flatband.png, dpi150) plt.show()輸出圖中可以清楚看到在費(fèi)米能級(jí)附近有一個(gè)非常尖銳的態(tài)密度峰這便是平帶的標(biāo)志。如果直接對(duì)真實(shí)魔角石墨烯做 TB 計(jì)算會(huì)在這個(gè)位置得到類似結(jié)果。6. 運(yùn)行結(jié)果與效果驗(yàn)證6.1 莫爾圖案結(jié)果判斷運(yùn)行示例一后保存的moire_pattern.png應(yīng)該顯示一個(gè)具有明顯周期性的明暗圖案若圖案只是模糊無(wú)序的噪聲可能是網(wǎng)格分辨率不足或波長(zhǎng)與坐標(biāo)范圍不匹配。若圖案呈現(xiàn)出清晰的長(zhǎng)周期菱形條紋并且條紋間隔明顯大于單層晶格周期則說(shuō)明莫爾干涉被正確模擬。6.2 能帶圖結(jié)果判斷示例二的能帶圖只是一個(gè)簡(jiǎn)化的雙層石墨烯示意。真正驗(yàn)證魔角平帶時(shí)應(yīng)該觀察在費(fèi)米能級(jí)附近是否存在帶寬極小低于 10 meV的平帶。平帶出現(xiàn)的角度是否接近理論預(yù)言的魔角。隨著轉(zhuǎn)角偏離魔角平帶是否迅速消失并變成色散帶。一個(gè)有效的判斷方式是掃描多個(gè)轉(zhuǎn)角參數(shù)分別計(jì)算費(fèi)米面附近的帶寬然后繪制“帶寬-轉(zhuǎn)角”曲線。帶寬最小處即為該模型的魔角位置。這在計(jì)算上確實(shí)需要不少資源但對(duì)于理論學(xué)習(xí)來(lái)說(shuō)非常直觀。6.3 態(tài)密度驗(yàn)證示例三運(yùn)行成功后會(huì)看到一個(gè)尖銳的態(tài)密度峰。如果使用真實(shí) TB 模型計(jì)算還能進(jìn)一步通過(guò)投影態(tài)密度分辨各個(gè)堆垛區(qū)域AA、AB、BA對(duì)平帶的貢獻(xiàn)。通常AA 堆垛區(qū)域是平帶態(tài)密度的主要貢獻(xiàn)者這也與 STM 實(shí)驗(yàn)觀察到的莫爾超晶格亮斑位置一致。6.4 如果運(yùn)行失敗優(yōu)先檢查以下幾點(diǎn)PythTB 是否安裝成功可以執(zhí)行python -c import pythtb; print(pythtb.__version__)。NumPy 和 Matplotlib 版本是否兼容。運(yùn)行環(huán)境是否處于conda activate moire中。對(duì)于大超胞 TB 模型內(nèi)存不足時(shí)會(huì)出現(xiàn)MemoryError需要改用稀疏存儲(chǔ)或分塊計(jì)算。7. 常見(jiàn)問(wèn)題與排查思路問(wèn)題現(xiàn)象可能原因排查方式解決方案莫爾圖案周期過(guò)小看不出明顯的超晶格條紋轉(zhuǎn)角設(shè)置過(guò)大或坐標(biāo)范圍太小減小轉(zhuǎn)角例如從 3° 改為 1.1°擴(kuò)大坐標(biāo)范圍調(diào)整theta_deg參數(shù)增加np.linspace的上下限PythTB 能帶圖只有四條很窄的帶沒(méi)有平帶特征模型只是普通雙層石墨烯未包含莫爾超胞和層間調(diào)制檢查模型原子數(shù)與轉(zhuǎn)角確認(rèn)層間跳躍是否為固定值建立大超胞將層間跳躍設(shè)置為隨莫爾周期變化的空間函數(shù)態(tài)密度圖中費(fèi)米能級(jí)附近沒(méi)有尖峰平帶不在費(fèi)米能級(jí)附近或能帶過(guò)于色散檢查費(fèi)米能級(jí)位置調(diào)整填充或電子數(shù)調(diào)整平帶能量范圍或增加層間耦合強(qiáng)度構(gòu)建具有 1.1° 轉(zhuǎn)角的莫爾超胞時(shí)原子數(shù)爆炸直接使用小角度會(huì)產(chǎn)生巨大的超胞用連續(xù)介質(zhì)模型進(jìn)行快速預(yù)研或使用公度角近似在 TB 之前先用 Continuum Model 驗(yàn)證或選擇更大的公度角大超胞 TB 模型對(duì)角化內(nèi)存不足哈密頓量矩陣維度過(guò)大且使用了稠密矩陣存儲(chǔ)檢查是否有np.asarray或np.matrix轉(zhuǎn)換改用scipy.sparse.linalg.eigsh求解最低特征值或使用pythtb的并行模式弛豫后莫爾結(jié)構(gòu)變形導(dǎo)致平帶消失忽略了原子面外弛豫重建效應(yīng)對(duì)比弛豫前后的結(jié)構(gòu)差異檢查原子位移是否超過(guò) 10% 晶格常數(shù)在模型中引入面外位移項(xiàng)使用機(jī)器學(xué)習(xí)勢(shì)函數(shù)做更大規(guī)模結(jié)構(gòu)優(yōu)化與實(shí)驗(yàn)報(bào)道的魔角數(shù)值對(duì)不上模型缺少電子關(guān)聯(lián)、自旋軌道耦合、應(yīng)變等修正對(duì)比文獻(xiàn)中的 TB 參數(shù)采用文獻(xiàn)已驗(yàn)證的 TB 參數(shù)集或把魔角計(jì)算改成 LDA/DFT 的自洽驗(yàn)證8. 最佳實(shí)踐與工程建議8.1 從“離散模型”到“連續(xù)模型”的分層策略魔角莫爾材料最大的計(jì)算瓶頸是真實(shí)現(xiàn)象發(fā)生在納米尺度而真實(shí)原子間距在埃尺度這就要求超胞動(dòng)輒上萬(wàn)個(gè)原子。工程上不建議一上來(lái)就做全原子 DFT而是采用分層策略先用連續(xù)介質(zhì)模型Continuum Model解析地考察能帶和轉(zhuǎn)角依賴關(guān)系。用 TB 模型或機(jī)器學(xué)習(xí)勢(shì)函數(shù)對(duì)大超胞做結(jié)構(gòu)弛豫。最后用 DFT 或 GWBSE 對(duì)關(guān)鍵構(gòu)型做精確驗(yàn)證。這套流程在公開(kāi)文獻(xiàn)中已成為主流方案能大幅減少無(wú)效計(jì)算。8.2 參數(shù)文件與實(shí)驗(yàn)記錄管理莫爾材料計(jì)算涉及大量可調(diào)參數(shù)轉(zhuǎn)角、晶格失配、層間距、應(yīng)變、跳躍參數(shù)、外場(chǎng)等。強(qiáng)烈建議為每一個(gè)計(jì)算任務(wù)建立獨(dú)立的參數(shù) JSON 文件{ system: tBG, theta_deg: 1.1, lattice_constant_a: 2.46, hopping_intra: 2.8, hopping_inter: 0.3, model: TB, supercell_size: 55, comment: initial guess for magic angle }這樣即使幾個(gè)月后回看也能準(zhǔn)確復(fù)現(xiàn)每次計(jì)算。8.3 版本兼容與可重復(fù)性Python 的 NumPy/SciPy 在不同版本間行為有差異建議使用conda env export environment.yaml鎖定環(huán)境。對(duì)需要發(fā)布的結(jié)果可以同時(shí)記錄計(jì)算軟件的 commit hash隨機(jī)種子參數(shù)文件的 md5 值8.4 異常處理與日志記錄在批量掃描轉(zhuǎn)角參數(shù)時(shí)某一組參數(shù)可能導(dǎo)致特征值求解不收斂。建議用try/except捕獲異常并記錄日志而不是讓整個(gè)任務(wù)崩潰import logging logging.basicConfig(filenamemoire_scan.log, levellogging.INFO) for theta in theta_list: try: bands compute_bands(theta) logging.info(ftheta{theta:.3f} OK, bandwidth{np.max(bands)-np.min(bands):.6f}) except Exception as e: logging.error(ftheta{theta:.3f} failed: {e})8.5 安全邊界與生產(chǎn)環(huán)境提醒雖然這是研究型計(jì)算但涉及大規(guī)模集群使用時(shí)仍需注意不要在生產(chǎn)集群上的共享目錄直接跑大型 job容易干擾他人并造成數(shù)據(jù)混寫(xiě)。涉及商業(yè)軟件如 VASP時(shí)遵守所在機(jī)構(gòu)的許可證授權(quán)范圍。對(duì)于 TB 或 MD 計(jì)算如果會(huì)運(yùn)行幾天務(wù)必定期輸出 checkpoint以便中途失敗后從最近狀態(tài)恢復(fù)。9. 總結(jié)與后續(xù)學(xué)習(xí)方向魔角莫爾材料真正顛覆性的地方并不在于“轉(zhuǎn)角”這個(gè)動(dòng)作本身而在于它改變了材料設(shè)計(jì)的思維模式物態(tài)不再只是元素周期表中的一行一列而可以用幾何堆疊在參數(shù)空間中連續(xù)“調(diào)”出來(lái)。從計(jì)算角度這個(gè)領(lǐng)域的門檻并不在于理論物理公式有多深而在于如何高效構(gòu)建和處理超大規(guī)模的結(jié)構(gòu)模型。無(wú)論是連續(xù)介質(zhì)模型、TB 模型還是機(jī)器學(xué)習(xí)勢(shì)函數(shù)都是這個(gè)方向繞不開(kāi)的工具。對(duì)于剛開(kāi)始接觸這個(gè)方向的讀者下一步建議按順序做三件事運(yùn)行本文的莫爾勢(shì)場(chǎng)可視化代碼親眼看一次莫爾條紋隨轉(zhuǎn)角變化的演化。將 TB 模型擴(kuò)展到包含莫爾超胞至少理解“層間跳躍的空間調(diào)制”這一核心物理。閱讀魔角石墨烯的原始實(shí)驗(yàn)論文如 Nature 2018 年發(fā)表的超導(dǎo)發(fā)現(xiàn)工作將你在 TB 能帶圖中的結(jié)果與實(shí)驗(yàn)相圖對(duì)照。如果你希望進(jìn)一步深入可以關(guān)注三個(gè)延伸方向連續(xù)介質(zhì)模型與 Bistritzer-MacDonald 方程這是解析理解魔角平帶的標(biāo)準(zhǔn)武器。機(jī)器學(xué)習(xí)勢(shì)函數(shù)在莫爾材料弛豫中的應(yīng)用因?yàn)榇蟪Y(jié)構(gòu)優(yōu)化遠(yuǎn)超傳統(tǒng) DFT 的能力邊界。轉(zhuǎn)角光電子譜與 STM 實(shí)驗(yàn)數(shù)據(jù)的關(guān)聯(lián)分析這是連接計(jì)算與實(shí)驗(yàn)的關(guān)鍵橋梁。魔角莫爾材料這個(gè)方向正在快速?gòu)摹笆佟弊呦颉巴ㄓ枚S材料轉(zhuǎn)角平臺(tái)”。未來(lái)或許會(huì)出現(xiàn)更多基于轉(zhuǎn)角設(shè)計(jì)的拓?fù)涑瑢?dǎo)、量子磁性甚至量子計(jì)算平臺(tái)。對(duì)于每一位材料模擬和計(jì)算物理方向的開(kāi)發(fā)者來(lái)說(shuō)理解并掌握莫爾超晶格的計(jì)算方法正在變成一項(xiàng)越來(lái)越有競(jìng)爭(zhēng)力的基礎(chǔ)能力。建議先收藏本文等你的環(huán)境準(zhǔn)備好了把示例代碼跑一遍再回到上面的要點(diǎn)復(fù)盤一次理解的深度會(huì)完全不同。 魔角莫爾材料：從莫爾條紋到平帶工程的計(jì)算實(shí)踐
2.2° 的魔法當(dāng)兩張?jiān)印皾O網(wǎng)”疊在一起材料科學(xué)被改寫(xiě)了如果你對(duì)材料科學(xué)的印象還停留在“摻雜”、“高溫高壓”和“元素周期表排列組合”那么這兩年最值得關(guān)注的一個(gè)技術(shù)方向可能會(huì)徹底顛覆你的認(rèn)知僅僅把一張?jiān)雍穸鹊牟牧闲D(zhuǎn)一個(gè)特定的角度疊在另一張同樣材料的上面就能讓絕緣體變成超導(dǎo)體讓普通金屬變成磁性材料。這不是科幻設(shè)定而是魔角莫爾材料Magic Moiré Materials研究領(lǐng)域正在發(fā)生的真實(shí)突破。很多人第一次聽(tīng)到“魔角石墨烯”時(shí)第一反應(yīng)是“這又是材料學(xué)家在炒概念”。但實(shí)際上這個(gè)方向解決的是一個(gè)長(zhǎng)期困擾凝聚態(tài)物理和材料科學(xué)的根本問(wèn)題我們能不能不用改變材料的化學(xué)成分僅僅通過(guò)幾何堆疊和轉(zhuǎn)角主動(dòng)“設(shè)計(jì)”出全新的量子物態(tài)傳統(tǒng)方法調(diào)整材料性質(zhì)靠的是“往鍋里加料”。而魔角莫爾材料提供了一條完全不同的路徑——“我什么都不加只是把兩張網(wǎng)錯(cuò)開(kāi)一點(diǎn)角度疊起來(lái)”。這個(gè)思路的優(yōu)雅之處在于它把材料設(shè)計(jì)從“化學(xué)配方”時(shí)代推進(jìn)到了“幾何工程”時(shí)代。這篇文章我會(huì)從一個(gè)開(kāi)發(fā)者比較容易理解的“干涉條紋”現(xiàn)象講起拆解莫爾超晶格的形成機(jī)制、魔角為什么“魔”以及為什么這個(gè)方向會(huì)成為當(dāng)前凝聚態(tài)物理和材料計(jì)算領(lǐng)域最熱門的賽道之一。同時(shí)我會(huì)給出一個(gè)基于 Python 的莫爾勢(shì)場(chǎng)模擬示例、TB 模型構(gòu)建思路和能帶計(jì)算流程幫助剛接觸這個(gè)領(lǐng)域的讀者快速建立從物理概念到計(jì)算實(shí)踐的完整鏈路。1. 這篇文章真正要解決的問(wèn)題如果你不是凝聚態(tài)物理專業(yè)出身接觸“魔角莫爾材料”時(shí)通常會(huì)有三個(gè)困惑第一個(gè)困惑莫爾條紋我見(jiàn)過(guò)但為什么它是材料科學(xué)的大事件你小時(shí)候可能見(jiàn)過(guò)兩張窗紗或蚊帳疊在一起產(chǎn)生的那種明暗相間的花紋。兩根絲線本來(lái)間距很小但重疊后出現(xiàn)了一組間距大得多的條紋。莫爾超晶格本質(zhì)上是同一個(gè)物理現(xiàn)象兩個(gè)周期接近但略有偏差的原子晶格疊加后會(huì)形成一個(gè)全新的、周期大得多的公共周期結(jié)構(gòu)。這個(gè)看似簡(jiǎn)單的幾何效應(yīng)在電子世界里會(huì)引發(fā)劇烈的物理性質(zhì)變化。第二個(gè)困惑為什么“轉(zhuǎn)一個(gè)角度”就能改變材料的導(dǎo)電性關(guān)鍵在于“平帶”flat band。當(dāng)兩層石墨烯的轉(zhuǎn)角接近某個(gè)特定值理論預(yù)言約為 1.1°時(shí)電子能帶會(huì)被急劇壓縮形成幾乎不隨動(dòng)量變化的“平帶”。平帶意味著電子的動(dòng)能幾乎被“凍結(jié)”了此時(shí)電子與電子之間的相互作用庫(kù)侖排斥在總能量中的占比極大提升。換句話說(shuō)關(guān)聯(lián)效應(yīng)被放大了。第三個(gè)困惑這個(gè)方向跟我一個(gè)寫(xiě)代碼的有什么關(guān)系如果你做的是材料計(jì)算、量子化學(xué)模擬、機(jī)器學(xué)習(xí)勢(shì)函數(shù)或者高性能計(jì)算魔角莫爾材料是一個(gè)極佳的研究與應(yīng)用場(chǎng)景。莫爾超晶格的原子數(shù)量從幾千到幾萬(wàn)甚至幾十萬(wàn)起步遠(yuǎn)超傳統(tǒng) DFT密度泛函理論能處理的規(guī)模。這意味著必須依賴 TB緊束縛模型、深度學(xué)習(xí)勢(shì)函數(shù)、連續(xù)介質(zhì)模型等新的計(jì)算范式。換句話說(shuō)這個(gè)領(lǐng)域天然就是計(jì)算驅(qū)動(dòng)和數(shù)據(jù)驅(qū)動(dòng)的計(jì)算工具鏈的開(kāi)發(fā)者、算法工程師在這里有大量發(fā)揮空間。本文的目標(biāo)是讓有 Linux 環(huán)境下 Python 使用經(jīng)驗(yàn)和一定線性代數(shù)基礎(chǔ)的讀者在閱讀完這篇文章后能夠理解魔角莫爾材料背后的核心概念并跑通一個(gè)最小化的莫爾超晶格模擬與能帶計(jì)算流程。2. 基礎(chǔ)概念與核心原理從莫爾條紋到魔角平帶2.1 莫爾條紋與莫爾超晶格莫爾條紋的物理本質(zhì)是空間頻率的“差頻”現(xiàn)象。想象兩把間距分別為 (a_1) 和 (a_2) 的梳子疊在一起當(dāng)它們的間距很接近時(shí)重疊區(qū)域會(huì)出現(xiàn)一個(gè)遠(yuǎn)大于二者間距的“拍頻”周期[ L \approx \frac{a_1 \cdot a_2}{|a_1 - a_2|} ]在二維材料中當(dāng)我們把兩層六角晶格例如石墨烯以一個(gè)小角度 (\theta) 疊加時(shí)會(huì)形成一個(gè)莫爾周期[ L_M \frac{a}{2 \sin(\theta/2)} ]其中 (a) 是晶格常數(shù)。當(dāng) (\theta \approx 1.1^\circ) 時(shí)(L_M) 可以達(dá)到晶格常數(shù)的約 55 倍即約 13.4 納米。這個(gè)納米級(jí)的周期結(jié)構(gòu)對(duì)于電子來(lái)說(shuō)如同一個(gè)“人工原子晶格”其上每個(gè)周期格點(diǎn)都像一個(gè)巨大的“人造原子”。這個(gè)概念之所以重要是因?yàn)樗锩缘馗淖兞恕罢{(diào)控材料”的方式以前要實(shí)現(xiàn)對(duì)電子的約束必須依靠物理刻蝕或摻雜形成勢(shì)壘現(xiàn)在只需“疊放 旋轉(zhuǎn)”就能讓電子感受到周期性勢(shì)場(chǎng)。2.2 魔角與平帶在莫爾超晶格中電子的能帶結(jié)構(gòu)與轉(zhuǎn)角 (\theta) 有著強(qiáng)依賴關(guān)系。當(dāng)轉(zhuǎn)角較大時(shí)例如 (5^\circ)兩層材料之間的層間耦合較弱電子行為還接近單層石墨烯的線性色散。但隨著轉(zhuǎn)角逐漸減小到某個(gè)“魔幻”的角度層間耦合會(huì)把能帶強(qiáng)烈地重整化導(dǎo)致費(fèi)米面附近的電子速度趨近于零——這就是平帶。需要特別強(qiáng)調(diào)的是魔角石墨烯中的“魔角”并不止一個(gè)。理論計(jì)算表明在約 (1.1^\circ) 附近有一系列魔角候選但 (1.1^\circ) 左右的第一個(gè)魔角效應(yīng)最顯著。在魔角位置平帶帶寬可以小于 10 meV這意味著電子的動(dòng)能尺度meV量級(jí)遠(yuǎn)低于相互作用能量尺度10-30 meV。此時(shí)平帶中的電子不能再被當(dāng)作“自由電子氣”處理而必須用強(qiáng)關(guān)聯(lián)物理的語(yǔ)言來(lái)描述。2.3 莫爾材料 vs 傳統(tǒng)摻雜半導(dǎo)體理解魔角莫爾材料為什么“顛覆認(rèn)知”最好的方式是與傳統(tǒng)半導(dǎo)體調(diào)控手段做對(duì)比。維度傳統(tǒng)摻雜半導(dǎo)體魔角莫爾材料調(diào)控方式化學(xué)摻雜、外加電場(chǎng)、應(yīng)力工程幾何堆疊、轉(zhuǎn)角控制、垂直電場(chǎng)低溫量子物態(tài)需要極低溫和復(fù)雜外場(chǎng)在轉(zhuǎn)角精確時(shí)可以自發(fā)涌現(xiàn)能帶尺度eV 量級(jí)meV 量級(jí)關(guān)聯(lián)效應(yīng)弱關(guān)聯(lián)為主單電子近似成立強(qiáng)關(guān)聯(lián)主導(dǎo)單電子近似失效規(guī)?；苽涑墒旃饪獭⑼庋庸に嚩S材料堆疊/干法轉(zhuǎn)移尚處實(shí)驗(yàn)室階段計(jì)算難度常規(guī) DFT 可很好處理需要超越 DFT 的模型和算法這個(gè)對(duì)比可以幫助理解魔角莫爾材料為什么讓凝聚態(tài)物理學(xué)家興奮它把原本只能在極低溫極高壓等極端條件下觀察到的量子現(xiàn)象變成了“桌面實(shí)驗(yàn)”中可以主動(dòng)調(diào)節(jié)的研究平臺(tái)。2.4 Simon Becker 團(tuán)隊(duì)的工作定位從公開(kāi)資料看Simon Becker 團(tuán)隊(duì)的工作主要集中在“魔角莫爾材料”的精確制備、表征和物性調(diào)控方面。這類研究的核心價(jià)值在于把“魔角”概念從石墨烯拓展到更廣泛的二維材料體系同時(shí)通過(guò)掃描隧道顯微鏡/輸運(yùn)測(cè)量等手段證實(shí)平帶中的強(qiáng)關(guān)聯(lián)電子行為。這類工作其實(shí)是近年來(lái)凝聚態(tài)物理諾貝爾獎(jiǎng)級(jí)發(fā)現(xiàn)的延續(xù)——從實(shí)驗(yàn)上證實(shí)了“轉(zhuǎn)角自由度”可以作為一種強(qiáng)大的材料設(shè)計(jì)維度。對(duì)于開(kāi)發(fā)者而言可以這樣理解過(guò)去的材料工程師是“配方師”現(xiàn)在這個(gè)領(lǐng)域的研究者是“建筑師”他們不改變材料的原子種類只改變?cè)訉拥呐帕蟹绞健?. 環(huán)境準(zhǔn)備與前置條件在進(jìn)行莫爾超晶格模擬前需要準(zhǔn)備好計(jì)算環(huán)境。本文沒(méi)有依賴特定的商業(yè)軟件思路是用 Python 生態(tài)完成“莫爾勢(shì)場(chǎng)可視化 緊束縛模型構(gòu)建 能帶計(jì)算”的最小閉環(huán)。3.1 操作系統(tǒng)與硬件建議推薦使用 Ubuntu 20.04 LTS 或更新版本因?yàn)楹罄m(xù)如果需要擴(kuò)展到 GPU 加速的深度學(xué)習(xí)勢(shì)函數(shù)計(jì)算Linux 環(huán)境兼容性最好。Windows 用戶可以通過(guò) WSL2 搭建同樣的環(huán)境。硬件方面僅跑本文的示例代碼4 核 CPU 8GB 內(nèi)存就足夠。但如果要構(gòu)建超過(guò) 1 萬(wàn)個(gè)原子的莫爾超晶格并進(jìn)行 TB 模型對(duì)角化建議至少 16GB 內(nèi)存并考慮使用稀疏矩陣求解器。如果跑機(jī)器學(xué)習(xí)勢(shì)函數(shù)分子動(dòng)力學(xué)則需要 NVIDIA GPU建議顯存 16GB 以上。3.2 Python 環(huán)境與依賴包推薦使用 Anaconda 或 Miniconda 管理環(huán)境。創(chuàng)建一個(gè)專用的虛擬環(huán)境conda create -n moire python3.10 -y conda activate moire核心依賴包如下pip install numpy scipy matplotlib pip install ase # Atomic Simulation Environment用于構(gòu)建原子結(jié)構(gòu)版本要求以實(shí)際安裝為準(zhǔn)本文代碼基于 Python 3.10 及以上版本編寫(xiě)NumPy 2.x 或 1.26 均可。3.3 輔助工具可視化結(jié)構(gòu)可以配合 VESTA 或 ASE GUI 查看生成的原子結(jié)構(gòu)。如果需要做第一性原理驗(yàn)證可以從 PWscfQuantum ESPRESSO或者 VASP 開(kāi)始但本文不涉及 DFT 實(shí)操。如果需要處理 TB 模型推薦 PythTB 包pip install pythtb或自行用 NumPy 實(shí)現(xiàn)。4. 核心流程拆解從物理問(wèn)題到計(jì)算實(shí)現(xiàn)一個(gè)完整的莫爾材料模擬流程包含以下幾個(gè)步驟。4.1 構(gòu)建雙層的莫爾超晶格模型這是最關(guān)鍵的步驟。構(gòu)建莫爾超晶格有兩種常見(jiàn)方案方案一直接構(gòu)建大超胞。根據(jù)層間轉(zhuǎn)角確定超胞大小一個(gè)能顯示清晰莫爾周期的超胞往往包含數(shù)千到數(shù)萬(wàn)原子適合后續(xù)用機(jī)器學(xué)習(xí)勢(shì)函數(shù)做結(jié)構(gòu)弛豫和分子動(dòng)力學(xué)。方案二使用連續(xù)介質(zhì)模型Continuum Model。不顯式構(gòu)建原子結(jié)構(gòu)而是用一個(gè)周期性勢(shì)場(chǎng)描述兩層之間的相互作用。這種方法計(jì)算量極小適合快速解析平帶和電子結(jié)構(gòu)。后面的示例代碼采用這個(gè)思路。4.2 確定轉(zhuǎn)角與失配轉(zhuǎn)角決定了莫爾波長(zhǎng)。對(duì)于石墨烯晶格常數(shù) 2.46 ?(1.1^\circ) 對(duì)應(yīng)的莫爾周期約為 13.4 nm。為了在有限尺寸的超胞內(nèi)精確實(shí)現(xiàn)目標(biāo)轉(zhuǎn)角通常使用公度角commensurate angle的概念即找到一個(gè)整數(shù) ((m, n)) 使得超胞能夠在兩層晶格上同時(shí)閉合。這套方法在文獻(xiàn)中已經(jīng)很成熟感興趣的讀者可以搜索“石墨烯莫爾超晶格公度角構(gòu)建”。4.3 構(gòu)建緊束縛TB模型TB 模型是連接原子尺度與電子結(jié)構(gòu)的橋梁。在莫爾體系中可以只考慮石墨烯的 (p_z) 軌道并引入層間跳躍interlayer hopping隨距離變化的近似層間跳躍在莫爾周期上會(huì)呈現(xiàn)出局域增強(qiáng)或減弱的圖案這恰好對(duì)應(yīng)了“AA 堆垛區(qū)域”和“AB/BA 堆垛區(qū)域”的周期性分布。4.4 計(jì)算能帶與態(tài)密度在 TB 模型下能帶計(jì)算轉(zhuǎn)化為對(duì)一個(gè)大稀疏矩陣的哈密頓量 (H(k)) 進(jìn)行對(duì)角化。因?yàn)槟獱柍艽髮?shí)空間哈密頓量維數(shù)很容易超過(guò)數(shù)萬(wàn)通常選擇使用稀疏矩陣特征值求解器scipy.sparse.linalg計(jì)算最低若干條能帶或者沿高對(duì)稱路徑掃描。4.5 驗(yàn)證與可視化能帶圖是否展現(xiàn)出“平帶”是判斷模型是否捕捉到魔角效應(yīng)的核心指標(biāo)。可視化部分要輸出能帶色散圖和某一能量窗口內(nèi)的局域態(tài)密度圖。5. 完整示例與代碼實(shí)現(xiàn)下面用三個(gè)代碼示例完整演示“莫爾超晶格勢(shì)場(chǎng)可視化 — TB 模型能帶計(jì)算 — 平帶態(tài)密度驗(yàn)證”的流程。5.1 示例一莫爾超晶格勢(shì)場(chǎng)可視化這個(gè)示例用 NumPy 生成一個(gè)二維勢(shì)場(chǎng)模擬兩層六角晶格在不同轉(zhuǎn)角下的莫爾干涉圖案。# 文件路徑moire_potential.py import numpy as np import matplotlib.pyplot as plt def hex_lattice_potential(x, y, wavelength1.0): 生成一個(gè)簡(jiǎn)化的六角晶格周期勢(shì)場(chǎng)。 六角晶格可以用兩組成 60 度夾角的格矢描述 我們疊加三個(gè)方向上的余弦波來(lái)近似。 k1 np.array([1.0, 0.0]) k2 np.array([0.5, np.sqrt(3) / 2.0]) k3 k1 - k2 phase1 2 * np.pi * (k1[0] * x k1[1] * y) / wavelength phase2 2 * np.pi * (k2[0] * x k2[1] * y) / wavelength phase3 2 * np.pi * (k3[0] * x k3[1] * y) / wavelength return np.cos(phase1) np.cos(phase2) np.cos(phase3) def rotate_coords(x, y, theta_deg): 將坐標(biāo)旋轉(zhuǎn) theta_deg 度角度制。 theta np.deg2rad(theta_deg) x_rot x * np.cos(theta) - y * np.sin(theta) y_rot x * np.sin(theta) y * np.cos(theta) return x_rot, y_rot # 構(gòu)建網(wǎng)格 nx, ny 1200, 1200 x np.linspace(-15, 15, nx) y np.linspace(-15, 15, ny) X, Y np.meshgrid(x, y) # 第一層晶格 V1 hex_lattice_potential(X, Y, wavelength1.0) # 第二層晶格旋轉(zhuǎn) theta_deg theta_deg 3.0 # 使用 3 度轉(zhuǎn)角莫爾條紋清晰可見(jiàn) X_rot, Y_rot rotate_coords(X, Y, theta_deg) V2 hex_lattice_potential(X_rot, Y_rot, wavelength1.0) # 莫爾總勢(shì)場(chǎng)干涉圖案 V_moire V1 V2 plt.figure(figsize(8, 8)) plt.pcolormesh(X, Y, V_moire, cmapviridis, shadingauto) plt.title(fMoiré Pattern at θ {theta_deg}°) plt.xlabel(x) plt.ylabel(y) plt.colorbar(labelPotential) plt.savefig(moire_pattern.png, dpi150) plt.show()運(yùn)行方式python moire_potential.py如果一切正常圖像中會(huì)看到明顯的菱形或多邊形莫爾超晶格條紋。調(diào)整theta_deg為更小的角度如 1.1°莫爾周期會(huì)變得更大。5.2 示例二基于 PythTB 的魔角石墨烯能帶計(jì)算PythTB 是一個(gè)輕量的緊束縛模型工具包適合新手快速驗(yàn)證物理模型。下面的示例構(gòu)建一個(gè)雙層石墨烯的簡(jiǎn)化莫爾模型只考慮一個(gè)魔角附近的格子。# 文件路徑tb_band_structure.py # 需要先安裝 pythtbpip install pythtb import pythtb as ptb import numpy as np import matplotlib.pyplot as plt # 定義晶格常數(shù) a 2.46 # 石墨烯原始晶格常數(shù)單位埃 # 構(gòu)建雙層石墨烯的簡(jiǎn)化模型每個(gè)原胞包含4個(gè)原子 # 這里僅演示沿高對(duì)稱路徑計(jì)算能帶的方法 lat [[a, 0.0], [0.0, a]] # 4個(gè)基矢坐標(biāo)下層A、下層B、上層A、上層B orb [[0.0, 0.0], [a/2, 0.0], [0.0, a/2], [a/2, a/2]] model ptb.tb_model(2, 2, lat, orb) # 最近鄰跳躍 model.set_hop(0, 0, 1, 2.8) # 層內(nèi)最近鄰跳躍 model.set_hop(1, 1, 0, 2.8) model.set_hop(0, 0, 3, 2.8) model.set_hop(1, 2, 3, 2.8) # 層間跳躍簡(jiǎn)化不考慮轉(zhuǎn)角的空間調(diào)制 model.set_hop(0, 0, 2, 0.3) # 主要層間跳躍 model.set_hop(0, 1, 3, 0.3) # 沿 K → Gamma → M 路徑計(jì)算 path [[0.0, 0.0], [1/3, 1/3], [0.0, 0.0], [0.5, 0.0]] (kpts, labels, paths) model.k_path(path, 201) evals model.solve_all(kpts) # 作圖 fig, ax plt.subplots(figsize(6, 5)) for band in evals: ax.plot(paths, band, colorblack, linewidth1.2) ax.set_ylabel(Energy (eV)) ax.set_xticks(labels) ax.set_xticklabels([r$\Gamma$, r$K$, r$\Gamma$, r$M$]) ax.set_title(TB band structure (schematic)) plt.savefig(tb_bands.png, dpi150) plt.show()注意這是一個(gè)簡(jiǎn)化的雙層石墨烯模型并沒(méi)有真正實(shí)現(xiàn) 1.1° 轉(zhuǎn)角下的莫爾超胞。如果要真實(shí)刻畫(huà)魔角需要建立數(shù)千原子的超胞并將層間跳躍設(shè)置為隨莫爾周期呈空間調(diào)制形式。這段代碼的作用是幫助讀者理解 PythTB 的基本使用流程。5.3 示例三魔角平帶態(tài)密度驗(yàn)證當(dāng)模型正確包含莫爾周期和層間耦合調(diào)制時(shí)能帶會(huì)在費(fèi)米能級(jí)附近出現(xiàn)電子速度接近零的平帶。此時(shí)態(tài)密度會(huì)出現(xiàn)一個(gè)尖銳的峰。下面的代碼演示如何基于 TB 模型輸出態(tài)密度。# 文件路徑dos_check.py import numpy as np import matplotlib.pyplot as plt def gaussian_broadening(energies, smearing0.005): 將一組能級(jí)用高斯展寬得到態(tài)密度。 energies: 所有 k 點(diǎn)所有能帶的能量數(shù)組 smearing: 高斯展寬寬度單位 eV e_min np.min(energies) - 5 * smearing e_max np.max(energies) 5 * smearing grid np.linspace(e_min, e_max, 800) dos np.zeros_like(grid) for e in energies.flatten(): dos np.exp(-0.5 * ((grid - e) / smearing) ** 2) dos * 1.0 / (np.sqrt(2 * np.pi) * smearing * len(energies.flatten())) return grid, dos # 模擬一個(gè)包含平帶的能量集合 # 假設(shè)在 -0.05 eV 附近存在大量平帶電子態(tài) np.random.seed(42) flat_band_energies np.random.normal(-0.05, 0.003, size3000) scattered_energies np.random.normal(0.15, 0.1, size5000) all_energies np.concatenate([flat_band_energies, scattered_energies]).reshape(-1, 1) grid, dos gaussian_broadening(all_energies, smearing0.005) plt.figure(figsize(7, 5)) plt.plot(grid, dos, colordarkred, linewidth2) plt.axvline(0, colorgray, linestyle--, linewidth1) plt.xlabel(Energy (eV)) plt.ylabel(DOS (arb. units)) plt.title(Density of States with Van Hove Singularity / Flat Band) plt.savefig(dos_flatband.png, dpi150) plt.show()輸出圖中可以清楚看到在費(fèi)米能級(jí)附近有一個(gè)非常尖銳的態(tài)密度峰這便是平帶的標(biāo)志。如果直接對(duì)真實(shí)魔角石墨烯做 TB 計(jì)算會(huì)在這個(gè)位置得到類似結(jié)果。6. 運(yùn)行結(jié)果與效果驗(yàn)證6.1 莫爾圖案結(jié)果判斷運(yùn)行示例一后保存的moire_pattern.png應(yīng)該顯示一個(gè)具有明顯周期性的明暗圖案若圖案只是模糊無(wú)序的噪聲可能是網(wǎng)格分辨率不足或波長(zhǎng)與坐標(biāo)范圍不匹配。若圖案呈現(xiàn)出清晰的長(zhǎng)周期菱形條紋并且條紋間隔明顯大于單層晶格周期則說(shuō)明莫爾干涉被正確模擬。6.2 能帶圖結(jié)果判斷示例二的能帶圖只是一個(gè)簡(jiǎn)化的雙層石墨烯示意。真正驗(yàn)證魔角平帶時(shí)應(yīng)該觀察在費(fèi)米能級(jí)附近是否存在帶寬極小低于 10 meV的平帶。平帶出現(xiàn)的角度是否接近理論預(yù)言的魔角。隨著轉(zhuǎn)角偏離魔角平帶是否迅速消失并變成色散帶。一個(gè)有效的判斷方式是掃描多個(gè)轉(zhuǎn)角參數(shù)分別計(jì)算費(fèi)米面附近的帶寬然后繪制“帶寬-轉(zhuǎn)角”曲線。帶寬最小處即為該模型的魔角位置。這在計(jì)算上確實(shí)需要不少資源但對(duì)于理論學(xué)習(xí)來(lái)說(shuō)非常直觀。6.3 態(tài)密度驗(yàn)證示例三運(yùn)行成功后會(huì)看到一個(gè)尖銳的態(tài)密度峰。如果使用真實(shí) TB 模型計(jì)算還能進(jìn)一步通過(guò)投影態(tài)密度分辨各個(gè)堆垛區(qū)域AA、AB、BA對(duì)平帶的貢獻(xiàn)。通常AA 堆垛區(qū)域是平帶態(tài)密度的主要貢獻(xiàn)者這也與 STM 實(shí)驗(yàn)觀察到的莫爾超晶格亮斑位置一致。6.4 如果運(yùn)行失敗優(yōu)先檢查以下幾點(diǎn)PythTB 是否安裝成功可以執(zhí)行python -c import pythtb; print(pythtb.__version__)。NumPy 和 Matplotlib 版本是否兼容。運(yùn)行環(huán)境是否處于conda activate moire中。對(duì)于大超胞 TB 模型內(nèi)存不足時(shí)會(huì)出現(xiàn)MemoryError需要改用稀疏存儲(chǔ)或分塊計(jì)算。7. 常見(jiàn)問(wèn)題與排查思路問(wèn)題現(xiàn)象可能原因排查方式解決方案莫爾圖案周期過(guò)小看不出明顯的超晶格條紋轉(zhuǎn)角設(shè)置過(guò)大或坐標(biāo)范圍太小減小轉(zhuǎn)角例如從 3° 改為 1.1°擴(kuò)大坐標(biāo)范圍調(diào)整theta_deg參數(shù)增加np.linspace的上下限PythTB 能帶圖只有四條很窄的帶沒(méi)有平帶特征模型只是普通雙層石墨烯未包含莫爾超胞和層間調(diào)制檢查模型原子數(shù)與轉(zhuǎn)角確認(rèn)層間跳躍是否為固定值建立大超胞將層間跳躍設(shè)置為隨莫爾周期變化的空間函數(shù)態(tài)密度圖中費(fèi)米能級(jí)附近沒(méi)有尖峰平帶不在費(fèi)米能級(jí)附近或能帶過(guò)于色散檢查費(fèi)米能級(jí)位置調(diào)整填充或電子數(shù)調(diào)整平帶能量范圍或增加層間耦合強(qiáng)度構(gòu)建具有 1.1° 轉(zhuǎn)角的莫爾超胞時(shí)原子數(shù)爆炸直接使用小角度會(huì)產(chǎn)生巨大的超胞用連續(xù)介質(zhì)模型進(jìn)行快速預(yù)研或使用公度角近似在 TB 之前先用 Continuum Model 驗(yàn)證或選擇更大的公度角大超胞 TB 模型對(duì)角化內(nèi)存不足哈密頓量矩陣維度過(guò)大且使用了稠密矩陣存儲(chǔ)檢查是否有np.asarray或np.matrix轉(zhuǎn)換改用scipy.sparse.linalg.eigsh求解最低特征值或使用pythtb的并行模式弛豫后莫爾結(jié)構(gòu)變形導(dǎo)致平帶消失忽略了原子面外弛豫重建效應(yīng)對(duì)比弛豫前后的結(jié)構(gòu)差異檢查原子位移是否超過(guò) 10% 晶格常數(shù)在模型中引入面外位移項(xiàng)使用機(jī)器學(xué)習(xí)勢(shì)函數(shù)做更大規(guī)模結(jié)構(gòu)優(yōu)化與實(shí)驗(yàn)報(bào)道的魔角數(shù)值對(duì)不上模型缺少電子關(guān)聯(lián)、自旋軌道耦合、應(yīng)變等修正對(duì)比文獻(xiàn)中的 TB 參數(shù)采用文獻(xiàn)已驗(yàn)證的 TB 參數(shù)集或把魔角計(jì)算改成 LDA/DFT 的自洽驗(yàn)證8. 最佳實(shí)踐與工程建議8.1 從“離散模型”到“連續(xù)模型”的分層策略魔角莫爾材料最大的計(jì)算瓶頸是真實(shí)現(xiàn)象發(fā)生在納米尺度而真實(shí)原子間距在埃尺度這就要求超胞動(dòng)輒上萬(wàn)個(gè)原子。工程上不建議一上來(lái)就做全原子 DFT而是采用分層策略先用連續(xù)介質(zhì)模型Continuum Model解析地考察能帶和轉(zhuǎn)角依賴關(guān)系。用 TB 模型或機(jī)器學(xué)習(xí)勢(shì)函數(shù)對(duì)大超胞做結(jié)構(gòu)弛豫。最后用 DFT 或 GWBSE 對(duì)關(guān)鍵構(gòu)型做精確驗(yàn)證。這套流程在公開(kāi)文獻(xiàn)中已成為主流方案能大幅減少無(wú)效計(jì)算。8.2 參數(shù)文件與實(shí)驗(yàn)記錄管理莫爾材料計(jì)算涉及大量可調(diào)參數(shù)轉(zhuǎn)角、晶格失配、層間距、應(yīng)變、跳躍參數(shù)、外場(chǎng)等。強(qiáng)烈建議為每一個(gè)計(jì)算任務(wù)建立獨(dú)立的參數(shù) JSON 文件{ system: tBG, theta_deg: 1.1, lattice_constant_a: 2.46, hopping_intra: 2.8, hopping_inter: 0.3, model: TB, supercell_size: 55, comment: initial guess for magic angle }這樣即使幾個(gè)月后回看也能準(zhǔn)確復(fù)現(xiàn)每次計(jì)算。8.3 版本兼容與可重復(fù)性Python 的 NumPy/SciPy 在不同版本間行為有差異建議使用conda env export environment.yaml鎖定環(huán)境。對(duì)需要發(fā)布的結(jié)果可以同時(shí)記錄計(jì)算軟件的 commit hash隨機(jī)種子參數(shù)文件的 md5 值8.4 異常處理與日志記錄在批量掃描轉(zhuǎn)角參數(shù)時(shí)某一組參數(shù)可能導(dǎo)致特征值求解不收斂。建議用try/except捕獲異常并記錄日志而不是讓整個(gè)任務(wù)崩潰import logging logging.basicConfig(filenamemoire_scan.log, levellogging.INFO) for theta in theta_list: try: bands compute_bands(theta) logging.info(ftheta{theta:.3f} OK, bandwidth{np.max(bands)-np.min(bands):.6f}) except Exception as e: logging.error(ftheta{theta:.3f} failed: {e})8.5 安全邊界與生產(chǎn)環(huán)境提醒雖然這是研究型計(jì)算但涉及大規(guī)模集群使用時(shí)仍需注意不要在生產(chǎn)集群上的共享目錄直接跑大型 job容易干擾他人并造成數(shù)據(jù)混寫(xiě)。涉及商業(yè)軟件如 VASP時(shí)遵守所在機(jī)構(gòu)的許可證授權(quán)范圍。對(duì)于 TB 或 MD 計(jì)算如果會(huì)運(yùn)行幾天務(wù)必定期輸出 checkpoint以便中途失敗后從最近狀態(tài)恢復(fù)。9. 總結(jié)與后續(xù)學(xué)習(xí)方向魔角莫爾材料真正顛覆性的地方并不在于“轉(zhuǎn)角”這個(gè)動(dòng)作本身而在于它改變了材料設(shè)計(jì)的思維模式物態(tài)不再只是元素周期表中的一行一列而可以用幾何堆疊在參數(shù)空間中連續(xù)“調(diào)”出來(lái)。從計(jì)算角度這個(gè)領(lǐng)域的門檻并不在于理論物理公式有多深而在于如何高效構(gòu)建和處理超大規(guī)模的結(jié)構(gòu)模型。無(wú)論是連續(xù)介質(zhì)模型、TB 模型還是機(jī)器學(xué)習(xí)勢(shì)函數(shù)都是這個(gè)方向繞不開(kāi)的工具。對(duì)于剛開(kāi)始接觸這個(gè)方向的讀者下一步建議按順序做三件事運(yùn)行本文的莫爾勢(shì)場(chǎng)可視化代碼親眼看一次莫爾條紋隨轉(zhuǎn)角變化的演化。將 TB 模型擴(kuò)展到包含莫爾超胞至少理解“層間跳躍的空間調(diào)制”這一核心物理。閱讀魔角石墨烯的原始實(shí)驗(yàn)論文如 Nature 2018 年發(fā)表的超導(dǎo)發(fā)現(xiàn)工作將你在 TB 能帶圖中的結(jié)果與實(shí)驗(yàn)相圖對(duì)照。如果你希望進(jìn)一步深入可以關(guān)注三個(gè)延伸方向連續(xù)介質(zhì)模型與 Bistritzer-MacDonald 方程這是解析理解魔角平帶的標(biāo)準(zhǔn)武器。機(jī)器學(xué)習(xí)勢(shì)函數(shù)在莫爾材料弛豫中的應(yīng)用因?yàn)榇蟪Y(jié)構(gòu)優(yōu)化遠(yuǎn)超傳統(tǒng) DFT 的能力邊界。轉(zhuǎn)角光電子譜與 STM 實(shí)驗(yàn)數(shù)據(jù)的關(guān)聯(lián)分析這是連接計(jì)算與實(shí)驗(yàn)的關(guān)鍵橋梁。魔角莫爾材料這個(gè)方向正在快速?gòu)摹笆佟弊呦颉巴ㄓ枚S材料轉(zhuǎn)角平臺(tái)”。未來(lái)或許會(huì)出現(xiàn)更多基于轉(zhuǎn)角設(shè)計(jì)的拓?fù)涑瑢?dǎo)、量子磁性甚至量子計(jì)算平臺(tái)。對(duì)于每一位材料模擬和計(jì)算物理方向的開(kāi)發(fā)者來(lái)說(shuō)理解并掌握莫爾超晶格的計(jì)算方法正在變成一項(xiàng)越來(lái)越有競(jìng)爭(zhēng)力的基礎(chǔ)能力。建議先收藏本文等你的環(huán)境準(zhǔn)備好了把示例代碼跑一遍再回到上面的要點(diǎn)復(fù)盤一次理解的深度會(huì)完全不同。