| 查看: 3197 | 回复: 14 | ||||
[交流]
【已解决】【求助】为何PWscf不能得到类似于CASTEP那样的收敛测试曲线呢？
|
准备转入PWscf，因为这个免费 于是例行公事一样的先对模型进行了收敛测试，使用PBE的赝势，同时使用了CASTEP在完全相同的结构、K点下进行单点能计算。结果……CASTEP下得到了很圆满的收敛测试曲线（抛物线）：  而使用PWscf，得到的居然是个直线……：  有谁也曾经碰到过这样的问题么……不确定收敛参数的PWscf的计算结果会有意义么…… [ Last edited by souledge on 2011-1-2 at 22:20 ] |
回复此楼» 收录本帖的淘帖专辑推荐
 第一原理资料汇编 |
» 猜你喜欢
售T0P一区SCI文章，我:8O5.51.O.54,科目齐全,可＋急
已经有3人回复
-
售T0P一区SCI文章，我:8O5.51.O.54,科目齐全,可＋急
已经有6人回复
-
售T0P一区SCI文章，我:8O5.51.O.54,科目齐全,可＋急
已经有4人回复
-
售T0P一区SCI文章，我:8O5.51.O.54,科目齐全,可＋急
已经有6人回复
-
售T0P一区SCI文章，我:8O5.51.O.54,科目齐全,可＋急
已经有6人回复
-
售T0P一区SCI文章，我:8O5.51.O.54,科目齐全,可＋急
已经有4人回复
-
国社科通讯评审
已经有6人回复
-
面上的水很深
已经有21人回复
-
小木虫没落了，除了祈祷帖子，几乎看不到有价值的帖子
已经有15人回复
-
5个%2F，非固定段22个字母
已经有9人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 还是决定了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/1892
-
重磅分享｜全链路AI材料研发闭环平台PULSE即将发布，重构新材料研发范式
+2/126
-
【求职】科研助理/生物技术相关岗位 意向工作地点 ：河南郑州
+1/85
-
不用懂统计，用初中数学锤爆这篇造假学位论文：10+处铁证全解析
+1/83
-
深圳大学黄妍斐特聘教授课题组专任教师（助理教授）招聘
+1/79
-
澳大利亚国立大学【全奖】：AI+钙钛矿光伏方向（长期有效）
+2/76
-
【重磅招聘】小巨人企业西顿高薪招聘水性树脂研发博士、硕士（广州）
+1/58
-
祈福，求顺利！
+1/33
-
西安交通大学补亚忠课题组招聘博士生
+1/28
-
祈福C口面上！！！！
+1/24
-
希望面上有好消息:D:D:D
+2/12
-
华中农业大学植物科学技术学院李国田老师课题组劳动聘用制科研助理招聘公告
+1/12
-
中国科学院宁波材料技术与工程研究所招收研究生
+1/8
-
人生失去与收获
+1/7
-
诚邀生物医用材料领域学者加入Biofunctional Materials青年编委（Scopus收录）
+1/5
-
Scientific Journal of Intelligent Systems Research 期刊简介
+1/5
-
出国留学的目的是什么？
+1/4
-
东北大学理学院力学李健课题组招收2027级硕士推免研究生
+1/2
-
福州大学吴衔誉教授课题组 诚招2027届博士生
+1/2
-
北理工集成电路杰青团队 | 诚招科助理
+1/1
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
