This work introduces a novel, fully robust and highly-scalable, $h$-adaptive aggregated unfitted finite element method for large-scale interface elliptic problems. The new is based on recent distributed-memory implementation of the atop highly-scalable Cartesian forest-of-trees mesh engine. It follows classical approach weakly coupling nonmatching discretisations at to model internal discontinu...