Density-functional study of the thermodynamic properties and the pressure-temperature phase diagram of Ti

The phonon spectra of α, β, and ω Ti were studied using the supercell approach. The lattice vibrational energy was calculated in the quasiharmonic approximation using both first-principles phonon density of state and Debye model. The thermal electronic contribution to the free energy was evaluated from the integration over the electronic density of states. The Helmholtz energy was thus obtained by combining them with the 0 K total energy calculated within the framework of all-electron projector-augmented-wave method. The thermodynamic properties of α and ω Ti calculated by phonon and Debye model are very close to each other. The predicted enthalpy, entropy, bulk modulus, thermal expansion coefficient and heat capacity of α are in good agreement with experiments. By comparing with the experimental enthalpy of β, we found that the 0 K total energy calculated from bcc Ti is incorrect. This problem can be solved by shifting the total energy of β down by 8 kJ mol-1 to match the experimental value. With the Gibbs energy calculated from the Debye model as a function of pressure and temperature, the phase transformation conditions of α→ω, α→β, and β→ω were identified. The predicted transition temperature between α and β at ambient pressure and the triple point are close to experiments. It was found that the entropy plays an important role in the ω→α and α→β transitions, and the thermal electronic contribution to the Gibbs energy cannot be neglected for studying the α→β transition. Our calculations also showed that zero-point energy is crucial to predict the transition pressure of Ti at low temperatures.