计算化学公社

标题: 求助：Pt111的HER氢吸附后构型问题 [打印本页]

作者
Author: lyyylyyyl 时间: 2026-1-28 15:32
标题: 求助：Pt111的HER氢吸附后构型问题
如图是我优化完的Pt111面后的CONTCAR图片，固定了最底两层原子，为何最顶层的Pt发生这么大的跳动，老师说这不正常，我还计算了Fe，Ni等相同结构的纯金属111面，优化后的结构就很正常，没有明显形变，不知是哪出了问题。K点是自动生成的331，附INCAR
# control
SYSTEM=Cu2O_Au
#SYMPREC=1E-10
ISTART=0
ICHARG=2
#IWAVE=1
LCHARG=.F.
LWAVE=.F.
LREAL=.F.
#IVDW=12
NCORE=8

#SCF
IALGO=48
EDIFF=1E-4
ISMEAR=0
SIGMA=0.1
PREC=normal
NELM=150
NELMIN=10
ENCUT=350
LMAXMIX=4

#OPT
EDIFFG=-0.2
NSW=500
IBRION=2 # 0: MD, 1: Newton, 2: CG
ISIF=2 # 3 for bulk, 2 for surface and others
#IDIPOL=3 # 1: x, 2:y, 3, z, 4: all

#DOS
#LORBIT=11
#NBANDS=
#EMIN=-12
#EMAX=12
#NEDOS=1001

#spin- polar
ISPIN=1
#MAGMOM=8*0 6*0 1*0 2*2

作者
Author: UW_0728. 时间: 2026-1-28 19:37
本帖最后由 UW_0728. 于 2026-1-28 19:42 编辑

金属体系显然不应该设置ISMEAR = 0，看VASP wiki上关于ISMEAR的说明；ISMEAR=0算金属体系这种错误的方法得到的结果基本没法用，优化成上面那样离谱的样子也算情理之中。还有，你的EDIFFG = -0.2，这么宽松的收敛限怎么可能保证收敛到预期的极小点结构？设成-0.02才行。

此外，针对结构优化而言，你的EDIFF太松（建议1E-5甚至1E-6），PREC也没有设置成Accurate，这也会对结构优化产生不良影响，虽然这些问题不至于得到这么离谱的结果

作者
Author: lyyylyyyl 时间: 2026-1-28 21:28

UW_0728. 发表于 2026-1-28 19:37
金属体系显然不应该设置ISMEAR = 0，看VASP wiki上关于ISMEAR的说明；ISMEAR=0算金属体系这种错误的方法得 ...

对吧，真的太离谱了，我直接在ms里替换成其他原子，都不会优化成这样。我试试改参数再计算一下。谢谢你！

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