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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » Abinit » 【求助】请问这是怎么回事，应该怎么设定参数，谢谢
61/1返回列表
查看: 562  |  回复: 5
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看
当前主题已经存档。

identation

金虫 (著名写手)

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

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

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

» 猜你喜欢
1楼 2009-12-25 20:58:37
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wenjlu

木虫 (正式写手)
★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
identation(金币+5):谢谢参与
identation(金币+10,VIP+0):input已贴，望指教 谢谢 12-26 09:25
贴你的input文件
赞一下(8人) 回复此楼
2楼2009-12-26 09:13:46
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

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

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。
赞一下(11人) 回复此楼
4楼2009-12-27 18:01:39
已阅   回复此楼   关注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问题大吗？谢谢！
第一次做磁性体系的计算，自选设定不太会呀，求指教！谢谢！
赞一下(6人) 回复此楼
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 ]
赞一下(12人) 回复此楼
6楼2009-12-28 15:19:11
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 identation 的主题更新
61/1返回列表
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
信息提示
关闭
请填处理意见
关闭