def _calc_temp(self, V_a, T0=None):
        if T0 is None:
            T0 = self.eos_mod.refstate.ref_temp()

        a1, a2 = self._calc_strain_coefs()
        fstr_a = self._calc_fstrain(V_a)
        T_a = T0*np.sqrt(1 + a1*fstr_a + 0.5*a2*fstr_a**2)
        return T_a