Lanthanide trihalide molecules LnX3 (X = F, Cl, Br, I) were quantum chemically investigated, in particular detail for Ln = Lu (lutetium). We applied density functional theory (DFT) at the nonrelativistic and scalar and SO-coupled relativistic levels, and also the ab initio coupled cluster approach. The chemically active electron shells of the lanthanide atoms comprise the 5d and 6s (and 6p) valence atomic orbitals (AO) and also the filled inner 4f semivalence and outer 5p semicore shells. Four different frozen-core approximations for Lu were compared: the (1s(2) -4d(10) ) [Pd] medium core, the [Pd+5s(2) 5p(6) = Xe] and [Pd+4f(14) ] large cores, and the [Pd+4f(14) +5s(2) 5p(6) ] very large core. The errors of LuX bonding are more serious on freezing the 5p(6) shell than the 4f(14) shell, more serious upon core-freezing than on the effective-core-potential approximation. The LnX distances correlate linearly with the AO radii of the ionic outer shells, Ln(3+) -5p(6) and X(-) -np(6) , characteristic for dominantly ionic Ln(3+) -X(-) binding. The heavier halogen atoms also bind covalently with the Ln-5d shell. Scalar relativistic effects contract and destabilize the LuX bonds, spin orbit coupling hardly affects the geometries but the bond energies, owing to SO effects in the free atoms. The relativistic changes of bond energy BE, bond length Re , bond force k, and bond stretching frequency vs do not follow the simple rules of Badger and Gordy (Re ∼BE∼k∼vs ). The so-called degeneracy-driven covalence, meaning strong mixing of accidentally near-degenerate, nearly nonoverlapping AOs without BE contribution is critically discussed.