I'm attempting to calculate a crystal containing praseodymium and have a few inquiries. Is it possible to compute this compound accurately? The program documentation regarding basis sets and pseudopotentials does not specify anything for lanthanides or actinides, and previous forum posts suggested that calculations for Pr and other elements were not supported. Are these elements now supported? With version 2023, I can perform calculations using GGA with LCAO without encountering any errors, and the input parameters seem good. However, attempting to use PAW shows me that no pseudopotentials are available for Pr, and LDA yields an error. I'm uncertain whether the parameters for describing Pr atoms in the GGA calculation are correct.

Additionally, I aim to incorporate spin-orbit coupling (SOC) and make some corrections to the bandgap in the calculation. The DFT-1/2 correction is not defined for Pr. Does this mean the parameter is not automatically defined, and I can assign it a value, or does it indicate that this correction is not applicable for Pr? Can those corrections be applied to systems with Pr?

Thank you for your assistance.