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

客户端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

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

» 猜你喜欢

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

求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的回帖

goldenfisher

金虫 (著名写手)

★
sunyang1988(金币+1):谢谢交流 2010-12-11 21:35:10

引用回帖:

Originally posted by wdlmm1214 at 2010-12-11 13:23:11:

你是说，计算反铁磁只需把同一种原子不同自旋的分成两类就可以了，也就是说如果我的原胞中有四个磁性原子，自选方向分别为 + - + -，我是分成两类，还是分成四类。

我记得这个在你提到的那个例子里面的README里面有说明。做反铁磁的时候，是要构造网络状的自旋取向的。如果不这么设置，虽然说也有可能在最后得到这个态，但是大多数时候你会得到一个不是反铁磁的亚稳态。

赞一下(1人)

7楼2010-12-11 21:21:11

已阅回复此楼关注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 个回答