本文整理汇总了Python中ase.utils.eos.EquationOfState.fit0方法的典型用法代码示例。如果您正苦于以下问题：Python EquationOfState.fit0方法的具体用法？Python EquationOfState.fit0怎么用？Python EquationOfState.fit0使用的例子？那么恭喜您, 这里精选的方法代码示例或许可以为您提供帮助。您也可以进一步了解该方法所在类ase.utils.eos.EquationOfState的用法示例。
在下文中一共展示了EquationOfState.fit0方法的1个代码示例，这些例子默认根据受欢迎程度排序。您可以为喜欢或者感觉有用的代码点赞，您的评价将有助于系统推荐出更棒的Python代码示例。
示例1: bulk_summary
# 需要导入模块: from ase.utils.eos import EquationOfState [as 别名]
# 或者: from ase.utils.eos.EquationOfState import fit0 [as 别名]
def bulk_summary(self, plot, a0):
natoms = len(self.atoms)
eos = EquationOfState(self.volumes, self.energies)
v, e, B = eos.fit()
x = (v / self.atoms.get_volume())**(1.0 / 3)
self.log('Fit using %d points:' % len(self.energies))
self.log('Volume per atom: %.3f Ang^3' % (v / natoms))
if a0:
a = a0 * x
self.log('Lattice constant: %.3f Ang' % a)
else:
a = None
self.log('Bulk modulus: %.1f GPa' % (B * 1e24 / units.kJ))
self.log('Total energy: %.3f eV (%d atom%s)' %
(e, natoms, ' s'[1:natoms]))
if plot:
import pylab as plt
plt.plot(self.volumes, self.energies, 'o')
x = np.linspace(self.volumes[0], self.volumes[-1], 50)
plt.plot(x, eos.fit0(x**-(1.0 / 3)), '-r')
plt.show()
bulk = self.atoms.copy()
bulk.set_cell(x * bulk.cell, scale_atoms=True)
self.write_optimized(bulk, e)
return e, v, B, a