Density-Functional Tight-Binding Parameters for Bulk Zirconium: A Case Study for Repulsive Potentials

Hybrid Density Functional Tight Binding (DFTB)─Molecular Mechanics Approach for a Low-Cost Expansion of DFTB Applicability

The density functional-based tight binding (DFTB) method has seen a rise in adoption for materials modeling, as it offers significant improvement in scalability with accuracy comparable to the density functional theory (DFT) when good parameterizations exist. The cost reduction in DFTB compared to DFT is achieved by the pre-parameterization of the elements of the Hamiltonian matrix as well as the repulsion potential between all pairs of atoms. Parameterization for new systems with accuracies competitive with DFT in specific applications requires specialized manpower and computational resources. This prevents the application of the DFTB method to systems for which it was not parameterized. In this work, we explore an approach to address the problem of missing parameters of DFTB by modeling the interactions with missing Slater-Koster parameters with an interatomic interaction potential. When the distance between two atoms modeled at the force-field level is sufficiently large, the approach results in accurate structural and electronic properties. The resulting calculation is therefore a hybrid between DFTB and molecular mechanics, a pure DFTB for atoms with a complete set of interatomic parameterizations, and a mix between DFTB and molecular mechanics for atoms with a missing interatomic parameterization. The approach is expected to be particularly useful for hybrid materials and interfaces. The method is tested on the examples of 2D materials, mixed oxides, and a large-scale calculation of an oxide-oxide interface.

Development of Density-Functional Tight-Binding Parameters for the Molecular Dynamics Simulation of Zirconia, Yttria, and Yttria-Stabilized Zirconia

In this work, a set of density-functional tight-binding (DFTB) parameters for the Zr-Zr, Zr-O, Y-Y, Y-O, and Zr-Y interactions was developed for bulk and surface simulations of ZrO₂ (zirconia), Y₂O₃ (yttria), and yttria-stabilized zirconia (YSZ) materials. The parameterization lays the ground work for realistic simulations of zirconia-, yttria-, and YSZ-based electrolytes in solid oxide fuel cells and YSZ-based catalysts on long timescales and relevant size scales. The parameterization was validated for the zirconia and yttria polymorphs observed under standard conditions based on density functional theory calculations and experimental data. Additionally, we performed DFTB-based molecular dynamics (MD) simulations to compute structural and vibrational properties of these materials. The results show that the parameters can give a qualitatively correct phase ordering of zirconia, where the tetragonal phase is more stable than the cubic phase at a lower temperature. The lattice parameters are only slightly overestimated by 0.05-0.1 Å (2% error), still within the typical accuracy of first-principles methods. Additionally, the MD results confirm that zirconia and yttria phases are stable against transformations under standard conditions. The parameterization also predicts that vibrational spectra are within the range of 100-1000 cm^-1 for zirconia and 100-800 cm^-1 for yttria, which is in good agreement with predictions both from full quantum mechanics and a recently developed classical force field. To further demonstrate the advantage of the developed DFTB parameters in terms of computational resources, we conducted DFTB/MD simulations of the YSZ4 and YS12 models containing approximately 750 atoms.