The nature of chemical bonds in heavy main-group diatomics is discussed from the viewpoint of effective bond orders, which are computed from spin-orbit wave functions resulting from contracted spin-orbit configuration interaction calculations. The reliability of the relativistic correlated wave functions obtained in such two-step spinorbit coupling frameworks is assessed by benchmark studies of the spectroscopic constants with respect to either experimental data, or state-of-the-art fully relativistic correlated calculations. The I 2 , At 2 , IO + , and AtO + species are considered, and differences and similarities between the astatine and iodine elements are highlighted.In particular, we demonstrate that spin-orbit coupling weakens the covalent character of the bond in At 2 even more than electron correlation, making the consideration of spin-orbit coupling compulsory for discussing chemical bonding in heavy (6p) main group element systems.