Bridging constrained random-phase approximation and linear response theory for computing Hubbard parameters (2505.03698v1)

Abstract: The predictive accuracy of popular extensions to density-functional theory (DFT) such as DFT+U and DFT plus dynamical mean-field theory (DFT+DMFT) hinges on using realistic values for the screened Coulomb interaction U. Here, we present a systematic comparison of the two most widely used approaches to compute this parameter, linear response theory (LRT) and the constrained random-phase approximation (cRPA), using a unified framework based on the use of maximally localized Wannier functions to define the underlying basis sets. We demonstrate good quantitative agreement between LRT and cRPA in cases where the strongly interacting subspace corresponds to an isolated set of bands. Differences can be assigned to neglecting the response of the exchange-correlation potential in cRPA and the presence of additional screening channels in LRT. Moreover, we show that in cases with strong hybridization between interacting and screening subspaces, the application of cRPA becomes ambiguous and can lead to unrealistically small U values, while LRT remains well-behaved. Our work clarifies the relation between both methods, sheds light on their strengths and limitations, and emphasizes the importance of using a consistent set of Wannier orbitals to ensure transferability of U values between different implementations.