We analyze a description of twisted graphene bilayers, that incorporates the deformation of the layers using state of the art interlayer atomic potentials, and a modification of the hopping parameters between layers in the light of the classic Slonczewski-Weiss-McClure parametrisation. We obtain narrow bands in all cases, but that their nature can be rather different. We will show how to describe the results by equivalent continuum models. Even though such models can be constructed, their complexity can vary, requiring many coupling parameters to be included, and the full in-layer dispersion must be taken into account. The combination of all these effects will have a large impact on the wave functions of the flat bands, and that modifications in details of the underlying models can lead to significant changes. A robust conclusion is that the natural strength of the interlayer couplings is higher than usually assumed, leading to shifts in the definition of the magic angles. The structure at the edges of the narrow bands, at the Γ point of the Brillouin Zone is also strongly dependent on parametrization. As a result, the existence, and size, of band gaps between the flat bands and the neighboring ones are changed. Hence, the definition of Wannier functions, and descriptions based on local interactions are strongly dependent on the description of the model at the atomic scale. * Francisco. Guinea@imdea.org † Niels.Walet@manchester.ac.uk; https://www.research.manchester.ac.uk/portal/niels.walet.html AA AB FIG. 1. Examples of a graphene bilayer in (approximate) AA and AB alignment. arXiv:1903.00364v2 [cond-mat.str-el]