1. Introduction
During the last years, the possibility of tailoring heterostructures by piling up several two-dimensional materials has become a reality [
1,
2]. Such structures are not only important for technological applications [
1] but also serve as a playground to observe new and exotic phases of matter [
3,
4]. Among such phases, we can cite the observation of the Hofstadter butterfly fractal spectrum and its associated topological states. Such experiment was made by using a moiré pattern of graphene over graphene [
5]. More recently, a phase-diagram akin to those seen in complex high-
superconductor materials has been observed in graphene over twisted graphene at certain magical rotation angles [
6]. Moreover, recent advances in controlling the rotation angle of graphene over graphene [
7] allows using it as an effective knob for tuning their electronic properties [
8]. Although these results are not fully understood at the moment, recent theoretical advances are based either on Hubbard Hamiltonians [
9] or Dirac Hamiltonians with an emergent effective potential for electrons [
10]. In any case, there are many open questions concerning the electronic and optical properties of 2D heterostructures which require adequate tools to describe their structure [
3,
11].
A typical example of the involved complexity is the actual determination of the lattice parameter of the moiré superlattice. For example, in graphene on Ir(1 1 1), obtained via a pyrolytic cleavage of ethylene [
12,
13], there is a slight lattice mismatch between the graphene and Ir(1 1 1) lattices, resulting in a moiré pattern with a repeat distance of 2.53 nm [
12]. However, refinements indicated that this was just an excellent approximation, since the actual structure comprises three beatings instead of two [
14].
Moreover, until recently most work has dealt with the commensurate case, that is, when two layers form a supercell, what in materials science is known as a coincidence site lattice. Consider that one of the layers
is hexagonal and spanned by
and the other
by
and
can be obtained from
by means of a rotation through an angle
. Then, a discrete
coincidence site lattice,
, exists only when
satisfies,
for integers
m and
n. In the materials science and crystallographic literature, an equivalent condition,
is used, where
p and
q are integers [
15]. Both can be derived using the methods described first by Ranganathan [
16], however the latter form involving half angles is in better agreement with the spinorial (Rodrigues vector, quaternions, and Cayley–Klein parameters) character of the rotation involved [
15,
17].
Researchers in the 2D matter community field frequently refer to moiré patterns, but no one has recalled that moiré patterns in (periodic or not) structures are related to Bollmann’s O-lattice concept [
18]. Here, we describe this tool. In fact, the need for such tools is becoming evident as, for example, very recently it has been claimed that there is an incommensurate heterostructure analog to a quasicrystal [
11], obtained by growing a quasicrystalline 30
twisted bilayer graphene (30
-tBLG), stabilized by a Pt(111) substrate [
11]. Beside this experiment, the 30
quasicrystalline twisted bilayer graphene was simultaneously discovered by Ahn, Sung Joon, et al. [
19]. There are now available some theoretical works to model such system [
20,
21]. The observed quasiperiodic order is predicted to exhibit quantum oscillations associated to exotic spiral Fermi surfaces [
22].
The aim of this paper is twofold: first, to clarify and present the tools necessary to describe such kind of heterostructures and, second, to develop a cut and projection approach to obtain the best fit lattice between two rotated graphene structures. One can claim that no such concepts are needed as a moiré pattern seems to be just a superposition of two lattices. However, the superposition of the lattices results in a modulation of the graphene sheet as strain needs to be minimized [
3,
10]. Thus, there is a deformation strain field which modulates the equilibrium position of atoms on each graphene layer [
3]. This follows the spirit of one of the oldest known quasiperiodic Hamiltonians: the Frenkel–Kontorova model [
23]. Moreover, all particles and quasiparticles, such as electrons or holes, feel an effective potential which is the sum of both deformed and modulated lattices [
10]. Such deformation field has a paramount importance to the observed superconductivity at magic angles [
8]. Thus, our approach here is in the spirit of the Frenkel–Kontorova scenario: an effective two-dimensional model.
The structure of this paper is as follows. In
Section 2, we provide some general remarks concerning the crystallography of graphene, including coincidence lattices and Bollmann’s O-lattice and point out its presence in the bilayers. The cut and projection method is applied to the generation of bilayers in
Section 3. A discussion and our conclusions are presented in
Section 4.
3. Cut and Projection Approach
A powerful tool to describe modulated crystals and quasicrystals is the cut and projection method. In this section, we adapt this method to describe rotated layers. For the sake of completeness, let us briefly review such method restricted to quasiperiodic tilings of the plane. Under certain conditions (see below), a quasiperiodic tiling of the plane generated by linear integer combinations of the vectors , , can be obtained by projecting points from a N-dimensional hypercubic lattice in the following way. Let be the hypercubic lattice generated by the standard basis and let be the unitary hypercube. Let be a two-dimensional subspace of that does not contain any point of the lattice except the origin. An open strip is defined and the projection onto of all the points inside S gives a quasiperiodic tiling of the plane.
This projection formalism works provided that there exists a projector such that
The existence of
is guaranteed if the set
forms a
eutactic star [
25,
26]. The orthogonal projection of a set of orthogonal vectors in
onto a two-dimensional space is called a eutactic star; if the vectors in
are orthonormal, the projection is called a
normalized eutactic star. In other words, if the set
forms a normalized eutactic star, then there exists a projector
such that Equation (
3) is fulfilled and actually:
A necessary and sufficient condition for a star to be eutactic is due to H. Hadwiger [
27]: a star
in
is a normalized eutactic if and only if
is fulfilled for all
.
A more practical form of the eutacticity criterion is obtained if the structure matrix
S is introduced. In this case,
S is the matrix whose columns are the components of the vectors
, with respect to a given fixed orthonormal basis of
. The matrix form of Hadwiger’s theorem sates that the star represented by
S is normalized eutactic if and only if
where
I is the
unit matrix.
To summarize, a quasiperiodic tiling of the plane generated by linear integer combinations of the vectors , , can be obtained by cut and projection from if the set of vectors form a normalized eutactic star.
To adapt the cut and projection method to the case of graphene bilayers, consider first the case of two lattices in the plane,
and
, rotated in between by a given angle
. Let
and
be generated by
and
, respectively, and given by:
where
and
.
As shown in
Figure 2, if
, both vector sets collapse and three vectors pointing to the vertices of a triangle are obtained. If
, six vectors pointing to six vectors of a dodecagon are obtained instead. Consequently, with linear integer combinations of the vectors in Equation (
6), two hexagonal lattices rotated by
can be generated.
With the vectors in Equation (
6), the structure matrix is built as:
and it turns out that
, that is, the star
is eutactic for every value of
.
Notice that the vectors in Equation (
6) depend on
hence, according to Equation (
4), the projector
also depends on
. As the column space of
is
, then varying
is the same as projecting onto a rotating
space and this, in turn, can be interpreted as the rotation of the two lattices
and
in between by
. The cut and projection method under this situation can be interpreted as discussed in what follows.
As described in Ref. [
17], for the case of square lattices, two rotated lattices can be viewed as a twin grain boundary and the GCSN (generalized coincidences sites networks) or good fit lattice can be obtained by means of the cut and projection method, using a tube instead of the standard strip. In the case of two hexagonal lattices, we should project from a six-dimensional hypercubic lattice equipped with the standard basis
. The projection matrix
has the property in Equation (
3) and its entries are given by Equation (
4), where the vectors of the star are Equation (
6).
The main idea behind the approach in Ref. [
17] is that, given two lattices rotated in between, a coincidence site lattice can be generated or not, depending on the rotation angle. Then, a generalization of the coincide lattice is proposed that consists in a geometrical array of points of good fit, defined through the concept of
neighbors. We say that two lattice points
and
are neighbors if and only if
where
is the Voronoi polygon of the point
(
). The generalized model postulates that the best fit structure is a set
G of points given by [
17]:
That is, if two points and are neighbors, a point of good fit is located at the middle between and .
It should be stressed that, in our case, the hexagonal lattice and the honeycomb are dual to each other, that is the Voronoi tesselation of the honeycomb is the hexagonal lattice and vice versa. In
Figure 3, two honeycomb structures rotated
in between are shown with black and green points, respectively; the Voronoi tessellation associated to each lattice is also displayed in green and black lines. According to the procedure described above, we should look for neighbouring points, that is, where two points (one black and the other green) are both contained in the intersection of their respective Voronoi polyhedra. Neighbors points are then replaced by a point of good fit located at the middle.
In Ref. [
17], it was shown that, given two lattices
and
rotated in between, the good fit lattice can be obtained by means of the cut and projection method, using a tube instead of the standard strip. Here, we adopt this approach but the way to determine the radius of the tube should be modified since the standard procedure will produce the best fit lattice associated to two rotated hexagonal lattices instead of two rotated honeycombs. As we shown below, the radius of the tube can be determined in order to obtain a honeycomb when
, instead of the hexagonal lattice.
Consider the case of
in Equation (
6), that is, the hexagonal lattice. In this case,
,
and
. The hexagonal lattice can be generated (in a redundant way) using, for instance,
or
:
By considering the vectors
,
and
, the honeycomb is generated as:
where
or 3. That is, the honeycomb is formed by an hexagonal lattice and a shifted copy of itself.
The space
is the row space of
, for
, which turns out to be:
and the distance
d of any lattice point
to
is given by:
From this equation, we can verify that the vectors
,
,
, and
belong to
(thus
), which together with the vectors in Equation (
11) form a central hexagon (large disks in
Figure 4). The distance from the vertices to the incenter of this hexagon is 2, which is
larger than the hexagon generated by the standard basis vectors of
that are project onto Equation (
6) and shown with small disks in
Figure 4. Therefore, linear integer combinations of the vectors
,
,
, and
generate the first hexagonal lattice, say
A, which lie in
, that is, the distance
d of all its points is 0. Now, we now should obtain a second hexagonal lattice
B, displaced by a vector
(dashed arrow in
Figure 4). As in the case of
in Equation (
10), the magnitude of
is obtained as
and it can be verified that
is obtained by the projection of
onto
; actually, two other displacement vectors can be obtained by the projection of either
or
. Any of these vectors are at a distance
from
, thus it seems that accepting points
such that
will produce the required
A and
B lattices. This is not the case, however; since the size of the vectors
in Equation (
6) is
, the projection method will include the basis vectors of
, which project onto an small hexagonal lattice, as shown with small disks in
Figure 4. These points can be excluded by a translation of
along the
direction. Thus, for instance with a translation by
, the basis vectors have distances to
given by
while vectors such as
have distances to
given by
. Consequently, by translating
by
and accepting points such that
will produce the required bipartite lattice.
The utility of this approach is that, once the honeycomb lattices are obtained by projection, the space can be rotated and the projected structure is transformed from periodic (honeycomb) to quasiperiodic with dodecagonal symmetry.
Some projection structures obtained when
is rotated are now shown. Consider a translation of
by
; in this case, we get
and
. In
Figure 5 and
Figure 6, we show projected structures obtained using a tube with radius
(
). The honeycomb for the case
and a dodecagonal quasiperiodic structure when
(
) are shown in
Figure 5, and two structures corresponding to
(the golden angle) and
are shown in
Figure 6 (compare with
Figure 1).
4. Discussion and Conclusions
As shown in
Figure 5 and
Figure 6, many types of lattices can be obtained and described by this approach. These coincidence patterns have paramount importance for the electronic properties. It has been experimentally observed that maximal coincidence regions represent maximal electron density probability regions [
8,
28]. Coincidence regions are where the flat bands, responsible for the strong electron–electron interaction, predominately localize [
8]. As anticipated in the Introduction, strain and disorder play important roles in observing correlated states [
8]. Despite this, our approach allows identifying coincidence points in a general scheme and therefore is very well suited to describe and produce effective renormalized Hamiltonians [
10,
29]. Moreover, the Fourier transform of structures obtained from the cut and projection method are given by the convolution of the window function with the higher-dimensional space reciprocal lattice. Thus, our method suggests a way to describe wavefunctions, Van Hove singularities, and band conductances, as has been made for topological phases in the Quantum Hall effect [
30,
31]. In particular, it is tempting to work on the known relationship between electron wavefunctions in twisted graphene bilayers and the Quantum Hall effect [
10], as, for this effect, it has been proved that topological phases are naturally described by the cut and projection method applied to the magnetic superlattice [
31]. The effects of phason disorder are also open questions [
32], as well as the relationship with flat bands in the Penrose lattice [
33], random binary alloys [
34] and impurity doped graphene [
35].
To produce effective Hamiltonians, care must be taken as hybridization alters the electronic band structure, especially when twisted by a small angle. As the p-orbital of the graphene is not spherical but very anisotropic, the hopping is not a simple function of the distance between the two atoms. However, early models relied in a tight-binding Hamiltonian approach with a hopping parameter
between layers which depends upon the distance
, where
R is the carbon–carbon distance in the horizontal direction of the layers and
z is the distance between layers [
36]. As
z is a constant and
, the hopping is determined by the ratio
. For graphene over graphene, the order of magnitude is such that
. Some models are available in the literature where effective 2D potentials are identified with with some kinds of coincidences defining an effective superlattice [
9].
In conclusion, we provide an approach to describe graphene bilayers using a method originally designed to treat quasicrystals: cut and projection. To achieve this goal, we first observed that twisted heterostructures can be viewed as twin grain boundaries. Then, the cut and projection method is used to generate the coincidence lattice or a best fit lattice. We show that the resulting window function is a tube. Our approach may be useful to describe the electronic properties of twisted graphene over graphene, in which it is known that site coincidences play a fundamental role in the complex quantum phase diagram [
10].