Abstract
We study spectral properties of perturbed discrete Laplacians on two-dimensional Archimedean tilings. The perturbation manifests itself in the introduction of non-trivial edge weights. We focus on the two lattices on which the unperturbed Laplacian exhibits flat bands, namely the \((3.6)^2\) Kagome lattice and the \((3.12^2)\) “Super-Kagome” lattice. We characterize all possible choices for edge weights which lead to flat bands. Furthermore, we discuss spectral consequences such as the emergence of new band gaps. Among our main findings is that flat bands are robust under physically reasonable assumptions on the perturbation, and we completely describe the perturbation-spectrum phase diagram. The two flat bands in the Super-Kagome lattice are shown to even exhibit an “all-or-nothing” phenomenon in the sense that there is no perturbation, which can destroy only one flat band while preserving the other.
Similar content being viewed by others
Avoid common mistakes on your manuscript.
1 Introduction
This paper is about discrete Schrödinger operators on Archimedean tilings, a class of periodic two-dimensional lattices that were already investigated by Johannes Kepler in 1619 [21]. They are natural candidates for the geometry of two-dimensional nanomaterials, and due to advances in this field, most prominently represented by graphene, they have increasingly become a focus of attention.
Much work has been devoted to understanding physical properties of such (new) materials [7, 42, 44]. Most importantly, it can be expected that the underlying geometry, that is the particular lattice, is a key feature determining physical properties of the system. In fact, in particular in the mathematical physics literature, investigations of the connection between the geometry (or topology) of a system and the spectral properties of the associated Hamiltonian have become ubiquitous. Classical examples in this context are so-called quantum waveguides [8,9,10] as well as quantum graphs [1, 3]; see also [26] for a relatively recent reference relevant in our context.
A closely related research direction is superconductivity: the existence of a boundary leads to boundary states in a superconductor with a higher critical temperature than the one of the bulk [17, 38, 40]. In this spirit, it seems very promising to also study the interplay of geometry and many-particle phenomena on Archimedean tilings. Yet another related investigation can be found in [19, 42] where another important quantum phenomenon, namely Bose–Einstein condensation, is examined. It turns out that so-called flat bands, that are infinitely degenerate eigenvalues of the Hamiltonian, play an important role in understanding such many-particle effects, and for other physical phenomena [23]. Flat bands have recently become a topic of increasing attention [2, 39]. One of the central motivations for this paper is to study robustness of flat bands under certain natural perturbations.
Two Archimedean tilings, the \((3.6)^2\) Kagome lattice and the \((3.12^2)\) tilingFootnote 1, which we shall dub Super-Kagome lattice for reasons that will become clear over the course of the article, stand out: they are the only Archimedean lattices on which the discrete, unweighted Laplacian has flat bands. In particular, the Kagome lattice is a prominent model in physics that has recently enjoyed increasing interest [4, 6, 32]. From a mathematical point of view, our paper is motivated by [36] where flat bands for the discrete, unweighted Laplacian on Archimedean tilings have been studied in great detail, in combination with an explicit calculation of the integrated density of states.
A priory, the flat-band phenomena on the Kagome and Super-Kagome lattice seem very sensitive to perturbations: if one replaces the adjacency matrix or the Laplacian by a variant with periodically chosen edge weights, one will generically destroy flat bands. However, the results of this paper suggest that, if one looks at proper, meaningful variants of the discrete Laplacian which respect certain, natural symmetries of the tiling (we call them monomeric Laplacians in Definition 2), then flat bands will persist. Since monomericity is a physically justifiable assumption, this makes a strong case that flat bands are a robust phenomenon, caused by the geometry of the lattice alone and specific to these two lattices, see Theorems 5, and 9.
Other questions of interest on periodic graphs concern existence, persistence and estimates on the width of spectral bands and the gaps between them [28, 28, 33]. We will completely identify the spectra as a function of the perturbation in these cases, see Theorems 7, and 10 as well as Figs. 3, and 5. This provides an exhaustive description of all nanomaterials based on Archimedean tilings on which discrete Laplacians can exhibit flat bands.
Our paper is organized as follows: Sections 2 and 3 are of introductory nature, introducing the notion of and arguing for the relevance of Archimedean tilings, and defining a proper notion of a discrete Laplace operator with non-uniform edge weights. Section 3 also introduces the notion of flat bands and argues why it suffices to restrict our attention to the \((3.6)^2\) Kagome and the \((3.12^2)\) Super-Kagome lattice. Sections 4 and 5 contain our main results on the Kagome and Super-Kagome lattice, respectively. The contributions of this paper are:
-
(i)
We identify the Kagome and Super-Kagome lattice as the only Archimed-ean lattices on which a natural class of periodic, weighted Laplacians can have flat bands (Theorem 4).
-
(ii)
We describe all periodic edge weights, which lead to the maximal possible number of bands on the Kagome and Super-Kagome lattice, and prove that this is equivalent to so-called monomericity of the edge weights (Theorems 5 and 9).
-
(iii)
We completely describe the spectrum in the monomeric Kagome and Super-Kagome lattice (Theorems 7 and 10). In particular, the monomeric Super-Kagome lattice has a surprisingly rich spectrum-perturbation phase diagram (Fig. 5), which might bear relevance for various applications.
-
(iv)
In the Super-Kagome lattice, under a weaker condition than monomericity, namely constant vertex weight, we explicitly describe all remaining “spurious” edge weights which have only one flat band. We describe the topology of this set within the parameter space and show in particular that it is disconnected from the monomeric two-band set (Theorem 11).
2 Archimedean Tilings
Archimedean, Keplerian or regular tilings are edge-to-edge tesselations of the Euclidean plane by regular convex polygons such that every vertex is surrounded by the same pattern of adjacent polygons. We will adopt the notation of [14] and use the (counterclockwise) order of polygons arranged around a vertex as a symbol for a tiling (this is unique up to cyclic permutations), see Fig. 1 for the \((3.6)^2\) Kagome lattice and the \((3.12^2)\) Super-Kagome lattice, which will be investigated in this paper.
The first systematic investigation from 1619 is due to Kepler who identified all 11 such tilings [21]Footnote 2. Most importantly, Archimedean tilings provide natural candidates for geometries of two-dimensional nanomaterials since they form natural, symmetric arrangements of a single building block, positioned at every vertex. And indeed, these lattices can be observed in many naturally occurring materials [11, 12, 24].
From a physical point of view, two-dimensional materials such as graphene are interesting since they feature so-called Dirac points, which are related to a specific behaviour of the electronic band structure of the material [13, 15, 31].
Also note that there are deep connections between Laplacians on these lattices, percolation, and self-avoiding walks, which have also been studied extensively [18, 20, 22, 34, 35, 41, 43, 46]. An important quantity in this context is the so-called connective constant, which is known only in few cases, for example on the hexagonal lattice [5].
3 Defining a Suitable Hamiltonian
Every Archimedean tiling can be regarded as an infinite discrete graph \(G=(V,E)\) with (countable) vertex set V and (countable) edge set E. We write \(v \sim w\) if the vertices v and w are joined by an edge and denote by
the vertex degree of v (which in the case of Archimedean lattice graphs is v-independent). Archimedean lattices are \(\mathbb {Z}^2\)-periodic, and there exists a cofinite \(\mathbb {Z}^2\)-action
that is a group of graph isomorphisms isomorphic to the group \(\mathbb {Z}^2\). It can be intuitively understood as a group of shifts, generated by two linearly independent vectors \(\omega _1\in \mathbb {R}^2\) and \(\omega _2 \in \mathbb {R}^2\), see Figs. 2 and 4 for illustrations. Let \(Q \subset V\) be a minimal (in particular finite) fundamental domain of this action, i.e. the quotient of V under the equivalence relation generated by the group of isomorphisms \((T_\beta )_{\beta \in \mathbb {Z}^2}\).
In the unweighted case, a natural, normalized choice for the Hamiltonian is the discrete Laplacian
as used for instance in [36]. It can be written as \(\Delta f = {\text {Id}}- \frac{1}{|v |} \Pi \) where \(\Pi \) is the adjacency matrix, that is \(\Pi (v,w) = 1\) if \(v \sim w\) and 0 else.
Introducing non-trivial edge weights, a natural candidate for a (normalized) Laplacian—similar to formula (2.11) in [27]—is:
where the edge weights \(\gamma _{vw}=\gamma _{wv} > 0\) and vertex weights \(\mu (v)\) satisfy the relation
As for (1), the spectrum of (2) is contained in [0, 2].
Remark 1
In the literature, one often finds the definition
as a normalized, discrete Laplacian. Note that, whenever \(\mu (v) \ne \mu (w)\) for some \(v \sim w\), then this will not lead to a self-adjoint operator, but it can be made self-adjoint on a suitably weighted \(\ell ^2(V)\)-space, cf. [25]. If all \(\mu (v)\) are the same, then this definition coincides with (2) and can be simplified to
Now, one can prescribe various degrees of the symmetry of the underlying Archimedean lattice to be respected by the Laplacian:
Definition 2
Consider an Archimedean tiling (V, E) with periodic edge weights \(\gamma _{vw} = \gamma _{wv} > 0\), that is \(\gamma _{vw} = \gamma _{T_\beta v T_{\beta } w}\) for all \(v,w \in V\) and \(\beta \in \mathbb {Z}^2\), and corresponding vertex weights \(\mu (v) = \sum _{w \sim v} \gamma _{vw}\). Define the Laplacian \(\Delta _\gamma \) as in (2). Then, we say that the Archimedean tiling with Laplacian \(\Delta _\gamma \)
-
(1)
has constant vertex weight, if there is \(\mu > 0\) such that \(\mu (v) = \mu \) for all \(v \in V\).
-
(2)
is monomeric if for all vertices \(v \in V\) the list of edge weights, arranged cyclically around v, coincides (up to cyclic permutations).
Clearly, (2) is stronger than (1). However, in either case, the Laplacian reduces to (4).
The term “monomeric” is inspired by the fact that the associated operators can be interpreted as describing properties of nanomaterials formed from one type of monomeric building block, positioned at every vertex of an Archimedean tiling. Clearly, monomeric Laplacians on Archimedean lattices have constant vertex weights, but the converse is not true in general. However, we will see in Theorems 5 and 9 that on the Kagome and Super-Kagome lattice, the validity of the converse implication is equivalent to existence (or persistence) of all flat bands. Also, monomericity seems a physically reasonable assumption for nanomaterials, which suggests that the emergence of flat bands, while a priori very sensitive to perturbations of coefficients in the operator, might nevertheless be robust within the class of physically relevant operators. Let us note that in the literature, one also finds investigations of other perturbations to periodic operators such as by potentials [4, 39] or magnetic fields [4].
Next, let \(\mathbb {T}^2 = \mathbb {R}^2 / \mathbb {Z}^2\) be the flat torus and define for every \(\theta \in \mathbb {T}^2\) the |Q|-dimensional Hilbert space
with inner product
Note that this inner product is independent of the choice of the fundamental cell. Given the Laplacian (4) on \(\ell ^2(V)\) with properties described in Definition 2, we define on \(\ell ^2(V)_{\theta }\) the operator
Clearly, (5) can be represented as a \(|Q |\)-dimensional Hermitian matrix. Due to Floquet theory, we have
and the following statement holds.
Proposition 3
Let \(E \in \mathbb {R}\). Then, the following are equivalent:
-
(i)
\(E \in \sigma ( \Delta _\gamma ^\theta )\) for all \(\theta \in \mathbb {T}^2\).
-
(ii)
\(E \in \sigma ( \Delta _\gamma ^\theta )\) for a positive measure subset of \(\theta \in \mathbb {T}^2\).
-
(iii)
There is an infinite orthonormal family eigenfunctions of \(\Delta _\gamma \) to the eigenvalue E. Each of them can be chosen to be supported on a finite number of vertices.
If any of (i) to (iii) is satisfied, we say that \(\Delta _\gamma \) has a flat band (at energy E).
The proof of Proposition 3 can be found in [29, 30]. Note that, in the \(\ell ^\infty (V)\) setting instead of the \(\ell ^2(V)\) setting, such infinitely degenerate eigenvalues are also referred to as “black hole eigenvalues” in [45]. Also, the existence of flat bands can be interpreted as a breakdown of the unique continuation principle [37].
In the Hilbert space setting working on \(\ell ^2(V)\), it is known that for constant edge weights, the discrete Laplacian has flat bands only on two of the 11 Archimedean lattices, namely the \((3.6)^2\) Kagome lattice and the \((3.12^2)\) Super-Kagome lattice [36]. Before turning to perturbed versions of those two lattices, one should verify that there won’t be any surprises on the other lattices:
Theorem 4
On the Archimedean lattices \((4^4)\), \((3^6)\), \((6^3)\), \((3^3.4^2)\), \((4.8^2)\), \((3^2.4.3.4)\), (3.4.6.4), (4.6.12), \((3^4.6)\), there is no choice of periodic (with respect to the fundamental cell on the lattice) edge weights \(\gamma _{vw} = \gamma _{wv} > 0\) which will make the weighted adjacency matrix
have a flat band. Consequently, also the Laplacian with constant or monomeric edge weights has no flat bands on these lattices.
Theorem 4 is proved in Appendix A by a series of straightforward but somewhat lengthy calculations in which one identifies certain terms in the associated characteristic polynomials and shows that there are no \(\theta \)-independent roots, employing Proposition 3 (this should be compared to the proofs of Theorems 5 and 9). Theorem 4 justifies to restrict our attention to the (perturbed) Kagome and Super-Kagome lattices from now on.
4 The Perturbed Kagome Lattice
In this section, we discuss the Kagome lattice with non-uniform (periodic) edge weights. The elementary cell of the Kagome lattice contains three vertices and six edges (one can think of the edges as arranged around a hexagon).
A priori, periodicity allows for six edge weights \(\gamma _1,...,\gamma _6 > 0\), and the Floquet Laplacian \(\Delta _\gamma ^\theta \) can be written as the Hermitian matrix
where \(w:=e^{i\theta _1}\) and \(z:=e^{i\theta _2}\). We denote the three real eigenvalues of \(\Delta ^{\theta }_{\gamma }\) by \(\lambda _1(\theta ,\gamma ) \le \lambda _2(\theta ,\gamma ) \le \lambda _3(\theta ,\gamma )\).
Note that the six degrees of freedom are to be further reduced, depending on the following symmetry conditions:
-
If we merely assume a constant vertex weight \(\mu > 0\), then identity (3) will impose the three additional linearly independent conditions
$$\begin{aligned} \begin{aligned} \gamma _1+\gamma _4&=\gamma _2+\gamma _5\ , \\ \gamma _3+\gamma _6&=\gamma _2+\gamma _5\ ,\\ \gamma _1+\gamma _3+\gamma _4+\gamma _6&= \mu \ , \end{aligned} \end{aligned}$$(7)and we end up with three degrees of freedom.
-
If we also assume monomericity, then it is easy to see that the only choice is the breathing Kagome lattice, cf. [16], with an edge weight \(\alpha > 0\) on all edges belonging to upwards pointing triangles and edge weight \(\beta > 0\) on all edges belonging to downwards pointing triangles, where \(2 (\alpha + \beta ) = \mu \). After fixing the vertex weight \(\mu \), this amounts to only one degree of freedom.
4.1 Flat Bands in the Perturbed Kagome Lattice
Theorem 5
Consider the perturbed Kagome lattice with Laplacian (4), fixed vertex weight \(\mu > 0\) and periodic edge weights \(\gamma _1,...,\gamma _6 > 0\), satisfying the condition (3) on vertex and edge weights. Then, the following are equivalent:
-
(i)
There exists a flat band.
-
(ii)
The vertex weights are monomeric. More explicitly, there are \(\alpha , \beta > 0\) with \(2 (\alpha + \beta ) = \mu \) such that
$$\begin{aligned}\begin{aligned} \gamma _{2}&=\gamma _{4}=\gamma _{6}:=\alpha , \\ \gamma _{1}&=\gamma _{3}=\gamma _{5}:=\beta . \end{aligned} \end{aligned}$$
The rest of this subsection is devoted to the proof of Theorem 5. We start by identifying flat bands using the weighted adjacency matrix
which is spectrally equivalent to \(\Delta _\gamma ^\theta \) up to scaling and shifting via the relation
In order to find flat bands, we will identify conditions for \(\theta \)-independent eigenvalues of \(\Pi _\gamma ^\theta \) and therefore calculate
where \(A:=\gamma _3+w\gamma _6\), \(B:=w\gamma _4+z\gamma _1\) and \(C:=\gamma _2+z\gamma _5\). Rearranging the terms yields
The prefactors
are linearly independent as measurable functions of \(\theta \) on \(\mathbb {T}^2\). Consequently, since all \(\gamma _i\) are positive, \(\theta \)-independent eigenvalues exist if and only if the w and z-independent terms in every line are zero. This is only possible for negative \(\lambda \), which (possibly after scaling the \(\gamma _i\) and \(\mu \) for the moment) can be assumed to equal \(-1\). Therefore, we obtain the conditions
and
Lemma 6
The only positive solutions (meaning all \(\gamma _i\) are nonzero) of (7), (9), (10) are
with \(x,y \in (0,1)\) and \(x+y=1\).
Proof
By a direct calculation (11) solves (7), (9), (10).
Conversely, assume that there are positive solutions \(\gamma _1,...,\gamma _6 > 0\). From (9) we obtain
and this implies \(\gamma _6 > \gamma _2 \gamma _4\) and \(\gamma _4 > \gamma _2\gamma _6\). Hence, combining both equations yields \(\gamma ^2_2\gamma _6 < \gamma _6\) which shows that \(\gamma _2 < 1\). In the same way, one proves \(\gamma _i < 1\) for every other i.
Next, let \(\gamma _2+\gamma _5:=\Lambda \). By (7) one immediately concludes \(\gamma _1+\gamma _4=\gamma _3+\gamma _6=\Lambda \). Now, we add (9) and (10) and rearrange the equations to obtain
By repeated factorization, the right-hand side simplifies to
and since for the left-hand side one has
we arrive at the polynomial \(\Lambda ^3-\frac{3\Lambda ^2}{2}+\frac{1}{2}=0\) the only positive solution of which is \(\Lambda =1\). Finally, adding the first the two equations of (9) yields
and this implies \(\gamma _5=\gamma _3\). Furthermore, adding the last two equations gives
giving \(\gamma _4=\gamma _2\). Conditions (7) hence give \(\gamma _1=\gamma _5\) and \(\gamma _6=\gamma _2\). This proves the statement. \(\square \)
We are now in the position to prove Theorem 5.
Proof of Theorem 5
Comparing \(\Pi ^{\theta }_\gamma \) with \(\Delta ^{\theta }_{\gamma }\) we conclude that \(\Delta ^{\theta }_\gamma \) has a flat band with edge weights \(\gamma _1,...,\gamma _6\) if and only if there exists \(\delta > 0\) such that \(\Pi ^{\theta }_{\gamma }\) has a flat band for edge weights \(\delta \gamma _1,...,\delta \gamma _6\). From this observation the statement follows directly taking Lemma 6 into account. \(\square \)
4.2 The Spectrum and Band Gaps in the Monomeric Kagome Lattice
In the case where the perturbed Kagome lattice has a flat band, we further study the structure of the rest of the spectrum. We reiterate that, due to Theorem 5, the existence of a flat band is equivalent to the weights being monomeric.
As shown, for instance, in [36], in the case where all edge weights are equal, the two other spectral bands, generated by the two other \(\theta \)-dependent eigenvalues of \(\Delta ^{\theta }_{\gamma }\), touch at \(E = 3/4\), and the derivative of the integrated density of states at \(E = 3/4\) vanishes – an indication that the spectral density at 3/4 is sufficiently thin for a gap to form under perturbation. And indeed, this is the statement of the next theorem, which also characterizes the width of the gap.
Theorem 7
(Band gaps in the perturbed Kagome lattice). Consider the perturbed Kagome lattice with fixed vertex weight \(\mu > 0\), and monomeric edge weights \(\alpha , \beta > 0\), satisfying \(2 (\alpha + \beta ) = \mu \) as characterized in Theorem 5. Then, the spectrum is given by
Furthermore, there is always a flat band at \(\frac{3}{2}\).
Remark 8
Theorem 7 states that, as soon as \(\alpha \ne \beta \), or alternatively, \(\alpha \ne \frac{\mu }{4}\), a spectral gap of width
will form around \(\frac{3}{4}\), see also Fig. 3. The flat band at \(\frac{3}{2}\) will always be connected to the energy band below it which means that the “touching” of the flat band at \(\frac{3}{2}\) is protected in the class of monomeric perturbations.
Proof
A calculation shows that the eigenvalues of \(\Delta _\gamma ^\theta \) with the choice \(2( \alpha + \beta )=\mu \) as in Theorem 5 are given by
and
where \(F(\theta ):=\cos (\theta _1) + \cos (\theta _2) + \cos (\theta _1 - \theta _2)\). The function \(\mathbb {T}^2 \ni \theta \mapsto F(\theta )\) takes all values in \([-3/2, 3]\), see Lemma 3.1 in [36], whence \(\lambda _1(\theta , \gamma )\) and \(\lambda _2(\theta , \gamma )\) take all values in the intervals
\(\square \)
5 The Perturbed Super-Kagome Lattice
In this section, we investigate the Archimedean tiling \((3.12^2)\) which we call Super-Kagome lattice. Its minimal elementary cell contains six vertices and nine edges: three edges on upwards pointing triangles, three edges on downwards pointing triangles, and three edges bordering two dodecagons, see Fig. 4.
Given a constant vertex weight \(\mu > 0\), the Floquet Laplacian (5) is a \(6\times 6\)-matrix given by
where \(w:=e^{i\theta _1}\), \(z:=e^{i\theta _2}\).
-
If we fix a constant vertex weight \(\mu > 0\), the condition \(\sum _{w \sim v} \gamma _{vw} = \mu \) for all \(v \in V\) leads to
$$\begin{aligned} \begin{aligned} \mu = \gamma _2+\gamma _3+\gamma _7 = \gamma _5+\gamma _6+\gamma _7&= \gamma _1+\gamma _2+\gamma _8 = \gamma _4+\gamma _5+\gamma _8 \\&= \gamma _1+\gamma _3+\gamma _9 = \gamma _4+\gamma _6+\gamma _9. \end{aligned} \end{aligned}$$(14)This can be seen to be a linear system of 6 linearly independent equations with 9 unknowns, so the solution space is 3-dimensional. More precisely, by appropriate additions, we infer the three identities
$$\begin{aligned} \begin{aligned} 2 \gamma _1 + \gamma _8 + \gamma _9&= 2 \gamma _7 + \gamma _2 + \gamma _3, \\ 2 \gamma _4 + \gamma _8 + \gamma _9&= 2 \gamma _7 + \gamma _5 + \gamma _6, \\ \gamma _2 + \gamma _3&= \gamma _5 + \gamma _6 \end{aligned} \end{aligned}$$(15)which imply \(\gamma _1 = \gamma _4\). The identities \(\gamma _2 = \gamma _5\), and \(\gamma _3 = \gamma _6\) follow by completely analogous calculations. This leaves us with 6 independent variables \(\gamma _1, \gamma _2, \gamma _3\), and \(\gamma _7, \gamma _8, \gamma _9\) which are, however, still subject to the three conditions
$$\begin{aligned} \gamma _2 + \gamma _3 + \gamma _7 = \gamma _1 + \gamma _2 + \gamma _8 = \gamma _1 + \gamma _3 + \gamma _9 = \mu \end{aligned}$$from (14). Therefore, we are left with three degrees of freedom.
-
If we additionally prescribe monomericity, it is easy to see that there is only one degree of freedom: All edges around triangles carry the weight \(\alpha > 0\), and all remaining edges (separating two dodecagons) carry the weight \(\beta > 0\) under the condition \(2 \alpha + \beta = \mu \).
5.1 Flat Bands in the Perturbed Super-Kagome Lattice
Theorem 9
Consider the perturbed Super-Kagome lattice with Laplacian (4), fixed vertex weight \(\mu > 0\), and periodic edge weights \(\gamma _1, \dots , \gamma _9 > 0\) satisfying the condition (3) on vertex and edge weights. Then, the following are equivalent:
-
(i)
There exist exactly two flat bands.
-
(ii)
The Super-Kagome lattice is monomeric. More explicitly, there are \(\alpha ,\beta > 0\) such that \( 2\alpha + \beta = \mu \) together with
$$\begin{aligned} \begin{aligned} \gamma _1=\gamma _2=\gamma _3=\gamma _4=\gamma _5=\gamma _6&=\alpha \ , \\ \gamma _7=\gamma _8=\gamma _9&=\beta \ . \end{aligned} \end{aligned}$$
Proof
Recall that in the constant vertex weight case, we have
and consider the weighted adjacency matrix
which is a shifted and scaled version of \(\Delta ^{\theta }_{\gamma }\). We calculate
Since \(w + {\overline{w}} = 2 \cos (\theta _1)\), \(z + {\overline{z}} = 2 \cos (\theta _2)\), and \(w {\overline{z}} + {\overline{w}} z = 2 \cos (\theta _1 - \theta _2)\) are linearly independent on \(\mathbb {T}^2\), \(\lambda \) is a \(\theta \)-independent eigenvalue if and only if the conditions
as well as
hold.Footnote 3 Conditions (17) imply that any \(\theta \)-independent eigenvalue of the matrix \(\Pi ^{\theta }_{\gamma }\) must satisfy
Since all \(\gamma _i\) are positive, the only way for these three equations to have the same set of solutions, that is for two flat bands to exist, is therefore
together with
This implies that the matrix \(\Pi _\gamma ^\theta \) can only have two \(\theta \)-independent eigenvalues if there are \(\alpha , \beta > 0\) with
that is the monomeric case, and the only candidates for these eigenvalues are \(-\beta \pm \alpha \). To see that they are indeed eigenvalues, one verifies by an explicit calculation that condition (18) is also fulfilled. This shows the stated equivalence.\(\square \)
Next, we further describe the spectrum of the monomeric Super-Kagome lattice.
Theorem 10
(Band gaps in the perturbed Super-Kagome lattice). Consider the perturbed Super-Kagome lattice with Laplacian (4) with fixed vertex weight \(\mu > 0\) and monomeric edge weights \(\alpha , \beta > 0\), satisfying \(2 \alpha + \beta = \mu \) as characterized in Theorem 9. Then, the spectrum is given by:
with flat bands at \(\frac{3 \alpha }{\mu }\) and \(2 - \frac{\alpha }{\mu }\).
The spectrum and the position of the flat bands are plotted in Fig. 5. The spectrum generically consists of two distinct intervals (bands) except for the case \(3 \alpha = 2 \beta \), that is \(\alpha = \frac{2 \mu }{7}\), in which the two bands touch and the spectrum consists of one interval with an embedded flat band in the middle as well as a flat band at its maximum. This case \(\alpha = \frac{2 \mu }{7}\) connects two regimes with different spectral pictures:
-
If \(\alpha > \frac{2 \mu }{7}\), the spectrum consists of two intervals the upper one of which has two flat bands at its endpoints. In the special case of uniform edge weights (that is \(\alpha = \frac{\mu }{3})\), this has already been observed, for instance in [36].
-
If \(\alpha < \frac{2 \mu }{7}\), the spectrum will again consist of two intervals each of which will have a flat band at its maximum. Somewhat surprisingly, the lower flat band has now attached itself to the lower interval \(I_2\) upon passing the critical parameter \(\alpha = \frac{2 \mu }{7}\).
Another noteworthy observation is that no gap opens within the intervals \(I_1\) and \(I_2\), despite them being generated by two distinct Floquet eigenvalues and the density of states measure vanishing at a point in the interior of the bands, see again [36] for plots of the integrated density of states in the case of constant edge weights. In particular, this distinguishes the monomeric Super-Kagome lattice from the monomeric Kagome lattice where such a gap indeed opens within the spectrum at points of zero spectral density.
Proof of Theorem 10
In the monomeric case, the characteristic polynomial \( \det (\lambda {\text {Id}}- \Pi _\gamma ^\theta )\) of the matrix \(\Pi _\gamma ^\theta \) simplifies to
where \(F(\theta _1, \theta _2) = \cos (\theta _1) + \cos (\theta _2) + \cos (\theta _1 + \theta _2)\). Its six roots are
whence the eigenvalues of \(\Delta _\gamma ^\theta \) are given by
Using that the map \(\mathbb {T}^2 \ni (\theta _1, \theta _2) \mapsto F(\theta _1, \theta _2)\) takes all values in the interval \((- 3/2, 3)\), we conclude that the bands, generated by \(\lambda _1(\theta , \gamma )\) and \(\lambda _2(\theta , \gamma )\), as well as the bands generated by \(\lambda _4(\theta , \gamma )\) and \(\lambda _5(\theta , \gamma )\) always touch, and the spectrum consists of the two intervals
\(\square \)
One might now wonder under which conditions only one flat band exists. The next theorem completely identifies all parameters for which one flat band exists:
Theorem 11
Consider the perturbed Super-Kagome lattice with Laplacian (4), fixed vertex weight \(\mu > 0\), and periodic edge weights \(\gamma _1, \dots , \gamma _9 > 0\) satisfying the condition (3) on vertex and edge weights. The set of \((\gamma _i)\) such that exactly one flat band exists consists of six connected components which have no mutual intersections and have no intersection with the two-flat-band parameter set, identified in Theorem 9.
The solution space is invariant under those permutations of the \(\gamma _i\), which correspond to rotations of the lattice by \(\frac{2\pi }{3}\), and \(\frac{4\pi }{3}\). Modulo these permutations, the two connected components can be described as follows
-
A one-dimensional submanifold, isomorphic to an interval, and explicitely descibed in equation (26),
-
Two one-dimensional submanifolds each isomorphic to an interval, explicitely described in (28), and (30), which intersect in a single point.
Proof of Theorem 11
Recall that due to the reductions made at the beginning of the section, after fixing the constant vertex weight \(\mu > 0\), the space of edge weights is a 3-dimensional manifold in the 6-dimensional parameter space \(\{ \gamma _1, \gamma _2, \gamma _3, \gamma _7, \gamma _8, \gamma _9 > 0 \}\), subject to the conditions
Furthermore, from the proof of Theorem 9 we infer that \(\Delta _\gamma \) has a flat band at \(\lambda \) if and only if the weighted adjacency matrix \(\Pi _\gamma ^\theta \) has the \(\theta \)-independent eigenvalue \(\tilde{\lambda }:= \mu (1-\lambda )\). This requires in particular that
holds with a certain combination of plus and minus signs. Now, if equality in (22) holds with all three signs positive or all three signs negative, respectively, then the argument in the proof of Theorem 9 shows that this already implies that the edge weights are monomeric, the identities also hold with the opposite sign, the additional condition (18) is fulfilled, and there are two flat bands. As a consequence, the only chance for the existence of exactly one flat band is (22) to hold with different signs in front of \( \gamma _7, \gamma _8, \gamma _9 \). Also, it is immediately clear that (22) with different signs does not allow for a monomeric and nonzero solution, and hence, the solution space consists of at most six mutually disjoint components which have no intersection with the two-flat-band manifold, identified in Theorem 9.
By symmetry, it suffices to investigate two out of these six cases:
and
To solve Case(- + +), combine the second identities in (21) and (23), to deduce
which, recalling \(\gamma _i > 0\), is only possible if \(\gamma _1 = \gamma _3\). But then, by (23), \(\gamma _7 = \gamma _8\). Calling \(\alpha ' := \gamma _2\), and \(\beta ' := \gamma _9\), we can use (21), to further express
Next, we eliminate \(\beta '\) by resolving the yet unused first identity in (23), which yields
This only has real solutions if \(\alpha '> \frac{8}{17} \mu > \frac{1}{3} \mu \), thus only
can be a positive solution. Furthermore, we need \(\beta ' \in (0, \mu )\), which is the case if and only if
We therefore find the one-parameter solution set
with energy
Finally, an explicit calculation shows that with these parameters, (18) is indeed fulfilled.
As for Case(+ - -), we combine the second identity in (21) with the second identity in (24) to deduce
Identity (27) has two types of solutions:
Case(+ - -)(a): \(\gamma _1 = \gamma _3\).
As before we find \(\gamma _7 = \gamma _8\). Let \(\alpha ' := \gamma _2\), \(\beta ' := \gamma _9\), and combine the remaining first identity in (24) with (25) to solve for \(\beta '\), finding
Only the solution
has a chance to be in \((0, \mu )\), and, indeed, this is the case if and only if
We obtain the one-parameter solution set
with energy
Again, an explicit calculation shows that (18) is fullfilled.
Case(+ - -)(b): The other solution of (27) is
We set \(\alpha '' := \gamma _1\), \(\beta '' := \gamma _3\), whence
and use (21) to infer
Plugging (29) into the yet unused first identity in (24), we arrive at
We observe that only the solution with a plus has a chance to be positive and it is easy to see that this solution takes values in \((0, \mu )\) for all \(\alpha '' \in (0, \mu )\). We obtain the one-parameter solution set
at energy
Again, an explicit calculation verifies that with these choices, (18) is fullfilled.
Finally, to conclude the claimed topological properties of the manifolds, we need to verify that the solution space (28) in Case(+ - -)(a) intersects the solution space (30) in Case(+ - -)(b) if and only if
\(\square \)
Remark 12
Theorems 9 and 11 imply that the six one-flat-band components and the two-flat-band component are mutually disjoint. However, a closer analysis of the extremal cases in Formulas (26), (28), and (30), as well as of the monomeric case, implies that when sending the parameters to their extremal values, the three one-dimensional manifolds corresponding to Case(+ - -) (a), and the two-flat-band-manifold of solutions converge to the two points
which themselves do no longer belong to the space of admissible parameters. Likewise, the limit of solutions of Case(+ - -) in (26) corresponding to \(\alpha ' = \frac{\mu }{2}\) corresponds to the point \(X_2\), see also Fig. 6.
Notes
We explain the notation for the lattices in Sect. 2.
All 11 Archimedean tilings are: the \((4^4)\) rectangular tiling, the \((3^6)\) triangular tiling, the \((6^3)\) hexagonal tiling, the \((3.6)^2\) Kagome lattice, the \((3.12^2)\) Super-Kagome lattice, the \((3^3.4^2)\) tiling, the \((4.8^2)\) tiling, the \((3^2.4.3.4)\) tiling, the (3.4.6.4) tiling, the (4.6.12) tiling, and the \((3^4.6)\) tiling.
As we will see later, despite its complexity, (18) will not impose further restrictions and holds in all relevant cases. This appears to be a consequence of symmetries of the lattice and the operator.
References
Baradaran, M., Exner, P.: Kagome network with vertex coupling of a preferred orientation. J. Math. Phys. 63(8), 083502 (2022)
Becker, S., Embree, M., Wittsten, J., Zworski, M.: Mathematics of magic angles in a model of twisted bilayer graphene. Prob. Math. Phys. 3, 69–103 (2022)
Berkolaiko, G., Kuchment, P.: Introduction to Quantum Graphs. American Mathematical Society, Providence (2013)
Bilitewski, T., Moessner, R.: Disordered flat bands on the kagome lattice. Phys. Rev. B 98, 235109 (2018)
Duminil-Copin, H., Smirnov, S.: The connective constant of the honeycomb lattice equals \(\sqrt{2+\sqrt{2}}\). Ann. Math. 175, 1653–1665 (2012)
Dias, D.P.: Topological properties of flat bands in generalized Kagome lattice materials. Thesis: https://kth.diva-portal.org/smash/get/diva2:1591789/FULLTEXT01.pdf
de Lima, F.C., Ferreira, G.J., Miwa, R.H.: Topological flat band, Dirac fermions and quantum spin Hall phase in 2d Archimedean lattices. Phys. Chem. Chem. Phys. 21(40), 22344–22350 (2019)
Exner, P., Kovařík, H.: Quantum Waveguides. Springer, New York (2015)
Exner, P.: Spectral properties of soft quantum waveguides. J. Phys. A Math. Theor. 53(35), 355302 (2020)
Exner, P.: Soft quantum waveguides in three dimensions. J. Math. Phys. 63(4), 042103 (2022)
Frank, F.C., Kasper, J.S.: Complex alloy structures regarded as sphere packings. I. Definitions and basic principles. Acta. Cryst 11, 184–190 (1958)
Frank, F.C., Kasper, J.S.: Complex alloy structures regarded as sphere packings. II. Analysis and classification of representative structures. Acta. Cryst 12, 483–499 (1959)
Fefferman, C.L., Weinstein, M.I.: Honeycomb lattice potentials and Dirac points. J. Am. Math. Soc. 25(4), 1169–1220 (2012)
Grünbaum, B., Shephard, G.C.: Tilings and Patterns. An Introduction. A Series of Books in the Mathematical Sciences. W. H. Freeman and Company, New York (1989)
Hou, J.-M., Chen, W.: Hidden symmetry and protection of Dirac points on the honeycomb lattice. Sci. Rep. 5, 17571 (2015)
Herrera, M.A.J., Kempkes, S.N., de Paz, M.B., García-Etxarri, A., Swart, I., Smith, C.M., Bercioux, D.: Corner modes of the breathing Kagome lattice: Origin and robustness. Phys. Rev. B 105, 085411 (2022)
Hainzl, C., Roos, B., Seiringer, R.: Boundary superconductivity in the BCS model. J. Spectr. Theory 12, 1507–1540 (2022)
Jacobsen, J.L.: High-precision percolation thresholds and Potts-model critical manifolds from graph polynomials. J. Phys. A Math. Theor. 47, 135001 (2014)
Julku, A., Bruun, G.M., Törmä, P.: Quantum geometry and flat band Bose–Einstein condensation. Phys. Rev. Lett. 127, 17 (2021)
Jacobsen, J.L., Scullard, C.R., Guttmann, A.J.: On the growth constant for square-lattice self-avoiding walks. J. Phys. A Math. Theor. 49(49), 494004 (2016)
Kepler, J.: Harmonices Mundi. Johann Plack (1619)
Kesten, H.: The critical probability of bond percolation on the square lattice equals \({\frac{1}{2}}\). Commun. Math. Phys. 74(1), 41–59 (1980)
Kollár, A.J., Fitzpatrick, M., Sarnak, P., Houck, A.A.: Line-graph lattices: Euclidean and non-Euclidean flat bands, and implementations in circuit quantum electrodynamics. Commun. Math. Phys. 376(3), 1909–1956 (2019)
Kulkarni, A.A., Hanson, E., Zhang, R., Thornton, K., Braun, P.V.: Archimedean lattices emerge in template-directed eutectic solidification. Nature 577(7790), 355–358 (2020)
Keller, M., Lenz, D., Wojciechowski, R.K.: Graphs and Discrete Dirichlet Spaces. Springer, Switzerland (2021)
Kuchment, P., Post, O.: On the spectra of carbon nano-structures. Commun. Math. Phys. 275(3), 805–826 (2007)
Korotyaev, E., Saburova, N.: Two-sided estimates of total bandwidth for Schrödinger operators on periodic graphs. Commun. Pure Appl. Anal. 21(5), 1691–1714 (2022)
Korotyaev, E., Saburova, N.: Spectral estimates for Schrödinger operators on periodic discrete graphs. St. Petersburg Math. J. 30(4), 667–698 (2019)
Kuchment, P.A.: On the Floquet theory of periodic difference equations. In: Geometrical and Algebraical Aspects in Several Complex Variables (Cetraro,: Sem. Conf.), vol. 8. EditEl, Rende , 1991, pp. 201–209 (1989)
Kuchment, P.: An overview of periodic elliptic operators. Bull. Am. Math. Soc. 53(3), 343–414 (2016)
Liu, Z., Wang, J., Li, J.: Dirac cones in two-dimensional systems: from hexagonal to square lattices. Phys. Chem. Chem. Phys. 15, 18855–18862 (2013)
Mielke, C., Das, D., Yin, J.X., et al.: Time-reversal symmetry-breaking charge order in a Kagome superconductor. Nature 602, 245–250 (2022)
Mohar, B., Woess, W.: A survey on spectra of infinite graphs. Bull. Lond. Math. Soc. 21(3), 209–234 (1989)
Nienhuis, B.: Exact critical point and critical exponents of \(\rm O (n)\) models in two dimensions. Phys. Rev. Lett. 49, 1062–1065 (1982)
Parviainen, R.: Estimation of bond percolation thresholds on the Archimedean lattices. J. Phys. A Math. Theor. 40, 9253–9258 (2007)
Peyerimhoff, N., Täufer, M.: Eigenfunctions and the integrated density of states on Archimedean tilings. J. Spectr. Theory 11, 461–488 (2021)
Peyerimhoff, N., Täufer, M., Veselić, I.: Unique continuation principles and their absence for Schrödinger eigenfunctions on combinatorial and quantum graphs and in continuum space. Nanosyst.: Phys. Chem. Math. 8(2), 216–230 (2017)
Samoilenka, A., Babaev, E.: Boundary states with elevated critical temperatures in Bardeen–Cooper–Schrieffer superconductors. Phys. Rev. B 101, 13 (2020)
Sabri, M., Youssef, P.: Flat bands of periodic graphs. J. Math. Phys. 64, 092101 (2023)
Samoilenka, A., Babaev, E.: Microscopic derivation of superconductor-insulator boundary conditions for Ginzburg-Landau theory revisited: Enhanced superconductivity at boundaries with and without magnetic field. Phys. Rev. B 103, 224516 (2021)
Sykes, M.F., Essam, J.W.: Exact critical percolation probabilities for site and bond problems in two dimensions. J. Math. Phys. 5, 1117–1127 (1964)
Shi, M., Yu, F., Yang, Y., et al.: A new class of bilayer Kagome lattice compounds with Dirac nodal lines and pressure-induced superconductivity. Nat. Commun. 13, 2773 (2022)
Suding, P.N., Ziff, R.M.: Site percolation thresholds for Archimedean lattices. Phys. Rev. E 60(1), 275–283 (1999)
Tarnai, T., Fowler, P.W., Guest, S.D., Kovács, F.: Equiauxetic hinged archimedean tilings. Symmetry 14, 232 (2022)
von Below, J., Lubary, J.A.: Isospectral infinite graphs and networks and infinite eigenvalue multiplicities. Netw. Heterog. Media 4(3), 453–468 (2009)
Veselić, I.: Spectral analysis of percolation Hamiltonians. Math. Ann. 331(4), 841–865 (2004)
Acknowledgements
JK would like to thank the Bergische Universität Wuppertal where parts of this project were done while being on leave from the FernUniversität in Hagen. This article is based on work from COST Action 18232 MAT-DYN-NET, supported by COST (European Cooperation in Science and Technology), https://www.cost.eu/actions/CA18232/. MT would like to thank the Mittag-Leffler Institute where parts of this work were initiated during the trimester Program “Spectral Methods in Mathematical Physics”.
Funding
Open Access funding enabled and organized by Projekt DEAL.
Author information
Authors and Affiliations
Corresponding author
Additional information
Communicated by Alain Joye.
Publisher's Note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Appendix A. No Flat Bands on Other Archimedean Lattices
Appendix A. No Flat Bands on Other Archimedean Lattices
In this appendix, we prove Theorem 4 stating that, beyond the Kagome and the Super-Kagome lattice, on no other Archimedean lattice any set of positive, periodic weights \(\gamma _j\) can make the weighted adjacency matrix \(\Pi _\gamma \) have flat bands.
Throughout the section, we write \(w := \exp (i \theta _1)\), and \(z := \exp (i \theta _2)\). Clearly, \({\overline{w}} = w^{-1}\), and \({\overline{z}} = z^{-1}\). In particular, the characteristic polynomial can be expanded in the form
where n is a finite number, depending on the particular lattice, and the \(P_{\alpha , \beta , \gamma } (\lambda )\) are polynomials in \(\lambda \) with coefficients themselves polynomials in the \(\gamma _i\). Since \(\Pi _\gamma \) is Hermitian,
The monomials \(w^\alpha z^\beta \) are trigonometric functions of \(\theta _1\) and \(\theta _2\). In particular, they are linearly independent as functions of \(\theta \). As a consequence, to see that there is no flat band, it suffices to see that one of the polynomials \(P_{\alpha , \beta , \gamma }\) has no root. Let us summarize this strategy:
Lemma 13
Let \(\Pi _\gamma \) be the weighted \(\theta \)-dependent adjacency matrix. Expand its characteristic polynomial in the form
Then, \(\lambda \) is a flat band if and only if \(P_{\alpha ,\beta ,\gamma }(\lambda ) = 0\) for all \(\alpha , \beta \).
In particular, if for some \(\alpha _0, \beta _0\), \(P_{\alpha _0, \beta _0, \gamma }\) is a nonzero monomial, only depending on the \(\gamma _j\), there cannot be a flat band.
For the \((4^4)\), \((6^3)\), \((3^6)\), and \((3^3.4^2)\) lattice, we will exclude the existence of flat bands by a direct calculations, whereas for the remaining lattices \((4.8^2)\), \((3^2.4.3.4)\), (3.4.6.4), (4.6.12), and \((3^4.6)\), we will resort to Lemma 13.
A.1. \((4^4)\) A fundamental cell is given by one vertex with two edge weights to choose. The weighted adjacency matrix \(\Pi _\gamma \) then has the only entry and hence the only eigenvalue
This is clearly \(\theta \)-dependent, so there cannot be a flat band.
A.2. \((3^6)\) A fundamental cell is given by one vertex with three edge weights to choose. The weighted adjacency matrix \(\Pi _\gamma \) then has the only entry
This is clearly \(\theta \)-dependent, so there cannot be a flat band.
A.3. \((6^3)\) A fundamental cell is given by two vertices with three edge weights to choose. The weighted adjacency \(\Pi _\gamma \) has the form
Its spectrum is
Both eigenvalues are clearly \(\theta \)-dependent, so there cannot be a flat band.
A.4. \((3^3.4^2)\) A fundamental cell is given by two vertices with five edge weights to choose. The weighted adjacency \(\Pi _\gamma \) has the form
Its spectrum is
Both eigenvalues are clearly \(\theta \)-dependent, so there cannot be a flat band.
From now on, we will use Lemma 13 to exclude existence of flat bands.
A.5. \((4.8^2)\) A fundamental cell is given by four vertices with six edge weights to choose. The weighted adjacency matrix \(\Pi _\gamma \) has the form
Arranging its characteristic polynomial according to powers of z and w, we see that the term \(P_{1,-1,\gamma }(\lambda )\), that is the polynomial with the prefactor \(w {\overline{z}}\), consists only of the monomial \(- \gamma _2 \gamma _3 \gamma _4 \gamma _5\). This can never be zero, so there cannot be a flat band.
A.6. \((3^2.4.3.4)\) A fundamental cell is given by four vertices with 10 edge weights to choose. The weighted adjacency \(\Pi _\gamma \) has the form
Arranging its characteristic polynomial according to powers of z and w, we see that the term \(P_{0,2,\gamma }(\lambda )\), corresponding to the prefactor \(z^2\), consists only of the monomial \(\gamma _1 \gamma _2 \gamma _9 \gamma _{10}\). This can never be zero, so there cannot be a flat band.
A.7. (3.4.6.4) A fundamental cell is given by six vertices with 12 edge weights to choose. The weighted adjacency \(\Pi _\gamma \) has the form
Arranging its characteristic polynomial according to powers of z and w, we see that the term \(P_{2,-2,\gamma }(\lambda )\), that is the polynomial with the prefactor \(w^2 {\overline{z}}^2\), consists only of the monomial \(\gamma _1 \gamma _2 \gamma _5 \gamma _{10} \gamma _{11} \gamma _{12}\). This can never be zero, so there cannot be a flat band.
A.8. (4.6.12) A fundamental cell consists of 12 vertices with 18 edge weights to choose. The weighted adjacency matrix \(\Pi _\gamma \) has the form
Arranging its characteristic polynomial according to powers of z and w, we see that the term \(P_{2,0,\gamma }(\lambda )\), that is the polynomial with the prefactor \(w^2\), consists only of the monomial \(\gamma _2 \gamma _3 \gamma _5 \gamma _6 \gamma _7 \gamma _8 \gamma _9 \gamma _{11} \gamma _{13} \gamma _{14} \gamma _{15} \gamma _{16}\). This can never be zero, so there cannot be a flat band.
A.9. \((3^4.6)\) A fundamental cell consists of 6 vertices with 15 edge weights to choose. The weighted adjacency matrix \(\Pi _\gamma \) has the form
Arranging its characteristic polynomial according to powers of z and w, we see that the term \(P_{2,0,\gamma }(\lambda )\), that is the polynomial with the prefactor \(w^2\), consists only of the monomial \(\gamma _1 \gamma _5 \gamma _9 \gamma _{10} \gamma _{11} \gamma _{13}\). This can never be zero, so there cannot be a flat band.
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.
About this article
Cite this article
Kerner, J., Täufer, M. & Wintermayr, J. Robustness of Flat Bands on the Perturbed Kagome and the Perturbed Super-Kagome Lattice. Ann. Henri Poincaré 25, 3831–3857 (2024). https://doi.org/10.1007/s00023-023-01399-7
Received:
Accepted:
Published:
Issue Date:
DOI: https://doi.org/10.1007/s00023-023-01399-7