the study of the spontaneous organization of microscopic objects such as colloids , amphiphiles or protein clusters into mesoscale structures ( mesoatoms ) is an active research area that offers challenges on many fronts @xcite .
one would like to create a rich warehouse of raw materials from which to engineer mesomolecules or bulk materials with novel mechanical , optical or electronic behavior .
a fertile source of potential mesoatoms is provided by the self - assembly of micron or smaller scale colloidal particles on _ curved _ interfaces .
the interplay between spatial curvature and condensed matter order in these systems has proven to be very rich .
spatial curvature leads to novel defect arrays in the ground state that have widespread implications for the fundamental and applied physics of curved phases of matter with crystalline , hexatic or nematic order .    in a recent paper kim _
et al _ reported the formation of toroidal micelles from the self - assembly of dumbbell - shaped amphiphilic molecules @xcite .
molecular dumbbells dissolved in a selective solvent self - assemble in an aggregate structure due to their amphiphilic character .
this process has been observed to yield coexisting spherical and open - ended cylindrical micelles .
these structures change slowly over the course of a week to toroidal micelles which thus appear more stable .
toroidal geometries also occur in microbiology in the viral capsid of the coronavirus _ torovirus _ @xcite .
the torovirus is an rna viral package of maximal diameter between @xmath0 and @xmath1 nm and is surrounded , as other coronaviridae , by a double wreath / ring of cladding proteins .
carbon nanotori form another fascinating and technologically promising class of toroidal crystals @xcite with remarkable magnetic and electronic properties .
the interplay between the ballistic motion of the @xmath2 electrons and the geometry of the embedding torus leads to a rich variety of quantum mechanical properties including pauli paramagnetism @xcite and aharonov - bohm oscillations in the magnetization @xcite .
ring closure of carbon nanotubes by chemical methods @xcite suggest that nanotubes may be more flexible than at first thought and provides another technique of constructing carbon tori .    a unified theoretical framework to describe the structure of toroidal crystals
is provided by the elastic theory of defects in a curved background @xcite .
this formalism has the advantage of far fewer degrees of freedom than a direct treatment of the microscopic interactions and allows one to explore the origin of the emergent symmetry observed in toroidal crystals as the result of the interplay between defects and geometry .
the latter is one of the fundamental hallmarks of two - dimensional non - euclidean crystals and leads to universal features observed in systems as different as viral capsids and carbon macromolecules .    in this paper
we provide a detailed analysis of the structural properties of toroidal crystals .
we show that the ground state has excess @xmath3fold disclination defects on the exterior of the torus and @xmath4fold defects on the interior .
the precise number of excess disclinations , as well as their arrangement , is primarily controlled by the aspect ratio of the torus . since defective regions are physically distinguished locations they are natural places for biological activity and chemical functionalization . a thorough understanding of the surface structure of crystalline assemblages could represent a significant step towards a first - principles design of entire libraries of nano and mesoscale components with precisely determined valence .
such mesoatoms could serve in turn as the building blocks for mesomolecules or bulk materials via spontaneous self - assembly or controlled fabrication .
the paper is organized as follows . in sec .
[ sec:2 ] we review the elastic theory of defects in two - dimensional curved geometries and derive the ground state energy of a toroidal crystal . in sec .
[ sec:3 ] we summarize some fundamentals of the geometry of triangulated tori and discuss how the intrinsically discrete problem of crystallography can be reconciled with the result of the continuous elastic theory presented in the previous section . sec .
[ sec:4 ] is devoted to the analysis of the crystalline structure arising from the solution of the elastic problem . in sec .
[ sec:5 ] we discuss the results of numerical minimization of the potential energy of a system of classical particles interacting via a short - range potential on the surface of a torus in the light of the results of sec .
[ sec:4 ] . finally , in sec .
[ sec:6 ] we specialize our analysis to the case of a `` fat torus '' of aspect ratio one and we show how the curvature singularity at the center of the torus is responsible for a remarkable curvature - driven transition to a disordered state .
some of the results presented here have been already announced in ref . .
a two - dimensional torus of revolution @xmath5 is described in parametric form by : @xmath6 where @xmath7 are the two radii of the torus . the metric tensor @xmath8 ( with determinant @xmath9 ) and the gaussian curvature @xmath10 can be written respectively as : @xmath11 k = \frac{\cos\psi}{r_{2}(r_{1}+r_{2}\cos\psi)}\,.\end{gathered}\ ] ] within the elastic theory of defects on curved surfaces the original interacting particle problem is mapped to a system of interacting disclination defects in a continuum elastic curved background .
disclinations are characterized by their topological or disclination charge , @xmath12 , representing the departure of a vertex from a the @xmath13fold coordination of a perfect triangular lattice .
thus @xmath14 , where @xmath15 is the coordination number of the @xmath16th vertex .
a classic theorem of euler requires the total disclination charge of any triangulation of a two - dimensional riemannian manifold @xmath17 to equal @xmath18 , where @xmath19 is the euler characteristic of @xmath17 . in the case of the torus @xmath20 , and
thus disclinations must appear in pairs of opposite disclination charge ( i.e. @xmath3fold and @xmath4fold vertices with @xmath21 and @xmath22 respectively ) in order to ensure disclination charge neutrality .
( color online ) the standard parametrization of a circular torus of radii @xmath23 and @xmath24 . ]
the total free energy of a toroidal crystal with @xmath25 disclinations can be expressed as @xmath26 where @xmath27 is the 2d young modulus and @xmath28 is the solution of the following poisson problem with periodic boundary conditions : @xmath29 where @xmath30 is the laplace - beltrami operator @xmath31 and @xmath32 is the total topological charge density @xmath33 of @xmath25 disclinations located at the sites @xmath34 together with the screening contribution due to the gaussian curvature @xmath35 of the embedding manifold .    the first term in eq
.   represents the long - range elastic distortion due to defects and curvature . the second term in eq .
is the defect core - energy representing the energy required to create a single disclination defect .
this quantity is related to the short - distance cut - off of the elastic theory and is proportional to the square of the topological charge times a constant @xmath36  @xcite .
finally @xmath37 is the free energy of a flat defect - free monolayer .
the gauss - bonnet theorem  @xcite requires the total topological charge to be zero on the torus : @xmath38 the function @xmath28 can be expressed in the green form : @xmath39 where @xmath40 , the laplacian green function on the torus , satisfies the equation : @xmath41 with @xmath42 the dirac delta function on the torus and @xmath43 is the surface area .
it is easy to prove that the solution of the traditional green - laplace equation does nt exist on closed manifolds like the torus , from which the extra term @xmath44 appearing in eq . .
@xmath44 is the eigenfunction of the laplacian associated with the null eigenvalue ( zero mode ) and is necessary because a pure isolated source ( giving rise to the delta - function on the right - hand side ) has no place to terminate the field on the closed torus .    as usual
the calculation of the green function can be remarkably simplified by conformally mapping the torus to a domain of the euclidean plane via a suitable system of isothermal coordinates .
intuitively the torus is conformally equivalent to a rectangular domain described by a system of cartesian coordinates . to make this explicit
, one can equate the metric of the torus in the coordinates @xmath45 to a conformally euclidean metric in the coordinates @xmath46 : @xmath47 where @xmath48 is a positive conformal factor .
taking @xmath49 and @xmath50 , the coordinate @xmath51 is determined by the differential equation : @xmath52 where @xmath53 , the aspect ratio of the torus , may be taken greater or equal to one without loss of generality . choosing the plus sign and integrating both sides of eq .
we find : @xmath54 taking @xmath55 $ ] , the integral yields : @xmath56 where @xmath57 in the transformed coordinate system @xmath46 the modified green - laplace equation reads : @xmath58 where @xmath59 and @xmath60 are now the euclidean laplacian and delta function .
the function @xmath40 can be expressed in the form : @xmath61 where @xmath62 is the laplacian green function on a periodic rectangle and the angular brackets stand for the normalized integral of the function @xmath62 with respect to the dotted variable : @xmath63 analogously the function @xmath64 is given by @xmath65 and ensures the neutrality property : @xmath66    using standard analysis the green function @xmath62 at the points @xmath67 and @xmath68 of a periodic rectangle of edges @xmath69 and @xmath70 can be calculated in the form : @xmath71 + \frac{1}{2\pi}\log\left|\frac{\vartheta_{1}(\frac{z - z'}{p_{1}/\pi}|\frac{ip_{2}}{p_{1 } } ) }   { \vartheta_{1}'^{\frac{1}{3}}(0|\frac{ip_{2}}{p_{1}})}\right|\,,\end{gathered}\ ] ] where @xmath72 , @xmath73 and @xmath74 is the jacobian theta function and reflects the double periodicity of the torus @xcite . a pedagogical derivation of the green function @xmath62
is reported in appendix [ app:1 ] .
substituting eq . in eq . with @xmath75 and @xmath76
we obtain : @xmath77 where @xmath78 represents the stress field due to the gaussian curvature of the torus and is given by : @xmath79+\frac{r-\sqrt{r^{2}-1}}{r}\,.\ ] ] the function @xmath80 is the stress field at the point @xmath81 arising from the elastic distortion due to a defect at @xmath34 and is given by @xmath82 + \frac{1}{4\pi^{2}r}\log(r+\cos\psi_{k})-\frac{\kappa}{4\pi^{2}}\operatorname{re}\{\operatorname{li_{2}}(\alpha e^{i\psi_{k}})\}\notag\\[7pt ] + \frac{1}{2\pi}\log\left|\vartheta_{1}\left(\frac{z - z_{k}}{\kappa}\bigg|\frac{2i}{\kappa}\right)\right|\ , , \label{eq : gamma_defects}\end{gathered}\ ] ] where @xmath83 is the usual eulerian dilogarithm and @xmath84 a derivation of the functions @xmath78 and @xmath80 is given in appendix [ app:2 ] . integrating the function @xmath28 on the torus gives the elastic energy of an arbitrary collection of disclinations . a detailed analysis of the crystalline structure arising from eqs .
, , and is carried out in sec .
[ sec:3 ] .
before analyzing the defect distribution arising from the elastic energy of eq .  , together with eq .
, it is necessary to understand the geometry of triangulated tori .
reconciling the predictions of a continuum elastic theory with the intrinsically discrete nature of crystallography requires an understanding of the possible lattices that can be embedded on the torus and the associated defects .
the problem of classifying the possible triangulations of the @xmath85torus has received considerable attention from mathematicians , physicists and chemists over the past twenty years .
lavrenchenko @xcite proved in 1984 that all the triangulations of the torus can be generated from 21 irreducible triangulations by certain sequences of operations called vertex splitting  .
after the discovery of carbon nanotubes in 1991 and the subsequent theoretical construction ( later followed by the experimental observation ) of graphitic tori , many possible tessellations of the circular torus have been proposed by the community @xcite . in this section we review the construction of a defect - free triangulated torus and we show how the most symmetric defective triangulations can be generally grouped into two fundamental classes corresponding to symmetry groups @xmath86 and @xmath87 respectively .
( color online ) construction of a defect - free triangulation of the torus . on top planar map of the triangulated torus corresponding to the @xmath88 configuration @xmath89 . on the bottom
@xmath90 chiral torus .
the edges of each one of the six tubular segments has been highlighted in red.,title="fig : " ] + ( color online ) construction of a defect - free triangulation of the torus .
on top planar map of the triangulated torus corresponding to the @xmath88 configuration @xmath89 . on the bottom
@xmath90 chiral torus .
the edges of each one of the six tubular segments has been highlighted in red.,title="fig : " ]    for the sake of consistency with the existing literature we adopt here the language developed to describe the structure of carbon nanotubes .
the structure of a triangulated cylinder can be specified by a pair of triangular lattice vectors @xmath91 and @xmath92 , called the _
chiral _ and _ translation _ vector respectively , which together define how the planar lattice is rolled up . in the canonical basis @xmath93 and @xmath94 ,
the vector @xmath91 has the form : @xmath95 the translation vector @xmath92 , on the other hand , can be expressed as an integer multiple @xmath96 of the shortest lattice vector @xmath97 perpendicular to @xmath91 .
the vector @xmath97 is readily found to be of the form : @xmath98 where @xmath99 denotes the greatest common divisor of @xmath100 and @xmath101 and enforces the minimal length .
the three - dimensional structure of the torus is obtained by connecting the edge @xmath102 of the rectangle in fig .
[ fig : triangulation ] ( top ) to @xmath103 and @xmath104 to @xmath105 .
the edge @xmath102 is then mapped to the external equator of the torus while the edge @xmath104 to the @xmath106 meridian .
the resultant toroidal lattice has characteristic chirality related to the initial choice of the vector @xmath91 . in the nanotubes literature _
armchair _ referes to the lattice obtained by choosing @xmath107 , _ zigzag _ to that obtained for @xmath108 and _ chiral _ to all other lattices .
an example of a @xmath88 chiral torus is shown in fig .
[ fig : triangulation ] ( bottom ) for the case @xmath109 , @xmath110 and @xmath111 .
the chirality is extremely important in graphitic carbon nanotube or nanotori , where it determines whether the electronic behavior of the system is metallic or semiconducting .    by euler
s theorem one can prove that the number of triangular faces @xmath112 and the number of vertices @xmath113 of a triangular toroidal lattice is given by : @xmath114 denoting @xmath115 the area of the rectangle with edges @xmath91 and @xmath92 and @xmath116 the area of a fundamental equilateral triangle , the number of vertices of a defect - free toroidal triangulation is then : @xmath117 the planar construction reviewed above allows only lattices with an even number of vertices .
defect - free toroidal deltahedra with an odd number of vertices are also possible and their construction is generally achieved by assembling congruent octahedral building blocks .
an example of this scheme will be briefly discussed in sec .
[ sec:5 ] for the case @xmath118 and @xmath119 .
we refer the reader to ref . for an comprehensive review of the topic .
the embedding of an equal number of pentagonal and heptagonal disclinations in the hexagonal network was first proposed by dunlap in 1992 as a possible way to incorporate positive and negative gaussian curvature into the cylindrical geometry of carbon tubules @xcite . according to the dunlap construction the necessary curvature
is incorporated by the insertion of `` knees '' ( straight cylindrical sections of the same diameter joined with a kink ) in correspondence with each pentagon - heptagon pair arising from the junction of tubular segments of different chirality ( see fig .  [
fig : knees ] ) .
( color online ) voronoi lattices of a tp@xmath120 prismatic ( top ) and ta@xmath120 antiprismatic ( bottom ) toroids with @xmath121 and @xmath122 . ]    in particular , a junction between a @xmath123 and a @xmath124 tube can be obtained by placing a @xmath4fold disclination along the internal equator of the torus and a @xmath3fold disclination along the external equator .
since the radii of the two segments of a junction are different by construction , the values of @xmath120 and @xmath125 are commonly chosen to minimize the ratio @xmath126 . by repeating the @xmath127 construction
periodically it is possible to construct an infinite number of toroidal lattices with an even number of disclinations pairs and dihedral symmetry group @xmath86 ( where @xmath128 is the total number of @xmath127 pairs , fig  [ fig : lucas - torus ] ) .
the structure of the lattice is described by the alternation of two motifs with crystalline axes mutually rotated by @xmath129 as a consequence of the connecting disclination .
one of the fundamental aspects of dunlap s construction is that all the disclinations are aligned along the two equators of the torus where the like - sign gaussian curvature is maximal .
as we will see below , this feature makes these arrangements optimal in releasing the elastic stress due to curvature .
( color online ) dunlap knees obtained by joining two straight tubular segments with @xmath123 and @xmath124 chirality . [ courtesy of a. a. lucas and a. fonseca , facults universitaries notre - dame de la paix , namur , belgium ] . ]
( color online ) five - fold polygonal torus obtained by joining @xmath130 and @xmath131 tubular segments via ten pairs of @xmath127 rings .
this structure was originally proposed by the authors of ref .   as a possible low - strain configuration for carbon nanotori .
[ courtesy of a. a. lucas and a. fonseca , facults universitaries notre - dame de la paix , namur , belgium ] . ]    another class of crystalline tori with dihedral antiprismatic symmetry @xmath87 was initially proposed by itoh _
et al _ @xcite shortly after dunlap .
aimed at reproducing a structure similar to the c@xmath132 fullerene , itoh s original construction implied ten disclination pairs and the point group @xmath133 .
in contrast to dunlap tori , disclinations are never aligned along the equators in antiprismatic tori , instead being staggered at some angular distance @xmath134 from the equatorial plane .
hereafter we will use the symbol ta@xmath120 to refer to toroidal deltahedra with @xmath128 disclination pairs and @xmath87 symmetry group .
a systematic construction of defected triangulations of the torus can be achieved in the context of planar graphs @xcite .
a topological embedding of a graph in a two - dimensional manifold corresponds to a triangulation of the manifold if each region of the graph is bounded by exactly three vertices and three edges , and any two regions have either one common vertex or one common edge or no common elements of the graph
. the simplest example of toroidal polyhedra with @xmath87 symmetry group , featuring only @xmath3fold and @xmath4fold vertices , can be constructed by repeating @xmath120 times the unit cell of fig .
[ fig : graph1]a .
these _ toroidal antiprisms _   have @xmath135 vertices and can be obtained equivalently from the edge skeleton of a @xmath136fold antiprism by attaching at each of the base edges a pentagonal pyramid and by closing the upper part of the polyhedron with @xmath120 additional triangles . by counting the faces one finds @xmath137 from which @xmath135 .
the simplest polyhedron of this family has @xmath138 and @xmath139 symmetry group ( see top left of fig . [
fig : toroidal - antiprisms ] ) and corresponds to the `` drilled icosahedron '' obtained by removing two parallel faces of an icosahedron and connecting the corresponding edges with the six lateral faces of an antiprism with triangular base ( i.e. a prolate octahedron ) . starting from this family of toroidal antiprisms a number of associated triangulations having the same defect structure
can be obtained by geometrical transformations such as the goldberg inclusion @xcite .
such transformations , popularized by caspar and klug for the construction of the icosadeltahedral structure of spherical viruses @xcite , consist in partitioning each triangular face of the original graph into smaller triangular faces in such a way that old vertices preserve their valence and new vertices have valence six .
the partition is obtained by specifying two integer numbers @xmath140 which define how the original vertices of each triangle are connected by the new edges so that the total number of vertices is increased by a factor @xmath141 .
( color online ) first six toroidal antiprisms obtained by repeating the unit cell of fig .  [ fig : graph1 ] . the first polyhedron on the left is the `` drilled icosahedron '' . ]
a general classification scheme for @xmath87 symmetric tori was provided by berger and avron @xcite in 1995 .
their scheme is based on the construction of unit graphs comprising triangular tiles of different _
generations_. in each generation , tiles are scaled in length by a factor @xmath142 with respect to the previous generation .
this rescaling approximates the non - uniformity of the metric of a circular torus .    ( a ) unit cell for toroidal antiprisms .
@xmath3fold vertices are circled and @xmath4fold vertices are boxed .
( b ) unit cell of a @xmath87 torus in the berger - avron construction .
the graph consists of four generation of tiles and the internal equator of the torus is mapped into the horizontal line passing to the mid - point between the 6th and the 7th vertex . ]
( color online ) unit cells for dunlap toroids of type @xmath143 and @xmath144 according to the classification scheme given here .
highlighted regions correspond to the central polygon . ]
dunlap toroids can be obtained from unit cells such as those shown in fig .
[ fig : graph2 ] . the geometrical properties of these graphs can be described in different ways . a particularly intuitive way , in the spirit of this paper , consists in specifying the distances between @xmath3 and @xmath4fold pairs .
one starts by drawing the smallest convex loop passing through defective sites .
this identifies a central polygon whose upper vertices ( @xmath145 and @xmath146 in fig.[fig : graph3 ] ) have degree five and lower vertices ( @xmath147 and @xmath148 in fig .  [
fig : graph3 ] ) have degree seven . then calling @xmath100 the distance between @xmath3fold vertices @xmath145 and @xmath146 , @xmath101 that between @xmath4fold vertices @xmath147 and @xmath148 and @xmath149 the length of the segment @xmath150 ( @xmath151 for a trapezoid ) , we can express the total number of triangles enclosed by the central polygon as : @xmath152 each @xmath4fold vertex sits at the apex of a diamond - shaped complex of @xmath153 triangles .
each @xmath3fold vertex , on the other hand , is at the apex of a triangular region of @xmath154 triangles .
the graph is completed by a rectangle of height @xmath155 and base of arbitrary length @xmath156 containing : @xmath157 triangles .
the total number of vertices is @xmath158 the final triangulation of the torus is obtained by repeating the prismatic unit cell @xmath159 times and therefore has @xmath160 vertices .
this scheme provides direct information on the arrangement of defective sites .
thus for instance an @xmath161 unit cell ( see top of fig .  [
fig : graph2 ] ) has @xmath3fold vertices separated by two lattice spacings and @xmath4fold vertices by one lattice spacing .
on the other hand the integers @xmath120 and @xmath125 giving the chirality of the two segments of the junction @xmath162 are given directly by as : @xmath163 thus the @xmath143 cell of fig .
[ fig : graph2 ] is obtained from the junction between a @xmath164 and a @xmath165 tubular segment .
( color online ) central polygon ( right ) and rectangular ( zig - zag ) region in our construction scheme of dunlap s toroids . in this example @xmath161 . ]    unit cells for tp(2)@xmath120 and tp(3)@xmath120 toroids . ]
dunlap s toroids are not the only examples of defective triangulations of the torus with dihedral prismatic symmetry group @xmath86 . with the help of numerical simulations
( see sec .
[ sec:5 ] ) we found two other classes whose unit cell is shown in fig .
[ fig : graph4 ] . unlike dunlap s toroids , the @xmath4fold vertices in these prismatic triangulations
are not aligned along the internal equator of the torus , but rather grouped in dimers normal to the equatorial plane .
@xmath3fold vertices are distributed along the external equator in the graph of fig .
[ fig : graph4]a or form a double ring above and below it in the case of the graph fig .
[ fig : graph4]b . toroidal deltahedra obtained by embedding the prismatic graphs of fig .
[ fig : graph4 ] on a circular torus are shown in fig .  [
fig : prismatic - deltahedra ] for the case of a @xmath3fold symmetric toroid with @xmath166 and a @xmath4fold symmetric toroid with @xmath167 . in the rest of the paper we will reserve the symbol tp@xmath120 for dunalp s toroids and refer with tp@xmath120a and tp@xmath120b to the other two classes of toroids with symmetry group @xmath86 and unit cell of as shown in fig .
[ fig : graph4]a and fig .
[ fig : graph4]b respectively .
( color online ) ( a ) tp5a and ( b ) tp7b toroids with @xmath166 and 49 obtained by repeating the unit cells of fig .
[ fig : graph4 ] .
@xmath4fold vertices form dimers normal to the equatorial plane while @xmath3fold vertices are ( a ) distributed along the external equator or ( b ) form a double ring above and below the equatorial plane . ]
all defective triangulations presented so far are characterized by an even number of disclination pairs .
regular tessellations of the torus comprising an odd number of defects pairs are also possible .
such tessellations are obtained by combining segments of prismatic and antiprismatic lattices with a consequent loss of dihedral symmetry .
[ fig : r_3.3_200 ] shows the voronoi diagram of a toroidal lattice , with @xmath168 and @xmath169 vertices , containing 11 disclination pairs . for an angular length of approximately @xmath170
the lattice is a prismatic @xmath133 toroid while in the remaining @xmath171 the local structure is that of an antiprismatic toroid .
the global structure has only bilateral symmetry about a sagittal plane dividing the lattice in two mirror halves and thus point group @xmath172 .
( color online ) voronoi diagram of a toroidal lattice with @xmath168 and @xmath169 vertices .
the lattice exhibits 11 disclination pairs and has @xmath172 symmetry group . ]    in the past few years , alternative constructions of triangulated tori have been proposed as well as novel geometrical and graph - theoretical methods to express the coordinates of their three - dimensional structures ( see for example kirby @xcite , lszl _ at al _
@xcite , diudea _ et al _ @xcite ) . here
we choose to focus on the defect structure associated with the two most important class tp@xmath120 and ta@xmath120 with groups @xmath86 and @xmath87 .
in this section we analyze the crystalline structures arising from the solution of the elastic problem and we show how the interplay between the geometry of the embedding torus , the topology of the lattice , and the mechanical properties of the microscopic units , here encoded in the young modulus @xmath27 and the core energy , lead to a rich variety of structures whose phase - diagram is presented at the end of the section . since the free energy of eq .   is minimized when disclinations best approximate the continuum gaussian curvature of the torus , it is clear that disclinations are most likely to be found in regions of like - sign gaussian curvature .
maximal curvature occurs along the external ( positive ) and internal ( negative ) equators , which thus constitute preferred regions for the appearance of disclinations in the ground state .
( color online ) the screening function @xmath173 for different values of the aspect ratio @xmath174 . ]    to analyze the elastic free energy we start by considering the energies of two opposite sign disclinations constrained to lie on the same meridian . the elastic free energy of this system is shown in fig .  [ fig : redbluetorus ] as a function of the angular separation between the two disclinations .
the energy is minimized for the positive ( @xmath3fold ) disclination on the external equator ( maximally positive gaussian curvature ) and the negative ( @xmath4fold ) disclination on the internal equator ( maximally negative gaussian curvature ) .
the picture emerging from this simple test case suggests that a good _ ansatz _ for an optimal defect pattern is a certain number @xmath175 of equally spaced @xmath176 disclinations on the external equator matched by the same number of equally spaced @xmath22 disclinations on the internal equator .
we name this configuration with the symbol @xmath177 , where @xmath175 stands for the total number of disclination pairs .
@xmath178 where the two pairs of numbers specify the @xmath45 coordinates of the positive and negative disclinations respectively .
a comparison of the energy of different @xmath177 configurations , as a function of aspect ratio and disclination core energy , is summarized in the phase diagram of fig .
[ fig : phase_diagram1 ] .
we stress here that only @xmath177 configurations with @xmath175 even have an embedding on the torus corresponding to lattices of the tp@xmath179 class .
nevertheless a comparison with @xmath180odd configurations can provide additional information on the stability of @xmath180even lattices . for small core energies , moreover , thermally excited configurations with a large number of defects and similar @xmath180polar distributions of topological charge
are expected to exhibit an elastic energy comparable in magnitude with that of these minimal constructions .
( color online ) elastic energy of a @xmath127 disclination dipole constrained to lie on the same meridian , as a function of the angular separation . in the inset ,
illustration of a circular torus of radii @xmath7 .
regions of positive and negative gaussian curvature have been shaded in red and blue respectively . ]
the defect core energy has been expressed here in the form : @xmath181 the core energy @xmath36 of a single disclination depends on the details of the crystal - forming material the corresponding microscopic interactions . a simple phenomenological argument
( see for example ref .  )
gives @xmath182 where @xmath100 the lattice spacing .
taking @xmath183 , with @xmath184 the area of the torus , yields : @xmath185 for a system of order @xmath186 subunits , then , the dimensionless core energy on the left hand side of eq .
is of order @xmath187 .
this estimate motivates our choice of the scale for @xmath188 in fig .
[ fig : phase_diagram1 ] .
( color online ) phase diagram for @xmath177 configurations in the plane @xmath189 . for @xmath190 $ ] and @xmath191 the structure
is given by a @xmath192 configuration with symmetry group @xmath193 . ]    for dimensionless core energies below @xmath194 and aspect ratios @xmath174 between 3.68 and 10.12 the ground state structure is the tp5 lattice corresponding to a double ring of @xmath176 and @xmath22 disclinations distributed on the external and internal equators of the torus as the vertices of a regular decagon ( the @xmath192 configuration ) .
the tp5 lattice has dihedral symmetry group @xmath193 . that this structure might represent a stable configuration for polygonal carbon toroids has been conjectured by the authors of ref .
, based on the argument that the 36@xmath195 angle arising from the insertion of ten pentagonal - heptagonal pairs into the lattice would optimize the geometry of a nanotorus consistently with the structure of the @xmath196 bonds of the carbon network ( unlike the 30@xmath195 angle of the @xmath13fold symmetric configuration originally proposed by dunlap ) . in later molecular dynamics simulations , han @xcite found that a @xmath3fold symmetric lattice , such as the one obtained from a ( 9,0)/(5,5 ) junction ( see fig .
[ fig : lucas - torus ] ) , is in fact stable for toroids with aspect ratio less then @xmath197 . the stability , in this case , results from the strain energy per atom being smaller than the binding energy of carbon atoms .
irrespective of the direct experimental observation of such disclinated toroidal crystals , which is still open , we have shown here , from continuum elasticity , that a @xmath3fold symmetric lattice indeed constitutes a minimum of the elastic energy for a broad range of aspect ratios and defect core energies .
( color online ) two examples of defect free `` crop circle '' toroids with @xmath198 and @xmath199 ( left ) and @xmath200 ( right ) . ]    for small aspect ratios the @xmath3fold symmetric configuration becomes unstable and is replaced by the @xmath201fold symmetric phase @xmath202 . as we mentioned , however , this configuration does nt correspond to a possible triangulation of the torus .
it is likely that the ground sate in this regime consists of ten skew disclination pairs as in the antiprismatic ta@xmath120 lattice .
the latter can be described by introducing a further degree of freedom @xmath134 representing the angular displacement of defects from the equatorial plane : @xmath203 \left.\left ( ( -1)^{2k}(\pi-\delta\psi),\frac{2\pi k}{n}\right)_{1 \le k \le n } \right\}\end{gathered}\ ] ] a comparison of the tp5 configuration and the ta5 configuration is shown in fig .
[ fig : phase_diagram2 ] for different values of @xmath134 .
the intersection points of the boundary curves with the @xmath204axis has been calculated by extrapolating the @xmath205 data points in the range @xmath206 $ ] with @xmath207 . for small @xmath134 and @xmath208
$ ] the prismatic tp5 configuration is energetically favored . for @xmath209 , however , the lattice undergoes a structural transition to the ta5 phase . for @xmath210
the prismatic symmetry of the tp5 configuration breaks down again . in this regime
, however , the elastic energy of both configurations rapidly rises because of the lower curvature and defects disappear .
( color online ) phase diagram of a @xmath3fold symmetric lattice in the plane @xmath205 . for small @xmath134 and @xmath174 in the range @xmath211
$ ] the prismatic tp5 configuration is energetically favored . for @xmath209
the system undergoes a structural transition to the antiprismatic phase ta5 . ]      for aspect ratio @xmath212 the tp5 phase is replaced by a defect free configuration ( @xmath213 in fig .
[ fig : phase_diagram1 ] ) so that configurations with defects are no longer energy minima .
any toroidal crystal with aspect ratio larger than @xmath214 is than energetically favored to be defect - free . in the thin torus limit the ground state structure
is directly related to the simple problem of finding the most efficient packing of congruent equilateral triangles on the torus of a given aspect ratio . given @xmath113 subunits ( vertices )
one seeks the densest packing of equilateral triangles of edge - length @xmath215 on the torus with aspect ratio @xmath174 , such that each vertex has valence six . using the planar construction described in sec .
[ sec:3 ] , the optimal choice of the indices @xmath88 , can be translated into the minimization of the following quantity : @xmath216 obtained by equating the magnitude of the chiral vector @xmath91 with that of the sectional circumference of the embedding torus , under the constraints : @xmath217 n,\,m,\,l\in\mathbb{z } \end{array } \right.\,.\ ] ] this construction successfully predicts the structure of the lattices of fig .  [
fig : cropcircle ] .
so far we have studied the elasticity of toroidal crystals exclusively in terms of interacting topological defects on a rigid toroidal substrate .
thus the elastic strain due to defects and curvature takes the form of pure stretching on the tangent plane of the torus and no out - of - plane deformation takes place . in a more realistic scenario ,
a crystalline torus would undergo both in - plane stretching and out - of - plane bending .
the latter implies an energy cost : @xmath218 with @xmath219 the mean curvature ( see ref .  ) .
the case of defect - free tori is simple enough to incorporate bending in the problem and see what the optimal aspect ratio of a defect - free torus would be as a function of the fppl - von krmn number @xmath220 representing the ratio of the stretching energy scale to the bending rigidity . in absence of defects the only source of stress
is given by the curvature .
thus @xmath221}{2r^{2}}\notag\\ & \hspace{0.47\columnwidth}+\operatorname{li_{2}}(\alpha^{2})-2\big\}\,.\label{eq : stretching}\end{aligned}\ ] ] summing eq .   and and taking the derivative with respect to @xmath174 ( assuming constant area ) , one obtains the following equation for the optimal value of @xmath174 : @xmath222=0\,.\end{gathered}\ ] ] the optimal aspect ratio @xmath174 as obtained form eq .
is shown in fig .
[ fig : optimalr ] as a function of @xmath223 . for @xmath224 ,
when the major contribution to the elastic energy is given by the bending , the optimal geometry is given by the clifford torus ( @xmath225 ) .
if , on the other hand , the in - plane stretching dominates , a `` skinny '' torus ( large @xmath174 ) is energetically favoured .
( color online ) optimal value of the aspect ratio @xmath174 as function of the fppl - von krmn number @xmath220 . for @xmath224
the clifford torus with @xmath225 is optimal .
larger values of @xmath223 favour instead a `` skinnier '' torus . ]      in this section we investigate the possibility of a regime of coexistence between isolated disclinations and grain boundary `` scars '' .
the existence of scars , first predicted in the context of spherical crystallography @xcite and later observed experimentally in spherical droplets coated with colloids @xcite , has become one of the fundamental signatures of dense geometrically frustrated systems .    in the regime of large particle numbers , the amount of curvature required to screen the stress field of an isolated disclination in units of lattice spacing becomes too large and disclinations are unstable to grain boundary `` scars '' consisting of a linear array of tightly bound @xmath127 pairs radiating from an unpaired disclination @xcite . in a manifold with variable gaussian curvature
this effects leads to a regime of coexistence of isolated disclinations ( in regions of large curvature ) and scars . in the case of the torus
the gaussian curvature inside ( @xmath226 ) is always larger in magnitude than that outside ( @xmath227 ) for any aspect ratio and so we may expect a regime in which the negative internal curvature is still large enough to support the existence of isolated @xmath4fold disclinations , while on the exterior of the torus disclinations are delocalized in the form of positively charged grain boundary scars .    to check this hypothesis we compare the energy of the tp5 lattice previously described with that of `` scarred '' configurations obtained by decorating the original toroid in such a way that each @xmath176 disclination on the external equator is replaced by a @xmath228 mini - scar .
the result of this comparison is summarized in the phase diagram of fig .
[ fig : phase_diagram3 ] in terms of @xmath174 and the number of vertices of the triangular lattice @xmath113 ( the corresponding hexagonal lattice has twice the number of vertices , i.e. @xmath229 ) .
@xmath113 can be derived from the angular separation of neighboring disclinations in the same scar by approximating @xmath230 , with @xmath231 the area of a hexagonal voronoi cell and @xmath100 the lattice spacing . when the aspect ratio is increased from 1 to 6.8
the range of the curvature screening becomes shorter and the number of subunits required to destroy the stability of the tp5 lattice decreases . for @xmath232
, however , the geodesic distance between the two equators of the torus becomes too small and the repulsion between like - sign defects takes over . thus the trend is inverted .
( color online ) isolated defects and scar phases in the @xmath233 plane .
when the number of vertices @xmath113 increases the range of the screening curvature becomes smaller than one lattice spacing and disclinations appear delocalized in the form of a @xmath228 grain boundary mini - scar . ]
in this section we report the result of a numerical minimization of a system of @xmath113 point - like particles constrained to lie on the surface of a torus and interacting via a pair potential of the form @xmath234 where @xmath235 denotes the euclidean distance in @xmath236 .
the problem of finding the minimal energy configuration of repulsively interacting points on a 2-manifold has become a standard problem of potential theory and has its paradigm in the classical thomson problem on the sphere .
the choice of the cubic potential is motivated here by the so called `` poppy seed bagel theorem '' @xcite , according to which the configuration of points that minimizes the riesz energy @xmath237 on a rectifiable manifold of hausdorff dimension @xmath238 is uniformly distributed on the manifold for @xmath239 . in the case of a torus of revolution
this implies that for small @xmath240 the points are mostly distributed on the exterior of the torus ( the interior becomes completely empty in the limit @xmath241 ) .
as @xmath240 is increased , however , the points cover a progressively larger portion of the surface .
the distribution becomes uniform for @xmath242 . on the other hand , since the number of local minima of the riesz energy increases with @xmath240 , it is practical to choose a value not much larger than two .
the choice @xmath243 has the further advantage of modelling a real physical system of neutral colloidal particles assembled at an interface @xcite and is therefore suitable for direct comparison with experiments on colloidal suspensions .    to construct low energy configurations we adopt an carefully designed hybrid optimization algorithm named tapping ( ta ) .
like other hybrid algorithms , ta consists in a combination of fast local optimizations and global stochastic moves designed to release the system from the local minimum to which it is confined at the end of a local minimization step .
a more detailed description of our algorithm is reported in appendix  [ app:3 ] .
we study four different aspect ratios : @xmath244 and @xmath245 . for each aspect ratio we consider several different particle numbers up to @xmath246 and each simulation is performed for @xmath247 to @xmath248 ta iterations .
the lowest energies found , as well as the number of defects in the corresponding configuration , are reported in table  [ tab : minima ] for a selected set of systems .
the corresponding lattices are shown in fig .
[ fig : numeric_tori ]    the lattices are best presented using a voronoi construction corresponding to the dual lattice of the delaunay triangulation . here
pentagonal faces are colored in red while heptagonal faces are colored in blue .
the complete set of data produced in our simulations together with a collection of interactive 3d graphics for each low energy configuration found is available on - line  @xcite . for fewer than @xmath249 particles the results of our numerical minimization
are in good agreement with the continuum elastic theory .
in particular for @xmath250 and @xmath251 and @xmath252 , we always find minimal energy configurations consisting of ten @xmath3fold disclinations on the outside of the torus and ten @xmath4fold disclinations in the inside as predicted by the elastic theory in the regime of @xmath253 . for @xmath198 and @xmath254
we also find the lowest energy configurations to be defect free .
for small numbers of particles we do nt expect the continuum approximation to accurately describe the lowest energy structure of the toroidal clusters presented in fig .
[ fig : numeric_tori ] . loosely speaking the limit of validity of the elastic theory can be quantified by requiring the average lattice spacing @xmath255^{1/2}$ ] to be much smaller than the radius @xmath24 of the torus .
this condition requires @xmath113 to be of order @xmath256 particles for a torus with aspect ratio @xmath257 .
remarkably , good agreement between the theory and simulations is found starting from much smaller values of @xmath113 and in some cases ( see the following discussion on the configuration with @xmath257 and @xmath258 ) , we already observe the onset of the ideal behavior predicted by theory for @xmath259 .
the occurrence of a ground state configuration with exact prismatic or antiprismatic symmetry , in particular , is only possible when the number of particles @xmath113 belongs to a specific sequence of `` magic numbers '' described in sec .
[ sec:3 ] .
nevertheless for @xmath113 outside such a sequence it is still possible to observe in the ground state a predominant prismatic or antiprismatic character depending on the aspect ratio .
.[tab : antiprism]maximum and minimum aspect ratio for which the toroidal antiprisms are a global minimum . [ cols="^,^,^",options="header " , ]     @xmath260 ) an interesting feature is observed at @xmath261 .
as in the case of @xmath257 we also find a minimum with @xmath262 symmetry group , but unlike the latter configuration , it belongs to the tp@xmath252a class and has @xmath3fold disclinations along the external equator .
a tp@xmath263b configuration is obtained again for @xmath167 . for @xmath264 and @xmath265 the global minimum is achieved by two spectacular antiprismatic configurations with @xmath266 and @xmath267 symmetry group respectively .
these toroids can be obtained from the toroidal antiprisms discussed in sec .
[ sec:3 ] by splitting   in two such that the total number of vertices is increased by one .
the two newly created vertices have coordination number @xmath268 , with @xmath149 is the coordination number @xmath269 , while two of the @xmath149 neighbors of @xmath269 gain a bond . ] one or more times the initial set of @xmath3fold vertices . thus starting from a @xmath270fold toroidal antiprism with @xmath271 vertices and splitting all @xmath272 @xmath3fold vertices one obtains @xmath273 vertices .
splitting twice all @xmath274 @xmath3fold vertices of a @xmath275fold toroidal antiprism with @xmath276 , on the other hand , we have @xmath277 . for @xmath278 the global minimum
is represented by a fascinating lattice of tp@xmath279a type .
lattices with @xmath280 and @xmath281 resemble very closely the structure of a tp@xmath279 graph while for @xmath282 the lattice has a more antiprismatic character with ten defect pairs .
@xmath119 ) three defect - free configurations are found at @xmath283 and @xmath284 .
the case @xmath118 is a particular example of a defect - free lattice that can not be obtained form the planar construction reviewed in sec .
[ sec:3 ] .
it consists of @xmath285 octahedra connected in the form of a chain .
since each octahedron is attached to other two , it contributes with six faces to the total face count .
thus @xmath286 and @xmath287 .
for @xmath288 we always find configurations with ten disclination pairs as expected from continuum elasticity .
for @xmath289 the regime of coexistence between isolated disclinations and scars described in sec .
[ sec : scars ] is observed .
the delocalization of isolated disclinations into scars , however , does nt take place at each defective site simultaneously and the regime of coexistence between positively charged scars and isolated @xmath4fold disclinations is preceded by a phase with isolated @xmath3 and @xmath4fold disclinations and scars .
for tori with aspect ratio as large as @xmath198 we find defect - free ground states every time it is possible to construct a purely @xmath13valent toroidal graph with the same number of vertices @xmath113 .
we have seen that disclination defects , forbidden in the lowest energy state of a planar crystal , may be energetically favored on a substrate of non - vanishing gaussian curvature .
it is therefore natural to ask whether large curvature can completely destroy crystalline order by driving the proliferation of a sufficiently high density of defects .
the resulting state would be amorphous .
the problem of generating amorphous structures by tiling a two - dimensional curved space with identical rigid subunits has drawn attention over the years , particularly through the connection to the structure of such disordered materials as supercooled liquids and metallic glasses .
since the work of frank @xcite the notion of geometrical frustration arises frequently in investigations of supercooled liquids and the glass transition .
a paradigmatic example is represented by the icosahedral order in metallic liquids and glasses which , although locally favored , can not propagate throughout all of three - dimensional euclidean space .
a two - dimensional analog , consisting of a liquid of monodisperse hard disks in a 2-manifold of constant negative gaussian curvature ( the hyperbolic plane ) was first proposed by nelson and coworkers in 1983 @xcite .
in such a system the impossibility of covering the entire manifold with a 6-fold coordinated array of disks mimics many aspects of the geometrical frustration of icosahedral order in three dimensions . in all these models of geometrical frustration , however , the origin of the disorder is primarily due to the short - range nature of the potential between the subunits . in a more realistic setting ,
part of the frustration is relieved by the fact that hexagonal unit cells can compress in order to match the underlying geometry .
the embedding of a triangular lattice on an axisymmetric torus , provides a particularly suitable playground to study curvature - driven disorder .
when @xmath290 the gaussian curvature on the inside of the torus grows like @xmath291 and diverges on the internal equator at @xmath292 .
we thus expect a high density of defects in the vicinity of the curvature singularity and a resultant loss of the local @xmath13fold bond orientational order .
in this regime the system will have crystalline regions on the outside of the torus and amorphous regions near the curvature singularity .
( color online ) top view of a defect free triangulation of a fat torus with @xmath293 and @xmath294 .
the corresponding elastic energy becomes very high in the interior of the torus where the triangles are more compressed to match the reduction of surface area . ]    in this section we substantiate this claim analytically based on the elastic theory of continuous distributions of edge dislocations on a `` fat '' torus .
our argument is based on the following construction .
as a consequence of the curvature singularity the surface area of an arbitrary wedge of angular width @xmath295 becomes smaller and smaller as the sectional angle @xmath296 increases and vanishes at @xmath292 .
if a defect - free lattice is embedded on such a wedge , bragg rows will become closer and closer as the singularity is approached with a consequent rise in the elastic energy ( see fig .
[ fig : fat_torus ] ) . an intuitive way to reduce
the distortion of the lattice is to recursively remove bragg rows as one approaches the point @xmath292 ( see fig .
[ fig : dislocation ] ) .
this is equivalent to introducing a growing density of edge dislocations .
this dislocation `` cloud '' will ultimately disorder the system by destroying the local @xmath13fold bond orientational order .
one might therefore view the curvature as playing the role of a local effective temperature which can drive `` melting '' by liberating disclinations and dislocations . in two - dimensional non - euclidean crystals at @xmath297 , however , the mechanism for dislocation proliferation is fundamentally different from the usual thermal melting .
while the latter is governed by an entropy gain due to unbinding of dislocation pairs , the amorphization at @xmath297 is due to the adjustment of the lattice to the geometry of the embedding manifold via the proliferation of defects and the consequent release of elastic stress .
a similar phenomenon occurs in the disorder - driven amorphization of vortex lattices in type - ii and high-@xmath298 superconductors @xcite .
( color online ) a schematic example of a dislocation pile - up on a square lattice resulting from the shrinking of the area on a regular wedge of a fat torus . ]
since the shrinking area per plaquette on the inside of the torus necessitates a high density of dislocations we may approximate the dislocation cloud in this region by a continuous distribution of burgers vector density @xmath299 .
minimizing the elastic energy with respect to @xmath299 yields a variational equation from which the optimal dislocation density can be calculated as a function of the ratio @xmath300 between the dislocation core energy @xmath301 and elastic energy scale @xmath302 with @xmath303 .    as a starting point , we calculate the green function @xmath40 in the fat torus limit @xmath290 ( i.e. @xmath304 and @xmath305 , see eq .  ) .
the conformal angle @xmath51 in this limit is : @xmath306 and to leading order of @xmath307 we have : @xmath308 \frac{\kappa}{4\pi^{2}}\operatorname{re}\{\operatorname{li_{2}}(\alpha e^{i\psi})\ } \rightarrow \frac{\kappa}{16\pi^{2}}\,\psi^{2}+\frac{1}{2\pi^{2}}\log\left(\cos\frac{\psi}{2}\right)\,.\end{gathered}\ ] ] to handle the limit of the jacobi theta function we can take @xmath309 , @xmath310 and calculate the limit @xmath311 .
this can be done by using the modular transformation properties of jacobi functions@xcite : @xmath312 thus @xmath313 , @xmath314 and @xmath315 , where : @xmath316 this is easily evaluated by means of the expansion : @xmath317 taking the logarithm and neglecting irrelevant constant terms , we obtain : @xmath318 which finally leads to : @xmath319 with @xmath320 . with the green function in hand
, we can calculate the effect of the curvature singularity at @xmath292 on the distribution of defects .
let @xmath299 be the burgers vector density of the dislocation cloud .
hereafter we work in a local frame , so that : @xmath321 with @xmath322 a basis vector in the tangent plane of the torus whose points are specified by the three - dimensional euclidean vector @xmath323 .
the quantity @xmath299 has to be such that : @xmath324 with @xmath325 the total burger s vector in a generic domain @xmath326 .
because on a closed manifold dislocation lines can not terminate on the boundary , extending the integration to the whole torus we have : @xmath327 since the basis vectors @xmath328 in eq .
have the dimension of length , contravariant coordinates @xmath329 have dimensions of an inverse area .
assuming all defects to be paired in the form of dislocations ( i.e. @xmath330 everywhere ) , the total energy of the crystal reads : @xmath331 where @xmath301 is the dislocation core energy and @xmath332 the function @xmath28 encoding the elastic stress due to the curvature and the screening contribution of the dislocation cloud obeys @xmath333 where @xmath334 is the levi - civita antisymmetric tensor on the torus : @xmath335 the stress function @xmath28 can be expressed in the form @xmath336 with    @xmath337 + 1\,,\label{eq : fat_gamma_screening}\\[7pt ] \frac{\gamma_{s}(\psi,\phi)}{y } = \int d^{2}y\,\epsilon_{k}^{i}\nabla_{i}b^{k}(\bm{y})\,g_{l}(\bm{x},\bm{y})\,.\label{eq : fat_gamma_defects}\end{gathered}\ ] ]    taking advantage of the closeness of the torus we can integrate eq .   by parts
so that : @xmath338 now we want reduce the integral term in eq . to a more friendly functional of @xmath299 , suitable for a variational approach .
given the azimuthal symmetry we assume that all dislocations are aligned along @xmath339 .
even though not necessarily true , we argue this to be a reasonable work hypothesis as well as a solid starting point to capture the essential physics of the fat limit . in this case
@xmath340 can be recast in the form @xmath341\,.\end{gathered}\ ] ] substituting eq .   and in eq .
and minimizing with respect to @xmath342 we can now write the variational equation : @xmath343 by inverting the order of integration in the integral on the right hand side , eq   can be expressed in the form of a fredholm equation of the second kind : @xmath344 where @xmath345 , @xmath346 and the kernel @xmath347 is given by : @xmath348 the function @xmath349 on the right hand side of eq .   is given by : @xmath350\sin\psi-2\operatorname{cl_{2}}(\psi+\pi)\bigr\}\end{gathered}\ ] ] where @xmath351 is the clausen function ( see ref .
1005 - 1006 ) defined as : @xmath352    ( color online ) the burgers vector component @xmath353 for different choices of @xmath354 . ]
as previously noted the dislocation core energy is @xmath301 is much smaller than the elastic energy scale @xmath302 .
eq   is then suitable to be solved in powers of the dimensionless number @xmath354 : @xmath355 the corrections to the zero - order term @xmath356 can be calculated recursively by solving a set of fredholm equations of the first kind : @xmath357 the function @xmath356 associated with the burgers vector density of the dislocation cloud in the limit @xmath358 , on the other hand , can be calculated directly from eq .   by setting the effective topological charge density on the right hand side to zero : @xmath359 for a torus of revolution the only nonzero christoffel symbols are @xmath360 \gamma_{\phi\phi}^{\psi } = r_{2}^{-1}\sin\psi(r_{1}+r_{2}\cos\psi)\,.\end{gathered}\ ] ] since @xmath361 by assumption , the first term in eq .   can be expressed as : @xmath362 & = ( 1+\cos\psi)\partial_{\psi}b^{\phi}-2\sin\psi\,b^{\phi}\,,\end{aligned}\ ] ] and eq .
becomes an ordinary differential equation @xmath363 whose solution is given by @xmath364 so that @xmath365 .
the burgers vector density @xmath342 obtained by a numerical solution of eq .
is shown in fig .
[ fig : burgers ] for different values of @xmath354 .
the burgers vector density is measured in units of @xmath366 the function @xmath342 has cubic singularities at @xmath367 and is approximately zero on the outside of the torus . the solid blue curve in fig .
[ fig : burgers ] represents the zeroth order solution of eq ..    now , in the theory of dislocation mediated melting a system at the solid liquid phase boundary is described as a crystalline solid saturated with dislocations . in three - dimensions , in particular , there is a strong experimental evidence of the existence of a critical dislocation density at the melting point @xmath368 where @xmath101 is the length of the length of the smallest perfect - dislocation burgers vector @xcite .
several theoretical works have motivated this evidence both for three - dimensional solids and vortex lattices in super conductors @xcite . on the other hand , given the existence of such a critical density , its value can be empirically used to determine whether a system is in a solid or liquid - like phase in the same spirit as the lindemann criterion . with this goal in mind
we can calculate the dislocation density by requiring @xmath369 with @xmath370 the density of single lattice spacing dislocations .
this yields : @xmath371 solving @xmath372 as a function of @xmath296 and @xmath113 we obtain the diagram of fig .
[ fig : melting - diagram ] .
as expected the inside of the torus contains an amorphous region whose angular size decreases with the number of vertices @xmath113 as a consequence of the reduction of the lattice spacing .
( color online ) phase diagram for curvature driven amorphization .
the inside of the torus contains an amorphous region whose angular size decreases with the number of vertices @xmath113 as a consequence of the reduction of the lattice spacing . ]
in this article we reported a comprehensive analysis of the ground state properties of torodial crystals . using the elastic theory of defects on curved substrates we identified the ground state structures of an arbitrary crystalline torus as a function of the aspect ratio @xmath174 and the ratio @xmath188 of the defect core energy to the elastic energy scale .
we showed that for a large range of aspect ratios and core energies the minimal energy structure of a toroidal crystal has ten disclinations pairs and symmetry group @xmath193 , as first conjectured by lambin _
et al _ more than ten years ago  @xcite . for large system sizes we proved isolated disclinations
are unstable to grain - boundary scars consisting of chains of tightly bound @xmath127 pairs radiating form an unpaired disclination . on a torus , where the gaussian curvature on the inside is always larger in magnitude than that on the outside , the occurrence of scars is marked by a state featuring isolated @xmath4fold disclinations in the inside of the torus together with scars on the outside . for a torus of aspect ratio close to one we showed how a diverging gaussian curvature on the internal equator is responsible for the remarkable occurrence of a curvature - driven transition of the system to a disordered , liquid - like , state .
the predictions of our elastic theory were compared with the results of a numerical study of a system of point - like particles constrained on the surface of a torus and interacting via a short range potential , with good agreement . from a purely geometrical point of view
we have introduced a number of novel toroidal polyhedra as well as a new construction and classification scheme for certain types of prismatic tori .
this project was supported by nsf grants dmr-0305407 and dmr-0705105 and by an allocation through the teragrid advanced support program .
lg is supported on a graduate fellowship by the syracuse biomaterials institute .
we acknowledge amand lucas and antonio fonseca for sharing with us the beautiful photos presented in sec.[sec:3 ] .
we also would like to thank eric j. west and m. cristina marchetti for inspiring discussions during the development of sec .
[ sec:6 ] .
lg is very grateful to dario giuliani who made possible the realization of the _ toroidal lattices database _
@xcite , the use of which has been crucial for obtaining many of the results and the ideas presented in this paper .
the modified laplacian green function on a periodic rectangle of edges @xmath69 and @xmath70 can be conveniently calculated in the form : @xmath373 where @xmath374 is the eigenfunction of the laplace operator with periodic boundary conditions : @xmath375 such that : @xmath376 u_{\lambda}(\xi,0 ) = u_{\lambda}(\xi , p_{2 } ) \end{array } \right.\,.\ ] ] in cartesian coordinates the eigenfunctions are simple plane waves of the form : @xmath377 where @xmath378 and @xmath379 are given by : @xmath380 and the eigenvalue @xmath354 is given by : @xmath381 calling for simplicity @xmath382 and @xmath383 , the function @xmath384 is given by : @xmath385 although is simple , it is very useful to rewrite it in terms of elliptic functions .
noting that the odd terms in cancel we have :    @xmath386 & = -\frac{2}{p_{1}p_{2 } } \left [    \sum_{m=1}^{\infty}\frac{\cos\frac{2\pi m}{p_{2}}(y-\eta)}{\left(\frac{2\pi m}{p_{2}}\right)^{2}}+    \sum_{n=1}^{\infty}\sum_{m=-\infty}^{\infty}\frac{\cos\frac{2\pi n}{p_{1}}(x-\xi)\cos\frac{2\pi m}{p_{2}}(y-\eta ) }     { \left(\frac{2\pi n}{p_{1}}\right)^{2}+\left(\frac{2\pi m}{p_{2}}\right)^{2 } } \right]\,.\end{aligned}\ ] ]    an equivalent expression can be obtained by isolating the @xmath108 contribution in the sum rather than the @xmath387 term . the first sum in eq .
can be evaluated easily by using : @xmath388 thus we have : @xmath389 & = -\frac{1}{2\,p_{1}p_{2}}\left(\frac{p_{2}^{2}}{6}-p_{2}|y-\eta|+|y-\eta|^{2}\right)\ , .
\end{aligned}\ ] ] the second sum in eq . can be evaluated with the help of the poisson summation formula : @xmath390 in particular , if we choose : @xmath391 we can write the second sum in eq .   as :    @xmath392 & = -\frac{2}{p_{1}p_{2 } } \sum_{n=1}^{\infty}\cos\frac{2\pi
n}{p_{1}}(x-\xi ) \sum_{k=-\infty}^{\infty}\int_{-\infty}^{\infty}dt\,\frac{\cos(2\pi k+p)t } { \left(\frac{2\pi n}{p_{1}}\right)^{2}+\left(\frac{2\pi t}{p_{2}}\right)^{2}}\end{aligned}\ ] ]    the integral can be easily calculated by considering : @xmath393 thus eq .   becomes : @xmath394 the sum in @xmath120 can be calculated by noting : @xmath395 where @xmath396 with @xmath397 and arbitrary choice of sign .
thus separating positive and negative values of @xmath398 in eq .   and using eq .
the function @xmath399 takes the final form : @xmath400 + \sum_{k=1}^{\infty}\log\left|1 - 2q^{2k}\cos\frac{2\pi}{p_{1}}(z-\zeta)+q^{4k}\right|\,,\end{gathered}\ ] ] where : @xmath401 the second term in eq .
can be expressed in terms of the jacobi theta function @xmath402 defined as : @xmath403 another useful relation can be obtain by taking the derivative of @xmath402 with respect to @xmath404 : @xmath405 thus we have : @xmath406 in this way we can express : @xmath407\,.\ ] ] then taking : @xmath408 and substituting in eq .
we obtain : @xmath409 -\frac{\pi}{p_{1}}|y-\eta| + \log\left|\frac{\vartheta_{1}(u , q)}{\vartheta_{1}'^{\frac{1}{3}}(0,q)}\right|\,.\end{gathered}\ ] ] combining eq .   with eq .
we conclude that : @xmath410 & = \frac{\log 2}{6\pi } - \frac{1}{2\,p_{1}p_{2}}|y-\eta|^{2 } + \frac{1}{2\pi}\log\left|\frac{\vartheta_{1}(u , q)}{\vartheta_{1}'^{\frac{1}{3}}(0,q)}\right|\,.\end{aligned}\ ] ] an alternative notation frequently used for the jacobi theta function is : @xmath411 with this choice we can write the green function in the final form : @xmath412 + \frac{1}{2\pi}\log\left|\frac{\vartheta_{1}(\frac{z-\zeta}{p_{1}/\pi}|\frac{ip_{2}}{p_{1 } } ) }   { \vartheta_{1}'^{\frac{1}{3}}(0|\frac{ip_{2}}{p_{1}})}\right|\,.\end{gathered}\ ] ]
in this appendix we derive the analytical expression for the stress functions @xmath413 and @xmath80 in eq .  .
the former is given by the integral : @xmath414\ ] ] it is convenient to keep the green function @xmath62 in the form of eq .  .
thus we have : @xmath415 which using eq .   and carrying out the integrals yields : @xmath416\ ] ] to calculate the second integral in eq .
it is convenient to invert the order of integration and use the result of eq .  :
@xmath417 & = \log\left[\frac{4(r^{2}-1)}{(r+\sqrt{r^{2}-1})^{2}}\right]-\frac{r-\sqrt{r^{2}-1}}{r}\,.\end{aligned}\ ] ] combining eq .   and
we obtain : @xmath418 + \frac{r-\sqrt{r^{2}-1}}{r}\,.\ ] ]    the elastic stress produced at the point @xmath81 by a disclination at @xmath34 is given by the function : @xmath419 to calculated the average @xmath420 one can again start from @xmath62 in the form of a series and use eq .  .
this yields @xmath421 to calculate the integrals in eq .
it is convenient to expand the conformal angle @xmath51 in fourier harmonics : @xmath422 with : @xmath423\ ] ] where @xmath424 .
thus we have : @xmath425 with @xmath426 where @xmath427 stands for the real part and @xmath428 is the usual euler s dilogarithm ( see ref
.   pp 1004 - 1005 ) .
the second integral in eq .
is give by : @xmath429\,.\end{gathered}\ ] ] to calculate the integral in eq .
one uses parseval s identity : @xmath430 where @xmath431 an arbitrary square - integrable function on the interval @xmath432 $ ] with fourier series : @xmath433 thus we have : @xmath434\,.\end{aligned}\ ] ] using eq .  , and we conclude that : @xmath435 & -\frac{\kappa}{8\pi}\operatorname{li_{2}}(\alpha^{2 } )   + \frac{\kappa}{4\pi^{2}}\operatorname{re}\{\operatorname{li_{2}}(\alpha e^{i\psi})\}\notag\\[5pt ] & + \frac{1}{4\pi^{2}r}\log\left[\frac{2(r^{2}-1)}{(r+\cos\psi)(r+\sqrt{r^{2}-1})}\right]\,.\end{aligned}\ ] ]
the tapping algorithm ( ta ) is a hybrid algorithm designed to find the optimal crystalline structure of particle systems constrained to lie on a curved surface and interacting with a long range potential of the form @xmath436 .
hybrid algorithms , such as basin - hopping @xcite and minima - hopping @xcite , have been successfully employed throughout the years to predict the crystalline structure of molecular clusters and proteins . in general
they combine fast local minimizations with global moves whose goal is to release the system from the local minimum it is confined at the end of a local minimization step .
a typical hybrid optimization routine can be summarized in the following two steps : _ 1 _ ) after all the independent variables have been randomly initialized , a local optimization is performed and a local minimum @xmath437 is determined ; _ 2 _ ) from @xmath437 a new configuration @xmath438 is constructed by applying a global ( generally stochastic ) move . the new configuration @xmath438 is then used as starting point for a new local optimization step .
the two steps are iterated until a stopping criterion is satisfied .
the goal is thus to explore the largest possible number of local minima and avoid visiting the same minimum too often .
the crucial point in designing an effective hybrid algorithm is clearly the choice of the global move .
there is no general rule to identify a successful global transformation @xmath439 and physical intuition and prior experience are typically the only guidelines . in the case of basin - hopping , for instance
, the global transformation consists in a monte carlo move in which all the particles of the system are randomly displaced in order to construct a new initial configuration from which a new trial minimum is obtained .
the step is accepted with probability @xmath440 , where @xmath441 is the energy difference between the new and previous minimum and @xmath442 is an inverse temperature adjusted to obtain a @xmath443 acceptance ratio . in the case of minima - hopping the escape step is performed by a short molecular dynamics simulation by assigning the particles a fixed kinetic energy .    the global move adopted in ta is inspired by the process of close packing of spherical objects by tapping and is motivated by the well established role of topological defects in determining the order of two - dimensional non - euclidean crystals as well as the picture of the potential energy surface ( pes ) of such systems as a multi - funnel landscape .
consider a system of say spherical objects confined in a two - dimensional box with an initial disordered configuration . a common way to bias the system toward a close - packed configuration is to provide it kinetic energy by gently tapping the box .
if the system is populated by locally ordered regions ( i.e. grains ) separated by clusters of defects , the primary effect of tapping is to produce a glide of defects inside the crystals with a subsequent rearrangement of grains .
this mechanism can be reproduced numerically in the following way .
the algorithm starts with a random distribution of particles and rapidly quenches the system by performing a fast local minimization .
once particles are trapped in a local minimum , defects are identified by a delaunay triangulation of the lattice .
then the system is tapped by adding to the defect positions a random displacement .
the magnitude of the displacement is given by the typical spacing associated with the particles number times a factor @xmath354 which represents the tapping strength .
this factor is initially set to @xmath444 .
after defects have been moved a new local minimization is performed in order to construct the trial configuration @xmath438 .
the energy of this configuration is compared with the energy of the previous minimum and the move is accepted if their difference is larger than some tolerance factor @xmath445 . if , on the other hand , the energy difference is smaller than @xmath445 , the system has relaxed again to the same minimum . in this case the tapping strength is increased of a factor 10 and the process is repeated until the system successfully hops to a new minimum . the tapping strength @xmath354 is then set to its initial value .
the process is iterated until the rate of discovery of new global minima drops below some threshold value or a maximum number of iterations is reached .    in the current implementation of the algorithm
, the local minimization step is performed using the fletcher - reeves conjugate gradient algorithm @xcite .
analytic expressions for the energy gradient and the hessian matrix are coded in the program in order to reduce the number of evaluations of the objective function during the relaxation step to one single event .
the delaunay triangulation is calculated via the dwyer s divide and conquer algorithm with alternate cuts  @xcite , which runs in @xmath446 time , making the identification of the defects particularly fast .
the main difference between ta and other hybrid algorithms ( including basin - hopping ) is that the escape move consists of adaptive displacements of defects only , rather than of the entire system . in the case of non - euclidean crystals , where the conformation of the energy landscape is subtly related to the arrangement of topological defects ,
this mechanism is believed to explore the pes more accurately . in systems as lennard - jones clusters or spin - glasses , the pes is characterized by an exponential number of local minima separated by energy barriers . for this reason the majority of the algorithms are specifically designed to allow the system to overcome a barrier by providing it a significant amount of energy . if the energy landscape , however , is characterized by the presence of multiple narrow funnels , as believed in this case , the previous methods become ineffective .
a funnel represents the basin of attraction of a given local minimum .
if the global minimum is also located at the bottom of a funnel , an algorithm that is attempting to locate it via a sequence of local minimization steps has a chance to find it exclusively by starting from a configuration already at the muzzle of the funnel .
such possibility , however , is ruled out if all the particles are displaced simultaneously during the escape move and the system is abruptly moved to a completely different place in the energy landscape . on the other hand , by adaptively tapping the defects it is possible to achieve a much finer inspection of the pes and possibly locate the funnel associated with the global minimum .
a copy of our code is available by request .
whitesides and m. boncheva , _ beyond molecules : self - assembly of mesoscopic and macroscopic components _ , proc .
usa * * 9**9 , 4769 - 4774 ( 2002 ) [ www.pnas.org/cgi/doi/10.1073/pnas.082065899 ] .
a. r. bausch , m. j. bowick , a. cacciuto , a. d. dinsmore , m. f. hsu , d. r. nelson , m. g. nikolaides , a. travesset and d. a. weitz , science * 299 * , 1716 ( 2003 ) .
t. einert , p. lipowsky , j. schilling , m. j. bowick and a. bausch , langmuir * 21 * , 12076 ( 2005 ) .
d. r. nelson , nano lett . * 2 * , 1125 ( 2002 ) .