原文地址:birch-m rnaghan作者:海角鱼 对al晶格常数进行优化过程中,根据侯老师的指南得到etvo.dat文件,该文件是关于体积与能量的对应关系,要得到优化了的晶格常数,需要利用birch-m rnaghan状态方程进行拟合.拟合过程可通过m rn.f小程序得以实现.m rn.f的编译过程如下: g77 -o m rn.x m rn.f 得到一个可执行文件m rn.x;该可执行文件的运行格式为:m rn.x inp tfile o tp tfile; inp tfile的格式如下: 文件格式:inp.m 内容大概为: 2 # 表示下面输入的能量单位;ry-1,ev-2,hartree-3 0.25 # 原胞体积晶胞体积 3.80 4.30 50 # 晶格常数的最小,最大值以及要拟合的点数目 10 # 计算了多少晶格常数用来拟合 3.80 -3. # 以下几行为计算得到的晶格常数与能量的对应值 . . .  # 不知道是否可以不要??待解.. 准备好inp.m后,运行m rn.x: $: .m rn.x inp.m o t.m 得到的o t.m包含所以拟合的结果信息,可通过vi编辑查看.