private double latticePar(double T, double par0, double TEC, double Tch) {
        if (temperature == 0) {
            return par0;
        } else {
            return par0 * (1 + TEC * Tch * Utils.debye(Tch / T));
        }
    }