Skip to main content

Full text of "The World's Colonisation and Trade Routes Formation as Imitated by Slime Mould"

See other formats


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