24小时热门版块排行榜    

查看: 402  |  回复: 5
当前主题已经存档。
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

identation

金虫 (著名写手)

[交流] 【求助】请问这是怎么回事,应该怎么设定参数,谢谢

铁磁性体系 设定自旋之后 收敛测试过程中出现这个warning,请问是怎么回事?求高手指教,感激不尽!

[ Last edited by identation on 2009-12-27 at 19:46 ]
回复此楼

» 猜你喜欢

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

xbelinda

银虫 (小有名气)

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
identation(金币+5):谢谢参与
identation(金币+10,VIP+0):谢谢,体系是半导体,不是金属;不是金属不是occopt就要选1吗,谢谢 12-27 19:43
ice_rain(金币+1,VIP+0):谢谢交流 12-27 20:06
是不是你的体系是金属,你用occopt=7,然后测试一下 tsmear。
你有两个ecut。
去掉diemac。
4楼2009-12-27 18:01:39
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 6 个回答

identation

金虫 (著名写手)

引用回帖:
Originally posted by wenjlu at 2009-12-26 09:13:
贴你的input文件

nsppol 2
nspden 2
prtden 1
spinat
0.0 0.0 0.0
0.0 0.0 0.0
0.0 0.0 0.0
0.0 0.0 0.0
0.0 2.0 0.0
0.0 2.0 0.0
0.0 2.0 0.0
0.0 2.0 0.0
0.0 0.0 0.0
0.0 0.0 0.0
0.0 0.0 0.0
0.0 0.0 0.0
0.0 0.0 0.0
0.0 0.0 0.0
0.0 0.0 0.0
0.0 0.0 0.0
0.0 0.0 0.0
0.0 0.0 0.0
0.0 0.0 0.0
0.0 0.0 0.0
ndtset 5
ecut: 20
ecut+ 3
nband  75            
kptopt 1 # Kpoint option ( default 0.0 )
ecut 20
ngkpt 10 10 6  # Number of k-points
nshiftk 1
shiftk 0.5 0.5 0.5
toldfe 1.0d-6 # Tolerance on the difference of total energy
iscf 5 # SCF cycle, CG based on the minim. of the energy
nstep 100 # Number of self-consistent field steps
diemac 12.0 # Macroscopic dielectric constant diemac of the system
# is a useful help to speed-up the SCF procedure
enunit 2 # Energy unit of output
ixc 1 # Exchange Correlation choice
znucl (此处略去)
ntypat 4
occopt 1
natom 20
typat 4*1 2*2 2*3 12*4
tsmear 0.01
acell (此处略去)
angdeg (此处略去)
rprim
1.0 0.0 0.0
0.0 1.0 0.0
0.0 0.0 1.0
xred
(此处略去)

谢谢!
3楼2009-12-26 09:24:46
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

identation

金虫 (著名写手)


TIGERYZZ(金币+1,VIP+0):谢谢回帖交流 12-27 21:23
引用回帖:
Originally posted by xbelinda at 2009-12-27 18:01:
是不是你的体系是金属,你用occopt=7,然后测试一下 tsmear。
你有两个ecut。
去掉diemac。

谢谢!体系不是金属 是半导体,半导体occopt=1,那其他的怎么设定呢?这个warning问题大吗?谢谢!
第一次做磁性体系的计算,自选设定不太会呀,求指教!谢谢!
5楼2009-12-27 19:44:41
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wenjlu

木虫 (正式写手)

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
identation(金币+10,VIP+0): 12-28 17:26
aylayl08(金币+1,VIP+0):谢谢解答 12-28 20:06
我猜测你做的是113双钙钛矿结构的计算。

出错,主要问题在你设定的spinat不对。
在做线性的磁性计算时候,只允许z方向有磁矩。其它方向为零。 即:
spinat=0.0   0. 0   2.0
关于非线性情况,就需要其它关键词来设定完成。

[ Last edited by wenjlu on 2009-12-28 at 15:20 ]
6楼2009-12-28 15:19:11
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通表情 高级回复 (可上传附件)
信息提示
请填处理意见