Orbital Magnetism in Ensembles of Ballistic Billiards

We calculate the magnetic response of ensembles of small two-dimensional structures at finite temperatures. Using semiclassical methods and numerical calculation we demonstrate that only short classical trajectories are relevant. The magnetic susceptibility is enhanced in regular systems, where these trajectories appear in families. For ensembles of squares we obtain a large paramagnetic susceptibility, in good agreement with recent measurements in the ballistic regime.