Design and characterization of all two-dimensional fragile topological bands

Designing topological materials with specific topological indices is a complex inverse problem, traditionally tackled through manual, intuition-driven methods that are neither scalable nor efficient for exploring the vast space of possible material configurations. In this work, we develop an algorithm that leverages the covariance matrix adaptation evolution strategy to optimize the Fourier representation of the periodic functions shaping the designer material's characteristics. This includes mass profiles or dielectric tensors for phononic and photonic crystals, respectively, as much as synthetic potentials applicable to electronic and ultra-cold atomic systems. We demonstrate our methodology with a detailed characterization of a class of topological bands known as "fragile topological", showcasing the algorithm's capability to address both topological characteristics and spectral quality. This automation not only streamlines the design process but also significantly expands the potential for identifying and constructing high quality designer topological materials across the wide range of platforms, and is readily extendable to other setups, including higher-dimensional and non-linear systems.