The article surveys the known results and conjectures about the analytic properties of dispersion relations and Fermi surfaces for periodic equations of mathematical physics and their spectral incarnations.
This brief survey mostly follows and extends the outlines given in the author’s lectures at the conferences “Learning from Insulators New Trends in the Study of Conductivity of Metals” at the Lorenz Center in Leiden, Netherlands (August 2021) and “Algebraic Geometry in Spectral Theory” at ICERM, Brown University (February 2023).
Dispersion relations and Fermi surfaces for periodic operators of mathematical physics are some of the most common and important notions in condensed matter physics, in particular in the lately immensely popular nano-materials (graphene, nanotubes, etc.), photonic crystal theory, and topological insulators. We attempt here to survey their basic analytic properties, their relations to the spectra of the corresponding operators, their physical meaning, and some related results and conjectures.
Since the readership of this text might come from diverse communities of mathematicians and physicists, the author apologizes in advance for providing some discussions that might seem obvious to some, but not necessarily to others. So, for a mathematician, who has not dealt with dispersion relations, we start with analogy of constant coefficients partial differential equations/operators (PDEs). We then switch to honest periodic operators and introduce the corresponding definitions. After that we start discussing the analytic properties of these varieties and their role for properties of solutions and the spectrum.
The word “algebraic” in the title is justified by the fact that in the discrete situations (e.g., the popular tight binding approximation) the arising problems of analytic geometry become the ones of algebraic geometry instead. This is certainly an “improvement,” since one can use powerful results of algebraic geometry, sometimes not available in the analytic case, on the other hand, the discrete (graph) case sometimes presents surprises impossible in the continuous case, due to non-trivial graph topology and absence of the uniqueness of continuation property. A significant progress in the discrete case has been occurring lately.
Regretfully, it would be impossible to provide any proofs here. So, the reader is referred to the literature for those. Many references can be found in the survey.52 However, since a lot has happened in the past several years, we’ll try to present many more citations. As always, even with the best of intentions, some important references might be missed, and the author offers his apologies in advance.
II. CONSTANT COEFFICIENT PDEs AND FOURIER TRANSFORM
A. “Dispersion relation” and “Fermi surface” for constant coefficient operators
The quotation marks here reflect the fact that no one uses these words in the constant coefficient case. Well, we are in the free world and will use them.
- The graph of L(ξ) (in or )i.e., the graph of the symbol, is “dispersion relation,” or “Bloch variety.”(2.6)
- The level surfaceis “Fermi surface” at the “energy level” λ.(2.7)
Note that using not only real, but also complex values of ξ and λ turns out to be important, both in the constant coefficient and periodic situations.
The following statement is clear, because both of these varieties are given as sets of zeros of polynomials.
Both Bloch and Fermi varieties are algebraic.
Notice that the “Fermi” surface is just the set of zeros of the symbol [L(ξ) − λ]. It is also called the characteristic set for the operator L(D) − λ and is known to carry many “secrets” about the operator (see, e.g., Ref. 42).
B. Systems of constant coefficient PDEs
- “Dispersion relation,” or “Bloch variety” in or isi.e., is the graph of the multiple-valued function ξ ↦ σ(L(ξ)).
- The “Fermi surface” at the energy level λ:
III. PERIODICITY ENTERS THE FRY
A. Group of periods
Our starting point in the constant coefficient case was using the Fourier transform. This was natural, since such operators commute with the group acting (transitively) on by shifts. Then L(ξ) is the restriction of L to the 1D space generated by the function eiξ·x, which is how any irreducible representation (irrep) of Γ looks like. This is just the common spectral direct integral decomposition into irreps (see, e.g., Ref. 14).
We fix choices of a fundamental domain W (Wigner–Seitz cell) for action of Γ on and of a fundamental domain B (Brillouin zone) for action of Γ*. Functions that are Γ-periodic can be considered as functions on the torus . Analogously, Γ*-periodic functions can be identified with functions on the “Brillouin torus” . These tori are supplied with the natural flat metrics pushed down from the space . We will denote the corresponding normalized to the volume 1 measures by dx and dk correspondingly.
L(k) = L(x, D + k) acting on function on the torus .
The advantage of the operators L(k) in comparison with L is that they act on a compact manifold (torus) and if they are elliptic, their spectra are discrete.
- Since the correspondence between quasimomenta k and irreps of Γ is not one-to-one, it is often convenient to introduce the Floquet multiplierIf k is real, z is in “Brillouin torus” .(3.5)
One can observe that in the discrete case, in the Floquet multipliers rather than quasimomenta representation, the Bloch and Fermi varieties become algebraic.
Before going further, the reader should be warned that the irrep expansions in the (often important) case of co-compact action of a discrete non-abelian (=non-commutative) group98 Γ do not seem too useful, due to complicated structure of the space of the irreducible representations, and in particular due to usage of non-unitary representations (corresponding to the complex quasi-momenta k) in the abelian case. It has been understood that the better language is the one of group algebras (not providing any complete reference, we just mention as token examples Refs. 38, 39, and 66).
IV. THE CLASS OF OPERATORS
Let us describe precisely the classes of operators under consideration. Here are main assumptions (although sometimes we will need to weaken or strengthen them):
- where the coefficients aα are Γ-periodic and smooth99(4.1)
- The operator is assumed to be elliptic, i.e., its principal symboldoes not vanish for any real non-zero vectors ξ.(4.2)
This is a crucial assumption, without which the tools and results of the whole analytic theory of Bloch and Fermi varieties essentially fall apart or have to be changed significantly.100 This applies for instance for time-periodic wave or Schrödinger type equations (see for some discussions, applications, and Refs. 22, 34, 50, 52, 80, 88, 90, and 92 and references therein).
It will be often assumed that the operator is self-adjoint in , although some basic facts, like for instance analyticity of the varieties, hold without this condition.
In some results one will need to require that the operator is of the second order m = 2. Indeed, some of the crucial statements do not hold for general periodic elliptic operators of higher order, due to possible absence of unique continuation property.50,52,75 The same issue is encountered in the periodic discrete or quantum graph cases.52
Some other (matrix) operators of mathematical physics, such as Maxwell or Dirac ones are also important and studied (see e.g., Refs. 7, 8, 11–13, 51, 53, and 68), but we will not address these in the current text.
V. AND FINALLY, BLOCH AND FERMI VARIETIES!
The reader should note that the above definitions of the Bloch and Fermi varieties do not specify whether the quasimomenta k and eigenvalues λ have to be real. One can allow them to be complex without a problem, to obtain the complex versions .
VI. ANALYTICITY OF BLOCH AND FERMI VARIETIES
We have seen in the case of constant coefficient equations and systems of those that the Bloch and Fermi varieties were “nice,” namely algebraic sets [i.e., were the sets of zeros of the polynomial with respect to (k, λ) function det(L(ξ) − λI)]. It would have be nice if we could have a similar description in the periodic case, something like “det(L(k) − λI) = 0.” However, our operators L(k) are unbounded for each k. Here is where ellipticity gets into the game. The operators L(k) now have discrete spectra, but a simple-minded attempt to define some “determinants” for them fails. One can use, however, the so called regularized determinants (see Refs. 17, 37, and 50), which save the day and enable one to prove the following result.50,52
Our analogy with the constant coefficient case is almost achieved, except that the varieties of interest happen to be analytic, rather than algebraic sets. In other words, they can be described as the sets of zeros of entire functions, rather than polynomials. In some questions this does not make that much of a difference, while some useful properties of algebraic sets just do not survive for the analytic ones. In particular, the following “almost all or almost nothing” property of algebraic sets comes handy: projection of an algebraic subset either is contained in a proper algebraic subset (i.e., is “almost nothing”), or contains the complement to such a subset (i.e., is almost “everything”). Regretfully, the zero sets of entire functions might not have this property. Fortunately, having estimates of the sort of (6.3) helps19,62 and raises some hopes.
VII. THE SPECTRUM
The spectrum σ(L) is the projection on the spectral (i.e., λ-) axis of the real Bloch variety .
The spectrum consists of the spectral bands Ij, which are the ranges of the band functions λj(k) and thus are closed finite intervals.
When j → ∞, both endpoints of Ij tend to infinity.
In dimension higher than 1, the spectral bands might overlap.
The spectral bands might leave some uncovered spectral gaps in between. The spectral gaps ends (edges) correspond to extrema of the dispersion relation λ(k).
VIII. BLOCH AND FERMI VARIETIES AND SPECTRUM
A. Geometry of the spectrum–bands and gaps
According to the statement (1) of Theorem 7.1, the Bloch variety (dispersion relation) determines the spectrum and tells us that it has the band-gap structure. A very important questions is whether the spectral gaps do exist. Fort instance, existence of semi-conductors and photonic crystals depends upon this property.1,43,51
The spectrum is usually purely absolutely continuous. At least, this expected to be true for all second order elliptic periodic operator. Absence of singular continuous spectrum is often a rather mundane conclusion of periodicity,52 and is in fact much more general - it holds for a class of analytically fibered operators,35 which our direct integrals provide. A different story is presence (or rather absence) of the pure point spectrum, i.e., of bound states.
The spectrum of any periodic elliptic second order operator with sufficiently “nice” coefficients is absolutely continuous, and thus no bound states exist.
This result has been proven for Schrödinger operators (4.3) almost half of a century ago in Ref. 89 (see also expositions in Refs. 50, 52, and 76). Proving this in presence of a periodic magnetic potential turned out to be much harder, took a long time, and was finally achieved in Ref. 87. Handling variable coefficients in the second order terms would not succumb to the previously used harmonic analysis tools. This has been achieved in Ref. 32, where a seemingly irrelevant symmetry conditions was imposed, but no one has succeeded to remove it.
Why the second order restriction? The reason is that there are examples of periodic elliptic operators of higher order (of elasticity type), which do have bound states. The reason is the absence of weak uniqueness of continuation property, and thus existence of compactly supported solutions (see discussion in Ref. 52). The role of the (weak) uniqueness continuation principle102 is emphasized by the second order operator example provided in Ref. 27, in which the coefficients in the leading term are just a notch below those guarantying uniqueness of continuation, and indeed, compactly supported solutions (and thus bound states) do appear. Graphs (discrete and quantum) notoriously lack uniqueness of continuation, and thus bound states can exist. Moreover, in this case if a bound state does exist, then the corresponding eigenspace is generated by compactly supported solutions.5,49,52 If such a result could be proven in the continuous case, it would have proven absolute continuity of the spectrum in one shot, due to uniqueness of continuation.
Another indication here is the following result.50,52
Existence of such a solution for second order periodic (or more general class) elliptic equations should be prohibited (at least in the self-adjoint case) by some kind of “uniqueness of continuation at infinity” (the Landis problem). Thus, the spectra must be absolutely continuous.
The only results in this direction known to the author are in Refs. 33 and 67. They apply only to Schrödinger operators with bounded electric potentials and also allow r to reach up to 4/3 (at least in the non-selfadjoint case). There is significant activity nowadays around the Landis problem, but so far none of the results produce progress in establishing the above conjecture.
An interesting and very useful observation from Ref. 89 (although it is not formulated there in such terms) is the following.
This allows one to prove absence of bound states by smartly choosing a quasimomentum value k0 with large imaginary part, so that the point k0 is not in .
B. Gaps existence and creation
Fortunately, the Nature has provided us with semi-conductors, but the story of creating photonic spectral gaps was not an easy one (both experimentally and theoretically).
What could be techniques for creating spectral gaps? For a Schrödinger operator spectral gaps can be created near the bottom of the spectrum using the following simple procedure: adding to the free operator (Laplacian) a localized deep potential well will produce an eigenvalue below the continuous spectrum. Now, repeating this well periodically at a sufficiently large distance will spread this eigenvalue into a small spectral band, thus still leaving a gap.
This procedure does not work for Maxwell operators, where there is no way to shift the bottom of the spectrum. This was exactly what the struggle for creating photonic crystals was about.
Another procedure is of creating periodically placed high contrast inclusions throughout the medium (a version of the above procedure). It is known that this can create gaps throughout the spectrum (see, e.g., Refs. 23–25, 41, 86, and 96). The idea is to spread throughout the medium small identical resonators with internal structure, so near their eigenfrequencies the propagation is suppressed. Such gaps are often called resonant. This idea goes back at least to Ref. 74, but its explicit implementations has been achieved so far only in the graph cases81 and quantum graph.5,16,69 The gap opening is obtained by “decorating the graph: either by” attaching to every graph vertex another graph (“kite”)81 or inserting an internal structure (“spider”) into each vertex,16 see Fig. 3 (compare with Ref. 77). The spider decorations look in some sense more practical, but they do not necessarily open gaps and thus is harder to handle.103
The old Bethe-Sommerfeld conjecture (BSC)6 claims that in dimension higher than 1, the number of spectral gaps is expected to be finite. After a long series of works of distinguished mathematicians, this was proven for the Schrödinger operator case in Ref. 71, where and in Ref. 73 one can find further references.
Simple counterexamples to BSC are easy to construct (by decorations) for periodic discrete graphs, although for the quantum graph case the conjecture is still open.
C. Gap edges location
For quite some time there has been a belief that the spectral gap edges must occur only at some “symmetric” values of quasimomenta in Brillouin zone. This assumption allows for significant dimension reduction when computing the spectrum as the set (while not resolving the density of the spectrum). This is true and has simple origin in the one-dimensional case. There is no such reason in higher dimension, and indeed, this has been disproved in Refs. 20 and 40. There are, however, cases when this does hold true.3
D. Wave packets behavior
It is well known that the propagation of narrow wave packets of Bloch waves, located at a point of the dispersion relation, is governed for some length of time by the constant coefficient operator whose symbol is the jet (a lower degree part of the Taylor expansion) of the dispersion relation at this point1,21 (kind of quasi-classical limit).
There are two especially important cases.
1. Spectral gap edge
The first one is at the edge of a spectral gap, i.e., at an extremum of the dispersion relation k ↦ λ.
Generically with respect to the parameters of the (elliptic periodic) operator, e.g., for “almost all” periodic potentials, the edges of the spectral gaps are isolated in the k-space and non-degenerate (i.e., the dispersion relation has a parabolic shape).
This condition claims that near a spectral edge the solutions behave like in the free (albeit anisotropic) space. This allows one to define the effective masses.1 It is also important for establishing asymptotics of the Green’s function near the gap edges, Liouville theorems, Anderson localization, and other issues.45,46,52,54–56
2. Dirac cones
Another very popular nowadays case is the appearance of (stable) Dirac cones (“diabolic points”) in the dispersion relation.
Existence of such a cone implies that (for some time) wave packets located near the cone behave like solutions of the Dirac’s equation for massless fermions.22 The famous graphene owes many of its exciting properties to existence of these cones.44 The mandatory appearance of them in the honeycomb geometry is established in Ref. 22, see also the nice group-theoretic justification in Ref. 4.
IX. AN ESOTERIC QUESTION: IRREDUCIBILITY
Let be an analytic (algebraic) variety, i.e., it can locally be represented as the set of all common zeros of finitely many analytic (polynomial) functions fj(z).
Such a set is reducible, if it can be represented as the union of two strictly smaller analytic (algebraic) subsets Z1, Z2 ⊂ Z. One should imagine the case when f(z) has a non-trivial factorization f1(z)f2(z) into the product of two analytic (polynomial) functions fj.
An analytic (algebraic) set is irreducible if it is not reducible. Another incarnation of the irreducibility is that the set of all its smooth points is connected.
Yet another view on irreducibility is that any small open patch of the variety uniquely determines the whole of it (by analytic continuation).
The following result was established in Ref. 48.
Let L = −Δ + V(x), where V(x) is continuous periodic potential in . Then the Bloch variety of L is irreducible.
This seems to be an esoteric statement. In fact, it is not. Indeed, it implies that
There are no bound states. I.e., the spectrum of L is purely absolutely continuous.
The Bloch variety shows an amazing rigidity: any small patch of a band function determines the dispersion relation. (Do not hope for doing this! It is a highly numerically unstable procedure.)
The proof of the first statement is simple. Indeed, as Theorem 8.4 claims, existence of a bound states implies existence of a flat Fermi surface . This component shows that the Bloch variety is reducible into two component: one is and another is the rest of the Bloch variety.
The following conjecture is believed, but still unproven (except the result mentioned above).
The Bloch variety of any periodic Schrödinger operator −Δ + V(x) (or even general second order elliptic periodic operator) is irreducible.
In discrete case, some instances of this conjecture have been proven in Ref. 36.
The analogous question of irreducibility has been asked about the Fermi surfaces (for all, or almost all values of the energy λ).
This irreducibility problem also seems esoteric. However, as it was discovered in Refs. 57–59, it is at the heart of proving the absence of embedded states, which is a long standing problem in dimensions 2 and higher. Namely, when a periodic medium is perturbed by a localized perturbation, the old important question is whether impurity eigenvalues can arise inside the spectral bands, rather than in the spectral gaps. This fact has been established long ago in 1D78,79 (see also Refs. 9, 29, 30, and 93–95). To the best of the author’s knowledge, there had been no progress in higher dimensions, till it was discovered in Refs. 57–59 that irreducibility of the Fermi surface for a given λ inside of the continuous spectrum prevents appearance of embedded impurity eigenvalues at this point. Examples of irreducibility were presented there as well. Even before these works, conditions of Fermi surface irreducibility in the 2D square lattice case were found in Ref. 36. Significant recent progress proving irreducibility for discrete models, as well as examples of reducibility, and constructions of the corresponding examples of appearance of embedded eigenvalues have been provided in Refs. 26, 31, 63–65, 70, 82, and 83.
P.K. acknowledges the support from the NSF DMS Grant No. 2007408. Thanks also go to many colleagues for their fruitful discussion and collaboration. This includes, but not limited to G. Berkolaiko, Ngoc Do, J. Fillman, L. Friedlander, W. Liu, B. Ong, R. Matos, S. Shipman, F. Sottile, and J. Zhao. Thanks also go to the reviewer.
Conflict of Interest
The author has no conflicts to disclose.
Peter Kuchment: Conceptualization (lead); Formal analysis (lead); Investigation (lead); Methodology (lead); Writing – original draft (lead); Writing – review & editing (lead).
Data sharing is not applicable to this article as no new data were created or analyzed in this study.