不同压强下ZnO的声子谱计算及其收敛性测试

通常对于一个经费不是那么充足，即没什么经费的计算模拟课题组来说（我没说我组，没说），究竟整几台服务器是个大问题。服务器的使用又经常会出现以组会为周期特征，以审稿意见为诱发因子的使用高峰期，和相应的低峰期，难以互补。故而超算就可以作为很好的润滑剂来使用。

北鲲云能够基本满足一切我等对超算的幻想，对于导师来说，最重要的肯定是便宜，便宜，和便宜；对于学生来说，最重要的则是简单，简单，傻瓜式的简单。前者我们可以在后续的计算中管窥一斑，后者，大家点开北鲲云的主页注册一下，登录控制台看一看就知道了。大部分常用的科研软件均皆预装，点开即用；所有过程都是点一点、拖一拖，无脑式完成。并且注册所赠送的试用金，真的够计算出不少东西的。这个刷新了我对超算所存有过的偏见。

在本文中，我来做一系列简单的声子谱和分析，来为大家演示在宝藏云中一站式科研的全过程。对于材料模拟来说，确定它的动力学稳定性是非常关键的一点。声子谱G点下的小小虚频，也是无数人的噩梦。但实际上即便是计算“出错”的声子谱，也是包含很多信息的，这一点我还尚在学习，将一点心得与此计算过程一切记录下来。我们来计算三维体材料ZnO在常压下和19Gpa静水压下的声子谱，每一个声子谱计算对他进行k点的测试，最后我们将声子谱对不同原子进行投影，进一步分析他的晶格振动情况。

相当的酷炫，也很流畅，对于我来说，MS只用来建模，所以装一个几个G的东西，很划不来，调用一个4核节点，一小时3毛2。

1，

a) 常压优化INCAR 与 晶格常数

晶格常数：3.165（a、b），5.106（c）

b) 190Gpa INCAR 与 晶格常数。

晶格常数：3.165（a、b），5.106（c）

我专门上vasp论坛看过，静水压的单位KB = 0.1Gpa。

2， 平平无奇的声子谱计算：

对两种晶格，按网传的10埃米原理，进行3-4-2的扩胞，共96个原子。对得到的超胞在1x1x1,2x2x2,3x3x3点网格下进行声子谱计算，以查看其收敛情况。声子谱输入文件如下：

在进行1x1x1（即单G点）计算时，我们可以采取单G点版本来计算，以缩短计算时间。只要简单的把提交任务脚本中的vasp_std 修改为 vasp_gam。(记得在计算更大k点的时候要调回来)，另外涂抹方案如-