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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » Vasp&MedeA » 【求助】如何优化单斜晶系反铁磁体系的晶格参数
131/2返回列表12下一页
查看: 1895  |  回复: 12
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看
当前主题已经存档。

hbl626

铁虫 (正式写手)

[交流] 【求助】如何优化单斜晶系反铁磁体系的晶格参数

我在优化单斜晶系反铁磁体系时，如果仅仅进行离子弛豫ISIF=2，是进行自旋极化，还是非自旋极化，结构都会很快的收敛，但如果要进行晶格参数优化，其它参数不变，仅仅将ISIF=3，体系往往优化一段时间后，没有收敛就停下来，这时的晶格参数于是试验值非常接近，但将ＣＯＮＴＣＡＲ　考到ＰＯＳＣＡＲ继续运行的话，晶格参数会变化很大，体积也有原来的７９变成８６.　晶格参数由试验值（ａ＝４．６５，　ｂ＝３．４１　，ｃ＝５．１０８　）变成（ａ＝４．３，　ｂ＝３．８，　ｃ＝５．１７）角度由９９.４８变成９１.７.
请问这是怎么回事？
回复此楼

» 猜你喜欢
1楼 2009-05-09 12:31:01
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

liumiao05

至尊木虫 (知名作家)

秦武汉风
★ ★
qasd(金币+2,VIP+0):谢谢~ 5-13 19:54
ISIF=3 一般就是定性估计，误差比较大
优化晶格常数一般用脚本，固定一个体积，对a,b,c进行优化，然后改变体积，也就是改变你的参考长度，重复上述过程，得到很多个E(V)函数值，拟合得到E-V曲线，能量最低点对应的晶格常数即为所求
赞一下(10人) 回复此楼
天下英雄出我辈，一入江湖岁月催。鸿图霸业谈笑间，不胜人生一场醉。提剑跨骑挥鬼雨，白骨如山鸟惊飞。尘世如潮人如水，只叹江湖几人回。
2楼2009-05-13 00:45:28
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

含草

金虫 (小有名气)
"优化晶格常数一般用脚本，固定一个体积，对a,b,c进行优化，然后改变体积，也就是改变你的参考长度，重复上述过程，得到很多个E(V)函数值，拟合得到E-V曲线，能量最低点对应的晶格常数即为所求" 意思是说通过一个脚本来产生一个POSCAR文件吧,可是如果固体一个体积的话,那里面应该怎么设呢,例如这个是设固定晶格常数

#!/bin/sh
rm WAVECAR
for i in 3.80 3.85 3.90 3.95 4.00 4.05 4.10 4.15 4.20 4.25 4.30
do
cat >POSCAR < Pd-fcc
$i
0.0 0.5 0.5
0.5 0.0 0.5
0.5 0.5 0.0
1
Direct
0.0 0.0 0.0
!
还有就是里面a,b,c的值都不一样,这在脚本文件里面应该如何处理,谢谢
赞一下 回复此楼
3楼2009-05-13 08:30:06
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

liumiao05

至尊木虫 (知名作家)

秦武汉风
★ ★
qasd(金币+2,VIP+0):O(∩_∩)O谢谢 5-13 21:07
可以在INCAR里面把ISIF设为4或者5，固定体积对晶胞进行形状弛豫(就是改变a,b,c优化基矢)，4还同时弛豫离子内坐标(如果不是原胞，内部还有原子)
如果是fcc或者bcc这类简单的结构，就不用对离子弛豫，对每个体积做静态计算即可
赞一下(10人) 回复此楼
天下英雄出我辈，一入江湖岁月催。鸿图霸业谈笑间，不胜人生一场醉。提剑跨骑挥鬼雨，白骨如山鸟惊飞。尘世如潮人如水，只叹江湖几人回。
4楼2009-05-13 20:43:00
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

含草

金虫 (小有名气)
谢谢,不过我还不是很明白,你的意思是说假如:for i in 3.80 3.85 3.90 3.95 4.00 4.05 4.10 4.15 4.20 4.25 4.30,这些值应该是体积的各个固定值,ISIF设为4,那么在脚本文件里面:
$i
0.0 0.5 0.5
0.5 0.0 0.5
0.5 0.5 0.0
应该怎么处理?a,b,c三个值可能都不等,我不知道这样处理才能得到晶格常数.
赞一下 回复此楼
5楼2009-05-13 21:32:37
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

liumiao05

至尊木虫 (知名作家)

秦武汉风
★ ★
wuchenwf(金币+2,VIP+0):谢谢 5-14 22:59
i给定了，初始体积确定了，比如 V=(i^3)*a*b*c，这里的a,b,c是初始值
优化之后变成a1,a2,a3 (CONTCAR里面可以看)，但是这个时候V=(i^3)*a1*b1*c1 ，没有变，因为体积固定了
弛豫内坐标的意思是说如果你的单元有多个原子，那么某些原子在晶胞内的相对位置也可能发生变化，已得到能量更低的结构，ISIF=4可以对这些内部自由度进行弛豫
比如原来有个原子在0.25，0.25，0.25，就可能弛豫成0.24，0.26，0.22
不知道我说清楚了没有
另外我也用了没多久，有不对的或者需要补充的地方请牛人指点
赞一下(10人) 回复此楼
天下英雄出我辈，一入江湖岁月催。鸿图霸业谈笑间，不胜人生一场醉。提剑跨骑挥鬼雨，白骨如山鸟惊飞。尘世如潮人如水，只叹江湖几人回。
6楼2009-05-14 01:39:28
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

hbl626

铁虫 (正式写手)
引用回帖:
Originally posted by liumiao05 at 2009-5-14 01:39:
i给定了，初始体积确定了，比如 V=(i^3)*a*b*c，这里的a,b,c是初始值
优化之后变成a1,a2,a3 (CONTCAR里面可以看)，但是这个时候V=(i^3)*a1*b1*c1 ，没有变，因为体积固定了
弛豫内坐标的意思是说如果你的单元有 ...

这样算出来结果只能是在试验值的基础上，在体积不变的情况下，优化了一下，a,b,c。下一步怎么办？怎样才能得到E~V关系式呢？谢谢
赞一下 回复此楼
7楼2009-05-14 10:01:57
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

liumiao05

至尊木虫 (知名作家)

秦武汉风
再改变i的值，V不就改变了吗，每次都可以得到一个优化能量E，在把得到的各个E对V进行拟合
赞一下 回复此楼
天下英雄出我辈，一入江湖岁月催。鸿图霸业谈笑间，不胜人生一场醉。提剑跨骑挥鬼雨，白骨如山鸟惊飞。尘世如潮人如水，只叹江湖几人回。
8楼2009-05-14 20:34:20
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

含草

金虫 (小有名气)
意思就是说得进行多次运算,在实验值的左右取一些体积的值,分别得到最低能量,最后找出E-V曲线,同时也可以得出晶格常数吧,不知道有没有理解错.那在POSCAR中
$i
0.0 0.5 0.5
0.5 0.0 0.5
0.5 0.5 0.0
那这三个基矢坐标里肯定跟a,b,c有关,是不是在坐标里先给个c/a,b/a的初始值呢
赞一下 回复此楼
9楼2009-05-14 22:03:39
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

liumiao05

至尊木虫 (知名作家)

秦武汉风
★
wuchenwf(金币+1,VIP+0):谢谢 5-14 22:59
是啊，就按实验值或者文献值给吧，给的差别太大了不利于优化，for里面只改变i就可以了
赞一下(10人) 回复此楼
天下英雄出我辈，一入江湖岁月催。鸿图霸业谈笑间，不胜人生一场醉。提剑跨骑挥鬼雨，白骨如山鸟惊飞。尘世如潮人如水，只叹江湖几人回。
10楼2009-05-14 22:07:31
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 hbl626 的主题更新
131/2返回列表12下一页
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 售SCI一区文章，我:8.O.551.O.5.4,科目全,可伽急 +3 6PrO1XyZ5URF 2026-11-13 3/150 2026-11-15 16:12 by 0aDbZAghdVq9
[考博] 售SCI文章，我:8O.5.5.1O.54,科目全,可十急 +3 SX5mJmPWDfOO 2026-11-13 3/150 2026-11-15 08:18 by xSrQNSSjcgoX
[硕博家园] 售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急 +4 SX5mJmPWDfOO 2026-11-13 4/200 2026-11-15 07:34 by xSrQNSSjcgoX
[论文投稿] 售SCI一区文章，我:8.O.551.O.5.4,科目全,可伽急 +4 XJXMuRjD7cK7 2026-11-13 4/200 2026-11-15 07:01 by xSrQNSSjcgoX
[考研] 售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急 +3 uKOSCocCb6Fm 2026-11-13 3/150 2026-11-15 05:11 by xSrQNSSjcgoX
[硕博家园] 售一区SCI文章T0P，我:8O.551.O54,科目全,可十急 +3 uKOSCocCb6Fm 2026-11-12 6/300 2026-11-15 02:59 by xSrQNSSjcgoX
[文学芳草园] 厉娜 +4 张丹枫1991 2026-11-09 6/300 2026-11-15 00:48 by Equinoxhua
[基金申请] 售SCI一区T0P文章，我:8.O.55.1.O.54,科目齐全,可+急 +6 FChKX8meu1ds 2026-11-09 8/400 2026-11-15 00:35 by Equinoxhua
[考博] 售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急 +4 G2rhSvqvhNzl 2026-11-10 6/300 2026-11-15 00:24 by Equinoxhua
[考研] 售SCI一区T0P文章，我:8O.55.1.O.54,科目全,可伽急 +5 YHiw1gou1Eda 2026-11-10 8/400 2026-11-14 22:51 by Equinoxhua
[考研] 售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急 +6 FChKX8meu1ds 2026-11-09 9/450 2026-11-14 22:39 by Equinoxhua
[考博] 27博士 +6 锦鲤送福至 2026-11-10 6/300 2026-11-14 21:23 by xupeii
[考研] 售一区SCI文章T0P，我:8O.551.O54,科目全,可十急 +8 zSbvVX9JFl9J 2026-11-09 9/450 2026-11-13 11:11 by Equinoxhua
[考博] 售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急 +4 TgvkbZZl9tqB 2026-11-11 6/300 2026-11-12 12:42 by ECLTXnJnRSXf
[基金申请] 售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急 +4 Z7qfCsdIiD1T 2026-11-10 7/350 2026-11-12 11:33 by HajghKOMynmz
[教师之家] 售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急 +4 czAiolzR0IJC 2026-11-10 7/350 2026-11-12 11:00 by HajghKOMynmz
[硕博家园] 售SCI文章，我:8O.5.5.1O.54,科目全,可十急 +4 czAiolzR0IJC 2026-11-10 8/400 2026-11-12 10:49 by HajghKOMynmz
[公派出国] 售SCI一区T0P文章，我:8O.55.1.O.54,科目全,可伽急 +6 FChKX8meu1ds 2026-11-09 11/550 2026-11-12 09:57 by HajghKOMynmz
[考研] 售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急 +8 cqv1zKnzd7hr 2026-11-09 12/600 2026-11-12 09:32 by HajghKOMynmz
[公派出国] 售SCI一区T0P文章，我:8O.55.1.O.5.4,科目齐全,可+急 +3 czAiolzR0IJC 2026-11-10 5/250 2026-11-12 06:36 by 9fOTqj074tIg
信息提示
关闭
请填处理意见
关闭