资源描述：
《石墨烯 graphene MS建模方法》由会员上传分享，免费在线阅读，更多相关内容在行业资料-天天文库。
1、1、打开materialstudio，新建一个工程，导入石墨graphite.msi（也可以自己build，然后添加原子）。2、build->makep1（目的是消除对称性，这样才能够删除一层原子）。3、删除一层原子（选中原子->delete）。4、修改晶格参数：build->crystal->rebuildcrystal，设置方位角，，5、构建supercell（方便掺杂，也为了好看）：build->symetry->supercell，构建一个5x5x1的超原胞。6、cleavesurface（为了能够添加真空层）：build->surface->cleavesurface,(h,k,
2、l)改为（0,0，-1）7、添加20埃真空层（添加真空层是为了减小层与层之间的影响，至少20埃，大点没关系，最多是计算时时间长一点）：build->srystal->buildvacuum。构建好后，模型如下：两种模型的建立方法：第一种，导入软件内置模型执行file–import–structure–ceramics–graphite.msi，获得双层石墨烯，层间距为0.34nm，将其扩充为6层，选定一层，将其移动到模型正中央，模型厚度为0.68*3nm；第二种方法，建立晶胞，选择模型为第183型，设置参数为2.46、2.46和3.4，然后将碳原子添加进去，设置坐标为0.333、0.667
3、和0.500，获得厚度为0.34nm的晶胞，将其扩充为6层，因此它的厚度与第一种一样。现在要确定两种模型的结点个数，为使体积接近，分别将其扩充为145和128个结点。如图，显而易见，第一种模型边沿布满结点，而第二种模型边沿没有结点。为使模型稳定，对它们初步先进行几何结构优化。优化以前，键角都是120°键长均为0.142nm。几何结构优化后，键长和键角均发生了一些轻微变化。(模型一)(模型二)导入石墨结构后，cleavesurface，取石墨的C方向(001)，选合适的thickness和position，使之只有一层原子（比如top:1.0，thickness0.6），得到表面后再buil
4、dvacuumslab，选thickness（比如20A），slabposition可以选负的（比如-10A），这样就得到了，你还可以重新定义二维晶胞的晶格参数(build->symmetry->redefinelattice,比如选B为-120，晶格就变成长方形的了，当然也可以在六方晶格的supercell上删掉一些原子得到长方形的supercell)，使之更适合你的需要。选择Import
5、Structures
6、ceramicsandimportgraphite.msi.然后选择,Build
7、Symmetry
8、MakeP1，删除其中一层,把剩下的一层移到中间,然后选择Build
9、Symm
10、etry
11、FindSymmetry…
12、FindSymmetry
13、ImposeSymmetry.接下来选择Build
14、Symmetry
15、Supercell，创建5*5*1的石墨烯超晶胞.接下来选择Build
16、Crystals
17、RebuildCrystal....把C值增加到15，点击rebuildcrystal.SampleText 优化步骤.点击CASTEP按钮Calculation打开CASTEPCalculation对话框.把Task从Energy改变为GeometryOptimization.点击More…打开CASTEPGeometryOptimization对话框.把Qualit
18、y改为Fine,然后勾选上theOptimizeCell复选框.接下来选择Electronic,把SCFtolerance改为Fine，然后点击More…打开CASTEPElectronicOption对话框.在Basis标签下勾选Usecustomenergycutoff并键入300.在k-point下点击theCustomgridparametersradio.在theGridparametersfield下,把a、b、c依次设置为3、3、1，并关闭对话框.最后,点击Run。布里渊区积分路径：diracpoint应该出现在K点。但是当构建超胞的时候，由于去掉了体系的对称性（makeP1
19、），K点不再是MS计算时默认的高对称点了。因此，在计算BandStructure时，MS的选项里没出现K点。有两种解决方法。一是用MS先确定模型的对称性，看看能不能使K点出现。（对于我选取的n×n的平行四边形超胞是可以的。）二是计算BandStructure时直接指定K点的左边。K点坐标是（-0.333，0.667，0）.我的n×n超胞G点在C原子上，能带计算是G-->K，能带在K点相交于fermi面
