Ab initio Phase Diagram of Ta2O5

Tantalum pentoxide (Ta2O5) is a polymorphic wide-bandgap semiconductor with outstanding dielectric properties and widespread use in optical and electronic technologies. Its rich structural diversity, arising from multiple polymorphs accessible under different synthesis conditions, has made Ta2O5 a long-standing subject of interest. However, a unified understanding of the thermodynamic stability and phase transitions of its polymorphs across pressure-temperature (P-T) space has remained elusive. Here, using first-principles calculations, we map the thermodynamic landscape of Ta2O5 and establish a comprehensive P-T phase diagram together with a phase-stability hierarchy. We find that Gamma-Ta2O5 and B-Ta2O5 dominate the phase diagram over a broad range of P-T conditions: Gamma-Ta2O5 is stabilized at low pressures, while B-Ta2O5 becomes thermodynamically favored at higher pressures up to ~ 60 GPa, beyond which Y-Ta2O5 emerges as the most stable phase. Crucially, nuclear quantum effects (NQEs) are shown to play a significant role in determining relative phase stability, contributing substantially to the Gibbs free energy and altering phase boundaries. A re-entrant phase transition between Gamma and B-Ta2O5 is predicted near ~ 2 GPa, revealing unexpected complexity in the phase behavior of this oxide. More generally, we identify a characteristic temperature (T_0), at which zero-point and thermal phonon contributions to the free energy become comparable, and show that T_0 is approximately one-third of the Debye temperature. This relationship provides a simple, physically transparent criterion for assessing the importance of NQEs in phase stability, with implications extending beyond Ta2O5 to a broad class of complex oxides.

💡 Research Summary

This paper presents a comprehensive first‑principles investigation of the thermodynamic stability and phase transitions of tantalum pentoxide (Ta₂O₅) across a wide pressure–temperature (P‑T) range. Using density‑functional theory (DFT) with the PBE‑GGA functional and projector‑augmented wave (PAW) potentials as implemented in VASP, the authors fully optimized ten known polymorphs of Ta₂O₅ (γ, γ₁, B, λ, LSR, δ, β_AL, β_R, Z, and Y). Structural relaxations were performed with high‑quality k‑point meshes and a plane‑wave cutoff of 600 eV. The resulting energy‑volume data were fitted to a Birch‑Murnaghan equation of state to obtain equilibrium volumes and bulk moduli as functions of pressure.

Phonon spectra and vibrational free energies were calculated via density‑functional perturbation theory (DFPT) and the PHONOPY package. By combining static electronic energies, phonon contributions, and the PV term, the Gibbs free energy G(P,T) for each polymorph was obtained. This allowed the construction of a full P‑T phase diagram by selecting, at each (P,T) point, the phase with the lowest G. The authors also evaluated the role of nuclear quantum effects (NQEs) by explicitly separating zero‑point energy (ZPE) from thermal phonon contributions.

Key findings include: (1) At ambient pressure, the γ‑phase is the most stable, followed closely by the newly identified γ₁‑phase. (2) A re‑entrant phase transition is predicted near 2 GPa, where the B‑phase becomes more stable than γ, but γ regains stability at slightly higher pressures before B dominates again. (3) Between roughly 6 GPa and 60 GPa, the B‑phase is the thermodynamic ground state, while the λ‑phase becomes the second‑most stable polymorph for pressures above ~4 GPa up to ~32 GPa. (4) Above ~60 GPa, the high‑pressure Y‑phase, characterized by TaO₁₀ polyhedra and the highest Ta–O coordination reported for this system, becomes the most stable structure. (5) Nuclear quantum effects contribute significantly to the Gibbs free energy, especially at low temperatures, shifting phase boundaries by up to ~1 GPa. The ZPE accounts for 5–10 % of the total free energy in the relevant pressure range.

A particularly insightful result is the identification of a characteristic temperature T₀ at which the zero‑point vibrational energy equals the thermal phonon energy. The authors find that T₀ is approximately one‑third of the Debye temperature (Θ_D) for each polymorph, i.e., T₀ ≈ Θ_D/3. This simple proportionality provides a practical criterion for estimating the importance of NQEs in complex oxides without performing full phonon calculations.

The paper also discusses the structural chemistry of the polymorphs: γ and γ₁ consist solely of distorted TaO₆ octahedra; B and λ contain edge‑sharing octahedra; LSR combines octahedra with pentagonal bipyramids; δ includes both octahedra and hexagonal bipyramids; Z is built entirely from pentagonal bipyramids; and Y features unprecedented TaO₁₀ polyhedra. The transition from octahedral‑dominated networks to higher‑coordination polyhedra under pressure explains the observed stability sequence.

Overall, the study delivers the first complete ab‑initio P‑T phase diagram for Ta₂O₅, elucidates the impact of nuclear quantum effects on phase stability, and proposes a universal T₀–Θ_D relationship that can be applied to other complex oxides. These results provide valuable guidance for high‑pressure synthesis, thin‑film deposition, and the design of Ta₂O₅‑based electronic and photonic devices, such as high‑k dielectrics, UV photodetectors, and resistive‑switching memory elements.