Shortcut model

An important question in statistical mechanics is the dependence of model behaviour on the dimension of the system. The shortcut model^[1][2] was introduced in the course of studying this dependence. The model interpolates between discrete regular lattices of integer dimension.

Introduction

The behaviour of different processes on discrete regular lattices have been studied quite extensively. They show a rich diversity of behaviour, including a non-trivial dependence on the dimension of the regular lattice.^{[3][4][5][6][7][8][9][10][11]} In recent years the study has been extended from regular lattices to complex networks. The shortcut model has been used in studying several processes and their dependence on dimension.

Dimension of complex network

Usually, dimension is defined based on the scaling exponent of some property in the appropriate limit. One property one could use ^[2] is the scaling of volume with distance. For regular lattices ${\displaystyle \textstyle \mathbf {Z} ^{d}}$ the number of nodes $\textstyle j$ within a distance ${\displaystyle \textstyle r(i,j)}$ of node $\textstyle i$ scales as ${\displaystyle \textstyle r(i,j)^{d}}$.

For systems which arise in physical problems one usually can identify some physical space relations among the vertices. Nodes which are linked directly will have more influence on each other than nodes which are separated by several links. Thus, one could define the distance ${\displaystyle \textstyle r(i,j)}$ between nodes $\textstyle i$ and $\textstyle j$ as the length of the shortest path connecting the nodes.

For complex networks one can define the volume as the number of nodes $\textstyle j$ within a distance ${\displaystyle \textstyle r(i,j)}$ of node $\textstyle i$, averaged over $\textstyle i$, and the dimension may be defined as the exponent which determines the scaling behaviour of the volume with distance. For a vector ${\displaystyle \textstyle {\vec {n}}=(n_{1},\dots ,n_{d})\in \mathbf {Z} ^{d}}$, where $\textstyle d$ is a positive integer, the Euclidean norm ${\displaystyle \textstyle \|{\vec {n}}\|}$ is defined as the Euclidean distance from the origin to ${\displaystyle \textstyle {\vec {n}}}$, i.e.,

${\displaystyle \|{\vec {n}}\|={\sqrt {n_{1}^{2}+\cdots +n_{d}^{2}}}.}$

However, the definition which generalises to complex networks is the ${\displaystyle \textstyle L^{1}}$ norm,

${\displaystyle \|{\vec {n}}\|_{1}=\|n_{1}\|+\cdots +\|n_{d}\|.}$

The scaling properties hold for both the Euclidean norm and the ${\displaystyle \textstyle L^{1}}$ norm. The scaling relation is

${\displaystyle V(r)=kr^{d},}$

where d is not necessarily an integer for complex networks. $\textstyle k$ is a geometric constant which depends on the complex network. If the scaling relation Eqn. holds, then one can also define the surface area ${\displaystyle \textstyle S(r)}$ as the number of nodes which are exactly at a distance $\textstyle r$ from a given node, and ${\displaystyle \textstyle S(r)}$ scales as

${\displaystyle S(r)=kdr^{d-1}.}$

A definition based on the complex network zeta function^[1] generalises the definition based on the scaling property of the volume with distance^[2] and puts it on a mathematically robust footing.

Shortcut model

The shortcut model starts with a network built on a one-dimensional regular lattice. One then adds edges to create shortcuts that join remote parts of the lattice to one another. The starting network is a one-dimensional lattice of $\textstyle N$ vertices with periodic boundary conditions. Each vertex is joined to its neighbors on either side, which results in a system with $\textstyle N$ edges. The network is extended by taking each node in turn and, with probability $\textstyle p$, adding an edge to a new location $\textstyle m$ nodes distant.

The rewiring process allows the model to interpolate between a one-dimensional regular lattice and a two-dimensional regular lattice. When the rewiring probability ${\displaystyle \textstyle p=0}$, we have a one-dimensional regular lattice of size $\textstyle N$. When ${\displaystyle \textstyle p=1}$, every node is connected to a new location and the graph is essentially a two-dimensional lattice with $\textstyle m$ and ${\displaystyle \textstyle N/m}$ nodes in each direction. For $\textstyle p$ between $\textstyle 0$ and $\textstyle 1$, we have a graph which interpolates between the one and two dimensional regular lattices. The graphs we study are parametrized by

${\displaystyle {\text{size}}=N,\,}$

${\displaystyle {\text{shortcut distance}}=m,\,}$

${\displaystyle {\text{rewiring probability}}=p.\,}$

Application to extensiveness of power law potential

One application using the above definition of dimension was to the extensiveness of statistical mechanics systems with a power law potential where the interaction varies with the distance $\textstyle r$ as ${\displaystyle \textstyle 1/r^{\alpha }}$. In one dimension the system properties like the free energy do not behave extensively when ${\displaystyle \textstyle 0\leq \alpha \leq 1}$, i.e., they increase faster than N as ${\displaystyle \textstyle N\rightarrow \infty }$, where N is the number of spins in the system.

Consider the Ising model with the Hamiltonian (with N spins)

${\displaystyle H=-{\frac {1}{2}}\sum _{i,j}J(r(i,j))s_{i}s_{j}}$

where ${\displaystyle \textstyle s_{i}}$ are the spin variables, ${\displaystyle \textstyle r(i,j)}$ is the distance between node $\textstyle i$ and node $\textstyle j$, and ${\displaystyle \textstyle J(r(i,j))}$ are the couplings between the spins. When the ${\displaystyle \textstyle J(r(i,j))}$ have the behaviour ${\displaystyle \textstyle 1/r^{\alpha }}$, we have the power law potential. For a general complex network the condition on the exponent $\textstyle \alpha$ which preserves extensivity of the Hamiltonian was studied. At zero temperature, the energy per spin is proportional to

${\displaystyle \rho =\sum _{i,j}J(r(i,j)),}$

and hence extensivity requires that ${\displaystyle \textstyle \rho }$ be finite. For a general complex network ${\displaystyle \textstyle \rho }$ is proportional to the Riemann zeta function ${\displaystyle \textstyle \zeta (\alpha -d+1)}$. Thus, for the potential to be extensive, one requires

${\displaystyle \alpha >d.\,}$

Other processes which have been studied are self-avoiding random walks, and the scaling of the mean path length with the network size. These studies lead to the interesting result that the dimension transitions sharply as the shortcut probability increases from zero.^[12] The sharp transition in the dimension has been explained in terms of the combinatorially large number of available paths for points separated by distances large compared to 1.^[13]

Conclusion

The shortcut model is useful for studying the dimension dependence of different processes. The processes studied include the behaviour of the power law potential as a function of the dimension, the behaviour of self-avoiding random walks, and the scaling of the mean path length. It may be useful to compare the shortcut model with the small-world network, since the definitions have a lot of similarity. In the small-world network also one starts with a regular lattice and adds shortcuts with probability $\textstyle p$. However, the shortcuts are not constrained to connect to a node a fixed distance ahead. Instead, the other end of the shortcut can connect to any randomly chosen node. As a result, the small world model tends to a random graph rather than a two-dimensional graph as the shortcut probability is increased.

References

1. O. Shanker (2007). "Graph Zeta Function and Dimension of Complex Network". Modern Physics Letters B. 21 (11): 639–644. Bibcode:2007MPLB...21..639S. doi:10.1142/S0217984907013146.
2. O. Shanker (2007). "Defining Dimension of a Complex Network". Modern Physics Letters B. 21 (6): 321–326. Bibcode:2007MPLB...21..321S. doi:10.1142/S0217984907012773.
3. O. Shanker (2006). "Long range 1-d potential at border of thermodynamic limit". Modern Physics Letters B. 20 (11): 649. Bibcode:2006MPLB...20..649S. doi:10.1142/S0217984906011128.
4. D. Ruelle (1968). "Statistical mechanics of a one-dimensional lattice gas". Communications in Mathematical Physics. 9 (4): 267. Bibcode:1968CMaPh...9..267R. doi:10.1007/BF01654281.
5. F. Dyson (1969). "Existence of a phase-transition in a one-dimensional Ising ferromagnet". Communications in Mathematical Physics. 12 (2): 91. Bibcode:1969CMaPh..12...91D. doi:10.1007/BF01645907.
6. J. Frohlich & T. Spencer (1982). "The phase transition in the one-dimensional Ising Model with 1/r² interaction energy". Communications in Mathematical Physics. 84: 87. Bibcode:1982CMaPh..84...87F. doi:10.1007/BF01208373.
7. M. Aizenman; J.T. Chayes; L. Chayes; C.M. Newman (1988). "Discontinuity of the magnetization in one-dimensional 1/|x−y|² Ising and Potts models". Journal of Statistical Physics. 50: 1. Bibcode:1988JSP....50....1A. doi:10.1007/BF01022985.
8. J.Z. Imbrie; C.M. Newman (1988). "An intermediate phase with slow decay of correlations in one dimensional 1/|x−y|² percolation, Ising and Potts models". Communications in Mathematical Physics. 118 (2): 303. Bibcode:1988CMaPh.118..303I. doi:10.1007/BF01218582.
9. E. Luijten & H.W.J. Blöte (1995). "Monte Carlo method for spin models with long-range interactions". International Journal of Modern Physics C. 6 (3): 359. Bibcode:1995IJMPC...6..359L. doi:10.1142/S0129183195000265.
10. R.H. Swendson & J.-S. Wang (1987). "Nonuniversal critical dynamics in Monte Carlo simulations". Physical Review Letters. 58 (2): 86. Bibcode:1987PhRvL..58...86S. doi:10.1103/PhysRevLett.58.86.
11. U. Wolff (1989). "Collective Monte Carlo Updating for Spin Systems". Physical Review Letters. 62 (4): 361–364. Bibcode:1989PhRvL..62..361W. doi:10.1103/PhysRevLett.62.361. PMID 10040213.
12. O. Shanker (2008). "Algorithms for Fractal Dimension Calculation". Modern Physics Letters B. 22 (07): 459–466. Bibcode:2008MPLB...22..459S. doi:10.1142/S0217984908015048.
13. O. Shanker (2008). "Sharp dimension transition in a shortcut model". J. Phys. A. 41 (28): 285001. Bibcode:2008JPhA...41B5001S. doi:10.1088/1751-8113/41/28/285001.