| 查看: 3440 | 回复: 14 | ||||
[交流]
【已解决】【求助】为何PWscf不能得到类似于CASTEP那样的收敛测试曲线呢？
|
准备转入PWscf，因为这个免费 于是例行公事一样的先对模型进行了收敛测试，使用PBE的赝势，同时使用了CASTEP在完全相同的结构、K点下进行单点能计算。结果……CASTEP下得到了很圆满的收敛测试曲线（抛物线）：  而使用PWscf，得到的居然是个直线……：  有谁也曾经碰到过这样的问题么……不确定收敛参数的PWscf的计算结果会有意义么…… [ Last edited by souledge on 2011-1-2 at 22:20 ] |
回复此楼» 收录本帖的淘帖专辑推荐
 第一原理资料汇编 |
» 猜你喜欢
国自然评审意见
已经有16人回复
-
售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急
已经有6人回复
-
售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急
已经有7人回复
-
售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急
已经有6人回复
-
有多少人是今天查系统知道结果的？
已经有18人回复
-
2026年叶企孙基金
已经有7人回复
-
我就是申请一个面上项目而已，这评审意见是按照杰青的条件评的吧？
已经有6人回复
-
售SCI文章，我:8O.5.5.1O.54,科目全,可十急
已经有4人回复
-
为什么资助数各大高校都创新高，自己申请怎么就这么难
已经有13人回复
-
国自然面上复盘~欢迎讨论
已经有14人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 还是决定了13年考研，希望得到点大家的经验，谢谢！
已经有229人回复
- 纤维素酶的菌液怎么处理之后得到纤维素酶粉末
已经有5人回复
- 求绳子的曲线方程
已经有6人回复
- 收敛精度
已经有7人回复
- 奇怪的磁热曲线
已经有10人回复
- 请教DSC 曲线问题，有两个尖锐的吸热峰，一个放热峰，核磁很干净！
已经有17人回复
- CASTEP模块中如何通过施加应变来得到应力，并形成应力应变曲线，详细实施步骤
已经有10人回复
- CASTEP 能量曲线上升 但是收敛了
已经有13人回复
- 【求助成功】castep 收敛性测试
已经有22人回复
- 【求助】CASTEP计算弹性常数得到应力应变曲线
已经有17人回复
- 【求助】castep收敛性测试问题
已经有20人回复
- 【求助】PWSCF计算介电极化时为何出现如下提示？
已经有19人回复
- 【求助】Pwscf计算结果与castep的比较
已经有6人回复
- 【求助】请教PWSCF晶体建模和自洽计算，不收敛
已经有21人回复
2楼2011-01-02 13:19:59
|
是使用的Bohr啊～，其中，参数XryParAU是从798976开始，以步长3779变化的变量 &CONTROL calculation = 'scf' , restart_mode = 'from_scratch' , wf_collect = .false. , outdir = '/home/vm/PW_Out/' , wfcdir = '/home/vm/PW_temp/' , pseudo_dir = '/home/vm/PWscf/pseudo/' , lkpoint_dir = .false. , disk_io = 'default' , verbosity = 'default' , tstress = .false. , tprnfor = .false. , tefield = .false. , dipfield = .false. , lelfield = .false. , lberry = .false. , / &SYSTEM ibrav = 2, celldm(1) = ${XryParAU}.14d-5, nat = 8, ntyp = 2, ecutwfc = 29.3994472 , / &ELECTRONS electron_maxstep = 150, startingpot = 'atomic' , startingwfc = 'atomic+random' , mixing_mode = 'plain' , diagonalization = 'david' , diago_full_acc = .false. , tqr = .false. , / ATOMIC_SPECIES Ti 42.90000 Ti.pbe-sp-van_ak.UPF C 12.01100 C.pbe-van_ak.UPF ATOMIC_POSITIONS crystal Ti 0.000000000 0.000000000 0.000000000 Ti 0.500000000 0.000000000 0.500000000 Ti 0.500000000 0.500000000 0.000000000 Ti 0.000000000 0.500000000 0.500000000 C 0.500000000 0.500000000 0.500000000 C 0.000000000 0.000000000 0.500000000 C 0.000000000 0.500000000 0.000000000 C 0.500000000 0.000000000 0.000000000 K_POINTS automatic 10 10 10 0 0 0 |
3楼2011-01-02 13:27:20
4楼2011-01-02 14:31:00
咋一开始用PWscf就这么纠结呢……5楼2011-01-02 15:38:37
6楼2011-01-02 15:51:22
7楼2011-01-02 15:59:00
8楼2011-01-02 16:18:37
9楼2011-01-02 20:59:27
10楼2011-01-02 21:03:04
11楼2011-01-02 21:39:46
|
我在想是不是建立晶体模型的时候出现了问题呢？ 看了一下PWscf给的例子，Si的结构应该如下图：  理论上，应该是FCC结构，Si的位置在(000)，(0.5,0.5,0)，(0,0.5,0.5)，(0.5,0,0.5)，(0.75,0.25,0.25)，(0.25,0.75,0.25)，(0.25,0.75,0.25)，(0.75,0.75,0.75)，同时，nat=8，ntyp=1。 但是例子中给的模型是nat=2，ntyp=1，Si的位置只出现在(000)和(0.25,0.25,0.25)这两个位置上。这是用的什么平移方法得到的呢？ 是不是也要用同样的方法来表示TiC呢？TiC结构如下，是否应该是TiC也是nat=2，ntyp=2，两个原子位置分别是Ti(000)，C(0.5,0.5,0.5)呢？结晶学学得比较糊涂，还望大家不吝指点……  |
12楼2011-01-02 21:40:58
13楼2011-01-02 21:44:04
确定了，的确是因为没有考虑FCC的平移条件，使用P1的原子坐标的所致。更改为nat=2后，计算的结果跟CASTEP相比差别也不大了，只有4.368 A那个点不知道为啥，突然有点走火，但是总体问题不大了，谢谢大家
14楼2011-01-02 22:18:36
15楼2013-03-29 16:33:19
投票:
