1
|
Ling X, Ju WB, Guo N, Zhu KJ, Wu CY, Hao QY. Effects of topological characteristics on rhythmic states of the D-dimensional Kuramoto model in complex networks. CHAOS (WOODBURY, N.Y.) 2022; 32:013118. [PMID: 35105134 DOI: 10.1063/5.0058747] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 06/02/2021] [Accepted: 11/24/2021] [Indexed: 06/14/2023]
Abstract
Synchronization is a ubiquitous phenomenon in engineering and natural ecosystems. While the dynamics of synchronization modeled by the Kuramoto model are commonly studied in two dimensions and the state of dynamic units is characterized by a scalar angle variable, we studied the Kuramoto model generalized to D dimensions in the framework of a complex network and utilized the local synchronous order parameter between the agent and its neighbors as the controllable variable to adjust the coupling strength. Here, we reported that average connectivity of networks affects the time-dependent, rhythmic, cyclic state. Importantly, we found that the level of heterogeneity of networks governs the rhythmic state in the transition process. The analytical treatment for observed scenarios in a D-dimensional Kuramoto model at D=3 was provided. These results offered a platform for a better understanding of time-dependent swarming and flocking dynamics in nature.
Collapse
Affiliation(s)
- Xiang Ling
- School of Automotive and Transportation Engineering, Hefei University of Technology, 230009 Hefei, People's Republic of China
| | - Wen-Bin Ju
- School of Automotive and Transportation Engineering, Hefei University of Technology, 230009 Hefei, People's Republic of China
| | - Ning Guo
- School of Automotive and Transportation Engineering, Hefei University of Technology, 230009 Hefei, People's Republic of China
| | - Kong-Jin Zhu
- School of Automotive and Transportation Engineering, Hefei University of Technology, 230009 Hefei, People's Republic of China
| | - Chao-Yun Wu
- School of Mathematics and Physics, Anqing Normal University, Anqing 246133, People's Republic of China
| | - Qing-Yi Hao
- School of Mathematics and Physics, Anqing Normal University, Anqing 246133, People's Republic of China
| |
Collapse
|
2
|
Bolfe M, Metz FL, Guzmán-González E, Castillo IP. Analytic solution of the two-star model with correlated degrees. Phys Rev E 2021; 104:014147. [PMID: 34412227 DOI: 10.1103/physreve.104.014147] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/19/2021] [Accepted: 07/07/2021] [Indexed: 11/06/2022]
Abstract
Exponential random graphs are important to model the structure of real-world complex networks. Here we solve the two-star model with degree-degree correlations in the sparse regime. The model constraints the average correlation between the degrees of adjacent nodes (nearest neighbors) and between the degrees at the end-points of two-stars (next nearest neighbors). We compute exactly the network free energy and show that this model undergoes a first-order transition to a condensed phase. For non-negative degree correlations between next nearest neighbors, the degree distribution inside the condensed phase has a single peak at the largest degree, while for negative degree correlations between next nearest neighbors the condensed phase is characterized by a bimodal degree distribution. We calculate the degree assortativities and show they are nonmonotonic functions of the model parameters, with a discontinuous behavior at the first-order transition. The first-order critical line terminates at a second-order critical point, whose location in the phase diagram can be accurately determined. Our results can help to develop more detailed models of complex networks with correlated degrees.
Collapse
Affiliation(s)
- Maíra Bolfe
- Physics Department, Federal University of Santa Maria, 97105-900 Santa Maria, Brazil
| | - Fernando L Metz
- Physics Institute, Federal University of Rio Grande do Sul, 91501-970 Porto Alegre, Brazil and London Mathematical Laboratory, 8 Margravine Gardens, London W6 8RH, United Kingdom
| | - Edgar Guzmán-González
- Departamento de Física, Universidad Autónoma Metropolitana-Iztapalapa, San Rafael Atlixco 186, Ciudad de México 09340, México
| | - Isaac Pérez Castillo
- Departamento de Física, Universidad Autónoma Metropolitana-Iztapalapa, San Rafael Atlixco 186, Ciudad de México 09340, México
| |
Collapse
|
3
|
Roy M, Poria S, Hens C. Assortativity-induced explosive synchronization in a complex neuronal network. Phys Rev E 2021; 103:062307. [PMID: 34271687 DOI: 10.1103/physreve.103.062307] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/25/2020] [Accepted: 05/25/2021] [Indexed: 11/07/2022]
Abstract
In this study, we consider a scale-free network of nonidentical Chialvo neurons, coupled through electrical synapses. For sufficiently strong coupling, the system undergoes a transition from completely out of phase synchronized to phase synchronized state. The principal focus of this study is to investigate the effect of the degree of assortativity over the synchronization transition process. It is observed that, depending on assortativity, bistability between two asymptotically stable states allows one to develop a hysteresis loop which transforms the phase transition from second order to first order. An expansion in the area of hysteresis loop is noticeable with increasing degree-degree correlation in the network. Our study also reveals that effective frequencies of nodes simultaneously go through a continuous or sudden transition to the synchronized state with the corresponding phases. Further, we examine the robustness of the results under the effect of network size and average degree, as well as diverse frequency setup. Finally, we investigate the dynamical mechanism in the process of generating explosive synchronization. We observe a significant impact of lower degree nodes behind such phenomena: in a positive assortative network the low degree nodes delay the synchronization transition.
Collapse
Affiliation(s)
- Mousumi Roy
- Department of Applied Mathematics, University of Calcutta, 92, A.P.C. Road, Kolkata 700009, India
| | - Swarup Poria
- Department of Applied Mathematics, University of Calcutta, 92, A.P.C. Road, Kolkata 700009, India
| | - Chittaranjan Hens
- Physics and Applied Mathematics Unit, Indian Statistical Institute, Kolkata 700108, India
| |
Collapse
|
4
|
Letellier C. Assessing synchronizability provided by coupling variable from the algebraic structure of dynamical systems. Phys Rev E 2020; 101:042215. [PMID: 32422746 DOI: 10.1103/physreve.101.042215] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/26/2019] [Accepted: 04/05/2020] [Indexed: 06/11/2023]
Abstract
Synchronization is a very generic phenomenon which can be encountered in a large variety of coupled dynamical systems. Being able to synchronize chaotic systems is strongly dependent on the nature of their coupling. Few attempts to explain such a dependency using observability and/or controllability were not fully satisfactory and synchronizability yet remained unexplained. Synchronizability can be defined as the range of coupling parameter values for which two nearly identical systems are fully synchronized. Our objective is here to investigate whether synchronizability can be related to the main rotation necessarily required for structuring any type of attractor, that is, whether synchronizability is significantly improved when the coupling variable is one of the variables involved in the main rotation. We thus propose a semianalytic procedure from a single isolated system to discard the worst variable for fully synchronizing two (nearly) identical copies of that system.
Collapse
Affiliation(s)
- Christophe Letellier
- Normandie Université-CORIA, Campus Universitaire du Madrillet, F-76800 Saint-Etienne du Rouvray, France
| |
Collapse
|
5
|
Gao J, Efstathiou K. Reduction of oscillator dynamics on complex networks to dynamics on complete graphs through virtual frequencies. Phys Rev E 2020; 101:022302. [PMID: 32168684 DOI: 10.1103/physreve.101.022302] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/22/2019] [Accepted: 01/22/2020] [Indexed: 11/06/2022]
Abstract
We consider the synchronization of oscillators in complex networks where there is an interplay between the oscillator dynamics and the network topology. Through a remarkable transformation in parameter space and the introduction of virtual frequencies we show that Kuramoto oscillators on annealed networks, with or without frequency-degree correlation, and Kuramoto oscillators on complete graphs with frequency-weighted coupling can be transformed to Kuramoto oscillators on complete graphs with a rearranged, virtual frequency distribution and uniform coupling. The virtual frequency distribution encodes both the natural frequency distribution (dynamics) and the degree distribution (topology). We apply this transformation to give direct explanations to a variety of phenomena that have been observed in complex networks, such as explosive synchronization and vanishing synchronization onset.
Collapse
Affiliation(s)
- Jian Gao
- Bernoulli Institute for Mathematics, Computer Science, and Artificial Intelligence, University of Groningen, P.O. Box 407, 9700 AK, Groningen, The Netherlands
| | - Konstantinos Efstathiou
- Division of Natural and Applied Sciences and Zu Chongzhi Center for Mathematics and Computational Science, Duke Kunshan University, No. 8 Duke Avenue, Kunshan 215316, China
| |
Collapse
|
6
|
Liu Y, Kurths J. Effects of network robustness on explosive synchronization. Phys Rev E 2019; 100:012312. [PMID: 31499821 DOI: 10.1103/physreve.100.012312] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/21/2019] [Indexed: 11/07/2022]
Abstract
Current studies have shown that there is a positive correlation between the network assortativity and robustness and that the assortativity also plays an important role in explosive synchronization. In this paper, taking the network robustness as a global property, we investigate its significance as well as the influence of its interaction with the assortativity on explosive synchronization. Our numerical results demonstrate that explosive synchronization is suppressed in extreme situations of both the robustness and assortativity. In addition, through appropriate adjustments of them, a maximum hysteresis area between the forward and backward transitions can be reached. Furthermore, our results might also provide reference for those who are interested in effects of network structure on synchronization, though this problem is still challenging as we show in the discussion.
Collapse
Affiliation(s)
- Yang Liu
- Potsdam Institute for Climate Impact Research, 14412 Potsdam, Germany.,Department of Computer Science, Technische Universität Berlin, 10587 Berlin, Germany
| | - Jürgen Kurths
- Potsdam Institute for Climate Impact Research, 14412 Potsdam, Germany.,Department of Physics, Humboldt University Berlin, 12489 Berlin, Germany.,Saratov State University, 410012 Saratov, Russia
| |
Collapse
|
7
|
Tian C, Bi H, Zhang X, Guan S, Liu Z. Asymmetric couplings enhance the transition from chimera state to synchronization. Phys Rev E 2017; 96:052209. [PMID: 29347748 DOI: 10.1103/physreve.96.052209] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/04/2017] [Indexed: 06/07/2023]
Abstract
Chimera state has been well studied recently, but little attention has been paid to its transition to synchronization. We study this topic here by considering two groups of adaptively coupled Kuramoto oscillators. By searching the final states of different initial conditions, we find that the system can easily show a chimera state with robustness to initial conditions, in contrast to the sensitive dependence of chimera state on initial conditions in previous studies. Further, we show that, in the case of symmetric couplings, the behaviors of the two groups are always complementary to each other, i.e., robustness of chimera state, except a small basin of synchronization. Interestingly, we reveal that the basin of synchronization will be significantly increased when either the coupling of inner groups or that of intergroups are asymmetric. This transition from the attractor of chimera state to the attractor of synchronization is closely related to both the phase delay and the asymmetric degree of coupling strengths, resulting in a diversity of attractor's patterns. A theory based on the Ott-Antonsen ansatz is given to explain the numerical simulations. This finding may be meaningful for the control of competition between two attractors in biological systems, such as the cardiac rhythm and ventricular fibrillation, etc.
Collapse
Affiliation(s)
- Changhai Tian
- Department of Physics, East China Normal University, Shanghai 200062, People's Republic of China
- School of Data Science, Tongren University, Tongren 554300, People's Republic of China
| | - Hongjie Bi
- Department of Physics, East China Normal University, Shanghai 200062, People's Republic of China
| | - Xiyun Zhang
- Department of Physics, East China Normal University, Shanghai 200062, People's Republic of China
| | - Shuguang Guan
- Department of Physics, East China Normal University, Shanghai 200062, People's Republic of China
| | - Zonghua Liu
- Department of Physics, East China Normal University, Shanghai 200062, People's Republic of China
| |
Collapse
|
8
|
Jalan S, Kumar A, Zaikin A, Kurths J. Interplay of degree correlations and cluster synchronization. Phys Rev E 2016; 94:062202. [PMID: 28085396 DOI: 10.1103/physreve.94.062202] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/16/2016] [Indexed: 06/06/2023]
Abstract
We study the evolution of coupled chaotic dynamics on networks and investigate the role of degree-degree correlation in the networks' cluster synchronizability. We find that an increase in the disassortativity can lead to an increase or a decrease in the cluster synchronizability depending on the degree distribution and average connectivity of the network. Networks with heterogeneous degree distribution exhibit significant changes in cluster synchronizability as well as in the phenomena behind cluster synchronization as compared to those of homogeneous networks. Interestingly, cluster synchronizability of a network may be very different from global synchronizability due to the presence of the driven phenomenon behind the cluster formation. Furthermore, we show how degeneracy at the zero eigenvalues provides an understanding of the occurrence of the driven phenomenon behind the synchronization in disassortative networks. The results demonstrate the importance of degree-degree correlations in determining cluster synchronization behavior of complex networks and hence have potential applications in understanding and predicting dynamical behavior of complex systems ranging from brain to social systems.
Collapse
Affiliation(s)
- Sarika Jalan
- Complex Systems Lab, Discipline of Physics, Indian Institute of Technology Indore, Simrol, Indore 453552, India
- Centre for Biosciences and Biomedical Engineering, Indian Institute of Technology Indore, Simrol, Indore 453552, India
| | - Anil Kumar
- Complex Systems Lab, Discipline of Physics, Indian Institute of Technology Indore, Simrol, Indore 453552, India
| | - Alexey Zaikin
- Institute for Women's Health and Department of Mathematics, University College London, London WC1H 0AY, United Kingdom
- Department of Applied Mathematics and Centre of Bioinformatics, Lobachevsky State University of Nizhni Novgorod, Nizhni Novgorod 603950, Russia
| | - Jürgen Kurths
- Potsdam Institute for Climate Impact Research, D-14412 Potsdam, Germany
| |
Collapse
|
9
|
Kim M, Mashour GA, Moraes SB, Vanini G, Tarnal V, Janke E, Hudetz AG, Lee U. Functional and Topological Conditions for Explosive Synchronization Develop in Human Brain Networks with the Onset of Anesthetic-Induced Unconsciousness. Front Comput Neurosci 2016; 10:1. [PMID: 26834616 PMCID: PMC4720783 DOI: 10.3389/fncom.2016.00001] [Citation(s) in RCA: 38] [Impact Index Per Article: 4.2] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/26/2015] [Accepted: 01/04/2016] [Indexed: 11/18/2022] Open
Abstract
Sleep, anesthesia, and coma share a number of neural features but the recovery profiles are radically different. To understand the mechanisms of reversibility of unconsciousness at the network level, we studied the conditions for gradual and abrupt transitions in conscious and anesthetized states. We hypothesized that the conditions for explosive synchronization (ES) in human brain networks would be present in the anesthetized brain just over the threshold of unconsciousness. To test this hypothesis, functional brain networks were constructed from multi-channel electroencephalogram (EEG) recordings in seven healthy subjects across conscious, unconscious, and recovery states. We analyzed four variables that are involved in facilitating ES in generic, non-biological networks: (1) correlation between node degree and frequency, (2) disassortativity (i.e., the tendency of highly-connected nodes to link with less-connected nodes, or vice versa), (3) frequency difference of coupled nodes, and (4) an inequality relationship between local and global network properties, which is referred to as the suppressive rule. We observed that the four network conditions for ES were satisfied in the unconscious state. Conditions for ES in the human brain suggest a potential mechanism for rapid recovery from the lightly-anesthetized state. This study demonstrates for the first time that the network conditions for ES, formerly shown in generic networks only, are present in empirically-derived functional brain networks. Further investigations with deep anesthesia, sleep, and coma could provide insight into the underlying causes of variability in recovery profiles of these unconscious states.
Collapse
Affiliation(s)
- Minkyung Kim
- Department of Anesthesiology, University of Michigan Medical SchoolAnn Arbor, MI, USA; Center for Consciousness Science, University of Michigan Medical SchoolAnn Arbor, MI, USA; Department of Physics, Pohang University of Science and TechnologyPohang, South Korea
| | - George A Mashour
- Department of Anesthesiology, University of Michigan Medical SchoolAnn Arbor, MI, USA; Center for Consciousness Science, University of Michigan Medical SchoolAnn Arbor, MI, USA; Neuroscience Graduate Program, University of Michigan Medical SchoolAnn Arbor, MI, USA
| | - Stefanie-Blain Moraes
- Department of Anesthesiology, University of Michigan Medical School Ann Arbor, MI, USA
| | - Giancarlo Vanini
- Department of Anesthesiology, University of Michigan Medical School Ann Arbor, MI, USA
| | - Vijay Tarnal
- Department of Anesthesiology, University of Michigan Medical School Ann Arbor, MI, USA
| | - Ellen Janke
- Department of Anesthesiology, University of Michigan Medical School Ann Arbor, MI, USA
| | - Anthony G Hudetz
- Department of Anesthesiology, University of Michigan Medical SchoolAnn Arbor, MI, USA; Center for Consciousness Science, University of Michigan Medical SchoolAnn Arbor, MI, USA; Neuroscience Graduate Program, University of Michigan Medical SchoolAnn Arbor, MI, USA
| | - Uncheol Lee
- Department of Anesthesiology, University of Michigan Medical SchoolAnn Arbor, MI, USA; Center for Consciousness Science, University of Michigan Medical SchoolAnn Arbor, MI, USA
| |
Collapse
|
10
|
Qu J, Wang SJ, Jusup M, Wang Z. Effects of random rewiring on the degree correlation of scale-free networks. Sci Rep 2015; 5:15450. [PMID: 26482005 PMCID: PMC4611853 DOI: 10.1038/srep15450] [Citation(s) in RCA: 7] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/30/2015] [Accepted: 09/08/2015] [Indexed: 01/01/2023] Open
Abstract
Random rewiring is used to generate null networks for the purpose of analyzing the topological properties of scale-free networks, yet the effects of random rewiring on the degree correlation are subject to contradicting interpretations in the literature. We comprehensively analyze the degree correlation of randomly rewired scale-free networks and show that random rewiring increases disassortativity by reducing the average degree of the nearest neighbors of high-degree nodes. The effect can be captured by the measures of the degree correlation that consider all links in the network, but not by analogous measures that consider only links between degree peers, hence the potential for contradicting interpretations. We furthermore find that random and directional rewiring affect the topology of a scale-free network differently, even if the degree correlation of the rewired networks is the same. Consequently, the network dynamics is changed, which is proven here by means of the biased random walk.
Collapse
Affiliation(s)
- Jing Qu
- School of Physics and Information Technology, Shaanxi Normal University, Xi’an 710119, China
| | - Sheng-Jun Wang
- School of Physics and Information Technology, Shaanxi Normal University, Xi’an 710119, China
| | - Marko Jusup
- Faculty of Sciences, Kyushu University, Fukuoka 819-0395, Japan
| | - Zhen Wang
- Interdisciplinary Graduate School of Engineering Sciences, Kyushu University, Fukuoka 816-8580, Japan
- School of Automation, Northwestern Polytechnical University, Xi’an 710072, China
| |
Collapse
|
11
|
Sadilek M, Thurner S. Physiologically motivated multiplex Kuramoto model describes phase diagram of cortical activity. Sci Rep 2015; 5:10015. [PMID: 25996547 PMCID: PMC4650820 DOI: 10.1038/srep10015] [Citation(s) in RCA: 28] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/03/2014] [Accepted: 03/18/2015] [Indexed: 12/04/2022] Open
Abstract
We derive a two-layer multiplex Kuramoto model from Wilson-Cowan type physiological equations that describe neural activity on a network of interconnected cortical regions. This is mathematically possible due to the existence of a unique, stable limit cycle, weak coupling, and inhibitory synaptic time delays. We study the phase diagram of this model numerically as a function of the inter-regional connection strength that is related to cerebral blood flow, and a phase shift parameter that is associated with synaptic GABA concentrations. We find three macroscopic phases of cortical activity: background activity (unsynchronized oscillations), epileptiform activity (highly synchronized oscillations) and resting-state activity (synchronized clusters/chaotic behaviour). Previous network models could hitherto not explain the existence of all three phases. We further observe a shift of the average oscillation frequency towards lower values together with the appearance of coherent slow oscillations at the transition from resting-state to epileptiform activity. This observation is fully in line with experimental data and could explain the influence of GABAergic drugs both on gamma oscillations and epileptic states. Compared to previous models for gamma oscillations and resting-state activity, the multiplex Kuramoto model not only provides a unifying framework, but also has a direct connection to measurable physiological parameters.
Collapse
Affiliation(s)
- Maximilian Sadilek
- Section for Science of Complex Systems, Medical University of Vienna, Spitalgasse 23, A-1090 Vienna, Austria
| | - Stefan Thurner
- Section for Science of Complex Systems, Medical University of Vienna, Spitalgasse 23, A-1090 Vienna, Austria
- Santa Fe Institute, 1399 Hyde Park Road, New Mexico 87501, USA
- IIASA, Schlossplatz 1, A-2361 Laxenburg, Austria
| |
Collapse
|
12
|
Peron TKD, Ji P, Rodrigues FA, Kurths J. Effects of assortative mixing in the second-order Kuramoto model. PHYSICAL REVIEW. E, STATISTICAL, NONLINEAR, AND SOFT MATTER PHYSICS 2015; 91:052805. [PMID: 26066210 DOI: 10.1103/physreve.91.052805] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 10/15/2014] [Indexed: 06/04/2023]
Abstract
In this paper we analyze the second-order Kuramoto model in the presence of a positive correlation between the heterogeneity of the connections and the natural frequencies in scale-free networks. We numerically show that discontinuous transitions emerge not just in disassortative but also in strongly assortative networks, in contrast with the first-order model. We also find that the effect of assortativity on network synchronization can be compensated by adjusting the phase damping. Our results show that it is possible to control collective behavior of damped Kuramoto oscillators by tuning the network structure or by adjusting the dissipation related to the phases' movement.
Collapse
Affiliation(s)
- Thomas K Dm Peron
- Instituto de Física de São Carlos, Universidade de São Paulo, São Carlos, São Paulo, Brazil
| | - Peng Ji
- Potsdam Institute for Climate Impact Research (PIK), 14473 Potsdam, Germany
- Department of Physics, Humboldt University, 12489 Berlin, Germany
| | - Francisco A Rodrigues
- Departamento de Matemática Aplicada e Estatística, Instituto de Ciências Matemáticas e de Computação, Universidade de São Paulo, Caixa Postal 668, 13560-970 São Carlos, São Paulo, Brazil
| | - Jürgen Kurths
- Potsdam Institute for Climate Impact Research (PIK), 14473 Potsdam, Germany
- Department of Physics, Humboldt University, 12489 Berlin, Germany
- Institute for Complex Systems and Mathematical Biology, University of Aberdeen, Aberdeen AB24 3UE, United Kingdom
| |
Collapse
|