Click here to login

Statistical Mechanics & Dynamical Systems Laboratory

Group leader:  Dr. Provata Astero

E-mail: a.provata “at”


PhD Students

Kasimatis Theodoros (t.kasimatis “at”

Tsigkri Nefeli-Dimitra (n.tsigkri “at”

MSc Students

Rontogiannis  Alexandros

External Collaborators

Prof. Andreou Christina  (Lubeck University, Germany)

Prof. Antonopoulos Christos  (Essex University, UK)

Prof. Mandilara Aikaterini  (Nazarbayev University, Kazakhstan)

Prof. Schoell Eckehard  (Technical Universitaet of Berlin, Germany)

Prof. Vlamos Takis  (Ionian University, Greece)

Previous Students

Argyropoulos George (MSc, 2017-2020)

Koulierakis Ioannis (MSc, 2017-2020)

Karakos George (MSc 2017-2019)

Labrakos Christos (MSc, 2017-2018)

Georgiou Antoine (BSc, 2015-2016)

Thanos Dimitrios (BSc, 2015-2018)

Kasimatis Theodoros (MA 2014-2016)

Tsigkri Nefeli-Dimitra (MA 2013-2015)

Dr. Hizanidis Johanne  (Postdoc, 2013-2017)

Flokas Lambros (BA 2015)

Mousa George (BSc 2016)

Theodoropoulos Spyros (BSc 2016)

Breki Christina-Marina (BA 2014-2015)

Dr. Bassett Jason (BSc, 2012-2013)

Vogiatzian Philippos-Arthuros  (BSc, 2012-2013)

Dr. Kovanis Michail  (BSc, 2012)

Dr. Katsaloulis Panayotis  (MSc, 2006-2007, PhD, 2007-2011)

Dr. Kouvaris Nikos  (PhD, 2007-2011)

Dr. Noussiou Vicky  (PhD, 2004-2009)

Prof. Oikonomou Thomas  (PhD, 2004-2008)

Dr. Rabias Ioannis  (Postdoc, 2002)

Dr. Sfyrakis Kostas  (Postdoc, 2002-2003)

Prof. Prassas Cyril Prassas (DEA, 2001)

Dr. Prassas Vassilios  (MSc, 1998)

Dr. Chantron Arnaud  (DEA, 1999)

Prof. Kalosakas George (postdoc, 2001)

Prof. Shabunin Alexey  (Postdoc, 1999,2001,2003,2007)

Prof. Chichigina Olga  (Postdoc, 2004)

Chondrou Presveia  (MSc 2003-2004)

Dr. Tsekouras Georgios-Artemios  (PhD, 2000-2004)


Group Photos

Group Photos 2017

Group Photos 2016




Group Photos 2015



The Laboratory of Statistical Mechanics and Dynamical Systems (STAT-DYN) was founded in 02/2004 as part of the Institute of Physical Chemistry, while from 2014 belongs to the Institute of Nanoscience and Nanotechnology. Its research focuses on the fields of :

    1. Non-linear Dynamics,
    2. Statistical Mechanics
    3. Networks of Interacting Neurons
    4. Computational Neuroscience
    5. Brain Dynamics
    6. Fractals
    7. Complex Systems and
    8. Network Science

      Images of health brain: Left: Frequency plot, Right: Chimera state

Reaction Diffusion process on 3D fractal substrate

Our aim is the development of methods and models for understanding the emergence and evolution of mesoscopic and macroscopic spatial patterns and temporal synchronization motifs due to the interactions between nonlinear elements (e.g., neuronal oscillators or population dynamics models). The spatiotemporal structures induced by the interacting units include chimera states in neural networks, aggregates, spirals and stripe formations in reaction diffusion systems, helices, fractals, synchronisation phenomena etc. which can be experimentally observed in material science, physics, chemistry and biology. Our studies in particular include research on fractal pattern formation and correlations near critical points but also research in open systems, dissipative in constant exchange with the environment. Away from the critical points in closed, isolated systems short range correlations and spatiotemporal patterns with well-defined length and time scales are studied (eg. spiral and stripe formations, helices etc.). At the critical points long range correlations are developed and information exchanges are studied at the long scales. In open, coupled dissipative systems the reduction of the phase space is studied, which leads to unexpected spatiotemporal phenomena such as the chimera states in coupled neuronal networks. The understanding of these structures at the micro-, meso- and macro scales and the interaction between these three levels of description has major technological impact in materials science and physical, chemical and biological processes.

For the study of complex systems in the lab we develop a) statistical methods and algorithms describing the evolution of complex morphologies and b) modelling tools for the dynamics of pattern formation and synchronization phenomena. Statistical methods include thermodynamic approaches, entropic (extensive and non-extensive) approaches, theory of complex networks, theory of long and short range distributions, Levi distributions, theory of random walks. For the study of the mechanisms creating complex patterns, non-linear dynamical systems of hierarchical complexity are used, together with mean-field theories, exact enumeration methods, real space renormalisation theory, theory of stochastic processes and numerical integration and kinetic Monte Carlo Methods.

Applications include, among others:

    • Spatial and temporal synchronization properties in neuronal networks

    • Complex brain dynamics

    • Analysis and modelling of the fractal architecture of the neuron axons spanning the human brain

    • Reaction-diffusion processes in Complex networks,

    • Studies of surface phenomena and aggregates,

    • Fractals: formation and evolution

    • Bioinformatics,

    • Statistical analysis and modelling of biological tissues and macromolecules,

Typical recent studies include:

1) Analysis and modelling of the complex fractal architecture of the neuron axons spanning the human brain

Based on the classical Magnetic Resonance Imaging (MRI) technique it is now possible to obtain high resolution brain images reaching scales as low as ~0.2 mm. Diffusion Tensor MRI tomography (DTI-MRI) allows for capturing the 3D motion of the water diffusion in the brain area by using the anisoptropic movement of water molecules to capture the positions of the axons connecting the neurons, see Fig. 1 for reconstruction of the neuron axons network from MRI recordings. Based on thorough analysis of MRI images from may healthy brains we have demonstrated self-similarity in the form of a hierarchical architecture of the connectivity matrix of the neuron axons network, which is proper to the healthy state of the brain, while its fractal dimensions were calculated to df~2.5. Using fractal and multifractal analyses we show that these complex connectivity patterns cannot be a result of simple or even chaotic processes because rare, unexpected configurations are present, indicative of the influence of long range exchange mechanisms. A further challenge is to extend the already existing models to cover the cases of pathological brain architectures (Alzheimer, Parkinson and Schizophrenia), a problem which has not been undertaken thus far.

Fig. 1. Typical neuron axons network representation from DTI-MRI images

In parallel with the analysis of the detailed structure of the neuron axons network, we use numerical simulations to illustrate brain dynamics, using precise network connectivity matrices as recorded by the MRI analysis. The hierarchical connectivity in the brain, we consider as a key factor for the development of the intricate spatiotemporal patterns covering multiple scales, necessary for the emergence of cognition. Coupling schemes between a large number of neurons follow three main connectivity patterns: all-to-all (global) coupling, nearest neighbours (or local ) coupling and non-local coupling where each neuron is linked with a finite number of neighbours. The last coupling pattern under special conditions produces the so-called ”chimera states”, where some of the oscillators (neurons) oscillate in phase while others are incoherent, see Fig. 2.

Fig. 2. Exemplary chimera state in 1D, 2D and 3D networks with nonlocal connectivity. We note that synchronous and asynchronous regions coexist simultaneously in the network.


The coexistence of coherent and incoherent states have been associated with partial brain functionality, as in simple brain tasks which engage only a few modules of the brain or in uni-hemispheric sleep. We believe that the presence of chimera states is far more important than proposed for the unihemispheric sleep and is directly related to brain functionality during, at least, simple tasks.

2) Reaction Diffusion Processes on Complex Networks

Abstract networks can be used for the description of the dynamics in reactive processes and they can give valuable information about the evolution and the products of each process. The abstract reactive networks are constructed as follows: The phase space of the variables (concentrations in reactive systems) is partitioned into a finite number of segments, which constitute the nodes of the abstract network. Transitions between the nodes are dictated by the dynamics of the reactive process and provide the links between the nodes. These are weighted networks, since each link weight reflects the transition rate between the corresponding states-nodes. With this construction the network properties mirror the dynamics of the underlying process and one can investigate the system properties by studying the corresponding abstract network.

As an example, the abstract network of the Lattice Limit Cycle (LLC) model was analysed and its transition matrix elements were computed via Kinetic (Dynamic) Monte Carlo simulations. For this model it was shown that the degree distribution follows a power law with exponent -1, while the average clustering coefficient c(N) scales with the network size N as with a power law exponent 1.46. The computed exponents classify the LLC abstract reactive network into the scale-free networks. This conclusion corroborates earlier investigations demonstrating the formation of fractal spatial patterns in LLC reactive dynamics due to stochasticity and to the clustering of homologous species. The present construction of abstract networks (based on the partition of the phase space) is generic and can be implemented with appropriate adjustments in many dynamical systems and in time series analysis.

Fig. 3: Representative snapshots of the LLC system from Kinetic Monte Carlo simulations. An initially homogeneous system (t=0 MCS) develops clusters and inhomogeneous distribution of species (t=10000 MCS).

3) Bioinformatics

The complexity of the primary structure of human DNA is explored using methods from nonequilibrium statistical mechanics, dynamical systems theory and information theory. A collection of statistical analyses are performed on the DNA data and the results are compared across species and with sequences derived from different stochastic processes. Although detailed balance seems to hold at the level of a binary alphabet in genomic data, it fails when all four basepairs are considered, suggesting spatial asymmetry and irreversibility. Furthermore, the block entropy does not increase linearly with the block size, reflecting the long range nature of the correlations in the human genomic sequences. To probe locally the spatial structure of the chain we study local quantities, such as the exit distances from a specific symbol, the distribution of recurrence distances and the Hurst exponent, all of which show power law tails and long range characteristics. These results suggest that human DNA can be viewed as a non-equilibrium structure maintained in its state through interactions with a constantly changing environment.

Fig. 4 : Exemplary Information Transfer Analysis between a sequence and its shift by n symbols in the 4-letter representation. Line with circles depicts the DNA sequence, line with crosses a random sequence and line with diamonds a model-generated sequence.

Complexity measures are also used to compare the genomic characteristics of different organisms belonging to distinct classes spanning the evolutionary tree: higher eukaryotes, amoebae, unicellular eukaryotes and bacteria. We have demonstrated that the conditional probability matrix for the four-letter and AT-CG alphabet is markedly asymmetric in eukaryotes while it is nearly symmetric in bacterial genomes. Overall, the conditional probability, the fluxes, the block entropy content and the exit distance distributions can be used as markers, discriminating between eukaryotic and prokaryotic DNA, allowing in many cases to discern details related to finer classes. In all cases the reduction from four letters to two mask some important statistical and spatial properties, with the AT-CG alphabet having higher ability of discrimination than the AG-CT one.