The aim is to find the maximum size of a set of mutually ske lines on a non
singular Hermitian surface in PG (3, q) for various values of q. For q = 9
such extremal sets are intricate combinatorial structures intimately connec
ted ith hemisystems, subreguli, and commuting null polarities. It turns out
they are also closely related to the classical quartic surface of Kummer.
Some bounds and examples are also given in the general case.