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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » QE(Pwscf) » 收敛测试，这种情况我该继续测试么,简直就是跳跃式的嘛
252/3返回列表上一页123下一页
查看: 5117  |  回复: 24
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

souledge

专家顾问 (著名写手)

【答案】应助回帖

内容已删除
赞一下 回复此楼
思想重于技巧，内涵重于表象
11楼2012-04-17 09:55:53
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

minmin_0082003

金虫 (正式写手)

莫愁莫愁
引用回帖:
10楼: Originally posted by gemucai at 2012-04-17 09:53:15:
截断能的收敛测试到了快收敛的时候是有振荡的，纵坐标可以采用前后两次能量之差，这样直观一些。

还要取绝对值对吧
赞一下 回复此楼
12楼2012-04-17 10:01:12
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

minmin_0082003

金虫 (正式写手)

莫愁莫愁
引用回帖:
11楼: Originally posted by souledge at 2012-04-17 09:55:53:
如果是孤立的分子，那每个收敛测试时用relax比较好，否则无法判断分子是否在最稳定态，对比起来没什么意义。
顺便，对于纯粹的孤立分子，如果能用商用软件最好，如Gaussian或Dmol3可以处理非周期分子体系，速度 ...

我想做这个分子吸附在表面上，所以我是不是该做表面的收敛测试，而不是分子的收敛测试。做好表面的收敛测试直接把分子加上就可以了吧。望指教
赞一下 回复此楼
13楼2012-04-17 10:05:32
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zhangguangping

木虫 (著名写手)
引用回帖:
9楼: Originally posted by souledge at 2012-04-17 09:52:19:
偷个懒，这里有详细的说明：http://en.wikipedia.org/wiki/Birch%E2%80%93Murnaghan_equation_of_state

谢谢专家的指导！
我说说我的理解，专家帮我看看对不对。

首先取一个比较大的k点密度，这样保证k点是足够的。然后保持这个k点密度不变，对cutoff值做测试。具体的测试办法是，对应一个cutoff值，不断的改变晶格常数的值，找到总能量和晶格常数的变化关系，然后利用下面的公式：

拟合出公式里面的V0和E0,这样就得到了这个cutoff值对应的V0和E0了。反复这样测试，最后得到cutoff变化时V0和E0的变化了。从而找到一个收敛的cutoff值。是这样的吧？


公式中的B0和B0导数可以从程序输出中得到是吗？


谢谢专家的指导！
赞一下 回复此楼
弘德明志博学笃行
14楼2012-04-17 10:11:53
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

buct2010

金虫 (正式写手)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★
感谢参与，应助指数 +1
minmin_0082003: 金币+8, ★★★很有帮助 2012-04-17 13:46:20
引用回帖:
13楼: Originally posted by minmin_0082003 at 2012-04-17 10:05:32:
我想做这个分子吸附在表面上，所以我是不是该做表面的收敛测试，而不是分子的收敛测试。做好表面的收敛测试直接把分子加上就可以了吧。望指教

截断能我一般都是按照参考文献取的，你找一下类似的体系，看一下别人的设置（赝势的类型也要和你一样），这个取大一些肯定没问题，关键是k点要测试到收敛（以我的经验，用参考文献的值能够确保收敛性，比如文献提到的截断能450，k点5*5*5。我会测试一下截断能400,450,500，此时k点取文献值。确定完截断能后，再试k点，比如4*4*4，5*5*5，6*6*6。你会发现确实是收敛的）。对于你算吸附，先做好表面的截断能、k点测试，把表面的结构弛豫好，再加上要吸附的分子，再弛豫。因为弛豫分子主要是为了得到分子能量，计算吸附能。所以我觉得你没必要用分子去测试截断能。在同一个计算中截断能要保持一致，一般氧的截断能比较高，当你的体系有氧的话，截断能一般是450或500。还有收敛性是看两个计算的能量差值，不是看绝对大小的，楼主没必要为此纠结，很容易解决的。
赞一下 回复此楼
15楼2012-04-17 10:44:32
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

minmin_0082003

金虫 (正式写手)

莫愁莫愁
引用回帖:
15楼: Originally posted by buct2010 at 2012-04-17 10:44:32:
截断能我一般都是按照参考文献取的，你找一下类似的体系，看一下别人的设置（赝势的类型也要和你一样），这个取大一些肯定没问题，关键是k点要测试到收敛（以我的经验，用参考文献的值能够确保收敛性，比如文献提 ...

那拿金表面来说，Au（111）的情况下，一个晶胞俩面只有一个Au原子。 我的表面是3*3 的晶胞四层。也就是有36个Au原子。那我该测谁的kpoint的收敛性呢，是那个含有一个Au原子的晶胞，还是测表面的呢。尽管我是觉得该测表面的，但是还是想确定一下。除了截断能跟kpoint，还有其他需要测试的么
赞一下 回复此楼
16楼2012-04-17 11:00:59
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

souledge

专家顾问 (著名写手)
引用回帖:
14楼: Originally posted by zhangguangping at 2012-04-17 10:11:53:
谢谢专家的指导！
我说说我的理解，专家帮我看看对不对。

首先取一个比较大的k点密度，这样保证k点是足够的。然后保持这个k点密度不变，对cutoff值做测试。具体的测试办法是，对应一个cutoff值，不断的改变晶 ...

张教主太客气也太牛气了，一眼就看穿了整个测试过程，的确就是这么个顺序，不愧是第一性区一号牛人啊！
赞一下 回复此楼
思想重于技巧，内涵重于表象
17楼2012-04-17 11:06:22
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zhangguangping

木虫 (著名写手)
引用回帖:
17楼: Originally posted by souledge at 2012-04-17 11:06:22:
张教主太客气也太牛气了，一眼就看穿了整个测试过程，的确就是这么个顺序，不愧是第一性区一号牛人啊！

确实不敢当，这个说明在第一性原理版块学习的东西多啊。亏了你们这些专家顾问，我们这些新手才能少走不少弯路，也同时长了不少见识。
赞一下 回复此楼
弘德明志博学笃行
18楼2012-04-17 11:10:46
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

minmin_0082003

金虫 (正式写手)

莫愁莫愁
引用回帖:
17楼: Originally posted by souledge at 2012-04-17 11:06:22:
张教主太客气也太牛气了，一眼就看穿了整个测试过程，的确就是这么个顺序，不愧是第一性区一号牛人啊！

你俩稍后再互相感谢呗，我问题还没解决呢。嘻嘻
赞一下 回复此楼
19楼2012-04-17 11:21:42
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

buct2010

金虫 (正式写手)

【答案】应助回帖

引用回帖:
16楼: Originally posted by minmin_0082003 at 2012-04-17 11:00:59:
那拿金表面来说，Au（111）的情况下，一个晶胞俩面只有一个Au原子。 我的表面是3*3 的晶胞四层。也就是有36个Au原子。那我该测谁的kpoint的收敛性呢，是那个含有一个Au原子的晶胞，还是测表面的呢。尽管我是觉得 ...

截断能当然是从开始就要测试的啊，先测试Au的原胞，得到晶格常数再去构建表面。开始是确定截断能，之后不同的的体系（表面、原胞等）都要测试k点的。你多看看侯老师的手册吧。
赞一下 回复此楼
20楼2012-04-17 11:24:42
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 minmin_0082003 的主题更新
252/3返回列表上一页123下一页
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 朋友圈看到的 +4 wangzilk 2026-08-18 5/250 2026-08-18 19:25 by wangzilk
[基金申请] 欢迎发来filecode的Mz6后的代码验证其规律 +51 医学老男孩 2026-08-13 110/5500 2026-08-18 19:07 by icerianxia
[基金申请] 明天放榜？ +5 Shxjjxjkx 2026-08-18 5/250 2026-08-18 18:14 by -大大大大大-
[基金申请] 时间戳变了，能看出什么问题？ +8 基诺咪客 2026-08-17 11/550 2026-08-18 17:13 by 悠然熊猫
[基金申请] 重要消息，中午系统在维护 +9 yuleib84 2026-08-18 10/500 2026-08-18 16:52 by archvillain
[教师之家] 两姐妹，一个辛辛苦苦30岁工作当大学教师，另一个23岁就工作当开车教练 +3 瞬息宇宙 2026-08-11 7/350 2026-08-18 16:42 by LNP@mRNA
[论文投稿] 投稿咨询 +4 wwm09 2026-08-17 6/300 2026-08-18 15:36 by wwm09
[基金申请] 今天维护系统维护 祝所有人 高中 +8 gjjjzhong 2026-08-18 9/450 2026-08-18 13:01 by 家与远方
[基金申请] 哪位老哥知道今年的国自然具体哪一天放榜？ +13 Ldrop2023 2026-08-13 16/800 2026-08-18 12:25 by 淀粉搬运工
[基金申请] 2027广东省杰青 +3 奶牛小黑 2026-08-15 8/400 2026-08-18 11:47 by gltch
[基金申请] 快农历七夕节了，轻松一下，男人悄悄话，女施主请不要进来。 +4 Tide man 2026-08-14 5/250 2026-08-18 11:35 by Tide man
[基金申请] 93BebMhtakh前后11位开头都是大写 +4 且听虎啸 2026-08-17 5/250 2026-08-18 00:49 by 蔡棒棒菂
[基金申请] filecode=后面第一个是大写字母 +8 wangze12014 2026-08-14 10/500 2026-08-17 17:05 by xter9665
[基金申请] 咱们一起用铁证分析2026国家社科基金中标与否 +7 启萌科技 2026-08-12 26/1300 2026-08-16 12:35 by 启萌科技
[精细化工] 招聘 金属平磨液，抛光液研发工程师 +3 小天0311 2026-08-14 3/150 2026-08-16 07:31 by H9PLUS
[基金申请] 各位道友，我要去昆明玩几天，回来见。 +7 Tide man 2026-08-14 8/400 2026-08-15 01:11 by arzu_hma
[基金申请] 重要来源：本周末出结果 +10 瞬息宇宙 2026-08-12 10/500 2026-08-13 15:46 by likettle
[基金申请] 不应该看fileCode +7 且听虎啸 2026-08-12 9/450 2026-08-13 14:27 by flydreamws
[基金申请] Filecode 又变了，巨变 +3 WH3796 2026-08-12 4/200 2026-08-13 14:13 by 小木虫6752397
[基金申请] 结合人工智能，周易传统文化，filecode打分制来了，3分以上希望很大。 +3 Tide man 2026-08-12 4/200 2026-08-13 08:35 by ZJTJZ
信息提示
关闭
请填处理意见
关闭