问题如下
1. CASTEP caculation对话框中的TAST中的\, Geometry Optimization, Dynamic, Elastic Constants\各有什么计算意义,尤其是Energy和Geometry的计算上的区别。 答:
Energy指的就是单点能的计算,或则说能量的计算。
Geometry optimization指的是结构优化,而优化的依据主要就是能量,所以,在Geometry optimization的过程中,CASTEP在调整结构,每得到一个结构就会计算它的能量,最终找到一个能量最低的结构。
CASTEP在执行Energy和Geometry optimization计算时都会计算体系的能量,而我们知道,体系的所有性能,包括电子结构、力学性质、热力学性质、光学性质等等,所有我们在properties里面可以勾选的性质选项都是在体系能量确定之后得到的,所以在CASTEP的操作中,当你选择Energy或者Geometry optimization时,都可以同时在properties里面勾选要计算的性质,不过,当选择Geometry optimization时,程序计算的是优化得到的稳定结构所对应的性质。CASTEP的TASK中还有一个properties的选项,这个选项是用来直接计算体系性质的,但前提是你已经做过Energy或Geometry optimization计算,体系的能量已经确定,这时候可以通过properties直接计算体系的各种性质,这个选项主要是为了方便使用者,不必每次计算性质都要从新计算体系的能量。
Dynamic是做动力学模拟,也就是基于牛顿运动方程研究体系中各个原子在指定热力学条件下如何运动,它与常说的分子动力学相比,最大的区别在于原子间相互作用是通过量子力学计算,或则说求解薛定谔方程确定的,而通常的分子动力学方法是通过基于经验参数的力场来描述原子间的相互作用。 Elastic constant是做体系的弹性力常数计算,属于力学性质范围。
2. Functional: LDA(CAPA); GGA(PW91,PBE, RPBE, PBSOL), HF, HF-LDA等各有什么优势,如何选择。 答
各种交换泛函的优势或则针对性,你可以参考每种泛函的文献,这一点MS的help文件中都有明确的标注,对于固体材料,特别是你们所研究的一些过渡金属氧化物材料,WC和PBESOL是首选。
3. 我建了BaTiO3这样一个晶胞,以晶格常数4.018nm作为初始参数,然后进行了几何优化,但计算完毕后,这个参数没有任何变化,是计算没有任何误差,还是我的设置有问题。计算文件在附件中。 答
如果你要优化这个参数,请在Geometry optimization中选择optimize cell(当选择Geometry optimization后,点击其后的MORE按钮,在跳出的面板中可以找到这个选项),不需要把对称性变为P1。
4. 还是BaTiO3这样一个晶胞 ,我把模型建好后,想选择,
Build | Symmetry | Primitive Cell,但“Primitive Cell”,为灰色,不可选,为什么 ?
答
BaTiO3的primitive cell与conventional cell相同,你可以看到conventional cell中只包含一个分子式的BaTiO3,也就是1个Ba;1个Ti;3个O。
5. Symmetry 下的unbuild crystal, Nonperiodic, Superstructure, Make P1, Redefine options各有什么作用? 答
Unbuild crystal: 得到最小非对称单元的结构
Nonperiodic: 去掉结构的周期性,形象地说就是把盒子去掉。
Superstructure: 构建超晶胞结构,也就是扩大最小重复单元(或则说晶胞) Make P1: 去掉晶体结构中的所有点对称操作,只保留其平移对称性 Redefine lattice: 重新定义晶胞中基矢的方向
6. 我还做了一个,BaTiO3和SrTiO3的对半掺杂,但是几何优化后得到的能带曲线一篇文献中的图(见附件)相变,我得到图曲线变化过于单调,不知是哪里的设置出了问题。计算文件也在附件中. 答
没有任何问题,你只需要调整一下显示即可。在快捷工具栏中找到并选择放大镜工具(3d viewer zoom mode),然后按住鼠标左键,在能带结构图上滑动,调整能带结构图的显示。
不过,你得注意,由于你把BaTiO3的对称性调整为P1,此时MS在计算能带时的默认路径与你给我们的图不同,不过确实不存在“曲线变化过于单调的问题”。