第一原理计算 CASTEP 实战守则.docVIP

CASTEP 实战守则 ? 模型选定 我们在进行材料物理模拟所需要做的第一步 (也是很重要的一步) 是模型的选定或建构。CASTEP 虽然内建了很多功能来预测晶胞参数 (边长，夹角) 与原子位置，但仍然仰赖使用者告诉它 要进行计算的系统是什么。在选定模型时，我们需切记如果系统内原子太多或是超晶胞体积太大，则计算量都会以平面波数的 3 次方增加到计算机难以负荷或使用者难以等待的程度，]此，我时时应考虑设计出一个足以表现出我们所想要研究的物理，而却又能使所有采用的超晶胞越少越好的模型。 如果研究上涉及到一系列原子数不同的大小系统需要做计算，最好能先计算小结构，不要一开始就送入大结构到计算机中。 如果要仿真的系统是含有杂质，则单位晶胞必须进一步放大成超晶胞以便使化学成分里的分数变成整数，因此晶胞会变得很大。 在某些特殊的情况，相互取代的元素种类是很类似的，(即在化学行为上类似)，则下一个版本的CASTEP会提供一种叫做虚拟晶体近似 (Virtural Crystal Approximation，VCA) 的方法，则模型里面的原子就可以指定成如0.3A元素加0.7B元素这种样子，因此总可以以最小晶胞来做计算的模型。但这种方法的精确度通常只适用于合金材料，故要小心使用。 模型的选定有许多人为抉择会含在里面，例如表面计算的层数，因此有些情况也需要进行所谓的收敛性测试。 ? ? ? Vps (poseudopotential) 选择 Vps选单： MS接口的设定是选用USP优先于NCP，USP有加速计算与减少内存使用的效果，其精确度也与NCP (norm-conserving potential) 相当。至于什么时候使用NPC呢 (就是在poseudopotential 选项中那些延伸文件名是 .recpot者)，使用到NCP的场合有： 1.某些CASTEP计算的功能尚未支持到USP，因此需要选用NCP。 2.为了要与已经发展的文献比较或进行验证 3.对计算的结果存疑者，能提供『多一种选择』 ( 注：至于延伸档名是*.psp者也是属于norm-conserving的一种，是TM potential，在文献上也常被使用，但所需的截止动能较高，因此计算代价较大。) Vps information： CASTEP/Materials Studio套件所包含的Vps都经过测试并纪录其在使用不同E-cut(截止功能)情况之下的总能收敛程度，其中后者会在使用者选定计算质量是coase，medium，fine时，依照收敛度的需求决定CASTEP进行计算时实际要用多高的平面波截止功能(及要用多少个平面波求作波函数及电荷密度的傅立业展开)。而前者以批注的方式纪录在Vps的文件头，可以使用一般的文字编辑器来阅读。我们只到在Vps数据库存放的地点 (C:\Program Files\Accelrys\MS Modeling 3.0\Data\Resources\Quantum\Castep\Potentials) 就可以直接打开阅读。 coarse medium fine： 使用coarse之E-cut所算出之总能，除以系统总原子数后，每个原子的误差在1.0 eV之内，同理：medium→0.3 eV之内，fine→0.1 eV之内。precise则是fine的E-cut值在加10% 同元素但不同的Vps： 同元素但不同的Vps，含有什么不同的选择？(1)近似半径不同 (如O_00.recpot用于氧化物，O_01.recpot用于分子) (2)价电子数不同(如13个价电子的Ga_00.recpot可用于GaN，但更简单的GaAs只需3个价电子的Ga_01即可) 最近 user 遇到的例子：以 norm-conserving potential 计算 InN ，选用了预设的In_00.recpot ，但它是为 金属态铟设计的 Vps。应选 In_01.recpot 才定合用于InN者。 ? ? ? ? 精确性测试 (accuracy test) 选择一个元素成份结构排列接近所要研究的系统，但是其单位晶胞成份是分子结构却是小得多，来进行物性的预测并与实验比较 (通常是晶胞参数或分子键长等简单的性质) 此一测试结构的选取上，往往也与什么结构恰好有实验值有关.自然在众多可能选取的已知实验所得的性质里，尽可能取与我们研究目标有关的物理量，则将来所获得结果的可信度也就越高。 ? ? ? ? ? 收敛性测试 举凡计算物理中利用数值方法的研究，其结果都不能具有数学上相当于定理证明的效力。任何计算物理的研究为了要确保其可靠性，那需要针对所有的人为选择参数进行收敛性测试，也就是说，要以要好的参数精细度精细度去计算同一个我们所要的物理量，直到所需的物理量在我们所需要的精确度范围之