结构优化的核心目标是通过调整原子位置或晶格参数,使体系的总能量达到最小值。在VASP中,这一过程通常通过迭代算法实现,例如共轭梯度法(Conjugate Gradient)或准牛顿法(Quasi-Newton)。
优化过程中,系统会根据当前的原子力和能量变化进行调整,直到满足预设的收敛条件(如能量变化小于设定值)。
结构优化的理论基础是密度泛函理论(DFT),其核心是通过交换-相关泛函(如GGA或LSDA)来描述电子结构。在VASP中,常用的交换-相关泛函包括Perdew-Wang 1991(PW91)和PBE等。这些泛函通过近似描述电子之间的相互作用,从而计算体系的总能量和结构。
VASP计算需要准备超算连接软件EASYCONNECT与SSH,建模软件VESTA,以及VASP的输入文件INCAR,KPOINTS,POTCAR
jp-minerals.org/vesta/en/download.html
EasyConnect下载-EasyConnect最新版下载V7.6.7.0
Downloading WinSCP-6.5.3-Setup.exe :: WinSCP
确认结构没有问题后,导出成VASP所需要的POSCAR格式
POSCAR文件可以从本地电脑上传至超算,内容如下
ISTART=0 #开始新的任务,随机产生初始波函数
ICHARG=2 #开始新的任务,从原子电荷密度产生体系初始电荷密度
PREC=M #计算精度,决定ENCUT,FFT的网格大小
EDIFF=1E-5 #相邻两步电子迭代的能量差收敛标准
EDIFFG=-0.1 #离子弛豫的force的收敛标准
ISMEAR=0 #费米能级附近电子占据数为高斯分布,适合金属、半导体、绝缘体
Automatic generation #注释行
判断计算是否完成,看到reach required accuracy代表计算已经完成。
第四步,用VESTA查看计算输出文件CONTCAR
在winscp软件中找到结构优化路径,将计算输出的CONTCAR文件下载
本文给大家介绍pt的晶体结构优化计算。首先需要一个pt的 CIF 文件,下载的位置可以去到网上的各种数据库里面,这是一个非常普通的晶体结构,这是它的打开之后的一个图像,pt是一个面心立方结构,一个晶胞中有四个原子,分别是在立方体的顶角和面心的位置。
大概了解之后,把它导出成计算所需要的 POSCAR 文件,比如就叫pt.vasp。接下来看一下它的结构。它是一个立方的结构,晶格常数接近4,四个原子,一个在顶角上,后面三个是在面心的位置。这个模型做好之后,接下来去到超算,计算位置还是 crystal 下面的 pt 文件夹。
第一步做结构优化,已经做好了一个opt的文件夹,把POSCAR 也放上来了,和刚才是完全一样的。再看一下其他的输入参数,KPOINTS 选择的是4x4x4,根据它的晶格常数来确定。 INCAR中ISPIN=2,因为它是过渡金属可能有磁性,然后 ISIF=3,这个晶格常数是要优化的,pt的POTCAR直接从这个赝势包里复制过来这一个文件就可以。
接下来提交计算,算出来它的磁距几乎为0,但是作为过渡金属,还是设置 ISPIN=2,否则有很大的风险。不是每个都会像pt那样没有磁性的。最后来看一下它的CONTCAR和POSCAR 晶格常数差多少?POSCAR 是3.977,前边是3.962,差了 0.015,这个是非常小的一个偏差,说明这个模型和参数没有什么大的问题。
🏅 华算科技提供专业的第一性原理、分子动力学、生物模拟、量子化学、机器学习、有限元仿真等代算服务。
🎯500+博士团队护航,累计助力5️⃣0️⃣0️⃣0️⃣0️⃣➕篇科研成果,计算数据已发表在Nature & Science正刊及大子刊、JACS、Angew、PNAS、AM系列等国际顶刊。 👏👏👏
声明:如需转载请注明出处(华算科技旗下资讯学习网站-学术资讯),并附有原文链接,谢谢!