| 查看: 3194 | 回复: 14 | ||||
[交流]
【已解决】【求助】为何PWscf不能得到类似于CASTEP那样的收敛测试曲线呢？
|
准备转入PWscf，因为这个免费 于是例行公事一样的先对模型进行了收敛测试，使用PBE的赝势，同时使用了CASTEP在完全相同的结构、K点下进行单点能计算。结果……CASTEP下得到了很圆满的收敛测试曲线（抛物线）：  而使用PWscf，得到的居然是个直线……：  有谁也曾经碰到过这样的问题么……不确定收敛参数的PWscf的计算结果会有意义么…… [ Last edited by souledge on 2011-1-2 at 22:20 ] |
回复此楼» 收录本帖的淘帖专辑推荐
 第一原理资料汇编 |
» 猜你喜欢
小木虫没落了，除了祈祷帖子，几乎看不到有价值的帖子
已经有15人回复
-
面上的水很深
已经有18人回复
-
国社科通讯评审
已经有5人回复
-
5个%2F，非固定段22个字母
已经有9人回复
-
生命口C13面上会评时间
已经有11人回复
-
内心匮乏
已经有16人回复
-
生命口面上会评结束了吗大家有知道结果吗？现在还没有结果正常不？
已经有9人回复
-
投英文期刊如何查重
已经有6人回复
-
不要再数国自然申请书的 filecode 的分隔符个数了
已经有21人回复
-
伯胺上甲基求助
已经有7人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 还是决定了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/1896
-
工材电工至今没任何消息，散金币以求心安吧~
+2/1136
-
北京工业大学材料科学与工程学院李建荣教授课题组诚聘电池相关领域应届毕业生和博士后
+1/578
-
散金币，祝福自己
+5/515
-
坐标山东，诚征女友
+1/411
-
【哈工大-材料学院-新材料及成形技术-招收博士、硕士研究生】长期有效
+2/118
-
【求职】科研助理/生物技术相关岗位 意向工作地点 ：河南郑州
+1/85
-
BBC儿童英语科普《Do you know》高清英字，培养思考习惯，非常好的启蒙科普片！
+1/85
-
不用懂统计，用初中数学锤爆这篇造假学位论文：10+处铁证全解析
+1/83
-
年薪40W+项目论文绩效，广州医科大学招收骨科、生物材料相关背景博士后2名
+1/81
-
复旦大学李鹏课题组/香港理工大学马凯凯联合培养博士
+1/65
-
【重磅招聘】小巨人企业西顿高薪招聘水性树脂研发博士、硕士（广州）
+1/58
-
北京理工大学杰青沈国震团队招聘柔性智能电子博士后
+1/33
-
基金申请
+1/7
-
7月15日-21日暨南大学2027年研究生招生线上宣讲会直播重磅开启！
+1/6
-
Scientific Journal of Intelligent Systems Research 期刊简介
+1/5
-
双一流高校招收博士（现代技术管理/环境技术经济管理）及硕士研究生（长期有效）
+1/5
-
Scientific Journal of Humanities and Social Sciences 期刊简介
+1/4
-
福州大学吴衔誉教授课题组 诚招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
