24小时热门版块排行榜    

查看: 2506  |  回复: 6

wenj674

金虫 (小有名气)

[交流] 几何优化的对称性约束问题 已有3人参与

各位大侠:
在晶体结构几何优化的时候,把其对称性去掉优化的效果好还是保留原来对称性约束优化的效果好呢?
请各位有经验的大侠不吝赐教啊!!!!
回复此楼

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

VASP

» 猜你喜欢

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

历史长河一条线,每段重现几百遍。你笑我疯,他说我贱,我都假装看不见。
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

buct2010

金虫 (正式写手)

★ ★ ★
小木虫(金币+0.5):给个红包,谢谢回帖
liliangfang(金币+2, 专家考核): 多谢指教 2011-12-30 10:00:34
在计算纯净晶体的时候建议保留对称性,因为纯净晶体是具有空间操作对称性的

掺杂、空位缺陷会改变晶胞对称性,所以计算时应该取消对称性,这样得到的结构更合理
"GaN1−xSbx ... the individual atomic positions have been fully optimized without any symmetry constraints.” 这篇文章计算掺杂时就是这样做的。Visible-light absorption and large band-gap bowing of GaN1−xSbx from first principles PHYSICAL REVIEW B 84, 075304 (2011)
2楼2011-12-30 09:24:15
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

信任开花

铁杆木虫 (知名作家)

★ ★ ★
小木虫(金币+0.5):给个红包,谢谢回帖
WDD880227(金币+2): 感谢交流~~~有帮助~~ 2011-12-30 15:34:52
引用回帖:
2楼: Originally posted by buct2010 at 2011-12-30 09:24:15:
在计算纯净晶体的时候建议保留对称性,因为纯净晶体是具有空间操作对称性的

掺杂、空位缺陷会改变晶胞对称性,所以计算时应该取消对称性,这样得到的结构更合理
"GaN1−xSbx ... the individual at ...

计算空位型缺陷的时候要先优化一下,步骤:
1. 完美超胞(2*2*2)先优化,取isif=3优化
2. 取上一步骤优化好的原子坐标文件CONTCAR,删除其中一个原子坐标,得到空位型缺陷的初始坐标文件,然后用isif=4优化,以期得到优化驰豫的空位型缺陷的结构

问题是isif优化后OUTCAR文件中的其中一行
external pressure =       -8.69 kB  Pullay stress =        0.00 kB

这个external pressure 有点大,是否是isif选择的问题

可是大家在计算缺陷的时候不是说缺陷的存在不影响晶格参数的大小,这个时候优化应该用isif=4,即固定体积优化么?

还有计算的时候ISYM参数取默认值,是否也要改为ISYM=0

期待专家解答,非常感谢
从不温柔,从不体贴,从不讲理
3楼2011-12-30 13:08:17
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

buct2010

金虫 (正式写手)

★ ★
小木虫(金币+0.5):给个红包,谢谢回帖
fzx2008(金币+1): 谢谢指教 2011-12-30 14:21:59
引用回帖:
3楼: Originally posted by 信任开花 at 2011-12-30 13:08:17:
计算空位型缺陷的时候要先优化一下,步骤:
1. 完美超胞(2*2*2)先优化,取isif=3优化
2. 取上一步骤优化好的原子坐标文件CONTCAR,删除其中一个原子坐标,得到空位型缺陷的初始坐标文件,然后用isif=4优化, ...

晶格常数优化采用ISIF=3
结构弛豫一般用ISIF=2,以保证基矢、倒格矢、布里渊区不变。你说的ISIF=4我没用过,不知道在什么场合适用。因为我关注能带结构,所以这样设置。
对称性是用ISYM来 控制
4楼2011-12-30 14:20:22
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

信任开花

铁杆木虫 (知名作家)

★
小木虫(金币+0.5):给个红包,谢谢回帖
引用回帖:
4楼: Originally posted by buct2010 at 2011-12-30 14:20:22:
晶格常数优化采用ISIF=3
结构弛豫一般用ISIF=2,以保证基矢、倒格矢、布里渊区不变。你说的ISIF=4我没用过,不知道在什么场合适用。因为我关注能带结构,所以这样设置。
对称性是用ISYM来 控制

非常感谢专家的回复

对于二楼你的回帖,是否可以做如下理解:

对于含空位的结构驰豫,ISYM=0才是更符合实际的情况
从不温柔,从不体贴,从不讲理
5楼2011-12-30 15:05:55
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

buct2010

金虫 (正式写手)

★ ★
小木虫(金币+0.5):给个红包,谢谢回帖
liliangfang(金币+1): 谢谢回复 2011-12-31 11:03:08
引用回帖:
5楼: Originally posted by 信任开花 at 2011-12-30 15:05:55:
非常感谢专家的回复

对于二楼你的回帖,是否可以做如下理解:

对于含空位的结构驰豫,ISYM=0才是更符合实际的情况

嗯,可以这么理解
6楼2011-12-31 10:38:45
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

09131305

新虫 (初入文坛)

★
小木虫: 金币+0.5, 给个红包,谢谢回帖
请教::优化2*3的CeO2(111)面时。。。对称性应该怎么设置??和2*2时一样么??
7楼2013-11-05 16:06:01
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 wenj674 的主题更新
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考博] 27博士申请 +3 白白白ly 2026-11-09 3/150 2026-11-09 21:39 by szq0819
[考博] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +3 zSbvVX9JFl9J 2026-11-09 3/150 2026-11-09 21:05 by FChKX8meu1ds
[文学芳草园] 厉娜 +3 张丹枫1991 2026-11-09 5/250 2026-11-09 16:53 by 张丹枫1991
[考博] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +5 kvewwGtKktqB 2026-11-04 12/600 2026-11-09 13:26 by Equinoxhua
[考研] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +5 ttH3f1oA6ace 2026-11-04 15/750 2026-11-09 13:19 by Equinoxhua
[硕博家园] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +5 5sJFDKp2kwuI 2026-11-06 16/800 2026-11-09 11:47 by GwrQVVR6Oenx
[考博] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +5 5sJFDKp2kwuI 2026-11-06 15/750 2026-11-09 11:41 by GwrQVVR6Oenx
[硕博家园] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +4 5sJFDKp2kwuI 2026-11-06 10/500 2026-11-09 08:00 by IhnEQauGmT3K
[考研] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +3 5sJFDKp2kwuI 2026-11-06 12/600 2026-11-09 06:32 by l48hFd8Nwae2
[公派出国] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +4 5sJFDKp2kwuI 2026-11-06 8/400 2026-11-09 00:14 by l48hFd8Nwae2
[教师之家] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +4 5sJFDKp2kwuI 2026-11-06 8/400 2026-11-08 23:01 by l48hFd8Nwae2
[教师之家] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +5 5sJFDKp2kwuI 2026-11-06 9/450 2026-11-08 22:46 by l48hFd8Nwae2
[教师之家] 延安大学回应“副教授留绝笔信”:确已离世,正核实举报内容 +3 yexuqing 2026-11-04 5/250 2026-11-07 13:22 by 苏东坡二世
[公派出国] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +5 ttH3f1oA6ace 2026-11-04 10/500 2026-11-06 16:12 by 5sJFDKp2kwuI
[硕博家园] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +4 ttH3f1oA6ace 2026-11-04 11/550 2026-11-06 16:10 by 5sJFDKp2kwuI
[公派出国] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +6 ttH3f1oA6ace 2026-11-04 12/600 2026-11-06 16:10 by 5sJFDKp2kwuI
[公派出国] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +5 kvewwGtKktqB 2026-11-04 14/700 2026-11-06 05:10 by qL1cU46YIdK0
[考研] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +4 ttH3f1oA6ace 2026-11-04 10/500 2026-11-06 02:04 by qL1cU46YIdK0
[硕博家园] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +3 ttH3f1oA6ace 2026-11-04 9/450 2026-11-06 02:02 by qL1cU46YIdK0
[硕博家园] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +4 hNpp4Y2CGvSC 2026-11-04 12/600 2026-11-06 01:53 by qL1cU46YIdK0
信息提示
请填处理意见