| 查看: 1506 | 回复: 2 | ||||
| 【悬赏金币】回答本帖问题，作者月亮市集将赠送您 35 个金币 | ||||
[求助]
对pw.x计算的一些问题，望解答！
|
||||
|
在网上看到有说“pw.x约定：在同一目录，并保持outdir、prefix一致时，先运行vc-relax、relax计算，在接着的scf、nscf、bands计算会读取之前弛豫后的结构，而忽略此时的结构设置。所以，此时做scf计算，修改的地方是：将calculation=’vc-relax’改成calculation=’scf’，其他部分与上一步输入文件相同。”对此我有以下几点疑问： 1.outdir、prefix一致时我做完自洽再做非自洽是不是就把自洽产生的文件覆盖了？如果我做完vc-relax、scf、nscf，发现我scf计算有问题，要重新修改scf.in，再次提交，覆盖的文件对scf计算会产生影响吗 2.比如说我在一个目录内做了vc-relax计算，但是scf计算我想测试不同参数会产生什么影响，然后我复制了这个做完vc-relax的文件夹产生一个新的文件夹，此时scf计算还会读取弛豫后的结构吗？ 3.计算bands时要加入nbnd="",那在scf计算的时候就要加入这个参数吗？这个加的早晚有什么影响？ 非常感谢，望您不吝赐教！ 发自小木虫Android客户端 |
回复此楼» 猜你喜欢
售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急
已经有3人回复
-
售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急
已经有6人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.54,科目齐全,可+急
已经有3人回复
-
售一区SCI文章T0P，我:8O.551.O54,科目全,可十急
已经有5人回复
-
售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急
已经有6人回复
-
售SCI一区T0P文章，我:8O.55.1.O.5.4,科目齐全,可+急
已经有5人回复
-
售SCI一区T0P文章，我:8O.55.1.O.5.4,科目齐全,可+急
已经有4人回复
-
售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急
已经有7人回复
-
售SCI文章，我:8O.5.5.1O.54,科目全,可十急
已经有6人回复
-
售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急
已经有4人回复
卡开发发
专家顾问 (著名写手)
Ab Initio Amateur
-
专家经验: +224 - 1ST强帖: 2
- 应助: 47 (小学生)
- 金币: 3575.8
- 散金: 14
- 红花: 223
- 帖子: 2477
- 在线: 1254小时
- 虫号: 1369423
- 注册: 2011-08-16
- 性别: GG
- 专业: 金属材料的磨损与磨蚀
- 管辖: 第一性原理
|
1、一般来说from scratch的情况不会对scf产生影响，但如果restart可能会有影响，或者读取波函数或势的情况会造成影响，我自己习惯上对nscf或者band计算之前会把scf目录备份一下防止出现一些乱七八糟的情况。 2、不会，你需要把out文件当中最后的结构手动更新到新的inp的结构部分，这样才是用的驰豫后的结构进行scf。 3、scf是否加入看场景。这里得先提到一个概念，就是平面波的基尺寸很大，所以一般不会求解全部的本征值，而是采用迭代对角求解较低的那些能级，nbnd就是控制求解能级数目的参数。然后存在如下三个场景： （1）对于绝缘体来说，基本上求解到价电子最高占据能级即可，即fix occupancy. （2）对金属或者窄带半导体或者磁性体系，一个是因为引入展宽（电子有限温度）来加速BZ积分收敛但会引入电子的分数占位，一个是自旋极化情况会导致alpha或者beta电子会占据到更高的能级，这样要的nbnd就会更大一些，程序尽管有预设，但也有不太正常的时候需要手动调。 （3）上述两种情况自洽后需要计算能带的时候，我们需要看到导带甚至更高能级的信息，这个时候会通过设置更大的nbnd求解那些较高的能级，为了提高效率采用非自洽（即一般nscf或band的任务）读取当前密度通过Harris-Foulkes泛函估计能级而不重新自洽，这在场景在(P)DOS和能带计算中比较常见。 |
» 本帖已获得的红花（最新10朵）
月亮市集
2楼2024-02-01 13:30:05
3楼2024-02-03 15:01:26
