networks of balls and springs or frames of nodes connected by compressible struts provide realistic models for physical systems from bridges to condensed solids .
their elastic properties depend on their coordination number @xmath0  the average number of nodes each node is connected to . if @xmath0 is large enough , the networks are elastic solids whose long - wavelength mechanical properties are described by a continuum elastic energy with non - vanishing elastic moduli . if @xmath0 is small enough , the networks have deformation modes of zero energy
they are floppy .
as @xmath0 is increased from the floppy side , a critical value , @xmath4 , is reached at which springs provide just enough constraints that the system has no zero - energy  floppy " modes @xcite ( or mechanisms @xcite in the engineering literature ) , and the system is _ isostatic_. the phenomenon of rigidity percolation @xcite
whereby a sample spanning rigid cluster develops upon the addition of springs is one version of this floppy - to - rigid transition .
the coordination numbers of whole classes of systems , including engineering structures @xcite ( bridges and buildings ) , randomly packed spheres near jamming @xcite , network glasses @xcite , cristobalites @xcite , zeolites @xcite , and biopolymer networks @xcite are close enough to @xmath4 that their elasticity and mode structure is strongly influenced by those of the isostatic lattice .
( a ) section of a kagome lattice with @xmath5 and @xmath6 three - site unit cells .
nearest - neighbor bonds , occupied by harmonic springs , are of length @xmath7 . the rotated row ( second row from the top )
represents a floppy mode .
next - nearest neighbor bonds are shown as dotted lines in the lower left hexagon .
the vectors @xmath8 , @xmath9 , and @xmath10 indicate symmetry directions of the lattice .
the numbers in the triangles indicate those that twist together under pbcs in zero modes along the three symmetry direction .
note that there are only @xmath11 of these modes .
( b ) section of a square lattice depicting a floppy mode in which all sites along a line are displaced uniformly .
( c ) twisted kagome lattice , with lattice constant @xmath12 , derived from the undistorted lattice by rigidly rotating triangles through an angle @xmath13 . a unit cell , bounded by dashed lines , is shown in violet .
arrows depict site displacements for the zone - center , i.e. , zero wavenumber , @xmath14 mode which has zero ( nonzero ) frequency under free ( periodic ) boundary conditions .
sites @xmath15 , @xmath16 , and @xmath17 undergo no collective rotation about their center of mass whereas sites @xmath15 , @xmath18 , and @xmath19 do .
( d ) superposed snapshots of the twisted lattice showing decreasing areas with increasing @xmath13 . ]
though the isostatic point always separates rigid from floppy behavior , the properties of isostatic lattices are not universal ; rather they depend on lattice architecture . here
we explore the the unusual properties of a particular class of periodic isostatic lattices derived from the two - dimensional kagome lattice by rigidly rotating triangles through an angle @xmath13 without changing bond lengths as shown in fig .
[ fig : periodic - isostatic ] .
the bulk modulus @xmath20 of these lattices is rigorously zero for all @xmath21 . as a result
, their poisson ratio acquires its limit value of @xmath22 ; when stretched in one direction , they _ expand _ by an equal amount in the orthogonal direction : they are maximally auxetic @xcite .
these modes represent collapse pathways @xcite of the kagome lattice .
modes of isostatic systems are generally very sensitive to boundary conditions @xcite , but the degree of sensitivity depends on the details of lattice structure . for reasons we will discuss more fully below , modes of the square lattice , which is isostatic , are in fact insensitive to changes from free boundary conditions ( fbcs ) to periodic boundary conditions ( pbcs ) , whereas those of the undistorted kagome lattice are only mildly so . the modes of both ,
however , change significantly when rigid boundary conditions ( rgbs ) are applied .
we show here that in all families of the twisted kagome lattice , modes depend sensitively on whether fbcs , pcbs or rgbs are applied : finite lattices with free boundaries have floppy surface modes that are not present in their periodic or rigid spectrum or in that of finite undistorted kagome lattices . in the long wavelength limit
, the surface floppy modes , which are present in any @xmath23 material with @xmath24 , reduce to surface rayleigh waves @xcite described by a conformally invariant energy whose analytic eigenfunctions are fully determined by boundary conditions . at shorter wavelengths ,
the surface waves become sensitive to lattice structure and remain confined to within a distance of the surface that diverges as the undistorted kagome lattice is approached . in the simplest twisted kagome lattice ,
all floppy modes are surface modes , but in more complicated lattices , including ones with uniaxial symmetry , we construct , there are both surface and bulk floppy modes .    arguments due to j.c .
maxwell @xcite provide a criterion for network stability : networks in @xmath1 dimensions consisting of @xmath25 nodes , each connected with central - force springs to an average of @xmath0 neighbors , have @xmath26 zero - energy modes when @xmath3 ( in the absence of redundant bonds - see below ) . of these a number , @xmath27 , which depends on boundary conditions , are trivial rigid translations and rotations , and the and the remainder are floppy modes of internal structural rearrangement . under fbcs
an pbcs , @xmath28 equals @xmath29 and @xmath1 , respectively . with increasing @xmath0 , mechanical stability
is reached at the isostatic point at which @xmath30 .
the maxwell argument is a global one ; it does not provide information about the nature of the floppy modes and does not distinguish between bulk or surface modes .
the kagome lattice of central force springs shown in fig .
[ fig : periodic - isostatic](a ) is one of many locally isostatic lattices , including the familiar square lattice lattice in two dimensions [ fig .
[ fig : periodic - isostatic](b ) ] and the cubic and pyrochlore lattices in three dimensions , with exactly @xmath2 nearest - neighbor ( @xmath31 ) bonds connected to each site not at a boundary . under pbcs , there are no boundaries , and every site has exactly @xmath23 neighbors .
finite , @xmath25-site sections of these lattices have surface sites with fewer than @xmath32 neighbors and of order @xmath33 zero modes . the free kagome lattice with @xmath34 and @xmath35 unit cells along its sides [ fig .  [ fig : periodic - isostatic](a ) ]
has @xmath36 sites , @xmath37 bonds , and @xmath38 zero modes , all but three of which are floppy modes . these modes , depicted in fig .
[ fig : periodic - isostatic](a ) , consist of coordinated counter rotations of pairs of triangles along the symmetry axes @xmath8 , @xmath9 and @xmath10 of the lattice . there are @xmath34 modes associated with lines parallel to @xmath8 , @xmath35 associated with lines parallel to @xmath10 , and @xmath39 modes associated with lines parallel to @xmath9 .    in spite of the large number of floppy modes in the kagome lattice , its longitudinal and shear lam coefficients , @xmath40 and @xmath41 , and its bulk modulus @xmath42 are nonzero and proportional to the nearest neighbor ( @xmath31 ) spring constant @xmath43 : @xmath44 , and @xmath45 .
the zero modes of this lattice can be used to generate an infinite number of distorted lattices with unstretched springs and thus zero energy @xcite .
we consider only periodic lattices , the simplest of which are the twisted kagome lattices obtained by rotating triangles of the kagome unit cell through an angle @xmath13 as shown in figs .
[ fig : periodic - isostatic](c ) and ( d ) @xcite .
these lattices have @xmath46 rather than @xmath47 symmetry and , like the undistorted kagome lattice , three sites per unit cell . as fig .
[ fig : periodic - isostatic](d ) shows , the lattice constant of these lattices is @xmath12 , and their area @xmath48 decreases as @xmath49 as @xmath13 increases .
the maximum value that @xmath13 can achieve without bond crossings is @xmath50 so that the maximum relative area change is @xmath51 .
since all springs maintain their rest length , there is no energy cost for changing @xmath13 , and as a result , @xmath20 is zero for every @xmath52 , whereas the shear modulus @xmath53 remains nonzero and unchanged .
thus , the poisson ratio @xmath54 attains its smallest possible value of @xmath22 . for any @xmath21 ,
the addition of next - nearest - neighbor ( @xmath55 ) springs , with spring constant @xmath56 ( or of bending forces between springs ) stabilizes zero - frequency modes and increases @xmath20 and @xmath57 .
nevertheless , for sufficiently small @xmath56 , @xmath57 remains negative .
figure [ fig : auxetic - phased ] shows the region in the @xmath58 plane with negative @xmath57 .
phase diagram in the @xmath59 plane showing region with negative poisson ratio @xmath57 . ]
we now turn to the linearized phonon spectrum of the kagome and twisted kagome lattices subjected to pbcs .
these conditions require displacements at opposite ends of the sample to be identical and thus prohibit distortions of the shape and size of the unit cell and rotations but not uniform translations , leaving two rather than three trivial zero modes .
the spectrum @xcite of the three lowest frequency modes along symmetry directions of the undistorted kagome lattice with and without @xmath55 springs is shown in fig .
[ fig : kagome - spectrum](a ) . when @xmath60 , there is a floppy mode for each wavenumber @xmath61 running along the entire length of the three symmetry - equivalent straight lines running from @xmath62 to @xmath63 to @xmath62 in the brillouin zone [ see inset to fig .
[ fig : kagome - spectrum ] ] .
when @xmath64 , there are exactly @xmath65 wavenumbers with @xmath61 along each of these lines for a total of @xmath66 floppy modes . in addition , there are three zero modes at @xmath67 corresponding to two rigid translations , and one floppy mode that changes unit cell area at second but not first order in displacements , yielding a total of @xmath68 zero modes rather than the @xmath69 modes expected from the maxwell count under fbcs .
this is our first indication of the importance of boundary conditions .
the addition of @xmath55 springs endows the floppy modes at @xmath60 with a characteristic frequency @xmath70 and causes them to hybridize with the acoustic phonon modes [ fig .
[ fig : kagome - spectrum](a ) ] @xcite .
the result is an isotropic phonon spectrum up to wavenumber @xmath71 and gaps at @xmath63 and @xmath62 of order @xmath72 .
remarkably , at nonzero @xmath13 and @xmath60 , the mode structure is almost identical to that at @xmath73 and @xmath74 with characteristic frequency @xmath75 and length @xmath76 . in other words ,
twisting the kagome lattice through an angle @xmath13 has essentially the same effect on the spectrum as adding @xmath55 springs with spring constant @xmath77 .
thus under pbcs , the twisted kagome lattice has no zero modes other than the trivial ones : it is  collectively " jammed in the language of references @xcite , but because it is not rigid with respect to changing the unit cell size , it is not strictly jammed .
( a)phonon spectrum for the undistorted kagome lattice .
dashed lines depict frequencies at @xmath60 and full lines at @xmath74 .
the inset shows the brillouin zone with symmetry points @xmath63 , @xmath62 , and @xmath78 .
note the line of zero modes along @xmath79 when @xmath60 , all of which develop nonzero frequencies for wavenumber @xmath80 when @xmath74 reaching @xmath81 on a plateau beginning at @xmath82 defining a length scale @xmath83 .
( b ) phonon spectrum for @xmath84 and @xmath60 .
the plateau along @xmath79 defines @xmath85 and its onset at @xmath86 defines a length @xmath87 . ]
to understand the origin of the differences in the zero - mode count for different boundary conditions , we turn to an elegant formulation @xcite of the maxwell rule that takes into account the existence of redundant bonds ( i.e. , bonds whose removal does not increase the number of floppy modes @xcite ) and states in which springs can be under states of self - stress . consider a ring network in two dimensions shown in fig .  [ fig : self - stress ] with @xmath88 nodes and @xmath89 springs with three springs of length @xmath7 and one spring of length @xmath90 .
the maxwell count yields @xmath91 zero modes : two rigid translations , one rigid rotation , and one internal floppy mode  all of which are  finite - amplitude " modes with zero energy even for finite - amplitude displacements .
when @xmath92 , the maxwell rule breaks down . in the zero - energy configuration , the long spring and the three short ones are colinear , and a prestressed state in which the @xmath90-spring is under compression and the three @xmath7-springs are under tension ( or vice versa ) but the total force on each node remains zero becomes possible .
this is called a state of self - stress .
the system still has three finite amplitude zero modes corresponding to arbitrary rigid translations and rotations , but the finite - amplitude floppy mode has disappeared . in the absence of prestress ,
it is replaced by two  infinitesimal " floppy modes of displacements of the two internal nodes perpendicular of the now linear network . in the presence of prestress
, these two modes have a frequency proportional to the square root of the tension in the springs .
thus , the system now has one state of self stress and one extra zero mode in the absence of prestress , implying @xmath93 , where @xmath94 is the number of states of self stress .
this simple count is more generally valid as can be shown with the aid of the equilibrium and compatibility matrices @xcite , denoted , respectively , as @xmath95 and @xmath96 .
@xmath95 relates the vector @xmath97 of @xmath98 spring tensions to the vector @xmath99 of @xmath100 forces at nodes via @xmath101 , and @xmath102 relates the the vector @xmath103 of @xmath100 node displacements to the vector @xmath104 of @xmath98 spring stretches via @xmath105 .
the dynamical matrix determining the phonon spectrum is @xmath106 .
vectors @xmath107 in the null space of @xmath95 , ( @xmath108 ) , describe states of self - stress whereas vectors @xmath109 in the null space of @xmath102 represent displacements with no stretch @xmath104 , i.e. , modes of zero energy . thus the nullspace dimensions of @xmath95 and @xmath102 are , respectively , @xmath94 and @xmath110 . the rank - nullity theorem of linear algebra @xcite states that the rank @xmath111 of a matrix plus the dimension of its null space equals its column number . since the rank of a matrix and its transpose are equal , the @xmath95 and @xmath102 matrices , respectively , yield the relations @xmath112 and @xmath113 , implying @xmath114 . under pbcs ,
locally isostatic lattices have @xmath2 exactly , and the maxwell rule yields
@xmath115 : there should be no zero modes at all .
but we have just seen that both the square and undistorted kagome lattices under pbcs have of order @xmath33 zero modes as calculated from the dynamical matrix , which , because it is derived from a harmonic theory , does not distinguish between infinitesimal and finite - amplitude zero modes .
thus , in order for there to be zero modes , there must be states of self - stress , in fact one state of self - stress for each zero mode .    in the square lattice under fbcs , @xmath116 and @xmath117 , there are no states of self stress , and @xmath118 zero modes depicted in fig .
[ fig : periodic - isostatic](b ) . under pbcs ,
the dimension of the nullspace of @xmath95 is @xmath119 , and there are also @xmath120 zero modes that are identical to those under fbcs .
we have already seen that there are @xmath121 zero modes in the free undistorted kagome lattice .
direct evaluations @xcite ( see text s1 ) of the dimension of the null spaces of @xmath95 and @xmath102 for the undistorted kagome lattice with pbcs yields @xmath122 when @xmath64 .
the zero modes under pbcs are identical to those under fbcs except that the @xmath123 modes associated with lines parallel to @xmath9 under fbcs get reduced to @xmath34 modes because of the identification of apposite sides of the lattice required by the pbcs as shown in fig .
[ fig : periodic - isostatic](a ) .
thus the modes of both the square and kagome lattices do not depend strongly on whether fcbs or pbcs are applied . under rbc s ,
however , the floppy modes of both disappear .
the situation for the twisted kagome lattice is different .
there are still @xmath124 zero modes under fbcs , but there are only two states of self stress under pbcs and thus only @xmath125 zero modes , as a direct evaluation of the null spaces of @xmath95 and @xmath102 verifies ( see appendix for details ) , in agreement with the results obtained via direct evaluation of the eigenvalues of the dynamical matrix @xcite .
all of the floppy modes under fbcs have disappeared .
( a ) ring - network with @xmath126 showing internal floppy mode . ( b )
ring - network with @xmath92 showing one of the two infinitesimal modes . ]
an effective long - wavelength energy @xmath127 for the low - energy acoustic phonons and nearly floppy distortions provides insight into the nature of the modes of the twisted kagome lattice .
the variables in this theory are the vector displacement field @xmath128 of nodes at undistorted positions @xmath129 and the scalar field @xmath130 describing nearly floppy distortions within a unit cell .
the detailed form of @xmath127 depends on which three lattice sites are assigned to a unit cell .
figure [ fig : periodic - isostatic](c ) depicts the lattice distortion @xmath14 for the nearly floppy mode at @xmath63 ( with energy proportional to @xmath131 ) along with a particular representation of a unit cell , consisting of a central asymmetric hexagon and two equilateral triangles , with @xmath132 sites on its boundary .
if sites @xmath15 , @xmath16 , and @xmath17 are assigned to the unit cell , then the distortion @xmath14 involves no rotations of these sites relative to their center of mass , and the harmonic limit of @xmath127 depends only on the symmetrized and linearized strain @xmath133 and on @xmath14 : @xmath134 , \label{eq : full - elastic}\end{aligned}\ ] ] where @xmath135 is the symmetric - traceless stain tensor , @xmath136 , @xmath137 , @xmath138 , @xmath139 , and @xmath140 . the last term
in which @xmath141 is a third - rank tensor , whose only non - vanishing components are @xmath142 , is invariant under operations of the group @xmath47 but not under arbitrary rotations .
the @xmath143 term is the only one that reflects the @xmath46 ( rather than @xmath47 ) symmetry of the lattice .
there are several comments to make about this energy .
the gauge - like coupling in which the isotropic strain @xmath144 appears only in the combination @xmath145 guarantees that the bulk modulus vanishes : @xmath14 will simply relax to @xmath146 to reduce to zero the energy of any state with nonvanishing @xmath144 .
the coefficient @xmath78 can be calculated directly from the observation that under @xmath14 alone , the length of every spring changes by @xmath147 , and this length change is reversed by a homogenous volume change @xmath148 . in the @xmath149 limit , @xmath150 , and
the energy reduces to that of an isotropic solid with bulk modulus @xmath151 if the @xmath152 and @xmath153 terms , which are higher order in gradients , are ignored .
the @xmath153 term gives rise to a term , singular in gradients of @xmath154 , when @xmath14 is integrated out that is responsible for the deviations of the finite - wavenumber elastic energy from isotropy . at small @xmath13
, the length scale @xmath155 appears in several places in this energy : in the length @xmath156 and in the ratios @xmath157 , @xmath158 , and @xmath159 . at length
scales much larger than @xmath155 , the @xmath152 and @xmath153 terms can be ignored , and @xmath14 relaxes to @xmath146 leaving only the shear elastic energy of an elastic solid proportional @xmath160 . at length scales shorter than @xmath155
, @xmath14 deviates from @xmath146 and contributes significantly to the form of the energy spectrum . if @xmath15 , @xmath18 and @xmath19 in fig .
[ fig : periodic - isostatic](d ) are assigned to the unit cell , then @xmath14 involves rotations relative to the lattice axes , and the energy develops a cosserat - like form @xcite that is a function of @xmath161 rather than @xmath14 .    the modes of our elastic energy in the long - wavelength limit ( @xmath162 )
are simply those of an elastic medium with @xmath24 . in this limit , there are transverse and longitudinal bulk sound modes with equal sound velocities @xmath163 and @xmath164 , where @xmath165 is the particle mass at each node and @xmath166 is the mass density .
in addition there are rayleigh surface waves @xcite in which there is a single decay length ( rather than the two at @xmath167 ) , and displacements are proportional to @xmath168 $ ] with @xmath169 for a semi - infinite sample in the right half plane so that the penetration depth into the interior is @xmath170 .
these waves have zero frequency in two dimensions when @xmath24 , and they do not appear in the spectrum with pbcs .
thus this simple continuum limit provides us with an explanation for the difference between the spectrum of the free and periodic twisted kagome lattices . under fbcs
, there are zero - frequency surface modes not present under pbcs .    further insight into how boundary conditions affect spectrum follows from the observation that the continuum elastic theory with @xmath24 depends only on @xmath171 .
the metric tensor @xmath172 of the distorted lattice is related to the strain @xmath173 via the simple relation @xmath174 ; and @xmath175/2 $ ] , which is zero for @xmath176 , is invariant , and thus remains equal to zero , under conformal transformations that take the metric tensor from its reference form @xmath177 to @xmath178 for any continuous function @xmath179 .
the zero modes of the theory thus correspond simply to conformal transformations , which in two dimensions are best represented by the complex position and displacement variables @xmath180 and @xmath181 .
all conformal transformations are described by an analytic displacement field @xmath182 . since by cauchy s theorem , analytic functions in the interior of a domain
are determined entirely by their values on the domain s boundary ( the  holographic " property @xcite ) , the zero modes of a given sample are simply those analytic functions that satisfy its boundary conditions .
for example , a disc with fixed edges ( @xmath183 ) has no zero modes because the only analytic function satisfying this fbc is the trivial one @xmath184 ; but a disc with free edges ( stress and thus strain equal to zero ) has one zero mode for each of the analytic functions @xmath185 for integer @xmath186 .
the boundary conditions @xmath187 and @xmath188 on a semi - infinite cylinder with axis along @xmath189 are satisfied by the function @xmath190 when @xmath191 , where @xmath192 is an integer .
this solution is identical to that for classical rayleigh waves on the same cylinder . like the rayleigh theory
, the conformal theory puts no restriction on the value of @xmath192 ( or equivalently @xmath193 ) .
both theories break down , however , at @xmath194 beyond which the full lattice theory , which yields a complex value of @xmath195 , is needed .
figure [ fig : ky - kx](a ) shows an example of a surface wave . at the bottom of this figure
, @xmath196 is an almost perfect sinusoid .
as @xmath197 decreases toward the surface , the amplitude grows , and in this picture reaches the nonlinear regime by the time the surface at @xmath198 is reached .
figure [ fig : ky - kx](b ) plots @xmath199 as a function of @xmath193 obtained both by direct numerical evaluation and by an analytic transfer matrix procedure @xcite for different values of @xmath13 ( text s1 ) .
the rayleigh limit @xmath200 is reached for all @xmath13 as @xmath201 .
interestingly the rayleigh limit remains a good approximation up to values of @xmath193 that increase with increasing @xmath13 .
the inset to fig .
[ fig : ky - kx ] , plots @xmath202 as a function of @xmath203 and shows that in the limit @xmath204 ( @xmath205 ) , @xmath199 obeys an @xmath13-independent scaling law of the form @xmath206 .
the full complex @xmath207 obeys a similar equation .
this type of behavior is familiar in critical phenomena where scaling occurs when correlation lengths become much larger than microscopic lengths .
the function @xmath208 approaches @xmath209 as @xmath210 and asymptotes to @xmath211 for @xmath212 .
thus for @xmath213 , @xmath214 and for @xmath215 , @xmath216 . as @xmath13 increases , @xmath217 is no longer much larger than one , and deviations from the scaling law result .
the situation for surfaces along different directions ( e.g. , along @xmath218 rather than @xmath198 ) is more complicated and will be treated in a future publication @xcite .
( a)lattice distortions for a surface wave on a cylinder , showing exponential decay of the surface displacements into the bulk .
this figure was constructed by specifying a small sinusoidal modulation on the bottom boundary and propagating lattice - site positions upward to the free boundary at the top under the constraint of of constant lengths and periodic boundary conditions around the cylinder .
distortions near and at the top boundary , which have become nonlinear , are not described by our linearized treatment .
( b ) @xmath219 as a function of @xmath220 for lattice rayleigh surface waves for @xmath221 , in order from bottom to top .
smooth curves are the analytic results from a transfer matrix calculation , and dots are from direct numerical calculations .
the dashed line is the continuum rayleigh limit @xmath200 .
curves at smaller @xmath13 break away from this curve at smaller values of @xmath207 than do those at large @xmath13 . at @xmath222 , @xmath223 diverges at @xmath224 . the inset plots @xmath202 as a function of @xmath225 for different @xmath13 .
the lower curve in the inset(black ) is the @xmath13-independent scaling function of @xmath226 reached in the @xmath227 limit .
the other curves from top to bottom are for @xmath228 , @xmath229 , @xmath230 , and @xmath231 ( chosen to best present results rather than to match the curves in the main figure ) .
curves for @xmath232 are essentially indistinguishable from the scaling limit .
the curve at @xmath233 stops because @xmath234 . ]
the @xmath46 twisted kagome lattice is the simplest of many lattices that can be formed from the kagome and other periodic isostatic lattices .
figures [ fig : otherauxetic ] ( a ) and ( b ) show two other examples of isostatic lattices constructed from the kagome lattice .
most intriguing is the lattice with pgg symmetry .
its geometry has uniaxial symmetry , yet its long - wavelength elastic energy is identical to that of the @xmath46 twisted kagome lattice , i.e , it is isotropic with a vanishing bulk modulus , and its mode structure near @xmath67 is isotropic as shown in fig .  [
fig : otherauxetic ] ( c ) .
thus , this system loses long - wavelength zero - frequency bulk modes of the undistorted kagome lattice to surface modes .
however , at large wavenumber , lattice anisotropy becomes apparent , and ( infinitesimal ) floppy bulk modes appear .
thus in this and related systems , a fraction of the zero modes of the under fbcs are bulk modes that are visible under pbcs , and a fraction are surface modes that are not .
randomly packed spheres above the jamming transition with average coordination number @xmath235 exhibit a characteristic frequency @xmath236 and length @xmath237 and a transition from a debye - like ( @xmath238 ) to a flat density of states at @xmath239 @xcite .
the square and kagome lattices with randomly added @xmath55 springs have the same properties @xcite .
a general  cutting " argument @xcite provides a procedure for perturbing away from the isostatic limit and an explanation for these properties .
however , it only applies provided a finite fraction of the of order @xmath240 floppy modes of a sample with sides of length @xmath241 cut from an isostatic lattice with pbcs are extended , i.e. , have wave functions that extend across the sample rather than remaining localized either in the interior or at the surface the sample . clearly the twisted kagome lattice , whose floppy modes are all surface modes , violates this criterion ; and indeed , the density of states of the lattice with @xmath242 shows debye - behavior crossing over to a flat plateau at @xmath243 . adding next nearest neighbor bonds gives rise to a length @xmath244 and crossover to the plateau at @xmath245 . the pgg lattice in fig .
[ fig : otherauxetic](a ) , however , has both extended and surface floppy modes , so its crossover to the a flat plateau occurs at @xmath246 rather than at @xmath247 or @xmath248 .
( a ) kagome - based lattice with pgg space group symmetry and uniaxial @xmath249 point group symmetry . ( b )
lattice with p6 space symmetry but global @xmath250 point - group symmetry .
( c ) density plot of the spectrum of the lowest frequency branch of the pgg uniaxial kagome lattice .
the spectrum is absolutely isotropic near the origin point @xmath63 , but it has a zero modes on two symmetry related continuous curves at large values of wave - number . ]
our study highlights the rich and remarkable variety of physical properties that isostatic systems can exhibit . under fbcs
, floppy modes can adopt a variety of forms , from all being extended to all being localized near surfaces to a mixture of the two . under pbc s ,
the presence of floppy modes depends on whether the lattice can or can not support states of self stress .
when a lattice exhibits a large number of zero - energy edge modes , its mechanical / dynamical properties become extremely sensitive to boundary conditions , much as do the electronic properties of the topological states of matter studied in    quantum systems @xcite .
the zero - energy edge modes observed in our isostatic lattices are collective modes whose amplitudes decay exponentially from the edge with a finite decay length , in direct contrast to the very localized and trivial floppy modes arising from dangling bonds .
we focussed primarily on high - symmetry lattices derived from the kagome lattice , but the properties they exhibit , namely a deficit of floppy modes in the bulk and the existence of floppy surface modes , are shared by any two - dimensional system with a vanishing bulk modulus ( or the equivalent in anisotropic systems ) .
three - dimensional analogs of the twisted kagome lattice @xcite can be constructed by rotating tetrahedra in pyrochlore and zeolite lattices @xcite and in cristobalites @xcite .
these lattices are anisotropic . with @xmath31 forces only , they exhibit a vanishing modulus for compression applied in particular planes rather than isotropically , but we expect them to exhibit many of the properties the two - dimensional lattices exhibit .
finally , we note that maxwell s ideas can be applied to spin systems such as the heisenberg anti - ferromagnet on the kagome lattice @xcite , and the possibility of unusual edge states in them is intriguing .
we are grateful for informative conversations with randall kamien , andea liu , and s.d . guest .
this work was supported in part by nsf under grants dmr- 0804900 ( t.c.l . and x.m . ) , mrsec dmr-0520020 ( t.c.l . and a.s . ) and jqi - nsf - pfc ( k. s. ) .
in this section , we will provide some of the details for calculating the equilibrium and compatibility matrices @xmath95 and @xmath102 and the dynamical matrix @xmath251 under periodic boundary conditions .
the equilibrium matrix @xmath95 relates the @xmath252 dimensional vector of bond tensions @xmath97 to the @xmath100 dimensional vector @xmath99 of forces at the nodes via @xmath101 .
we label unit cells by the their position vectors @xmath253 , where @xmath254 and @xmath255 are integers and @xmath256 are the primitive lattice vectors .
( here we take the lattice spacing to be one , i.e. , @xmath257 . )
each unit cell has three sites , labeled @xmath15 , @xmath16 , and @xmath17 as shown in fig .
[ fig : lattice1 ] , so we label the @xmath258 sites by @xmath259 , where @xmath260 and @xmath261 is the number of unit cells
. forces at the nodes are denoted by @xmath262 and their @xmath263 components as @xmath264 , where @xmath265 .
they have to be balanced by stretching forces in the springs located on the bonds .
there are six bonds per cell , which we can take to be the six bonds in the distorted hexagon shown in fig .
[ fig : lattice1 ] , oriented parallel to the unit vectors @xmath266 for @xmath267 , where @xmath13 is the twist angle of the twisted kagome lattice . when @xmath73 , these vector reduce to the edge vectors of a symmetric hexagon with @xmath268 .
the @xmath263 bonds , labeled with @xmath269 , are occupied by springs under tension @xmath270 that exert pulling ( for @xmath271 ) forces , direct along the bond vectors @xmath272 , on nodes .
we define vectors @xmath273 . because we have periodic boundary conditions , we can express both @xmath274 and @xmath275 in terms of their fourier components , @xmath276 the equilibrium matrix is diagonal in @xmath277 , so it breaks up into independent @xmath278 blocks for each of the @xmath6 independent vectors @xmath277 . for simplicity , we consider only the case @xmath64 for which @xmath279 and for which the number of sites is @xmath280 .
there are four bonds incident on each site .
the equations relating @xmath262 to @xmath275 follow from fig .
[ fig : lattice1 ] : @xmath281 we map @xmath282 to a single index @xmath165 with @xmath283 and @xmath284 and then fourier transform .
the result is that @xmath95 breaks up into independent @xmath278 blocks @xmath285 for each @xmath277 .
this leads to @xmath286 where @xmath287 the null space of this matrix is easily calculated with the aid of mathematica .    *
case 1 , @xmath73 * : this is the untwisted kagome lattice . when @xmath61 , the null space of @xmath288 is empty unless @xmath277 is perpendicular to one of the primitive lattice vectors : @xmath289 and @xmath290 , in which case there is a single vector , @xmath291 where @xmath292 is the number of cells , for each @xmath192 and @xmath293 in the null space .
when @xmath67 , there are three vectors , @xmath294 , @xmath295 , and @xmath296 in the null space .
these are the @xmath297 limits of @xmath298 , @xmath299 , and @xmath300 , respectively .
we can , therefore , take the set of null - space vectors of @xmath95 to be @xmath298 , @xmath299 , and @xmath299 for the @xmath34 values of @xmath301 , which include @xmath297 .
linear combinations of these vectors are also null - space vectors , and we can use them to construct vectors confined to a single line of bonds .
for example the vector @xmath302 corresponds to a state in which springs on all @xmath15-bond in cells with centers at @xmath303 for any @xmath254 and on all @xmath11-bonds in cells with centers at @xmath304 are all under the same tension .
this is a state of self - stress in which all bonds along a particular straight line parallel to the @xmath8-axis are under tension .
similar states of self - stress for any of the @xmath305 lines of bonds can be constructed .    * case 2 , @xmath306 * : this is the twisted lattice .
the null space of @xmath288 is empty for _ all _ @xmath307 , and it contains only two vectors when @xmath308 .
these correspond to two states of self - stress in which the stress the stress on the six bonds take on both positive and negative values that are identical in all unit cells .
the compatibility matrix @xmath102 relates the @xmath100 displacements @xmath103 to the @xmath309 stretches @xmath310 .
the displacements are labeled the same way as the forces , i.e. as @xmath311 and the stretches in the same way as the bond tensions , i.e. as @xmath312 .
the equations relating @xmath311 and @xmath313 are @xmath314 fourier transforming and using the same @xmath315 notation as for the equilibrium matrix , we obtain @xmath316 where @xmath317 as required by the constraint @xmath318 .
finally , the dynamical matrix is @xmath319 , whose components are @xmath320 \nonumber\\ d_{12}({\mathbf{q } } ) & = & = -\frac{1}{2 } \sqrt{3 } \cos ( 2 \alpha ) \nonumber\\ d_{13}({\mathbf{q } } ) & = & = d_{31}^*({\mathbf{q } } ) = - e^{-i \left(\frac{q_x}{2}+\frac{\sqrt{3 } q_y}{2}\right ) } \sin ^2\left(\frac{\pi } { 6}-\alpha \right ) -e^{-i q_x}\sin ^2\left(\alpha + \frac{\pi}{6}\right ) \nonumber\\ d_{14}({\mathbf{q } } ) & = & d_{41}^*({\mathbf{q } } ) = \frac{1}{2 } e^{-i q_x}\cos \left(\frac{\pi } { 6}-2 \alpha \right ) + \frac{1}{2}e^{-i \left(\frac{q_x}{2}+\frac{\sqrt{3 } q_y}{2}\right ) } \cos \left(2 \alpha + \frac{\pi } { 6}\right ) \nonumber\\ d_{15}({\mathbf{q } } ) & = & d_{51}^*({\mathbf{q } } ) = -2e^{-\frac{1}{2 } i \sqrt{3 } q_y } \cos ^2(\alpha ) \cos \left(\frac{q_x}{2}\right ) \nonumber\\ d_{16}({\mathbf{q } } ) & = & d_{61}^*({\mathbf{q } } ) = i e^{-\frac{1}{2 } i \sqrt{3 } q_y } \sin ( 2 \alpha ) \sin \left(\frac{q_x}{2}\right ) \nonumber\\ d_{22}({\mathbf{q } } ) & = & 2 \sin ^2(\alpha ) + \cos ^2\left(\frac{\pi } { 6}-\alpha \right ) + \cos ^2\left(\alpha + \frac{\pi } { 6}\right ) \nonumber\\ d_{23}({\mathbf{q } } ) & = & d_{32}^*({\mathbf{q } } ) = \frac{1}{2}e^{-i q_x } \cos \left(\frac{\pi } { 6}-2 \alpha \right ) + \frac{1}{2 } e^{-i\left(\frac{q_x}{2}+\frac{\sqrt{3 } q_y}{2}\right ) } \cos \left(2 \alpha + \frac{\pi } { 6}\right ) \nonumber\\ d_{24}({\mathbf{q}})&= & d_{42}^*({\mathbf{q } } ) = -e^{-i \left(\frac{q_x}{2}+\frac{\sqrt{3 } q_y}{2}\right ) } \cos ^2\left(\frac{\pi } { 6}-\alpha \right ) -e^{-i q_x}\cos ^2\left(\alpha + \frac{\pi}{6}\right ) \nonumber\\ d_{2,5}({\mathbf{q } } ) & = & d_{5,2}^*({\mathbf{q } } ) = ie^{-\frac{1}{2 } i \sqrt{3 } q_y } \sin ( 2 \alpha )   \sin \left(\frac{q_x}{2}\right ) \nonumber\\ d_{2,6}({\mathbf{q } } ) & = & d_{6,2}^*({\mathbf{q } } ) = -2e^{-\frac{1}{2 } i \sqrt{3 } q_y } \sin ^2(\alpha )   \cos \left(\frac{q_x}{2}\right ) \nonumber\\ d_{33}({\mathbf{q } } ) & = & 2 \sin ^2\left(\frac{\pi } { 6}-\alpha \right)+ 2 \sin ^2\left(\alpha + \frac{\pi } { 6}\right ) \nonumber\\ d_{34}({\mathbf{q } } ) & = & d_{43}^ * ( { \mathbf{q } } ) = 0 \nonumber\\ d_{3,5}({\mathbf{q } } ) & = & d_{5,3}({\mathbf{q } } ) = -e^{-i \left(\frac{\sqrt{3 } q_y}{2}-\frac{q_x}{2}\right ) } \sin ^2\left(\frac{\pi } { 6}-\alpha \right ) -e^{i q_x}\sin ^2\left(\alpha + \frac{\pi}{6}\right ) \nonumber\\ d_{3,6}({\mathbf{q } } ) & = & d_{6,3}^*({\mathbf{q } } ) = -\frac{1}{2 } e^{i q_x}\cos \left(\frac{\pi } { 6}-2 \alpha \right ) -\frac{1}{2}e^{i\left(\frac{q_x}{2}-\frac{\sqrt{3 } q_y}{2}\right ) } \cos \left(2 \alpha + \frac{\pi } { 6}\right ) \nonumber\\ d_{4,4}({\mathbf{q } } ) & = & 2 \cos ^2\left(\frac{\pi } { 6}-\alpha
\right)+ 2 \cos ^2\left(\alpha + \frac{\pi } { 6}\right ) \nonumber\\ d_{4,5}({\mathbf{q } } ) & = & d_{5,4}^*({\mathbf{q } } ) = -\frac{1}{2 } e^{iq_x}\cos \left(\frac{\pi } { 6}-2 \alpha \right ) -\frac{1}{2}e^{i\left(\frac{q_x}{2}-\frac{\sqrt{3 } q_y}{2}\right ) } \cos \left(2 \alpha + \frac{\pi } { 6}\right ) \nonumber\\ d_{4,6}({\mathbf{q } } ) & = & d_{6,4}^*({\mathbf{q } } ) = -e^{i \left(\frac{q_x}{2}-\frac{\sqrt{3 } q_y}{2}\right ) } \cos ^2\left(\frac{\pi } { 6}-\alpha \right ) -e^{i q_x}\cos ^2\left(\alpha + \frac{\pi}{6}\right ) \nonumber\\ d_{5,5}({\mathbf{q } } ) & = & \sin ^2\left(\frac{\pi } { 6}-\alpha \right)+ \sin ^2\left(\alpha + \frac{\pi } { 6}\right)+2 \cos ^2(\alpha ) \nonumber\\ d_{5,6}({\mathbf{q } } ) & = & d_{6,5}^*({\mathbf{q } } ) = \frac{1}{2}\cos \left(\frac{\pi } { 6}-2 \alpha \right ) + \frac{1}{2}\cos \left(2 \alpha + \frac{\pi } { 6}\right ) \nonumber\\ d_{6,6}({\mathbf{q } } ) & = & 2 \sin ^2(\alpha ) + \cos ^2\left(\frac{\pi } { 6}-\alpha \right ) + \cos ^2\left(\alpha + \frac{\pi } { 6}\right)\end{aligned}\ ] ] its eigenvalue spectrum can easily be calculated with the aid of mathematica .
the results are shown in fig .  3 of the main text .
surface modes decay exponentially into the bulk .
they are characterized by a wavevector @xmath321 parallel to the plane of the surface , their frequency @xmath322 , their decay length @xmath323 perpendicular to the surface . in lattices with nearest - neighbor forces only ,
@xmath322 and @xmath323 can be determined by setting @xmath324 , where @xmath277 is the wavevector that appears in the dynamical matrix and @xmath325 is the component of @xmath277 perpendicular to the surface ( @xmath326 is the unit vector perpendicular to the surface ) , setting @xmath327 , requiring @xmath328 = 0,\ ] ] where @xmath329 is the unit matrix , and that the equation of motion of the surface layer be satisfied . in the case of surface modes of with zero frequency ,
the relation @xmath330=0 $ ] determines @xmath331 as a function of @xmath321 .
the surface modes of the twisted kagome lattice have zero frequency , so we need only solve this equation .
we consider only the case in which the surface is perpendicular to the @xmath197-direction and set @xmath332 and @xmath333 .
we find @xmath334 = \frac{9}{64 } a(q)\left[r+r^{-1 } - 2g(q)\right ] \left[r+r^{-1 } - 2g^*(q)\right ] , \label{eq : d(q , r)}\ ] ] where @xmath335 \\
g(q ) & = & \frac{\cos(q/2)[3 - 2 \cos p + \cos p \cos ( 4 \alpha ) ] + 2 i \sqrt{3 } \sin^3 ( 3 q/2 ) \sin ( 4 \alpha)}{2 - \cos q + ( 1 - 2 \cos q ) \cos ( 4 \alpha)}\end{aligned}\ ] ] solving for @xmath111 in eq .
( [ eq : d(q , r ) ] ) and to a lattice spacing of @xmath336 , we obtain @xmath337 .\ ] ] this is the function that is plotted in fig .
5(b ) in the text . note that @xmath338 has both a real and an imaginary part , but that in the limit @xmath339 , @xmath340 .