24小时热门版块排行榜    

查看: 2659  |  回复: 22

tiaozaoququ

木虫 (著名写手)


[交流] siesta计算结果,让我崩溃

以前一直使用VASP计算材料性质,最近学习了下siesta,发现,对于同一种材料考虑自旋计算,有时候vasp与siesta的结果完全不同,比如一个算出材料为金属,另一个得出是半导体。另外,siesta中,采用LDA与GGA算出来的结果有时候也有本质性的区别。
    想问问各位虫友,你们有没有碰到这样的问题呢?本人新手,不知道怎么取舍这些结果?大侠们给个意见吧

[ Last edited by tiaozaoququ on 2011-12-7 at 18:16 ]
回复此楼

» 猜你喜欢

» 本主题相关商家推荐: (我也要在这里推广)

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

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

查看全部散金贴

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

tiaozaoququ(金币+1):谢谢参与
一般情况下vasp给出的结果比较靠谱吧。
siesta的计算结果如vasp不容是不是赝势的原因?
2楼2011-12-07 17:23:18
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

tiaozaoququ

木虫 (著名写手)


引用回帖:
2楼: Originally posted by uuv2010 at 2011-12-07 17:23:18:
一般情况下vasp给出的结果比较靠谱吧。
siesta的计算结果如vasp不容是不是赝势的原因?

如果不计自旋计算,VASP与siesta的结果大部分还是比较一致,但若考虑自旋,就出现不一致的问题了。赝势选取的不同,结果会有不同。
3楼2011-12-07 18:01:00
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

vv_china

铜虫 (小有名气)


★ ★
tiaozaoququ(金币+1):谢谢参与
uuv2010(金币+1): 欢迎参加讨论 2011-12-12 22:17:26
呵呵 有时候同一个计算软件,版本不同,计算结果都会有较大差异。
4楼2011-12-08 12:54:43
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

jixia1981

新虫 (初入文坛)


★ ★
tiaozaoququ(金币+1):谢谢参与
uuv2010(金币+1): 欢迎参加讨论 2011-12-12 22:17:41
siesta计算,基组的影响也可能会比较大,换一个精确的基组试试
5楼2011-12-08 21:00:28
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

tiaozaoququ

木虫 (著名写手)


引用回帖:
5楼: Originally posted by jixia1981 at 2011-12-08 21:00:28:
siesta计算,基组的影响也可能会比较大,换一个精确的基组试试

我现在怀疑我选用的赝势文件有问题
由于自己不会产生赝势,所以计算中所用的赝势都是从siesta官网上下的。不知道你们使用的是不是自己产生赝势?
6楼2011-12-09 13:09:22
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

tiaozaoququ

木虫 (著名写手)


引用回帖:
4楼: Originally posted by vv_china at 2011-12-08 12:54:43:
呵呵 有时候同一个计算软件,版本不同,计算结果都会有较大差异。

我全选用的是DZP基组
7楼2011-12-09 13:10:07
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

jixia1981

新虫 (初入文坛)


你使用的是那些元素的赝势?化合物的名称?
早期siesta是要自己产生赝势,比较费时。现在有了赝势可以下载,但是使用之前,要针对自己的体系进行测试,一个赝势的平移性毕竟是有限的。
8楼2011-12-09 13:51:29
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
呵呵~~~我也突然之间才发现~~~~~
楼主加油啊!
13楼2011-12-09 18:29:03
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

tiaozaoququ

木虫 (著名写手)


引用回帖:
13楼: Originally posted by imdiablo1 at 2011-12-09 18:29:03:
呵呵~~~我也突然之间才发现~~~~~
楼主加油啊!

我也当你听众算了,嘿嘿
14楼2011-12-09 20:02:17
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

tiaozaoququ

木虫 (著名写手)


引用回帖:
8楼: Originally posted by jixia1981 at 2011-12-09 13:51:29:
你使用的是那些元素的赝势?化合物的名称?
早期siesta是要自己产生赝势,比较费时。现在有了赝势可以下载,但是使用之前,要针对自己的体系进行测试,一个赝势的平移性毕竟是有限的。

如果测试发现赝势不好,就需要自己去产生了吗?关键是那个atom模块的介绍材料有点单薄呀,看完还是丈二摸不着头脑Ing。
15楼2011-12-09 20:15:18
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
★ ★
tiaozaoququ(金币+1):谢谢参与
uuv2010(金币+1): 多谢指导! 2011-12-16 09:20:55
uuv2010: 2011-12-16 09:21:00
如果你用LDA赝势的CA-PZ方案,那么你不必给出XC.functional这个选项,因为系统默认为LDA;但是,如果你用GGA的赝势,假若你没指定XC.functional为GGA,那么计算仍会采用默认的LDA赝势计算而不报错,DOS曲线就大不一样了,如果你用GGA赝势,你要指定XC.functional为GGA,同时指定XC.author为PBE,PW等等 参照手册 你看看是不是这个原因
20楼2011-12-12 20:34:15
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

tiaozaoququ

木虫 (著名写手)


引用回帖:
20楼: Originally posted by identation at 2011-12-12 20:34:15:
如果你用LDA赝势的CA-PZ方案,那么你不必给出XC.functional这个选项,因为系统默认为LDA;但是,如果你用GGA的赝势,假若你没指定XC.functional为GGA,那么计算仍会采用默认的LDA赝势计算而不报错,DOS曲线就大不 ...

不是这个原因啦
我发现Lda 与gga结果发生不一致时,都是在计及自旋时。
21楼2011-12-12 21:45:06
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

fengshiquan

铁杆木虫 (正式写手)


★ ★
tiaozaoququ(金币+1):谢谢参与
uuv2010(金币+1): 欢迎参加讨论 2011-12-16 09:21:20
siesta结果没有vasp的可信,毕竟vasp是商业性软件,里面的赝势都是调试过的
22楼2011-12-16 09:12:16
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

tiaozaoququ

木虫 (著名写手)


引用回帖:
22楼: Originally posted by fengshiquan at 2011-12-16 09:12:16:
siesta结果没有vasp的可信,毕竟vasp是商业性软件,里面的赝势都是调试过的

应该是的
可能用ATK算出的结果会比siesta可信
23楼2011-12-16 12:13:15
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
简单回复
2011-12-09 14:14   回复  
tiaozaoququ(金币+1):谢谢参与
hjduandian10楼
2011-12-09 14:41   回复  
tiaozaoququ(金币+1):谢谢参与
shotao11楼
2011-12-09 15:00   回复  
tiaozaoququ(金币+1):谢谢参与
imdiablo112楼
2011-12-09 15:22   回复  
tiaozaoququ(金币+1):谢谢参与
tiaozaoququ(金币+1): 你的人气好高呀,听众这么多 2011-12-09 15:45:28
panyuyu16楼
2011-12-09 20:20   回复  
tiaozaoququ(金币+1):谢谢参与
2011-12-09 20:24   回复  
tiaozaoququ(金币+1):谢谢参与
lijian197918楼
2011-12-09 20:25   回复  
tiaozaoququ(金币+1):谢谢参与
suhuishui19楼
2011-12-09 20:41   回复  
tiaozaoququ(金币+1):谢谢参与
相关版块跳转 我要订阅楼主 tiaozaoququ 的主题更新
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[论文投稿] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +3 jRl3mE6ddZGq 2026-08-03 7/350 2026-08-04 04:50 by k7OM8YghWbkC
[硕博家园] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +3 jRl3mE6ddZGq 2026-08-03 6/300 2026-08-04 04:47 by k7OM8YghWbkC
[找工作] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +3 cu9Nq1xK233Z 2026-08-03 8/400 2026-08-04 04:47 by k7OM8YghWbkC
[博后之家] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +3 cu9Nq1xK233Z 2026-08-03 8/400 2026-08-04 04:46 by k7OM8YghWbkC
[公派出国] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +3 cu9Nq1xK233Z 2026-08-03 7/350 2026-08-04 04:46 by k7OM8YghWbkC
[考博] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +3 cu9Nq1xK233Z 2026-08-03 5/250 2026-08-04 04:44 by k7OM8YghWbkC
[教师之家] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +4 cu9Nq1xK233Z 2026-08-03 8/400 2026-08-04 04:41 by k7OM8YghWbkC
[有机交流] 一个有机合成实验室都需要哪些设备? 50+3 kf2781974 2026-07-31 9/450 2026-08-04 03:56 by kf2781974
[教师之家] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +3 jRl3mE6ddZGq 2026-08-03 7/350 2026-08-04 03:35 by k7OM8YghWbkC
[考研] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +3 cu9Nq1xK233Z 2026-08-03 8/400 2026-08-04 03:21 by k7OM8YghWbkC
[硕博家园] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +4 cu9Nq1xK233Z 2026-08-03 12/600 2026-08-04 03:17 by k7OM8YghWbkC
[高分子] UV压敏胶开发 +3 ichall 2026-07-30 5/250 2026-08-03 14:40 by Sunrisepay
[基金申请] 面上提前没消息,有中的吗 +14 archvillain 2026-08-02 16/800 2026-08-03 12:23 by fuweiguochen
[基金申请] 面上再次挂了,太难了,躺也躺不了,倦也卷不过,小学校之殇! +19 低垂的野花 2026-07-31 25/1250 2026-08-03 09:25 by gy116024
[基金申请] 2026年国自然面上资助率 +16 布布和一二 2026-07-30 22/1100 2026-08-03 09:20 by gy116024
[教师之家] 基础研究怎么拉横向,学校到款任务越来越多,难以完成 拉横向,都有哪些途径啊 +8 锦衣卫寒战 2026-07-28 12/600 2026-08-03 08:22 by syl200707
[基金申请] 微信指数没变化,科研之友没阅读 +15 wangze12014 2026-07-28 19/950 2026-08-02 20:05 by 蔡棒棒菂
[考博] 2027年申博 50+3 射雕英雄胜 2026-07-30 3/150 2026-08-02 09:36 by lfy8008
[高分子] HXDI做水性聚氨酯乳液,是不是特别容易出渣 15+3 yuyusuv 2026-07-29 3/150 2026-07-31 09:21 by huizingga
[基金申请] 系统今天又提示维护了,估计离放榜不远了 +11 winnerche 2026-07-29 15/750 2026-07-30 23:08 by jnhyjjm
信息提示
请填处理意见