“…microscopic) parameters for our calculations is selected as (in cm −1 ): ∆ 1 = ∆ 2 = 22500, ∆ 3 = 25000, ∆ 4 = 21500, the SOC constant λ in the range 50 to 110, and the SSC constant ρ = 0.18 [28]. Additionally, for tetragonal symmetry the mixing coefficients [6,15] satisfy the relations: q = 1, s = 0. Note that the SOC constant λ in crystal is reduced from the free ion value λ 0 : λ = k × λ 0 , k ≤ 1, due to the orbital (covalency) reduction [16,17,19].…”