In this study we report about the aromaticity of the prototypical [(HtAc)(3)(mu(2)-H)(6)], [(HtTh)(3)(mu(2)-H)(6)](+), and [(HtPa)(3)(mu(2)-H)(6)] clusters via two magnetic criteria: nucleus-independent chemical shifts (NICS) and the magnetically induced current density. All-electron density functional theory calculations were carried out using the two-component zeroth-order regular approach and the four-component Dirac-Coulomb Hamiltonian, including scalar and spin-orbit relativistic effects. Four-component current density maps and the integration of induced ring-current susceptibilities clearly show that the clusters [(HtAc)(3)(mu(2)-H)(6)] and [(HtTh)(3)(mu(2)-H)(6)](+) are non-aromatic whereas [(HtPa)(3)(mu(2)-H)(6)] is anti-aromatic. However, for the thorium cluster we find a discrepancy between the current density plots and the classification through the NICS index. Our results also demonstrate the increasing influence of f orbitals, on bonding and magnetic properties, with increasing atomic number in these clusters. We think that the enhanced electron mobility in [(HtPa)(3)(mu(2)-H)(6)] is due the significant 5f character of its valence shell. Also the participation off orbitals in bonding is the reason why the protactinium cluster has the shortest bond lengths of the three clusters. This study provides another example showing that the magnetically induced current density approach can give more reliable results than the NICS index. (C) 2011 American Institute of Physics. [doi:10.1063/1.3632056]