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

回复此楼

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

QE学习资料

» 猜你喜欢

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

求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

引用回帖:

Originally posted by goldenfisher at 2010-12-13 21:17:23:

自洽计算时你也可以选点。不选点的目的是为了全空间采样。能带计算时你不需要全空间的能带结构，只需要特殊方向的能带就可以了。这是因为固体的能带只在分析某些特定的高对称点取向（不一定是高对称点）有意义 ...

原来是这样…………大概了解了…………这一部分在固体物理的能带论吗？我回去再查一下，这学期代本科生的固体物理，主要也是为了在巩固一遍，暑假自学了基泰尔的《固体物理学导论》一些基础知识有时的确不太清楚…………