The rich landscape of iterated simplicity

Cellular automata and network automata

Michiel Rollier

Public PhD defence ยท Friday 2 October 2026

Supervisors: Prof. Dr. Ir. Jan M. Baetens ยท Prof. Dr. Odemir M. Bruno ยท Dr. Aisling Daly

Joint PhD ยท Ghent University ยท University of Sรฃo Paulo

Ghent UniversityFaculty of Bioscience EngineeringUniversity of Sรฃo PauloBionamiX

The title of the thesis, over a landscape in the style of its cover. It is the result of a five-year joint PhD at BionamiX (Faculty of Bioscience Engineering, Ghent University, under Jan M. Baetens) and the Scientific Computing Group (Sรฃo Carlos Institute of Physics, University of Sรฃo Paulo, under Odemir M. Bruno). More on the research page and on ORCID.

The picture is not explained yet: within ten minutes the room will know what it is, and it comes back on the last slide.

A small experiment

  1. 1 Pick two people near you
  2. 2 Ignore your own hand
  3. 3 Wait for the signal

0 up

1 up

2 up

The talk starts cold, with a small experiment in participatory science: nothing is explained yet. Everyone picks two people near them and keeps them for the whole game. Each round you look only at their hands: if exactly one of them has a hand up, yours goes up; if none or both do, yours goes down. Your own hand does not count, and neither does which of the two is which: only how many are up.

Everybody first looks and decides, and only moves on the signal. That keeps the whole room in step.

Hands up if you were born in the first six months of the year

January to June

The start is random, but in a way everyone can check for themselves: hands up if you were born between January and June.

Practice 1

0 up

1 up

2 up

Two practice rounds, then back to the start and five real rounds. Each click starts a three-second countdown, which the room says out loud: three, two, one, switch. The rule stays on screen the whole time. Nothing on screen shows the pattern the room makes: the room itself is the computer.

You just simulated the parity rule!

Nobody saw the whole. Everybody followed the same simple, local rule.

That was the parity rule: your hand goes up when an odd number of your neighbours have theirs up. This is what a row of cells following it does, drawn round by round: space runs to the right, time runs down, each round one row lower. From a random start, like tonightโ€™s, it makes a jumble of triangles; from a single raised hand in the middle it draws a Sierpiล„ski triangle. The reset button draws a new random start.

Nobody saw the whole pattern; everyone only followed a simple, local rule. That is a cellular automaton, and that was the game: from here on, the cells are squares on a screen.

Elementary cellular automata (256)

0
1
0
1
1
0
1
0

01011010 in binary = rule 90

rule 168 dies out
rule 1 repeats
rule 30 looks random
rule 110 builds structures

The parity rule is one of many. A cell looks at itself and its two neighbours, which makes eight possible situations; a rule says, for each of them, whether the cell is on or off next. Read those eight answers as a binary number and you get the ruleโ€™s name: 01011010 is 90, so the parity rule is rule 90. Eight yes-or-no answers make 256 rules in all, the elementary cellular automata. Click an answer in the table to flip it and see which rule you get.

Four of them, started one by one from the same random row: rule 168 dies out, rule 1 repeats itself, blinking on and off, rule 30 looks random, and rule 110 builds structures that move and collide. Same recipe, very different outcomes. Each reset draws a new random row.

The Game of Life

Moving to two spatial dimensions

Eight neighbours

Moving to two spatial dimensions: John Conwayโ€™s Game of Life, from 1970. Every cell now has eight neighbours, the Moore neighbourhood, and counts how many of them are alive. A dead cell with exactly three live neighbours is born; a live cell with two or three survives; with fewer it dies of loneliness, with more of overpopulation.

From a random soup the grid churns for a while and then calms down into small blocks, blinkers and the odd traveller; a reset draws a new soup. These four rules are enough to build machines: this is Bill Gosperโ€™s glider gun, which fires a glider every thirty steps, forever, into an eater that swallows each one.

Still lifes

Patterns that the rule leaves exactly as they are.

Every cell here obeys the Game of Life. Nothing changes, until you disturb it.

Our result Create artwork for an international conference

Some patterns are left exactly as they are by the rule: every live cell has two or three live neighbours, and no dead cell has exactly three. These are still lifes. The five here are tiles of the mosaics of the cover, of five sizes (levels 2 to 6).

Fit many of them together and you can make a picture in which every cell obeys the Game of Life and nothing ever changes. The cover of the thesis is made this way. The colours are only decoration: drain them away and every cell is simply alive or dead, and still nothing changes. Take away one cell and the picture starts to unravel. A reset packs a new landscape of tiles; on this page you can also click or drag on the picture to disturb it yourself.

Life-like cellular automata (262โ€ฏ144)

Just as rule 90 is one of 256 elementary rules, the Game of Life is one of many Life-like rules: a dead cell is born with some numbers of live neighbours (B), a live cell survives with others (S). The Game of Life is B3/S23; there are 218=262โ€‰1442^{18} = 262\,144 such rules, and again they behave very differently.

Three of them, each started with its own button: Maze (B3/S12345) freezes into corridors within a few steps, the Game of Life keeps something going that is hard to name, and Replicator (B1357/S1357), the parity rule on eight neighbours, looks like noise. Each reset draws a new random start.

screen

Generalising even further

  1. Elementary cellular automata
  2. Life-like cellular automata
  3. Life-like network automata

What if the room played the game differently? First, from here on the people become nodes, and a hand up becomes a node that is on.

In the game you looked at your two neighbours, so the room was really one long row: link everybody to the person before and after them, snaking up from the back of the room to the front and round again, and everyone has two neighbours. That is a one-dimensional cellular automaton, an elementary one, like the rows of cells we just saw.

Now link everyone to the eight people around them: a grid, like the one the Game of Life and the other Life-like cellular automata live on. The room wraps round at its edges, so the people at a wall have eight too; the stubs at the edges are links to the far side.

Then take some of those links and swing them across the room to anyone at all. A few long links like these make a small-world network: everyone is suddenly only a few steps from everyone else, and running the same Life-like rule for everyone on it gives a Life-like network automaton. On this page, drag the slider: at 0 the grid comes back, at 100% the links are completely random.

But now some have four contacts and some thirteen, so a rule like โ€œexactly three neighboursโ€ no longer makes sense for everybody.

Life-like network automata (โˆž)

Being born

Surviving

ฯ•6,115\phi^{\htmlClass{phi-r}{5}}_{\htmlClass{phi-b}{6},\htmlClass{phi-s}{11}}

Our result Defining LLNAs to allow symmetry

A node with ten neighbours and a node with three cannot follow a rule that counts neighbours. So count the share instead: the fraction of a nodeโ€™s neighbours that are on, which runs round each ring from 0 at the top, clockwise, to 1, from light to dark. Cut it into intervals, five here. A node that is off (left) is born when the share falls in one of its chosen intervals; a node that is on (right) survives when the share falls in one of its chosen intervals. The chosen intervals are the outlined ones. An end drawn solid belongs to the interval, a dotted one does not: in ฯ•6,115\phi^5_{6,11} a node that is off is born from 15\tfrac15 up to and including 35\tfrac35. On this page, click an interval to choose it or not. The ruleโ€™s name reads the chosen intervals as two binary numbers, as rule 90 did on the elementary table: born is 6, the second and third intervals, and survive is 11, the first, second and fourth.

Each click cuts the share more finely, into 7, 9 and then 11 intervals: the resolution rr. The rule keeps the regions it had, but its name changes, and every finer resolution can write rules that a coarser one cannot: there is no end to them. At r=9r = 9, on a grid where everyone has eight neighbours, each interval holds exactly one number of live neighbours, and the Game of Life is ฯ•8,129\phi^9_{8,12}. The slider on this page sets any resolution from 2 to 13.

Defining these rules properly, for any network and at any resolution, is the first contribution of the thesis (Chapter 5, Physica D 2025), and the point of the definition is symmetry. Which ends each interval includes is chosen so that every rule has a cousin that behaves exactly the same with on and off swapped: the same patterns, in the other colours. The button on this page swaps a rule for its cousin: born becomes the complement of survive read backwards, and survive the complement of born read backwards; some rules are their own cousin. At an odd resolution the intervals below 12\tfrac12 include their left end, those above it their right end, and the middle one both. An even resolution has no middle interval, so there the two rings give 12\tfrac12 to different sides: the born ring to the interval just left of it, the survive ring to the one just right of it (Appendix C of the thesis).

The network matters too

ฯ•168,4869\phi^{9}_{168,486}

random
small-world ring
scale-free

Life-like network automata can run any rule on any network, and the network matters as much as the rule. Here one rule, ฯ•168,4869\phi^9_{168,486}, on the three networks of the rest of the talk, 150 nodes each, from a random start with three in ten on; the click starts all three. On the random network it never settles; on the small-world ring it fills up and freezes; on the scale-free network it all but dies out. Why? On the scale-free network 80 of the 150 nodes, more than half, have exactly four links, and a node with four links only ever sees none, a quarter, half, three quarters or all of its neighbours on: shares in which this rule never lets a node be born. So they can never switch on. On the ring every node has about eight neighbours, and on the random network anything from two to thirteen. Same rule, three different behaviours: a great many combinations to study. (On the page, โ†บ draws new starts.)

A word on the names. The random network is an Erdล‘sโ€“Rรฉnyi network: every pair of nodes is linked with the same small probability. The ring is a Wattsโ€“Strogatz network (WSr in the thesis): a ring of 150 nodes, each linked to its four nearest neighbours on either side, with a tenth of the links rewired at random. โ€œSmall-worldโ€ is a shorthand: depending on how much is rewired, such a network is close to the lattice, a small world, or close to a random network, and most of the rewired networks in the thesis lean towards random. The third is grown by nonlinear preferential attachment (NPA): every newcomer links to earlier nodes with a probability that grows with their degree kk, as kฮฑ+1k^\alpha + 1. With ฮฑ=1\alpha = 1 it resembles a Barabรกsiโ€“Albert scale-free network without being quite one; here ฮฑ=1.8\alpha = 1.8, so one hub is linked to 106 of the other nodes and the next ones to about fifty, on the way to a star (ฮฑ=0\alpha = 0 gives no hubs at all, a large ฮฑ\alpha one giant hub).

So what? Emergence in nature

Boids and your mouse

Starlings and a falcon

Nature is full of the same recipe: simple local rules, on a network of who reacts to whom. Every bird in a flock follows a few local rules: align with your neighbours, stay close, donโ€™t collide. Two birds are easy to understand; the flock as a whole, turning and splitting and closing up again around a threat, is not. The flock on the left is made of boids, the simulated birds of Craig Reynolds (1987), each following just those rules; the click starts it together with the film beside it. Move the mouse over them: they scatter from it and close up behind it. The same happens in the sky over Rome: a starling flock splits around a peregrine falcon (the small dark bird near the gap) and closes up again behind it. Then, in one long shot, the falcon chases the flock, and dark waves run through it as the birds swerve away, one after the other. Nobody in the flock sees the whole; each bird reacts only to the birds nearest to it. (Video: C. Hemelrijk et al. 2019, from the videos to โ€œDamping of waves of agitation in starling flocksโ€, CC BY 4.0; two clips joined and cropped.)

So what? Emergence in nature

Neurons firing

Experiencing and steering reality

A single neuron is well understood: it fires, or not, depending on its inputs. How billions of them together experience reality and steer through it, and even dream, is one of the hardest questions there is. The brain on the left is only an illustration: a neuron fires when enough of its inputs fired, then rests for a moment; the click starts it together with the film beside it. On the right, an octopus fast asleep. Every so often a sleeping octopus slips into a short phase of active sleep, and its skin comes alive: it darkens, breaks into blotches and pales again, running through the same patterns it uses when it is awake, to hide or to signal. That skin is a screen of tiny colour cells, each opened and closed by nerves from the brain, so what we are watching is the brain at work, with nothing outside to react to. The researchers found that its brain activity in these moments looks like waking activity, and compare it to our own REM sleep, the sleep in which we dream. Nobody knows what an octopus experiences; but simple cells, wired together, produce something that looks remarkably like a dream. (Video: an octopus in active sleep, A. Pophale et al., Nature 2023, CC BY 4.0; cropped and trimmed.)

Genotype: the rule

local, simple
what each node does

ฯ•170,1709\phi^{\htmlClass{phi-r}{9}}_{\htmlClass{phi-b}{170},\htmlClass{phi-s}{170}}

Phenotype: the behaviour

global, rich
what the whole does

The research question

What links a simple local rule to the rich behaviour it produces, and can we put that link to use?

Borrowing two words from biology (after S. A. Kauffman, 1990): the rule is the genotype, local and simple, what each node does; the behaviour is the phenotype, global and rich, what the whole does. Here the rule is Replicator, ฯ•170,1709\phi^{9}_{170,170}, on the random network of 150 nodes from before, started with a click: it never settles. The rule fits in two rings, and yet you cannot read what the network does off those rings. On this page, click the rings to change the rule while it runs.

That gap is what the thesis is about, and a Life-like network automaton is a good place to study it: the rule is as simple as a rule gets, the network can be anything, and every step can be watched. Its research question has two halves: what links a simple local rule to the rich behaviour it produces, and can we put that link to use? The rule never tells the whole story, but it tells more than nothing, and the behaviour in turn tells us something about the rule and about the network it runs on.

The thesis in one slide

Part I

Setting the stage

  • 2 cellular automataCNSNS 2025
  • 3 networks

Part II

Rules and their metrics

  • 5 metricsPhysica D 2025
  • 6 synchronisationPhysica D 2025

What links rule and behaviour?

Part III

Machine learning

Part IV

Applications in network science

  • 9 network classificationunder review
  • 10 impact analysiswork in progress

Can we put that link to use?

The thesis has four parts. Part I sets the stage: a taxonomy of the many kinds of cellular automata (A comprehensive taxonomy of cellular automata, CNSNS 2025) and the networks they can live on. Part II is about rules and their metrics: the Lyapunov spectrum, exact growth rates of defects for the simplest rules (under review, and too technical for tonight), and the metrics and the synchronisers of Essential metrics for Life on graphs (Physica D 2025). Part III, skipped tonight, is about cellular automata and neural networks (ACRI 2024 and a Springer book chapter, 2025). Part IV uses rules as instruments in network science: network classification (under review in Information Sciences, first author Lucas C. de Oliveira), and impact analysis, which node matters (unpublished, work in progress).

The two halves of the research question fall on two parts: Part II asks what links a rule and its behaviour, Part IV puts that link to use.

Four chapters make up the rest of the talk, two for each half: getting to know the rules (metrics, synchronisation) and applying them (network classification, impact analysis). They drop to the bottom of the screen and stay there, with the current one filled, so you can always see where we are.

How often does a node switch on?

ฯ•8,129\phi^{\htmlClass{phi-r}{9}}_{\htmlClass{phi-b}{8},\htmlClass{phi-s}{12}}

Hamming weightHW8=0.27\mathrm{HW}_8 = \htmlClass{avg-v}{0.27}

Can we say something about a rule before we run it? Start with one number. Take a node and its eight neighbours, each on or off at random, as in a network that is half on: 512 possible inputs. The ruleโ€™s two rings say which of them turn the node on: born, round a node that is off, and survive, round one that is on. The share of the 512 inputs that turn it on is the Hamming weight, HW8=129โˆ‘sโˆˆ{0,1}โˆ‘q=08(8q)โ€‰ฯ•(s,q/8)\mathrm{HW}_8 = \frac{1}{2^9}\sum_{s\in\{0,1\}}\sum_{q=0}^{8}\binom{8}{q}\,\phi(s, q/8): each slice counts as often as there are neighbourhoods in it, so the slices at the bottom of the rings, around half on, count the most. It is the average next density from half on. For the Game of Life it is 0.27: from half on, about a quarter of the nodes is on one timestep later. (From any other start the same question gives a curve, the average next density or mean-field curve, which Chapter 5 works out too; the talk shows those curves only briefly, in the race of the next chapter.) On this page the buttons switch between the Game of Life and a consensus-seeking rule; the two slides of this chapter share the rule, and its name is under the buttons.

One click brings in the consensus-seeking rule, the kind of rule the next chapter goes looking for: 0.50. It treats on and off alike, so from half on it stays at half on, on average.

How much does one neighbour matter?

ฯ•8,129\phi^{\htmlClass{phi-r}{9}}_{\htmlClass{phi-b}{8},\htmlClass{phi-s}{12}}

Our result Metrics for LLNAs, read off the rule before it runs

Boolean sensitivityBS8=3.17\mathrm{BS}_8 = \htmlClass{avg-v}{3.17}

The second number is about differences. Take the same node and its eight neighbours, and flip one of its nine inputs. Flip a neighbour, and the share of neighbours on moves one slice along the ring, one way or the other: the red arrows. Flip the node itself, and it moves to the other ring. Does the nodeโ€™s next state change? Count the inputs whose flip changes it, averaged over all 512 inputs: that is the Boolean sensitivity, BS8=129โˆ‘sโˆˆ{0,1}โˆ‘q=08(8q)(IS(s,q,8)+NS(s,q,8))\mathrm{BS}_8 = \frac{1}{2^9}\sum_{s\in\{0,1\}}\sum_{q=0}^{8}\binom{8}{q}\big(\mathrm{IS}(s,q,8)+\mathrm{NS}(s,q,8)\big), split into the node itself (IS) and its neighbours (NS), at most 9. For the Game of Life it is 3.17: of a nodeโ€™s nine inputs, about three matter. It is also how many nodes a single difference reaches in one timestep, on average: above one, a difference spreads. (Chapter 5 follows a difference of any size too, the average next defect or Derrida curve; its slope at the origin is this number.)

The consensus-seeking rule: 6.40, twice as sensitive as the Game of Life. A single difference spreads fast under this rule, and that will matter in the next chapter.

Our result: we worked out these numbers, and the two curves behind them, for Life-like network automata, so that they can be read off any of the 262 144 rules before running it. They are what the filter of the next slide is made of.

The mission: looking for consensus

262 144
512
27

Strong filter:
metrics inform behaviour

symmetry

simple thresholdthe winner

Can we really use what a rule looks like to predict what it will do? A test, in the form of a mission: looking for consensus. From a random start, after a while every node should have the same colour, all on or all off. The obvious candidates are threshold rules: do what most of your neighbours do. On a network they get stuck, in patches that each agree with themselves but never with each other.

There are 262โ€ฏ144 rules at resolution 9: every way of colouring the two rings, here the Game of Life. First filter: we do not care which colour wins, only that both are equally likely, so the rule must treat on and off alike: symmetry, it must be its own cousin, like the simple threshold rule shown here (the button says so). That leaves 512. Second filter, on the average next density: from half on and half off, the network must move away, towards all on or all off. That leaves 27, about one rule in ten thousand, found in a fraction of a second without a single simulation. Here each of them stands where it will run, with the two curves behind the numbers of before, just for a glance (the thesis has the details): in blue the average next density from any start, in red the average next defect, each with its tangent. Every blue curve leaves one half more steeply than the diagonal, so from half and half the network is pushed away; the red tangent is the Boolean sensitivity.

Now we can afford to run every one of them. One click puts them in place, the next runs them: all 27 on the same grid of 900 nodes, a fifth of its links rewired, from the same random start, for at most 100 timesteps, as in the thesis. Most of them freeze in patches within a few rounds, or keep churning, and fade; the simple threshold rule is stuck from timestep 10. Three reach agreement, and they are the only ones: of all 512 rules that treat on and off alike, no other reaches a lasting agreement in this race. A strong filter: the metrics inform the behaviour, and a run decides which rule is best. Run over hundreds of networks, one of them clearly is: phi488,4649phi^9_{488,464}, which brought more than 96% of the rewired grids of the thesis to agreement. It is among the most sensitive of the 27, but sensitivity alone is not enough (Chapter 6, Physica D 2025).

Watch it work

Our result Filter for rules that do what we want

The winner at work, on the grid of the race (900 nodes, eight neighbours each, a fifth of the links rewired), from the same random start with half of the nodes on. Patches form and fight, and from timestep 31 every node is on. The curve on the right is the share of nodes that are on, timestep by timestep. The buttons switch to the simple threshold rule, on the same network from the same start: within a dozen timesteps it freezes into patches, each content with its own majority, a few nodes blinking at the borders.

The slider under the grid rewires it: none at the left end, then on a logarithmic scale from one link in a thousand to all of them; back at 20% it is the raceโ€™s grid again. Moved while the rule runs, it only swaps the links: every node keeps its colour and the run goes on (one that was stuck may come loose). The talk stays there: on rewired grids like this one, the winner brought more than 96% of the networks of the thesis to agreement at a fifth rewired, but only 21% of the grids without rewiring and 37% of the fully rewired ones. On this grid and from this start it agrees within the thesisโ€™s 100 timesteps from about 2% to 50% rewired, and not below or above. The thesis ran 900 samples for every setting, on rewired rings, rewired grids and scale-free networks; these grids are a small illustration.

Our result: with the metrics as a filter, we find rules that do what we want before running a single one. (To say out loud, not on the slide: which raises the question the other way round. Can we start from rules that already do something we want, such as telling networks apart, and find out what makes them work?)

Does behaviour predict topology?

Hunting for distinct fingerprints

random
small-world ring
scale-free

The first application turns the question around: instead of asking what a rule does on a network, ask what it tells us about the network. Here one rule runs on the three networks of before, 150 nodes each, from a random start with three nodes in ten on: first the consensus-seeking rule of the mission. Under each network, three empty lines. The click starts it, for at most 100 timesteps, as in the thesis.

Under each network, three histograms sum up what its nodes did over the run: how evenly each node split its time between on and off, how irregular its row of ons and offs was, and how many of its neighbours were on (the Shannon entropy and the Lempelโ€“Ziv complexity of Miranda et al., and the neighbourhood density of Chapter 9). Together they are the networkโ€™s fingerprint, drawn over them as a real one, and a simple classifier can learn to tell such fingerprints apart. Under the consensus-seeking rule every network agrees on all off within about ten timesteps, and from then on every node does the same: three fingerprints alike, and the run stops a second later. A bad detective.

The next click switches to a good detective, ฯ•170,489\phi^{9}_{170,48}, from a new start, with the lines empty again; the click after runs it. Now the histograms differ, and so do the fingerprints: a whorl, a loop and a tented arch. On the ring every node keeps switching, irregularly, on about half of the time; on the random network about one node in five goes off for good (a node with two or four links can never be born under this rule); on the scale-free network more than half of the nodes never switch on again, and the rest switch in a more regular rhythm. So we look for rules whose fingerprints are clearly different. (G. H. B. Miranda and colleagues started classifying networks this way in 2016; the thesis did it on 23 datasets, with histograms split by degree and a support-vector machine. On the page the buttons switch the rule and โ†บ draws a new start; about one start in fifteen leaves the ring stuck under the consensus-seeking rule.)

What makes a good detective?

A consensus-seeking rule ฯ•488,4649\phi^{9}_{488,464}

6 of 16

A good detective ฯ•170,489\phi^{9}_{170,48}

10 of 16

Our result Good detectives change their mind often, but not always

What makes the difference can be read off the rule itself. Back to the rings, and to the red arrow of the sensitivity slide: one neighbour more or fewer moves a node one slice along its ring. The next click draws that arrow wherever it changes the answer, between every two neighbouring slices where the rule says something else. The consensus-seeking rule changes its mind at 6 of the 16 steps, the good detective at 10: a rule that changes its mind often is a sensitive rule, since a small change in the neighbourhood often changes what a node does.

This is our contribution to network classification, Chapter 9, work led by Lucas C. de Oliveira (with M. Rollier, J. M. Baetens and O. M. Bruno, under review). Over thousands of rules and 23 datasets, the best detectives change their mind often, but not always: at between half and nine in ten of the steps. That makes sense: such rules are sensitive enough to pick up how networks differ, but not so sensitive that it drowns in chaos; they behave in a complex way (a reading the thesis proposes but does not test). Our detective, at 10 of 16, sits inside that band. On 13 of the 23 datasets, more often tends to classify better (correlations between about 0.2 and 0.46), and searching among the rules that change their mind often, but not always, finds a rule within 5% of the best while trying only 10 to 20% of the rules. On social and bioinformatics networks there is no such pattern, but the method itself is mediocre there as well. (The count is what the thesis calls the jaggedness: the neighbour part of the Boolean sensitivity with every slice weighed the same, whatever the network. The Boolean sensitivity itself shows the same pattern, slightly weaker.)

What could this be good for, when nobody has the network? Think of a contact-tracing app such as Coronalert, launched in Belgium six years ago this week: phones sense which phones come near and, for privacy, nobody ever holds the network of who met whom. If every phone ran one simple rule with the phones around it and reported only its own ons and offs, their fingerprint could still say what kind of network we live in, for instance whether our contacts run through a few superspreaders. We cannot try thousands of rules on real people, so the rule would be chosen in advance, from simulations alone: one that changes its mind often, but not always. (A thought experiment; a backup slide shows it.) And the consensus-seeking rules, such bad detectives, turn out to be good at something else: finding who matters.

The stubborn node

Which individuals steer the group?

  1. Find a converging system
  2. Force a node to be stubborn

    Some are not strong enough

    Others convince the entire group

The second application uses the consensus-seeking rules that were such bad detectives. Here a toy village of 100 people, about eight contacts each: which individuals steer the group?

First, find a converging system: everyone follows a consensus-seeking rule, from a start with half of them on. The click runs it: with nobody stubborn, the village ends with every node on.

Then force a node to be stubborn: it stays off, no matter what. Some are not strong enough: when a quiet person at the edge is stubborn (ringed), nobody else follows. Others convince the entire group: when an important person is stubborn, the whole village flips to all off. Holding one node fixed like this is called clamping, or pinning: it is an intervention, not just an observation. (An illustration, from the same start every time.)

importance: low high

Importance score

Analysing the impact of individuals

Example: a health campaign

Budget: three ambassadors

Real rule unknown: use a stand-in

our ambassadorsmost contacts

Our result Rank who steers the group, with a stand-in rule

Now make every node stubborn in turn, one at a time, and record how much of the village ends up following it: that is each nodeโ€™s dynamical importance score, the darker, the more (in the thesis, holding a node fixed is called clamping). In this village most people sway only a small part of it, and the most important about half of it, on average. (Chapter 10, work in progress; an illustration over 100 starts that reach agreement, with each node stubborn in each.)

Why would we want to know who matters? Take an example: a health campaign with a small budget. Marketing campaigns and political propaganda ask the same question; the health campaign is the friendlier example. In 32 villages in Honduras, researchers started a water-purification and a vitamin programme with a few villagers each, to see how far it would spread, and compared ways of picking those few (Kim and colleagues, The Lancet, 2015). So which three villagers should be our ambassadors? We do not know how villagers really persuade each other, so we use a stand-in: the consensus-seeking rule, not to model the villagers, but to measure them. The three most important, ringed in red, become our ambassadors: an ambassador is a stubborn node, who keeps promoting the programme, no matter what. The dashed rings are the three villagers with the most contacts, 21 each: not the same people. Importance goes with the number of contacts (a rank correlation of 0.82 here), but not one-to-one: the many contacts of those three also listen to the centre of the village, and follow it. In the Honduran trial too, picking the villagers with the most friends did no better than picking at random, while picking friends that villagers named did better, at least for the vitamins. (An example, not a campaign I studied.)

Our result: a ranking of who steers the group that needs no knowledge of the real rule, only a stand-in rule and a simulator. Whether the stand-in picks the same people as the real rule would is the most important next step of Chapter 10: comparing it with, say, an epidemic model on the same networks. Taking the top three one by one also ignores that they may reach the same people; choosing them together is a next step. This part of the thesis is unpublished and still in progress.

Summary

Getting to know the rules

metrics: define and quantify LLNAs
synchronisation: filter for rules that do what we want

Applying them

network classification: understand rules that generate distinct fingerprints
impact analysis: assign a dynamical importance score

To sum up, along the two halves of the research question, and the chapters along the bottom of the slides.

Getting to know the rules. Metrics: we defined Life-like network automata so that every rule has its cousin, and worked out numbers and curves that quantify a rule before it runs, such as the Hamming weight and the Boolean sensitivity, here for one rule after another (Chapter 5).

Synchronisation: with those metrics as a filter, we find rules that do what we want, here a consensus-seeking rule that brings a rewired grid to agreement, from one random start after another (Chapter 6).

Applying them. Network classification: we understand which rules give networks distinct fingerprints, good detectives that change their mind often, but not always, here on the three networks (Chapter 9).

And impact analysis: making every node stubborn in turn gives each a dynamical importance score. Here the best-connected villager, the most important one, a quiet one and another important one, in turn: twice the village follows, twice it does not, each time what that person does in most starts (Chapter 10, work in progress).

Beyond this thesis

BionamiXSpatiotemporal modelling

Groundwater level Jan

  • lowhigh
  • well

The MACCA-GW model ยท ACRI 2026, accepted

Wildfires

  • tree
  • on fire
  • burnt
Influence maximisationMany open questions, many relevant applications
Inside artificial mindsThe same gap, with very high stakes

Where does this go next? First, at BionamiX: spatiotemporal modelling, where cellular automata model processes in space and time. Those need the extensions from the taxonomy of Chapter 2: rules that differ from cell to cell, more states, larger neighbourhoods.

Groundwater level: the MACCA-GW model, a cellular automaton for groundwater flow, applied to De Wilde Landen in the Netherlands, is accepted for the ACRI 2026 proceedings (with L. Landuyt and J. M. Baetens). It is not part of the thesis: a model in space, not a network automaton. Every cell holds the height of the water table; water flows to it from its four neighbours, rain falls, the land cover draws some back, and a line of wells pumps, more in summer. The map is a simplified version on a made-up reserve, not the calibrated model, started with this click: a year in ten seconds, the water table highest in spring and lowest in late summer, deepest along the wells. Next to established models such as MODFLOW it offers simplicity, and many runs at once for ensembles; in the paper, pumping more in winter and less in summer raises the spring water table by about a third of a metre, at no cost.

And wildfires, close to home: in August a fire in the Hautes Fagnes burnt some 3,000 hectares, one of the largest in Belgian history, and villages had to be evacuated. The picture is a wildfire as a simple cellular automaton, an illustration and not the BionamiX model, started with this click: green trees, red fire, charcoal where it has passed, a wind from the left; a burning tree burns for several timesteps and can set each of its eight neighbours alight, most easily downwind. Models like this ask where a fire goes next, and where a firebreak would help most.

Second: influence maximisation. I wish I had started on node importance earlier: it holds many open questions and many relevant applications. Which few people, chosen together, reach the most? Does a stand-in rule pick the same people as the real processes on a network, an epidemic, an opinion, a product spreading? The claims of the impact chapter still have to be checked against such processes, and on real networks. The picture, started with this click, is the village of before with our three ambassadors (ringed) holding on to the programme: from one start after another, the village follows them, nearly always.

And next for me: what goes on inside artificial minds. A neural network is also simple operations, iterated: every node adds up what comes in, each input with its weight, adds a bias, and passes the sum on only if it is positive (a ReLU). Here a small one with random weights and biases, started with this click, processes one input after another: each link lights up with what it passes on, a node fills as it fires, and the largest output is ringed. Nobody writes down what a trained network does: the same gap between rule and behaviour, with very high stakes. The results of this thesis do not carry over directly, but the way of thinking, from the parts to the whole and back, can help.

The rich landscapewhat we see

of iterated simplicitywhat we write down

We brought these two worlds closer together.

Thank you

michielrollier.be

every demo of tonight, and the thesis chapter by chapter

Ghent UniversityFaculty of Bioscience EngineeringUniversity of Sรฃo PauloBionamiX

The title once more, in two halves: the rich landscape is what we see, the behaviour; iterated simplicity is what we write down, the rule. This thesis brought those two worlds a little closer together.

Thank you. The colours drain out of the landscape, as on the still-life slide: every cell simply alive or dead, and nothing changes. The QR code leads to this site, with every demo of tonight and the thesis chapter by chapter.

The last click takes away one cell, and the still life breaks down, as on the still-life slide: the debris spreads over the landscape and keeps going while the slide is up. Where it comes near the words and the logos it is faded, so they stay readable.

What single question unifies the thesis?

  • How much does a simple local rule tell you about the global behaviour it produces, and what does the behaviour tell you about the rule and the network: genotype against phenotype.
  • Three objectives: rigorous foundations for these systems; how network topology shapes the dynamics; the systems turned into instruments, for neural networks and for network science.
  • In the thesisโ€™s words: to map, understand and apply the link between iterating simple local rules and the rich landscape of emergent outcomes.
  • Chapters 4 to 10 each answer one part; Chapter 11 puts the answers side by side.

Ch. 1 ยง1.3, Ch. 11 ยท back to the questions

What are the main limitations?

  • The theory is shallow: exact for a small corner of rules, one step deep elsewhere; the rest is simulation.
  • Narrow settings, broad conclusions: one resolution, one network family, one starting density, synthetic networks only.
  • Evaluation lagged behind generation: point estimates without error bars, no rule held out, an importance score whose reliability went unmeasured.
  • Mechanisms are proposed, not shown: that sensitive rules synchronise and jagged rules discriminate is a reading of a pattern.

Ch. 11 ยง11.3 ยท back to the questions

Why these networks and these parameter ranges?

  • A Life-like rule only reads the share of a nodeโ€™s neighbours that are on, so what matters is how the degrees are spread and how the links are arranged.
  • Each family is a dial for one of those: rewiring pp runs a ring or grid from clustered and long to random and short; the power ฮฑ\alpha runs preferential attachment from no hubs to one hub.
  • At p=0p = 0 the rings and grids are the neighbourhoods of classical cellular automata, which anchors each dial in known territory.
  • Synthetic networks only: empirical ones are the next step.

Ch. 3 ยท back to the questions

Is your Lyapunov exponent really an entropy?

  • It is a tangent-space quantity: the growth rate of a count of walks along which a defect can travel, which makes it entropy-like by construction.
  • For the parity rule on a graph the exponents are lnโกโˆฃฮปk(A)+aโˆ˜โˆฃ\ln\lvert\lambda_k(A) + a_\circ\rvert, the maximal one lnโกฯ(A+aโˆ˜I)\ln \rho(A + a_\circ I); for the rule without self, the logarithms of the adjacency eigenvalues.
  • The thesis itself calls it closer to an entropy than to a rate of separation.
  • Wolfram (1984) introduced both notions for cellular automata side by side.

Ch. 4 ยท back to the questions

Without a chaos claim, why keep the name โ€œLyapunov spectrumโ€?

  • For continuity with the literature; the thesis says plainly that it counts walks and claims no chaos.
  • What the chapter adds: exact spectra as benchmarks, which correct spectra published earlier.
  • And the Boolean Jacobian, whose expected row sum is the sensitivity that Chs. 5 and 6 use to read rules.
  • Still, the chapter stands somewhat apart from the rest of the thesis.

Ch. 4 ยง4.4 ยท back to the questions

Why odd resolutions and mirrored intervals?

  • Complementation symmetry: swapping on and off must map the rule space onto itself.
  • That needs the cut of the density axis mirrored around 12\tfrac12, with a middle interval that contains 12\tfrac12, hence an odd number of intervals.

Ch. 5 ยง5.2, App. C ยท back to the questions

How good is the mean-field approximation?

  • It is one step deep and annealed: every round, contacts are redrawn at random.
  • It ignores correlations and clustering, which the real network builds up at once.
  • Hence simulation remains necessary; the genotype narrows the search, it does not replace it.

Ch. 5 ยง5.3 ยท back to the questions

Is the genotype filter really cheaper than a genetic algorithm?

  • The genetic algorithm of Das et al. spends some 10710^7 to 10810^8 simulations on the search.
  • The filter simulates nothing: a bitwise condition and 512 mean-field curves, in a fraction of a second.
  • Verification costs both the same order: the saving is in the search.
  • The evolved rules of Das et al. synchronise better, on lattices: ฯ•23,479\phi^9_{23,47} is not the best synchroniser available, and the rule spaces differ.

Ch. 6 ยง6.1.3 ยท back to the questions

Does small-world structure cause synchronisation?

  • No causal claim. The rule fails for pโ‰ค0.04p \le 0.04, peaks near pโ‰ˆ0.3p \approx 0.3 (94%), beyond the small-world band, and drops to 39% at p=1p = 1.
  • Short paths look necessary but not sufficient, for this rule.
  • Rewired grids behave differently: 82% unrewired, above 95% for 0.02โ‰คpโ‰ค0.40.02 \le p \le 0.4.
  • At p=0.2p = 0.2 the thesis reports 83% from the screening and 78% from the sweep, both at half the nodes on: two estimates of one point.
  • pp is not the same amount of rewiring in both families: a ring changes 2pโˆ’p22p - p^2 of its links, a grid pp. The text now says so; the simulations stand, only the comparison between the two pp axes shifts.

Ch. 3 ยง3.3.2, Ch. 6 ยง6.1.3 ยท back to the questions

Can the networkโ€™s spectrum predict where the synchroniser works, as for oscillators?

  • Not studied.
  • The spectrum governs the dynamics exactly only for affine rules such as parity; ฯ•23,479\phi^9_{23,47} is not affine.
  • The closest bridge is pinning control, which Ch. 11 names as a direction.

Chs. 4, 6, 11 ยท back to the questions

Can a classifier discover behaviour classes, or only transfer known ones?

  • Its output only transfers the labels it was trained on.
  • Its internal representation can support exploration without labels, but that is not discovery.
  • A self-supervised attempt turned out hard.

Ch. 8 ยท back to the questions

How much confidence do the 98% and 61% deserve?

  • The test set was kept apart for the final numbers, but the validation set chose the augmentations.
  • Each entry is the best of a few runs, as its caption says; there are no error bars.
  • No rule was held out, so 61% measures recognising rules the network has seen.
  • The lesson does not depend on the decimals: first check whether scanning the rule table already solves the task.

Ch. 8 ยท back to the questions

What makes the best classification rules special? The edge of chaos?

  • The analysis works at the level of partitions of the rule space by jaggedness, not rule by rule; the best individual rules were not examined one by one.
  • Calling the rules in the sweet spot โ€œcomplexโ€, or placing them at an edge of chaos, is a hypothesis; no measure of their behaviour was computed to test it.
  • The edge of chaos is contested anyway: Mitchell, Hraber and Crutchfield (1993) could not reproduce Packardโ€™s 1988 result.
  • Part of it may be counting: partitions are largest near a jaggedness of 12\tfrac12, and the best of many rules looks better than the best of few.

Ch. 9 ยง9.4.2, ยง9.6.3, ยง9.7 ยท back to the questions

Why not graph neural networks or graph kernels for classification?

  • Beating them was not the aim; on social and chemical benchmarks they do better.
  • The rule-based pipeline needs no node attributes, handles networks of any size, and comes with one interpretable parameter: the rule.

Ch. 9 ยง9.7 ยท back to the questions

What did the joint doctorate with Sรฃo Paulo bring?

  • Network automata as instruments: the pipeline that classifies networks by the patterns a rule draws on them comes from the Sรฃo Paulo group.
  • In Ch. 9, first-authored by Lucas de Oliveira, my part is the link from jaggedness to the genotype measures of Ch. 5, and review and editing.
  • Ch. 10, on node importance, grew out of the same collaboration.

Chs. 9, 10 ยท back to the questions

Isnโ€™t the stubborn person just a zealot from the voter model?

  • Related, and the thesis says so: a zealot never changes its mind, and neither does a clamped node.
  • But the aim differs: not a model of opinions, a probe of how much each node matters.
  • Every node is clamped in turn: an intervention, not one zealot observed.
  • The rule is deterministic, non-monotone and chosen by its genotype. A quantitative comparison with zealot results is not in the thesis; it would be the natural test.

Ch. 3 ยง3.5, Ch. 10 ยท back to the questions

Clamping is an intervention; what does merely swapping states do?

  • Swap the most central person who is up with the least central one who is down, again and again.
  • On preferential-attachment networks with 8 contacts: at most 10 swaps by closeness flip the outcome 90% of the time.
  • Random swaps only push the outcome towards a coin toss: the 90% mark is never reached.
  • That random reference covers one setting only, ฮฑ=0\alpha = 0, reconstructed from saved, binned output: thin, and a full rerun is the next step.
  • Clamping needs no such control: it tests every node, blind to any centrality.

Ch. 10 ยง10.3.1 ยท back to the questions

Is dynamical importance confounded by conditioning on consensus?

  • ฮท\eta pools clampings that change a nodeโ€™s start (its hand was up) with ones that only confirm it (already down).
  • Which starts reach consensus is itself linked to centrality; the thesis does not separate the two.
  • On the toy network of tonight, the Spearman correlation of ฮท\eta and degree is 0.820.82 overall, 0.820.82 from the changing clampings only, 0.570.57 from the confirming ones only.
  • So the ranking survives the split here, on one toy network; that is an illustration, not a result.
  • Why not split the thesisโ€™s own scores? The raw clamping outputs are no longer available, so it needs a rerun; it is listed as the next step.

Ch. 10 ยง10.2.2, ยง10.4, Tab. 11.1 ยท back to the questions

Could the score find keystone species, or whom to vaccinate?

  • It needs a known network and a dynamics that can be simulated on it.
  • The automaton is a generic probe, not a model of the real process: whether its ranking carries over is exactly the untested step.
  • Choosing a set of nodes rather than ranking single ones is influence maximisation, the direction Ch. 11 proposes.

Chs. 10, 11 ยท back to the questions

Why discrete automata, and what is lost against continuous models?

  • Gained: exact computation with no discretisation error, and networks, where there is no continuum to discretise.
  • Lost: calibration to physical parameters and built-in conservation laws.
  • The groundwater model sits right on that trade-off.

Ch. 1, MACCA-GW ยท back to the questions

Why is the groundwater work not in the thesis?

  • It is a spatially explicit cellular-automaton model of groundwater flow, not a network automaton.
  • It does not fit the thesisโ€™s story from rules to networks, so it is published separately, in the ACRI 2026 proceedings.
  • Next to MODFLOW it offers simplicity and many runs at once for ensembles; not parallelism as such, since parallel MODFLOW 6 exists.

MACCA-GW ยท back to the questions

With one more year, what would you do, and what differently?

  • Finish Ch. 10 properly: faster simulation, the two kinds of clamping split, the scoreโ€™s reliability measured, a comparison with an epidemic or Ising dynamics, empirical networks.
  • Differently from day one: held-out test sets, and every raw output stored.

Ch. 10 ยง10.4, Ch. 11 ยท back to the questions

What do cellular automata offer interpretability research?

  • Circuits whose weights are known exactly: a cellular automaton written into a neural network by hand (Ch. 7).
  • The shortcut lesson: high accuracy that came from reading the rule off the first row, not from the behaviour (Ch. 8).
  • An idea for future work, not a result of the thesis.
  • It is also a deliberate change of field, not only a bridge.

Chs. 7, 8 ยท back to the questions

Who meets whom?

Identifying a network nobody can see

a few hubs:
superspreaders

The rule to load: one that changes its mind often, but not always.

Ch. 9 ยท back to the questions

What could this be good for, when nobody has the network? Six years ago this week, on 30 September 2020, Belgium launched its contact-tracing app, Coronalert: every phone sensed over Bluetooth which phones came near it, and to protect our privacy nobody ever held the network of who met whom. Now imagine every phone running one simple rule with the phones around it, and reporting only its own ons and offs, never whom it met. The next click starts them: a small network of 48 phones, hidden.

From those rows of ons and offs alone, the fingerprint says what kind of network we live in: whether our contacts run through a few hubs, superspreaders, or are well mixed. Here the hidden network is scale-free, with one phone linked to 34 others (its links appear on this click), and the phones run the good detective, so the fingerprint is the scale-free one of before. A thought experiment, of course: real contacts change from hour to hour.

But which rule do we load into millions of phones? We cannot try thousands of rules on real people. This is where knowing what makes a good detective pays off: a rule that changes its mind often, but not always, chosen in advance from simulations alone. And the consensus-seeking rules, such bad detectives, turn out to be good at something else: finding who matters.

fullscreen ยท or press the clicker to start in this window