Track synthesis, materials and catalysis as it is published.
Search papers, labs, and topics across Lattice.
This paper addresses the limitations of standard density-functional tight-binding (DFTB) for vanadium oxides by introducing a parameter-free reference Hamiltonian derived from superposed neutral atoms and the full PBE exchange-correlation potential. The proposed model accurately reproduces PBE energies and charge responses for VO2, achieving sub-meV accuracy for energy differences and capturing the self-consistent potential shift on the V 3d shell. This approach overcomes significant deviations in band-index centroid separation caused by DFTB approximations, offering a more reliable description of these materials.
A parameter-free tight-binding model achieves sub-meV accuracy for VO2 energy differences and captures V 3d shell self-consistent shifts, surpassing standard DFTB limitations.
Density-functional tight binding has yet to yield a transferable description of vanadium oxides; for the standard second-order construction and the tested cells, we demonstrate that two of its approximations deviate by electron volts and offset each other in the band-index centroid separation between the V 3d and O 2p bands: the neutral two-centre reference omits electron-volt environmental terms, whereas the Mulliken charge model assigns the compact V 3d shell less than half of the self-consistent potential shift observed in PBE. Instead, we build a tight-binding model where the reference Hamiltonian derives from superposed neutral atoms, the full PBE exchange-correlation potential and frozen-core orthogonality, without adjustable terms, and where the orbital-resolved charge response employs three screening factors obtained from PBE linear response and an unscreened Coulomb range. Against the primary PAW reference, the reference energy at fixed geometry and the stated numerical settings matches the plane-wave Harris-Foulkes energy differences of VO2, exhibiting signed deviations of +2.6 meV/atom for M1 relative to rutile (an energy difference of 92 meV/atom) and +6.2 and -5.2 meV/atom under 3% volume compression and expansion; meanwhile, the charge response captures 21 of 22 fitted and all 8 held-out matched PBE responses within tolerance, providing the V 3d shell with nearly the same self-consistent shift as PBE (-3.89 versus -3.92 eV in rutile). The convexity of the spinless model at fixed geometry prevents it from spontaneously differentiating symmetry-equivalent sites under the symmetries respected by its implemented operators and sampling, so within this model such differentiation in the Magnéli phases must arise from electron-lattice coupling or an explicit spin or correlation term.