资源描述:
《石墨烯 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面