| 查看: 1196 | 回复: 5 | |||
enola捐助贵宾 (正式写手)
|
[交流]
【分享】pwscf模拟水分子 已有5人参与
|
|
看见有人在MS版面上提到水分子的模拟，今天有空自己做了下，可是出来结果不对（输入时，特意把H－O－H键角是105 deg,但pw优化后成了97.5deg ）不知道是不是水中的H得考虑自旋阿,还是其他有忽略的地方，输入文件如下，请大家批评阿： &CONTROL calculation = 'relax', prefix = 'H2O', pseudo_dir = '/home/enola/pwscf/espresso-3.2.3/pseudo', outdir = '/home/enola/tmp', tprnfor=.ture forc_conv_thr=1.0D-4 nstep=100 / &SYSTEM ibrav = 0, nat = 3, ntyp = 2, ecutwfc = 35.D0, ecutrho = 150.D0, / &ELECTRONS conv_thr = 1.D-7, mixing_beta = 0.7D0, / &IONS ion_dynamics='bfgs', pot_extrapolation = 'second_order', wfc_extrapolation = 'second_order', / CELL_PARAMETERS cubic 20.0 0.0 0.0 0.0 20.0 0.0 0.0 0.0 20.0 ATOMIC_SPECIES H 1.00d0 H.fpmd.UPF O 16.0d0 O.BLYP.UPF ATOMIC_POSITIONS (bohr) O 0.0099 0.0099 0.0000 H 1.4325 -0.2243 -0.0001 H -0.2243 1.8325 0.0002 K_POINTS {Gamma} |
回复此楼» 猜你喜欢
2026年机械制造与材料应用国际会议 （ICMMMA 2026）
已经有4人回复
-
磺酰氟产物，毕不了业了！
已经有6人回复
-
求助:我三月中下旬出站，青基依托单位怎么办？
已经有9人回复
-
Cas 72-43-5需要30g，定制合成，能接单的留言
已经有8人回复
-
北京211副教授，35岁，想重新出发，去国外做博后，怎么样？
已经有8人回复
-
论文终于录用啦！满足毕业条件了
已经有25人回复
-
自荐读博
已经有3人回复
-
不自信的我
已经有5人回复
-
投稿Elsevier的杂志（返修），总是在选择OA和subscription界面被踢皮球
已经有8人回复
» 本主题相关商家推荐: (我也要在这里推广)
» 本主题相关价值贴推荐，对您同样有帮助:
- PWSCF中如何填写基本的参数
已经有9人回复
- PWscf计算homo lumo或者valence conduction band，问题已解决
已经有14人回复
- 【已解决】【求助】为何PWscf不能得到类似于CASTEP那样的收敛测试曲线呢？
已经有14人回复
- 【求助】pwscf错误提示diagonalization(ZHEGV*) falied
已经有8人回复
- 【求助】PWSCF计算介电极化时为何出现如下提示？
已经有19人回复
- 【求助】请教PWSCF晶体建模和自洽计算，不收敛
已经有21人回复
2楼2010-06-05 15:39:32
y1ding
铁杆木虫 (著名写手)
- 1ST强帖: 1
- 应助: 61 (初中生)
- 贵宾: 0.33
- 金币: 5959.3
- 散金: 1
- 红花: 21
- 帖子: 1884
- 在线: 491.1小时
- 虫号: 142265
- 注册: 2005-12-21
- 专业: 凝聚态物性 II ：电子结构
3楼2010-06-05 18:19:42
goldenfisher
金虫 (著名写手)
- 1ST强帖: 3
- 应助: 55 (初中生)
- 贵宾: 0.153
- 金币: 5779.1
- 散金: 2908
- 红花: 21
- 帖子: 1552
- 在线: 351.5小时
- 虫号: 117764
- 注册: 2005-11-26
- 性别: GG
- 专业: 原子和分子物理
★ ★ ★
小木虫(金币+0.5):给个红包，谢谢回帖交流
enola(金币+2): 2010-06-05 19:55:33
aylayl08(金币+2):谢谢指点 2010-06-05 20:28:50
小木虫(金币+0.5):给个红包，谢谢回帖交流
enola(金币+2): 2010-06-05 19:55:33
aylayl08(金币+2):谢谢指点 2010-06-05 20:28:50
|
赝势用的是模守恒吧，我想ecutrho = 150.D0应该差不多了。但是ecutwfc太小，有H在这里，你必须用大一些。建议至少60Ry. 此外，我想很可能是初始构型的问题。我对水分子进行过无数次优化，都是对的。 你可以考虑用别的赝势试一下。 版本太低了啊。呵呵 我总觉得整个输入文件缺点什么。 [ Last edited by goldenfisher on 2010-6-5 at 19:49 ] |
4楼2010-06-05 19:47:24
enola
捐助贵宾 (正式写手)
- 1ST强帖: 2
- 应助: 13 (小学生)
- 金币: 1485.9
- 散金: 1549
- 红花: 14
- 帖子: 847
- 在线: 102.5小时
- 虫号: 947564
- 注册: 2010-01-24
- 专业: 半导体光电子器件
5楼2010-06-06 10:10:57
6楼2010-06-10 23:07:25