pub fn num_xtb_basis_functions(z: u8) -> usize
Count basis functions: 1 (s-only for H), 4 (s+p), or 9 (s+p+d for TM).