  We report the initial results of deep eROSITA monitoring of the "magnificent
seven" isolated neutron stars (INSs). Thanks to a combination of high count
statistics and good energy resolution, the eROSITA datasets unveil the
increasingly complex energy distribution of these presumably simple thermal
emitters. For three targets, we report the detection of multiple (in some
cases, phase-dependent) spectral absorption features and deviations from the
dominant thermal continuum. Unexpected long-term changes of spectral state and
timing behaviour have additionally been observed for two INSs. The results pose
challenging theoretical questions on the nature of the variations and
absorption features and ultimately impact the modeling of the atmosphere and
cooling of highly magnetised neutron stars.
