Computation of finite temperature mechanical properties were performed as detailed in our recent study. Two approaches were used in this work: the finite strain difference approach and the strain-fluctuation method. We focused on the example of the bulk modulus K, defined as K = -V (∂P /∂V ) T . Values of K reported in this work correspond to the equilibrium volume.