Norm-Conserving 4f-in-Core Pseudopotentials and Basis Sets Optimized for Trivalent Lanthanides (Ln = Ce–Lu)

-containing molecular complexes, aqueous solutions, and bulk solids to validate the accuracy, reliability, and efficiency of the optimized 4f-in-core GTH pseudopotentials and basis sets. The 4f-in-core GTH pseudopotentials successfully replicate the main features of lanthanide structural chemistry and reaction energetics, particularly for nonredox reactions. The chemical bonding features and solvation shells, hydrolysis energetics, acidity constants, and solid-state properties of selected lanthanide systems are also discussed in detail by utilizing these new 4f-in-core GTH pseudopotentials. This work bridges the idea of keeping highly localized 4f electrons in the atomic core and efficient pseudopotential formalism of GTH, thus providing a highly efficient approach for studying lanthanide chemistry in multi-scale modeling of constituent-wise and structurally complicated systems, including electronic structures of the condensed phase and first-principles molecular dynamics simulations.