| 查看: 2719 | 回复: 14 | ||||
[交流]
【已解决】【求助】为何PWscf不能得到类似于CASTEP那样的收敛测试曲线呢？
|
准备转入PWscf，因为这个免费 于是例行公事一样的先对模型进行了收敛测试，使用PBE的赝势，同时使用了CASTEP在完全相同的结构、K点下进行单点能计算。结果……CASTEP下得到了很圆满的收敛测试曲线（抛物线）：  而使用PWscf，得到的居然是个直线……：  有谁也曾经碰到过这样的问题么……不确定收敛参数的PWscf的计算结果会有意义么…… [ Last edited by souledge on 2011-1-2 at 22:20 ] |
回复此楼» 收录本帖的淘帖专辑推荐
 第一原理资料汇编 |
» 猜你喜欢
296求调剂
已经有14人回复
-
生物学调剂
已经有6人回复
-
300分求调剂 （085501机械专硕，本科扬大）
已经有3人回复
-
283求调剂 086004考英二数二
已经有6人回复
-
农学0904 312求调剂
已经有5人回复
-
085600材料与化工329分求调剂
已经有10人回复
-
复试调剂
已经有9人回复
-
材料与化工调剂
已经有35人回复
-
调剂求收留
已经有10人回复
-
一志愿哈工大 085600 277 12材科基求调剂
已经有18人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 还是决定了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人回复
» 抢金币啦！回帖就可以得到:
-
2007年一起玩耍的老虫子们，想死你们啦~~~出来冒个泡吧……
+2/188
-
郑州人在北京，真诚找对象，非诚勿扰
+1/168
-
长江大学化工学院2026年硕士研究生调剂公告（理学、工学均可）
+1/89
-
大湾区大学李红庚与中国科学技术大学王柳教授联合招收博士后
+1/87
-
汕头大学硕士研究生调剂
+2/50
-
浙师大化与材学院2026年有大量硕士研究生调剂名额 （线上 调剂）
+1/41
-
中医药广东省实验生物成像团队调剂信息
+1/41
-
科研助理（RA）招聘｜适合想做芯片/架构/硬件安全方向的同学（地点：深圳理工大学）
+1/39
-
抗病毒性传染病创新药物全国重点实验室—接收化学、生物、药学、生物与医药调剂考生
+1/16
-
山西大同大学物理学专业还有调剂名额，欢迎调剂！
+1/12
-
哈尔滨工业大学（深圳）赵怡潞课题组诚招博士后、博士生
+1/12
-
齐鲁工业大学-材料方向招收两名专硕调剂
+1/11
-
大连大学-六盘水师范学院联合培养接收调剂
+1/10
-
深圳理工大学博士后招聘-先进拉曼技术课题组
+1/10
-
深圳理工大学研究助理教授招聘-先进拉曼技术课题组
+1/8
-
福建理工大学材料学院学硕专业代码08开头还有调剂名额
+1/3
-
线上复试：河南理工大学 化工学院 接收调剂硕士研究生
+1/2
-
长安大学智能爆破课题组招收2026级博士和硕士研究生
+1/2
-
中山大学课题组招实习生本硕皆可，有细胞培养，分子生物，动物实验经验优先
+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
