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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » Vasp&MedeA » 【求助】如何优化单斜晶系反铁磁体系的晶格参数
北京石油化工学院2026年研究生招生接收调剂公告
131/2返回列表12下一页
查看: 1677  |  回复: 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下一页
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 材料调剂 +10 一样YWY 2026-04-06 10/500 2026-04-06 21:05 by lbsjt
[考研] 0703化学调剂 348分 +7 唉我超真没招了 2026-04-06 7/350 2026-04-06 19:50 by lijunpoly
[考研] 286求调剂 +12 Faune 2026-04-06 12/600 2026-04-06 17:40 by dongzh2009
[考研] 求调剂 +4 熊二想上岸 2026-04-06 4/200 2026-04-06 17:17 by dongzh2009
[考研] 求调剂 +10 chenxrlkx 2026-04-05 10/500 2026-04-06 11:31 by 猪会飞
[考研] 一志愿武汉理工大学-085601材料工程（专硕）-总分353求调剂 +3 2626262626li 2026-04-02 3/150 2026-04-06 09:08 by 无际的草原
[考研] 085600，320分求调剂 +7 大馋小子 2026-04-01 8/400 2026-04-05 21:19 by 学员8dgXkO
[考研] 296求调剂 +3 汪！？！ 2026-04-05 4/200 2026-04-05 20:13 by 啵啵啵0119
[考研] 考研调剂生寻找导师 +3 顾瞻考研啊 2026-04-05 3/150 2026-04-05 18:18 by 啵啵啵0119
[考研] 081200-11408-276学硕求调剂 +4 崔wj 2026-04-04 5/250 2026-04-05 14:06 by imissbao
[考研] 294求调剂 +6 Grey_Ey 2026-04-02 9/450 2026-04-04 22:07 by hemengdong
[考研] 一志愿沪9，求生物学调剂，326分 +6 刘墨墨 2026-04-04 6/300 2026-04-04 19:44 by 唐沐儿
[考研] 22408，264求调剂 +3 ywh729 2026-04-03 4/200 2026-04-04 11:04 by ywh729
[考研] 274求调剂 +9 顺理成张 2026-04-03 10/500 2026-04-03 15:10 by 啊俊！
[考研] 319求调剂 +18 太容易1018 2026-04-01 18/900 2026-04-03 11:18 by linyelide
[考研] 一志愿复旦材料，英一专硕，总分357调剂 +4 1050389037 2026-04-02 5/250 2026-04-02 21:40 by dongzh2009
[考研] 材料340分调剂 +7 夏夜晚风_long 2026-04-02 9/450 2026-04-02 21:20 by dongzh2009
[考研] 362求调剂 +14 西南交材料专硕3 2026-03-31 14/700 2026-04-02 17:50 by yunlongyang
[考研] 270调剂 +7 maxjxbsk 2026-04-02 7/350 2026-04-02 09:50 by yulian1987
[考研] 一志愿北京科技大学085601材料工程英一数二初试总分335求调剂 +5 双马尾痞老板2 2026-03-31 5/250 2026-04-01 09:04 by oooqiao
信息提示
关闭
请填处理意见
关闭