Noncoplanar and chiral spin states on the way towards Nel ordering in fullerene Heisenberg models

Using a high -accuracy variational Monte Carlo approach based on group -convolutional neural networks, we obtain the symmetry -resolved low -energy spectrum of the spin -1/2 Heisenberg model on several highly symmetric fullerene geometries, including the famous C60 buckminsterfullerene. We argue that as the degree of frustration is lowered in large fullerenes, they display characteristic features of incipient magnetic ordering: Correlation functions show high -intensity Bragg peaks consistent with Neel-like ordering, while the low -energy spectrum is organized into a tower of states. Competition with frustration, however, turns the simple Neel order into a noncoplanar one. Remarkably, we find and predict chiral incipient ordering in a large number of fullerene structures.