<p>Phonon frequencies at the &Gamma;, X, M, and R points of the Brillouin zone (BZ) in the cubic phase of Ba(Ti,Zr,Hf)₃ were calculated by the frozen phonon method using density functional theory with the hybrid exchange correlation functional PBE0. The calculations use linear combination of atomic orbitals basis functions as implemented in the crystal09 computer code. The Powell algorithm was applied for basis-set optimization. In agreement with experimental observations, the structural instability via soft modes was found only in BaTiO_3. A good quantitative agreement was found between the theoretical and experimental phonon frequency predictions in BaTiO₃ and BaZrO₃. It is concluded that the hybrid PBE0 functional is able to predict correctly the structural stability and phonon properties in perovskites under consideration. The symmetry of phonons for nonzero k vector (X, M, and R points of the BZ) is given for each calculated phonon frequency, and it will be useful for future comparison of the calculated results with the experimental data on neutron scattering (such data are absent at present).</p>