版块导航: 正在加载中...

客户端APP下载

登录注册

应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜

当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

wdlmm1214

新虫 (初入文坛)

应助: 0 (幼儿园)
金币: 195.8
帖子: 8
在线: 5.9小时
虫号: 1076207

[交流] 【求助】求助pwscf一些参数的物理意义和设定经验

Sample Text
本人刚学pwscf，有很多参数不知道什么意思看说明书没看懂，希望高手指点。
我先通过例子去问问题。# self-consistent calculation
cat > feo.scf.in << EOF
&control
pseudo_dir = '$PSEUDO_DIR',
outdir='$TMP_DIR'
prefix='feo',
/
&system
ibrav=  0, celldm(1)=8.19, nat=  4, ntyp= 3,
ecutwfc = 30.0, ecutrho = 240.0,
occupations='smearing', smearing='mp', degauss=0.02,
nspin=2,
starting_magnetization(2)= 0.5, starting_magnetization(3)=-0.5,
lda_plus_u=.true.
Hubbard_U(2)=4.3, Hubbard_U(3)=4.3,
starting_ns_eigenvalue(3,2,2)=1.d0, starting_ns_eigenvalue(3,1,3)=1.d0
/
&electrons
mixing_beta = 0.3
/
CELL_PARAMETERS
0.50 0.50 1.00
0.50 1.00 0.50
1.00 0.50 0.50
ATOMIC_SPECIES
O 1.  O.LDA.US.RRKJ3.UPF
Fe1  1.  Fe.pz-nd-rrkjus.UPF
Fe2  1.  Fe.pz-nd-rrkjus.UPF
ATOMIC_POSITIONS crystal
O  0.25 0.25 0.25
O  0.75 0.75 0.75
Fe1 0.0  0.0  0.0
Fe2 0.5  0.5  0.5
K_POINTS automatic
4 4 4 0 0 0
这是Feo的例子，相信很多人都看过了。
1.FeO中明明只有Fe,O,Fe也应该全是二价，两种原子为什么ntyp= 3，2个Fe还用了两个一样的赝势，为什么不把两个氧原子看成不同类原子，这个类是怎么分的，VASP中算磁性，是对每个原子分别设定磁距很方便，但在这里说明书里说针对每一类原子去设定，我有点不明白，是不是为了表现反铁磁的磁结构，把不同自旋的原子看成不同类啊？那要是螺旋磁结构又该如何呢，希望高手帮忙解答。
2.是针对LDA+U的问题有两个参数一个是starting_ns_eigenvalue(3,2,2)=1.d0，我怎么也没看懂这个参数的意思，还有一个是和Hubbard_alpha(I)，在VASP利用U和J，来计算LDA+U，但这个Hubbard_alpha，我不明白啥意思，怎么去设定，这个最关键，希望高手能给你一些这个参数的设定经验。
3.mixing_beta 这个参数，在软件包的例子中也总变，希望高手能告诉一些设定的经验。
4.例子中有些KPOINTS是手动设定的例如
K_POINTS
28
0.0 0.0 0.0 1.0
0.0 0.0 0.1 1.0
0.0 0.0 0.2 1.0
0.0 0.0 0.3 1.0
0.0 0.0 0.4 1.0
0.0 0.0 0.5 1.0
0.0 0.0 0.6 1.0
0.0 0.0 0.7 1.0
0.0 0.0 0.8 1.0
0.0 0.0 0.9 1.0
0.0 0.0 1.0 1.0
0.0 0.0 0.0 1.0
0.0 0.1 0.1 1.0
0.0 0.2 0.2 1.0
0.0 0.3 0.3 1.0
0.0 0.4 0.4 1.0
0.0 0.5 0.5 1.0
0.0 0.6 0.6 1.0
0.0 0.7 0.7 1.0
0.0 0.8 0.8 1.0
0.0 0.9 0.9 1.0
0.0 1.0 1.0 1.0
0.0 0.0 0.0 1.0
0.1 0.1 0.1 1.0
0.2 0.2 0.2 1.0
0.3 0.3 0.3 1.0
0.4 0.4 0.4 1.0
0.5 0.5 0.5 1.0
我想知道怎么手动设定K点，尤其后面的权重，怎么设定。Sample Text

» 收录本帖的淘帖专辑推荐

» 猜你喜欢

面上没中，邀请各位路过的虫友分析一下分数已经有13人回复
售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急已经有8人回复
两块石头已经有8人回复
初秋的晨风已经有6人回复
学科评审组评审是指会评吗？已经有9人回复
售SCI一区T0P文章，我:8O.55.1.O.54,科目全,可伽急已经有14人回复
广西大学-广州大学招聘博士后欢迎广大优秀人才！！！已经有7人回复
科研人应该花精力去思考如何解决问题，而不是去凝练问题已经有17人回复
Ei源刊怎么投已经有4人回复
哈尔滨工业大学韩晓军教授课题组招收2027年硕士推免生及博士研究生已经有3人回复

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

求PWSCF使用说明书已经有3人回复
[求助]PWscf安装已经有11人回复
【已解决】【求助】为何PWscf不能得到类似于CASTEP那样的收敛测试曲线呢？已经有14人回复
【求助】pwscf错误提示diagonalization(ZHEGV*) falied 已经有8人回复
【求助】利用abinit或pwscf软件在单机上可以计算半导体材料的哪些性质? 已经有8人回复
【求助】关于PWSCF介电常数和声子谱的计算已经有7人回复
【求助】PWSCF计算介电极化时为何出现如下提示？已经有19人回复
【求助】用PWSCF中的LDA+U做反铁磁绝缘体系的弛豫计算已经有23人回复
【求助】请教PWSCF晶体建模和自洽计算，不收敛已经有21人回复

» 抢金币啦！回帖就可以得到:

查看全部散金贴

1楼 2010-12-10 23:03:37

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

改王之王

金虫 (正式写手)

1ST强帖: 1
应助: 0 (幼儿园)
金币: 1120.7
帖子: 817
在线: 113.5小时
虫号: 968169

★
zhang668(金币+1):多谢交流 2010-12-30 20:23:07

我也有一点不太懂…………就是计算能带一般选高对称的K点…………但是计算自洽时自己不自动选K点时，按什么规律选取呢？看了看例子好像和算能带取的不一样……

赞一下(1人)

9楼2010-12-12 08:37:30

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

查看全部 14 个回答

锐利的碎片

木虫 (正式写手)

1ST强帖: 1
应助: 136 (高中生)
金币: 3637.1
帖子: 988
在线: 1224.9小时
虫号: 961933

★
ym23(金币+1):谢谢帮助 2010-12-10 23:56:05
wdlmm1214(金币+5): 2010-12-11 02:24:26

alpha不用管，那个是用来算U值的。mixing beta一般默认，如果计算不收敛可以适当改一下。

赞一下(2人)

2楼2010-12-10 23:27:19

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

hbhehe

木虫 (著名写手)

应助: 0 (幼儿园)
金币: 2626
帖子: 1827
在线: 136.5小时
虫号: 1120429

★
wdlmm1214(金币+1):谢谢参与

同问，帮顶一下...

3楼2010-12-11 09:20:31

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

暖秋

新虫 (小有名气)

应助: 0 (幼儿园)
金币: 12.5
帖子: 70
在线: 36小时
虫号: 797547

★ ★
wdlmm1214(金币+1):谢谢参与
wdlmm1214(金币+5): 2010-12-11 13:24:01
sunyang1988(金币+1):谢谢交流 2010-12-11 13:58:24

计算磁序的时候，确实是把自旋不同的原子看成不同类的原子，至于你说的螺旋结构，我不清楚，也没做过。。。抱歉

赞一下(3人)

4楼2010-12-11 10:21:06

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

简单回复

liusha339910楼

2010-12-12 09:09 回复

祝福祝福祝福祝福祝福祝福祝福祝福

查看全部 14 个回答