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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » QE(Pwscf) » 收敛测试，这种情况我该继续测试么,简直就是跳跃式的嘛
252/3返回列表上一页123下一页
查看: 4661  |  回复: 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下一页
普通表情 龙 兔 虎 猫
信息提示
关闭
请填处理意见
关闭