计算化学公社

标题: 求助vasp结构优化，高压下晶格膨胀 [打印本页]

作者
Author: 嘿嘿 时间: 2025-9-17 12:51
标题: 求助vasp结构优化，高压下晶格膨胀
我设置pstress=100，施加压力后，晶格还变大了，不知道什么原因，原子之间也不形成键了，这怎么办啊，求助各位老师。 (, 下载次数 Times of downloads: 0)

作者
Author: 嘿嘿 时间: 2025-9-17 12:53
这怎么办呀，不知道是什么原因，我是掺杂进一个Zn，zn替换了gd的位置，我
来做结构优化

作者
Author: kikiw 时间: 2025-9-17 18:25
10GPa这么大的压力。
我觉得你首先考虑下晶胞有没有变化，暂时看来可能是晶胞替换位置后abc扩大三倍
另外Zn和Gd的配位数也不一样

作者
Author: 嘿嘿 时间: 2025-9-17 19:28

kikiw 发表于 2025-9-17 18:25
10GPa这么大的压力。
我觉得你首先考虑下晶胞有没有变化，暂时看来可能是晶胞替换位置后abc扩大三倍
另 ...

是的，但是我实验已经成功合成了，zn掺杂的。但是现在遇到个情况再不加压时，只做结构优化，他也会晶格膨胀，跟之前一样。老师有没有什么方法，可以不让他晶格膨胀呢。现在我是这些参数

ISTART=0
ICHARG=2
ISPIN=1
PREC=Accurate

ALGO=N
NELM=100
EDIFF=1E-5
ENCUT=600

IBRION=2
NSW=400
ISIF=3
EDIFFG=-0.02

ISMEAR=1
SIGMA=0.2

ISYM=0

作者
Author: get-it 时间: 2025-9-17 19:39
估计零压时你就没做晶胞优化吧

作者
Author: 嘿嘿 时间: 2025-9-17 19:53

get-it 发表于 2025-9-17 19:39
估计零压时你就没做晶胞优化吧

对，我0压力时，做晶格优化，就一直膨胀。现在我开启自旋结构优化一下。这些参数一优化，晶胞又膨胀了3倍
1.0
      7.7020001411       0.0000000000       0.0000000000
   -3.8510000706       6.6701277822       0.0000000000
      0.0000000000       0.0000000000       6.6640000343
O Gd Zn S
8 7 1 4
Direct
   0.166669995       0.333330005       0.372000009
   0.333330005       0.166669995       0.628000021
   0.666670024       0.333330005       0.372000009
   0.833329976       0.166669995       0.628000021
   0.166669995       0.833329976       0.372000009
   0.333330005       0.666670024       0.628000021
   0.666670024       0.833329976       0.372000009
   0.833329976       0.666670024       0.628000021
   0.333330005       0.166669995       0.268000007
   0.833329976       0.166669995       0.268000007
   0.666670024       0.333330005       0.731999993
   0.333330005       0.666670024       0.268000007
   0.166669995       0.833329976       0.731999993
   0.833329976       0.666670024       0.268000007
   0.666670024       0.833329976       0.731999993
   0.166669995       0.333330005       0.731999993
   0.000000000       0.000000000       0.000000000
   0.500000000       0.000000000       0.000000000
   0.000000000       0.500000000       0.000000000
   0.500000000       0.500000000       0.000000000
incar文件
ISTART=0
ICHARG=2
PREC=Accurate

ALGO=N
NELM=80
EDIFF=1E-5
ENCUT=600

IBRION=1
NSW=400
ISIF=3
EDIFFG=-0.02

ISMEAR=1
SIGMA=0.2

ISYM=0

kpoints
Auto
0
Gamma
3 4 4
0.0 0.0 0.0

作者
Author: 嘿嘿 时间: 2025-9-17 20:01

嘿嘿发表于 2025-9-17 19:53
对，我0压力时，做晶格优化，就一直膨胀。现在我开启自旋结构优化一下。这些参数一优化，晶胞又膨胀了3倍 ...

这是之前的参数，现在加了自旋，磁矩，来优化

作者
Author: 嘿嘿 时间: 2025-9-17 21:28
加个自旋和磁矩，就好了

作者
Author: 乐平 时间: 2025-9-17 22:05

嘿嘿发表于 2025-9-17 21:28
加个自旋和磁矩，就好了

“加个自旋和磁矩，就好了” 是什么意思? 什么东西好了？
好了是什么意思？晶胞参数不增大了？

作者
Author: 嘿嘿 时间: 2025-9-17 23:14

乐平发表于 2025-9-17 22:05
“加个自旋和磁矩，就好了” 是什么意思? 什么东西好了？
好了是什么意思？晶胞参数不增大了？

对价格自旋和磁矩，晶胞就正常优化了，还收敛了

作者
Author: sobereva 时间: 2025-9-18 02:25

嘿嘿发表于 2025-9-17 12:53
这怎么办呀，不知道是什么原因，我是掺杂进一个Zn，zn替换了gd的位置，我
来做结构优化

有别人回复之前若需要对帖子进行修改、补充，应直接编辑原帖，不要通过回帖进行补充，导致信息零零碎碎，这点在置顶的新社员必读贴里明确说了。

明明是VASP的问题，发帖的时候必须选择VASP分类，别选“其它”，以后注意

作者
Author: 嘿嘿 时间: 2025-9-18 09:41

sobereva 发表于 2025-9-18 02:25
有别人回复之前若需要对帖子进行修改、补充，应直接编辑原帖，不要通过回帖进行补充，导致信息零零碎碎， ...

好的，下次注意

欢迎光临计算化学公社 (http://bbs.keinsci.com/)