<p>The calculations based on the linear combination of atomic orbitals have been performed for the low-temperature phase of BaTiO₃ crystal. Structural and electronic properties, as well as phonon frequencies were obtained using hybrid PBE0 exchange-correlation functional. The calculated frequencies and total energies at different volumes have been used to determine the equation of state and thermal contribution to the Helmholtz free energy within the quasiharmonic approximation. For the first time, the bulk modulus, volume thermal expansion coefficient, heat capacity, and Gr&uuml;neisen parameters in BaTiO₃ rhombohedral phase have been estimated at zero pressure and temperatures form 0 to 200 K, based on the results of first-principles calculations. Empirical equation has been proposed to reproduce the temperature dependence of the calculated quantities. The agreement between the theoretical and experimental thermodynamic properties was found to be satisfactory.</p>