The World's Colonisation and Trade Routes
Formation as Imitated by Slime Mould
Andrew Adamatzky
University of the West of England, Bristol, United Kingdom
This is unedited preprint with low-resolution photographs.
Final and edited version of this paper is published in
Int. J. Bifurcation Chaos, 22, 1230028 (2012) [26 pages]
DOI: 10.1142/S0218127412300285
Abstract
The Plasmodium of Physarum polycephalum is renowned for spanning sources of nutrients with
networks of protoplasmic tubes. The networks transport nutrients and metabolites across the Plas-
modium's body. To imitate a hypothetical colonisation of the world and formation of major trans-
portation routes we cut continents from agar plates arranged in Petri dishes or on the surface of
a three-dimensional globe, represent positions of selected metropolitan areas with oat flakes and
inoculate the Plasmodium in one of the metropolitan areas. The Plasmodium propagates towards
the sources of nutrients, spans them with its network of protoplasmic tubes and even crosses bare
substrate between the continents. From the laboratory experiments we derive weighted Physarum
graphs, analyse their structure, compare them with the basic proximity graphs and generalised graphs
derived from the Silk Road and the Asia Highway networks.
Keywords: biological transport networks, unconventional computing, slime mould
1 Introduction
Nature-inspired computing paradigms and experimental laboratory prototypes are demonstrated reason-
able success in approximation of shortest, and often collision-free, paths between two given points in an
Euclidean space or a graph. Examples include ant-based optimisation of communication networks |15) .
approximation of a shortest path in experimental reaction-diffusion chemical systems [T], gas-discharge
analog systems [35] , spatially extended crystallisation systems [5] , fungi mycelia networks (25] , and maze
solving by Physarum polycephalum |29) . Amongst all explored so far experimental prototypes of path-
computing devices slime mould P. polycephalum is the most user-friendly and easiest to cultivate and
observe biological substrate. In its plasmodium stage, P. polycephalum spans scattered sources of nutri-
ents with a network of protoplasmic tubes. Nakagaki and colleagues demonstrated that the plasmodium
optimises its protoplasmic network to cover all sources of nutrients and to assure quick and fault-tolerant
distribution of nutrients in the Plasmodium's body [281 [23 HOI ELL] • There are experimental evidences [2]
that in the course of the Plasmodium's foraging behaviour the protoplasmic network follows the Tous-
1
Adamatzky A.
Int. J. Bifurcation Chaos, 22, 1230028 (2012) [26 pages]
saint [43, 25, 21 hierarchy of proximity graphs however neither of the known proximity graphs is exactly
matched by the protoplasmic networks [7].
The slime mould's protoplasmic network is responsible for transportation of nutrients and metabolites,
and intra-cellular communication. Therefore it is a matter of the natural scientific curiosity to compare the
protoplasmic transport networks with man-made transport network. Two first comparisons were made
by Adamatzky and Jones [3] — the slime mould imitation of the UK motorways, and Tero and colleagues
[4"0] — the slime mould imitation of the Tokyo rail roads. These pioneering results were followed by a
series of laboratory experiments on evaluation of motorway/highway networks in Mexico 4 , Spain and
Portugal [BJ, Brazil [S], Australia 9 and Canada [10] . These papers demonstrated that P. polycephalum
approximates existing road networks up to some degree, especially with regards to the core spanning
structure, but the slime mould also allows for a wide range of structural deviations in the transport
networks and solutions are not always optimal or predictable. The shape of any particular country and
the spatial configuration of major urban areas, represented by sources of nutrients, might substantially
determine the resultant topology of the protoplasmic networks developed. In certain countries, the
road network designs were influenced not only by logic but also by political, economic and somewhat
idealogical decisions. Therefore, more experiments are required to provide a viable generalisation, to
undertake a comparative analyses of slime mould transport networks constructed in different countries
and to built a theory of P. polycephaum based road planning and associated urban development. In
present paper we decided to undertake an experimental laboratory exercise on playing a hypothetical
scenario of the colonisation of the World by large-scale amorphous substrate and formation of the world-
wide transportation routes crossing countries and linking continents. The resultant protoplasmic networks
were compared, at the abstract level of generalised graphs, with ancient road network — the Silk Road,
and future road network — the Asian Highways.
2 Materials and methods
The plasmodium of P. polycephalum is cultivated in a plastic container, on paper kitchen towels moistened
with still water, and fed with oat flakes. For experiments we use 220 x 220 mm polystyrene square Petri
dishes and 2% agar gel (Select agar, by Sigma Aldrich) as a substrate and also the three-dimensional
globe 15 cm diameter (Stellanova, Germany). When experiments are conducted in the Petri dishes the
hot agar gel is poured in the dishes to fill the dishes by 2-3 mm. When experimenting with the globe we
poured hot gel onto the globe while rotating the globe, so it is covered by agar uniformly. The globe was
covered by the agar gel layer by layer. When one layer cooled down and settled, next layer of hot get
applied. When agar cooled down, the shapes corresponding to oceans and seas are cut out and removed,
and only the dry land is represented by the agar substrate.
We considered 24 metropolitan areas as follows (Fig. Ilk):
1.
Beijing
G.
Ho Chi Minn
11.
Karachi
16.
Lagos
21.
Lima
2.
Seoul
7.
Jakarta
12.
Tehran
17.
Kinshasa
22.
Sao Paulo
3.
Tokyo
8.
Kolkata
13.
Moscow
18.
New York
23.
Canberra
4.
Hong Kong
9.
Mumbai
14.
Istanbul
19.
Mexico City
5.
Ha Noi
10.
Delhi
15.
London
20.
Bogota
24.
Wellington
The elements of U were selected based on their population size, proximity to other urban/metropolitan
2
Adamatzky A.
Int. J. Bifurcation Chaos, 22, 1230028 (2012) [26 pages]
Adamatzky A.
Int. J. Bifurcation Chaos, 22, 1230028 (2012) [26 pages]
Adamatzky A.
Int. J. Bifurcation Chaos, 22, 1230028 (2012) [26 pages]
(a) (b)
Figure 3: Plasmodium crosses bare space between the continent-shaped agar plates.
areas, representativeness of a continent/country. The following area were not chosen for experiments,
despite being amongst largest metropolitan areas: Manila is not included due to its proximity to Jakarata;
Shanghai due to its proximity to Beijing; Buenos Aires due to its proximity to San Paolo; Osaka-Kobe-
Kyoto due their proximity to Tokyo; Dhaka due to its proximity to Kalkutta. The capital cities Canberra
and Wellington are by no means in the list of largest cities or metropolitan areas, however they were
included to give slime mould a change to move to these countries. To represent areas of U we placed oat
flakes (each flake weights 9-13 mg and is 5-7 mm in diameter) in the positions of agar plate corresponding
to the areas (Figs. [T]d and[2|.
At the beginning of each experiment an oat flake colonised by plasmodium (25-30 mg plasmodial
weight) is placed in Beijing area. We have chosen Beijing as a starting point for the slime mould because
Beijing is one of the most populated cities in the world [13] and also because the Chinese civilisation is
amongst earliest civilisations [2"r?l 126) .
We undertook 38 experiments in total: eight experiments with the globe and 30 experiments with the
Petri dish. The Petri dishes and the globe with plasmodium were kept in darkness (the globe was placed
in a glass container to prevent drying of agar gel), at temperature 22-25°C, except for observation and
image recording. Periodically, the dishes were scanned with an Epson Perfection 4490 scanner and the
globe was photographed with FujiFilm FinePix. As exemplified in Fig. [3] absence of a humid growing
substrate prevented plasmodium from spreading 'uncontrollably' into parts of substrate corresponding to
oceans and seas yet allowed the plasmodium to migrate between the continents when necessary.
To generalise our experimental results we constructed a Physarum graph with weighted-edges. A
Physarum graph is a tuple P = (U,E, w), where U is a set of urban areas, E is a set edges, and
w : E — > [0, 1] associates each edge of E with a probability (or weights). For every two regions a and b
from U there is an edge connecting a and b if a Plasmodium's protoplasmic link is recorded at least in
one of k experiments, and the edge (a, b) has a probability calculated as a ratio of experiments where
protoplasmic link (a, b) occurred in the total number of experiments k. We do not take into account the
exact configuration of the protoplasmic tubes but merely their existence. Further we will be dealing with
threshold Physarum graphs P(9) = (U, T(E), io, 0). The threshold Physarum graph is obtained from
Physarum graph by the transformation: T(E) = {e G E : w(e) > &}. That is all edges with weights less
than 6 are removed.
5
Adamatzky A.
Int. J. Bifurcation Chaos, 22, 1230028 (2012) [26 pages]
3 Results
3.1 Scenarios of colonisation
Scenarios of the plasmodium development in Petri dishes and on the globe are strikingly similar. Here
we consider one example of colonisation of U in the Petri dish and three examples of the colonisation
on the globe. In first day after inoculation in the Petri dish, the plasmodium propagates from Beijing
to Seoul and Tokyo in the east; from Beijing to Hong Kong, Ha Noi and Ho Chi Minh in the south;
and from Beijing to Kolkata, Mumbai, Delhi and Karachi in the west (Fig. [4^,). On second day the
plasmodium explores the space north of Beijing, links Karachi and Tehran, and Tehran and Istanbul
with protoplasmic tubes. It also growth from Istanbul to Moscow and from Moscow to London, and
from Tehran to Lagos and Kinshasa (Fig. [4]b) . The plasmodium grows from London to Iceland, then
to Greenland. It propagates from Greenland to Canada and eventually reaches New York and Mexico
City on the fifth day of the experiment (Fig. [4ji). The final stage of spanning U takes place on the sixth
day of the experiment: the plasmodium propagates from Mexico City to Bogota and from Bogota to
Lima and Sao Paulo. At the same time the plasmodium grows from Ho Chi Minh to Jakarta and from
Jakarta to Canberra and Wellington (Fig. [4ji) . As soon as all sources of nutrients, representing the areas
of U, are spanned by the plasmodium, the plasmodium remains in its configuration for a couple of days
(Fig. |4j;) . By that time the nutrients become depleted, and the substrate contaminated with products of
metabolism and the plasmodium tries to migrate to other areas and/or form sclerotium.
Two examples shown in Fig. [5j demonstrate that colonisation of East and South Asia can go by dif-
ferent scenarios however propagations towards Western Europe and Americas have matching trajectories.
Let us consider scenario in Fig. [5^,. In the second day after inoculation plasmodium propagates from
Beijing to Seoul and then to Tokyo, and from Beijing to Delhi and to Kolkata, and from Delhi to Karachi
and to Mumbai. On the third day, plasmodium propagates from Karachi to Tehran to Istanbul to Moscow
(Fig. |6(a)| ) . It then continuous to explore space north-east of Moscow. On the fourth day plasmodium
develops protoplasmic links from Delhi to Kolkata, from Kolkata to Ha Noi and Ho Chi Minh, and from
Ho Chi Minh to Jakarta; and, completes the route from Moscow to Beijing. It also explores space east
of Jakarta (Fig. 6(b)). On the fifth, sixth and seventh days slime mould propagates from Jakarta to
Canberra (Fig. 6(c)[ ), and from Istanbul to London. The plasmodium reaches New York from London via
Iceland, Greenland and Canada on ninth day after inoculation in Beijing. Then it propagates from New
York to Mexico City, Bogota, Lima and Sao Paulo (Fig. [5^) .
Early stages of colonisation (Fig. [5^,) show the following priorities of the Plasmodium's development:
east then west then south. The scenario with early development: east then south then west, is shown in
Fig. [5]3 and illustrated in Figs. [7(a)| ff(cj\ [7(b)] |7(d")j
In first two days after inoculation in Beijing the plasmodium propagates to Seoul and Tokyo in the
east, to Delhi in the south-west and Hong Kong in the south-east. It also ventures north to Siberia. On
the third day the plasmodium grows from Ha Noi to Kolkata and Ho Chi Minh, from Ho Chi Minh to
Jakarta and Canberra (Fig. 7(a) ). On the forth day the plasmodium propagates from Delhi to Karachi,
Tehran, Istanbul, Moscow, London, and ventures from London to Iceland (Fig. 7(c)). Protoplasmic
links between London and Lagos and Kinshasa are developed by the slime mould at the fifth day of the
experiment. On the sixth and seventh days plasmodium crosses Greenland and Canada towards USA. It
firstly reaches New York and then propagates towards Mexico City, from where it spans Bogota, Lima
and Sao Paul o (Fig . [7(b)T ). This scenarios also displays a pronounced transport link between Beijing and
London (Fig. [7(d)).
Example of incomplete spanning is shown in Fig. [8] Transport links between Beijing, Seoul and Tokyo,
and Beijing and Hong Kong are established in two days after inoculation (Fig. |9(d)[ ). On the third day
the plasmodium connects the following areas by the chain of protoplasmic tubes Hong Kong- Ha Noi-
6
Adamatzky A. Int. J. Bifurcation Chaos, 22, 1230028 (2012) [26 pages]
(e) 5 days
Figure 4: Example of plasmodium development on a configuration continent-shaped agar plates.
7
Adamatzky A. Int. J. Bifurcation Chaos, 22, 1230028 (2012) [26 pages]
2 nd day 3 ,d day - -* 4* 1 day - -*- 5 th ' - 7" days ■ ■ ■ ► 9 ,h day
10 th day
Mi
■
YV^"",-' V *** 1 ' . -AusujliJ \
(a)
2" d day »*■ 3 ,d day -
lively t 3
. . • v
>■ 5 ,h day ■ ■ ■► 6 ,h day ■♦ 7'" day
(b)
Figure 5: Two scenarios of plasmodial active zones propagation on the globe.
8
Adamatzky A.
Int. J. Bifurcation Chaos, 22, 1230028 (2012) [26 pages]
Figure 6: Experimental examples of scenario Fig. [5p. (a) Propagation of slime mould from Karachi to
Tehran to Istanbul to Moscow on the third day after inoculation, (b) On the fourth day of experiment
Plasmodium links Ha Noi to Ho Chi Minh to Jakarta with its protoplasmic tubes and explores space east
of Jakarta, (c) Plasmodium's active zone enters Australia at the port Derby.
9
Adamatzky A. Int. J. Bifurcation Chaos, 22, 1230028 (2012) [26 pages]
(c)
(d)
Figure 7: Illustrations of the scenario Fig. |5p. |(a)| Plasmodium propagates from South-East Asia to
Australia, (c) Plasmodium propagates west and north-west. |(b)| Plasmodium reaches USA via Iceland,
Greenland and Canada and then spans Latin America. |(d)| Transport link between Beijing and London.
10
Adamatzky A.
Int. J. Bifurcation Chaos, 22, 1230028 (2012) [26 pages]
2 nd day
»• 3 rd day - - • 4 ,h day
>. 5 th day
■ ► 6 th day
re? r v
/ ' FlnLii-dV^/-'.
© ©
Australia
Figure 8: Examples of incomplete spanning of urban areas. Trajectories of plasmodial active zone prop-
agation in ten days of experiment.
Kolkata- Delhi- Karachi- Tehran- Istanbul (Fig. |9(b)[ ). Moscow and Lagos and Kinshasa are spanned
by the protoplasmic network on the fourth day of colonisation (Fig. 9(c) ). It takes slime mould one more
day to cross South Atlantic to make the link between Lagos and Sao Paulo (Fig. 9(a) ). Protoplasmic links
(Sao Paulo, Lima), (Sao Paulo, Bogota), (Bogota, Mexico City), (Mexico City, New York) are grown on
the sixth day of the experiment (Fig. [8}) .
Generalised Physarum graphs derived from experiments in Petri dishes (let us call them P-graphs) and
the globe (G-graphs) are shown in Fig. [ToJ P-graphs represent G-graphs sufficiently. Only three edges
(albeit represented in one experiment each) of P-graphs are not represented by G-graphs: (Kinshasa,
Sao Paulo), (Moscow, Istanbul) and (Beijing, Hong Kong). In contrast, eight edges of P-graphs are not
represented by G-graphs: (Hong Kong, Jakarta), (Ha Noi, Ho Chi Minh), (Jakarta, Kolkata), (Tehran,
Lagos), (Istanbul, Kinshasa), (London, Mexico City), (Mexico City, Bogota) and (Lagos, Sao Paulo).
Further in the paper we consider generalised Physarum graphs as a union of P- and G-graphs.
Generalised Physarum graphs for critical values of
10
Generalised Physarum
graphs P(9) are sub-graphs of P- and G-graphs for 9 > J|. Physarum graphs are planar for 8 >
however with acquiring planarity the Physarum graph P(^) loses connectivity and Wellington becomes
isolated (Fig. flOp). The highest value of 8 for which graph P(8) — {Wellington } remains connected
are shown in Fig.
14
is || (Fig. |10[:). When 8 increases to || Americas become disconnected from Eurasia and Africa
Ame rica n transport pathways become a single chain New York- Mexico City- Bogota-
and
Lima- Sao Paulo
(Fig. 10(1). Further increase of 8 to || leads to isolation of Sao Paulo and Canberra (Fig.
For 8
22
:S8
10
0-
all Americas areas of U become isolated and the transport network at this continent
virtually disappears. The graph P(||) has two components of connectivity:
the chain Beijing- Seoul- Tokyo
11
Adamatzky A.
Int. J. Bifurcation Chaos, 22, 1230028 (2012) [26 pages]
"1
Figure 9: Illustration of the scenario Fig. |8] [(a)] Plasmod ium crosses South Atlantic from West Africa
to South America, (c) Protoplasmic tube connecting Istanbul to Lagos develops on the fourth day of
experiment. |(d)| Transport links Beijing, Seoul and Tokyo, and Beijing and Hong Kong are established in
two days after inoculation. |(b)| Chain Hong Kong- Ha Noi- Kolkata- Delhi- Karachi- Tehran- Istanbul
is developed on the third day of experiment.
12
Adamatzky A.
Int. J. Bifurcation Chaos, 22, 1230028 (2012) [26 pages]
(a) (b)
Figure 10: Physarum graphs obtained from experiments in (a) the Petri dishes, P-graph, and (b) the
globe, G-graph. Thickness of a line is proportional to the line's weight.
• the component consisting of the chain Jakarta- Ho Chi Minh- Kolkata- Delhi- Tehran- Istanbul-
London- Lagos- Kinshasa with two branches Ha Noi- Hong Kong and Delhi- Mumbai attached
and the cycle Tehran- Moscow- Istanbul- Tehran (Fig. 10').
Generalised Physarum graph is transformed to the acyclic graph at 9 = || when the edge (Tehran,
Moscow) becomes cut off (Fig. 10
3.2 Proximity graphs
A planar graph consists of nodes which are points of the Euclidean plane and edges which are straight
segments connecting the points. A planar proximity graph is a planar graph where two points are
connected by an edge if they are close in some sense. A pair of points is assigned a certain neighbourhood,
and points of the pair are connected by an edge if their neighbourhood is empty. Here we consider the
most common proximity graphs as follows.
MST: The Euclidean minimum spanning tree [32] is a connected acyclic graph which has minimum
possible sum of edges' lengths (Fig. [12^,).
RNG: Points a and b are connected by an edge in the Relative Neighbourhood Graph if no other
point c is closer to a and b than dist{a, b) j3Hj (Fig- [l2p).
GG: Points a and b are connected by an edge in the Gabriel Graph if disc with diameter dist{a, b)
centred in middle of the segment ab is empty [Ml [25] (Fig. 12 x).
In general, the graphs relate as MST C RNG C GG 43, 25, 2l|; this is called Toussaint hierarchy.
Why do we need to compare Physarum graphs with the proximity graphs MST, RNG and GG? The
minimum spanning tree helps to evaluate optimality of the protoplasmic networks: minimal distances of
nutrient transportation yet complete covering of the sources of nutrients (sites of U) . Being an acyclic
graph the spanning tree is sensitive to a structural damage. Removal of a single edge might transform
the spanning tree to two disconnected trees. The relative neighbourhood graph and the Gabriel graph
show higher degree of fault-tolerance (depending on exact configuration of nodes). It also considered
to be optimal in terms of total edge length and travel distance. The graphs RNG and GG are used
in geographical variational analysis [TBI [25] , simulation of epidemics [12] ; and design of ad hoc wireless
13
Adamatzky A. Int. J. Bifurcation Chaos, 22, 1230028 (2012) [26 pages]
Figure 11: Generalizei^Physarum graphs P(#) for selected values of 9.
Adamatzky A.
Int. J. Bifurcation Chaos, 22, 1230028 (2012) [26 pages]
Figure 12: Relative neighbourhood graph (a) and Gabriel graph (b) constructed on urban areas of U.
15
Adamatzky A.
Int. J. Bifurcation Chaos, 22, 1230028 (2012) [26 pages]
Figure 13: Intersections of generalised Physarum graphs (abc) P(^) and (def) P(§§) with proximity
graphs: (ad) minimum spanning tree MST, (be) relative neighbourhood graph RNG, (cf) Gabriel graph
GG.
16
Adamatzky A.
Int. J. Bifurcation Chaos, 22, 1230028 (2012) [26 pages]
Figure 14: Graphs of the Silk Road and the Asian Highways, (a) the Silk Road graph, comprised of the
Silk Road (thick line) and associated trade routes (thin lines); the graph is derived from the scheme of
Trans- Asian trade 500 BC - AD 750 [35] • (b) the Asian Highway network graph; the graph is derived
from the Asian Highway Network Map [34] .
networks [33] E3 EES [HI SS]- The proximity graphs, especially RNG, are invaluable in simulation
of human-made, road networks; these graphs are validated in studies of Tsukuba central district road
networks 46, 47J. Gabriel graphs, particularly their relaxed versions [13] are an ideal tool for path finding
and online routing on planar graphs [12] .
Intersections of generalised Physarum graphs with the proximity graphs is shown in (Fig. 13).
Finding 1. The minimum spanning tree and the relative neighbourhood graph are sub- graphs o/P(gg).
Gabriel graph would be a sub-graph o/P(^) if the Physarum graph had the edges (Hong Kong, Seoul)
and (New York, Bogota).
The MST differs from RNG only by absence of the single edge (Mexico City, Bogota), the rest
is the same and is matched by edges of the generalised Physarum graph P(gg) (Fig. 13 lb). The fact
that GG C P(l|)U(Hong Kong, Seoul) U (Mexico City, Bogota) (Fig. 13p) plays against the geographical
validity of the Gabriel graph. The only physically possible routes from New York to Bogota and from
Hong Kong to Seoul are maritime routes; overland route from New York to Bogota is passing via Mexico
City and from Hong Kong to Seoul via Beijing.
Finding 2. The
Kolkata- Delhi- Karachi
longest chain in MST fl P(§§) * s C\ — (Canberra
Tehran- Istanbul- London- New York-
Jakarta- Ho
Mexico City)
Chi Minh- Ha Noi-
Two fragments of the chain C\\ Jakarta- Ho Chi Minh- Ha Noi and Kolkata to London (Fig. 131)
match the Silk Road, see Sect. 3.3 The chain C\ is elongated by segment (Bogota- Lima- Sao Paulo) in
RNGnP(if) (Fig. Ill;) and GGnP(Af) (Fig.pt).
3.3 Protoplasmic Networks, the Silk Road and the Asian Highways
To evaluate how well plasmodial networks match large-scale transportation systems we compare gener-
alised Physarum graphs with graphs derived from the Silk Road and the Asia Highway network.
17
Adamatzky A.
Int. J. Bifurcation Chaos, 22, 1230028 (2012) [26 pages]
The Silk Road is a long-distance trade route across Asia developed in the first millennium to transport
goods between China and Western Asia and Northern Europe [UJ . Commonly the Silk Road is comprised
of the main overland route from Luoyang to Seleucia/Clesiphon and associated overland and maritime
routes reaching as far Japan in the east and as far as Colchester and Holborough in the west [351 HI] •
A graph of the Silk Road is shown in Fig. [14^ ,. When constructing the Silk road graph we tracked
geographical positions of modern cities along the historical Silk Road and associated trade routes, with
minor assumptions. Thus the Silk Road is represented by edge (Beijing, Tehran). The road actually
started at Luoyang and passed via Hecatompylos and ended in Ctesiphon. We assumed Hecatompylos
was in proximity of Tehran and Beijing is in proximity of Luoyang. Similarly, we assume that edge
Delhi to Kolkata corresponds to the trade route connecting Mathura (proximity of Delhi) and Tampluk
via Pataliputra (proximity of Kolkata). Link (Tokyo, Hong Kong) corresponds to a maritime route
from south of Japan to Oc Eo via Guangzhou and passing between the continental China and Hong
Kong island. The edge (Tehran, Istanbul) represents land route from Hecatompylos to Dura Europos
and Antioch combined with maritime trade route from Antioch to Constantinple. The edge (Istanbul,
London) represents maritime route Constantinople - Athens - Rome - Massilia - Colchester.
The Asian Highways is a network of roads selected for regional transport cooperation "initiative aimed
at enhancing efficiency and development of the road transport infrastructure in Asia" |34j . The Asian
Highway network consists of 141,000 km of roads running across 32 member States. When constructing
Asian Highway graph (Fig. 14 d) we adopted the following match between the graph's edges and the
highways (see map of the Asian Highway network in |34j ) :
• (Beijing, Seoul): potential route AH1 Beijing to Shenyang and existing route AH1 Shenyang to
Seoul:
• (Beijing, Hong Kong): route AH1 Beijing to Hong Kong via Zhengzhou, Xinyang, Xianglan;
• (Beijing, Delhi): route Al Beijing to Zhengzhou, potential routes AH34 Zhengzhou to Xi'an, AH5
Xi'an to Langzhou, AH42 Langzhou to Lhasa, and final existing routes AH42 Lhasa to Zhangmu
and AH2 Zhangmu to Delhi;
• (Beijing, Tehran): comprised of segments of the following highways (ordered from Beijing to Tehran)
AH1, AH34, AH5, AH4, AH65, AH62, AH76, AH1;
• (Beijing, Moscow): AH3 Beijing to Ulan Ude, AH5 from Ulan Ude to Moscow via Irkutsk, Novosi-
birsk, Omsk, Petropavlosk, Chelyabinsk, Samar;
• (Seoul, Tokyo): maritime route Pusan to Fukuoka and overland route AH1 Fukuoka to Tokyo;
• (Hong Kong, Ha Noi): potential route AH1 Guangzhou to Nanning, AH1 Nanning to Hanoi;
• (Ha Noi, Ho Chi Minh): AH1;
• (Ha Noi, Kolkata): AH14 Hanoi - Kummig - Mandalay and AH1 Mandalay - Dispur - Kolkata;
• (Ho Chi Minh, Jakarta): AH1 Ho Chi Minh to Kabin Bun, AH19 Kabin Bun to Bangkok, AH2
Bangkok to Hat Yai, AH18 Hat Yai to Singapore, and maritime Singapore to Jakarta;
• (Kolkata, Mumbai): AH46 and AH47;
• (Kolkata, Delhi): AH1;
• (Delhi, Karachi): AH1 Delhi to Lahore, AH2 Lahore to Rohn, AH4 Rohn to Karachi;
18
Adamatzky A.
Int. J. Bifurcation Chaos, 22, 1230028 (2012) [26 pages]
• (Karachi, Tehran): AH7 Karachi to Kandahar and AH1 Kandahar - Dilaram - Herat - Zhabzevar
- Tehran;
• (Tehran, Moscow): AH8 Tehran - Baku - Astahan - Volgograd - Moscow;
• (Tehran, Istanbul): AH1 Tehran - Yerevan - Ankara and AH5 Ankara to Istanbul.
The Asian Highway graph has the same number of edges as the Silk Road graph however they are
not isomorphic (Fig. 14 1.
Finding 3. Slime mould P. polycephalum approximates over 76% of the Silk Road routes and the Asian
Highway routes.
Physarum graph P(<§j) (Fig. nib) approximates 13 of 17 edges of the Silk Road graph(Fig. 14 1) and
also 13 of 17 edges of the Asian Highway graph (Fig.[l4p). Physarum graph P(g§) (Fig- H : ) approximates
11 of 17 edges of the Silk Road graph and also 11 oFl7 edges of the Asian Highway graph.
Finding 4. Transport links (Beijing, Tehran) and (Beijing, Hong Kong) of the Silk Road and the Asian
Highway network are never approximated by P. polycephalum.
The following routes of the Silk Road are never approximated by any Physarum graph: (Delhi,
Tehran), (Beijing, Hong Kong), (Tokyo, Hong Kong), (Hong Kong, Jakarta). Routes (Beijing, Tehran)
and (Beijing, Kolkata) are approximated by P(J|) but not P(§§).
The following routes of the Asian highways are not approximated by any Physarum graph: (Beijing,
Tehran), (Beijing, Hong Kong), (Beijing, Delhi) and (Jakarta, Kolkata). Routes (Beijing, Moscow) and
(Ho Chi Minh, Kolkata) are approximated by P(gg) but not P(||).
4 Conclusions
To imitate a hypothetical scenario of the World's colonisation and emergence of principal transport
networks we represented the continents with geometrically-shaped plates of non-nutrient agar and the
major metropolitan areas with sources of nutrients and inoculated plasmodium of P. polycephalum in
Beijing. We analysed scenarios of the plasmodium propagation and colonisation of the metropolitan
areas. We derived weighted generalised Physarum graphs from the protoplasmic networks recorded in
the laboratory experiments.
We found that the Physarum graphs include basic proximity graphs — the minimum spanning tree,
the relative neighbourhood graph and the Gabriel graph constructed on the metropolitan areas — as
their sub-graphs. The longest chain of transport links presented in the minimum spanning tree and in
the Physarum graph with edge weights exceeding 0.36 is Canberra- Jakarta- Ho Chi Minh- Ha Noi-
Kolkata- Delhi- Karachi- Tehran- Istanbul- London- New York- Mexico City. We found that slime
mould P. polycephalum approximates over 76% of the Silk Road routes and the Asian Highway network.
Transport links (Beijing, Tehran) and (Beijing, Hong Kong) of the Silk Road and the Asian Highway
network are never approximated by P. polycephalum.
We believe our experimental results will inspire further thoughts, paradigms and approaches for
re-evaluation of historical findings on the emergence of ancient roads and will help to design future trans-
continental pathways. The results could be also applied in analysis of pandemics' dynamics |111 I18|
and in shaping unorthodox approaches to prediction of international military conflicts and battlefield
operations [19].
19
Adamatzky A.
Int. J. Bifurcation Chaos, 22, 1230028 (2012) [26 pages]
References
[11 Adamatzky A., De Lacy Costello B., Asai T. Reaction-Diffusion Computers, Elsevier, Amsterdam,
2005.
Adamatzky A. Developing proximity graphs by Physarum Polycephalum: Does the plasmodium
follow Toussaint hierarchy? Parallel Processing Letters 19 (2008) 105-127.
Adamatzky A. and Jones J. Road planning with slime mould: If Physarum built motorways it
would route M6/M74 through Newcastle [arXiv:0912.3967 vl [nlin.PS] (Dec 2009). Int J Bifurcaton
and Chaos 20 (2010) 3065-3084.
Adamatzky A., Martinez G. J., Chapa-Vergara S. V., Asomoza-Palacio R., Stephens C. R. Approx-
imating Mexican highways with slime mould. Natural Computing 10 (2010) 1195-1214.
Adamatzky A. Hot ice computer. Physics Lett A 347 (2009) 264-271.
Adamatzky A. and Alonso-Sanz R. Rebuilding Iberian motorways with slime mould. Biosystems 105
(2010) 89-100.
Adamatzky A. Physarum Machines: Making Computers from Slime Mould (World Scientific, 2010).
Adamatzky A. and de Oliveira P.P.B. Brazilian highways from slime mold's point of view. Kybernetes
11 (2011).
Adamatzky A. and Prokopenko M. Slime mould evaluation of Australian motorways, Int J Parallel,
Emergent and Distributed Systems (2011) DOL10.1080/17445760.2011. 616204
Adamatzky A. and Akl S.G. Trans-Canada Slimeways: Slime mould imitates the Canadian transport
network (2011) |arXTv: 1105.5084^ 1 [nlin.PS]
Beck E.J. Mays N., Whiteside AW. (Eds.) The HIV Pandemic: Local and Global Implications.
(Oxford University Press, 2008).
Bose P. and Morin P. Competitive online routing in geometric graphs. Theoretical Computer Science
324 (2004) 273-288.
Bose P., Cardinal J., Collette S., Demaine E. D., Polop B., Taslakian P., Zehk N. Relaxed Gabriel
Graphs. In: 21st Canadian Conference on Computational Geometry 2009, Vancouver, BC, August
1719, 2009.
Beijing Municipal Bureau of Statistics (2011). pittp: //www. bj stats .gov, cn/esite/
Dorigo M. and Stutzle T. Ant Colony Optimization MIT Press, 2004.
Gabriel K. R. and R. R. Sokal. A new statistical approach to geographic variation analysis. Systematic
Zoology, 18 (1969) 259-278.
Gernet J. A History of Chinese Civilization. (Cambridge University Press, 1996).
Haour-Knipe M. and Rector R. (Editor) Crossing Borders: Migration, Ethnicity and AIDS (Social
Aspects of AIDS) (Taylor & Francis, 1996.)
[19] Ilachinski A. Artificial War: Multiagent-Based Simulation of Combat (World Scientific, 2004).
20
Adamatzky A. Int. J. Bifurcation Chaos, 22, 1230028 (2012) [26 pages]
[20] Jaromczyk J.W. and Kowaluk M. A note on relative neighbourhood graphs Proc. 3rd Ann. Symp.
Computational Geometry, 1987, 233-241.
[21] Jaromczyk J. W. and G. T. Toussaint, Relative neighborhood graphs and their relatives. Proc. IEEE
80 (1992) 1502-1517.
[22] Jarrett T. C, Ashton D. J., Fricker M., Johnson N. F. Interplay between function and structure in
complex networks Phys. Rev. E 74 (2006) , 026116.
[23] Li X.-Y. Application of computation geometry in wireless networks. In: Cheng X., Huang X., Du D.-
Z. (Eds.) Ad Hoc Wireless Networking (Kluwer Academic Publishers, 2004) 197-264.
[24] Liu X.s The Silk Road in World History (Oxford University Press, 2010).
[25] Matula D. W. and Sokal R. R. Properties of Gabriel graphs relevant to geographical variation research
and the clustering of points in the same plane. Geographical Analysis 12 (1984) 205-222.
[26] Makeham J. and Buell P. D. China: The World's Oldest Living Civilization Revealed (Thames &
Hudson, 2008).
[27] Muhammad R. B. A distributed graph algorithm for geometric routing in ad hoc wireless networks.
J Networks 2 (2007) 49-57.
[28] Nakagaki T., Yamada H., and Toth A. Maze-solving by an amoeboid organism. Nature 407 (2000)
470.
[29] Nakagaki T., Yamada H., and Toth A. Path-finding by tube morphogenesis in an amoeboid organism.
Biophys. Chem. 92 (2001) 47-52.
[30] Nakagaki T., lima M., Ucda T., Nishiura Y., Saigusa T., Tero A., Kobayashi R., and Showalter K.
Minimum-risk path finding by an adaptive amoeba network. Phys. Rev. Lett. 99 (2007) 068104.
[31] Nakagaki T., Makoto I., Ueda T., Nishiura T., Saigusa T., Tero A., Kobayashi R., and Showalter K.
Minimum-risk path finding by an adaptive amoebal network. Phys. Rev. Lett. 99 (2007) 68104.
[32] Nesetril J., Milkova E., Nesctrilova H., Otakar Boruvka on minimum spanning tree problem, Discrete
Mathematics 233 (2001) 3-36.
[33] Santi P. Topology Control in Wireless Ad Hoc and Sensor Networks (Wiley, 2005).
[34] Priority Investment Needs for the Development of the Asian Highway Network. Economic and Social
Commission for Asia and the Pacific. (United Nations, New York, 2006), p. 2.
[35] Reyes D. R., Ghanem M. G., George M. Glow discharge in micro fluidic chips for visible analog
computing. Lab on a Chip 1 (2002) 113-116.
[36] Scarre C. (General Editor). Past Worlds. The Times Atlas of Archaeology. (Times Books, 1996), p.
190-191.
[37] Song W.-Z., Wang Y., Li X.-Y. Localized algorithms for energy efficient topology in wireless ad hoc
networks. In: Proc. MobiHoc 2004 (May 24-26, 2004, Roppongi, Japan).
[38] Stephenson S. L. and Stempen H. Myxomycetes: A Handbook of Slime Molds. (Timber Press, 2000).
21
Adamatzky A.
Int. J. Bifurcation Chaos, 22, 1230028 (2012) [26 pages]
Supowit K.J. The relative neighbourhood graph, with application to minimum spanning tree J. ACM
30 (1988) 428-448.
Tero A., Takagi S., Saigusa T., Ito K., Bebber D.P., Fricker M.D., Yumiki K., Kobayashi R., Nakagaki
T. Rule for biologically inspired adaptive network design. Science 327 (2010) 439-442.
Thubron C. Shadow of the Silk Road (Harper Perennial, 2008).
arXiv:
Toroczkai Z. and Guclu H. Proximity networks and epidemics. Physica A 378 (2007) 68.
phys ics/0701255vi"1
Toussaint G. T., The relative neighborhood graph of a finite planar set, Pattern Recognition 12
(1980) 261-268.
Tsuda S., Aono M., Gunji Y.-P. Robust and emergent Physarum logical-computing. Biosystems 73
(2004) 45-55.
Wan P.-J., Yi C.-W. On the longest edge of Gabriel Graphs in wireless ad hoc networks. IEEE Trans,
on Parallel and Distributed Systems 18 (2007) 111-125.
Watanabe D. A study on analyzing the road network pattern using proximity graphs. J of the City
Planning Institute of Japan 40 (2005) 133-138.
Watanabe D. Evaluating the configuration and the travel efficiency on proximity graphs as trans-
portation networks. Forma 23 (2008) 81-87.
22