The underestimation of high pressure in DFT+$U$ simulation for the wide range cold-pressure of lanthanide metals

Abstract

Density functional theory plus $U$ (DFT+$U$) is one of the most efficient first-principles methods to simulate the cold pressure properties of strongly-correlated materials. However, the applicability of DFT+$U$ at ultra-high pressure is not sufficiently studied, especially in the widely-used augmented schemes [such as projector augmented wave (PAW) and linearized augmented plane wave (LAPW)]. This work has systematically investigated the performance of DFT+$U$ in PAW and LAPW at the pressure up to several hundred GPa for the lanthanide metals, which is a typical strongly-correlated series. We found DFT+$U$ simulation in PAW exhibits an unphysical underestimating of force at high pressure. By delicate analysis and comparison with local-orbital-independent hybrid functional results, we have demonstrated that this unphysical behavior is related to a normalization problem on the local density matrix caused by the overlap of local orbitals in PAW under high pressure. Additionally, we observed a slight softening of force in DFT+$U$ in heavy lanthanides (Tm, Yb and Lu) at high pressure comparing with the DFT results without the influence of local orbital overlap, and it might be related to the enhancement of bonding effect in correlation correction methods at high pressure. Our work reveals the underestimating of high pressure in DFT+$U$ simulation, analyses two sources of this unusual behavior and proposes their mechanism. Most importantly, our investigation highlights the breakdown of DFT+$U$ for high pressure simulation in VASP package based on PAW framework.

arXiv comment20 pages, 8 figures

References in corpus (8)

• Implementation of the Projector Augmented Wave LDA+U Method: Application to the Electronic Structure of NiO
• Implementation of the LDA+U method using the full potential linearized augmented plane wave basis
• The f-electron challenge: localized and itinerant states in lanthanide oxides united by GW@LDA+U
• Ground state properties and high pressure behavior of plutonium dioxide: Systematic density functional calculations
• Structural stability and thermodynamics of CrN magnetic phases from ab initio and experiment
• Structural behavior of uranium dioxide under pressure by LSDA+U calculations
• Equation of state of paramagnetic CrN from ab initio molecular dynamics
• Projector self-consistent DFT+U using non-orthogonal generalized Wannier functions