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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » MS » 【求助】k-piont 怎么确定的啊？
172/2返回列表上一页12
查看: 1664  |  回复: 16
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

luxia0210

铁杆木虫 (著名写手)

GG
★
zmrright(金币+1):谢谢参与
，在赚个金币哈
回复此楼
潜心百年，为求一败！
11楼2010-05-04 18:50:39
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

acridine

木虫 (正式写手)
★ ★ ★ ★
小木虫(金币+0.5):给个红包，谢谢回帖交流
xiaohunhun(金币+3):xiexie 2010-05-04 20:55:00
楼主的语序有点乱啊，似乎很匆忙。
>1. 我在算LiBH4,但是我的k-point 总是弄不好，我想应该是与晶体的晶格参数有关，是不是啊
k的设置不仅仅和晶格参数有关，与你计算的材料、要讨论性质等等都有关系。准确的设置k-piont就如同楼上各位提到的“收敛性测试”

>2. 我在 run的时候总是报错，我想应该是k-point 出错了。我改变为小的就好了。。。但是咋确定呢。。。？
这里给出的信息太少，在给出错误信息之前，确实无法确定是否是k点设置的问题。比如内存不足也会出现这种现象：高精度出错，降低精度才可以计算下去。
赞一下(2人) 回复此楼
12楼2010-05-04 19:36:27
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

awmc2008

至尊木虫 (文坛精英)

学习者
★
小木虫(金币+0.5):给个红包，谢谢回帖交流
学习学习，谢谢高手。
赞一下(1人) 回复此楼
xuexizhe
13楼2010-05-04 22:48:14
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zmrright

金虫 (著名写手)

谢谢，我搜一下。。。

引用回帖:
Originally posted by guifan at 2010-05-03 22:13:50:
忘了提醒一下，可以关注金属所的胡青苗教授，直接检索他的文章就能发现

谢谢你好。。。
回复此楼
做好自我，不枉此生
14楼2010-05-04 23:21:44
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zmrright

金虫 (著名写手)

谢谢

引用回帖:
Originally posted by acridine at 2010-05-04 19:36:27:
楼主的语序有点乱啊，似乎很匆忙。
>1. 我在算LiBH4,但是我的k-point 总是弄不好，我想应该是与晶体的晶格参数有关，是不是啊
k的设置不仅仅和晶格参数有关，与你计算的材料、要讨论性质等等都有关系。准确的 ...

我尝试一下。。。谢谢
回复此楼
做好自我，不枉此生
15楼2010-05-04 23:23:22
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zmrright

金虫 (著名写手)

刚搜集的，虽然没太懂，但是大家商量看。。。

★ ★
aylayl08(金币+2):鼓励一下讨论 2010-05-07 19:55:52
大家好：

我根据自己在计算过程中遇到的收敛性的问题，做了一点关于CASTEP计算参数调试的总结。当然，这些只是我的个人心得体会，大家如果觉得有什么问题请提出指正，或是如果有自己独到的见解，请跟帖加上来，供大家一起学习提高。

如果计算在Max. cycle全部用完之后仍然不能收敛，要先查看一下.castep文件，如果能量比刚开始计算时有明显下降，而且震动幅度也明显减小，只是还没有降到收敛标准以下，那说明Max. cycle还不够，要增加这个数值。

但多数情况下，我们会发现，能量并没有比刚开始计算时的能量有明显下降，而是围绕着某个平均值在反复地震荡，而且震动幅度也没有减小。那说明是参数设定的问题，此时就算再怎么增加Max. cycle的数值，计算最终还是不会收敛的。

这个时候，最好是增加energy cutoff的数值，降低K point set的数值，增加smearing width数值（当然，如果通过调试前两个数值就可以使计算收敛的话，就最好不要动这个数值）。

再有，如果计算是使用的density mixing的话，那要将mix_charge_amp，mix_spin_amp的数值降低，将mix_energy_cutoff的数值增加到energy cutoff数值的3到4倍，增加mix_history_length的数值。

当然，如果调试上面的参数仍然发现收敛很辛苦的话，就将mixing scheme转成All Band/EDFT,这会解决收敛问题，但是计算时间将会是使用density mixing计算所用时间的3倍以上。


总结的挺好的，但是指出几点不足，即：

1. 参数调试过于经验化，如果楼主能结合对程序原理的认识， 或许能够把每个参数调节的意义给阐述的更为清晰些； 从优化算法， 体系特性， 物理过程等方面着手来进行分析， 或许可以得到更加全面的认识（如果楼主能够结合上述内容对CASTEP收敛进行详细的分析和总结的话，我会考虑接受楼主的文章，并且发表在小木虫计算与模拟研究RCS期刊上）。

2. 建议楼主在进行平面波DFT优化之前，通过其他的办法， 比如经典力学优化以及DMol中局域波函数进行优化， 这样的话，可以帮助CASTEP在最短的时间内找到合适的最小值， 也就是通常意义上说的加快收敛。


great!
其实在做收敛速度无非是两个目的：１，计算精度要高；２计算速度要快．
一般平面波程序，在计算精度上，主要需要测试ecut(倒空间的最大平面波截断，理论上要求无穷大，因为那样基矢才是完备的．但实际上不需要，也做不到)．还有就是ｋ点取样．理论上要计算无穷个ｋ点，但实际也不可能．我觉得应该先测试ecut,因为你首先要保证每个ｋ点能量是正确的．然后再测试ｋ点．其实到了，最后即使ｋ点取的很大的时候，还是有很大的误差的，这主要来自于系统误差．所以有时候，会造成震荡的情况，这主要还是由于ｅｃｕｔ太小，或者是有的ｋ点根本就难以算准．这个时候，如无特殊需要，基本就认为它收敛了．除非你想要在０．１ｍｅＶ的误差．　
但是根据我的想法，我就不认为它能将计算精度控制在０．１ｍｅＶ．即使收敛测试达到了．（当然,这是个人看法！）
第二，就是计算速度．因为目前ＤＦＴ主要还是用自洽的方法，在理论上应该是用新的电荷密度等于旧的电荷密度，但实际上如果这样做，一是收敛很慢，二是容易形成震荡情形。一般的做法是rho_in=mix*tho_old+(1-mix)*rho_old。（具体到每个程序，肯定是做法不同了。）一般认为mix=0.3比较合理。这样收敛的非常快。　　所以很多程序在测试各种mix的时候，就是在测试这个。
最后，选择合适的参数。在你所需要的精度范围内，最快的计算。　我们的目的是要解决问题，有时候不需要那么精确。



引用：“围绕着某个平均值在反复地震荡，而且震动幅度也没有减小。
这个时候，最好是增加energy cutoff的数值，降低K point set的数值，增加smearing width数值（当然，如果通过调试前两个数值就可以使计算收敛的话，就最好不要动这个数值）。”

我觉得不是减小Ｋ，而应当增加Ｋ
如果不是金属体系，调节ｓｍｅａｒｉｎｇ没有意义
如有不妥，还望不吝赐教：Ｐ
赞一下(1人) 回复此楼
做好自我，不枉此生
16楼2010-05-05 00:03:57
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wzh上善若水

木虫 (著名写手)
★
小木虫(金币+0.5):给个红包，谢谢回帖交流
需要进行进行收敛测试，可查看相关文章学习收敛性测试
赞一下(1人) 回复此楼
知 止 而 后 有 定 ； 定 而 后 能 靜 ； 靜 而 后 能 安 ； 安 而 后 能 慮 ； 慮 而 后 能 得
17楼2010-05-09 20:38:02
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 zmrright 的主题更新
172/2返回列表上一页12
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
信息提示
关闭
请填处理意见
关闭