1、打開materialstudio,新建一個工程，導入石墨（也可以自己build，然后添加原子）。2、 build->makepl（目的是消除對稱性，這樣才能夠刪除一層原子）。3、 刪除一層原子（選中原子->delete）。4、 修改晶格參數：build->crystal->rebuildcrystal,設置方位角，，5、 構建supercell（方便摻雜，也為了好看）：build->symetry->supercell，構建一個5x5x1的超原胞。6、 cleavesurface（為了能夠添加真空層）：build->surface->cleavesurface,（h,k,l改為（0,0，-1）7、 添加20埃真空層（添加真空層是為了減小層與層之間的影響，至少20埃，大點沒關系，最多是計算時時間長一點）：build->srystal->buildvacuum。構建好后，模型如下：兩種模型的建立方法：第一種，導入軟件內置模型執行file-import-structure-ceramics-，獲得雙層石墨烯，層間距為，將其擴充為6層，選定一層，將其移動到模型正中央，模型厚度為*3nm；第二種方法，建立晶胞，選擇模型為第183型，設置參數為、和，然后將碳原子添加進去，設置坐標為、和，獲得厚度為的晶胞，將其擴充為層，因此它的厚度與第一種一樣。現在要確定兩種模型的結點個數，為使體積接近，分別將其擴充為145和128個結點。如圖，顯而易見，第一種模型邊沿布滿結點，而第二種模型邊沿沒有結點。為使模型穩定，對它們初步先進行幾何結構優化。優化以前，鍵角都是120°鍵長均為。幾何結構優化后，鍵長和鍵角均發生了一些輕微變化。（模型一）（模型二）幾何結構優化后，鍵長和鍵角均發生了一些輕微變化。（模型一）（模型二）導入石墨結構后，cleavesurface，取石墨的C方向（001），選合適的thickness和position，使之只有一層原子（比如top:，thickness），得到表面后再buildvacuumslab，選thickness（比如20A），slabposition可以選負的（比如-10A），這樣就得到了，你還可以重新定義二維晶胞的晶格參數（build->symmetry->redefinelattice,比如選B為-120，晶格就變成長方形的了，當然也可以在六方晶格的supercell上刪掉一些原子得到長方形的supercell），使之更適合你的需要。選擇Import|Structures|ceramicsandimport.然后選擇,Build|Symmetry|MakeP1,刪除其中一層,把剩下的一層移到中間，然后選擇Build|Symmetry|FindSymmetry...|FindSymmetry|ImposeSymmetry.接下來選擇Build|Symmetry|Supercell,創建5*5*1的石墨烯超晶胞.接下來選擇Build|Crystals|RebuildCrystal....把C值增加到15,點擊rebuildcrystal.SampleText優化步驟.點擊CASTEP按鈕Calculation打開CASTEPCalculation對話框.把Task從Energy改變為GeometryOptimization.點擊More…打開CASTEPGeometryOptimization對話框.把Quality改為Fine,然后勾選上theOptimizeCell復選框.接下來選擇Electronic,把SCFtolerance改為Fine，然后點擊More…打開CASTEPElectronicOption對話框.在Basis標簽下勾選Usecustomenergycutoff并鍵入300.在k-point下點擊theCustomgridparametersradio在theGridparametersfield下，把a、b、c依次設置為3、3、1，并關閉對話框.最后，點擊Run。布里淵區積分路徑:臣］BriLLgXijJl.如涂鋌h二 一'一ResetBrillouinzonepathFromTo臨］日.日前6.-0000.000.-0J3-30.6670.000拎-0.3330-8670.000F:0.0000.5000.000F:0.0000.5000.000g：0.000o.aooo.ooo廠pispla^jreciprocal 疾；X?￡?Helipdiracpoint應該出現在K點。但是當構建超胞的時候，由于去掉了體系的對稱性（makeP1），K點不再是MS計算時默認的高對稱點了。因此，在計算BandStructure時，MS的選項里沒出現K點。有兩種解決方法。一是用MS先確定模型的對稱性，看看能