  We develop an efficient numerical algorithm for the identification of a large
number of saddle points of the potential energy function of Lennard- Jones
clusters. Knowledge of the saddle points allows us to find many thousand
adjacent minima of clusters containing up to 80 argon atoms and to locate many
pairs of minima with the right characteristics to form two-level systems (TLS).
The true TLS are singled out by calculating the ground-state tunneling
splitting. The entropic contribution to all barriers is evaluated and
discussed.
