找到ordering后能量最低的结构

vasp计算完后生成100个结构,即100个文件夹,我们需要对文件夹中的vasprun.xml文件进行操作,这个文件几乎包括了所有的信息。
首先我们需要找到能量最低的结构:

from pymatgen import Structure
from pymatgen.io.vasp.sets import MPRelaxSet
from pymatgen.io.vasp import Vasprun
all_energy = [] #python中min()函数的用法是返回列表元素中的最小值,我们先构建一个空的列表
for i in range(100):
    v = Vasprun(r"D:\\one-hundred structures\\" + str(i) + "_relax\\vasprun.xml")  #这里路径中的用变量写法在另一篇博文里面提到过,注意斜杠和反斜杠的区别
    #print(v.final_energy)
    all_energy.append(v.final_energy)  #这里得到的v.final_energy是浮点数类型,可放入列表中比较大小
n = all_energy.index(min(all_energy)) #遍历all_energy列表中的所有元素,找到最小的能量和其在列表中对应的位置
print(n , min(all_energy))

#输出结果:
1 -154.8229777 eV  #这个在列表中的位置与文件名一致;如这里位置是1,就是1_relax这个结构能量最低,结构最稳定
  • 0
    点赞
  • 1
    收藏
    觉得还不错? 一键收藏
  • 0
    评论

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值