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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » MS » castep 计算问题
272/3返回列表上一页123下一页
查看: 5345  |  回复: 26
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

零下1℃

铁虫 (正式写手)
★ ★ ★
franch: 金币+3, 谢谢回帖交流。。。 2014-07-29 13:08:45
引用回帖:
8楼: Originally posted by lzn_t at 2014-07-17 22:09:19
很感谢你的回复！在3楼贴出的.castep文件中，可以看到系统默认对3个E-cut值进行了计算，分别是540eV、545eV、550eV，比较后面给出的三者的Final Energy，发现其差值约在0.0003eV。这个值是否足够呢？...

我是这样理解的，精度够不够，不能但看我们的差值是小数点后的第几位，还要看你所计算的东西。。打个比方，你现在计算的是硅硅键的键能，键能远大于0.0003eV时和在0.0003eV附近时，我们对0.0003eV这样的误差理解是不一样的。所以，我们不能单看数值大小，还是要横向比较的。由于我不知道你具体计算的是什么，所以不能妄下结论，所以还是LZ仔细比对一下较好。
赞一下 回复此楼
11楼2014-07-18 08:37:09
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lzn_t

木虫 (小有名气)
引用回帖:
7楼: Originally posted by 卡开发发 at 2014-07-17 12:38:27
6L的说法有道理，并且从结果来看迭代一步的速度相当慢，大概1000s走一个scf，假设一个构型算20steps，几何优化跳一个点出来大约5小时。计算32核心的机器上采用了28核心，我不确定这个是32核心是32线程，主频如何我 ...

谢谢您的指点！计算的体系可以确定是磁性体系，但是我在结构优化时没有勾选 spin-polarized，这主要是考虑到在后面的Energy计算时，需要分别进行铁磁性和反铁磁性计算。如果在结构优化时就选择磁性，那么优化所得到的结构就会是铁磁性的（这里假定铁磁性系统能量更低）；这时再使用此铁磁性结构来计算反铁磁性，是不是不太合理呢？

[ 发自手机版 http://muchong.com/3g ]
赞一下 回复此楼
12楼2014-07-21 19:08:39
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lzn_t

木虫 (小有名气)
引用回帖:
11楼: Originally posted by 零下1℃ at 2014-07-18 08:37:09
我是这样理解的，精度够不够，不能但看我们的差值是小数点后的第几位，还要看你所计算的东西。。打个比方，你现在计算的是硅硅键的键能，键能远大于0.0003eV时和在0.0003eV附近时，我们对0.0003eV这样的误差理解是 ...

非常感谢您的指教！

[ 发自手机版 http://muchong.com/3g ]
赞一下 回复此楼
13楼2014-07-21 19:10:02
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

卡开发发

专家顾问 (著名写手)

Ab Initio Amateur
★ ★ ★
franch: 金币+3, 谢谢回帖交流。。。 2014-07-29 13:09:01
引用回帖:
12楼: Originally posted by lzn_t at 2014-07-21 19:08:39
谢谢您的指点！计算的体系可以确定是磁性体系，但是我在结构优化时没有勾选 spin-polarized，这主要是考虑到在后面的Energy计算时，需要分别进行铁磁性和反铁磁性计算。如果在结构优化时就选择磁性，那么优化所得到 ...

不是，体系只要有磁性（铁磁性反铁磁性）无论什么时候都要考虑自旋极化。优化后的结构的结果未必一样。
既然不考虑自旋极化能够scf部分计算收敛，这意味着你的体系可能存在铁磁性和无磁性，反铁磁性应该不必考虑了。
上面的问题我多嘴一句，0.0003这个差距太小，很难比较，这个0.003是为了体结构优化对基函数作调整才在不同的截断能算三次，如果固定晶格优化则不会出现这样。一般建议50~150这个区间做测试。
赞一下 回复此楼
不一定挂在论坛，计算问题问题欢迎留言。
14楼2014-07-21 19:16:28
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lzn_t

木虫 (小有名气)
引用回帖:
14楼: Originally posted by 卡开发发 at 2014-07-21 19:16:28
不是，体系只要有磁性（铁磁性反铁磁性）无论什么时候都要考虑自旋极化。优化后的结构的结果未必一样。
既然不考虑自旋极化能够scf部分计算收敛，这意味着你的体系可能存在铁磁性和无磁性，反铁磁性应该不必考虑了 ...

谢谢！在最早的计算时，我并没有考虑反铁磁性，并且从结构优化到能量计算都是考虑了自旋极化；从计算结果来看，体系也全部是铁磁性的。但是在审稿人的意见中，除了计算精度外，还建议对反铁磁性进行分析。您这里说不必考虑反铁磁性，能说一下原因吗？
从这几天的收敛性测试来看，E_cut在500eV时，Final energy 与收敛值（600eV时的值）的差小于0.0002eV；k点为4x4x4时，Final energy与收敛值（5x5x5时的值）的差小于0.00002eV，那么这样的设置（500eV，4x4x4）是否足够呢？

[ 发自手机版 http://muchong.com/3g ]
赞一下 回复此楼
15楼2014-07-21 19:59:59
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

卡开发发

专家顾问 (著名写手)

Ab Initio Amateur
引用回帖:
15楼: Originally posted by lzn_t at 2014-07-21 19:59:59
谢谢！在最早的计算时，我并没有考虑反铁磁性，并且从结构优化到能量计算都是考虑了自旋极化；从计算结果来看，体系也全部是铁磁性的。但是在审稿人的意见中，除了计算精度外，还建议对反铁磁性进行分析。您这里说 ...

第一个问题，如果无磁性计算能够成功我觉得有很大可能性会被优化到无磁性的情况，当然这点我不是完全确定（所以我用了应该这个词），有条件的话可以验证一下才会稳妥。

第二个问题，大体来说是可以了，如果没采用展宽并且空带数目充裕。如果采用了展宽，那么采用外推至0K的能量，但是这样计算能量和态的时候仍然要去掉展宽。另外建议再比较600+的趋势及k的趋势，有时候可能有跳点。
赞一下(1人) 回复此楼
不一定挂在论坛，计算问题问题欢迎留言。
16楼2014-07-21 20:24:25
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

零下1℃

铁虫 (正式写手)
引用回帖:
13楼: Originally posted by lzn_t at 2014-07-21 19:10:02
非常感谢您的指教！
...

客气了，相互学习。
回复此楼
17楼2014-07-22 08:38:01
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lzn_t

木虫 (小有名气)
引用回帖:
16楼: Originally posted by 卡开发发 at 2014-07-21 20:24:25
第一个问题，如果无磁性计算能够成功我觉得有很大可能性会被优化到无磁性的情况，当然这点我不是完全确定（所以我用了应该这个词），有条件的话可以验证一下才会稳妥。

第二个问题，大体来说是可以了，如果没采 ...

很感谢您的帮助！
从我后来的计算来看，对于同一个体系，无论是无磁性还是磁性计算，都是可以优化成功的。而对这个体系的磁性计算，无论是铁磁性还是反铁磁性，也都可以得到SCF收敛的结果。从我计算的模型来说，铁磁性体系能量是最低的，其次为反铁磁性，能量最高的是无磁性计算得到的结果。

还有一点需要请教一下：在审稿人的意见中提出“It has to be said whether the semicore sp states of Mn pseudopotential were included in the valence band.”我在计算时使用的是系统的默认自带赝势，在计算的castep文件中有这样的描述“Pseudo atomic calculation performed for Mn 3d5 4s2”。那么，是否可以回答：“The semicore sp states of Mn pseudopotential were not included in the valence band.”？
赞一下 回复此楼
18楼2014-08-14 16:04:44
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

卡开发发

专家顾问 (著名写手)

Ab Initio Amateur
引用回帖:
18楼: Originally posted by lzn_t at 2014-08-14 16:04:44
很感谢您的帮助！
从我后来的计算来看，对于同一个体系，无论是无磁性还是磁性计算，都是可以优化成功的。而对这个体系的磁性计算，无论是铁磁性还是反铁磁性，也都可以得到SCF收敛的结果。从我计算的模型来说，铁 ...

对，这种赝势确实没有考虑半芯态的情形。
赞一下(1人) 回复此楼
不一定挂在论坛，计算问题问题欢迎留言。
19楼2014-08-14 18:16:49
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

西北偏北1092

铜虫 (正式写手)
请问是如何判断体系收敛了呢？本人新手，希望解答！谢谢！
赞一下 回复此楼
20楼2014-08-31 09:09:00
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 lzn_t 的主题更新
272/3返回列表上一页123下一页
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 化工335求调剂 +3 摸摸猫猫头 2026-03-02 3/150 2026-03-02 20:55 by 平凡的七月
[考研] 材料学硕318求调剂 +10 February_Feb 2026-03-01 10/500 2026-03-02 20:18 by hypershenger
[考研] 0856求调剂285 +11 吕仔龙 2026-02-28 11/550 2026-03-02 20:15 by hypershenger
[考研] 0805总分292，求调剂 +10 幻想之殇 2026-03-01 10/500 2026-03-02 19:51 by zhukairuo
[考研] 261求调剂 +3 陆lh 2026-03-01 3/150 2026-03-02 19:32 by zhukairuo
[考研] 275求调剂 +7 明远求学 2026-03-01 7/350 2026-03-02 19:22 by zhukairuo
[考博] 26超级电容器申博 +3 dhdjdjend 2026-02-25 3/150 2026-03-02 19:21 by 轻松不少随
[考研] 一志愿中石油（华东）本科齐鲁工业大学 +3 石能伟 2026-03-02 3/150 2026-03-02 18:54 by caszguilin
[考研] 284求调剂 +5 天下熯 2026-03-02 5/250 2026-03-02 18:38 by caszguilin
[考研] 291 求调剂 +3 化工2026届毕业� 2026-03-02 3/150 2026-03-02 12:55 by houyaoxu
[考研] 26考研报考西工大材料308分求调剂 +4 weizhong123 2026-03-01 4/200 2026-03-02 12:46 by 无际的草原
[基金申请] 面上模板改不了页边距吧？ +6 ieewxg 2026-02-25 7/350 2026-03-02 12:44 by stidwellNK
[考研] 化工专硕342，一志愿大连理工大学，求调剂 +6 kyf化工 2026-02-28 7/350 2026-03-02 10:56 by 无际的草原
[考研] 调剂 +3 13853210211 2026-03-02 4/200 2026-03-02 10:16 by 13853210211
[考研] 材料类求调剂 +11 wana_kiko 2026-02-28 14/700 2026-03-02 08:46 by 聪明的大松鼠
[考研] 328求调剂 +3 aaadim 2026-03-01 5/250 2026-03-01 17:29 by njzyff
[考研] 307求调剂 +5 wyyyqx 2026-03-01 5/250 2026-03-01 15:21 by Fff-1
[考研] 调剂 +3 简木ChuFront 2026-02-28 3/150 2026-03-01 11:46 by 王伟要上岸啊
[考研] 寻找调剂 +4 LYidhsjabdj 2026-02-28 4/200 2026-03-01 10:56 by sunny81
[高分子] 求环氧树脂研发1名 +3 孙xc 2026-02-25 11/550 2026-02-28 16:57 by ichall
信息提示
关闭
请填处理意见
关闭