def energyForce(self):
        self.getVaspRun_vasp()
        forces = parseVasprun('forces')
        stress = parseVasprun('stress')
        c = shell_exec("grep TOTEN OUTCAR|tail -1")
        energy = sscanf(c, "free  energy   TOTEN  =      %f eV")[0]
        return forces, stress, energy