This paper was converted on www.awesomepapers.org from LaTeX by an anonymous user.
Want to know more? Visit the Converter page.

Infinite Bar-Joint Frameworks, Crystals and Operator Theory

J.C. Owen D-Cubed, Siemens PLM Software, Park House,
Castle Park, Cambridge UK
owen.john.ext@siemens.com
 and  S.C. Power Dept. Math. Stats.
Lancaster University
Lancaster LA1 4YF
U.K.
s.power@lancaster.ac.uk
Abstract.

A theory of flexibility and rigidity is developed for general infinite bar-joint frameworks (G,p)(G,p). Determinations of nondeformability through vanishing flexibility are obtained as well as sufficient conditions for deformability. Forms of infinitesimal flexibility are defined in terms of the operator theory of the associated infinite rigidity matrix R(G,p)R(G,p). The matricial symbol function of an abstract crystal framework is introduced, being the matrix-valued function on the dd-torus representing R(G,p)R(G,p) as a Hilbert space operator. The symbol function is related to infinitesimal flexibility, deformability and isostaticity. Various generic abstract crystal frameworks which are in Maxwellian equilibrium, such as certain 44-regular planar frameworks, are proven to be square-summably infinitesimally rigid as well as smoothly deformable in infinitely many ways. The symbol function of a three-dimensional crystal framework determines the infinitesimal wave flexes in models for the low energy vibrational modes (RUMs) in material crystals. For crystal frameworks with inversion symmetry it is shown that the RUMS appear in surfaces, generalising a result of F. Wegner [32] for tetrahedral crystals.

1991 Mathematics Subject Classification:
52C75, 46T20.
Key words and phrases. infinite bar-joint framework, vanishing flexibility, rigidity operator.

1. Introduction

Infinite bar-joint frameworks appear frequently as idealised models in the analysis of deformations and vibration modes of amorphous and crystalline materials. See [11], [18], [6], [14], [32] and [36] for example and the comments below. Despite these connections there has been no extended mathematical analysis of such models. Notions of rigidity, flexibility, deformability, constrainedness, independence and isostaticity, for example, are usually employed either in the sense of their usage for a finite approximating framework or in a manner drawn from experience and empirical fact in the light of the application at hand. It seems that a deeper understanding of the models is of considerable interest in its own right and that a mathematical development may prove useful in certain applications. In what follows we shall provide formal definitions of the terms above in quite a wide variety of forms and we examine some of their inter-relationships and manifestations.

Suppose that one starts with a flexible square bar-joint framework in two dimensions and that this is then extended periodically to create an infinite periodic bar-joint network. Is the resulting assemblage, with inextendible bars, continuously flexible in two dimensional space? A moment’s reflection reveals a proliferation of flexibility, such as sheering motions with one half of the network fixed.

However such movement is dramatically infinite and a natural second question is whether for periodic frameworks such as these there are flexes for which the total joint movement is finite. The less obvious negative answer in these cases offers some satisfaction in that it is consistent with so-called Maxwell counting in the case where the average number of degrees of freedom of the joints matches the average number of constraints per joint. (See Theorem 5.2 and Corollary 5.3.)

On the other hand if, for the grid example, one rigidifies alternate squares by adding diagonal bars, as shown, then the resulting structure of corner linked rigid squares remains properly flexible, although now uniquely so. In fact the unique flex has an affinely contracting character (See Definition 4.3) with alternating rotation of the squares.

More generally the flexibility of polytope networks in two and three dimensions continues to be of interest in the modeling of crystals and amorphous materials, especially with regard to their low frequency vibrational modes. Such modes appear, for example, in higher order symmetry phases of tetrahedral crystals and are referred to as rigid unit modes (RUMs). Indeed in the paper of Giddy et al [11] the alternating flex of the squares framework above has been associated with vibrational modes in perovskite. See also Hammond et al [18], Wegner [32], as well as Goodwin et al [14] for a useful overview. At the same time, in the modeling of amorphous materials, such as glasses, there is interest in understanding the critical probabilities that guarantee flexibility and rigidity for classes of randomly constructed frameworks. See, for example, Chubynsky and Thorpe [6] for the recent determination of such probabilities in simulation experiments.

Refer to caption
Figure 1. The grid framework in the plane, 𝒢2{\mathcal{G}}_{{\mathbb{Z}}^{2}}.
Refer to caption
Figure 2. The corner-joined squares framework, 𝒢sq{\mathcal{G}}_{sq}.

Formally, a framework in d{\mathbb{R}}^{d} (or bar-joint framework, or distance-constraint framework) is a pair 𝒢=(G,p){\mathcal{G}}=(G,p) where G=(V,E)G=(V,E) is a simple connected graph and p=(p1,p2,)p=(p_{1},p_{2},\ldots) is a framework vector made up of framework points pip_{i} in d{\mathbb{R}}^{d} associated with the vertices v1,v2,v_{1},v_{2},\dots of VV. The framework edges are the (closed) line segments [pi,pj][p_{i},p_{j}] associated with the edges EE of the graph G=(V,E)G=(V,E). As the ellipsis suggest, we allow GG to be a countable graph. We shall also define a crystal framework 𝒞{\mathcal{C}} as a framework with translational symmetry which is generated by a connected finite motif of edges and vertices. (See Definition 4.2.)

When GG is finite and the framework points are generically located in 2\mathbb{R}^{2} then a celebrated theorem of Laman [21], well-known in structural engineering and in the discrete mathematics of rigidity matroids [15], gives a simple combinatorial criterion for the minimal infinitesimal rigidity of the framework; the graph itself satisfies Maxwell’s counting rule 2|V||E|=32|V|-|E|=3, and subgraphs G=(V,E)G^{\prime}=(V^{\prime},E^{\prime}) must comply with 2|V||E|32|V^{\prime}|-|E^{\prime}|\geq 3. This is a beautiful result since the rigidity here is the noncombinatorial requirement that the kernel of an associated rigidity matrix R(G,p)R(G,p) has the smallest dimension (namely three) for some (and hence all) generic framework. On the other hand frameworks with global symmetries, or even with ”symmetric elements” (such as parallel edges) are not generic, that is, algebraic dependencies do exist between the framework point coordinates. Such frameworks arise in classical crystallography on the one hand and in mathematical models in structural engineering and in materials science on the other. See, for example, Donev and Torquato [7], Hutchinson and Fleck [19], Guest and Hutchinson [16] and various papers in the conference proceedings [31].

The present paper develops two themes. The first concerns a mathematical theory of deformability and rigidity for general infinite frameworks, with frequent attention to the case of periodic frameworks. There is, unsurprisingly, a great diversity of infinite framework flexing phenomena and we introduce strict terminology and some methods from functional analysis to capture some of this. In the second theme we propose an operator theory perspective for the infinitesimal (first order) flexibility of infinite frameworks.

Particularly interesting classes of infinite frameworks, from the point of view of flexibility, are those in the plane whose graphs are 44-regular and those in three dimensions whose graphs are 66-regular. In this case the graphs are in Maxwell equilibrium, so to speak, and so in a generic framework realisation any flex must activate countably many vertices. This is also the case for various periodic realisations such as the kagome framework, 𝒢kag{\mathcal{G}}_{kag}, formed by corner-joined triangles in regular hexagonal arrangement, and frameworks in three dimensions formed by pairwise corner-joined tetrahedra. Despite being internally rigid in this way (Definition 2.18 (vi)) these frameworks admit diverse deformations. For example we note that the kagome framework admits uncountably many distinct deformations and in Theorem 4.4 we note that d{\mathbb{Z}}^{d}-periodic cell-generic grid frameworks in d{\mathbb{R}}^{d} admit deformations associated with affine transformations.

A significant phenomenon in the infinite setting is the appearance of vanishing flexibility. This means, roughly speaking, that the framework is a union of finite flexible subframeworks but the extent of flexibility diminishes to zero as the size of these subframeworks increases, so that the infinite assemblage is inflexible. Elementary examples were indicated in [25] but we give more subtle examples here which are due to flex amplification at second order distances through concatenation effects. In particular there are bounded infinitesimal flexes in periodic frameworks that admit no continuous extensions and which do not arise as the derivative of a smooth deformation. We also note that there are 2{\mathbb{Z}}^{2}-periodic crystal frameworks which are somewhat paradoxical, being indeformable despite the flexibility of all supercell subframeworks. On the other hand, in the positive direction, in Theorem 2.20 we give a general result which identifies a uniform principle for the existence of a deformation. The proof uses the Ascoli-Arzela theorem on the precompactness of equicontinuous families of local flexes. It remains an interesting open problem to determine necessary and sufficient conditions for the rigidity and bounded rigidity of periodic planar frameworks.

The operator theory perspective for frameworks was suggested in [25] as an approach to a wider understanding of infinitesimal flexibility and rigidity. In this consideration the rigidity matrix is infinite and determines operators between various normed sequence spaces associated with nodes and with edges. Also, in [27] we have given a direct proof of the Fowler-Guest formula [9] for symmetric finite frameworks which is based on the commutation properties of the rigidity matrix as a linear transformation and this adapts readily to the infinite case and the rigidity operators of crystal frameworks. Indeed, translational symmetry ensures that the rigidity matrix R(G,p)R(G,p) intertwines the coordinate shift operations. We consider square summable flexes and stresses and for distance regular locally finite frameworks R(G,p)R(G,p) is interpreted as a bounded operator between Hilbert spaces. Also, enlarging to complex Hilbert spaces the Fourier transform R(G,p)1{\mathcal{F}}R(G,p){\mathcal{F}}^{-1} is identified as a multiplication operator

MΦ:L2(𝕋d)nL2(𝕋d)mM_{\Phi}:L^{2}({\mathbb{T}}^{d})\otimes{\mathbb{C}}^{n}\to L^{2}({\mathbb{T}}^{d})\otimes{\mathbb{C}}^{m}

given by an m×nm\times n matrix-valued function Φ(z)\Phi(z) on the dd-torus. The function Φ\Phi for 𝒞{\mathcal{C}} is referred to as the matricial symbol function associated with the particular generating motif. The terminology and notation is borrowed from standard usage for Toeplitz operators and multiplication operators (see [3] for example). Many aspects of infinitesimal flexibility and isostaticity are expressible and analysable in terms of the matricial symbol function and its associated operator theory. For example a straightforward consequence of the operator theoretic approach is the square summable isostaticity of various nondegenerate regular frameworks that satisfy Maxwell counting, such as grid frameworks and the kagome framework.

An explicit motif-to-matrix function algorithm is given for the progression

𝒞=(G,p)R(G,p)Φ(z).{\mathcal{C}}=(G,p)\to R(G,p)\to\Phi(z).

Furthermore the identification of infinitesimal periodic-modulo-phase flexes and their multiplicities is determined by the degeneracies of Φ(z)\Phi(z) as zz ranges on the dd-torus. In particular, the function

μ(z):=dimkerΦ(z):𝕋d.\mu(z):=\dim\ker\Phi(z):{\mathbb{T}}^{d}\to{\mathbb{Z}}_{.}

gives a determination of the mode multiplicity of periodic-modulo-phase infinitesimal flexes.

In the discussions below we are mainly concerned with properties of mathematical bar-joint frameworks. (The framework bars are indestructibly inextensible, the joints are located deterministically, they maintain perfect, frictionless fit and may even coincide.) Nevertheless, analysis of matricial symbol functions and their degeneracies are particularly relevant to the description and analysis of Rigid Unit Modes in material crystals. We show that for crystal frameworks with inversion symmetry the set of RUMS is a union of surfaces. This generalises and provides an alternative perspective for a recent result of Wegner [32] for tetrahedral crystals.

Operator theory methods have proven beneficial in many areas of mathematics and applications. This is evidently the case for multi-linear systems theory and in control theory for example. Infinite rigidity matrix analysis seems to possess some similitudes with these areas and it seems to us that here too the operator turn will be a useful one.

The development is as follows. Section 2 gives a self-contained account of continuous flexibility and rigidity, and vanishing flexibility and one-sided flexibility is proven for various periodic infinite frameworks. Forms of flexibility, such as bounded flexes, square summable flexes, summable flexes and vanishing flexes are defined and determined for some specific examples. Sufficient conditions are obtained for the existence of a smooth flex and a flex extension problem for generic finite frameworks is posed. A positive resolution of this problem would provide a natural extension of Laman’s theorem to infinite frameworks. In Section 3 we consider infinitesimal theory for general infinite frameworks and determine a number of rigidity operators and their flex and stress spaces. (The topic is taken up in more detail for crystal frameworks in Section 5.) In Section 4 we consider (abstract) crystal frameworks in two or three dimensions. These are generated by a motif and a discrete translation group. Various forms of deformations are considered, such as strict periodic flexibility, flow-periodic flexibility and flexes with reduced periodicity and symmetry. Also we indicate the flat torus model for crystal frameworks and recent results of Ross and Whiteley in this direction, including a periodic analogue of Laman’s theorem. In the final section we consider the matricial symbol function approach and various examples. In particular we determine the (unit cell) infinitesimal wave flex multiplicities for the kagome net framework by factoring the determinant of the matricial symbol function.

Acknowledgments. Some of the developments here have benefited from discussions and communications with Robert Connelley, Patrick Fowler, Simon Guest, Elissa Ross and Walter Whiteley during and following the Summer Research Workshop on ”Volume Inequalities and Rigidity”, organised by Károly Bezdek, Robert Connelly, Balázs Csikös and Tibor Jordán, at the Department of Geometry, Institute of Mathematics, Eotvos Lorand University in July 2009.

2. Infinite Bar-joint Frameworks

In this section we give a self-contained rigourous development of infinite frameworks and examine the nonlinear aspects of their flexibility by continuous deformations and their associated rigidity. In the next section we consider infinitesimal flexibility and rigidity in a variety of forms.

2.1. Continuous flexibility and rigidity

We first define continuous flexes and continuous rigidity. The latter means, roughly speaking, that the framework admits no proper deformations that preserve the edge lengths. The definition below gives straightforward generalisations of terms used for finite frameworks. In that case we note that a continuous flex is often referred to as a finite flex while in engineering models it is referred to as a finite mechanism.

Unless we specify otherwise we shall assume that the frameworks under consideration are proper in that the framework points do not lie on a hyperplane in the ambient space d{\mathbb{R}}^{d} and that the framework edges [pi,pj][p_{i},p_{j}] have nonzero lengths |pipj||p_{i}-p_{j}|.

Definition 2.1.

Let (G,p)(G,p) be an infinite framework in 2{\mathbb{R}}^{2}, with connected abstract graph G=(V,E)G=(V,E), V={v1,v2,}V=\{v_{1},v_{2},\dots\} and p=(p1,p2,)p=(p_{1},p_{2},\dots).

(a) A base-fixed continuous flex, or, simply, a flex of (G,p)(G,p), is a function
p(t)=(p1(t),p2(t),)p(t)=(p_{1}(t),p_{2}(t),\dots) from [0,1][0,1] to V2\prod_{V}{\mathbb{R}}^{2} with the following properties;

(i) p(0)=pp(0)=p,

(ii) each coordinate function pi:[0,1]2p_{i}:[0,1]\to{\mathbb{R}}^{2} is continuous,

(iii) for some base edge (va,vb)(v_{a},v_{b}) with |papb|0|p_{a}-p_{b}|\neq 0, pa(t)=pa(0)p_{a}(t)=p_{a}(0) and pb(t)=pb(0)p_{b}(t)=p_{b}(0) for all tt,

(iv) each edge distance is conserved; |pi(t)pj(t)|=|pi(0)pj(0)||p_{i}(t)-p_{j}(t)|=|p_{i}(0)-p_{j}(0)| for all edges (vi,vj)(v_{i},v_{j}), and all tt, and

(v) p(t)pp(t)\neq p for some t(0,1]t\in(0,1].

(b) The framework (G,p)(G,p) is flexible, or more precisely, continuously flexible, if it possesses a base-fixed continuous flex.

(c) The framework (G,p)(G,p) is rigid, or continuously rigid, if it is not flexible.

Similarly one defines base-fixed continuous flexes and continuous rigidity for proper frameworks in d{\mathbb{R}}^{d} by replacing a base edge by an appropriate set of framework points with maximal affine span.

The simplest kind of continuously rigid framework in the plane is one which is a union of continuously rigid finite frameworks. In particular the following theorem follows simply from the theorem of Laman indicated in the introduction.

Theorem 2.2.

Let GG be a connected graph which is the union of a sequence of finite Laman graphs. Then every generic realisation (G,p)(G,p) in the plane is continuously rigid.

By generic, or, more precisely, algebraically generic, we mean, as is usual, that the coordinates of any finite set of framework points is algebraically independent over the rational numbers. Unlike the case of finite frameworks it is possible to construct two generic frameworks with the same abstract graph one of which is flexible and one of which is rigid. Accordingly it seems appropriate to formulate the following definition to extend the usual usage.

Definition 2.3.

An infinite simple connected graph GG is said to be rigid, or generically rigid, for two dimensions, if every generic framework (G,p)(G,p) in the plane is rigid.

Note that if GG is rigid and if HH is a containing graph for GG with the same vertex set then every generic framework (H,p)(H,p) in the plane is rigid.

It seems likely that the converse to the theorem above holds. That is, if HH does not contain a sequentially Laman graph (in the sense below) with the same vertex set, then HH is not generically rigid. We comment more on this later in Section 2.6.

Rigidity and flexibility are properties of the entire framework and it is such entire features and their inter-relationships that are of primary interest in what follows. One would like to understand the relationship with small scale or local structure, such as local counting conditions and local connectivity. Additionally, as above, one would like to relate entire properties to sequential features that pertain to an exhausting chain of finite subframeworks and for this the following definition is helpful.

Definition 2.4.

If PP is a property for a class of finite, simple, connected graphs (resp. frameworks) then a graph GG (resp. framework 𝒢=(G,p){\mathcal{G}}=(G,p)) is sequentially PP or σ\sigma-PP if GG is the union of graphs in some increasing sequence of vertex induced finite subgraphs G1G2,G_{1}\subseteq G_{2}\subseteq\dots, and each graph GkG_{k} (resp. framework (Gk,p)(G_{k},p)) has property PP.

For example, we may refer to an infinite graph as being σ\sigma-Laman, or σ\sigma-Laman1-1 and an infinite framework as being σ\sigma-rigid. To say that an infinite framework is σ\sigma-flexible, or sequentially flexible, is rather vacuous since it usually prevails for trivial reasons. (One can construct countably infinite edge complete frameworks that may fail to be so but all our examples are sequentially flexible.) The more interesting property is the failure of sequentially rigidity as in the following definition.

Definition 2.5.

A framework 𝒢{\mathcal{G}} is said to have vanishing flexibility if it is continuously rigid but not σ\sigma-rigid.

An important topic in finite rigidity frameworks is that of global rigidity, also called unique rigidity. This holds, in the two dimensional setting, if there is, up to congruency, only one embedding of the framework in 2{\mathbb{R}}^{2}. This implies, for example, that these frameworks admit no foldings, and indeed have unique diagrams up to rotations and reflection. One can extend the term to infinite frameworks but we do not consider this issue here at all. One might be tempted to say that a rigid infinite framework, especially one with vanishing flexibility, is globally rigid, but we refrain from doing so because of conflict with this usage.

The term ”global” for global rigidity is natural since rigidity for finite frameworks is equivalent to the ”local” property that ”nearby equivalent frameworks are congruent”. That is, if there exists ϵ>0\epsilon>0 such that if (G,p)(G,p^{\prime}) is a finite framework with |pipi|<ϵ|p_{i}-p_{i}^{\prime}|<\epsilon, for all ii, and if (G,p)(G,p^{\prime}) is equivalent to (G,p)(G,p), in the sense that corresponding edges have the same length, then (G,p)(G,p^{\prime}) and (G,p)(G,p) are congruent. See, for example, Gluck [12] and Asimow and Roth [1].

An infinite simple graph GG is locally finite if for every vertex vv there are finitely many incident edges. Amongst such graphs are those for which there is an upper bound to the degree of the vertices, as in the case of the graphs of crystal frameworks. Within this class a graph GG is said to be rr-regular if every vertex has degree rr. We remark that the theory of tilings provides a wealth of examples of planar frameworks which are 44-regular.

Remark 2.6.

In what follows we consider only locally finite frameworks. Without this assumption it is possible to construct quite wildly flexing planar linkages. In fact, given a continuous function f:[0,1]2f:[0,1]\to{\mathbb{R}}^{2} one can construct an infinite linkage, in the sense of the definition below, and a base-fixed flex p(t)p(t) with a motion pv(t)p_{v}(t) for a particular vertex vv that is equal to f(t)f(t). This includes the possibility of space filling curves. This is a consequence of a continuous analogue of a well known theorem of Kempe which asserts that any finite algebraic curve in the plane can be simulated by a finite linkage. For more details see Owen and Power [26].

2.2. Linkages

The removal of a framework edge from a rigid framework may result in flexibility which is, roughly speaking, of a one-dimensional nature. We reserve the term linkage for such a mathematical object, which we formally specify in the next definition. We remark that finite frameworks are also referred to as linkages, particularly when they are flexible, perhaps with several degrees of freedom, but this should not cause confusion.

A two-sided continuous flex p(t)p(t) of (G,p)(G,p) is defined as above but for the replacement of [0,1][0,1] by [1,1][-1,1]. The following formal definition of an infinite linkage reflects the fact that the initial motion of a base-fixed linkage is uniquely determined by the angle change at any flexible joint.

Definition 2.7.

A linkage in 2{\mathbb{R}}^{2} is a finite or infinite connected framework 𝒢=(G,p){\mathcal{G}}=(G,p) in 2{\mathbb{R}}^{2} for which there exists a continuous two-sided base-fixed flex p(t)p(t) with framework edges [pi,pj],[pj,pk][p_{i},p_{j}],[p_{j},p_{k}] such that the cosine angle function

g(t)=pi(t)pj(t),pk(t)pj(t)g(t)=\langle p_{i}(t)-p_{j}(t),p_{k}(t)-p_{j}(t)\rangle

is strictly increasing on [1,1][-1,1], and such that p(t)p(t) is the unique two-sided flex q(t)q(t) of 𝒢\mathcal{G} with ql(t)=pl(t)q_{l}(t)=p_{l}(t), for l=i,j,kl=i,j,k.

Many interesting finite linkages were considered in the nineteenth century in connection with mechanical linkages. See, for example, Kempe [20]. Note however that the definition is liberal in allowing coincident joints and self-intersecting flexes. Also the definition refers to local deformation behaviour and this does not rule out the possibility of bifurcations occurring in a parameter extension of the given flex.

It is a simple matter to construct diverse infinite linkages by tower constructions or progressive assembly. (See, for example, the Cantor tree frameworks of [25].) However, some such constructions lead to frameworks with vanishing flexibility and so are not linkages in this case. An elementary illustration is given in Figure 4 wherein a two-way infinite rectangular strip linkage is augmented by adding flex-restricting cross braces in an alternating fashion. If the brace lengths tend to the diagonal length from above then the infinite framework is rigid. Evidently in this case the triangle inequality is playing a role in isolating one real solution to the solution set V(G,p)V(G,p) defined below. One can also construct examples in which this isolation is less evident, with all joint angles bounded away from zero and π/2\pi/2 for example.

A more interesting and subtle form of vanishing flexibility is due to progressive flex amplification rather than local flex restrictions. Roughly speaking, if a small flex is initiated at a particular joint and the flex propagates in some amplifying manner, then the triangle inequality at some far remove may prohibit any further increase. If the framework is infinite then no local joint flex may be possible at all. The strip framework of concatenated levers in Figure 5 gives an example where the amplification is evident, while the rigidity of the strip framework in Figure 6 and the trapezium strip in Figure 7 is less evident. The lever framework has a natural infintesimal flex, in the sense of Section 3, which is unbounded. The corresponding flexes of Figures 6 and 7 however, are bounded with amplification unfolding as a second order effect. This is proven in the next subsection.

Refer to caption
Figure 3. An infinite rectangle strip linkage.
Refer to caption
Figure 4. A restricted rectangle strip.
Refer to caption
Figure 5. Rigid but not σ\sigma-rigid.
Refer to caption
Figure 6. Rigid but not σ\sigma-rigid.
Refer to caption
Figure 7. Rigid but not σ\sigma-rigid.
Refer to caption
Figure 8. A periodic half-strip which is only right-flexible.

It is a straightforward matter to incorporate the vanishing flexibility of the strips above as subframeworks of a 2{\mathbb{Z}}^{2}-periodic framework. This process is indicated in Figures 9, 10 and 11 below, where the infinite frameworks are determined as the periodic extensions of the given unit cell. Figure 9 shows a linkage formed as a ”fence lattice” composed of infinite horizontal and vertical σ\sigma-rigid bands. Figure 10 shows an analogue where the infinite bands have been replaced by rigid strip frameworks. Figure 11 is an elaboration of this in which cross braces have been introduced to remove the flexibility. (Only one edge is needed for this whereas the example given is periodic.) The additional degree 2 vertex in the cell ensures that the framework is not σ\sigma-rigid, while the infinite bands remain vanishingly rigid.

Let us note that for the framework in Figure 11, with its curious mixture of rigidity and flexibility, one can add any finite number of additional degree 2 vertices without changing the rigidity of the framework. In particular we have a construction that proves the following proposition.

Refer to caption
Figure 9. Unit cell for a ”fence lattice” linkage.
Refer to caption
Figure 10. Unit cell for a modified fence lattice linkage.
Refer to caption
Figure 11. Unit cell for a rigid periodic framework which is not σ\sigma-rigid.
Proposition 2.8.

Let c>2c>2. Then there is a 2{\mathbb{Z}}^{2}-periodic framework in 2{\mathbb{R}}^{2} which is rigid, which is not σ\sigma-rigid and for which the average vertex degree is less than cc.

One can readily extend this fanciful idea in various ways to obtain such structures in higher dimensions. For example, start with a one-dimensionally periodic σ\sigma-rigid girder in 3D and augment it with trapezium ”tents” of alternating height to creates vanishingly rigid girders. Also periodically interpolate any number of degree two vertices into the tent top edges without removing the vanishing flexibility. Join infinitely many such component girders periodically at appropriate tent-top edges to create a fence framework and add linear jointed cross braces to create, finally, a 2D periodic grid which is continuously rigid in 3D, which is not σ\sigma-rigid and which has average coordination number arbitrarily close to two.

2.3. Relative rigidity and the extension of flexes.

For a finite or infinite framework 𝒢=(G,p){\mathcal{G}}=(G,p) in 2{\mathbb{R}}^{2} define the function

fG:V2E,fG(q)=(|qiqj|2)e=(vi,vj).f_{G}:\prod_{V}{\mathbb{R}}^{2}\to\prod_{E}{\mathbb{R}},\quad f_{G}(q)=(|q_{i}-q_{j}|^{2})_{e=(v_{i},v_{j})}.

This is the usual edge function of the framework and depends only on the abstract graph GG.

Definition 2.9.

The solution set, or configuration space, of a framework 𝒢=(G,p){\mathcal{G}}=(G,p), denoted V(G,p)V(G,p), is the set fG1(fG(p)).f_{G}^{-1}(f_{G}(p)). This is the set of all framework vectors qq for GG that satisfy the distance constraints equations

|qiqj|2=|pipj|2, for all edges e=(vi,vj).|q_{i}-q_{j}|^{2}=|p_{i}-p_{j}|^{2},\mbox{ for all edges }{e=(v_{i},v_{j})}.

In general the solution set of an infinite framework need not be a real algebraic variety even when it is ”finitely parametrised”. In less wild situations it can be useful to relate V(G,p)V(G,p) to the algebraic variety V(H,πH(p))V(H,\pi_{H}(p)) associated with a finite subgraph HH of GG, or with an elementary subgraph such as a tree, or even a set of vertices.

Definition 2.10.

An infinite bar-joint framework (G,p)(G,p) in d{\mathbb{R}}^{d} is rigid over a subframework (H,πH(p))(H,\pi_{H}(p)) if every continuous flex of (G,p)(G,p) which is constant valued on (H,πH(p))(H,\pi_{H}(p)) is constant. Similarly, if HH is a subgraph of a countable connected simple graph GG then GG is rigid over HH, or generically rigid over HH if, for every generic frameworks (G,p)(G,p), every continuous flex of (G,p)(G,p) which is constant-valued on (H,πH(p))(H,\pi_{H}(p)) is constant.

We may also form the following associated notions.

Definition 2.11.

An infinite framework (G,p)(G,p) in d{\mathbb{R}}^{d} is finitely determined if it is rigid over (H,p)(H,p) for some finite subgraph HH and is finitely flexible if it is flexible and finitely determined.

Finite flexibility in the sense above is a strong property in which paths from pp in the solution set V(G,p)V(G,p) are determined near pp by the finite algebraic variety V(H,p)V(H,p). Note that the term ”infinitely flexible” is not appropriate to describe a flexible framework which is not finitely flexible since it is possible to construct linkages, in our formal sense, which are not finitely flexible. This is the case for the periodic framework in Figure 12 which is a linkage because of partial vanishing flexibility.

Refer to caption
Figure 12. A linkage which is not finitely determined.

An important class of frameworks which appear in mathematical models are those that are distance regular.

Definition 2.12.

A framework 𝒢=(G,p){\mathcal{G}}=(G,p) in d{\mathbb{R}}^{d} is distance-regular if there exist 0<m<M0<m<M such that for all edges (i,j)(i,j),

m<|pipj|M.m<|p_{i}-p_{j}|\leq M.

For such a framework (G,p)(G,p) it is natural to consider the nearby frameworks with the same graphs but with slightly perturbed framework points (and therefore edge lengths). If a property holds for all such frameworks, for some perurbation distance ϵ\epsilon then we call such a property a stable property for the the framework.

Formally, an ϵ\epsilon-perturbation of a distance regular framework 𝒢=(G,p){\mathcal{G}}=(G,p) is a framework 𝒢=(G,p){\mathcal{G}}^{\prime}=(G,p^{\prime}) for which |pipi|<ϵ|p_{i}-p_{i}^{\prime}|<\epsilon, for all corresponding framework points. Recall that a finite framework in d{\mathbb{R}}^{d} is said to be ϵ\epsilon-rigid if it is congruent to every equivalent ϵ\epsilon-perturbation. Let us say that a general framework is perturbationally rigid if it is ϵ\epsilon-rigid for some ϵ\epsilon. It is a well-known fact that perturbational rigidity and rigidity are equivalent in the case of algebraically generic finite frameworks [1], [12]. However, it is straightforward to see that this equivalence thoroughly fails for general infinite frameworks (see [25]).

Definition 2.13.

Let 𝒢{\mathcal{G}} be a distance-regular framework. Then 𝒢{\mathcal{G}} is stably rigid (resp stably flexible) if it is rigid (resp. flexible) and for sufficiently small ϵ>0\epsilon>0 every ϵ\epsilon-perturbation of 𝒢{\mathcal{G}} is rigid (resp. flexible).

Likewise, if PP is any particular property of a distance-regular infinite framework then we may say that 𝒢{\mathcal{G}} is stably PP if, for some ϵ>0\epsilon>0, the property PP holds for all ϵ\epsilon-perturbations.

Proposition 2.14.

The periodic trapezium strip frameworks, with alternating unequal heights a>b>oa>b>o, are rigid. In particular the rectangle strip linkage (of Figure 7) is not stably flexible.

Refer to caption
Figure 13.
Proof.

Figure 13 shows the displacement of a double trapezium to the right. Let the three vertical bar lengths be a,ba,b and aa units with a>b>0a>b>0. The displaced position has angles A,B,CA,B,C at the base line and angles D,E,F,GD,E,F,G occurring relative to the trajectory tangents of the displaced vertices. For a subsequent incremental change δA\delta A, with resulting incremental changes δB,δC,δD,δE,δF,δG\delta B,\delta C,\delta D,\delta E,\delta F,\delta G we can see from simple geometry that to first order

a(δAcosD)=b(δBcosE).a(\delta A\cos D)=b(\delta B\cos E).

Suppose now that AA is regarded as a specialisation of the input angle α\alpha with resulting output angle β=β(α)\beta=\beta(\alpha), so that at α=A\alpha=A we have β(α)=B\beta(\alpha)=B. Then

dβdα|α=A=limδA0δBδA=abcosDcosE.\frac{d\beta}{d\alpha}|_{\alpha=A}=\lim_{\delta A\to 0}\frac{\delta B}{\delta A}=\frac{a}{b}\frac{\cos D}{\cos E}.

Similarly, with angle CC regarded as the output angle γ(B)\gamma(B) for the angle transmission function γ=γ(β)\gamma=\gamma(\beta), we have

dγdβ|β=B=bacosFcosG\frac{d\gamma}{d\beta}|_{\beta=B}=\frac{b}{a}\frac{\cos F}{\cos G}

and so

dγdα|α=A=dγdβ|β=Bdβdα|α=A=cosDcosEcosFcosG.\frac{d\gamma}{d\alpha}|_{\alpha=A}=\frac{d\gamma}{d\beta}|_{\beta=B}\frac{d\beta}{d\alpha}|_{\alpha=A}=\frac{\cos D}{\cos E}\frac{\cos F}{\cos G}.

Note that since B>AB>A we have also D>ED>E, and since B>CB>C we have F>GF>G, from which it follows that both ratios above are less than one. Thus certainly 0<γ(α)<10<\gamma^{\prime}(\alpha)<1 for 0<α<α10<\alpha<\alpha_{1} where α1\alpha_{1} is the first positive angle for which γ(α1)=0.\gamma^{\prime}(\alpha_{1})=0.

It follows, from the mean value theorem, that the double trapezium angle transmission function is an increasing differentiable function with

γ(0)=0,0<γ(α)<α,for0<α<α1.\gamma(0)=0,\quad 0<\gamma(\alpha)<\alpha,\quad\mbox{for}\quad 0<\alpha<\alpha_{1}.

It follows immediately that the right-semi-infinite trapezium strip is right flexible.

We let λ=γ(α1)\lambda=\gamma(\alpha_{1}), which we refer to as the locking angle. Note the second trapezium of the double admits no increase of this angle. In view of the above we have λ<α1\lambda<\alpha_{1}.

Suppose now that 𝒢{\mathcal{G}} is the two-way-infinite trapezium strip, with aba\neq b. Let p(t)p(t) be a flex and suppose that for a fixed framework edge with length aa the angle AA is greater than zero for some time t1>0t_{1}>0 and that t1t_{1} is the first such time. Then certainly 0<A<λ0<A<\lambda. Note that An=γn(A),n=1,2,A_{-n}=\gamma^{-n}(A),n=1,2,\dots are the angles of the edges of length aa, counted off to the left. In view of the function dominance 0<γ(α)<α0<\gamma(\alpha)<\alpha it follows that An>α1A_{-n}>\alpha_{1} for some nn, which is a contradiction. ∎

The argument above also shows that the semi-infinite trapezium strip framework of Figure 7 has a continuous flex but no two-sided flex.

Remark 2.15.

It seems to be of interest to analyse strip frameworks in further detail. For example, a trapezium strip framework is not stably rigid, despite the apparent ”robustness” of the argument above. To see this use surgery in the following way. Remove one cross bar, then push the rightmost semifinite strip to the right, by an angle perturbation A=ϵ>0A=\epsilon>0. Now insert a replacement bar of the required length. One can flex the resulting structure towards the left to restore the position of the right hand strip. Indeed, this is all the flexibility the framework has. The possible flex of an ϵ\epsilon perturbation, such as the one described, seems to be of order ϵ\epsilon and so there does seem to be ”approximate rigidity”.

Remark 2.16.

Consider a periodic trapezium grid framework 𝒢trap{\mathcal{G}}_{trap} obtained by perturbing 𝒢2{\mathcal{G}}_{{\mathbb{Z}}^{2}} by adding a fixed small positive value to the yy coordinate of the framework points pijp_{ij} for the odd values of ii and jj. It can be shown that this framework is rigid over any linear subframework. This contrasts with the grid framework itself which is freely flexible over its xx and yy axes in the following specific local sense: all sufficiently small flexes of the subframework extends to a flex of the whole framework. On the other hand note that Theorem 4.4 shows that 𝒢trap{\mathcal{G}}_{trap} is deformable.

2.4. Forms of flexibility

It seems to be a fundamental and interesting issue to determine the ways in which infinite bar-joint frameworks are rigid or continuously flexible. In this section we give some further definitions, we give sufficient conditions for the existence of a proper flex and we contemplate a plausible infinite framework version of Laman’s theorem.

Flexes are often infinitely differentiable or smooth in the sense of the following formal definition. This is the case for example, for the ”alternation” flexes of 𝒢sq{\mathcal{G}}_{sq} and 𝒢kag{\mathcal{G}}_{kag}.

Definition 2.17.

A continuous base-fixed two-sided flex p(t):t[1,1]p(t):t\in[-1,1] of a framework (G,p)(G,p) in d{\mathbb{R}}^{d} is a smooth flex if each coordinate function pi(t)p_{i}(t) is infinitely differentiable.

The smoothness of a flex is a local requirement whereas the following terms impose various increasing forms of global constraint. In particular rotational flexes of infinite frameworks with unbounded diameter are not bounded flexes, while a translational flex of an infinite framework is a bounded flex but is not a vanishing flex. Adopting a term that has been used in applications [14] we refer to flexes which are not bounded as colossal flexes.

Definition 2.18.

A continuous flex p(t)=(pk(t))k=1,(t[0,1])p(t)=(p_{k}(t))_{k=1}^{\infty},(t\in[0,1]) of an infinite framework (G,p)(G,p) in d{\mathbb{R}}^{d} is said to be

(i) a bounded flex if for some M>0M>0 and every kk and tt,

|pk(t)pk(0)|M,|p_{k}(t)-p_{k}(0)|\leq M,

(ii) a colossal flex if it is not bounded,

(iii) a vanishing flex if p(t)p(t) is a bounded flex and if the maximal displacement

pkpk(0)=supt[0,1]|pk(t)pk(0)|\|p_{k}-p_{k}(0)\|_{\infty}=\sup_{t\in[0,1]}|p_{k}(t)-p_{k}(0)|

tends to zero as kk\to\infty,

(iv) a square-summable flex if

k=1pkpk(0)2<,\sum_{k=1}^{\infty}\|p_{k}-p_{k}(0)\|_{\infty}^{2}<\infty,

(v) a summable flex if

k=1pkpk(0)<,\sum_{k=1}^{\infty}\|p_{k}-p_{k}(0)\|_{\infty}<\infty,

(vi) an internal flex if for all but finitely many kk the function pk(t)p_{k}(t) is constant.

Also we say that (G,p)(G,p) has a deformation (resp. bounded or vanishing deformation) if it has a base-fixed flex p(t)p(t) (which is bounded or vanishing).

Definition 2.19.

A connected infinite locally finite proper framework in two or three dimensions is boundedly rigid (resp. summably rigid, square-summably rigid, smoothly rigid, internally rigid) if there is no deformation, that is, no base-fixed proper continuous flex, which is bounded (resp. summable, square-summable, smooth, internal).

2.5. Sufficient conditions for flexibility

There is a sense in which vanishing flexibility is almost the only obstacle to the existence of a flex of a framework all of whose finite subframeworks are flexible. More precisely, in the hypotheses of the next theorem we assume that there are two distinguished framework vertices, p1,p2,p_{1},p_{2}, such that any finite subframework (H,p)(H,p) containing p1,p2p_{1},p_{2} has a flex which properly separates this pair in the sense of condition (ii) below. The additional requirement needed is that there is a family of flexes of the finite subframeworks whose restrictions to any given subframework (H,p)(H,p) are uniformly smooth in the sense of condition (i). Note that the constant here depends only on (H,p)(H,p) and indeed the resulting flex may of necessity be a colossal flex.

Theorem 2.20.

Let (G,p)(G,p) be an infinite locally finite framework in d{\mathbb{R}}^{d} with a connected graph, let

(G1,p)(G2,p),(G_{1},p)\subseteq(G_{2},p)\subseteq\dots,

be subframeworks, determined by finite subgraphs Gr=(Vr,Er)G_{r}=(V_{r},E_{r}) with union equal to GG and let v1,v2v_{1},v_{2} be vertices in G1G_{1}. Suppose moreover that for each r=1,2,,r=1,2,\dots, there is a base-fixed smooth flex pr(t)=(pkr(t))k=1|V(Gr)|p^{r}(t)=(p^{r}_{k}(t))_{k=1}^{|V{(G_{r})}|} of 𝒢r=(Gr,p){\mathcal{G}}_{r}=(G_{r},p) such that

(i) for each finite framework 𝒢l{\mathcal{G}}_{l} the set l{\mathcal{F}}_{l} of restriction flexes

{pr(t)|𝒢l:rl}\{p^{r}(t)|{\mathcal{G}}_{l}:r\geq l\}

have uniformly bounded derivatives, that is, there are constants Ml,l=1,2,,M_{l},l=1,2,\dots, such that

|ddtpkr(t)|Ml for rl,vkVl,|\frac{d}{dt}p_{k}^{r}(t)|\leq M_{l}\quad\mbox{ for }r\geq l,v_{k}\in V_{l},

(ii) the framework points p1,p2p_{1},p_{2} are uniformly separated by each flex pr(t)p^{r}(t) in the sense that

|p1r(1)p2r(1)||p1r(0)p2r(0)|c|p^{r}_{1}(1)-p^{r}_{2}(1)|-|p^{r}_{1}(0)-p^{r}_{2}(0)|\geq c

for some positive constant cc.

Then (G,p)(G,p) has a deformation.

Note that it is essential that the separated vertices of condition (ii) are the same for each subgraph. To see this note that the two-way infinite trapezium strip framework considered in Figure 7 has smooth deformations on each of its finite strip subframeworks, each of which ”separates” some two vertices (at the end of the strip) by a fixed positive distance. Nevertheless the infinite strip fails to have a deformation.

Proof.

For l=1,2,,l=1,2,\dots, let XlX_{l} be the space of continuous functions from [0,1][0,1] to d|Vl|{\mathbb{R}}^{d|V_{l}|} and note that the family l{\mathcal{F}}_{l}, by the hypotheses, is an equicontinuous family in XlX_{l}. Moreover with respect to the supremum norm l{\mathcal{F}}_{l} is a bounded set. By the Ascoli-Arzela theorem (see [30] for example) l{\mathcal{F}}_{l} is precompact and in particular for 1{\mathcal{F}}_{1} there is a subsequence r1,r2,r_{1},r_{2},\dots such that the restrictions prk|𝒢1p^{r_{k}}|{\mathcal{G}}_{1} converge uniformly in their (finitely many) coordinates to a flex q1q^{1} of 𝒢1{\mathcal{G}}_{1}. Relabel the sequence (prk)(p^{r_{k}}) as (p(k,1))(p^{(k,1)}). The restrictions of these flexes to 𝒢2{\mathcal{G}}_{2} similarly have a convergent subsequence, say (p(k,2))(p^{(k,2)}), and so on. From this construction select the diagonal subsequence (p(k,k))(p^{(k,k)}). This converges coordinatewise uniformly to a coordinatewise continuous function

q:[0,1]d×d×.q:[0,1]\to{\mathbb{R}}^{d}\times{\mathbb{R}}^{d}\times....

Since the restriction of qq to every finite subframework is a flex, qq satisfies the requirements of a flex of (G,p)(G,p), except possibly the properness requirement (v) of Definition 2.1. In view of (ii) however, q(0)q(1)q(0)\neq q(1) and so qq is a continuous flex of (G,p)(G,p). ∎

Remark 2.21.

Computer simulations provide evidence for the fact that small random perturbations of 𝒢2{\mathcal{G}}_{{\mathbb{Z}}^{2}} yield frameworks that are flexible. That is it seems that 𝒢2{\mathcal{G}}_{{\mathbb{Z}}^{2}} is stably flexible. It would be interesting if the theorem above could assist in a proof of this.

2.6. Flex extensions and generic rigidity

Let us recall a version of Laman’s theorem.

Theorem 2.22.

Let (G,p)(G,p) be an algebraically generic finite framework. Then (G,p)(G,p) is infinitesimally rigid if and only if the graph GG has a vertex induced subgraph HH, with V(H)=V(G)V(H)=V(G), which is maximally independent in the sense that 2|V(H)|=|E(H)|+32|V(H)|=|E(H)|+3 and 2|V(H)||E(H)|+32|V(H^{\prime})|\geq|E(H^{\prime})|+3 for every subgraph HH^{\prime} of HH.

For convenience we refer to a maximally independent finite graph as a Laman graph. We remark that any Laman graph can be obtained from a triangle graph by a sequence of moves known as Henneberg moves. The first of these adds a new vertex with two connecting edges while the second breaks an edge into two at a new vertex which is then connected by a new edge to another point of the graph.

Now let GG be an infinite graph which contains a subgraph HH on all the vertices of GG and suppose that HH is σ\sigma-Laman. In view of Laman’s theorem every algebraically generic realisation of HH (and hence GG) in the plane is σ\sigma-rigid and so continuously rigid. Is the converse true ? That is, if every generic realisation of an infinite graph GG is rigid does GG necessarily contain a σ\sigma-Laman subgraph with the same vertex set.

To see this one needs to show is that if GG is σ\sigma-(Laman-1) and not σ\sigma-Laman, then there exists a vertex generic realisation (G,p)(G,p) in the plane which has a continuous flex. That is, we want to build up a flex of the infinite structure by adding new vertices and edges, in the least handicapping way, to allow all, or most of the flex of an initial finite subgraph to be extended.

Alternatively, and more explicitly, suppose that one starts with a generic connected Laman-1 framework (G1,p)(G_{1},p) with nn vertices and an infinite sequence of Henneberg move ”instructions”. These instructions yield a unique infinite graph. Is it possible to choose associated framework points pn,pn+1,..p_{n},p_{n+1},.. so judiciously that some (perhaps small) flex of (G1,p)(G_{1},p) extends fully to each successive finite extension framework (and hence the infinite framework) ?

3. Rigidity Operators and Infinitesimal Rigidity

In previous sections we have considered some variety in the nature of continuous flexes p(t)p(t) and how they might distinguished. A companion consideration is the analysis of various spaces of infinitesimal flexes. This gives insight into continuous flexes since the derivative p(0)p^{\prime}(0) of a differentiable flex p(t)p(t) is an infinitesimal flex.

Here we give an operator theory perspective for an infinitesimal theory of infinite frameworks in which the rigidity matrix R(G,p)R(G,p) is viewed as a linear transformation or linear operator between various spaces. The domain space contains a space of infinitesimal flexes, which lie in the kernel of the rigidity operator, while the range space contains a space of self-stress vectors namely those in the kernel of the transpose of R(G,p)R(G,p).

3.1. Infinitesimal rigidity and the rigidity matrix.

Recall that for a finite framework (G,p)(G,p) in d{\mathbb{R}}^{d} with n=|V|n=|V| an infinitesimal flex is a vector u=(u1,,un)u=(u_{1},\dots,u_{n}) in the vector space v=dd{\mathcal{H}}_{v}={\mathbb{R}}^{d}\oplus\dots\oplus{\mathbb{R}}^{d} such that the orthogonality relation pipj,uiuj=0\langle p_{i}-p_{j},u_{i}-u_{j}\rangle=0 holds for each edge (vi,vj)(v_{i},v_{j}). This condition ensures that if each pip_{i} is perturbed to pi(t)=pi+tuip_{i}(t)=p_{i}+tu_{i}, with tt small, then the edge length perturbations are of second order only as tt tends to zero. That is, for all edges,

|pi(t)pj(t)||pipj|=O(t2).|p_{i}(t)-p_{j}(t)|-|p_{i}-p_{j}|=O(t^{2}).

If q(t):[1,1]vq(t):[-1,1]\to{\mathcal{H}}_{v} is a two-sided smooth flex of the finite framework (G,p)(G,p) then q(0)q^{\prime}(0) is an infinitesimal flex and for a generic finite framework every infinitesimal flex arises in this way. See Asimow and Roth [1] for example.

Associate with an infinite framework (G,p)(G,p) the product vector space

v=Vd=dd{\mathcal{H}}_{v}=\prod_{V}{\mathbb{R}}^{d}={\mathbb{R}}^{d}\oplus{\mathbb{R}}^{d}\oplus\dots

consisting of all sequences u=(u1,u2,)u=(u_{1},u_{2},\dots). Conceptually such a vector corresponds to a specification of instantaneous velocities, or to a perturbation sequence, applied to the framework joints. Define an infinitesimal flex of (G,p)(G,p) to be a vector uu in v{\mathcal{H}}_{v} for which, as above, pipj,uiuj=0\langle p_{i}-p_{j},u_{i}-u_{j}\rangle=0 holds for each edge (vi,vj)(v_{i},v_{j}), and let fl{\mathcal{H}}_{fl} denote the linear space of all these vectors. In the planar case fl{\mathcal{H}}_{fl} contains the three-dimensional linear subspace (assuming GG has at least one edge) of the infinitesimal flexes that arise from the isometries of 2{\mathbb{R}}^{2}. Note that a nonzero rotation infinitesimal flex uu is an unbounded sequence if and only if (G,p)(G,p) is an unbounded framework. We denote the space of rigid body motion infinitesimal flexes as rig{\mathcal{H}}_{rig}.

Definition 3.1.

An infinite framework (G,p)(G,p) is infinitesimally rigid if every infinitesimal flex is a rigid body motion infinitesimal flex.

If GG is infinite then v{\mathcal{H}}_{v} contains properly the direct sum space v00=V2{\mathcal{H}}_{v}^{00}=\sum_{V}\oplus{\mathbb{R}}^{2} consisting of vectors whose coordinates are finitely supported, in the sense of being finitely nonzero. The following definition is convenient and evocative.

Definition 3.2.

An infinite framework (G,p)(G,p) is internally infinitesimally rigid if every finitely supported infinitesimal flex is the zero flex.

Refer to caption
Figure 14. Unit cell for an internally infinitesimally flexible periodic framework.

We now give the usual direct definition of the rigidity matrix R(G,p)R(G,p) of a framework (G,p)(G,p), allowing GG to be infinite. This matrix could also be introduced via the Jacobian of the equation system that defines V(G,p)V(G,p) since 2R(G,p)2R(G,p) is the Jacobian evaluated at pp.

Write pi=(xi,yi)p_{i}=(x_{i},y_{i}), ui=(uix,uiy),i=1,2,u_{i}=(u^{x}_{i},u^{y}_{i}),i=1,2,\dots, and denote the coordinate difference xixjx_{i}-x_{j} by xijx_{ij}. The rigidity matrix is an infinite matrix R(G,p)R(G,p) with rows indexed by edges e1,e2,e_{1},e_{2},\dots and columns labeled by vertices but with multiplicity two, namely v1x,v1y,v2x,v2y,v_{1}^{x},v_{1}^{y},v_{2}^{x},v_{2}^{y},\dots. Note that any matrix of this shape, with finitely many nonzero entries in each row, provides a linear transformation from v{\mathcal{H}}_{v} to e=E{\mathcal{H}}_{e}=\prod_{E}{\mathbb{R}}.

Definition 3.3.

The rigidity matrix of the infinite framework (G,p)(G,p), with p=(pi)=(xi,yi)p=(p_{i})=(x_{i},y_{i}), is the matrix R(G,p)R(G,p) with entries xij,xji,yij,yjix_{ij},x_{ji},y_{ij},y_{ji} occurring in the row with label e=(vi,vj)e=(v_{i},v_{j}) with the respective column labels vix,vjx,viy,vjyv^{x}_{i},v_{j}^{x},v^{y}_{i},v^{y}_{j}, and with zero entries elsewhere.

In particular a vector uu in v{\mathcal{H}}_{v} is an infinitesimal flex if and only if R(G,p)u=0R(G,p)u=0.

Definition 3.4.

(i) A stress (or, more properly, a self-stress) of a finite or infinite framework (G,p)(G,p) is a vector w=(we)w=(w_{e}) in e=E{\mathcal{H}}_{e}=\prod_{E}{\mathbb{R}} such that R(G,p)tw=0.R(G,p)^{t}w=0.

(ii) A finite or infinite framework (G,p)(G,p) is isostatic, or absolutely isostatic, if it is infinitesimally rigid and has no nonzero self-stresses.

Since it is understood here that GG is a locally finite graph a rigidity matrix hase finitely many entries in each column and so its transpose corresponds to a linear transformations from e{\mathcal{H}}_{e} to v{\mathcal{H}}_{v}.

In the finite case a self-stress represents a finite linear dependence between the rows of the rigidity matrix, which one might abbreviate, with language abuse, by saying that the corresponding edges of the framework are linearly dependent. A self-stress vector w=(we)eEw=(w_{e})_{e\in E} can be simply related to a vector b=(be)b=(b_{e}) conceived of as a sequence b=(be)b=(b_{e}) of bar tension forces with a resolution, or balance, at each node. Indeed, for such a force vector bb the vector ww for which we=|pipj|1bew_{e}=|p_{i}-p_{j}|^{-1}b_{e} (e=(vi,vj)e=(v_{i},v_{j})) is a stress vector. Thus there is a simple linear relationship between the space of internal stresses and the space of resolving bar tensions. We shall not consider here the more general stress vectors, important in engineering applications, that arise from an external loading vectors.

Let e00{\mathcal{H}}_{e}^{00} be the space of finitely supported vectors in e{\mathcal{H}}_{e}. We say that an infinite framework (G,p)(G,p) is finitely isostatic if it is internally infinitesimally rigid and if the finite support stress space str00:=stre00{\mathcal{H}}_{str}^{00}:={\mathcal{H}}_{str}\cap{\mathcal{H}}_{e}^{00} is equal to {0}.\{0\}. It is straightforward to see that the grid frameworks 𝒢d{\mathcal{G}}_{{\mathbb{Z}}^{d}}, in their ambient spaces, are finitely isostatic, as is the kagome framework.

Between the extremes of infinitesimal rigidity and internal rigidity there are other natural forms of rigidity such as those given in the following definition. Write \ell^{\infty}, 2\ell^{2} and c0c_{0} to indicate the usual Banach sequence spaces for countable coordinates, and write e,v,,v0{\mathcal{H}}_{e}^{\infty},{\mathcal{H}}_{v}^{\infty},\dots,{\mathcal{H}}_{v}^{0} for the corresponding subspaces of e{\mathcal{H}}_{e} and v{\mathcal{H}}_{v}.

Definition 3.5.

An infinite framework (G,p)(G,p) is

(i) square-summably infinitesimally rigid (or infinitesimally 2\ell^{2}-rigid) if

v2kerR(G,p)={0},{\mathcal{H}}_{v}^{2}\cap\ker R(G,p)=\{0\},

(ii) boundedly infinitesimally rigid (or infinitesimally \ell^{\infty}-rigid) if

vkerR(G,p)=vrig,{\mathcal{H}}_{v}^{\infty}\cap\ker R(G,p)={\mathcal{H}}_{v}^{\infty}\cap{\mathcal{H}}_{rig},

(iii) vanishingly infinitesimally rigid (or infinitesially c0c_{0}-rigid) if

v0kerR(G,p)={0},{\mathcal{H}}_{v}^{0}\cap\ker R(G,p)=\{0\},

(iv) square-summably isostatic (or 2\ell^{2}-isostatic) if it is infinitesimally 2\ell^{2}-rigid and

stre2={0},{\mathcal{H}}_{str}\cap{\mathcal{H}}_{e}^{2}=\{0\},

(v) boundedly isostatic if it is boundedly infinitesimally rigid and stre={0}{\mathcal{H}}_{str}\cap{\mathcal{H}}_{e}^{\infty}=\{0\},

(vi) vanishingly isostatic if it is vanishingly infinitesimally rigid and stre0={0}{\mathcal{H}}_{str}\cap{\mathcal{H}}_{e}^{0}=\{0\}.

There is companion terminology for flexes and stresses. Thus we refer to vectors in v2kerR(G,p){\mathcal{H}}_{v}^{2}\cap\ker R(G,p) as square summable infinitesimal flexes and so on.

Example 3.6.

Let us use the shorthand (,p)({\mathbb{N}},p) to denote a semi-infinite framework in 2{\mathbb{R}}^{2} whose abstract graph is a tree with a single branch, with edges (v1,v2),(v2,v3),(v_{1},v_{2}),(v_{2},v_{3}),\dots and where p=(pi),pi=(xi,yi),i=1,2,p=(p_{i}),p_{i}=(x_{i},y_{i}),i=1,2,\dots . Then, writing xijx_{ij} and yijy_{ij} for the differences xixjx_{i}-x_{j} and yi=yjy_{i}=y_{j}, as before, the rigidity matrix with respect to the natural ordered bases takes the form

R(,p)=[x12y12x21y21000x23y23x32y3200000].R({\mathbb{N}},p)=\begin{bmatrix}x_{12}&y_{12}&x_{21}&y_{21}&0&\dots&&\\ 0&0&x_{23}&y_{23}&x_{32}&y_{32}&0&\dots\\ 0&0&0&0&*&*&*&\dots\\ \vdots&&&&&&&\end{bmatrix}.

With respect to the coordinate decomposition v=xy{\mathcal{H}}_{v}={\mathcal{H}}_{x}\oplus{\mathcal{H}}_{y} we have

R(,p)=[RxRy]=[DxDy][TT]R({\mathbb{N}},p)=\begin{bmatrix}R_{x}&R_{y}\end{bmatrix}=\begin{bmatrix}D_{x}&D_{y}\end{bmatrix}\begin{bmatrix}T&\\ &T\end{bmatrix}

where Rx=DxT,Ry=DyTR_{x}=D_{x}T,R_{y}=D_{y}T, where DxD_{x} and DyD_{y} are the diagonal matrices

Dx=[x12000x23000x340]Dy=[y12000y23000y340]D_{x}=\begin{bmatrix}x_{12}&0&0&&\dots\\ 0&x_{23}&0&&\dots\\ 0&0&x_{34}&0&\\ \vdots&&&\ddots&\end{bmatrix}\quad D_{y}=\begin{bmatrix}y_{12}&0&0&&\dots\\ 0&y_{23}&0&&\dots\\ 0&0&y_{34}&0&\\ \vdots&&&\ddots&\end{bmatrix}

and

T=[1100110011].T=\begin{bmatrix}1&-1&0&&\dots\\ 0&1&-1&&\dots\\ 0&0&1&-1&\\ \vdots&&&\ddots&\end{bmatrix}.

If we now identify the domain and range spaces in the natural way for these coordinates then we have T=IUtT=I-U^{t} where UtU^{t} is the transpose of the forward unilateral shift operator on the linear space of real sequences.

The analogous framework (,p)({\mathbb{Z}},p) has a similar matrix structure in all respects except that in place of the Toeplitz matrix TT one has the corresponding two-way infinite Laurent matrix IW1I-W^{-1} where WW is the forward bilateral shift. In both cases, the two-dimensional subspace spanned by the translation flexes is evident, being spanned by the constant vectors in x{\mathcal{H}}_{x} and y{\mathcal{H}}_{y}. Evidently there are infinitely many finitely supported flexes and in fact it is possible to identify kerR(G,p)\ker R(G,p) as a direct product space.

One can use operator formalism to examine the space of stresses. In the case of the simple framework (,p)({\mathbb{Z}},p) note that Wt=W1W^{t}=W^{-1} and

ker(IW)=ker(IW1)=e\ker(I-W)=\ker(I-W^{-1})={\mathbb{R}}e

where ee is the vector with every entry equal to 11. Since

R(,p)t=[IWIW][DxDy]R({\mathbb{Z}},p)^{t}=\begin{bmatrix}I-W&\\ &I-W\end{bmatrix}\begin{bmatrix}D_{x}\\ D_{y}\end{bmatrix}

it follows that a vector ww is a stress vector if and only if DxweD_{x}w\in{\mathbb{R}}e and DyweD_{y}w\in{\mathbb{R}}e. Thus for some constants α,β\alpha,\beta we have xi,i+1wi=α,yi,i+1wi=βx_{i,i+1}w_{i}=\alpha,y_{i,i+1}w_{i}=\beta, and so for all ii

yiyi+1xixi+1=βα.\frac{y_{i}-y_{i+1}}{x_{i}-x_{i+1}}=\frac{\beta}{\alpha}.

This colinearity condition shows that the space of stresses is trivial unless the framework points pip_{i}, ii\in{\mathbb{Z}} are colinear in which case str{\mathcal{H}}_{str} is one dimensional. This includes the colinear cases in which pp is a bounded sequence and the framework lies in a finite line segment in 2{\mathbb{R}}^{2}.

Example 3.7.

With similar notational economy write (r,p)({\mathbb{Z}}^{r},p), (resp. (r,p)({\mathbb{N}}^{r},p)) for frameworks associated with the grid graph with vertex set labeled by rr-tuples of integers (resp. positive integers) n=(n1,,nr)n=(n_{1},\dots,n_{r}) where the edges correspond to vertex pairs (n,n±ej)(n,n\pm e_{j}), where e1,,ere_{1},\dots,e_{r} are the usual basis elements. The ambient space for the framework is either understood or revealed by the entries of the vector pp.

Again one can use operator formalism to analyse the space of stresses as a vector subspace of {\mathcal{H}}. In the special case of the regular grid framework 𝒢2{\mathcal{G}}_{{\mathbb{Z}}^{2}} in 2{\mathbb{R}}^{2} one can see that the vector subspace str{\mathcal{H}}_{str} is a direct product vector space (like {\mathcal{H}} itself) whose product basis is indexed by (two-way infinite) linear subframeworks parallel to the coordinate axes. This is also true for a general orthogonal grid framework such as the bounded grid framework determined by the framework points (±(1(1/2)i),±(1(1/2)j)).(\pm(1-(1/2)^{i}),\pm(1-(1/2)^{j})).

Note that we have defined an infinitesimal flex in a local way, being the verbatim counterpart of the usage for finite frameworks. In particular the notion takes no account of the possibility of (second order) amplification or vanishing flexibility.

Our examples above indicate the importance of shift operators and in the next section we see that the bilateral shift operators, in their Fourier transform realisation as multiplication operators, play a central role in the discussion of periodic frameworks.

Let us note here that approximate infinitesimal flexes are natural for infinite frameworks and the operator theoretic perspective.

Definition 3.8.

An approximate square-summable flex of an infinite framework (G,p)(G,p) is a sequence of finitely supported unit vectors u1,u2,u_{1},u_{2},\dots in v2{\mathcal{H}}_{v}^{2} (or 𝒦v2{\mathcal{K}}_{v}^{2}) such that R(G,p)un20\|R(G,p)u_{n}\|_{2}\to 0 as nn\to\infty.

Let (G,p)(G,p) be a distance regular framework. Then it is straightforward to show that the rigidity matrix determines a bounded Hilbert space operator RR. It is the metrical and geometric properties of the action of RR and its transpose that have relevance to rigidity theory rather than the spectral theory of RR. However we do have the almost vacuous statement that the existence of approximate (square-summable) flexes corresponds to the point 0 belonging to the approximate point spectrum of RR. For if 0 lies in the approximate point spectrum then (by definition) RvnRv_{n} is a null sequence for some sequence of unit vectors vnv_{n}, and approximation of these unit vectors by vectors with finite support yields, after normalisations, an approximate square-summable flex sequence (un)(u_{n}).

It is implicit in the matricial function association below that the rigid unit modes of translationally periodic frameworks are tied to the existence of approximate flexes.

4. Crystal Frameworks and Flexibility

In previous sections we have constructed frameworks to illustrate various definitions and properties. It is perhaps of wider interest to understand, on the other hand, how extant infinite frameworks, such as those suggested by crystals or repetitive structures, may be flexible. Accordingly we now define crystal frameworks and investigate various forms of flexibility and rigidity.

4.1. Periodic and crystal frameworks.

We have already observed some properties of the basic examples of the grid framework 𝒢2{\mathcal{G}}_{{\mathbb{Z}}^{2}}, the squares framework 𝒢sq{\mathcal{G}}_{sq} and the kagome framework, 𝒢kag{\mathcal{G}}_{kag}. In 3{\mathbb{R}}^{3} we also have analogues, such as the cube framework 𝒢cube{\mathcal{G}}_{cube}, the octagon framework 𝒢oct{\mathcal{G}}_{oct} and the kagome net framework 𝒢knet{\mathcal{G}}_{knet}, which consist, respectively, of vertex-joined cubes, octahedra and tetrahedra, with no shared edges or faces, each in a natural periodic arrangement. The frameworks 𝒢kag{\mathcal{G}}_{kag}, 𝒢cube,𝒢oct{\mathcal{G}}_{cube},{\mathcal{G}}_{oct} and 𝒢knet{\mathcal{G}}_{knet} are polytope body-pin frameworks but we consider the polytope rigid units as bar-joint subframeworks formed by adding some, or perhaps all, internal edges. As such these frameworks are examples of crystal frameworks in the sense of the formal definition below. We first comment on a wider notion of periodicity.

Definition 4.1.

An affinely periodic framework in d{\mathbb{R}}^{d} is a framework 𝒢=(G,p){\mathcal{G}}=(G,p) for which there exists a non-trivial discrete group of affine transformations Tg,g𝒟T_{g},g\in{\mathcal{D}}, where each TgT_{g} acts on framework points and framework edges.

For example, the two-way infinite dyadic cobweb framework of Figure 15 is affinely periodic for the dilation doubling map and the four-fold dihedral group D4D_{4}.

Refer to caption
Figure 15. An affinely periodic cobweb framework.

For another example we may take the infinite {\mathbb{Z}}-periodic framework in three dimensions for which Figure 16 forms a perspective view down a central axis, with framework vertices (±1,±1,m),m(\pm 1,\pm 1,m),m\in{\mathbb{Z}}. Here the affine group is an isometry group isomorphic, as a group, to ×C2×D4{\mathbb{Z}}\times C_{2}\times D_{4}.

To illustrate the following definition observe in Figure 16 a template of six edges and three vertices which generates the kagome framework by the translations associated with the parallelogram unit cell. Borrowing crystallographic terminology we refer to such a template as a motif for the framework and the chosen translation group. The following formal definition gives a convenient way of specifying abstract crystal frameworks.

Refer to caption
Figure 16. A motif and unit cell for the kagome framework.
Definition 4.2.

A crystal framework 𝒞=(G,p){\mathcal{C}}=(G,p) in d{\mathbb{R}}^{d} is a connected bar-joint framework for which there is a discrete group of translation isometries 𝒯={Tg:g𝒟}{\mathcal{T}}=\{T_{g}:g\in{\mathcal{D}}\}, a finite connected set FeF_{e} of frameworks edges and a finite set FvF_{v} of framework vertices (being a subset of the vertex set of FeF_{e}) such that

(i) the unions

g𝒟Tg(Fe),g𝒟Tg(Fv)\cup_{g\in{\mathcal{D}}}T_{g}(F_{e}),\quad\quad\cup_{g\in{\mathcal{D}}}T_{g}(F_{v})

coincide with the sets of framework edges and vertices, and

(ii) these unions are unions of disjoint sets.

We denote a crystal framework 𝒞{\mathcal{C}} by the triple (Fv,Fe,𝒯)(F_{v},F_{e},{\mathcal{T}}) or by the triple (Fv,Fe,d)(F_{v},F_{e},{\mathbb{Z}}^{d}) in the case of integer translations. An associated unit cell for 𝒞{\mathcal{C}} may be defined as a set which contains FvF_{v} and no other framework points and for which the translates under the translation isometries are disjoint and partition the ambient space. For example in the case of 𝒢2{\mathcal{G}}_{{\mathbb{Z}}^{2}} we may take the semiopen set [0,1)2[0,1)^{2} or the set [0,1)×[1/2,5/2)[0,1)\times[1/2,5/2) as unit cells. Such parallelepiped unit cells are useful for us for torus models for crystal frameworks. Voronoi cells (Brillouin zones) also play a unit cell role in applications but we shall not need such geometric detail here.

In many applied settings the appropriate framework models have ”short” edges, spanning no more than two adjacent unit cells. Here we allow general edges which may span a chain of adjacent cells.

Recall from elementary crystallography that, modulo orthogonal transformations, there are 14 different forms (or symmetry types) in which a countable set of isolated points can be arranged with translational symmetry throughout three-dimensional space. These arrangements are called the Bravais lattices and the translation group 𝒯{\mathcal{T}} above corresponds to such a lattice. Thus each point of the framework lies in the Bravais lattice generated by the orbit of its unique corresponding motif vertex under the translational group.

4.2. Deformability and flow flexibility.

Recall that a general (countable, locally finite, connected) framework 𝒢{\mathcal{G}} is rigid if there is no base-fixed continuous flex and is boundedly rigid, or boundedly nondeformable, if there is no bounded base-fixed continuous flex p(t)p(t). Recall, from Section 2, that bounded flexes are those for which there is an absolute constant MM such that for every vertex vv the time separation |pv(t)pv(0)||p_{v}(t)-p_{v}(0)| is bounded by MM for all tt and all vv.

We first describe a context for the standard ”alternation” flexes of 𝒢sq{\mathcal{G}}_{sq} and 𝒢kag{\mathcal{G}}_{kag} and certain periodic flexes of 𝒢kag{\mathcal{G}}_{kag} with reduced symmetry. For these nonbounded flexes translational periodicity is maintained but relative to an affine flow of the ambient space. By an affine flow we mean simply a path tAtt\to A_{t} of affine transformations of d{\mathbb{R}}^{d} which is pointwise continuous. The simplest such flow in two dimensions is a contracting flow such as At(x,y)=((1tc)x,(1tc))yA_{t}(x,y)=((1-tc)x,(1-tc))y, 0<c<10<c<1. The alternating flexes of 𝒢sq{\mathcal{G}}_{sq} and 𝒢kag{\mathcal{G}}_{kag} are associated with such a flow.

Definition 4.3.

Let d=2,3d=2,3, let tAtt\to A_{t} be a flow of d{\mathbb{R}}^{d} and let 𝒞=(G,p){\mathcal{C}}=(G,p) be a crystal framework for the translation group 𝒯={Tg:g𝒟}{\mathcal{T}}=\{T_{g}:g\in{\mathcal{D}}\}. A flow-periodic flex of 𝒞{\mathcal{C}}, relative to the flow and the translation group 𝒯{\mathcal{T}}, is a continuous flex p(t)p(t) such that for each tt the framework (G,p(t))(G,p(t)) is 𝒯t{\mathcal{T}}^{t}- periodic, where 𝒯t={AtTgAt1:g𝒟}{\mathcal{T}}^{t}=\{A_{t}T_{g}A_{t}^{-1}:g\in{\mathcal{D}}\}.

Let us say simply that 𝒞{\mathcal{C}} is affinely periodically deformable if there exists a non trivial flow-periodic flex.

A flow-periodic flex p(t)p(t) for 𝒞{\mathcal{C}} can be defined for a given flow if and only if for each tt one can continuously solve the distance constraint equations for the vertex positions of the motif, with the periodicity constraint, in the AtA_{t}-deformed unit cell. Equivalently, the motif and unit cell define a finite framework on a torus with the noninterior edges of the motif providing reentrant (”locally geodesic”) edges on the torus. For example consider the torus framework in Figure 17.

Refer to caption
Figure 17. A torus framework.

A horizontal affine contraction Ht:(x,y)((1t)x,y)H_{t}:(x,y)\to((1-t)x,y) leads to a continuous flex of the torus framework and hence to a colossal flex of the associate crystal framework. The same is true for the vertical affine contraction VtV_{t} and for the skew affine transformation

St:(x,y)(x+(sint)y,y+(1cost)x).S_{t}:(x,y)\to(x+(\sin t)y,y+(1-\cos t)x).

Note that this particular transformation preserves both the cyclic width and cyclic height of the torus.

For an illustration of the method we note the following affine deformation result for what might be termed periodic cell-generic grid frameworks in d{\mathbb{R}}^{d}. For more general results of this nature see also Borcea and Streinu [2].

Theorem 4.4.

Let 𝒞=(d,p){\mathcal{C}}=({\mathbb{Z}}^{d},p) be a grid framework in d{\mathbb{R}}^{d} which is (n1,,nd)(n_{1},\dots,n_{d})-periodic and which is a 1/31/3-perturbation of 𝒢d{\mathcal{G}}_{{\mathbb{Z}}^{d}} in the sense that

|p(k1,,kd)(k1,,kd)|<1/3|p(k_{1},\dots,k_{d})-(k_{1},\dots,k_{d})|<1/3

for 0kini1,1id0\leq k_{i}\leq n_{i}-1,1\leq i\leq d. Then 𝒞{\mathcal{C}} is affinely periodically deformable.

Proof.

We sketch the proof in case d=2d=2 with (n1,n2)=(n,m)(n_{1},n_{2})=(n,m). For convenience re-scale the framework so that the large cell [0,n1)×[0,m1)[0,n-1)\times[0,m-1) of 𝒞{\mathcal{C}} becomes to the usual unit cell. Let (Fv,Fe,2)(F_{v},F_{e},{\mathbb{Z}}^{2}) be the motif representation of (2,p)({\mathbb{Z}}^{2},p). Suppose first that n=m=2n=m=2 and the motif has four noninterior edges (corresponding to the four reentrant edges on the associated torus). Let FeF_{e}^{\prime} be FeF_{e} with one vertical and one horizontal reentrant edge removed and let 𝒞{\mathcal{C}}^{\prime} be the associated framework. (See Figure 17.) For tt taking positive values in some finite interval there are affine deformations of 𝒞{\mathcal{C}}^{\prime} associated with each of the flows Ht,Vt,StH_{t},V_{t},S_{t}. Moreover there is an affine deformation associated with any composition At=Vβ(t)Hα(t)StA_{t}=V_{\beta(t)}H_{\alpha(t)}S_{t}, where the functions α\alpha and β\beta are any continuous functions with α(0)=β(0)=0\alpha(0)=\beta(0)=0, where tt takes values in some finite interval. Now note that we may chose α(t)\alpha(t) so that the separation distance corresponding to the omitted horizontal edge is constant, and we may choose β(t)\beta(t) similarly for the omitted vertical edge. Note also the essential fact that Vβ(t)V_{\beta(t)} does not change the cyclic width. In this way an affine deformation for 𝒞{\mathcal{C}} is determined. (In the case of generic points, is determined uniquely on some finite interval).

The same principle operates for general (n,m)(n,m). The subframework 𝒞{\mathcal{C}}^{\prime}, with two deleted noninterior edges (one for each coordinate direction) once again has a flow periodic deformation associated with a composition At=Vβ(t)Hα(t)StA_{t}=V_{\beta(t)}H_{\alpha(t)}S_{t}, for tt in a sufficiently small interval. This follows from a simple induction argument. Once again the functions α\alpha and β\beta are chosen to provide a deformation of 𝒞{\mathcal{C}}^{\prime} in which the separations for the omitted edges is constant. In this way a flow periodic deformation of 𝒞{\mathcal{C}} itself is determined. ∎

The deformations obtained in the last proposition are colossal deformations and this seems to be a necessary condition if some form of periodicity is to be maintained. It would be interesting to determine when such frameworks possess bounded ”unstructured” deformations.

One may also consider affine deformations relative to a subgroup of 𝒟{\mathcal{D}} associated with a supercell. By a supercell we mean a finite union (or possibly an infinite union) of adjacent cells which tile the ambient space by translations from the subgroup. In this case the deformations maintain only a longer period of translational symmetry and for an infinite linear supercell one may even forgo translational symmetry in one direction.

To illustrate this let us note a class of interesting deformations of the kagome framework (G,p)(G,p) which have this form.

Start with an alternation flex p(t)p(t) which leaves fixed a particular vertex pp_{*} and leaves invariant a line of hexagon diameters. Note that this flex (which is not base-fixed) has bilateral symmetry with respect to the line (viewed as a mirror line) but that the inversion symmetry of (G,p)(G,p) about the fixed vertex is broken for t>0t>0. Perform identical surgeries on the frameworks (G,p(t))(G,p(t)) as follows:

(i) cut the frameworks along the fixed line,

(ii) effect a reflection of one of the resulting half-planes frameworks in the orthogonal line through pp_{*} and rejoin the half-planes to create a new flex of (G,p(0))(G,p(0)).

We might view this flex as one with a symmetry transition line. The hexagons in this line maintain inversion symmetry while elsewhere the hexagons maintain bilateral symmetry.

One may perform such surgery on several parallel surgery lines simultaneously. Performing surgery on countably many lines leads to the following theorem. This shows, roughly speaking, that the kagome framework sits in its configuration space as an infinitely singular point in the sense that it is the starting point for uncountably many distinct flexes. (Here we use the term distinct in the sense of Definition 2.7.)

Theorem 4.5.

There are uncountably many distinct flow periodic flexes of the kagome framework.

4.3. Crystal frameworks and periodic infinitesimal flexibility.

Let 𝒞=(G,p)=(Fv,Fe,2){\mathcal{C}}=(G,p)=(F_{v},F_{e},{\mathbb{Z}}^{2}) be a crystal framework in the plane for the integer translation isometry group. A natural form of infinitesimal flex for 𝒞{\mathcal{C}} is that of a 11-cell periodic flex in the following sense.

Definition 4.6.

A vector u=(uv)vVu=(u_{v})_{v\in V} in the real vector space v{\mathcal{H}}_{v}, or in the complex vector space 𝒦v{\mathcal{K}}_{v}, is a 11-cell-periodic infinitesimal flex for 𝒞{\mathcal{C}} (or simply a periodic flex if there is no ambiguity) if ukerR(G,p)u\in\ker R(G,p) and uκ+n=uκu_{\kappa+n}=u_{\kappa} for all κFv\kappa\in F_{v} and ndn\in{\mathbb{Z}}^{d}.

Such an infinitesimal flex uu is a bounded sequence of vectors determined by periodic extension of what we may call the motif flex vector umotif=(uκ)κFv.u_{motif}=(u_{\kappa})_{\kappa\in F_{v}}. There is a one-to-one correspondence between these periodic flexes and the finite vectors that are in the kernel of the motif rigidity matrix Rm(G,p)R_{m}(G,p) which we may define as the natural ”periodic completion” of the |Fe|×d|Fv||F_{e}|\times d|F_{v}| submatrix of R(G,p)R(G,p). This is the natural representing matrix for R(G,p)R(G,p) viewed as a linear transformation between the finite dimensional subspaces of 2{\mathbb{Z}}^{2}-periodic vectors. If there are no reflexive edges in the motif then

Similarly we define periodic infinitesimal stresses as those which correspond to periodic extensions of vectors in the cokernel of the motif rigidity matrix.

The crystal framework 𝒞{\mathcal{C}}, with given motif and discrete translation group indexed by 𝒟{\mathcal{D}}, is said to be 11-cell-periodically isostatic (or simply periodically isostatic) if the only 11-cell periodic flexes are translation flexes and if there are no nontrivial 11-cell periodic stresses.

Once again it can be helpful to consider a flat torus model for such frameworks. Note for example that the generic periodic framework defined by Figure 17 is periodically rigid and indeed periodically isostatic.

The following interesting periodic variant of Laman’s theorem has been obtained recently by E. Ross [28]. Let us say that the planar periodic framework 𝒞{\mathcal{C}} is topologically proper if 𝒞{\mathcal{C}} is connected and for every torus subframework motif Fv,FeF_{v}^{\prime},F_{e}^{\prime} with 2|Fv||Fe|=22|F_{v}^{\prime}|-|F_{e}^{\prime}|=2 there is an edge cycle from FeF_{e}^{\prime} which properly wraps around the torus in the sense that the associated homotopy class is nonzero.

Theorem 4.7.

Let 𝒞=(Fv,Fe,2){\mathcal{C}}=(F_{v},F_{e},{\mathbb{Z}}^{2}) be an infinite framework in the plane which is periodic for the integer translation group and is topologically proper. Then following are equivalent.

(i) 2|Fv||Fe|=22|F_{v}|-|F_{e}|=2 and for all edge induced submotifs Fv,FeF_{v}^{\prime},F_{e}^{\prime} we have 2|Fv||Fe|22|F_{v}^{\prime}|-|F_{e}^{\prime}|\geq 2.

(ii) 𝒞{\mathcal{C}} is periodically isostatic in the sense that the periodic vectors in kerR(G,p)\ker R(G,p) are spanned by the two translations (periodic rigidity) and the periodic vectors in cokerR(G,p)\operatorname{coker}R(G,p) are zero (periodic stress free).

As we have seen, a periodically infinitesimally rigid framework in the sense above may nevertheless correspond to a flat torus framework which can flex infinitesimally (and even deform) relative to contraction or expansion of the unit cell. Allowing such vertical and horizontal freedoms adds two further degrees of freedom to the constraint equations. It would be of interest to obtain a similar characterisation in this case as well as for the triply flexible affine torus.

5. The Matricial Symbol Function of a Crystal Framework

We now derive matrix function operators for general (abstract) crystal frameworks. Our approach is decidedly Hilbert space theoretic and allows for an extended conceptual framework for rigidity analysis. First we outline the standard identification of operators commuting with shift operators and multiplication operators in a Fourier transform space. We adopt complex scalars, replacing the v{\mathcal{H}}_{v} spaces of Section 3 by their complex counterparts, 𝒦v,𝒦e,𝒦v,𝒦v2,{\mathcal{K}}_{v},{\mathcal{K}}_{e},{\mathcal{K}}_{v}^{\infty},{\mathcal{K}}_{v}^{2}, and so forth.

5.1. Matrix function multiplication operators.

Let ϕ\phi be a continuous complex-valued function on the unit circle 𝕋{\mathbb{T}}. It defines a multiplication operator TT on the usual complex Hilbert space L2(𝕋)L^{2}({\mathbb{T}}). Its representing matrix with respect to the standard orthonormal basis {zn:n}\{z^{n}:n\in{\mathbb{Z}}\} is the ×{\mathbb{Z}}\times{\mathbb{Z}} indexed matrix with entries

Ti,j=ϕzj,zi=ϕzji,1=ϕ^(ij),T_{i,j}=\langle\phi z^{j},z^{i}\rangle=\langle\phi z^{j-i},1\rangle=\hat{\phi}(i-j),

the (ij)th(i-j)^{th} Fourier coefficient of ϕ\phi.

Similarly, let Φ(z)\Phi(z) be a continuous matrix-valued function on the two or three dimensional torus 𝕋d{\mathbb{T}}^{d} taking values in the space of n×mn\times m complex matrices. One can specify such a function Φ(z)=Φ(z1,,zd)\Phi(z)=\Phi(z_{1},\dots,z_{d}) in terms of a matrix of scalar functions Φ(z)=[ϕk,l(z)]k=1,l=1n,m\Phi(z)=[\phi_{k,l}(z)]_{k=1,l=1}^{n,m}. Given in this way Φ(z)\Phi(z) defines a multiplication operator from L2(𝕋)mL^{2}({\mathbb{T}})\otimes{\mathbb{C}}^{m} to L2(𝕋)nL^{2}({\mathbb{T}})\otimes{\mathbb{C}}^{n}. Indeed, choose a basis {ξl}\{\xi_{l}\} for m{\mathbb{C}}^{m} and an associated basis {ξlzp}\{\xi_{l}\otimes z^{p}\} for the domain space. Similarly let {ηkzp}\{\eta_{k}\otimes z^{p}\} be a basis for the codomain. Then the operator TT of multiplication by Φ(z)\Phi(z) may be defined by specifying, for each continuous function f(z)f(z) on 𝕋d{\mathbb{T}}^{d},

T(ξlf(z))=k=1mηkϕk,l(z)f(z),T(\xi_{l}\otimes f(z))=\sum_{k=1}^{m}\eta_{k}\otimes\phi_{k,l}(z)f(z),

and extending by linearity and continuity. Thus we can see that the representing matrix for TT with respect to these bases is determined by the Fourier coefficients of the matrix entries for Φ\Phi:

T(ξlzp),ηkzqmL2(𝕋d)=ϕk,l(z)zpq,1L2(𝕋d)=ϕ^k,l(qp).\langle T(\xi_{l}\otimes z^{p}),\eta_{k}\otimes z^{q}\rangle_{{\mathbb{C}}^{m}\otimes L^{2}({\mathbb{T}}^{d})}=\langle\phi_{k,l}(z)z^{p-q},1\rangle_{L^{2}({\mathbb{T}}^{d})}=\hat{\phi}_{k,l}(q-p).

Viewing an element of L2(𝕋d)mL^{2}({\mathbb{T}}^{d})\otimes{\mathbb{C}}^{m} as a function F(z)F(z) taking values in m{\mathbb{C}}^{m} (strictly speaking, taking values almost everywhere), and similarly for vectors in the codomain, one also considers TFTF as the function Φ(z)F(z)\Phi(z)F(z). This operator is usually denoted as MΦM_{\Phi}.

Suppose now, starting afresh, we have orthonormal bases

{ξl,p:1lm,pd}{ηk,p:1kn,pd}\{\xi_{l,p}:1\leq l\leq m,p\in{\mathbb{Z}}^{d}\}\quad\{\eta_{k,p}:1\leq k\leq n,p\in{\mathbb{Z}}^{d}\}

for Hilbert spaces such as 𝒦v2{\mathcal{K}}_{v}^{2} and 𝒦e2{\mathcal{K}}_{e}^{2} respectively. Suppose also that we are given an operator TT from 𝒦v2{\mathcal{K}}_{v}^{2} to 𝒦e2{\mathcal{K}}_{e}^{2} by means of its representing matrix with respect to these bases and suppose moreover that this matrix has the translational symmetry above in the sense that the matrix entries

Tξl,p,ηk,q𝒦e2\langle T\xi_{l,p},\eta_{k,q}\rangle_{{\mathcal{K}}_{e}^{2}}

are independent of pqp-q. Let v:nL2(𝕋d)𝒦v2{\mathcal{F}}_{v}:{\mathbb{C}}^{n}\otimes L^{2}({\mathbb{T}}^{d})\to{\mathcal{K}}^{2}_{v} be the Fourier transform and e{\mathcal{F}}_{e} the Fourier transform for 𝒦e2{\mathcal{K}}^{2}_{e}. These are simply the unitary operators determined by the natural bijection of basis elements, namely ξl,pξlzp\xi_{l,p}\to\xi_{l}\otimes z^{p} and ηk,pηkzp\eta_{k,p}\to\eta_{k}\otimes z^{p}, respectively. Thus, from the scalar matrix entries the functions ϕk,l\phi_{k,l} are determined and hence the matrix function Φ(z)\Phi(z). In this way one identifies e1Tv{\mathcal{F}}_{e}^{-1}T{\mathcal{F}}_{v} as the associated multiplication operator Φ{\mathcal{M}}_{\Phi}.

5.2. The rigidity matrix as a multiplication operator

Let us keep in view the special case of the 2{\mathbb{Z}}^{2}-periodic grid framework determined by an algebraically generic quadrilateral in the unit cell.

Let p1,,p4p^{1},\dots,p^{4} be four framework vertices in the unit cell [0,1)2[0,1)^{2} constituting the motif set FvF_{v} of 2{\mathbb{Z}}^{2}-periodic quadrilateral grid (2,p)({\mathbb{Z}}^{2},p). Let pijt=pt+(i,j)p^{t}_{ij}=p^{t}+(i,j), for (i,j)2(i,j)\in{\mathbb{Z}}^{2}, be the general framework points and let e1,,e8e^{1},\dots,e^{8} denote the eight edges which form the motif set FeF_{e} given by

e1=[p1,p2],e2=[p2,p3],e3=[p3,p4],e4=[p4,p1],e^{1}=[p^{1},p^{2}],e^{2}=[p^{2},p^{3}],e^{3}=[p^{3},p^{4}],e^{4}=[p^{4},p^{1}],
e5=[p1,02,p1],e6=[p0,13,p2],e7=[p1,04,p3],e8=[p0,11,p4],e^{5}=[p^{2}_{-1,0},p^{1}],e^{6}=[p^{3}_{0,-1},p^{2}],e^{7}=[p^{4}_{1,0},p^{3}],e^{8}=[p^{1}_{0,1},p^{4}],

and let ei,jt=et+(i,j)e^{t}_{i,j}=e^{t}+(i,j) be a typical edge. Write the corresponding basis for the vector space 𝒦e0{\mathcal{K}}^{0}_{e} as

{ηt,i,j:1t8,(i,j)2}.\{\eta_{t,i,j}:1\leq t\leq 8,(i,j)\in{\mathbb{Z}}^{2}\}.

Also, label the natural basis for 𝒦v0{\mathcal{K}}^{0}_{v} as the set of vectors {ξs,i,jx,ξs,i,jy}\{\xi^{x}_{s,i,j},\xi^{y}_{s,i,j}\}, where ss ranges from 11 to 44 and (i,j)(i,j) range through 2{\mathbb{Z}}^{2}.

The rigidity matrix R(G,p)R(G,p) satisfies the following symmetry equations,

Wi,jR(G,p)Ui,j=R(G,p)W_{i,j}R(G,p)U_{i,j}=R(G,p)

with respect to the translation shift operators Ui,jU_{i,j} on 𝒦v00{\mathcal{K}}^{00}_{v} and Wi,jW_{i,j} on 𝒦v00{\mathcal{K}}^{00}_{v}. By the discussion above the rigidity matrix for this framework gives rise to an 8×88\times 8 matrix function Φ(z1,z2)=(ϕi,j(z1,z2)\Phi(z_{1},z_{2})=(\phi_{i,j}(z_{1},z_{2}) on the torus 𝕋2{\mathbb{T}}^{2} in 2{\mathbb{C}}^{2} which is determined by the equations

ϕ^k,{x,s}((i,j))=(R(G,p)ξs,i,jx,ηk,0,0,\hat{\phi}_{k,\{x,s\}}(-(i,j))=\langle(R(G,p)\xi^{x}_{s,i,j},\eta_{k,0,0}\rangle,

together with a companion set of equations for the yy-labeled basis elements. Furthermore there is a simple algorithm for computing the symbol matrix function from the motif (Fv,Fe)(F_{v},F_{e}) and the motif rigidity matrix.

For a general crystal frameworks in the plane we have the following recipe for identifying Φ(z1,z2)\Phi(z_{1},z_{2}). There is an entirely similar identification of the matrix function of a crystal framework in higher dimensions.

Theorem 5.1.

Let Φ(z)\Phi(z) be the matricial symbol function of the crystal framework 𝒞=(Fv,Fe,2){\mathcal{C}}=(F_{v},F_{e},{\mathbb{Z}}^{2}). The entry of the rigidity matrix determined by the edge eke_{k} of FeF_{e} and the column labeled vs,i,jxv^{x}_{s,i,j} (resp. vs,i,jyv^{y}_{s,i,j}) provides the (i,j)(-i,-j)-th Fourier coefficient of ϕk,{x,s}{\phi}_{k,\{x,s\}} (resp. ϕk,{y,s}{\phi}_{k,\{y,s\}}).

Note that for the motif for the periodic quadrilateral grid each edge has vertices with different ss index. This is rather typical and in such cases it follows that ϕk,{x,l}{\phi}_{k,\{x,l\}} is either zero or has one nonzero Fourier coefficient. In the case of a ”reflexive” edge in FeF_{e}, of the form [pκ,0,pκ,δ][p_{\kappa,0},p_{\kappa,\delta}] there are entries in the columns for κ\kappa corresponding to (1z¯δ)ve(1-\overline{z}^{\delta})v_{e} where vev_{e} is the usual vector of coordinate differences appearing in the rigidity matrix for the edge ee. That is, ve=pκ,0pκ,δv_{e}=p_{\kappa,0}-p_{\kappa,\delta}.

Thus we may obtain the following identification where the direct construction of Φ\Phi from 𝒞{\mathcal{C}} is given above.

Theorem 5.2.

Let 𝒞=(Fv,Fe,d){\mathcal{C}}=(F_{v},F_{e},{\mathbb{Z}}^{d}) be a crystal framework with m=|Fe|m=|F_{e}|, n=|Fv|n=|F_{v}| and with rigidity operator RR from 𝒦v2{\mathcal{K}}_{v}^{2} to 𝒦e2{\mathcal{K}}_{e}^{2}.

(i) The motif (Fv,Fe)(F_{v},F_{e}) determines a matrix-valued trigonometric function Φ:𝕋dMm,dn()\Phi:{\mathbb{T}}^{d}\to M_{m,dn}({\mathbb{C}}) for which there is a unitary equivalence e1Rv=MΦ{\mathcal{F}}_{e}^{-1}R{\mathcal{F}}_{v}=M_{\Phi} where MΦM_{\Phi} is the multiplication operator

MΦ:ndL2(𝕋d)mL2(𝕋d).M_{\Phi}:{\mathbb{C}}^{nd}\otimes L^{2}{({\mathbb{T}}^{d})}\to{\mathbb{C}}^{m}\otimes L^{2}({\mathbb{T}}^{d}).

(ii) The framework 𝒞{\mathcal{C}} is square-summably rigid (resp. square-summably stress-free) if the column rank (resp. row rank) of Φ(z)\Phi(z) is maximal for almost every z𝕋dz\in{\mathbb{T}}^{d}.

The following corollary shows that the existence of a square summable flex is a rather strong condition for a crystal framework and in fact ”typical” frameworks in Maxwell counting equilibrium have no such flexes. What is of particular significance for such frameworks however is the presence of approximate flexes, in the sense noted in Section 3, as these are related to the local rank degeneracies of Φ(z)\Phi(z) and to periodic flexes.

Corollary 5.3.

The following are equivalent for a crystal framework with Maxwell counting equilibrium.

(i) 𝒞{\mathcal{C}} has an nonzero internal (finitely supported) infinitesimal flex.

(ii) 𝒞{\mathcal{C}} has an nonzero summable infinitesimal flex.

(iii) 𝒞{\mathcal{C}} has an nonzero square-summable infinitesimal flex.

Proof.

Note that (i) implies (ii) and (ii) implies (iii) so it remains to show (iii) implies (i).

Suppose that the matrix Φ(z)\Phi(z) has a square-summable vector f(z)f(z) in its kernel. Then det(Φ(z))=0\det(\Phi(z))=0 for almost all points in the support of f(z)f(z). Since detΦ(z)\det\Phi(z) is a mutivariable trigonometric polynomial it necessarily vanished identically, that is the polynomial detΦ(z)\det\Phi(z) is the zero polynomial.

Recall that for any square matrix XX over a ring we have
XX~=det(X)InX\tilde{X}=\det(X)I_{n} where X~\tilde{X} is given by the usual formula
X~i,j=(1)i+jXi,jdet(Xij)\tilde{X}_{i,j}=(-1)^{i+j}X_{i,j}\det(X^{ij}) involving the cofactors (Xij)(X^{ij}). In particular if X~\tilde{X} is not zero and det(X)=0\det(X)=0 then for a nonzero column vector f(z)f(z) of X~\tilde{X} we have Xf(z)=0Xf(z)=0. Applying this to X=ΦX=\Phi provides the desired polynomial (finitely supported) vector in this case.

If X~\tilde{X} happens to be the zero matrix we make use of the minimal polynomial lemma below. To apply the lemma let qq be the minimum polynomial of XX so that q(X)=0q(X)=0. By the lemma we have q(X)=Xq1(X)q(X)=Xq_{1}(X) and q1(X)0q_{1}(X)\neq 0 by minimality. Thus there is a nonzero vector f(z)f(z) in the range of q1(X)q_{1}(X) and this provides the desired local flex. ∎

Lemma 5.4.

Let q(λ)q(\lambda) in R[λ]R[\lambda] be the minimal polynomial of an nn by nn matrix XX whose entries lie in an integral domain RR and suppose that detX=0\det X=0. Then q(0)=0q(0)=0.

Proof.

Let p(λ)p(\lambda) be the characteristic polynomial det(λInX)\det(\lambda I_{n}-X). Recall from the Cayley-Hamilton theorem that p(X)=0p(X)=0. In the algebraic closure of the field of fractions of RR the linear factors of q(λ)q(\lambda) and p(λ)p(\lambda) agree. See Lange [22] for example. Thus

q(λ)=(λαi)si,p(λ)=(λαi)riq(\lambda)=\prod(\lambda-\alpha_{i})^{s_{i}},\quad p(\lambda)=\prod(\lambda-\alpha_{i})^{r_{i}}

with siris_{i}\leq r_{i} for each ii. In particular for suitably large nn the polynomial q(0)nq(0)^{n} is divisible by p(0)p(0). Since detX=0\det X=0 it follows that p(0)=0p(0)=0 and so q(0)n=0q(0)^{n}=0. Since q(0)q(0) is in RR it is equal to zero and the proof is complete. ∎

5.3. From motifs to matrix functions.

We now consider some examples. Adjusting notation, write z,wz,w for the coordinates variables of 𝕋2{\mathbb{T}}^{2} so that a trigonometric polynomial takes the form a finite sum

ϕ(z,w)=i,jϕ^(i,j)ziwj.\phi(z,w)=\sum_{i,j}\hat{\phi}(i,j)z^{i}w^{j}.

Examining the recipe of Theorem 5.1 above in the case of the periodic quadrilateral grid we see that the internal edges e1,,e4e^{1},\dots,e^{4} of the motif provide four rows for Φ(z,w)\Phi(z,w) whose entries are constant functions with constants corresponding to those of the rows for the rigidity matrix R(G,p)R(G,p). The other edges provide rows according to the following simplified monomial rule. As we have intimated above, this rule applies generally when all the ”external” edges of the motif have vertices with distinct ss index:

the entries of the ethe^{th} row of Φ(z,w)\Phi(z,w) that appear in the columns for the external (non-motif) vertex are the corresponding entries for the motif rigidity matrix multiplied by the monomial z¯iw¯j\overline{z}^{i}\overline{w}^{j}, where (i,j)(i,j) is the shift index for the cell occupied by the external vertex.

Thus we obtain the matrix identification in the next proposition. We write pi=(xi,yi)p_{i}=(x_{i},y_{i}), for i=1,,4i=1,\dots,4. Note, for example, that the entry x12+1=x1(x21)x_{12}+1=x_{1}-(x_{2}-1) denotes the constant function corresponding to the entry of R(G,p)R(G,p) for row e5e^{5} and column x1x_{1}, while (x12+1)z-(x_{12}+1){z} is the entry corresponding to the xx-coordinate of the external vertex of e5e^{5}.

Proposition 5.5.

The matricial symbol function Φ(z,w)\Phi(z,w) of the periodic quadrilateral grid framework (2,p)({\mathbb{Z}}^{2},p) for the four vertex motif is the matrix function on 𝕋2{\mathbb{T}}^{2} given by

[x12y12x12y12000000x23y23x23y23000000x34y34x34y34x41y410000x41y41x12+1y12(x12+1)zy12z000000x23y23+1x23w(y23+1)w000000x341y34(1x34)z¯y34z¯x41w¯(1y41)w¯0000x41y411]\left[\begin{array}[]{cccccccc}{x_{12}}&{y_{12}}&-{x_{12}}&-{y_{12}}&0&0&0&0\\ \vskip 6.0pt plus 2.0pt minus 2.0pt\cr 0&0&{x_{23}}&{y_{23}}&-{x_{23}}&-{y_{23}}&0&0\\ \vskip 6.0pt plus 2.0pt minus 2.0pt\cr 0&0&0&0&{x_{34}}&{y_{34}}&-{x_{34}}&-{y_{34}}\\ \vskip 6.0pt plus 2.0pt minus 2.0pt\cr-{x_{41}}&-{y_{41}}&0&0&0&0&{x_{41}}&{y_{41}}\\ \vskip 6.0pt plus 2.0pt minus 2.0pt\cr{x_{12}}+1&{y_{12}}&-{(x_{12}}+1){z}&-{y_{12}}{z}&0&0&0&0\\ \vskip 6.0pt plus 2.0pt minus 2.0pt\cr 0&0&{x_{23}}&{y_{23}}+1&-{x_{23}}{w}&-{(y_{23}}+1){w}&0&0\\ \vskip 6.0pt plus 2.0pt minus 2.0pt\cr 0&0&0&0&{x_{34}}-1&{y_{34}}&(1-{x_{34}})\overline{z}&-{y_{34}}\overline{z}\\ \vskip 6.0pt plus 2.0pt minus 2.0pt\cr-{x_{41}}\overline{w}&(1-{y_{41}})\overline{w}&0&0&0&0&{x_{41}}&{y_{41}}-1\end{array}\right]

Furthermore, if the four points in the unit cell have algebraically independent coordinates then (2,p)({\mathbb{Z}}^{2},p) is square-summably isostatic.

Proof.

The first four rows are linearly independent since the points p1,,p4p_{1},\dots,p_{4} have algebraically independent coordinates. It can be shown from elementary linear algebra for almost every z,wz,w the full matrix has rank 88. Thus Theorem 5.2 applies. ∎

As a special case we obtain the matrix symbol function for 𝒢2{\mathcal{G}}_{{\mathbb{Z}}^{2}} determined by the eight-edged motif, namely

Φ(z,w)=1/2[1010000000010100000010100100000110z0000000010w00000010z¯00w¯000001]\Phi(z,w)=1/2\left[\begin{array}[]{cccccccc}-1&0&1&0&0&0&0&0\\ \vskip 6.0pt plus 2.0pt minus 2.0pt\cr 0&0&0&-1&0&1&0&0\\ \vskip 6.0pt plus 2.0pt minus 2.0pt\cr 0&0&0&0&1&0&-1&0\\ \vskip 6.0pt plus 2.0pt minus 2.0pt\cr 0&-1&0&0&0&0&0&1\\ \vskip 6.0pt plus 2.0pt minus 2.0pt\cr 1&0&-{z}&0&0&0&0&0\\ \vskip 6.0pt plus 2.0pt minus 2.0pt\cr 0&0&0&1&0&-{w}&0&0\\ \vskip 6.0pt plus 2.0pt minus 2.0pt\cr 0&0&0&0&-1&0&\overline{z}&0\\ \vskip 6.0pt plus 2.0pt minus 2.0pt\cr 0&\overline{w}&0&0&0&0&0&-1\end{array}\right]

The determinant in this case is

1256zw(1z¯)2(1w¯)2.-{\frac{1}{256}}\,zw{{\left(1-\overline{z}\right)^{2}\left(1-\overline{w}\right)^{2}}}.

On the other hand for this symmetric quadrilateral framework 𝒢2{\mathcal{G}}_{{\mathbb{Z}}^{2}} the formalism above can be applied to a smaller unit cell with the simple two-edged motif of Figure 18.

Refer to caption
Figure 18. A motif and unit cell for 𝒢2{\mathcal{G}}_{{\mathbb{Z}}^{2}}.

Note that the index ss for the recipe takes the value 11 only and both motif edges have vertices with the same ss index. It follows that the rigidity operator of 𝒢2{\mathcal{G}}_{{\mathbb{Z}}^{2}} associated with the motif is unitarily equivalent to the multiplication operator on the Hilbert space 2L2(𝕋2){\mathbb{C}}^{2}\otimes L^{2}({\mathbb{T}}^{2}) determined by the matricial symbol function

Ψ(z,w)=[z¯100w¯1].\Psi(z,w)=\begin{bmatrix}\overline{z}-1&0\\ 0&\overline{w}-1\end{bmatrix}.

One can similarly verify the following.

Proposition 5.6.

(i) The matricial symbol function of the kagome framework determined by the six-edge motif is given by

Φkag(z,w)=14[202000001313130013202z0000013z¯w3z¯ww¯3w¯0013]\Phi_{kag}({z},{w})={\frac{1}{4}\left[\begin{array}[]{cccccc}-2&0&2&0&0&0\\ \vskip 6.0pt plus 2.0pt minus 2.0pt\cr 0&0&1&-\sqrt{3}&-1&\sqrt{3}\\ \vskip 6.0pt plus 2.0pt minus 2.0pt\cr-1&-\,\sqrt{3}&0&0&1&\sqrt{3}\\ \vskip 6.0pt plus 2.0pt minus 2.0pt\cr 2&0&-2\,{z}&0&0&0\\ \vskip 6.0pt plus 2.0pt minus 2.0pt\cr 0&0&-1&\sqrt{3}&\overline{z}{w}&-\sqrt{3}\overline{z}{w}\\ \vskip 6.0pt plus 2.0pt minus 2.0pt\cr\overline{w}&\sqrt{3}\overline{w}&0&0&-1&-\sqrt{3}\end{array}\right]}

(ii) The determinant of Φkag(z,w)\Phi_{kag}(z,w) is a multiple of

zw(z¯1)(w¯1)(z¯w¯)zw(\overline{z}-1)(\overline{w}-1)(\overline{z}-\overline{w})

(iii) The kagome framework is square-summably isostatic.

We remark that the function matrix association above has also arisen in engineering in the analysis of Hutchinson and Fleck [19] of the stresses and rigidity of the kagome repetitive truss framework. This is derived from the crystallographic perspective of wave periodic flexes and Bloch’s theorem.

5.4. Symmetry equations for infinite frameworks.

Let 𝒞{\mathcal{C}} be a crystal framework in d{\mathbb{R}}^{d} with complex Hilbert space rigidity operator RR. Also let 𝒦v2,𝒦e2{\mathcal{K}}_{v}^{2},{\mathcal{K}}^{2}_{e} be as before, let 𝒦fl2{\mathcal{K}}^{2}_{fl} be the space of square-summable infinitesimal flexes, let 𝒦str2{\mathcal{K}}^{2}_{str} be the space of square-summable infinitesimal stresses and let =𝒦v2𝒦fl2{\mathcal{M}}={\mathcal{K}}_{v}^{2}\ominus{\mathcal{K}}^{2}_{fl}, 𝒩=𝒦e2𝒦str2{\mathcal{N}}={\mathcal{K}}_{e}^{2}\ominus{\mathcal{K}}^{2}_{str} be their complementary spaces.

As well as commuting with the coordinate shift operators the rigidity operator satisfies commutation relations for every isometric symmetry of 𝒞{\mathcal{C}}. For example, let SS be an isometric (not necessarily linear) operator on d{\mathbb{R}}^{d} which effects a symmetry of 𝒞{\mathcal{C}}. There is an induced unitary operator SeS_{e} on the complex Hilbert space 2(E)\ell^{2}(E) and an analogous operator SvS_{v} on 2(V)\ell^{2}(V) which permute the standard basis elements. Write S~v\tilde{S}_{v} for the isometric (not necessarily linear) operator SvSS_{v}\otimes S on v=2(V)d{\mathcal{H}}_{v}=\ell^{2}(V)\otimes{\mathbb{C}}^{d} where SS is also viewed as an isometric operator on d{\mathbb{C}}^{d}. Then we have the fundamental symmetry equation SeR=RS~vS_{e}R=R\tilde{S}_{v}. (See [27].)

Moreover, if 𝒢{\mathcal{G}} is the symmetry group of 𝒞{\mathcal{C}} arising from isometric transformations of d{\mathbb{R}}^{d} and if ρe\rho_{e} and ρ~v\tilde{\rho}_{v} are the associated representations of 𝒢{\mathcal{G}} on 𝒦e2{\mathcal{K}}_{e}^{2} and 𝒦v2{\mathcal{K}}_{v}^{2} then we have the symmetry equations

ρe(g)R=Rρ~v(g) for g𝒢.\rho_{e}(g)R=R\tilde{\rho}_{v}(g)\mbox{ for }g\in{\mathcal{G}}.

We have shown in [27] how such symmetry equations may be used to obtain a simple proof of a unitary equivalence which implies the Fowler-Guest formula [9] together with various generalisations. In the present setting we have the following analogue which also leads to counting conditions for isostatic and rigid frameworks.

Theorem 5.7.

Let 𝒞{\mathcal{C}} be a crystal framework with isometry symmetry group 𝒢{\mathcal{G}} and let

ρ~v=ρρfl\tilde{\rho}_{v}=\rho_{\mathcal{M}}\oplus\rho_{fl}
ρe=ρ𝒩ρstr\rho_{e}=\rho_{\mathcal{N}}\oplus\rho_{str}

be the decompositions associated with the spaces of square-summable infinitesimal flexes and stresses. Then ρ\rho_{\mathcal{M}} and ρ𝒩\rho_{\mathcal{N}} are unitarily equivalent representations. In particular if 𝒞{\mathcal{C}} is square summably isostatic then ρ~v\tilde{\rho}_{v} and ρe\rho_{e} are unitarily equivalent.

Proof.

In the square-summably isostatic case RR is a bounded operator with trivial kernel and trivial cokernel. The partially isometric part UU of the polar decomposition R=U(RR)1/2R=U(R^{*}R)^{1/2} is therefore unitary and it is a standard verification that this unitary also intertwines the representations. In general the symmetry equations show that the space 𝒦fl2{\mathcal{K}}^{2}_{fl} is reducing for ρ~v\tilde{\rho}_{v} and 𝒦str2{\mathcal{K}}^{2}_{str} is reducing for ρe\rho_{e} and so the asserted direct sum decompositions do exist. Now the restriction of RR to {\mathcal{M}} maps to 𝒩{\mathcal{N}} with trivial kernel and cokernel and so as before the restriction representations are unitarily equivalent. ∎

5.5. From matrix function to wave modes

For a crystal framework 𝒞=(Fv,Fe,d){\mathcal{C}}=(F_{v},F_{e},{\mathbb{Z}}^{d}) with a given motif let Φ:𝕋dMm,nd\Phi:{\mathbb{T}}^{d}\to M_{m,nd} be the associated matricial symbol function.

Definition 5.8.

The mode multiplicity function of 𝒞{\mathcal{C}} associated with the given motif and translation group is the function μ:𝕋d+\mu:{\mathbb{T}}^{d}\to{\mathbb{Z}}_{+} given by μ(z)=dimkerΦ(z).\mu(z)=\dim\ker\Phi(z).

For d=2d=2 we also consider μ\mu as being parametrised by coordinates s,ts,t in [0,1)×[0,1)[0,1)\times[0,1) so that (z,w)𝕋2(z,w)\in{\mathbb{T}}^{2} corresponds to (e2πis,e2πit)(e^{2\pi is},e^{2\pi it}). From the determinant calculations above we obtain the following.

Proposition 5.9.

(i) For the grid framework 𝒢2{\mathcal{G}}_{{\mathbb{Z}}^{2}} and the 8-edged motif the mode multiplicity function has values μ(0,0)=2\mu(0,0)=2,

μ(s,0)=μ(0,t)=1,\mu(s,0)=\mu(0,t)=1,

if ss and tt are nonzero, and is zero otherwise.

(ii) For the kagome framework 𝒢kag{\mathcal{G}}_{kag} and the 6-edged motif the mode multiplicity function has values μ(0,0)=2,\mu(0,0)=2, at the origin while μ(s,s)=μ(s,0)=μ(0,t)=1\mu(s,s)=\mu(s,0)=\mu(0,t)=1 if s0s\neq 0 and tt\neq and is zero otherwise.

Consider now the wave flexes of 𝒞{\mathcal{C}} which we define as the infinitesimal flexes which are 11-cell-periodic modulo a phase factor.

Definition 5.10.

Let 𝒞=(Fv,Ev,d){\mathcal{C}}=(F_{v},E_{v},{\mathbb{Z}}^{d}) be a crystal framework in d{\mathbb{R}}^{d}. An (infinitesimal) wave flex of 𝒞{\mathcal{C}} is a complex infinitesimal flex u=(uv)vVu=(u_{v})_{v\in V} which is wave periodic (or more precisely 11-cell wave periodic) in the sense that for some vector qq in d{\mathbb{R}}^{d}

uκ+n=e2πiq,nuκu_{\kappa+n}=e^{2\pi i\langle q,n\rangle}u_{\kappa}

for each vertex κ\kappa in the motif set FvF_{v} and each dd-tuple n=(n1,,nd)n=(n_{1},\dots,n_{d}).

The values of the mode multiplicity function corresponds to the dimension of the spaces of wave flexes. To see this suppose that Φ(z)\Phi(z) is the matricial symbol function for 𝒞=(G,p)=(Fv,Fe,3){\mathcal{C}}=(G,p)=(F_{v},F_{e},{\mathbb{Z}}^{3}), that w𝕋3w\in{\mathbb{T}}^{3} and that detΦ(w)=0\det\Phi(w)=0. Then Φ(w)um=0\Phi(w)u_{m}=0 for some nonzero complex motif vector umu_{m}. For the Dirac delta function δw(z)\delta_{w}(z) on 𝕋3{\mathbb{T}}^{3} we have, informally, Φ(z)F(z)=0\Phi(z)F(z)=0 for all zz on the 33-torus where F(z)F(z) is the function δw(z)um\delta_{w}(z)u_{m}. Thus, taking Fourier transforms it follows that the wave periodic vector u=(F)u={\mathcal{F}}(F) in 𝒦v{\mathcal{K}}_{v} satisfies R(G,p)u=0R(G,p)u=0. This shows that the bounded (and non square-summable) vector

u=(uκ,n)κFv,n3=(F)(κ,n)=wnuκu=(u_{\kappa,n})_{\kappa\in F_{v},n\in{\mathbb{Z}}^{3}}={\mathcal{F}}(F)(\kappa,n)=w^{n}u_{\kappa}

is a wave flex. The Dirac delta argument can be rigourised in the usual manner and so we obtain the following theorem.

Theorem 5.11.

Let 𝒞=(Fv,Fv,d){\mathcal{C}}=(F_{v},F_{v},{\mathbb{Z}}^{d}) be a crystal framework in d{\mathbb{R}}^{d} with associated symbol function Φ(z)\Phi(z) and suppose that 𝒞{\mathcal{C}} has no internal (finitely supported) self-stresses and no internal (finitely supported) flexes. Then

(i) d|Fv|=|Fe|d|F_{v}|=|F_{e}|, the associated matricial symbol function Φ(z)\Phi(z) on the dd-torus is square with full rank almost everywhere and 𝒞{\mathcal{C}} is square-summably isostatic.

(ii) infinitesimal wave flexes (that are 1-cell periodic-modulo-phase) exist, with phase factor w𝕋dw\in{\mathbb{T}}^{d}, if and only if detΦ(w)=0~\det\Phi(w)=0. In this case the dimension of the corresponding space of wave flexes is dimkerΦ(w)\dim\ker\Phi(w).

It follows also from this that the zero set of the determinant of the matricial symbol function associated with the motif of an appropriate supercell coincides with the phases of supercell-periodic wave flexes.

Of particular computational and theoretical interest is what one might refer to as the wave flex acquisition when a crystal framework deforms under a colossal flex to a framework with higher symmetry. This phenomenon serves as a model for the appearance of so-called Rigid Unit Modes (RUMs) in higher symmetry phases of various material crystals.

The following theorem generalises an interesting result of Wegner [32] for tetrahedral crystals. It may be viewed as an expression of the simple principle that additional symmetry often entails additional flexibility. Our proof applies to arbitrary crystal frameworks and is quite direct, benefitting somewhat from the economy of operator theory formalism. We remark that there are also natural operator algebra perspectives that are relevant to symmetry considerations. The result shows in particular for d=3d=3 that the RUM set is typically a union of surfaces, being the zero set of a single real-valued polynomial, rather than the intersection of the zero sets of the real and imaginary part of a complex polynomial.

For a given crystal framework and motif write Ω\Omega for the subset of 𝕋d{\mathbb{T}}^{d} formed by the phases ω\omega of the 11-cell periodic wave flexes. We also refer to this informally as the RUM set or the RUMs of 𝒞{\mathcal{C}}.

Theorem 5.12.

Let 𝒞=(Fv,Fe,d){\mathcal{C}}=(F_{v},F_{e},{\mathbb{Z}}^{d}) be a crystal framework in d{\mathbb{R}}^{d} with d|Fv|=|Fe|d|F_{v}|=|F_{e}| and suppose that 𝒞{\mathcal{C}} possesses inversion symmetry. Then the set Ω\Omega has the form

Ω=𝕋dV(p)\Omega={\mathbb{T}}^{d}\cap V(p)

where V(p)V(p) is the zero set of a complex polynomial p(z1,,zd,z¯1,z¯d)p(z_{1},\dots,z_{d},\overline{z}_{1}\dots,\overline{z}_{d}) which is real-valued on 𝕋d{\mathbb{T}}^{d}.

Proof.

Let us denote the set of framework vertices and edges as

𝒱={κ+n:κFv,nd}{\mathcal{V}}=\{\kappa+n:\kappa\in F_{v},n\in{\mathbb{Z}}^{d}\}
={e+n:eFv,nd}{\mathcal{E}}=\{e+n:e\in F_{v},n\in{\mathbb{Z}}^{d}\}

Effecting a translation, if necessary, we may assume that the inversion is σ:xx\sigma:x\to-x, that the unit cell has inversion symmetry and σ(Fv)=Fv\sigma(F_{v})=F_{v}. It may or may not be possible to re-choose the noninternal edges of FeF_{e} so that that σ(Fe)=Fe\sigma(F_{e})=F_{e}. Suppose first that this is the case. We show that detΦ(z)\det\Phi(z) itself is either real-valued or purely imaginary, from which the stated form for Ω\Omega follows.

With the notation of Section 5.4 we have the symmetry equation

ρe(σ)R=Rρ~v(σ).\rho_{e}(\sigma)R=R\tilde{\rho}_{v}(\sigma).

Recall that ρe(σ)\rho_{e}(\sigma) and ρ~v(σ)\tilde{\rho}_{v}(\sigma) are the isometric operators induced by the (isometric) symmetry element σ\sigma. Taking Fourier transforms this equation takes the form

U(σ)MΦ(z)=MΦ(z)V(σ)U(\sigma)M_{\Phi(z)}=M_{\Phi(z)}V(\sigma)

where U(σ)U(\sigma) and V(σ)V(\sigma) are the unitary operators determined by their action on the distinguished orthonormal bases. In view of the assumption we have

U(σ)(ηkzl)=ησe(k)zl,k=1,,|Fe|U(\sigma)(\eta_{k}\otimes z^{l})=\eta_{\sigma_{e}(k)}\otimes z^{-l},\quad k=1,\dots,|F_{e}|

(for certain induced permutations σe,σv\sigma_{e},\sigma_{v}) and

V(σ)(ξkxizl)=ξσv(k)xizl,k=1,,|Fv|,i=1,,d.V(\sigma)(\xi_{k}^{x_{i}}\otimes z^{l})=\xi_{\sigma_{v}(k)}^{x_{i}}\otimes z^{-l},\quad k=1,\dots,|F_{v}|,i=1,\dots,d.

That is, these operators have the form

U(σ)=EσJ,V(σ)=VσJU(\sigma)=E_{\sigma}\otimes J,\quad V(\sigma)=V_{\sigma}\otimes J

where JJ is the inversion unitary operator on L2(𝕋d)L^{2}({\mathbb{T}}^{d}) given by (Jf)(z)=f(z¯)(Jf)(z)=f(\overline{z}), and where EσE_{\sigma} and VσV_{\sigma} are scalar permutation matrices induced by σ\sigma. Substituting these forms we see that

Φ(z)=(EσJ)1Φ(z)(VσJ)=Eσ1Φ(z¯)Vσ.\Phi(z)=(E_{\sigma}\otimes J)^{-1}\Phi(z)(V_{\sigma}\otimes J)=E_{\sigma}^{-1}\Phi(\overline{z})V_{\sigma}.

It follows that

detΦ(z)=detΦ(z¯)det(Eσ)1det(Vσ)=(1)τdet(Φ(z))¯\det\Phi(z)=\det\Phi(\overline{z})\det(E_{\sigma})^{-1}\det(V_{\sigma})=(-1)^{\tau}\overline{\det(\Phi(z))}

for some integer τ\tau, since the determinant polynomial has real coefficients and so the determinant is either real or purely imaginary.

In the general case because of edges lying on inversion axes (lines through the origin) σ\sigma cannot be made to act freely on the motif edges. However for each motif edge index kk there is a monomial shift factor zp(k)z^{p(k)} such that with DD the diagonal matrix function

D=diag(zp(1),,zp(d))D=diag(z^{p(1)},\dots,z^{p(d)})

we have

U(σ)=D(EσJ).U(\sigma)=D(E_{\sigma}\otimes J).

Now we see that detΦ(z)=(1)τzpdet(Φ(z))¯\det\Phi(z)=(-1)^{\tau}z^{p}\overline{\det(\Phi(z))} for some multi-index pp. Write Φ=F1+iF2\Phi=F_{1}+iF_{2} with the FiF_{i} real-valued and (1)τzp=eiγ(z)(-1)^{\tau}z^{p}=e^{i\gamma(z)} with γ=γ(z)\gamma=\gamma(z) real. Equating real and imaginary parts it follows that either cosγ=1\cos\gamma=1 in which case detΦ(z)=det(Φ(z))¯\det\Phi(z)=\overline{\det(\Phi(z))} and Φ\Phi is real-valued, or F1(z)=(sinγ)(1cosγ)1F2(z)F_{1}(z)=(\sin\gamma)(1-\cos\gamma)^{-1}F_{2}(z). In all cases, for zz on the torus, F1(z)=0F_{1}(z)=0 if and only if and F2(z)=0F_{2}(z)=0, as required. ∎

5.6. Honeycomb frameworks and the kagome net.

The explicit formulation of the matricial symbol function makes clear a simple additive principle which is useful for calculation : adding internal edges to a unit cell results in adding the same number of rows to the symbol function to form the new symbol function, and each of these rows has the standard form with constant entries.

Consider the following honeycomb frameworks based on regular
hexagons.

Refer to caption
Figure 19. A 33-regular honeycomb framework.
Refer to caption
Figure 20. A 44-regular honeycomb framework.

The 33-regular (symmetric) honeycomb framework has a 99 by 1212 matricial symbol function. There are vectors in the kernel of the multiplication operator providing internal infinitesimal flexes and there are evident bounded infinitesimal flexes associated with ”parallel deformations”, for example.

A 44-regular framework may be constructed by rigidifying the hexagon internal to the unit cell, by adding three cross edges. The matricial symbol function is correspondingly enlarged by three scalar rows, appearing in rows 99 to 1212 in the following matrix.

[1/301/3000000000001/61/631/61/6300000000001/61/631/61/6300000000001/301/3000000000001/61/631/61/631/61/63000000001/61/6300001/3000001/3z¯01/6w¯1/63w¯00001/61/630000001/6w¯1/63w¯00001/6z¯1/63z¯001/31/3300001/31/330000001/31/3300001/31/330000002/3000002/30]\left[\begin{array}[]{cccccccccccc}-1/3&0&1/3&0&0&0&0&0&0&0&0&0\\ \vskip 6.0pt plus 2.0pt minus 2.0pt\cr 0&0&-1/6&-1/6\,\sqrt{3}&1/6&1/6\,\sqrt{3}&0&0&0&0&0&0\\ \vskip 6.0pt plus 2.0pt minus 2.0pt\cr 0&0&0&0&1/6&-1/6\,\sqrt{3}&-1/6&1/6\,\sqrt{3}&0&0&0&0\\ \vskip 6.0pt plus 2.0pt minus 2.0pt\cr 0&0&0&0&0&0&1/3&0&-1/3&0&0&0\\ \vskip 6.0pt plus 2.0pt minus 2.0pt\cr 0&0&0&0&0&0&0&0&1/6&1/6\,\sqrt{3}&-1/6&-1/6\,\sqrt{3}\\ \vskip 6.0pt plus 2.0pt minus 2.0pt\cr 1/6&-1/6\,\sqrt{3}&0&0&0&0&0&0&0&0&-1/6&1/6\,\sqrt{3}\\ \vskip 6.0pt plus 2.0pt minus 2.0pt\cr 0&0&0&0&-1/3&0&0&0&0&0&1/3\,\overline{z}&0\\ \vskip 6.0pt plus 2.0pt minus 2.0pt\cr 1/6\,\overline{w}&1/6\,\sqrt{3}\overline{w}&0&0&0&0&-1/6&-1/6\,\sqrt{3}&0&0&0&0\\ \vskip 6.0pt plus 2.0pt minus 2.0pt\cr 0&0&-1/6\,\overline{w}&1/6\,\sqrt{3}\overline{w}&0&0&0&0&1/6\,\overline{z}&-1/6\,\sqrt{3}\overline{z}&0&0\\ \vskip 6.0pt plus 2.0pt minus 2.0pt\cr-1/3&-1/3\,\sqrt{3}&0&0&0&0&1/3&1/3\,\sqrt{3}&0&0&0&0\\ \vskip 6.0pt plus 2.0pt minus 2.0pt\cr 0&0&1/3&-1/3\,\sqrt{3}&0&0&0&0&-1/3&1/3\,\sqrt{3}&0&0\\ \vskip 6.0pt plus 2.0pt minus 2.0pt\cr 0&0&0&0&2/3&0&0&0&0&0&-2/3&0\end{array}\right]

The cross-barred hexagon is continuously rigid but is infinitesimally flexible. We refer to this flex as the in-out flex. (See Whiteley [34] for related discussions.) It is straightforward to show that this flex extends periodically as a 1-cell-periodic infinitesimal flex of the framework.

The motif matrix function is obtained by evaluating at z=w=1z=w=1. This has rank 88 reflecting the four independent 11-cell-periodic flexes corresponding to horizontal translation, to vertical translation, to local rotation of the cross-barred hexagon, and to the periodic ”in-out” flex.

We define the kagome net framework as the three dimensional crystal framework 𝒢knet=(Fv,Fe,3){\mathcal{G}}_{knet}=(F_{v},F_{e},{\mathbb{Z}}^{3}) determined by the motif diagram in Figure 21. This figure shows a tetrahedron internal to a parallelepiped unit cell together with six additional non-internal edges, each of which extends a tetrahedron edge as indicated.

Refer to caption
Figure 21. A twelve-edged motif for the kagome net, 𝒢knet{\mathcal{G}}_{knet}.

For the given motif the motif rigidity matrix Rm(𝒢knet)R_{m}({\mathcal{G}}_{knet}) takes the form of a 1212 by 1212 matrix valued function on the 33-torus. Different motifs, for the same translation group, give matrix functions that are related by permutations on columns and monomial multiplication on rows. It can be shown that the determinant of such a matrix function is equal to a scalar and monomial multiple of the polynomial

(z1)(w1)(u1)(zw)(wu)(uz).\left(z-1\right)\left(w-1\right)\left(u-1\right)\left(z-w\right)\left(w-u\right)\left(u-z\right).

This is in complete analogy with the kagome framework and in the light of earlier discussions may draw the following conclusions.

Theorem 5.13.

(i) The kagome net 𝒢knet{\mathcal{G}}_{knet} is square-summably isostatic and possesses no internal infinitesimal flexes.

(ii) The mode multiplicity function μ:𝕋3+\mu:{\mathbb{T}}^{3}\to{\mathbb{Z}}_{+} of the kagome net framework has support equal to the intersection of 𝕋3{\mathbb{T}}^{3} with the six planes

z=1,w=1,u=1,z=w,w=u,u=z.z=1,w=1,u=1,z=w,w=u,u=z.

5.7. Crystallography and Rigid Unit Modes.

Rigidity theory and a Hilbert space operator viewpoint have led us to determine, ab initio, the matricial symbol function Φ(z)\Phi(z) of an abstract crystal framework and motif. The zeros of detΦ(z)\det\Phi(z) correspond to the phases of infinitesimal periodic-modulo-phase wave flexes and the mode multiplicity function μ(w)\mu(w) of Φ\Phi detects the multiplicities of independent wave flexes for ww. There are close connections between these mathematical observations and the Rigid Unit Modes (RUMs) that are observed in certain material crystals through diffraction experiments. This connection is part of the motivation for the formulation of the function μ(z)\mu(z). It seems to us that the explicit algorithm for the passage from crystal motif to matrix function will provide a useful computational and theoretical tool for identifying RUMs and their relationships.

Infinitesimal wave flexes for abstract crystal frameworks appear in the classical crystallography of Born von Karman theory. For example in the case of variants of quartz the low energy excitation modes result from rigid unit motions of SiO4SiO_{4} molecular tetrahedra within a tetrahedral net. These modes are therefore modelled by the first order periodic-modulo-phase infinitesimal rigid motions of abstract tetrahedral net frameworks. For this connection see, for example, Chapter 3 of [35]. Comparisons of extensive simulations and experimental results have shown that the RUM modes of simulated crystals are closely correlated with observed RUMs. For this see the seminal paper of Giddy et al [11], and also, for example, Dove et al [8] and Goodwin et al [14].

An interesting advance has been obtained recently by Wegner [32] who has derived mathematically, rather than through simulation, the RUM sets for various idealised tetrahedral crystals that model the geometry of material crystals, including β\beta-cristobalite, HP tridymite, β\beta-quartz, α\alpha-cristobalite and α\alpha-quartz. These models are standard framework models that coincide with crystal bar-joint frameworks in our terminology and the results are obtained by determining vanishing determinants (zero sets of detΦ\det\Phi in our formalism) by means of computer algebra and symmetry reductions. In particular Wegner obtains an analytic derivation of some of the curious surfaces observed in experiments of Dove et al.

References

  • [1] L. Asimow and B. Roth, The rigidity of graphs, Trans. Amer. Math. Soc., 245 (1978), 279-289.
  • [2] C.S. Borcea and I. Streinu, Periodic frameworks and flexibility, Proc. R. Soc. A, doi:10.1098/rspa.2009.0676.
  • [3] A. Böttcher, B. Silbermann, Analysis of Toeplitz operators, Springer-Verlag, Berlin, 1990.
  • [4] R.Connelly, P.W.Fowler, S.D.Guest, B.Schulze, W.J.Whiteley, ”When is a pin-jointed framework isostatic?”, International Journal of Solids and Structures, 46 (2009), 762 - 773.
  • [5] D. Cox, J. Little, D. O’Shea, Ideals, Varieties and Algorithms, Springer-Verlag, 1992.
  • [6] M. V. Chubynsky and M. F. Thorpe, Algorithms for three-dimensional rigidity analysis and a first-order percolation transition, Physical Review E 76, 2007.
  • [7] A. Donev and S. Torquato, Energy-efficient actuation in infinite lattice structures, J. Mech Phys. Solids, 51 (2003) 1459-1475.
  • [8] M. T. Dove et al, Floppy Modes in Crystalline and Amorphous Silicates; in Rigidity Theory and Applications, Edited by M. F. Thorpe and P. M. Duxbury, Kluwer Academic/Plenum Publishers, New York (1999).
  • [9] P.W. Fowler and S.D. Guest, A symmetry extension of Maxwell’s rule for rigidity of frames, International Journal of Solids and Structures 37 (2000) 1793 - 1804.
  • [10] M. G. Frost and C. Storey, Equivalence of a matrix over [s,z]{\mathcal{R}}[s,z] with its Smith form, Int. J. Control, 28 (1978), 665-671.
  • [11] A. P. Giddy, M.T. Dove, G.S. Pawley, V. Heine, The determination of rigid unit modes as potential soft modes for displacive phase transitions in framework crystal structures. Acta Crystallogr., A49 (1993), 697 - 703.
  • [12] H. Gluck, Almost all simply connected closed surfaces are rigid, in Geometric Topology, Lecture Notes in Math., no. 438, Springer-Verlag, Berlin, 1975, pp. 225-239.
  • [13] I. Gohberg, P. Lancaster, L. Rodman, Matrix Polynomials, Classics in Applied Mathematics, SIAM, 2009.
  • [14] A. L. Goodwin, M. Calleja, M. J. Conterio, M. T. Dove, J. S. O. Evans, D. A. Keen, L. Peters, M. G. Tucker, Colossal Positive and Negative Thermal Expansion in the Framework Material Ag3[Co(CN)6], Science 8 February 2008: Vol. 319. no. 5864, pp. 794 - 797.
  • [15] J. Graver, B. Servatius and H. Servatius, Combinatorial rigidity, Graduate Texts in Mathematics, vol 2, Amer. Math. Soc., 1993.
  • [16] S. D. Guest and J. W. Hutchinson, On the determinacy of repetitive structures, Journal of the Mechanics and Physics of Solids, 51 (2003), 383 - 391.
  • [17] B. Grunbaum, G.C. Shephard, Tilings and patterns, W.H. Freeman, New York, 1987.
  • [18] K.D. Hammonds, M.T. Dove, A.P. Giddy, V. Heine, B. Winkler, Rigid-unit phonon modes and structural phase transitions in framework silicates, American Mineralogist, 81 (1996), 1057-1079.
  • [19] R.G. Hutchinson and N.A. Fleck, The structural performance of the periodic truss, Journal of the Mechanics and Physics of Solids, 54 (2006) 756-782.
  • [20] A. B. Kempe, On a general method of describing plane curves of the nth degree by linkwork, Proc. London Math. Soc., 7 (1876), 213-216.
  • [21] G. Laman, On graphs and the rigidity of plane skeletal structures, J. Engineering Mathematics, 4 (1970), 331-340.
  • [22] S. Lang, Algebra, Addison Wesley, 1993.
  • [23] J.C. Maxwell, On the calculation of the equilibrium and stiffness of frames. Philosophical Magazine 27, (1864), 294 - 299.
  • [24] J.C. Owen and S.C. Power, The non-solvability by radicals of generic 3-connected planar Laman graphs, Trans. Amer. Math. Soc., 359 (2007), 2269-2303.
  • [25] J.C. Owen and S.C. Power, Infinite bar-joint frameworks, Proceedings of the Symposium in Applied Computing, (SAC 2009) March 8-12, 2009, Honolulu, Hawaii, U.S.A.
  • [26] J.C. Owen and S.C. Power, Continuous curves from infinite Kempe linkages, Bull. London Math. Soc., 2009; doi: 10.1112/blms/bdp087
  • [27] J.C. Owen and S.C. Power, Frameworks, symmetry and rigidity, to appear in the Inter. J. of Computational Geometry and Applications, arXiv:0812.3785.
  • [28] E. Ross, private communication, July 2009.
  • [29] B. Roth, Rigid and flexible frameworks, American Math. Monthly, Vol. 88 (1981), 6-21.
  • [30] W. Rudin, Functional Analysis, McGraw and Hill, NY, 1991.
  • [31] M. F. Thorpe and P. M. Duxbury, Rigidity Theory and Applications, Edited by Kluwer Academic/Plenum Publishers, New York (1999).
  • [32] F. Wegner, Rigid-unit modes in tetrahedral crystals, Journal of Physics: Condensed Matter, Volume 19, Issue 40, (2007).
  • [33] S.A. Wells, M. T. Dove, M. G. Tucker and K. Trachenko, Real-space rigid-unit-mode analysis of dynamic disorder in quartz, cristobalite and amorphous silica, J. Phys.: Condens. Matter 14 (2002), 4645 4657.
  • [34] W. Whiteley, Matroids and rigid structures, in ”Matroid Applications” ed. N. White, Encyclodedia of Mathematics and its applications 40 (1992), 1-51.
  • [35] B.T.M.Willis and A.W. Pryor, Thermal vibrations in crystallography, Cambridge University Press, 1975.
  • [36] M. Wyart, Rigidity-based approach to the boson peak in amorphous solids: from sphere packing to amorphous silica, preprint 2008, arXiv:0806.4596v1.