[Paper Review] Hybrid density functional theory LCAO calculations on phonons in Ba (Ti,Zr,Hf) O3

Robert A Evaestov|arXiv (Cornell University)|Oct 6, 2010

Ferroelectric and Piezoelectric Materials3 citations

TL;DR

This study employs hybrid density functional theory with LCAO basis sets in the CRYSTAL09 code to calculate phonon frequencies in cubic Ba(Ti,Zr,Hf)O3 using the frozen phonon method. It finds that only BaTiO3 exhibits structural instability via soft modes, with excellent agreement between theoretical and experimental phonon frequencies in BaTiO3 and BaZrO3, confirming the PBE0 functional's accuracy for perovskite stability and vibrational properties.

ABSTRACT

Phonon frequencies at Γ,X,M,R-points of Brilloin zone in cubic phase of Ba(Ti,Zr,Hf)O3 were first time calculated by frozen phonon method using density functional theory (DFT) with hybrid exchange correlation functional PBE0. The calculations use linear combination of atomic orbitals (LCAO) basis functions as implemented in CRYSTAL09 computer code. The Powell algorithm was applied for basis set optimization. In agreement with the experimental observations the structural instability via soft modes was found only in BaTiO3. A good quantitative agreement was found between the theoretical and experimental phonon frequency predictions in BaTiO3 and BaZrO3. It is concluded that the hybrid PBE0 functional is able to predict correctly the structural stability and phonon properties in perovskites.

Motivation & Objective

• To investigate the phonon behavior and structural stability of cubic Ba(Ti,Zr,Hf)O3 perovskites using advanced quantum mechanical methods.
• To assess the predictive capability of the hybrid PBE0 functional in describing phonon frequencies and soft modes in complex perovskite oxides.
• To evaluate the performance of the LCAO method with CRYSTAL09 in modeling vibrational properties in transition metal oxide perovskites.
• To compare theoretical phonon frequencies with experimental data to validate the computational approach for complex oxides.

Proposed method

• The frozen phonon method was employed to calculate phonon frequencies at high-symmetry points (Γ, X, M, R) in the Brillouin zone of cubic Ba(Ti,Zr,Hf)O3.
• Density functional theory (DFT) with the hybrid exchange-correlation functional PBE0 was used to improve accuracy over standard GGA functionals.
• Linear combination of atomic orbitals (LCAO) basis sets were implemented via the CRYSTAL09 computational code for electronic structure calculations.
• The Powell algorithm was applied to optimize the basis set, enhancing the convergence and reliability of the electronic and vibrational properties.
• Phonon frequencies were computed by finite displacement of atoms in the supercell, and the Hessian matrix was derived from energy derivatives.
• The calculations were performed for all three end-member compositions: BaTiO3, BaZrO3, and BaHfO3 to compare stability and vibrational responses.

Experimental results

Research questions

• RQ1 Does the PBE0 functional correctly predict the structural stability of cubic Ba(Ti,Zr,Hf)O3 perovskites through phonon analysis?
• RQ2 How accurately do LCAO-based DFT calculations reproduce experimental phonon frequencies in BaTiO3 and BaZrO3?
• RQ3 Which perovskite in the Ba(Ti,Zr,Hf)O3 system exhibits soft mode-driven structural instability?
• RQ4 To what extent does the choice of exchange-correlation functional (PBE0) improve phonon frequency predictions compared to standard DFT functionals?
• RQ5 Can the LCAO method in CRYSTAL09 reliably model vibrational modes in complex perovskite oxides with multiple cations?

Key findings

• Only BaTiO3 exhibits a structural instability via soft modes, consistent with experimental observations of ferroelectric transition.
• Theoretical phonon frequencies in BaTiO3 show excellent quantitative agreement with experimental data, validating the PBE0 functional for this system.
• A similarly good agreement was found between theory and experiment for phonon frequencies in BaZrO3, indicating robust predictive power.
• The hybrid PBE0 functional successfully captures the correct vibrational and stability trends across the Ba(Ti,Zr,Hf)O3 solid solution series.
• The LCAO method with basis set optimization via the Powell algorithm yields reliable and converged results for phonon dispersion in perovskites.
• The study confirms that hybrid DFT with PBE0 and LCAO is a suitable and accurate approach for predicting phonon properties and structural stability in complex perovskite oxides.