24小时热门版块排行榜    

查看: 2051  |  回复: 4
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

repent03

金虫 (小有名气)

[求助] Pwscf带电计算

请问该如何在Pwscf上实现带电(charged cell)计算?

是需要在input档中加上什么参数嘛?

我搜到的资讯都是VASP的

请各位虫友不吝解惑

感谢!!
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐,对您同样有帮助:

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

repent03

金虫 (小有名气)

引用回帖:
Originally posted by souledge at 2011-05-16 15:57:04:
LZ不看软件的帮助文件是不对滴……

抱歉,我才接触没多久

导师忽然要我试试quantum-espresso

您提供的资讯是加电场的计算

我想问:像一个体系如果吸付了OH负离子

总电子数应该是体系电子数+OH电子数+1

这个在VASP可透过NELECT来达到

那在Pwscf呢?

感谢!!
3楼2011-05-16 17:19:24
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 5 个回答

souledge

专家顾问 (著名写手)

【答案】应助回帖

repent03(金币+3): 感谢,我之后发问会先去看userguide的 2011-05-16 17:20:00
引用回帖:
lelfield   LOGICAL
Default: .FALSE.  
If .TRUE. a homogeneous finite electric field described
through the modern theory of the polarization is applied.
This is different from "tefield=.true." !

引用回帖:
efield REAL
Default: 0.D0  
Amplitude of the finite electric field (in Ry a.u.;
1 a.u. = 36.3609*10^10 V/m). Used only if lelfield=.TRUE.
and if k-points (K_POINTS card) are not automatic.

引用回帖:
efield_cart(i), i=1,3 REAL
Default: (0.D0, 0.D0, 0.D0)  
Finite electric field (in Ry a.u.=36.3609*10^10 V/m) in
cartesian axis. Used only if lelfield=.TRUE. and if
k-points (K_POINTS card) are automatic.

LZ不看软件的帮助文件是不对滴……
思想重于技巧,内涵重于表象
2楼2011-05-16 15:57:04
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

souledge

专家顾问 (著名写手)

【答案】应助回帖


franch(金币+1): 谢谢回帖交流 2011-05-16 20:24:16
repent03(金币+3): 感謝 2011-05-17 17:20:26
引用回帖:
Originally posted by repent03 at 2011-05-16 17:19:24:
抱歉,我才接触没多久

导师忽然要我试试quantum-espresso

您提供的资讯是加电场的计算

我想问:像一个体系如果吸付了OH负离子

总电子数应该是体系电子数+OH电子数+1

这个在VASP可透过NELECT来 ...

不好意思,刚刚没看明白问题,不过,答案还是在软件自带的文档中啊:
引用回帖:
tot_charge         REAL
Default:         0.0

    total charge of the system. Useful for simulations with charged cells.
    By default the unit cell is assumed to be neutral (tot_charge=0).
    tot_charge=+1 means one electron missing from the system,
    tot_charge=-1 means one additional electron, and so on.

    In a periodic calculation a compensating jellium background is
    inserted to remove divergences if the cell is not neutral.

思想重于技巧,内涵重于表象
4楼2011-05-16 19:27:54
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

goldenfisher

金虫 (著名写手)

【答案】应助回帖

★ ★ ★
franch(金币+3): 谢谢指点 2011-05-16 21:20:52
repent03(金币+4): 感謝 2011-05-17 17:20:43
设置tot_charge=+N或者-N, N为带电数目。
但是此处你需要注意一些问题,带点数目不能太多,有时候还需要和对应的赝势结合起来使用。
5楼2011-05-16 21:12:21
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[教师之家] 跳槽后在研项目怎么办? +5 简单化xn 2026-08-22 9/450 2026-08-23 10:59 by Tide man
[基金申请] 2026国自然函评费到账 +15 羊腰板 2026-08-21 16/800 2026-08-23 10:45 by process2012
[硕博家园] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +3 2JOx3r2CYEgw 2026-08-21 9/450 2026-08-23 10:44 by IXZuIJ2Q7OVy
[基金申请] 今天放榜吗? +15 布布和一二 2026-08-19 16/800 2026-08-23 09:55 by 张春生
[论文投稿] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +3 2JOx3r2CYEgw 2026-08-22 4/200 2026-08-23 04:04 by OEbVnUOu01ol
[硕博家园] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +3 2JOx3r2CYEgw 2026-08-22 7/350 2026-08-23 03:55 by OEbVnUOu01ol
[论文投稿] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +3 2JOx3r2CYEgw 2026-08-22 5/250 2026-08-23 03:52 by OEbVnUOu01ol
[论文投稿] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +3 QTy3jDtz1uLt 2026-08-21 9/450 2026-08-23 03:04 by OEbVnUOu01ol
[基金申请] 93BebMhtakh前后11位开头都是大写 +7 且听虎啸 2026-08-17 8/400 2026-08-22 21:55 by 医学老男孩
[基金申请] 范进中举一文的中心思想 +3 zjuhero 2026-08-22 4/200 2026-08-22 20:02 by mycsru
[基金申请] filecode,4个jtjc了 +13 ziyangfang 2026-08-19 16/800 2026-08-22 17:08 by WH3796
[基金申请] 今天基金会出结果吗?20260819 +16 kkkl_v 2026-08-19 17/850 2026-08-22 16:12 by 阿布Abu
[基金申请] 时间戳今天,20号变了 +5 archvillain 2026-08-20 5/250 2026-08-22 06:12 by hui_daxiao
[基金申请] 放榜前的不淡定 40+4 snowwithsea 2026-08-19 14/700 2026-08-21 23:51 by cratir
[基金申请] 建议基金发布提前给出明确的时间点 +10 kulium 2026-08-21 13/650 2026-08-21 21:58 by alongwaytogo
[基金申请] 今日不放榜?网传国自然预计 8 月 27 日可查结果 +16 医学老男孩 2026-08-20 20/1000 2026-08-21 21:13 by Ldrop2023
[基金申请] 看来今天不会放榜了? +8 chengyan1220 2026-08-21 11/550 2026-08-21 17:52 by dcqxinyang
[基金申请] 应该是下周三26日公布了吧? +4 哈哈蛤? 2026-08-21 4/200 2026-08-21 10:58 by Vivilian
[基金申请] 今天放榜没戏了吧 +9 yuleib84 2026-08-19 11/550 2026-08-21 10:06 by gltch
[基金申请] 重要消息,中午系统在维护 +11 yuleib84 2026-08-18 12/600 2026-08-20 11:09 by xskun
信息提示
请填处理意见