24小时热门版块排行榜    

查看: 2919  |  回复: 9

niaofei

铁虫 (初入文坛)

[求助] 怎么优化有两个晶格常数的系统 已有3人参与

请教各位大牛,
比如alpha石英这种六方晶系,有两个晶格常数a,c.
在POSCAR里我是把c的值直接用a的倍数表示,优化a,再优化c/a这个比例,好呢? 还是,先固定住c,优化a,再优化c?
回复此楼

» 收录本帖的淘帖专辑推荐

答疑解惑

» 猜你喜欢

» 本主题相关价值贴推荐,对您同样有帮助:

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

小虫虫LG

木虫 (正式写手)

【答案】应助回帖


感谢参与,应助指数 +1
liliangfang: 金币+1, 谢谢交流 2015-07-03 07:53:38
直接ISIF=3优化不是很好的吗。
2楼2015-07-02 16:18:32
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

hs0218

金虫 (正式写手)

DDr.

【答案】应助回帖


感谢参与,应助指数 +1
liliangfang: 金币+1, 谢谢交流 2015-07-03 07:53:44
类似这种长方体的结构 a a c 结构。
可以先用ISIF=2优化一次原子,再用ISIF=7优化一次格子,之后就可以得到合理的a c 值。
路漫漫其修远兮
3楼2015-07-02 20:47:33
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

niaofei

铁虫 (初入文坛)

引用回帖:
3楼: Originally posted by hs0218 at 2015-07-02 08:47:33
类似这种长方体的结构 a a c 结构。
可以先用ISIF=2优化一次原子,再用ISIF=7优化一次格子,之后就可以得到合理的a c 值。

你的意思是没有必要做循环,找最小能量. 直接跑两次就行乐?
4楼2015-07-03 00:00:41
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

hs0218

金虫 (正式写手)

DDr.

【答案】应助回帖


liliangfang: 金币+1, 谢谢交流 2015-07-04 08:15:43
引用回帖:
4楼: Originally posted by niaofei at 2015-07-03 00:00:41
你的意思是没有必要做循环,找最小能量. 直接跑两次就行乐?...

在文献或者实验里找到a c的值。在自己的设置条件下用 2 和7优化两次得到晶格常数
路漫漫其修远兮
5楼2015-07-03 19:53:36
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xyj880506

银虫 (小有名气)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与,应助指数 +1
liliangfang: 金币+1, 谢谢交流 2015-07-04 08:15:50
niaofei: 金币+15, ★★★很有帮助 2015-07-29 05:31:27
优化三次,分别用ISIF=2,7,2
用isif=7时必须用能量收敛标准
6楼2015-07-03 21:28:22
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

李玥玥

铜虫 (正式写手)

引用回帖:
2楼: Originally posted by 小虫虫LG at 2015-07-02 16:18:32
直接ISIF=3优化不是很好的吗。

想问一下:
1.直接用ISIF=3不会造成晶格变形吗?我在计算缺陷结构的时候,先用ISIF=3优化得到一个结构,然后拿这个结构再优化一次,输入参数都一样,它会继续优化,给出一个新的结构,能量也相差较多。这个该怎么解释呢?我的猜想是因为ISIF=3优化造成的晶格变化不稳定。
2. 如果用ISIF=3,其他参数有没有需要注意特别设置的,使得计算出来的结构更加准确?

碰到过很多直接用ISIF=3的同学,也有老师不推荐ISIF=3,推荐使用ISIF=2727优化的。已晕!!!
7楼2017-07-19 09:35:45
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

李玥玥

铜虫 (正式写手)

引用回帖:
3楼: Originally posted by hs0218 at 2015-07-02 20:47:33
类似这种长方体的结构 a a c 结构。
可以先用ISIF=2优化一次原子,再用ISIF=7优化一次格子,之后就可以得到合理的a c 值。

请问你用ISIF=7优化的时候,有没有什么参数需要注意的呢?我用ISIF=7优化晶格参数的时候总是报错,

ZBRENT:fatal erroe in bracketing
Please rerun with smaller EDIFF, or CONTCAR to POSCAR.
8楼2017-07-19 09:38:54
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

李玥玥

铜虫 (正式写手)

引用回帖:
6楼: Originally posted by xyj880506 at 2015-07-03 21:28:22
优化三次,分别用ISIF=2,7,2
用isif=7时必须用能量收敛标准

请问,”ISIF=7必须设置能量收敛标准“这一说法能不能给解释一下为什么呀?我在用ISIF=7优化时,EDIFFG=-0.01总是不收敛,报错。看到你这个说法,找到一丝希望呀,不过想要知道这个说法的理论方面该怎么解释呢?为何不能在这个时候设置力收敛呢?

看到了希望大神回复一下哈,感谢感谢
9楼2017-07-19 09:42:29
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

啊洪江

新虫 (初入文坛)

遇到和楼主一样的问题,优化四方结构,用ISIF=3,基体优化没问题,但是一掺杂就很难收敛,即使收敛了结构常数也与常理相悖。不知楼主的问题是如何解决了?

发自小木虫Android客户端
10楼2017-08-04 15:48:23
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 niaofei 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 欢迎发来filecode的Mz6后的代码验证其规律 +31 医学老男孩 2026-08-13 68/3400 2026-08-15 23:19 by 医学老男孩
[基金申请] filecode=后面第一个是大写字母 +6 wangze12014 2026-08-14 7/350 2026-08-15 20:12 by gltch
[基金申请] filecode +8 cratir 2026-08-14 12/600 2026-08-15 18:08 by zyfgau
[基金申请] 哪位老哥知道今年的国自然具体哪一天放榜? +7 Ldrop2023 2026-08-13 7/350 2026-08-15 15:57 by mzhh2000
[基金申请] 各位道友,我要去昆明玩几天,回来见。 +7 Tide man 2026-08-14 8/400 2026-08-15 01:11 by arzu_hma
[基金申请] 咱们一起用铁证分析2026国家社科基金中标与否 +7 启萌科技 2026-08-12 22/1100 2026-08-14 23:45 by Noways
[基金申请] 奇怪,两个人的filecode固定段从头到尾一模一样 +8 布布和一二 2026-08-10 11/550 2026-08-14 14:58 by Equinoxhua
[基金申请] 应该是93bebmhtak前后十一个字符比较关键 +23 Lanmanbaby 2026-08-09 37/1850 2026-08-14 13:40 by Equinoxhua
[硕博家园] 请教兼职经验 +3 是阿文鸭 2026-08-09 3/150 2026-08-14 12:07 by HER12025
[硕博家园] 读博的好处 +4 lnee 2026-08-11 4/200 2026-08-14 10:20 by ahsoarli
[基金申请] filecode +15 documentary 2026-08-10 17/850 2026-08-14 10:08 by kissu88
[基金申请] FileCode能看出啥? +10 要乐观耀哥 2026-08-10 32/1600 2026-08-14 09:37 by 要乐观耀哥
[文学芳草园] 阿姨 +4 汪汪锅 2026-08-09 4/200 2026-08-13 19:43 by arzu_hma
[基金申请] 静等基金结果 +8 gjjjzhong 2026-08-10 21/1050 2026-08-13 17:56 by 且听虎啸
[基金申请] 重要来源:本周末出结果 +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
[基金申请] 分享一下我之前已中青C的计划书的filecode +4 布布和一二 2026-08-11 5/250 2026-08-13 12:56 by cratir
[基金申请] 国自然结果 +4 Vierhys 2026-08-10 8/400 2026-08-10 15:06 by Vierhys
[基金申请] 面上项目filecode邪修 +5 西山十月 2026-08-09 7/350 2026-08-10 07:32 by 仁砚薪传
信息提示
请填处理意见