“…According to theḡ[= (g + 2g ⊥ )/3 ≈ 2.34] of the experimental g factors of Nd 3+ in LiYF 4 [6], it can be attributed to the lowest doublet Γ 6 , whose average valueḡ would be about 2.67 for a 4f 3 ion [7,8]. In [7,8] only the contributions to the EPR parameters from the first-order perturbation terms were included. However, besides the lowest Γ 6 , the other 10 irreducible representations Γ x (i.e., five Γ 6 and five Γ 7 ) due to the tetragonal splitting of the ground 4 I 9/2 and the first excited 4 I 11/2 levels would mix with the lowest Γ 6 via crystal-fieldĤ CF and orbital angular momentumĴ (or hyperfine structure equivalent operatorN) interactions and lead to the second-order perturbation contributions to g factors (or hyperfine structure constants), as pointed out in [9,10].…”