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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » QE(Pwscf) » 收敛测试，这种情况我该继续测试么,简直就是跳跃式的嘛
252/3返回列表上一页123下一页
查看: 4734  |  回复: 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下一页
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 化学工程321分求调剂 +9 大米饭！ 2026-03-15 12/600 2026-03-16 21:59 by 大米饭！
[考研] 药学383 求调剂 +3 药学chy 2026-03-15 4/200 2026-03-16 20:51 by 元子^0^
[基金申请] 今年的国基金是打分制吗？ 50+3 zhanghaozhu 2026-03-14 3/150 2026-03-16 17:07 by 北京莱茵润色
[教师之家] 焦虑 +7 水冰月月野兔 2026-03-13 9/450 2026-03-16 10:00 by Quakerbird
[考研] 求老师收留调剂 +4 jiang姜66 2026-03-14 5/250 2026-03-15 20:11 by Winj1e
[考研] 本科南京大学一志愿川大药学327 +3 麦田耕者 2026-03-14 3/150 2026-03-14 20:04 by 外星文明
[基金申请] 面上和青基一样限30页不合理 +5 wowsunflower 2026-03-10 7/350 2026-03-14 17:21 by kingkocxr
[考研] 255求调剂 +3 李嘉慧， 2026-03-12 4/200 2026-03-14 16:58 by 有只狸奴
[考研] 材料080500调剂求收留 +3 一颗meteor 2026-03-13 3/150 2026-03-14 10:54 by peike
[考研] 材料工程专硕，一志愿中国矿业大学，总分314，求调剂 +5 无懈可击的巨人 2026-03-10 5/250 2026-03-14 00:37 by JourneyLucky
[考研] 2026考研调剂+本科延边大学+山东大学+生物化学与分子生物学+有项目经验 +3 ccdsscjy 2026-03-10 3/150 2026-03-14 00:12 by JourneyLucky
[考研] 341求调剂 +3 番茄头--- 2026-03-10 3/150 2026-03-13 23:07 by JourneyLucky
[考研] 285 求调剂 资源与环境 一志愿北京化工大学 +3 未名考生 2026-03-10 3/150 2026-03-13 23:04 by JourneyLucky
[考研] 求材料调剂 +5 隔壁陈先生 2026-03-12 5/250 2026-03-13 22:03 by 星空星月
[考研] 0856材料与化工301求调剂 +5 奕束光 2026-03-13 5/250 2026-03-13 22:00 by 星空星月
[考研] 一志愿西南交大，材料专硕317求调剂 +5 lx8568 2026-03-11 5/250 2026-03-13 21:43 by peike
[考研] 考研调剂 +4 芬达46 2026-03-12 4/200 2026-03-13 16:04 by ruiyingmiao
[考研] 308求调剂 +3 是Lupa啊 2026-03-12 3/150 2026-03-13 14:30 by 求调剂zz
[考研] 工科0856专硕化学工程269能调剂吗 +10 我想读研11 2026-03-10 10/500 2026-03-13 10:14 by Yuyi.
[考研] 0703化学调剂 +3 三dd. 2026-03-10 3/150 2026-03-10 15:45 by peike
信息提示
关闭
请填处理意见
关闭