Ab initio calculations for the investigation of crystal field in Ln complexes offers advantages over other approaches while permitting the application of crystal field analysis in a more insightful way than earlier phenomenological approaches. However, computational demands of ab initio approaches arise from the complex multiplet structure of many-electron atoms and the possibility of errors ˃ 20 % in comparison to experimentally determined crystal field energies.