zyj19831206 发表于 2016-7-13 22:37
是不是还需要编程什么的？那个程序只是得到平均值？

里面的z应该不是变量，否则Varvg(z)就不是个定值了。最简单的处理就是直接对z∈a or b进行平均：

c=a or b
Vavrg(z∈c)=1/c*∫_(z-c/2~z+c/2) Vavrg(z')dz={∫_(z-c/2~z+c/2) Vavrg(z')dz'}/{∫_(z-c/2~z+c/2)dz'}
如果按照△z等间隔的矩形方法来数值积分，离散化之后
{∫_(z-c/2~z+c/2) Vavrg(z')dz'}→∑_(z-c/2~z+c/2) Vavrg(z')△z
{∫_(z-c/2~z+c/2)dz'}→∑_(z-c/2~z+c/2)△z=Nc△z，Nc是z-c/2~z+c/2的格点数目
Vavrg(z∈c)=∑_(z-c/2~z+c/2) Vavrg(z')△z/Nc△z=∑_(z-c/2~z+c/2) Vavrg(z')/Nc，其实也就是第一次平均得到的曲线在z∈c取平均，间隔△z不太大的情况下，误差应该不会很大。要想结果好点的话，可以写个程序，或者excel也能实现（以前见到过有用excel做RK的）。

卡开发发 发表于 2016-7-13 18:46
哦，这样就很清楚了，划分是按照原有bulk的尺寸来估计的（也就是a和b的值）。知道划分的话只要对a和b的两 ...

是不是还需要编程什么的？那个程序只是得到平均值？

zyj19831206 发表于 2016-7-13 17:36
李明宪的书里运用的是moving-slab averaging这个工具来实现的，不知道具体怎么搞。

哦，这样就很清楚了，划分是按照原有bulk的尺寸来估计的（也就是a和b的值）。知道划分的话只要对a和b的两个区域分别对Vavrg_a(z)和Vavrg_b(z)对z数值积分除以高度就行，当然，粗略一点的话直接数值平均也成。

zyj19831206 发表于 2016-7-13 16:30
那在异质结中，不同的物质的平均势能怎么获得？貌似运用这个程序只能得到整体的平均势能。

这我不清楚，主要不知道两个物质的分界面怎么进行划分电势，如果知道划分的方式，分界面隔开的两部分各自做平均应该就行了。

卡开发发 发表于 2016-1-24 16:31
需要将电势文件在Analysis导出到xsd文件当中，然后建立一个Slice，名称得叫"Slice1"，然后再运行脚本。你可 ...

那在异质结中，不同的物质的平均势能怎么获得？貌似运用这个程序只能得到整体的平均势能。

需要将电势文件在Analysis导出到xsd文件当中，然后建立一个Slice，名称得叫"Slice1"，然后再运行脚本。你可能没有把静电势导出并建立Slice。这个是以前回答别人的http://emuch.net/bbs/viewthread. ... =1&target=blank可以做参考，最多再多对z做个平均。

1、xsd文件处于激活状态下，在castep的calculation选项卡中设置好计算参数以及在job control中勾选retain server files；
2、点击calculation选项卡中的files-save files，产生计算工作目录，激活*.param文件，在其中加入WRITE_FORMATTED_POTENTIAL: TRUE；
3、激活*.param文件的状态下在castep的calculation选项卡中files-run files，计算开始进行；
4、计算完成后在MS中job explorer中右键计算完成的任务-remote control，会自动激活网页浏览器，里面有计算产生的所有文件，找到*..pot_fmt右键另存为到工作目录即可。