Experiments on bilayer graphene unveiled a fascinating realization of stacking disorder where triangular domains with well-defined Bernal stacking are delimited by a hexagonal network of strain solitons. Here we show by means of numerical simulations that this is a consequence of a structural transformation of the moiré pattern inherent of twisted bilayer graphene taking place at twist angles θ below a crossover angle θ = 1.2 • . The transformation is governed by the interplay between the interlayer van der Waals interaction and the in-plane strain field, and is revealed by a change in the functional form of the twist energy density. This transformation unveils an electronic regime characteristic of vanishing twist angles in which the charge density converges, though not uniformly, to that of ideal bilayer graphene with Bernal stacking. On the other hand, the stacking domain boundaries form a distinct charge density pattern that provides the STM signature of the hexagonal solitonic network. * oleg.yazyev@epfl.ch arXiv:1711.08647v2 [cond-mat.mes-hall] 11 Apr 2018 Bilayer graphene (BLG) shares many of the properties of monolayer graphene while also showing a number of pronounced differences. For instance, its equilibrium structural configuration reveals the massive nature of its charge carriers [1], the possibility of inducing a tunable band gap by applying a transverse electric field [2-4] and quantum Hall valley ferromagnetism [5]. These properties are a result of the coupling between the two layers.In order to describe the atomic structure of bilayer graphene the relative position of the two layers has to be defined. In many situations it is sufficient to specify a unique interlayer displacement vector that defines the stacking configuration. As a general property of graphitic structures, the low-energy configuration is represented by the Bernal stacking [6,7]. However, the stacking configuration is not immune to disorder which can manifest, for example, in boundaries that connect two domains with energetically degenerate yet topologically inequivalent stacking configurations, AB and BA [8][9][10][11]. Such stacking domain boundaries are realized by strain solitons, which are segments with a characteristic width where the strain that arises from interfacing two inequivalent stacking domains is confined.Recent studies have shown that strain solitons can be displaced by the action of a scanning tunneling microscope tip, but do not vanish due to their topological nature [8,12]. From the theoretical point of view, the two-dimensional extension of the Frenkel-Kontorova model predicts the emergence of strain solitons with a typical width of a few nanometers [13] while their density is defined by the twist angle.In other situations the stacking configuration cannot be uniquely defined on the whole surface of the sample since the two layers cannot be superimposed by a rigid in-plane shift. This is the case of twisted bilayer graphene where one layer is rotated relative to another, a system that has been widely reported ...