Magic angles are a hot topic in condensed matter physics: when two sheets of graphene are twisted by those angles, the resulting material is superconducting and the so-called energy bands are flat. In 2011, Bistritzer and MacDonald proposed a model that is experimentally very accurate in predicting magic angles. In this talk, I will introduce some recent mathematical progress on the Bistritzer--MacDonald's model, including the mathematical characterization of magic angles and flat bands, and generic existence of Dirac cones. I will also discuss some new mathematical discoveries in twisted multilayer graphene.