24小时热门版块排行榜    

查看: 3613  |  回复: 21
【有奖交流】积极回复本帖子,参与交流,就有机会分得作者 identation 的 49 个金币 ,回帖就立即获得 2 个金币,每人有 1 次机会
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

identation

金虫 (著名写手)

[交流] 【求助】请教PWSCF晶体建模和自洽计算,不收敛已有6人参与

大家好,想请教一下晶体结构建模的事。BTO是I4/m空间群,a=10.36666Bohr,c=1.413066*10.36666Bohr;超晶胞有20个原子,原胞为10个原子;超晶胞的原子坐标如下(从MS中导出):
为什么我计算的结构不收敛呢?

[ Last edited by identation on 2010-4-11 at 11:48 ]
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

goldenfisher

金虫 (著名写手)

★ ★ ★ ★ ★
gavinliu7390:谢谢解答! 2010-04-08 15:24
aylayl08(金币+5):谢谢解答 2010-04-08 18:31
identation(金币+10): 2010-04-11 20:31
引用回帖:
Originally posted by enola at 2010-04-08 11:41:33:

但是我有个问题啊?这个结构由应力的,计算出来的结构会不会有影响呢?
特别是超胞?
每次计算,如果对自己的结构不确定结构有没有到最优,我们有没有必要都先做relax呢?
谢谢

最好是做一下优化。因为程序的处理总会出现一些精度的问题,而且原子数目越多,这种对结构的预测的误差就越大,每个程序、每个参数对应的最优结构只有一个。有时候一点小小的误差也会导致结果的不同。
对于只包含原胞的体系,优化出来的结果差不多,应该做的是体积优化,也就是说找到最优的晶格常数。
13楼2010-04-08 15:10:30
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 22 个回答

goldenfisher

金虫 (著名写手)

★ ★ ★ ★
identation(金币+5):谢谢参与
gavinliu7390(金币+3):谢谢参与! 2010-04-07 17:08
identation(金币+25): 2010-04-07 17:37
第一: ecutwfc = 30, 小了点。你的Ba原子应该大一些
第二: occupations='fixed', 不知道你为什么要选用fixed,可以试一试smearing。
第三:仔细检查你的原子坐标,感觉有点怪,你可以用xcrysden把你的构型画出来看看到底对不对。
建议,verbosity = 'low' ,去掉
具体的要下来看看才知道。
2楼2010-04-07 16:21:42
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

gavinliu7390

木虫 (著名写手)

叶落鹰飞


identation(金币+5):谢谢参与
identation(金币+5): 2010-04-07 17:37
mixing_beta = 0.3 ,   一般这个都用0.7 , 我估计这就是 你不收敛的原因.
ecutwfc = 30 ,   这个肯定是小.
我觉得K点也好像有点大.不过不好说.
收敛测试是一定要的.
真理是一点点接近的!
3楼2010-04-07 17:10:28
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lono75

铁杆木虫 (著名写手)

★ ★ ★ ★
identation(金币+5):谢谢参与
gavinliu7390(金币+3):谢谢测试! 2010-04-07 18:33
identation(金币+15):但是ibrv=7是体心四方呀,我取得是20个原子的supercell,堆起来不是应该是简单四方点阵吗?谢谢 2010-04-07 22:46
若体系是绝缘体,设置occupations='fixed'
ibrav =7
原子位置没问题
提高赝势能量截断
计算时混合因子 0.7 为默认值
增大mixing_ndim 及减小  mixing_beta 可使体系易收敛
otherelse,+U
if not convergance is still,
then the structure maybe have problem calculated by the method.

[ Last edited by lono75 on 2010-4-11 at 19:23 ]
由天命,成人事!
4楼2010-04-07 18:23:05
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
提示: 如果您在30分钟内回复过其他散金贴,则可能无法领取此贴金币
普通表情 高级回复(可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 科研之友最近访客那出现来自北京的未知用户,什么情况? +9 lihangshi 2024-06-28 13/650 2024-06-29 12:37 by zhengxianzb
[基金申请] 能看出是否上会了吗 +11 articlefan 2024-06-23 16/800 2024-06-29 12:34 by sunbetty
[硕博家园] 有机合成硕士毕业8年,想读博士,但又无法全职 +5 青灯孤案 2024-06-26 7/350 2024-06-29 12:16 by wod河畔
[访问学者] 公派访问学者申请还是审核中吗? +3 326lhpqk 2024-06-28 4/200 2024-06-29 12:12 by 加菲的树袋熊
[硕博家园] 博士该不该读? +13 L1009225316 2024-06-25 15/750 2024-06-29 09:58 by 自强2019
[基金申请] 我劝基金委要善良+一个大胆的设想 +13 我是大王 2024-06-28 14/700 2024-06-29 09:32 by Higraduate
[基金申请] 2024基金必中 +5 fzy721521 2024-06-28 5/250 2024-06-29 09:20 by 德尚中行
[基金申请] 2020面地青已有会评专家接到电话了 +10 dream200606 2024-06-28 12/600 2024-06-29 08:25 by llhljsy
[找工作] offer选择 +13 羡渔樵 2024-06-26 13/650 2024-06-28 23:59 by 高敖曹
[药学] 大佬们来讨论下,为什么分子里加入了有毒性的官能团分子反而没毒了 +4 HF111001 2024-06-24 5/250 2024-06-28 22:45 by 渣男2019
[基金申请] 请问三个字的中文名字,在文中到底怎么表示? 5+3 zhongyuan764 2024-06-22 8/400 2024-06-28 21:56 by zhongyuan764
[博后之家] 75批博士后面上项目公布了 +4 ishow1989 2024-06-28 4/200 2024-06-28 20:24 by sizhouyi
[基金申请] 有当年结题的在研国基对申请影响大不大 1+4 sunjc 2024-06-28 13/650 2024-06-28 17:09 by lijian7338
[育儿交流] 暑假孩子一个人在家可以干什么 +11 唠叨的映 2024-06-25 12/600 2024-06-28 12:53 by George9278
[基金申请] 评审专家会不会很在意申请人的单位啊 +11 lancet0903 2024-06-24 14/700 2024-06-28 09:45 by xingmoon0869
[论文投稿] 论文润色 +3 自强不息a?a 2024-06-26 3/150 2024-06-28 01:23 by Cherry_1215
[有机交流] 有机物的酸酐如何准确表征 10+3 方酱 2024-06-24 7/350 2024-06-26 17:53 by 宁静远行
[找工作] 老博士选高校 5+5 woodhead8386 2024-06-24 8/400 2024-06-26 16:35 by 大大熊呀
[基金申请] 2024安徽省哲社立项名单公示 +3 robin_work 2024-06-24 6/300 2024-06-25 14:37 by Pickfoot
[论文投稿] OSA期刊审稿逾期 +3 Thomas_Squid 2024-06-22 3/150 2024-06-23 15:20 by wspglt
信息提示
请填处理意见