%The paper contains a large number of figures that are available upon
%request.
BODY:
\documentstyle[12pt]{article}
\newcounter{num}
%\baselineskip 20pt
%\magnification 1200
\renewcommand{\baselinestretch}{1.8}
%\font\tenrm=cmr10
%\font\ninerm=cmr9
%\font\twelvermb=cmbx12
%\def\today{\ifcase\month\or
% January\or February\or March\or April\or May\or June\or
% July\or August\or September\or October\or November\or
% December\fi\space\number\day, \number\year}
%{\nopagenumbers
%\headline={\ifnum \count0=1 { \ninerm\hfil\today}\else\hfil\fi }
%\footline={\ifnum \count0=1 {~\hfil}\else \hfil\fi}
%{~} \bigskip
%\twelvermb
\title{A Classification Scheme for Toroidal Molecules}
%\tenrm \vskip 2.0cm
\author{Jorge Berger and Joseph E.~Avron \\
Department of Physics,
Technion - Israel Institute
of Technology \\
32000 Haifa, Israel\\
e-mail: PHR76JB@VMSA.TECHNION.AC.IL }
\begin{document}
\addvspace{-25 mm}
\maketitle
\noindent We construct a class of periodic tilings of the plane,
which corresponds to toroidal arrangements of trivalent atoms, with
pentagonal, hexagonal and heptagonal rings. Each tiling is characterized by a
set of four integers and gives rise to an
infinite family of toroidal molecules. These may be chosen to fit a desired
shape, size or number of atoms in the molecule. The configurational
energy and the delocalisation energy of several molecules
obtained in this way were computed
for Tersoff and H\"uckel models. The results indicate that many of these
molecules would be stable. We studied the influence of size on
the H\"{u}ckel spectrum: it bears both similarities and differences as
compared with the case of tubules.
\newpage
\section{Introduction}
The last decade witnessed the discovery of diverse \cite{Dov} hollow
graphitic structures. These may be classified into sphere-like
\cite{Fullerenes} and cylinder-like \cite{Tubules}. A natural question is
whether additional topologies, such as tori, also occur. It is
reasonable to expect that molecules that are multiply connected and have
negative curvature (i.e. have regions that approximate a saddle) would have
distinct physical and chemical properties. For example,
certain quantum and transport phenomena depend on multiconnectivity
\cite{Rieman}, and it
is known that negative curvature has important consequences for the
dynamics on the surface \cite{Chaos}. Toroidal molecules may have
applications in host-guest chemistry \cite{hostguest}, where
the ``hole'' in the torus accomodates the ``guest". It
has also been suggested that toroidal molecules would have applications
as components in nanotechnology \cite{torus} and as catalizers
\cite{Ijima}.
Structures that are actually half tori have been observed
\cite{Ijima} when two concentric cylinders inter-connect.
In the present study we shall focus on structures in the shape of a
doughnut. Miscellaneous proposals for toroidal molecules have been
suggested \cite{torus,hexagons,Chern}.
The purpose of the present study is to present a classification
scheme for a wide class of toroidal molecules . We shall consider
symmetric (mainly ${\rm D}_{nd}$) molecules.
Our tool is a set of tiling rules \cite{Tiling}, adapted to the
description of symmetric noncrumpled toroidal molecules. These
tiling rules reflect the energy penalty for
stretched bonds and deformed bond angles. Our method is
amenable for computarized search (as the ``orange-peel" method \cite{Peel}
in the case of spheroids). For simple properties, such as the spectrum
in H\"{u}ckel
model, it suffices to give the indices which characterize a molecule as
the input to a short program and the results are automatically obtained.
The tiling provides also metric information.
An announcement of some of the results in the present study is found in
\cite{A-B}.
As an application we describe results of a numerical study of several
toroidal carbon molecules. The study involved:
\begin{enumerate}
\item The relaxation of model carbon atoms represented by tiling to
a stable molecular configuration, through minimization of the
$\sigma-$bonds geometric potential energy. We have used a
phenomenological potential due to Tersoff \cite{Tersof}. From here,
average energies and binding energies of the least bound
atoms were obtained.
\item Homo-lumo and other properties of the $\pi$-electrons in
H\"uckel model.
\end{enumerate}
The results (Table I) support the conjecture that many
toroidal carbon molecules would be stable.
\section{Tilings}
A toroidal carbon molecule can be represented by a bi-periodic
planar graph, where atoms sit at vertices and each vertex is trivalent.
Since carbon bonds prefer angles near $120^\circ$, graphs
with pentagonal, hexagonal and heptagonal rings are favored and we allow
only these. In hexagons-only tori, as considered in \cite{hexagons},
the hexagons must be deformed in order
to reproduce the curvature of the toroidal surfaces, leading to
prohibitively large energies.
We shall represent a molecule by its {\it dual graph}, i.e. a
periodic planar triangular graph with vertices of valences $5,\ 6$ and
$7$. These triangles are called {\it tiles}, and the atoms are
located at the centers of each tile.
We are interested in smooth tori that approximate the surface of
revolution with cartesian representation
\begin{eqnarray}
R\left\{\left(1-\frac{\cos\theta}{\eta}\right)\cos\phi \ ,
\left(1-\frac{\cos\theta}{\eta}\right)\sin\phi \ ,
\frac{\sin\theta}{\eta}\right\} %1
\end{eqnarray}
in the toroidal coordinates $\{\theta,\phi\}$. The parameter
$R>0$ determines the size, and the parameter $\eta > 1$ determines
the shape. (Large $\eta$ corresponds to thin tori.) Eq.~(1)
induces the metric and area
\begin{eqnarray}
&&\eta^2 R^{-2}ds^2=d\theta^2+(d\phi)^2 \
(\eta-\cos\theta)^2 \ ; \nonumber \\ %2
&&\eta^2 R^{-2}dA = d\theta d\phi (\eta-\cos\theta) \ .
\end{eqnarray}
We shall refer to lines of fixed $\theta$ (resp. $\phi$) as latitudes
(resp. longitudes) and draw them horizontal (resp. vertical).
Our set of tiles includes various sizes,
reflecting the non uniform metric in the $\{\theta,\phi\}$
plane. (Large tiles are placed where $\cos\theta$ is large and vice
versa.) For simplicity, we shall take all our tiles as similar
isosceles triangles. In order to accomodate
vertices with valence 5, 6 and 7, these triangles have to be right-angled.
Due to the privileged role of latitudes and longitudes, we shall (initially)
require that each tile have at least one side along one of these directions.
More general structures will be discussed in Section 4.
It follows by Pythagoras that the different tiles are related by
scaling by powers of $\sqrt{2}$. The largest tiles will be called the
{\em first generation}; the tiles in the {\em $i$-th generation}
will have an area which is $2^{i-1}$ times smaller.
In the first generation we shall have the two (big) isosceles right-angled
tiles of area one in the tiling plane, shown in Fig.~1a. We
allow only the two orientations (related by a $\pi$ rotation)
with the hypotenuse along a
latitudinal line; this will be required from every odd generation,
giving rise to an array of equidistand horizontal lines. (This
feature agrees with metric (2), which
depends on $\theta$, but not on $\phi$. In this way the
first generation tiles give the shortest latitudinal lengths.)
At the second generation we add the four
tiles of area $\frac{1}{2}$ in Fig.~1b. At generation $p+2$, we add
the tiles obtained from those in generation $p$ by scaling by a factor
of $\frac{1}{2}$. Clearly, tiles that share a common edge
differ by at most one generation.
Since there is only one way of matching
odd generation tiles, the tiling of an odd-generation region
is uniquely determined by its border; on the other hand,
since a square can be tiled in two ways (see Fig.~2), the tiling
of an even-generation region
involves a choice of orientation for each square.
Since the linear dimension of the last generation is smaller than
that of the first by a factor $2^{(m-1)/2}$, and since the first generation
has the hypotenuse along a latitudinal line (what may introduce
an additional factor $2^{1/2}$ in the metric along the latitudinal
direction), Eq.~(2) can be used to relate $\eta$ to the number of generations
$m$:
\begin{eqnarray}
\frac{\eta+1}{\eta-1} \sim 2^{[m/2]} \ , %3
\end{eqnarray}
where $[m/2]=$ integral part of $m/2$.
Therefore, the
number of generations plays a role analogous to $1/\eta$.
A cylinder can be built with one generation, skinny tori
are associated with first and second generation tiles and, the fatter
the tori, the more generations are needed.
Unlike the sphere, which is unique up to scaling, the tori of Eq.~(1)
depend on the shape parameter $\eta$. {\it A priori}, it is not clear
whether this family has a distinguished member. It turns out that there
is such a torus and the case $\eta = \sqrt{2}$ is special. In the theory
of surfaces the torus of revolution with $\eta = \sqrt{2}$ is known as the
Clifford torus and is distinguished by being the minimizer of the
Willmore functional, i.~e., it minimizes the square of the mean curvature
\cite{Bensimon}.
The fattest torus we shall study has five generations, inspired by
the Clifford torus. As we shall see, it does have a marked stability.
For fatter tori,
axially symmetric structures are not expected \cite{Bensimon}.
Valence plays the role of
curvature: valence 5 (resp. 6, 7) vertices carry positive (resp. zero,
negative) curvature.
We want no 5 and 7 valence vertices embedded in single generation regions.
The pair of tiles of the odd generations have indeed
this behavior built into them: only sixfold vertices can be made
with them alone.
On the other hand, Fig.~2 shows a tiling by even
generation tiles with valences 5 and 7.
These pentagon-heptagon pairs do not contribute to the global curvature
required by metric (2), but merely produce crumpling, and should be
avoided by an additional rule; pentagons or heptagons will be allowed
only at the interfaces between different generations.
The connection between the tiling and the Riemanian curvature can be
made more stringent. The predominance of pentagons over heptagons
should be proportional to the curvature. In our case, it should be
a decreasing function of $\cos\theta$ and, also, of the size of the
tile. However, since there is an energy penalty for every
additional pentagon-heptagon pair, and since we are interested in
molecules with numbers of atoms which are not too large, we shall
compromise the smoothness of the curvature on behalf of a limited
number of pentagons and heptagons. Accordingly, pentagons (resp. heptagons)
will be located only where the
tiles are smallest (resp. largest). In summary, we impose the
following:\\
\noindent {\em 5-6-7 Rule}: For a torus with $m$ generations, all the
heptagons lie at the interfaces between the 1st and the 2nd generation
and all the pentagons lie at the interfaces between the $m$-th and the
$(m-1$)-th generation. (If there are only two generations,
both the heptagons and the pentagons lie at the interfaces between
them.) Furthermore, we require that around every heptagon (resp.
pentagon) there be a majority of tiles of generation 1 (resp. $m$).\\
Figs. 3-5 are examples of tilings which obey the 5-6-7 rule.
\noindent {\em Length}: We define the length $L(\theta)$ of a latitudinal
line as the number of edges it contains. Clearly,
$L(\theta)$ is larger when the line passes through high generation regions.
In Fig.~5, $L(\theta)$ varies from 1 to 3.
Likewise, we define the length of a longitudinal line as the number
of horizontal lines it intersects. Our tiling rules ensure that the
total length of longitudinal lines is independent of $\phi$, as it
should. This length will be dubbed the ``girth'' and denoted by
$g$. We also define $g_i(\phi)$ as the length within the $i$-th
generation regions. Note that if a longitudinal line crosses
several regions of the same generation $i$, then $g_i(\phi)$
is defined as the {\em total} length across all these regions.
In Fig.~5, $g_1 =g_3 =2$ and $g_2 =g_4 =1$ (independent of $\phi$).
Clearly, $\/ \sum^m_{i=1} g_i(\phi)=g$.\\
\section{Skeletons and Relationships}
Given a set of interfaces between generations, the entire tiling is
determined uniquely by the 5-6-7 rule (although not for every set of
interfaces does a legal tiling exist). The set of interfaces will be called
the {\it skeleton}. The ($2 p-1,\ 2p)$ interfaces are zigzag lines (e.g.
Fig.~6a) and the $(2p, \ 2p+1$) interfaces are horizontal lines
(e.g. Fig.~6b). Therefore, regions that contain $2p-1$ and $2p$
tiles lie within a horizontal stripe, which will be called a
{\em $p$-th stripe}. (If there are only 2 generations, there is just
one stripe; if the number $m$ of generations is odd, then the
last stripe has tiles of the $m$-th generation only.)
Figs.~3 and 4 have one stipe; Fig.~5 has two.
Since according to metric (2) the length of the latitudinal
lines has only one minimum and one maximum, we require that the
regions of generation $i$, $1\leq i\leq m$, be located in monotonic
and consecutive order, both when going from region 1 to $m$ by
increasing or by decreasing $\theta$. The 5-6-7 rule imposes a very
strong constraint on the class of possible skeletons with two or more
stripes: with the exception of
the vertices at which there is a majority of generation 1, and, if
$m$ is even, the vertices with a majority of generation
$m$, all the vertices of the zigzag interfaces lie at the stripe
borders and all the zigzag interfaces are parallel to each other (see
Fig.~7). This has
two consequences: First, the tiling obeys the symmetry ${\rm D}_{nd}$, which
is the highest symmetry of a torus in three dimensions
that could be expected from a discrete construction. The mirrors of the
symmetry are at the planes of constant $\phi$ that contain the
pentagons and the heptagons. The centers of inversion lie midway
between the heptagons (or pentagons) which are located at different
zigzag lines. Second, the $g_i$'s are independent of $\phi$, in
agreement with the desired metric (2).
If the skeleton consists of just one stripe, then the
symmetry ${\rm D}_{nd}$, the constancy of the $g_i$'s and the zigging
between two lines of constant latitude don't follow from the 5-6-7
rule (e. g. Figs.~8-9), but we impose them in order to mimic the
metric~(2).
The independence of the $g_i$'s on $\phi$ rules out the tori
considered in Ref. \cite{Chern}, which are cylindrical tubes connected by
elbows.
We are thus left with a class of tilings such that the skeleton
(and the entire tiling) is completely determined by four numbers:
$m$ (the number of generations), $g_1$ and $g_m$ (the contributions
of the first and last generation to the girth), and $z$ (the number of
edges in a zig of the (1,2) interface). If $i \ne 1,m$, then
$g_i=2^{p-1}z$, where $p$ is the number of the stripe.
The girth of the torus is given by
\setcounter{num}{4}
\setcounter{equation}{0}
\def\theequation{\thenum\alph{equation}}
\begin{eqnarray}
g=g_1+g_{\rm mid} + g_{\rm end} + g_m \ , %4a
\end{eqnarray}
with
\begin{eqnarray}
g_{\rm mid} = \left\{
\begin{array}{ll}
0 & m \leq 2 \\
z & 2> z$. Since
experimentally found ``turn around" edges have a small difference between
the external and internal radius (relative to the radii themselves),
elongation as in Fig.~8 is necessary. The outer cylinder contains tiles
of generation 2, connections contain the zigzagging boundaries, and
the inner cylinder contains tiles of generation 1 and separating
fringes of generation 2. From geometric arguments and
comparison with Eqs.~(9) and (10), we estimate $g_1$ and $g_2$ as $ 0.4
\times ($length of the cylinders), $nz$ as $2.5 \times ($difference in
radii of cylinders) and the total horizontal witdth of the separating
fringes as $5 \times ($the internal radius) minus $2.5 \times ($the
external radius), with all lengths in \AA. Note that in this case
Eq.~(1) is not a good approximation
and that interlayer interaction has not been taken into account.
The primitive tori ${\rm C}_{120}$ and ${\rm C}_{240}$ in \cite{torus}
may be traced to be $g_1 =2, \ g_4 =z=1$ and the leapfrogged of
$g_1 =g_3 =2, \ z=1$, respectively.
For the purpose of illustration, we have focused on small values for
$g_1,g_m$ and $z$. Molecules with longitudinal
perimeters of the order of that of C$_{60}$ imply $g \sim 10$.
For tori with roughly circular cross section, we may also expect that the
difference
in the latitudes of the pentagon and the heptagon lying at the same
longitude be at least about $2\pi/3$.
Fig.~11 shows three dimensional views of the molecule $m=2, \ g_1=g_2=z=4$
after relaxation to the minimal energy, and Fig.~12 shows $m=5,\ g_1 =
g_5=2, \ z = 1$. This is the five-generations molecule with the smallest
possible tiling. The difference in $\eta$ is apparent, with Fig.~12 being
close to a Clifford torus \cite{Bensimon}.
Comparison of different levels of inflation of $m=g_1=g_2=z=2$ shows that
the configurational energy per atom relative to graphite decreases
roughly as $1/N$. The same result was obtained in \cite{torus}.
However, the energy of the least bound atoms saturates.
This effect is probably due to the fact that the number of pentagons and
heptagons is kept unchanged, and they have to bear the entire burden of
curvature. In Fig.~11a it is clearly seen that pentagons stick out as
protuberances. It is also interesting to note in Fig.~11b how the
pentagons deviate from the latitude $\theta = 0$ which they have in the
tiling.
Fig.~13 shows a chiral molecule, as obtained by shifting by $j=1$ the tiling
of the molecule in Fig.~12 at $\theta^* = \pi $.
The influence of inflation on the H\"{u}ckel spectrum was studied for
$g_1=g_2=2$ and $1 \leq z \leq 3$. The results are summarized in Figs.~14
and 15. Fig.~14 describes the homo-lumo gap. Since there are $N$ energy
levels, and since they are bound in the interval $[3\beta,-3\beta]$, the gap
might naively be expected to be inversely proportional to $N$, and thus to
$\ell^2$. The trend is, however, that the gap is inversely proportional to
$\ell$. This result resembles the case of tubules, where the gap is
inversely proportional to the diameter \cite{diameter}. There are some
inflation values for which the regular trend is not obeyed. These are a
reminiscent of gap closure for tubules whose width is divisible by 3.
Fig. 15 shows how the delocalisation energy per atom approaches that of
graphite. The difference from graphite is inversely proportional to
$\ell^2$.
\section{Conclusion}
Tori are closed surfaces, like spheres,but unlike spheres are not
simply connected. They may provide an arena for a variety of physical and
chemical properties.
Table I shows that, from the energetic point
of view, there are toroidal molecules which are roughly as stable as
C$_{60}$. One could say that the fact that toroidal
carbon molecules have not been observed is as significant as the fact
that C$_{60}$ molecules were not observed until recently. The energies and
geometries we have found are similar to those in \cite{torus}, in spite
of having used a different potential and minimization technique.
We have developed a classification scheme which can be used to ``design''
toroidal molecules with prescribed dimensions. As long as these molecules
are not found, it may help to suggest candidates; if they are found, it
may help to analyze them.
{\bf Acknowledgments:}
We are grateful to M.~Kaftory for his help with Figs.~11-13 and
to M.~Dresselhaus, L.A.~Chernozatonskii and T.~Schlick for correspondence. The
research is supported by GIF and DFG through SFB288, and the Fund for the
Promotion of Research at the Technion.
\newpage
\begin{thebibliography}{article}
\bibitem{Dov} E. $\bar{\rm O}$sawa, M. Yoshida, and M. Fujita, MRS Bulletin
XIX(11) 33 (1994).
\bibitem{Fullerenes} H.W. Kroto, J.R. Heath, S.C. O'Brien, R.F. Curl and
R.E. Smalley, Nature {\bf 318} 162 (1985); H.W. Kroto, Nature
{\bf 329} 529 (1987); W.~Kr\"atschmer, L.D.~Lamb, K.~Fostiropoulos, and
D.R.~Huffman, Nature {\bf 347} 354-7 (1990); H.W.~Kroto, A.W. Allaf and
S.P. Balm, Chem. Rev. {\bf 91} 1213-35 (1991); R.F.~Curl and R.E.~Smalley,
Sci.~Am.~October 1991,pp 54-61.
\bibitem{Tubules} P.W.~Fowler, J. Chem. Soc. Faraday Trans., {\bf
86} 2073-7 (1990); S.~Ijima, Nature {\bf 354} 56-8 (1991); T.W.~Ebbesen
and P.M.~Ajayan, Nature {\bf 358} 220-2 (1992); Z.Ya.~Kosakovskaya {\it et
al.}, JETP Lett. {\bf 56} 26-30 (1992); N.~Hamada, S.~Sawada, and
A.~Oshiyama, Phys. Rev. Lett. {\bf 68} 1579-81 (1992);
A.A.~Lucas,P.H.~Lambin, and R.E.~Smalley, J.~Phys.
Chem. Solids {\bf 54} 587-93 (1993); S.~Iijima and T.~Ichihashi, Nature
{\bf 363}, 603 (1993); D.~S.~Bethune {\it et al.}, Nature {\bf 363}, 605
(1993).
%\bibitem{Heptagons} In D.~Vanderbilt and J.~Tersoff,
%Phys. Rev. Lett. {\bf 68} 511-3 (1992) an intriguing proposal for a
%negatively curved periodic surface with hexagonal and heptagonal carbon
%rings was studied.
\bibitem{Rieman} A. Shapere and F. Wilczek, {\it Geometric Phases in
Physics},
World Scientific, Singapore, (1989); J.~E.~Avron, M.~Klein, A.~Pnueli and
L.~Sadun, Phys.~Rev.~Lett.~{\bf 69}, 128-131, (1992).
\bibitem{Chaos} N.~Balasz and A.~Voros, Phys.~Rep. {\bf 143} 103, (1986);
see also {\em Proceedings of the Houches Summer School}, Eds. M.J.
Gianoni and A.~Voros, North Holland, Amsterdam (1989).
\bibitem{hostguest}
See e.g. {\it Host-Guest Molecular Interactions: From Chemistry to Biology},
Ciba Foundation Symposia, Wiley, (1991).
\bibitem{torus} S. Itoh, S. Ihara
and J. Kitakami, Phys. Rev. B {\bf 47} 1703 (1993);
S.~Ihara and S.~Itoh, Phys. Rev. B {\bf 47}, 12908 (1993);
S.~Itoh and S.~Ihara, Phys. Rev. B {\bf 48}, 8323 (1993).
\bibitem{Ijima} S.~Iijima P.~M.~Ajayan and T.~Ichihashi, Phys. Rev. Lett.
{\bf 69}, 3100 (1992); S.~Iijima, Mater. Sci. Engineering B {\bf 19},
172 (1993).
\bibitem{hexagons} T.G.~Schmalz, W.A.~Seitz, D.J.~Klein, and G.E.~Hite,
J.~Am.~Chem.~Soc. {\bf 110}, 1113-27 (1988); E.~C.~Kirby, R.~B.~Mallion and
P.~Pollak, J.~Chem. Soc. Faraday Trans., {\bf 89}, 1945 (1993).
\bibitem{Chern} L.A.~Chernozatonskii, Phys.~Lett.~A {\bf 170}, 37
(1992); B.I. Dunlap, Phys. Rev. B, {\bf 46} 1933 (1992);
B.~I.~Dunlap, Phys. Rev. B, {\bf 49} 5643 (1994);
E.G.~Gal'pern, I.V.~Stankevich, A.L.~Chistyakov and L.A.~Chernozatonskii,
Fullerene Sience \& Technology, {\bf 2}, 1-11, (1994).
\bibitem{Tiling} B.~Gr\"unbaum and G.C.~Shepard, {\it Tiling and
Patterns}, W.H. Freeman, N.Y. (1987).
\bibitem{Peel} D.~E.~Manolopoulos, J.~C.~May and S.~E.~Down,
Chem. Phys. Lett. {\bf 181}, 105 (1991).
\bibitem{A-B} J.E.Avron and J.Berger, Phys. Rev. A, {\bf 51}, xxxx (1995).
\bibitem{Tersof} J.~Tersoff, Phys.~Rev.~B {\bf 37} 6991-7000 (1988);
J.~Tersoff, Phys.~Rev.~Lett. {\bf 61} 2879 (1988).
\bibitem{Bensimon} U. Pinkall and I. Sterling, {\em Willmore Surfaces},
The Mathematical Intelligencer,
{\bf 9}(2)(1987); U.~Seifert, Phys. Rev. Lett. {\bf 66},2404 (1991);
T.~J.~Willmore, {\em Riemmannian Geometry} (Clarendon,Oxford,1993);
X.~Michalet, F.~J\"{u}licher,B.~Fourcade, U.~Seifert et
D.~Bensimon, La Recherche {\bf 25} 1012-1018 (1994).
\bibitem{Gold} M.Goldberg, Tohoku Math. J. {\bf 43},104 (1937).
%\bibitem{hepta} There is evidence thatthe necks in certain microtubules
%have heptagonal rings.
%\bibitem{Sah} P.W. Fowler and Woolrich, Chem.~Phys.~Lett. {\bf 127} 78-83
%(1986); D. Babic, D.J. Klein and C.H. Sah, Chem. Phys. Let. {\bf 211}
%235-41 (1993).
\bibitem{Leapfrog} P.W.~Fowler, Chem. Phys. Lett. {\bf 131} 444-50 (1986);
P.W. Fowler and J.I. Steer, J. Chem. Soc., Chem. Commun. 1403-5 (1987).
\bibitem{MWF} D.E. Manolopoulos, D.R. Woodwall and P.W.
Fowler, J. Chem Soc., Farad. Trans. {\bf 88}, 2427-2435 (1992).
\bibitem{Dress}
R.~Saito, M.~Fujita,G.~Dresselhaus, and M.S.~Dresselhaus, Appl.
Phys. Lett. {\bf 60} 2204-6 (1992).
%R. Saito, G.~Dresselhaus and M.S.~Dresselhaus, Chem.
%Phys. Lett. {\bf 195} 537 (1992).
\bibitem{helices} S.~Ihara, S.~Itoh and J. Kitakami,
Phys. Rev. B {\bf 48}, 5643 (1993); X.~B.~Zhang {\it et al.},
Europhys. Lett., {\bf 27}, 141 (1994).
\bibitem{Powel} M.J.D.Powell, Mathematical Programming {\bf 12} 241-254
(1977).
\bibitem{Graph} {\em Chemical applications of Graph theory}, A.T. Balaban,
ed., Academic Press, N.Y. (1976); D.M. Cvetkovic, M.Doob and H. Sachs, {\em
Spectra of Graph theory and applications}, Academic Press, N.Y. (1979).
\bibitem{diameter} C.T. White, D.H. Robertson and J.W.Mintmire, Phys. Rev. B
{\bf 47}, 5485 (1993); R.~Saito,G.~Dresselhaus and M.~S.~Dresselhaus,
J.~Appl. Phys. {\bf 73}, 494 (1993); C.~H.~Olk and J.~P.~Heremans,
J.~ Mater. Res. {\bf 9},259 (1994).
\end{thebibliography}
\newpage
%\bigskip
\begin{center}
Table I
\end{center}
\begin{tabular}{lcccccccc}
$m$ & 2 & 2 & 2 & 2 & 3 & 4 & 5 & $5^1$ \\
$g_1$ & 2 & 4 & 6 & 2 &4 & 2 & 2 & 2\\
$g_m$ & 2 & 4 & 6 & 2 &2 & 1 & 2 & 2\\
$z$ & 2 & 4 & 6 & 3 &6 & 1 & 1 & 1\\
$N$ & 120 & 480 & 1080 & 180 &1200 & 120 & 240
& 240 \\
$\epsilon_h(\beta)$ & 0.43 & 0.16 & 0.08 &0.45 &0.07 & 0.04
& 0.19 & 0.20\\
$\epsilon_\ell(\beta)$ & -0.27~ & -0.22~ & -0.02~ &-0.18~
&-0.02~ & 0.04 & 0.01 & -0.04~ \\
$10^3(\bar{\epsilon}-\bar{\epsilon}_G) (|\beta |)$ & 16 & 4 & 2 &11
&2 & 25 &11 & 9\\
$10^3(\bar{E} - E_G)/|E_G|$ & 123 & 32 & 16 & 118 & 16 & 76
& 39 & 43\\
$10^3(E_{\rm worst} - E_G)/|E_G|$ & 181 & 114 & 108 & 170 & 115 &
122 & 101 & 100\\
$N(\bar{E} - E_G)/|E_G|$ & 14.8 & 15.6 & 17.5 & 21.2 & 19.5 & 9.1
& 9.4 & 10.4\\
inner radius ($\rm \AA$) & 4.05 & 8.33 & 12.57~ & 6.87 &
14.15~ & 2.02 & 2.09 & 2.12\\
outer radius ($\rm \AA$) & 7.29 & 14.21~ & 21.17~ & 10.14~
&23.42~ & 6.02 & 8.49 & 8.61\\
height ($\rm \AA)$ & 2.93 & 5.68 & 8.00 & 3.37 & 8.13 & 4.57 & 6.89 & 6.96
\end{tabular}
\pagebreak
%\centerline{Table II}
%
%\begin{tabular}{cccc}
%& graphite&C$_{60}$&dodecahedron \\
%\\
%atoms per molecule&$\infty$&60&20 \\
% $\epsilon_h$($\beta$) & 0 & 0.62 & 0 \\
% $\epsilon_\ell$($\beta$) & 0 & $-0.14$~ & 0 \\
%$\frac{N(\epsilon_\ell-\epsilon_n)}{6|\beta|}$
%& - & 7.6~ & 0 \\
%$\bar{\epsilon}(\beta)$ & 1.57 & 1.55 & 1.47 \\
%$\bar{E} - E_G$ & 0 & 0.67 & 1.62 \\
%radius ($\AA$) & $\infty$ & 3.69 & 1.93
%\end{tabular}
%\pagebreak
\begin{center} Table II \end{center}
\begin{tabular}{|cc|ccc|ccc|} \hline
\verb+\+ & $g_2$ & 1 & 2 & 3 & 1 & 2 & 3 \\
$z$ & \verb+\+ & &$g_1=2$ & & & $g_1=4$ & \\ \hline
1 & & 35 & 17 & 22 & 36 & ~0 & ~1 \\
2 & & 61 & 70 & ~8 & 11 & 51 & 13 \\
3 & & 71 & 62 & 22 & 44 & 27 & 21 \\ \hline
\end{tabular}
\pagebreak
\begin{center}
Table Captions
\end{center}
\noindent {\bf Table I}
Energetic and geometric values for several toroidal molecules, according to
H\"uckel's theory and to Tersoff's potential. All the molecules have 5 unit
cells. The second and third columns are inflations
of the first. $N$ is the number of atoms in the molecule,
$\epsilon_{h,\ell}$ is the energy of the highest
occupied/lowest unoccupied molecular orbital, $\bar{\epsilon}$ is the
average energy of delocalized electrons, $\bar{E}$ is the average geometric
energy per atom, $E_{\rm worst}$ is the geometric energy
associated with the least bound atom, and
$\bar{\epsilon}_G=1.575 \beta$ and
$E_G=-7.40{\rm eV}$ are the delocalisation and
the geometric energy per atom of an ideal graphitic plane.
The molecule in the last column (denoted by $m=5^1$) is
obtained by a chiral shift $j=1$ of $(m,g_1,g_m,z,n)=(5,2,2,1,5)$
\noindent {\bf Table II} $10^2 \times
(\epsilon_{\ell}-\epsilon_h)/|\beta|$ for molecules with 2 generations
and 5 unit cells, for various values of $g_1, \ g_2$ and $z$.
\newpage
\centerline{Figure Captions}
\noindent {\bf Fig. 1} (a) Odd generation tiles. (b) Even generation tiles.
\noindent {\bf Fig. 2} Even generation tiling with unnecessary crumpling.
\noindent {\bf Fig. 3} The tiling $m=g_1=g_2=z=2$.
\noindent {\bf Fig. 4} The tiling $m=g_1=g_2=2, \ z=3$.
\noindent {\bf Fig. 5} The tiling $m=4, \ g_1=2, \ g_4=z=1$.
\noindent {\bf Fig. 6} (a) Interface between even generation higher than
the odd generation. (b) Interface between odd generation higher than
the even generation.
\noindent {\bf Fig. 7} Skeleton of tiling that obeys the 5-6-7 rule and has
more than one stripe. (In the figure $m=6$ and there are 3 stripes.)
The numbers denote the generation of the tiles in the
region. $g_1$, $g_m$ and $z$ determine completely the unit cell of the
tiling.
\noindent {\bf Fig. 8} The tiling $g_1=2, \ g_2=z=1$, after elongation by 1.
\noindent {\bf Fig. 9} Three unit cells of the tiling $g_1=g_2=z_1=2,
\ z_2=1$.
\noindent {\bf Fig. 10} Unit cell of the tiling $g_1=2, \ g_2=z=1$, after
a Goldberg inclusion with indices $\ell_1=1, \ \ell_2=2$. The circles signal
the original vertices before the inclusion.
\noindent {\bf Fig. 11} Three-dimensional view of the molecule
$(m,g_1,g_2,z,n)=(2,4,4,4,5)$, which minimizes Tersoff's potential. For
clarity, only atoms in the foreground are shown. (a) View from ``above".
(b) Side view.
\noindent {\bf Fig. 12} Same as Fig. 11, for $(m,g_1,g_2,z,n)=(5,2,1,1,5)$.
\noindent {\bf Fig. 13} Same as Fig. 12, with a chiral shift $j=1$.
\noindent {\bf Fig. 14} Gap in the H\"uckel spectrum for molecules
$(m,g_1,g_2,z,n)=(2,2\ell,2\ell,\ell,5)$, $(2,2\ell,2\ell,\ell,6)$,
$(2,2\ell,2\ell,2\ell,5)$, $(2,2\ell,2\ell,2\ell,6)$,
$(2,2\ell,2\ell,3\ell,5)$, $(2,2\ell,2\ell,3\ell,6)$. For $z/\ell=1$
(resp. 2, 3), calculations were carried up to $\ell=9$ (resp. 8, 5).
The numbers in the graph stand for $z/\ell$. Calculations for $n=6$ were
carried up to $\ell=5$; for the values of $\ell$ for which a number appears
only once, the values of the gap for $n=5$ and $n=6$ are either identical
or indistinguishable within the resolution of the graph. For every $z/\ell$,
$\ell(\epsilon_{\ell}-\epsilon_h)$ exhibits approximate periodicity as a
function of $\ell$, with period 3.
\noindent {\bf Fig. 15} Delocalisation energy for the same molecules studied
in Fig. 14. $\bar{\epsilon}_G=1.575 \beta$
is the delocalisation energy per atom of an ideal graphitic plane.
\end{document}