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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » QE(Pwscf) » 【求助】求助pwscf一些参数的物理意义和设定经验
141/1返回列表
查看: 3579  |  回复: 13
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

wdlmm1214

新虫 (初入文坛)

[交流] 【求助】求助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学习资料

» 猜你喜欢

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

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

查看全部散金贴
1楼 2010-12-10 23:03:37
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
回帖支持 ( 显示支持度最高的前 50 名 )

锐利的碎片

木虫 (正式写手)
★
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的回帖

暖秋

新虫 (小有名气)
★ ★
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的回帖

暖秋

新虫 (小有名气)
★
sunyang1988(金币+1):谢谢交流 2010-12-11 19:24:37
wdlmm1214(金币+5): 2010-12-13 00:03:46
引用回帖:
Originally posted by wdlmm1214 at 2010-12-11 13:23:11:

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

分成两类，这两类的磁矩设置为相反，就OK了。。。
赞一下(2人) 回复此楼
6楼2010-12-11 18:11:04
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

whw19850730

铁杆木虫 (小有名气)
★
zhang668(金币+1):多谢交流 2010-12-23 15:28:21
[quote]Originally posted by wdlmm1214 at 2010-12-11 13:23:11:

starting_ns_eigenvalue(3,2,2)=1.d0
这个是调整3d轨道占据数的，因为不调整，自动优化的占据数可能不是基态。
赞一下(2人) 回复此楼
13楼2010-12-23 11:52:35
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通回帖

hbhehe

木虫 (著名写手)
★
wdlmm1214(金币+1):谢谢参与
同问，帮顶一下...
回复此楼
3楼2010-12-11 09:20:31
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wdlmm1214

新虫 (初入文坛)
引用回帖:
Originally posted by 暖秋 at 2010-12-11 10:21:06:
计算磁序的时候，确实是把自旋不同的原子看成不同类的原子，至于你说的螺旋结构，我不清楚，也没做过。。。抱歉

你是说，计算反铁磁只需把同一种原子不同自旋的分成两类就可以了，也就是说如果我的原胞中有四个磁性原子，自选方向分别为 +  -  +  -，我是分成两类，还是分成四类。
赞一下 回复此楼
5楼2010-12-11 13:23:11
已阅   回复此楼   关注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的回帖

goldenfisher

金虫 (著名写手)
★
sunyang1988(金币+1):谢谢交流 2010-12-11 21:35:24
说实话，你要问的问题真的很多，而且都太基础了，都是可以直接在帮助文件里找到的。你不懂是因为你没有去理解它们的物理意义。等你看了相关文献你就明白了。
赞一下(1人) 回复此楼
8楼2010-12-11 21:24:34
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

改王之王

金虫 (正式写手)
★
zhang668(金币+1):多谢交流 2010-12-30 20:23:07
我也有一点不太懂…………就是计算能带一般选高对称的K点…………但是计算自洽时自己不自动选K点时，按什么规律选取呢？看了看例子好像和算能带取的不一样……
赞一下(1人) 回复此楼
9楼2010-12-12 08:37:30
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

goldenfisher

金虫 (著名写手)
★ ★ ★
zzy870720z(金币+3):谢谢指教 2010-12-13 23:00:08
引用回帖:
Originally posted by 改王之王 at 2010-12-12 08:37:30:
我也有一点不太懂…………就是计算能带一般选高对称的K点…………但是计算自洽时自己不自动选K点时，按什么规律选取呢？看了看例子好像和算能带取的不一样……

自洽计算时你也可以选点。不选点的目的是为了全空间采样。能带计算时你不需要全空间的能带结构，只需要特殊方向的能带就可以了。这是因为固体的能带只在分析某些特定的高对称点取向（不一定是高对称点）有意义，其他方向就可能重复了。
这些道理应该在固体物理书上讲得很清楚啊
赞一下(1人) 回复此楼
11楼2010-12-13 21:17:23
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

改王之王

金虫 (正式写手)
引用回帖:
Originally posted by goldenfisher at 2010-12-13 21:17:23:


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

原来是这样…………大概了解了…………这一部分在固体物理的能带论吗？我回去再查一下，这学期代本科生的固体物理，主要也是为了在巩固一遍，暑假自学了基泰尔的《固体物理学导论》一些基础知识有时的确不太清楚…………
赞一下 回复此楼
12楼2010-12-13 21:42:42
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

菲顶顶

铜虫 (正式写手)
★
小木虫: 金币+0.5, 给个红包，谢谢回帖
楼主知道K点的权重怎么弄了吗，
赞一下 回复此楼
14楼2016-03-31 21:09:58
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
简单回复
liusha339910楼
2010-12-12 09:09   回复  
 祝福祝福祝福祝福祝福祝福祝福祝福
相关版块跳转 我要订阅楼主 wdlmm1214 的主题更新
141/1返回列表
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考博] 【2027博士申请】纳米药物递送方向 20+3 13586093586 2026-08-03 5/250 2026-08-08 21:54 by 13586093586
[基金申请] 好奇怪的filecode +3 布布和一二 2026-08-08 4/200 2026-08-08 17:46 by zhanghaozhu
[基金申请] 售SCI一区T0P文章，我:8.O.55.1.O.54,科目齐全,可+急 +3 qTvzQBRAHjCi 2026-08-07 5/250 2026-08-08 17:02 by oEVWOejN9taj
[基金申请] 售SCI文章，我:8O.5.5.1O.54,科目全,可十急 +3 Vi50GxzrFcSG 2026-08-07 5/250 2026-08-08 16:27 by oEVWOejN9taj
[公派出国] 售SCI一区T0P文章，我:8O.55.1.O.5.4,科目齐全,可+急 +3 Vi50GxzrFcSG 2026-08-07 5/250 2026-08-08 15:59 by oEVWOejN9taj
[论文投稿] 售SCI一区T0P文章，我:8.O.55.1.O.54,科目齐全,可+急 +3 WQyTGMbfH7Cx 2026-08-07 4/200 2026-08-08 15:02 by oEVWOejN9taj
[博后之家] 售SCI一区T0P文章，我:8O.55.1.O.54,科目全,可伽急 +3 HEQlVqMTIA7d 2026-08-07 6/300 2026-08-08 14:42 by oEVWOejN9taj
[教师之家] 售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急 +4 HEQlVqMTIA7d 2026-08-07 5/250 2026-08-08 14:27 by oEVWOejN9taj
[论文投稿] 售SCI一区T0P文章，我:8O.55.1.O.5.4,科目齐全,可+急 +3 HEQlVqMTIA7d 2026-08-07 5/250 2026-08-08 14:22 by oEVWOejN9taj
[论文投稿] 售SCI文章，我:8O.5.5.1O.54,科目全,可十急 +3 KXLV3nuBVBY7 2026-08-07 6/300 2026-08-08 14:01 by oEVWOejN9taj
[硕博家园] 售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急 +3 Vi50GxzrFcSG 2026-08-07 6/300 2026-08-08 10:19 by 3OOjAIS77qg2
[考博] 售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急 +3 HEQlVqMTIA7d 2026-08-07 3/150 2026-08-08 03:47 by 6vVgjDL4CnGu
[基金申请] 娱乐 +6 Tide man 2026-08-03 6/300 2026-08-07 22:40 by 铁帽子农民
[基金申请] 化学口download_prp&amp;fileCode的固定段好像这几天一直没变，有变的大神么？ +3 Tide man 2026-08-07 4/200 2026-08-07 22:39 by Tide man
[基金申请] 听说今天filecode变了 +24 布布和一二 2026-08-06 47/2350 2026-08-07 16:02 by zhiyanjiang
[基金申请] filecode变化情况 +6 布布和一二 2026-08-07 22/1100 2026-08-07 14:45 by 且听虎啸
[基金申请] 8月时间戳变的，举个手。玩一下，释放压力 +9 archvillain 2026-08-04 11/550 2026-08-05 20:06 by wlwhappy
[基金申请] 好消息？这个有何含义？？？ +8 Tide man 2026-08-05 10/500 2026-08-05 16:14 by xmuxiaoyu
[基金申请] 纯娱乐，不喜欢勿喷 +7 Tide man 2026-08-04 10/500 2026-08-04 15:10 by loufangrui
[基金申请] 什么时候能放榜呀？ +3 Jacob678 2026-08-03 3/150 2026-08-03 16:14 by gltch
信息提示
关闭
请填处理意见
关闭