We model the PCS using the LFT Hamiltonian and the point-dipole approximation (PDA). We employ parameters obtained from ab initio ligand field theory (AILFT). For the finite-field approach, we diagonalize the exact LFT Hamiltonian in the presence of the external field B 0 for each orientation and compute the electronically-averaged induced field according to Eq. ( ) by using the operator