| 查看: 1444 | 回复: 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} |
回复此楼» 猜你喜欢
可靠消息，8月21号(明天)放榜
已经有3人回复
-
只有每年这种时候来逛逛小木虫
已经有20人回复
-
欢迎发来filecode的Mz6后的代码验证其规律
已经有125人回复
-
今天基金会出结果吗？20260819
已经有15人回复
-
时间戳变了，能看出什么问题？
已经有16人回复
-
filecode，4个jtjc了
已经有13人回复
-
为什么余额宝的年化利率越来越低？主要原因有哪些？
已经有5人回复
-
今天放榜没戏了吧
已经有5人回复
-
今天放榜吗？
已经有15人回复
-
filecode=后面第一个是大写字母
已经有12人回复
» 本主题相关商家推荐: (我也要在这里推广)
» 本主题相关价值贴推荐，对您同样有帮助:
- 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
