We have investigated the thermoelasticity of body-centered cubic (bcc) tantalum from first principles by using the linearized augmented plane wave (LAPW) and mixed--basis pseudopotential methods for pressures up to 400 GPa and temperatures up to 10000 K. Electronic excitation contributions to the free energy were included from the band structures, and phonon contributions were included using the particle-in-a-cell (PIC) model. The computed elastic constants agree well with available ultrasonic and diamond anvil cell data at low pressures, and shock data at high pressures. The shear modulus $c_{44}$ and the anisotropy change behavior with increasing pressure around 150 GPa because of an electronic topological transition. We find that the main contribution of temperature to the elastic constants is from the thermal expansivity. The PIC model in conjunction with fast self-consistent techniques is shown to be a tractable approach to studying thermoelasticity.

To be appear in Physical Review B