Search

The $f$-electron materials have many unique properties under pressure, thus of great interest in high-pressure physics and related industrial fields. However, the $f$-electrons pose a substantial challenge to simulations since the electron correlation effects. In this work, we present a first-principles calculation scheme for the equations of state (EoS) of $f$-electron materials. The self-consistent pressure-dependent on-site Coulomb correction is performed based on our recently developed doubly screened Coulomb correction approach. We investigate the zero-temperature EoS over a wide range of pressures and the phase stabilities of four prototypical lanthanide and actinide metals, Pr, Eu, Th and U. The simulated compressive properties are in better agreement with the experimental data than those obtained by conventional density functional theory (DFT) and fixed-parameter DFT+$U$ approaches. The pressure-induced phase transitions can also be well described.

Strongly correlated systems containing d/f-electrons present a challenge to conventional density functional theory (DFT), such as the widely used local density approximation (LDA) or generalized gradient approximation (GGA). In this work, we developed a doubly screened Coulomb correction (DSCC) approach to perform on-site Coulomb interaction correction for strongly correlated materials. The on-site Coulomb interaction between localized d/f-electrons is determined from a model dielectric function that includes both the static dielectric and the Thomas-Fermi screening. All parameters of the dielectric model are efficiently obtained from self-consistent calculations. We applied DSCC to simulate the electronic and magnetic properties of typical 3d, 4f and 5f strongly correlated systems. The results show that the accuracy of DSCC is comparable to hybrid functionals, but an order of magnitude faster. In addition, DSCC can reflect the difference in the Coulomb interaction of the same element between metallic and insulating situations, similar to the popular but computationally expensive constrained random phase approximation (cRPA) approach. This feature suggests that DSCC is also a promising method for simulating Coulomb interaction parameters., Comment: 11 pages, 5 figures, 2 tables

As correlation strength has a key influence on the simulation of strongly correlated materials, many approaches have been proposed to obtain the parameter using first-principles calculations. However, the comparison of the different Coulomb strengths obtained using these approaches and an investigation of the mechanisms behind them are still needed. Taking lanthanide metals as an example, we research the factors that affect the effective Coulomb interaction strength, $U_{\rm eff}$, by local screened Coulomb correction (LSCC), linear response (LR) and constrained random-phase approximation (cRPA) in VASP. The $U^{\rm LSCC}_{\rm eff}$ value increases from 4.75 eV to 7.78 eV, $U^{\rm LR}_{\rm eff}$ is almost stable at about 6.0 eV (except for Eu, Er and Lu), and $U^{\rm cRPA}_{\rm eff}$ shows a two-stage decreasing trend in both light and heavy lanthanides. To investigate these differences, we established a scheme to analyze coexistence and competition between the orbital localization and the screening effect. We find that LSCC and cRPA are dominated by the orbital localization and the screening effect, respectively, whereas LR shows a balance of the competition between the two factors. Additionally, the performance of these approaches is influenced by different starting points from PBE and PBE+$U$, especially for cRPA. Our results provide useful knowledge for understanding the $U_{\rm eff}$ of lanthanide materials, and similar analyses can also be used in the research of other correlation strength simulation approaches., Comment: 10 pages, 8 figures, 1 table

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., Comment: 20 pages, 8 figures