Mind Frontier · No. 32

Dynamical systems for cognition: attractors, bifurcations and metastability

Rajesh Ranjan Mahapatra

Attractors, hysteresis, critical slowing down, metastability and criticality explain conscious state transitions only when each predicts a signature that a rival account does not, and each signature comes with a control that can defeat it.

© 2026 Rajesh Ranjan Mahapatra. 74 min read

Rajesh Ranjan Mahapatra

Dynamical systems for cognition: attractors, bifurcations and metastability

Attractors, hysteresis, critical slowing down, metastability and criticality explain conscious state transitions only when each predicts a signature that a rival account does not, and each signature comes with a control that can defeat it.

TLDR

  1. A dynamical concept explains a conscious state transition only when it predicts a measurable signature that a rival account does not, such as a hysteresis loop that survives slow change.
  2. Lag mimics hysteresis: in my sweep-rate test a slow single-state system showed a wider loop than a true two-state switch at fast sweeps, and its loop fell to zero when the sweep slowed.
  3. Rising autocorrelation is a weak warning: changing noise produced a stronger warning than a real fold in every one of 40 runs, while large-noise jumps tipped with no warning.
  4. Metastability and criticality need controls: an uncoupled population gives half the peak metastability index, and subsampling to one per cent made a near-critical process read 0.20 instead of 0.98.
  5. Report each dynamical claim with its signature and its control: a slower sweep, a steady-noise check, an uncoupled population or a subsampling-corrected estimator.
Opening scene

The patient who woke at a lower dose#

Invented scenario Everything in this story is invented. It is here to give intuition, never evidence.

Dr Kabir Menon is three months into his anaesthesia training at a teaching hospital in Bengaluru. His patient this morning is Joseph D'Souza, a retired bus conductor of 58, who is having a knee replaced. The drug is propofol, given by a pump that estimates how much of it has reached the brain. The pump shows that estimate as a number, the effect-site concentration, in micrograms per millilitre.

Kabir raises the target slowly, in small steps, and asks Mr D'Souza to squeeze his hand after each one. At 2.5 he squeezes. At 2.8 he squeezes, more slowly. At 3.0 he does not. Kabir writes the number down.

Two hours later the operation is over, and Kabir lets the concentration fall. He asks again for a squeeze at each step. At 2.5 nothing happens. At 2.0 nothing. At 1.4 the hand closes, and Mr D'Souza opens his eyes and asks whether they have started.

Kabir shows the two numbers to his consultant, Dr Anjali Rao. He lost response at 3.0 and found it again at 1.4. Going down, he stayed under well past the point where he went under on the way up.

"His brain has two stable states," Kabir says. "Awake and anaesthetised. It is an attractor system with hysteresis. I read about it last week. It is called neural inertia."

Anjali has run this hospital's theatre lists for eleven years. She likes the idea and she does not trust the word. "Or the pump's estimate of what is in his brain is wrong," she says. "The number on the screen is a model of an average patient. If the drug reaches his brain more slowly than the model thinks, you would see the same gap with no attractor at all. What would tell the two stories apart?"

That question is the subject of this paper. Dynamical systems theory gives mind science a rich vocabulary: attractors, bifurcations, hysteresis, critical slowing down, metastability and criticality. Each word names something exact in mathematics. In papers about the brain and consciousness, each word is also used loosely, as if naming a pattern explained it.

Here is the answer in plain words, before the argument. A dynamical concept carries explanatory weight only when it predicts a signature that a rival explanation does not. An attractor predicts that a perturbed state returns, and that time spent in a state follows a particular law. A bifurcation with hysteresis predicts a gap between the way in and the way out that does not shrink when the change is made slowly. Critical slowing down predicts that recovery from small disturbances gets slower before a transition, but the same slowing can appear for other reasons. Metastability and criticality are the hardest to test, because the measurements that seem to show them are easy to produce in other ways.

The consequence for Kabir is practical. Neural inertia is a real hypothesis with real evidence behind it, but his two numbers do not test it. The consequence for mind science is broader. The vocabulary of dynamics should be used as a set of predictions, and each prediction comes with a control that could defeat it.

[AUTHOR NOTE: if you have watched an induction or emergence from anaesthesia, or have been anaesthetised yourself, a two-sentence first-person account of what the transition looked or felt like would ground this scene.]

The paper builds the case in steps. Part 1 says what an attractor is and what it predicts. Part 2 turns to bifurcations and hysteresis, and builds an instrument that separates a true bistable switch from a slow lag. Part 3 tests the early-warning signals that are meant to precede a transition. Parts 4 and 5 take on metastability and criticality, the two concepts most often linked with consciousness. Part 6 asks what the Kashmiri Śaiva doctrine of spanda, a pulse at the root of every state, can and cannot add.

The running example: one gap, two explanationsDiagram of the invented running example: a patient lost response at 3.0 and regained it at 1.4; a trainee proposes two stable states, a consultant proposes drug lag.Theatre four: a patient who woke at a lower dose than he went underproposesproposesJoseph D'SouzaKnee replacement under propofol.Lost response at 3.0, regained itat 1.4Dr Kabir MenonTrainee. "Two attractors: neuralinertia"Dr Anjali RaoConsultant. "Or the pump's estimateis wrong"Story 1: two stable statesThe gap stays when the drug is changed slowlyStory 2: drug lagThe gap shrinks to zero when the drug is changedslowlyNumbers are effect-site concentrations estimated by the pump, in micrograms per millilitre.
Figure 1. The running example: one gap, two explanations. The invented scenario that runs through the paper. The patient sits between the trainee and the consultant; their two stories make different predictions about slow change. What to notice: the two numbers alone fit both stories. Illustrative (invented scenario)
Text description of figure 1

Top row: Dr Kabir Menon, trainee, says two attractors and neural inertia; Joseph D'Souza, the patient, lost response at an estimated 3.0 and regained it at 1.4; Dr Anjali Rao, consultant, says the pump's estimate may be wrong. Bottom row: story 1, two stable states, predicts the gap stays under slow change; story 2, drug lag, predicts it shrinks to zero.

Part 1

What an attractor is, and what it predicts#

Start with a ball in a bowl. Push it and it rolls back to the bottom. Push it harder, over the rim, and it rolls into the next bowl instead. That picture is most of what people mean when they say a brain state is an attractor.

The mathematics makes the picture exact. A dynamical system is a rule that says how a state changes from one moment to the next. The state can be anything that can be written as a list of numbers: the positions of a ball, the firing rates of a thousand neurons or the concentrations of chemicals in a cell. The set of all possible states is called the state space. The rule turns every point in that space into a direction of motion, so the whole space is filled with arrows, like wind on a weather map.

An attractor is a region of the state space that nearby states flow into and do not leave. The simplest is a fixed point, a single state where all motion stops, like the bottom of the bowl. The set of states that flow into a given attractor is its basin. A system with two attractors has two basins, divided by a boundary. On the boundary sits a balance point, unstable in the way a pencil on its tip is unstable.

The double well#

The standard model of a two-state system is called the double well. Its state is one number, x. The rule says x moves downhill on a curve with two valleys, one at x = −1 and one at x = +1, separated by a hill at x = 0. Add random kicks, noise, and the state jiggles around in one valley. Occasionally a run of kicks in the same direction carries it over the hill, and it settles in the other valley.

Two predictions follow. The first concerns how long the state stays in one valley before it hops. In 1940 the Dutch physicist Hendrik Kramers worked out the answer for a particle jostled by heat (Kramers 1940). The mean time spent in a valley grows exponentially with the height of the hill divided by the strength of the noise. Halve the noise and the wait does not double. It squares, roughly.

The second prediction concerns the shape of the waiting times. Each hop is set off by a lucky run of kicks, and the system has no memory of how long it has waited. That makes the waiting times follow an exponential distribution, the same law as the time between radioactive decays. Short waits are the most common and long ones fall off steadily. One mark of an exponential is that the standard deviation of the waits equals their mean.

An attractor pair is two valleys and a quick crossingTwo panels: a curve with two valleys at minus one and plus one and a hill at zero, and a time series that hops between the two valleys.A noisy double well: two valleys, and a state that hops between themvalley at −1valley at +1hill: height ¼State xPotential-10+1Time (6,000 steps of 0.1 units, noise 0.45)Simulated with a fixed seed. The state spends long spells near −1 or +1 and crosses quickly.
Figure 2. An attractor pair is two valleys and a quick crossing. Left: the shape of the double-well potential used throughout Parts 1 to 3. Right: one simulated run. Notice that the state lingers near each valley and spends little time on the hill, which is the third signature in Part 1. Simulated (fixed seed)
Text description of figure 2

Left panel: potential curve with valleys at x = −1 and x = +1, separated by a hill of height one quarter at x = 0; a ball sits in the left valley. Right panel: a simulated trajectory over 6,000 steps at noise 0.45, with long spells near −1 and +1 and quick transitions.

A minimal worked case: measuring the hill#

I ran the double well with six levels of noise, using fixed seeds so the runs can be repeated. For each level I recorded how long the state stayed in a valley before crossing to the other. At a noise strength of 0.35 the mean stay was 283 time units. At 0.5 it was 36. At 0.6 it was 19.

Kramers' law says the logarithm of the mean stay should rise in a straight line against one over the noise strength squared. The slope of that line should be twice the height of the hill, which in this model is a quarter. So the slope should be 0.5. The fitted slope from my runs was 0.49. At the noise level of 0.45, I also checked the shape of the waits. Over 124 stays the standard deviation was 1.17 times the mean, close to the value of one that an exponential predicts.

Waiting times follow the law the attractor predictsLine chart: the log of mean dwell time rises linearly with one over noise squared, with fitted slope 0.49 against a predicted 0.5.Kramers' law: log of mean time in a valley against one over noise squared34567834561 / (noise strength)²Natural log of mean staySimulated mean staySlope 0.5 predictedby theoryfitted slope 0.49
Figure 3. Waiting times follow the law the attractor predicts. My run, five noise levels from 0.35 to 0.6, fixed seeds. Mean stays ranged from 283 time units at noise 0.35 to 19 at noise 0.6. Notice that the points lie on a straight line whose slope, 0.49, matches the 0.5 fixed in advance by the height of the hill. Measured on synthetic data
Text description of figure 3

Five points: at 1/s² of 2.8, 4.0, 4.9, 6.25 and 8.2 the log of mean stay is about 2.9, 3.6, 4.0, 4.6 and 5.6. A dashed line of slope 0.5 runs through them. The fitted slope is 0.49.

Stays in a valley are memorylessHistogram of 124 dwell times in bins of 20 time units, highest in the first bin and falling steadily.Distribution of 124 stays in a valley at noise 0.45Number of stays390–202320–401340–601660–80880–10010100–1202120–1404140–1604160–1800180–2001200–2200220–240Mean stay 64 time units; standard deviation 1.17 times the mean. An exponential law gives exactly 1.
Figure 4. Stays in a valley are memoryless. My run, fixed seed, 800,000 steps. Notice that the shortest stays are the most common and the counts fall steadily, the shape of an exponential law; a timetable would give a peak instead. Measured on synthetic data
Text description of figure 4

Twelve bins of 20 time units from 0 to 240. Counts: 39, 23, 13, 16, 8, 10, 2, 4, 4, 0, 1, 0. Mean stay 64 time units; the standard deviation is 1.17 times the mean, close to the value of 1 for an exponential law.

This is a small result, and that is the point. The concept of an attractor has done work here. It predicted two things in advance: a straight line on a particular plot, with a slope fixed by the hill, and a particular shape of waiting times. Both could have failed. A system that switched on a timetable would give the wrong shape. A system with no hill would give no straight line.

Explorer

Turn up the noise, watch the stays shrink

Pick a noise strength for the double well. The panel shows the mean time the state stayed in a valley in my runs, and how many crossings that average is based on.

Noise 0.45: mean stay 54 time units, from 73 crossings. At 0.35 it was 283; at 0.6 it was 19.

What this is for: feeling how steeply waiting times depend on noise when there is a real barrier. What it is not: a model of any brain; it is the simplest system with two attractors.

What the brain data would need to show#

For a brain state to count as an attractor in this exact sense, three things need to be seen. Perturbations should decay: a brief push, from a sound or a magnetic pulse, should be followed by a return to the same pattern. Time in the state should follow a law that depends on noise and on the depth of the basin, and changing either should change it in the predicted direction. And the transition should be quick compared with the time spent in each state, because the system spends little time on the hilltop.

Each of these can be measured. Something close to the first has been tried: in 2005 Marcello Massimini and colleagues gave the brain brief magnetic pulses during waking and deep sleep, and found that in sleep the response stayed local and short-lived (Massimini et al. 2005). Whether that reflects a return to an attractor is a matter of interpretation. The second and third are rarely tested at all. Most papers that call a brain state an attractor show only that the state is stable for a while. Stability for a while is what any slowly changing system shows, and it is not evidence of an attractor.

Part 2

Bifurcations, hysteresis and the instrument#

A bifurcation is a change in the number or kind of attractors as some setting of the system is turned. The setting is called a control parameter. For the brain under propofol, the obvious candidate is the drug concentration. For a pair of moving fingers, it might be the speed of a metronome.

The simplest kind is called a saddle-node or fold bifurcation. Take the double well and tilt it, by adding a constant push c to the rule. As c grows, the left valley gets shallower and the right one deeper. At a critical tilt, the left valley and the hilltop meet and vanish together. A state that was sitting in the left valley now has nowhere to rest, and it runs to the right. For the double well this happens at a tilt of about 0.385, or 2 divided by 3 times the square root of 3.

Now reverse. Tilt the other way, slowly. The state sits in the right valley, which is now the only valley at first. Even when the left valley reappears, the state stays where it is, because the right valley is still there and deep. It drops back to the left only when the right valley vanishes, at a tilt of −0.385.

So the system's state depends on the path taken. Going up, it switches at +0.385. Going down, it switches at −0.385. Between those values, two stable states exist side by side, and which one you see depends on where you came from. This is hysteresis, from the Greek for lagging behind.

A fold makes the way in differ from the way outLine chart of a hysteresis loop: the upward sweep stays low until a tilt of about 0.44, the downward sweep stays high until about minus 0.46.Tilted double well swept slowly up, then down-0.8-0.400.40.8-1.5-1011.5Tilt c (the control parameter)State xSweeping upSweeping downjumps up near +0.44drops near −0.46
Figure 5. A fold makes the way in differ from the way out. My run, small noise, 60 steps at each of 81 tilts. The theory puts the folds at ±0.385; the sweep crossed at +0.44 and -0.46, a little late because it was not infinitely slow. Notice the middle range, where both states exist and history decides which you see. Simulated (fixed seed)
Text description of figure 5

Two curves of state against tilt from −0.8 to 0.8. Sweeping up, the state stays near −1 until the tilt reaches about +0.44, then jumps to about +1. Sweeping down, it stays near +1 until about −0.46, then drops. Between the folds both states exist.

Eighty years of dynamical ideas, and when each was testedTimeline of eleven milestones from Kramers in 1940 to Toker and colleagues in 2022.From escape rates to consciousness: eleven steps1940KramersEscape over a barrier1975KuramotoCoupled oscillators1984KelsoFinger switching1985HKB modelTwo valleys for twopatterns1987Scholz et al.Slowing before theswitch2003Beggs, PlenzAvalanches in cortex2009Scheffer et al.Early warnings2010Friedman et al.Neural inertia2017Kim et al.Ring attractor in afly2018Wilting, PriesemannSubsampling fix2022Toker et al.Near-criticalconsciousness
Figure 6. Eighty years of dynamical ideas, and when each was tested. Dates are publication dates from the sources cited in this paper. Notice that the tests of a prediction often came years after the idea: the finger model in 1985, the measured slowing in 1987. Sourced (publication dates)
Text description of figure 6

Kramers 1940, escape over a barrier; Kuramoto 1975, coupled oscillators; Kelso 1984, finger switching; HKB model 1985; Scholz and colleagues 1987, slowing before the switch; Beggs and Plenz 2003, avalanches; Scheffer and colleagues 2009, early warnings; Friedman and colleagues 2010, neural inertia; Kim and colleagues 2017, ring attractor; Wilting and Priesemann 2018, subsampling; Toker and colleagues 2022, near-critical consciousness.

The confound: lag looks like hysteresis#

Back to Anjali's objection. A gap between the way in and the way out is the signature of hysteresis. But a gap also appears in any system that responds slowly to a change, if the change is made faster than the system can follow.

Picture a heavy door on a stiff spring, pushed open by a rising wind. If the wind rises quickly, the door lags behind. When the wind falls, the door lags again, on the other side. Plot door against wind and you get a loop, even though the door has only one resting position for every wind speed. There is no second valley, only slowness.

For propofol, the slowness has a name. The drug has to travel from the blood to the brain, and the pump estimates that delay from a model of an average patient. If Mr D'Souza's brain takes up the drug more slowly than the model assumes, the estimate on the screen will run ahead of his brain on the way up and behind it on the way down. A loop appears, from the pump's error alone.

The two explanations make different predictions about speed. A loop caused by lag shrinks as the change is made more slowly, because a slow change gives the system time to catch up. A loop caused by bistability does not shrink to zero, because the two valleys are there however slowly you go. Its width settles at the distance between the two fold points.

The instrument: a sweep-rate test#

That difference is the basis of the instrument for this paper. It is a test anyone with a model or an experiment can run. Sweep the control parameter up and back down at several speeds. Measure the width of the loop at each speed. Plot width against slowness.

I ran the test on two systems. The first is the tilted double well from above, which is truly bistable. The second is a monostable system with a single resting state that moves smoothly with the control parameter, but which relaxes slowly, like the heavy door. Both were swept from −0.8 to +0.8 and back, with the same small noise and fixed seeds. The speed was set by how many time steps the system was given at each of 81 settings, from 10 to 3,000.

Slowing the sweep separates a switch from a lagLine chart: the bistable loop width falls from 1.14 to about 0.77 and levels off; the monostable control falls from 1.32 to zero.Loop width against sweep slowness, for a two-state switch and a slow one-state system10301003001000300000.20.40.60.81.01.21.4Steps at each setting (log scale; more steps is a slower sweep)Loop widthTwo-state switch(bistable)Slow one-statesystem (control)Theory for theswitch: 0.77
Figure 7. Slowing the sweep separates a switch from a lag. My instrument run, fixed seeds, same noise for both systems. Notice the crossing: at the fastest sweep the control has the wider loop, 1.32 against 1.14, so one fast sweep would mislead. Slowed down, the control falls to zero and the switch settles at the theoretical width. Measured (instrument output)
Text description of figure 7

Six sweep speeds from 10 to 3,000 steps per setting. Bistable widths: 1.14, 0.96, 0.84, 0.78, 0.78, 0.76, levelling at the theoretical 0.77. Monostable control widths: 1.32, 0.52, 0.16, 0.02, 0.00, 0.02.

The results separate cleanly. At the fastest sweep, both systems show wide loops: 1.14 for the bistable system and 1.32 for the slow monostable one. The monostable loop has the larger width, so a single fast sweep would point you the wrong way. As the sweep slows, the monostable loop collapses: 0.52, then 0.16, then 0.02 and finally zero. The bistable loop narrows only a little and then stops, at 0.76 to 0.78. The theory predicts a width of 0.77, twice the fold point.

The sweep-rate test, with its actual outputsTable of loop widths at six sweep speeds for the bistable system and the monostable control.Instrument output: loop width at six sweep speedsSteps per settingTwo-state switchOne-state controlReading101.141.32Both look bistable; control looks stronger300.960.52Control loop halves1000.840.16Control loop nearly gone3000.780.02Switch at theory; control gone1,0000.780.00Unchanged3,0000.760.02Unchanged
Figure 8. The sweep-rate test, with its actual outputs. Every number in the middle columns is what the instrument returned. Notice the rule the table supports: a loop that keeps its width as the sweep slows is evidence of two states; a loop that shrinks towards zero is evidence of lag. Measured (instrument output)
Text description of figure 8

Rows for 10, 30, 100, 300, 1,000 and 3,000 steps per setting. Two-state switch: 1.14, 0.96, 0.84, 0.78, 0.78, 0.76. One-state control: 1.32, 0.52, 0.16, 0.02, 0.00, 0.02. Readings note that the control looks stronger at first and then vanishes.

The negative control. The monostable system is the control. It has no attractor structure at all, beyond its single resting state. A test that called it bistable would be broken. At fast sweeps, a test based on one loop would call it bistable, and that is exactly the mistake Kabir is at risk of making. Slowing down is what exposes it.

A loop with no second stateLine chart: two curves lag either side of a straight dashed resting line, opening a loop without any second state.The control system at a fast sweep: one resting state, slow relaxation-0.8-0.400.40.8-2-1012Control parameter cState xSweeping upSweeping downResting state
Figure 9. A loop with no second state. The negative control at 30 steps per setting. Its only resting state is the dashed line. Notice that the two sweeps still open a loop of width 0.52, because the state cannot keep up; a single sweep like this is what Kabir has. Simulated (negative control)
Text description of figure 9

The resting state is a straight dashed line from −2 at c = −0.8 to +2 at c = 0.8. The upward sweep lags below it and the downward sweep lags above it, forming a loop of width 0.52 at this speed.

Instrument

Slow the sweep yourself

Choose how many steps each system gets at each setting. More steps means a slower sweep. Compare the loop widths for the two-state switch and the slow one-state control.

10 steps: switch 1.14, control 1.32. At 300 steps: switch 0.78, control 0.02.

What this is for: seeing why one fast sweep cannot tell a switch from a lag. What it is not: a protocol for patients; safe rates in people are a separate question.

Where the instrument fails. The test has limits, and they matter for the brain. It needs sweeps slow enough for a lag to die out, and in a patient or an animal there may be no such speed that is also safe and practical. It assumes that noise is small. With large noise, the state can hop between the two valleys before either fold is reached, and the loop shrinks for that reason too, which mimics a lag. And it assumes that the system itself does not change during the sweep, which a living brain under a drug for hours may well do.

The sweep-rate test works only in one cornerSchematic heat map with sweep slowness on the x-axis and noise on the y-axis; a safe zone at bottom right, a lag confound at left and a noise confound at top.When a hysteresis loop can be read as two statesLag confound: a loop appears withor without two statesSafe: a loop that survives heremeans two statesNoise confound: hops between valleysshrink the loopSweep slowness compared with the system's relaxation time →Noise compared with barrier height →Zones set by hand from the instrument runs and the theory; boundaries are approximate.
Figure 10. The sweep-rate test works only in one corner. A schematic failure surface for the instrument. Notice that the safe zone needs both conditions at once: sweeps slow compared with the system and noise small compared with the barrier. A patient on one list sits in the coral zone. Illustrative (hand-set zones)
Text description of figure 10

A shaded grid. Left, coral: lag confound, where a loop appears with or without two states. Bottom right, teal: safe, where a loop that survives means two states. Top, indigo: noise confound, where hops between valleys shrink the loop. Zones are hand-set.

Back to theatre four#

The sweep-rate test shows what Kabir would need. A patient cannot be swept a dozen times. But the same logic was applied, with more care, in animals. In 2010 Eliot Friedman, Max Kelz and colleagues measured the drug doses at which mice and fruit flies lost and regained the righting reflex, which is a standard sign of anaesthesia in animals (Friedman et al. 2010). They found that the doses differed on the way in and the way out. In mice they measured drug concentrations in the brain itself, and the brain concentration at loss of the reflex was higher than at its return. That is the measurement that rules out the pump error.

They added a second test. Flies with a mutation in a gene called Shaker, which affects sleep, showed a collapsed loop, without obvious changes in how they took up or cleared the drug. In mice unable to make the signalling chemical noradrenaline, the loop was altered, and restoring that signalling in the nervous system brought it back to normal. The authors proposed a barrier between the waking and anaesthetised states, a tendency of the nervous system to resist a change of state, and called it neural inertia.

This is a case where the hysteresis concept carries weight. It predicted a loop that survives the removal of drug lag, and the loop survived. It predicted that changing the barrier, here by a gene, would change the loop, and it did.

Kabir's two numbers are not evidence of this. They are what both stories predict. Anjali's suggestion, when he brings it to her, is to look at the patient's brain activity on the monitor at both crossings, and to see whether the pattern in each state is the same whichever way it was entered. That is a second signature of two attractors, and it costs nothing to check.

Part 3

Early warnings and the false alarm#

The finger experiment found that recovery slows before a switch. That observation has a general form. As a valley grows shallow, the force pulling the state back to the bottom grows weak. A small push takes longer to die away. The system becomes sluggish, and it also wanders further, because the weak pull lets noise carry it further from the bottom.

This is called critical slowing down. It can be measured without pushing the system at all, from its natural fluctuations. Take a record of the state over time. Measure how similar each moment is to the one before, a quantity called the lag-one autocorrelation. Near a stable bottom, the state forgets quickly and the autocorrelation is low. As the valley flattens, the state forgets slowly and the autocorrelation rises towards one. The variance, the spread of the fluctuations, rises too.

In 2009 Marten Scheffer and colleagues gathered evidence that this rising autocorrelation and variance appear before sudden shifts in many systems: lakes turning murky, climate changes in the past, fish populations collapsing and epileptic seizures (Scheffer et al. 2009). The idea spread quickly to neuroscience and psychiatry, as a possible warning before a seizure, a relapse of depression or a change in conscious state.

A rising warning with two possible causesFour-panel storyboard of the lake thought experiment: a clear lake, a rising warning, a change in the weather as a second cause and a storm-driven flip with no warning.Thought experiment 3: the lake that slowed1clear lakeA clear lakeFertiliser run-off riseseach year; a flip to murkywater is feared.2murky dayscome in runsA warning appearsDay-to-day clarity becomesmore strongly correlatedover five years.3weather spellsgot longerA second causeThe weather that stirs thelake now comes in longerspells.4one storm:flip, no warningA flip with no warningOne great storm kicks adeep-valley lake over thehill.
Figure 11. A rising warning with two possible causes. Four panels of the invented thought experiment, read left to right. Notice that panels 2 and 3 show the same measurement with two different causes, and panel 4 shows a transition the measurement would miss. Illustrative (thought experiment)
Text description of figure 11

Panel 1: a clear lake with rising run-off. Panel 2: clarity becomes more autocorrelated. Panel 3: the weather spells have lengthened, which would produce the same rise. Panel 4: a single storm flips the lake with no prior warning.

A worked case: three ways to see a warning#

I ran three versions of the double well, each for 6,000 steps with fixed seeds. In the first, the tilt rises slowly towards the fold point, so that the valley really does flatten. In the second, the tilt stays fixed, far from the fold, but the noise becomes more correlated over time, like the changing weather. In the third, the tilt stays fixed and the noise is large, so that the state can be kicked over the hill without warning.

For each run I measured the lag-one autocorrelation and the variance in a moving window, up to the moment the state crossed the hill or the run ended. I summarised the trend with a standard rank measure called Kendall's tau, which is +1 if the measure only ever rises and 0 if it has no trend. Ecologists usually treat a tau above about 0.5 as a warning.

Two warnings, one false, and a transition with noneLine chart of three autocorrelation traces: the fold case and the changing-noise case both rise; the large-noise case stays flat and still tips.Lag-one autocorrelation in a moving window, three versions of the double well00.250.50.7510.50.60.70.80.91.0Fraction of the run elapsedLag-one autocorrelationApproaching a fold(tipped)Changing noise, nofold (did not tip)Large noise, farfrom fold (tipped)
Figure 12. Two warnings, one false, and a transition with none. My runs, fixed seeds, window 400 steps, each curve ending when the state crossed the hill or the run ended. Trend scores (Kendall's tau): approach 0.69, changing noise 0.89, large noise 0.11. Notice that the false alarm rises highest. Measured on synthetic data
Text description of figure 12

Three traces of lag-one autocorrelation over the run. Approaching a fold: rises, trend 0.69, and the state tips. Changing noise with no fold: rises most, trend 0.89, and the state never tips. Large noise far from the fold: flat, trend 0.11, and the state tips anyway.

The approaching fold gave a clear warning. In the run shown, the autocorrelation trend was 0.69 and the variance trend 0.68, and the state tipped. That is the textbook case.

The changing noise gave a stronger warning and no transition. In the run shown, the autocorrelation trend was 0.89, higher than for the real fold, and the state never tipped. This is the false alarm. The measure cannot tell a flattening valley from a change in the noise that drives the system.

The large-noise run tipped with no warning. Its autocorrelation trend was 0.11. This is the counterexample: a real transition between two attractors with no slowing down beforehand. Peter Ditlevsen and Sigfus Johnsen made the same point for climate records in 2010 (Ditlevsen and Johnsen 2010). Slowing down is expected only when a transition is caused by a valley vanishing. Noise-driven jumps between stable valleys give no warning, and many brain transitions may be of that kind.

To check that these runs were typical, I repeated each case 40 times with different seeds. The approaching fold gave a trend above 0.5 in 95 per cent of runs, and tipped in all of them. The changing noise gave a trend above 0.5 in every run, and tipped in only 12 per cent. The large noise gave a trend above 0.5 in 5 per cent of runs, and tipped in 53 per cent.

Over 40 runs, the warning tracks the noise as much as the foldBar chart in three pairs: fold warned 95 per cent and tipped 100; changing noise warned 100 and tipped 12; large noise warned 5 and tipped 53.Share of 40 runs with a warning (trend above 0.5) and share that tippedPer cent of runs95%Fold: warned100%Fold: tipped100%Changing noise:warned12%Changing noise:tipped5%Large noise: warned53%Large noise: tippedWarning: Kendall's tau of the autocorrelation trend above 0.5. Tipped: the state crossed the hill within 6,000 steps.
Figure 13. Over 40 runs, the warning tracks the noise as much as the fold. My replicate runs, 40 seeds per case. Notice the middle pair: changing noise warned in every run and tipped in 12 per cent, and the last pair: large noise tipped in 53 per cent with a warning in only 5 per cent. Measured on synthetic data
Text description of figure 13

Six bars. Approaching a fold: warning in 95 per cent of runs, tipped in 100 per cent. Changing noise: warning in 100 per cent, tipped in 12 per cent. Large noise: warning in 5 per cent, tipped in 53 per cent.

In 2012 Carl Boettiger and Alan Hastings pointed to a further trap (Boettiger and Hastings 2012). If you pick out, from many records, the ones that ended in a transition, and then look for warnings before them, you will find warnings more often than chance, even in systems that give none. The selection itself manufactures the signal. Their fix is to state in advance what you will measure and on which records.

Predict first

Which run gave the strongest warning?

Three versions of the double well were run 40 times each. In one the valley flattened towards a fold; in one the noise became more correlated with no fold; in one the noise was large. Which gave a rising-autocorrelation warning most often?

The changing noise: it warned in every one of 40 runs and tipped in only 12 per cent. The real fold warned in 95 per cent. The large-noise runs tipped in 53 per cent with a warning in only 5 per cent.

What this is for: testing your intuition about early warnings before reading the numbers. What it is not: a claim about any particular brain or lake.

What this means for anaesthesia and sleep#

The question for consciousness research is which kind of transition the brain makes. If the loss of consciousness under a drug is a fold, warnings should appear as the drug rises. It could instead be a noise-driven hop between valleys that both exist across a wide range of doses. Then warnings should be absent, and the moment of the transition should vary from one induction to the next in a way that looks random.

Neither of these is known in general. Studies report changes in the brain's rhythms before loss of consciousness, but few are designed to separate a flattening valley from a change in the drive. For the drug case the drive is the drug itself, so a change in the noise is hard to rule out. The early-warning concept carries weight only when the noise is known to be steady, or when the warning is checked against a record where no transition followed.

The active test: push and time the recovery#

There is a better version of the warning, and the finger experiments used it. Instead of reading slowness from the system's own fluctuations, push it and time how long it takes to recover. Scholz, Kelso and Schöner did this by briefly disturbing one finger and measuring how long the pattern took to return (Scholz, Kelso and Schöner 1987). As the metronome sped up towards the switch, recovery took longer.

The active test has an advantage over the passive one. Its result does not depend on the colour of the background noise, because the push is a known, separate event. A lake whose weather has changed will show rising autocorrelation, but a known disturbance to it will still be absorbed at the old speed if its valley has not flattened. Recovery time after a deliberate push measures the depth of the valley directly, where autocorrelation measures it only through the noise.

In the brain, the push can be a sound, a flash or a magnetic pulse, and the recovery is the brain's response returning to its baseline. Massimini's pulses were a step in that direction. A study that delivered the same push at several drug doses on the way into anaesthesia, and timed the recovery, would test for a fold far more cleanly than any passive measure.

Part 4

Metastability: states that never quite settle#

An attractor is a state the system goes to and stays in. Much of what a brain does seems different. Thoughts do not settle. Patterns of activity form, hold for a fraction of a second, and give way to others. Kelso, who began with the fingers, later gave this a name in the brain: metastability, a regime in which the parts of a system tend towards coordination without becoming locked in it (Tognoli and Kelso 2014).

The word has two meanings, and this causes trouble. In physics, a metastable state is a shallow valley, like the supercooled water: stable for a while, but not the deepest state. In brain dynamics, it often means something else: a system with no stable states at all, which passes near a series of almost-stable patterns, lingering at each one before moving on. Mikhail Rabinovich and colleagues described one version of this, a chain of saddle points, each of which the state approaches and then leaves along a particular route (Rabinovich, Huerta and Laurent 2008).

A saddle point is a balance point that attracts from some directions and repels in others, like a mountain pass. A state can travel towards a pass along the valley floor, slow down near it, and then leave along a different valley. A chain of such passes, each leading to the next, gives a sequence of patterns that repeats in a fixed order with dwell times at each one. Rabinovich's group proposed this as a model for sequences in the insect sense of smell and in other neural processes that unfold in a fixed order.

The two meanings make different predictions. Shallow valleys predict memoryless, exponential dwell times, as in Part 1, and a random order of visits. A chain of passes predicts a fixed order and dwell times that grow longer, run by run, as noise falls. A study that calls the brain metastable should say which kind it means, because the data can tell them apart.

A worked case: oscillators that half agree#

The standard model is a set of oscillators, each with its own natural rhythm, pulling on one another. It was introduced by the Japanese physicist Yoshiki Kuramoto in 1975, and Steven Strogatz later wrote a clear history of it (Strogatz 2000). Each oscillator is a point going round a circle at its own speed. The coupling strength, K, sets how strongly each one is pulled towards the average of the others.

With weak coupling, the oscillators run independently and their average is close to zero. With strong coupling, they lock together and move as one. Between these, they half agree. Clusters form and dissolve, and the degree of agreement rises and falls over time.

The degree of agreement is measured by a number called the order parameter, R. It is 0 when the oscillators are scattered round the circle and 1 when they move as one. A common measure of metastability is how much R varies over time. A locked system has a steady R near 1. A scattered one has a steady R near 0. A metastable one has an R that swings.

I ran 64 oscillators with natural speeds drawn from a bell curve, at eleven coupling strengths. Metastability, the standard deviation of R over time, peaked at 0.11 near K = 1.4, and fell to zero once the oscillators locked at K = 3 and above. The theory for an infinite population puts the onset of synchrony at K = 1.6 for this spread of speeds. The peak sits just below it, where clusters keep forming and breaking.

Metastability peaks just below the onset of synchronyLine chart: the metastability index rises to a peak of 0.11 at coupling 1.4 and falls to zero above 3; an uncoupled control sits at 0.05.Kuramoto model, 64 oscillators: swing in agreement against coupling strength0.511.6234500.020.040.060.080.10.12Coupling strength KStandard deviation of RMetastability:spread of R overtimeUncoupled control:0.051theory: synchrony begins at 1.6locked: R above 0.93
Figure 14. Metastability peaks just below the onset of synchrony. My run, fixed seed, natural speeds from a bell curve of width 1. The peak, 0.11 at K = 1.4, sits just below the onset of synchrony. Notice the dashed control: with no coupling at all, finite size alone gives 0.05, almost half the peak. Measured on synthetic data
Text description of figure 14

Metastability index at eleven coupling strengths from 0.5 to 5: 0.071, 0.089, 0.112, 0.106, 0.042, 0.029, 0.011, 0.010, then 0 from K = 3. The mean order parameter R rises from 0.14 to 0.98 over the same range. A dashed line marks the uncoupled control at 0.051.

The control. With no coupling at all, R still varied, with a standard deviation of 0.05. That variation comes only from having a finite number of oscillators, which drift into and out of chance alignment. So half of the peak value is what you would get from a system that does nothing. A metastability index means little until it is compared with the value for the same number of uncoupled parts.

A true story of an attractor seen directly#

Most claims of attractors in the brain are inferred from behaviour or from large-scale recordings. One is close to direct. In 2015 Johannes Seelig and Vivek Jayaraman recorded from a ring of neurons in the brain of the fruit fly, in a structure called the ellipsoid body (Seelig and Jayaraman 2015). They found a single bump of activity that moved round the ring as the fly turned, tracking its heading like the needle of a compass. The bump persisted in darkness, when the fly had no view to follow.

Theorists had long proposed a network for exactly this, called a ring attractor. Neurons arranged in a ring excite their neighbours and inhibit everything else, so that one bump of activity can sit anywhere on the ring and stay there. In 2017 Sung Soo Kim, Hervé Rouault, Shaul Druckmann and Jayaraman tested the prediction directly (Kim et al. 2017). They used light to switch on neurons at a different place on the ring, writing in a false heading. The circuit took the new bump and held it, with the old one gone. The network enforced one bump, as a ring attractor must.

A perturbation moved the bump, and the circuit held the new oneTwo rings of sixteen neurons: before, a bump of activity near the top right; after optogenetic stimulation, a bump at the bottom left and none at the old place.A ring attractor for heading: one bump, held in place, moved by a pushBeforebump marks the headingAfternew bump held; old onegonelight writes a new bumpSchematic of the design in Kim et al. 2017: local excitation and global inhibition allow one bump only.
Figure 15. A perturbation moved the bump, and the circuit held the new one. A schematic of the fly experiment in Part 4, drawn from the published description, not from its data. Notice what makes it strong evidence: the state was pushed to a new place, and the network kept it there with one bump, as a ring attractor must. Schematic (after Kim et al. 2017)
Text description of figure 15

Left ring: sixteen neurons with a bump of activity at one position, marking the fly's heading. Arrow: light switches on neurons elsewhere. Right ring: the bump now sits at the new position and the old bump has gone. A note says the design relies on local excitation and global inhibition.

This is the strongest kind of evidence for an attractor: a perturbation that moves the state, followed by the predicted persistence. It concerns heading, not consciousness. But it shows that the concept can be tested in neural tissue and pass.

What metastability can and cannot explain#

Gustavo Deco, Viktor Jirsa and Anthony McIntosh argued in 2010 that the brain's resting activity sits near the edge of instability, where it can explore many patterns (Deco, Jirsa and McIntosh 2010). Michael Breakspear's 2017 review set out how models of large-scale brain dynamics produce such behaviour, and how they can be fitted to recordings (Breakspear 2017). Francisco Varela and colleagues had earlier proposed that brief synchrony between distant regions binds them into a single moment of experience, with the synchrony constantly forming and dissolving (Varela et al. 2001).

These ideas share a claim: that a flexible, conscious brain lives between too little coordination and too much. That claim is testable, and some of it has support. Measures of how varied the brain's patterns are tend to drop under deep anaesthesia and in dreamless sleep. But a drop in variety is predicted by many accounts, including ones with no special role for metastability. The concept carries weight only if it predicts something more specific, such as the dwell times of each pattern, or how they change as coupling is altered by a drug, and if those predictions beat the deaf choir.

Part 5

Criticality and the subsampling trap#

The last concept is the most ambitious. A system is critical when it sits exactly at the boundary between two kinds of behaviour, such as order and disorder. At that boundary, events of every size occur. A small disturbance can die out at once or spread through the whole system, and no size is typical. In 1987 Per Bak, Chao Tang and Kurt Wiesenfeld proposed that some systems tune themselves to this boundary without any outside help, which they called self-organised criticality, using a pile of sand as the example (Bak, Tang and Wiesenfeld 1987).

The simplest model for the brain is a branching process. Each active neuron triggers, on average, m others in the next moment. If m is below 1, activity dies out. If it is above 1, activity grows until something stops it. At m = 1 the system is critical, and bursts of activity, called avalanches, come in all sizes.

In 2003 John Beggs and Dietmar Plenz recorded from slices of rat cortex on a grid of electrodes and found such avalanches (Beggs and Plenz 2003). Their sizes followed a power law, a straight line on a log-log plot, with the slope expected for a critical branching process. The idea that the brain is critical took off, and so did the claim that criticality is good for computation, because a critical system responds to the widest range of inputs (Chialvo 2010).

The trap: what you do not record#

A recording sees only a tiny fraction of the neurons in a brain. That matters more than it seems. In 2018 Johannes Wilting and Viola Priesemann showed that estimating m from a subsample gives a badly biased answer, and that the bias makes systems look less critical than they are (Wilting and Priesemann 2018). They proposed a fix, which they called a multistep regression estimator.

I ran the test myself. A branching process ran for 100,000 steps with a true m of 0.98, close to critical. I then recorded only a random fraction of its events: all of them, a tenth, a hundredth and a thousandth. For each, I estimated m in two ways. The naive way regresses the activity at each step on the activity at the step before. The multistep way regresses across many lags and reads m from how the relationship decays with the lag.

Recording a hundredth of the events hides a near-critical processBar chart: naive estimates fall from 0.98 to 0.72, 0.20 and 0.02 as recording thins; multistep estimates stay at 0.98.Estimated branching ratio of a process with true value 0.98, by share of events recordedEstimated branching ratio0.98Naive, all0.98Multistep, all0.72Naive, 1 in 100.98Multistep, 1in 100.20Naive, 1 in1000.98Multistep, 1in 1000.02Naive, 1 in1,0000.98Multistep, 1in 1,000Naive: regression of each step on the one before. Multistep: Wilting and Priesemann's regression across many lags.
Figure 16. Recording a hundredth of the events hides a near-critical process. My run, 100,000 steps, fixed seeds. Notice the coral bars falling to 0.20 at one in a hundred while the indigo bars stay at 0.98. As a negative control, a process with no branching gave a naive estimate of −0.004. Measured on synthetic data
Text description of figure 16

Eight bars in four pairs. All events recorded: naive 0.98, multistep 0.98. One in ten: naive 0.72, multistep 0.98. One in a hundred: naive 0.20, multistep 0.98. One in a thousand: naive 0.02, multistep 0.98.

With every event recorded, both methods gave 0.98. With a tenth recorded, the naive estimate fell to 0.72. With a hundredth, it fell to 0.20, and with a thousandth to 0.02. The multistep estimate stayed at 0.98 throughout, and 0.977 at a thousandth. The same process, recorded at a hundredth, looks to the naive method like one far from criticality. As a negative control, I ran the estimator on activity with no branching at all, subsampled to a hundredth. It gave −0.004, correctly close to zero.

The other trap: power laws without criticality#

The second trap runs the other way. A power law in avalanche sizes is not proof of criticality. In 2017 Jonathan Touboul and Alain Destexhe showed that power laws can arise in models that are not critical at all, including some driven only by noise (Touboul and Destexhe 2017). So can the scaling relationships that were thought to be a stronger sign. A straight line on a log-log plot is what many processes produce over a limited range of sizes.

Criticality and consciousness#

The bold version of the claim links criticality to consciousness directly. In 2014 Robin Carhart-Harris and colleagues proposed that the brain under psychedelic drugs moves towards criticality and greater disorder, and that normal waking consciousness sits just below it (Carhart-Harris et al. 2014). In 2022 Daniel Toker and colleagues analysed cortical recordings from a range of conscious and unconscious states. They argued that during conscious states the brain's slow electrical rhythms sit near the boundary between stability and chaos, and that they move away from it during unconsciousness (Toker et al. 2022). In 2019 Fontenele and colleagues reported that in anaesthetised rats, signs of criticality appeared only in a narrow range of cortical states, between more synchronised and less synchronised regimes (Fontenele et al. 2019).

These are serious studies. A careful review in 2022 by Jordan O'Byrne and Karim Jerbi set out how much remains unsettled (O'Byrne and Jerbi 2022). Different measures of criticality disagree. The subsampling and power-law traps are common. And "near critical" often comes without a stated distance. Criticality carries weight in consciousness research only when three things are in place. The estimator must survive subsampling. The claim must say how far from critical each state is, in numbers. And moving the brain along that axis must change consciousness in the direction predicted, not merely go along with it.

What a good criticality study would report#

It is easy to say what such a study would contain. It would estimate the branching ratio with a method that survives subsampling, and report the naive estimate beside it so that readers can see the size of the correction. It would report the distance from critical as a number with an uncertainty, for each state. It would include a control recording in which the same analysis should find no criticality, such as activity shuffled in time. And it would test a prediction that only a critical account makes, such as a sharp change in the response to small pushes as the drug crosses a particular dose.

Some recent work meets several of these conditions. Few studies meet all of them. Until more do, the claim that consciousness requires criticality is a promising hypothesis, and it should be graded as one.

Part 6

Spanda: a pulse at the seam between states#

The Kashmiri Śaiva tradition has a text devoted to change itself. The Spandakārikā, the "Stanzas on Vibration", is a short work of the ninth century, attributed to Vasugupta or to his pupil Bhaṭṭa Kallaṭa (Spandakārikā, GRETIL; Spandakārikā with Kṣemarāja, GRETIL; Dyczkowski 1992). Its subject is spanda, a word that means a slight movement, a throb or a pulse. In the text, spanda is the activity of consciousness itself, the movement by which every state of mind arises and passes.

This paper's working theory, the Recognition Theory of Mind (RTM), borrows the word for its first layer. RTM is a hypothesis under test, with five layers: L1 spanda, the pulse of activity; L2 prakāśa, the shining or showing of something; L3 vimarśa, reflective grasp; L4 pratyabhijñā, recognition; and L5 cit, awareness as such. If any part of RTM should meet dynamical systems theory, it is the first layer.

What the stanzas say#

Three verses bear directly on this paper. The translations are my own close renderings of the Sanskrit on GRETIL, whose numbering I follow, checked against Mark Dyczkowski's translation.

The opening verse praises the god "by whose opening and closing of the eyes the world dissolves and arises" (yasyonmeṣanimeṣābhyāṃ jagataḥ pralayodayau, 1.1). The pair unmeṣa and nimeṣa, opening and closing, is the text's basic image. States of the world, and of mind, come and go as a rhythm of opening and closing.

The third verse turns to states of consciousness. "Even in the differences of waking and the other states, it spreads through them without being divided from them. It never departs from its own nature as the perceiver" (jāgradādivibhede 'pi tadabhinne prasarpati, nivartate nijān naiva svabhāvād upalabdhṛtaḥ, 1.3). The states change. Something does not.

The third section gives a practice: "When one is absorbed in a single thought and another arises, that is unmeṣa; one should observe it for oneself" (ekacintāprasaktasya yataḥ syād aparodayaḥ, unmeṣaḥ sa tu vijñeyaḥ svayaṃ tam upalakṣayet, 3.9). An earlier verse makes the same move with stronger states. In great anger, in delight, in the moment of asking "what shall I do?", or while running, "the state one reaches, there spanda is established" (1.22).

Where the stanzas and the mathematics agree#

There are three points of contact, and I state them as correspondences, not identities.

First, both treat change as primary. Dynamics describes a system by its rule of motion, not by a list of its states. The stanzas describe consciousness by its pulse, and its states as what the pulse produces. Second, both direct attention to the seam between states. The informative moments in a dynamical system are its transitions: the dwell times, the slowing before a fold, the path of escape. Verses 1.22 and 3.9 place spanda at exactly such moments. Third, both distinguish what changes from what persists. A dynamical system has states that change and a rule that does not. Verse 1.3 has states that change and a perceiver that does not.

The stanzas and the mathematics agree on where to look, and part on who looksTable comparing six features of the Spandakārikā with dynamical-systems ideas, with a judgement of how closely each corresponds.Where the Spandakārikā and dynamical systems correspondIn the SpandakārikāIn dynamical systemsCorrespondenceUnmeṣa and nimeṣa: states arise and dissolve asan opening and closing (1.1)A system is described by its rule ofmotion, not by a list of statesHolds as an orientationSpanda is found at peaks: anger, delight, "whatshall I do?" (1.22)Informative events are transitions: dwelltimes, slowing, escapeHolds closelyWatch the arising of a new thought (3.9)Measure transitions between metastablepatternsHolds; a first-personmethodThe perceiver persists through waking and theother states (1.3)The rule persists while states changeBreaks: a rule has nopoint of viewAgitation (kṣobha) dissolving reveals thehighest state (1.9)Lower noise deepens and steadies anattractorHolds looselyThe aim is liberationThe aim is predictionBreaks
Figure 17. The stanzas and the mathematics agree on where to look, and part on who looks. A break-point map for Part 6. Verse numbers follow the GRETIL text; the left column is my paraphrase and the right column my judgement. Notice the fourth row: the persisting perceiver is where the bridge breaks. Textual reading and judgement
Text description of figure 17

Six rows. Opening and closing (1.1) matches description by a rule of motion, holds as an orientation. Spanda at peaks (1.22) matches transitions as informative events, holds closely. Watching a new thought arise (3.9) matches measuring transitions, holds as a first-person method. The persisting perceiver (1.3) against a persisting rule breaks. Agitation dissolving (1.9) matches lower noise, holds loosely. Liberation against prediction breaks.

Where they part#

The third correspondence is also where the bridge breaks. In dynamics, what persists is the rule, a mathematical object with no point of view. In the stanzas, what persists is the perceiver, upalabdhṛ, the one who knows each state. These are different kinds of thing. A rule of motion does not experience the states it generates. The stanzas say the perceiver does.

What this means for RTM#

The tradition earns a place in this paper on one point: it tells you where to look. If RTM's first layer is the brain's activity and its rhythms, the tradition's advice is to study the transitions, not the states, and to study them from the inside as well as from the outside. That advice matches the paper's thesis. The signatures that make dynamical concepts carry weight are signatures of transitions.

It also sets a test for RTM itself. If L1 spanda is the layer at which dynamics applies, then the measured statistics of transitions should constrain the higher layers. A change in dwell times, under a drug or with practice, should show up in what is shown (L2) and in what is grasped (L3). If transition statistics could change freely without any effect on the higher layers, the first layer would be doing no work in the theory.

[AUTHOR NOTE: if your own practice includes the kind of attention to the arising of thoughts described in Spandakārikā 3.9, a short account of what you noticed, and how reliably, would help readers judge the thought experiment.]

Switch the lens

One induction, four readers

The same moment, when Mr D'Souza stops squeezing the hand, read by four people with different questions.

I need to know the dose at which he goes under and the dose at which he comes back, for this patient, today. Whether the gap is inertia or lag changes how I time the end of the case.

I want the loop at three speeds, and the waiting time to response at a fixed dose. A loop that keeps its width and memoryless waits would tell me there are two attractors.

Before any warning signal is believed, I want a record with no transition analysed the same way, and the analysis fixed before anyone looks.

The stanzas would point to the moment itself, the seam between one state and the next, and ask what persists across it. That is a question about the perceiver, which no curve answers.

What this is for: seeing that one transition raises four different questions. What it is not: a claim that the four readings are equally testable.

The strongest case against

Everything is a dynamical system, so the words explain nothing#

Here is the strongest case against my thesis, put as its best defenders would put it. It comes from two directions at once.

The first says I am too generous. Every physical system can be written as a dynamical system. Every system that stays in a state for a while has, in some description, an attractor. Every system that changes suddenly has, in some description, a bifurcation. If the vocabulary fits everything, it explains nothing about consciousness in particular. The signatures I demand, such as exponential dwell times and loops that survive slow sweeps, are generic properties of noisy nonlinear systems. Finding them in a brain tells you that the brain is a noisy nonlinear system, which no one doubted.

The second says I am too strict. Tim van Gelder argued in 1998 that cognition is best understood as the behaviour of a dynamical system, and that asking dynamics to prove itself against rival explanations misreads the situation (van Gelder 1998). On his view, the rivals, such as computational models with stored symbols, are themselves dynamical systems described at a coarser level. There is no neutral test between them. The question is which description is most useful, and dynamics wins on usefulness wherever timing matters.

I take the first objection seriously, and it sharpens my claim. The signatures are generic, and that is why they are useful. A generic signature tells you which class of mechanism you are in, and the class is what the rival explanations deny. The drug-lag story says there is one valley. The two-attractor story says there are two. A loop that survives slow sweeps decides between them, whatever else is true. The vocabulary says nothing about consciousness in particular, I agree. It says something about the shape of the transitions between conscious and unconscious states, and that is a question about consciousness.

The second objection I accept in part. Van Gelder is right that dynamics and computation are not opposites, and Part 4's ring attractor is both: a network that computes heading by being an attractor. He is wrong, I think, that there is no neutral test. In each worked case of this paper there was a rival with different predictions. The rival was not always computational. It was a lag, a change in the noise, a deaf choir or a subsampled recording. The test is not dynamics against computation. It is one dynamical story against another.

Finally, there is a version of the case against that I cannot answer. It says that the brain is so high in dimension, and so far from the simple models, that the clean signatures may never be measurable in it. A double well has one variable. A brain has billions. If that is right, the concepts may be true of the brain and still never carry weight in explaining it, because their predictions cannot be checked. I hold the opposite view on the evidence of the ring attractor and of Friedman's loop, which were checked. But two cases do not settle a field.

I should also say what I concede. The vocabulary of dynamics has shaped good experiments even where the concepts were used loosely. The idea that the brain might sit near a critical point sent people to record avalanches, and those recordings are useful whatever the final verdict on criticality. A loose concept can be a good guide to where to look. My claim is narrower: it cannot be an explanation of what is found until it has passed a test that a rival could fail.

Rival views of dynamics in mind science, on two axesScatter plot of seven views of dynamics in mind science, placed by explanatory weight given to dynamics and by strictness about tests.Comparison atlas: how much each view credits dynamics, and how much it demands of testsStrict, scepticalStrict, dynamicalNeitherDynamical, permissiveHow much explanatory weight the view gives dynamical concepts →How strict the view is about tests →lowhighlowhighEach concept must earn itsplace by a signature (thispaper)Cognition is dynamics; noneutral test (van Gelder)Coordination dynamics withtested signatures (Kelso)Consciousness nearcriticality (Toker;Carhart-Harris)Criticality claims needstricter tests (O'Byrne,Jerbi)Early warnings are oftenwishful (Ditlevsen, Johnsen)Brain states as computation,dynamics as detail
Figure 18. Rival views of dynamics in mind science, on two axes. Positions are my judgement of each view's central commitments, not measurements. Notice that my position sits with Kelso's, in the top-right: dynamics carries weight, but only through tested signatures. Judgement (positions are interpretive)
Text description of figure 18

This paper, centre right and high on strictness. Van Gelder, far right and low. Kelso, right and high. Toker and Carhart-Harris, right and middle. O'Byrne and Jerbi, left and high. Ditlevsen and Johnsen, left and fairly high. Computation-first views, bottom left.

Refutation conditions

What would change my mind#

Concrete observations would weaken the thesis, and I list them so they can be checked.

  • If a dynamical concept, used without any distinguishing signature, were shown to predict new results about conscious state transitions that rival explanations missed, my demand for signatures would be too strict.
  • If the sweep-rate test were run on an anaesthetic in animals with brain concentrations measured, and the loop collapsed to zero at the slowest safe rate, neural inertia would lose its main support and the hysteresis concept would carry no weight there.
  • If early-warning signals before loss of consciousness were found to be reliable, with a matched false-alarm rate below one in five in records with no transition, critical slowing down would carry more weight than I give it.
  • If subsampling-corrected estimates of the branching ratio were the same in waking, sleep and anaesthesia, the link between criticality and consciousness would fail at the level of that measure.
  • For RTM, if the statistics of transitions in brain activity could be changed widely, by drugs or by practice, with no change in reported experience, the first layer would be doing no work in the theory.

What I still don't know#

I don't know whether the loss of consciousness under a drug is a fold or a noise-driven hop, and I don't know of a human study designed to tell them apart. I don't know whether any brain-wide metastability index beats its own uncoupled control by a wide margin in real data, because few studies report the control. I don't know how far the subsampling correction carries over from a simple branching process to brain recordings, where the activity is not a single branching chain. And I don't know whether first-person reports of the arising of thoughts can be made reliable enough to count as data about dwell times.

There is one more gap I should own. My simulations use one-variable and simple many-variable models. They show what each signature looks like and how it can be faked. They do not show that any of these signatures is present in a human brain.

The open questions, placed by how answerable they areScatter plot of five open questions placed by importance and by how answerable they are.What I don't know, by importance and tractabilityAnswerable, minorDo these firstPark for nowImportant but hardImportance to the thesisHow answerable with current toolslowhighlowhighIs loss of consciousness afold or a noise-driven hop?Does a brain metastabilityindex beat its uncoupledcontrol?Does the subsampling fixcarry over to brainrecordings?Can first-person reports ofthought onsets be reliable?Can a sweep-rate test be runsafely in people?
Figure 19. The open questions, placed by how answerable they are. My placement of the open questions from the refutation section. Notice the top-right pair: both can be tested by reanalysing existing recordings with the controls described in Parts 4 and 5. Judgement
Text description of figure 19

Fold or hop: very important, middling tractability. Metastability against control: important, answerable. Subsampling fix in brain data: important, answerable. First-person reports: middling importance, hard. Sweep-rate test in people: important, hard.

Claim Ledger

Eleven claims, graded

Each main claim, how sure I am, what kind of evidence carries it and what observation would refute it. Confidence uses three grades: established, likely and speculative.

IDClaimConfidenceEvidence typeWhat would refute it
C13.1A dynamical concept carries explanatory weight in a theory of state transitions only when it predicts a signature that a rival account does not.LikelyArgument; worked cases in this paperA dynamical concept, used without a distinguishing signature, that predicts new results its rivals missed.
C13.2For a noisy two-valley system, mean time in a valley follows Kramers' law and stays are close to exponential.EstablishedPublished theory (Kramers 1940); measured here (slope 0.49 against 0.5)None plausible for this model class; the open question is whether brain states belong to it.
C13.3A hysteresis loop from a single sweep cannot distinguish two stable states from a lag.EstablishedMeasured in this paper (control loop 1.32 against 1.14 at fast sweep)None plausible; this follows from the dynamics of lag.
C13.4A loop that keeps its width as the sweep slows is evidence of bistability; one that shrinks to zero is evidence of lag.EstablishedMeasured in this paper; standard bifurcation theorySystems where noise-driven hopping shrinks a true bistable loop at slow sweeps, which the instrument would misread.
C13.5Anaesthetic hysteresis in mice and flies is not fully explained by drug lag, and supports a barrier between states.LikelyPublished evidence (Friedman et al. 2010)Replications with brain concentrations finding the loop collapses at slow enough rates.
C13.6Rising autocorrelation can be produced by changing noise with no approaching transition.EstablishedMeasured here (warning in 40 of 40 runs, 12 per cent tipped); published (Ditlevsen and Johnsen)None plausible; this is a property of the measure.
C13.7Noise-driven transitions between stable states give no early warning.EstablishedMeasured here (5 per cent warned, 53 per cent tipped); published theoryNone plausible for noise-driven hops.
C13.8A metastability index is uninterpretable without the value for an uncoupled system of the same size.LikelyMeasured here (control 0.05 against peak 0.11)Brain data where the uncoupled value is negligible for every practical size.
C13.9Naive branching-ratio estimates from subsampled recordings are biased towards zero; multistep regression corrects this.EstablishedPublished (Wilting and Priesemann 2018); measured here (0.20 against 0.98 at 1 in 100)None plausible for branching processes; transfer to brain data is open.
C13.10Conscious states sit nearer a critical point than unconscious ones.SpeculativePublished evidence (Toker et al. 2022; Fontenele et al. 2019) with open methodological issuesSubsampling-corrected estimates showing no difference between conscious and unconscious states.
C13.11RTM's first layer, read as spanda, predicts that changes in transition statistics constrain what is shown and grasped at higher layers.SpeculativeArgument; textual reading of the SpandakārikāWide changes in transition statistics with no change in reported experience.

Dated predictions

IDDated predictionJudged byHow it will be judged
P13.1A published study that changes an anaesthetic at three or more rates, with brain or effect-site concentration measured, will find an induction–emergence gap that shrinks with slower change but stays above zero at the slowest rate.31 December 2028Pass if one such study reports a gap that levels off above zero; fail if two such studies find it falls to within measurement error of zero.
P13.2When naive and subsampling-corrected branching-ratio estimates are compared on the same human intracranial recordings in waking, the corrected estimate will be closer to 1 by at least 0.05.30 June 2029Pass if a published comparison shows the difference; fail if a comparison on at least ten patients finds differences under 0.02.
The mechanisms are settled; the claims about consciousness are notDot chart of eleven claims by confidence grade.Claim Ledger: confidence by claimRefutedSpeculativeLikelyEstablishedC13.1 A dynamical concept carries explanatoryweight in a theory…C13.2 For a noisy two-valley system, mean time ina valley follo…C13.3 A hysteresis loop from a single sweepcannot distinguish t…C13.4 A loop that keeps its width as the sweepslows is evidence…C13.5 Anaesthetic hysteresis in mice and flies isnot fully expl…C13.6 Rising autocorrelation can be produced bychanging noise w…C13.7 Noise-driven transitions between stablestates give no ear…C13.8 A metastability index is uninterpretablewithout the value…C13.9 Naive branching-ratio estimates fromsubsampled recordings…C13.10 Conscious states sit nearer a criticalpoint than unconsci…C13.11 RTM's first layer, read as spanda,predicts that changes i…
Figure 20. The mechanisms are settled; the claims about consciousness are not. Each claim from the ledger placed by my confidence grade. Notice that confidence is highest for claims about the measures themselves and lowest for claims about consciousness and RTM. Judgement
Text description of figure 20

Eleven rows, one per claim C13.1 to C13.11. C13.2, C13.3, C13.4, C13.6, C13.7 and C13.9 are established; C13.1, C13.5 and C13.8 likely; C13.10 and C13.11 speculative.

Glossary

Words used in this paper, in plain English

State space
The set of all possible states of a system, each written as a list of numbers.
Attractor
A region of state space that nearby states flow into and do not leave, such as a fixed point.
Basin
The set of states that flow into a particular attractor.
Bifurcation
A change in the number or kind of attractors as a control parameter is turned.
Fold (saddle-node) bifurcation
The point where a valley and a hill meet and vanish together, forcing a jump.
Hysteresis
Dependence of the current state on the path taken, shown as a loop when a control is swept up and down.
Neural inertia
Friedman and colleagues' term for the nervous system's tendency to resist a change between conscious and anaesthetised states.
Critical slowing down
Slower recovery from small disturbances as a valley flattens before a fold, seen as rising autocorrelation and variance.
Lag-one autocorrelation
How strongly each value in a time series predicts the next one.
Metastability
A regime in which parts of a system tend towards coordination without locking into it.
Order parameter (R)
In the Kuramoto model, the degree of agreement among oscillators, from 0 (scattered) to 1 (locked).
Criticality
Sitting at the boundary between two kinds of behaviour, where events of all sizes occur.
Branching ratio (m)
The average number of events each event triggers in the next moment; 1 is critical.
Subsampling
Recording only a fraction of the units in a system, which biases naive estimates of its dynamics.
Spanda (स्पन्द)
A slight movement or pulse; in the Spandakārikā, the activity of consciousness by which states arise and pass.
Unmeṣa (उन्मेष)
Opening, as of the eyes; the arising of a state or a new thought.
Nimeṣa (निमेष)
Closing, as of the eyes; the dissolving of a state.
Kṣobha (क्षोभ)
Agitation or disturbance of the mind.
Upalabdhṛ (उपलब्धृ)
The perceiver, the one who knows each state.
Prakāśa (प्रकाश)
Shining or showing; in RTM, the second layer, something appearing at all.
Conclusion

Kabir slows down#

Dynamical systems theory gives mind science a set of exact words, and each one makes a prediction. An attractor predicts that a pushed state returns, and that time in a state follows a law set by noise and by the depth of the basin. A fold with hysteresis predicts a loop that survives slow change. Critical slowing down predicts rising sluggishness before a fold, and it can be faked by changing noise. Metastability predicts swings in coordination that should beat an uncoupled control. Criticality predicts a branching ratio near one, and a naive estimate from a subsample will hide it.

Use each dynamical word as a prediction with a stated control: a slower sweep, a steady-noise check, an uncoupled population or a subsampling correction. Where the prediction is tested and passes, the word explains something. Where it is not tested, it is a description.

The Spandakārikā points to the same place from the other side. It tells the practitioner to watch the seam between one state and the next. The mathematics says that the seam is where the informative events are. Neither tells you, by itself, who is watching. The curve and the verse each stop short of that question, in different places. That is a question for the higher layers of RTM, and for later papers in this series.

Kabir does not get a second patient swept up and down a dozen times. What he gets, from Anjali, is a habit. At his next induction he raises the drug more slowly, and he writes down the brain monitor's pattern at both crossings, as well as the hand squeeze. At emergence he does the same. When the gap shrinks with the slower induction, he says nothing about attractors. When it does not, on a later list, he brings Anjali both sets of numbers.

"Now it might be inertia," she says. "Now you have something to test."

The habit generalises beyond theatre four. Whenever a paper about the brain uses a dynamical word, ask three questions. What signature does the word predict? What rival would produce the same data without it? And was the control run that would tell them apart? If all three have answers, the word is doing work. If not, it is a promising name for something not yet explained, and it should be read as one.

The whole argument on one pageSummary flow: four concepts above, each with a predicted signature, and four controls below, with the thesis in a band at the bottom.Four concepts, four signatures, four controlsAttractorPredicts return after a pushand exponential staysFold with hysteresisPredicts a loop thatsurvives slow changeCritical slowingPredicts rising sluggishnessbefore a foldMetastability,criticalityPredict swings and abranching ratio near 1Control: a timetablePeaked stays; no response tonoiseControl: a lagLoop shrinks to zeroControl: changing noiseWarning with no foldControl: uncoupled,subsampledFinite-size swings; biasedestimatesA dynamical word explains a transition only when its signature beats its control
Figure 21. The whole argument on one page. A one-page summary of the paper. Top row: what each concept predicts. Middle row: the control that could fake the prediction. What to notice: the bottom band is the thesis in one sentence. Summary (conceptual)
Text description of figure 21

Top row: attractor, fold with hysteresis, critical slowing, metastability and criticality, each with its signature. Middle row: the controls, a timetable, a lag, changing noise and an uncoupled or subsampled system. Bottom band: a dynamical word explains a transition only when its signature beats its control.

Further reading

Further reading, in a sensible order#

  1. Strogatz's account of the Kuramoto model, for the clearest introduction to coupled oscillators (Strogatz 2000).
  2. Haken, Kelso and Bunz, and the follow-up tests of critical slowing in finger movements, for a dynamical model that predicted its own signatures (Haken, Kelso and Bunz 1985; Scholz, Kelso and Schöner 1987).
  3. Scheffer and colleagues on early-warning signals, read with Ditlevsen and Johnsen on their limits (Scheffer et al. 2009; Ditlevsen and Johnsen 2010).
  4. Friedman and colleagues on neural inertia, for hysteresis in anaesthesia with drug lag ruled out (Friedman et al. 2010).
  5. Wilting and Priesemann on subsampling, and O'Byrne and Jerbi on the state of the criticality debate (Wilting and Priesemann 2018; O'Byrne and Jerbi 2022).
  6. Dyczkowski's translation of the Spandakārikā with its commentaries, for spanda in its own terms (Dyczkowski 1992).
Sources

Where every claim comes from

Each source was opened and checked against its publisher page, DOI record, arXiv record or a primary-text edition before it was cited. Anything I could not confirm is marked “to verify”.

Primary research

  1. Kramers, H. A. (1940). Brownian motion in a field of force and the diffusion model of chemical reactions. Physica 7(4), 284–304. https://doi.org/10.1016/S0031-8914(40)90098-2
  2. Kelso, J. A. S. (1984). Phase transitions and critical behavior in human bimanual coordination. American Journal of Physiology 246(6), R1000–R1004. https://doi.org/10.1152/ajpregu.1984.246.6.R1000
  3. Haken, H., Kelso, J. A. S., and Bunz, H. (1985). A theoretical model of phase transitions in human hand movements. Biological Cybernetics 51, 347–356. https://doi.org/10.1007/BF00336922
  4. Kelso, J. A. S., Scholz, J. P., and Schöner, G. (1986). Nonequilibrium phase transitions in coordinated biological motion: critical fluctuations. Physics Letters A 118(6), 279–284. https://doi.org/10.1016/0375-9601(86)90359-2
  5. Scholz, J. P., Kelso, J. A. S., and Schöner, G. (1987). Nonequilibrium phase transitions in coordinated biological motion: critical slowing down and switching time. Physics Letters A 123(8), 390–394. https://doi.org/10.1016/0375-9601(87)90038-7
  6. Friedman, E. B., Sun, Y., Moore, J. T., et al. (2010). A Conserved Behavioral State Barrier Impedes Transitions between Anesthetic-Induced Unconsciousness and Wakefulness: Evidence for Neural Inertia. PLoS ONE 5(7), e11903. https://doi.org/10.1371/journal.pone.0011903
  7. Massimini, M., Ferrarelli, F., Huber, R., Esser, S. K., Singh, H., and Tononi, G. (2005). Breakdown of Cortical Effective Connectivity During Sleep. Science 309, 2228–2232. https://doi.org/10.1126/science.1117256
  8. Scheffer, M., Bascompte, J., Brock, W. A., et al. (2009). Early-warning signals for critical transitions. Nature 461, 53–59. https://doi.org/10.1038/nature08227
  9. Ditlevsen, P. D., and Johnsen, S. J. (2010). Tipping points: Early warning and wishful thinking. Geophysical Research Letters 37, L19703. https://doi.org/10.1029/2010GL044486
  10. Boettiger, C., and Hastings, A. (2012). Quantifying limits to detection of early warning for critical transitions. Journal of the Royal Society Interface 9, 2527–2539. https://doi.org/10.1098/rsif.2012.0125
  11. Tognoli, E., and Kelso, J. A. S. (2014). The Metastable Brain. Neuron 81(1), 35–48. https://doi.org/10.1016/j.neuron.2013.12.022
  12. Rabinovich, M., Huerta, R., and Laurent, G. (2008). Transient Dynamics for Neural Processing. Science 321, 48–50. https://doi.org/10.1126/science.1155564
  13. Strogatz, S. H. (2000). From Kuramoto to Crawford: exploring the onset of synchronization in populations of coupled oscillators. Physica D 143, 1–20. https://doi.org/10.1016/S0167-2789(00)00094-4
  14. Seelig, J. D., and Jayaraman, V. (2015). Neural dynamics for landmark orientation and angular path integration. Nature 521, 186–191. https://doi.org/10.1038/nature14446
  15. Kim, S. S., Rouault, H., Druckmann, S., and Jayaraman, V. (2017). Ring attractor dynamics in the Drosophila central brain. Science 356, 849–853. https://doi.org/10.1126/science.aal4835
  16. Deco, G., Jirsa, V. K., and McIntosh, A. R. (2011). Emerging concepts for the dynamical organization of resting-state activity in the brain. Nature Reviews Neuroscience 12, 43–56 (online 2010). https://doi.org/10.1038/nrn2961
  17. Breakspear, M. (2017). Dynamic models of large-scale brain activity. Nature Neuroscience 20, 340–352. https://doi.org/10.1038/nn.4497
  18. Varela, F., Lachaux, J.-P., Rodriguez, E., and Martinerie, J. (2001). The brainweb: Phase synchronization and large-scale integration. Nature Reviews Neuroscience 2, 229–239. https://doi.org/10.1038/35067550
  19. Bak, P., Tang, C., and Wiesenfeld, K. (1987). Self-organized criticality: An explanation of the 1/f noise. Physical Review Letters 59, 381–384. https://doi.org/10.1103/PhysRevLett.59.381
  20. Beggs, J. M., and Plenz, D. (2003). Neuronal Avalanches in Neocortical Circuits. Journal of Neuroscience 23(35), 11167–11177. https://doi.org/10.1523/JNEUROSCI.23-35-11167.2003
  21. Chialvo, D. R. (2010). Emergent complex neural dynamics. Nature Physics 6, 744–750. https://doi.org/10.1038/nphys1803
  22. Wilting, J., and Priesemann, V. (2018). Inferring collective dynamical states from widely unobserved systems. Nature Communications 9, 2325. https://doi.org/10.1038/s41467-018-04725-4
  23. Touboul, J., and Destexhe, A. (2017). Power-law statistics and universal scaling in the absence of criticality. Physical Review E 95, 012413. https://doi.org/10.1103/PhysRevE.95.012413
  24. Carhart-Harris, R. L., Leech, R., Hellyer, P. J., et al. (2014). The entropic brain: a theory of conscious states informed by neuroimaging research with psychedelic drugs. Frontiers in Human Neuroscience 8, 20. https://doi.org/10.3389/fnhum.2014.00020
  25. Toker, D., Pappas, I., Lendner, J. D., et al. (2022). Consciousness is supported by near-critical slow cortical electrodynamics. Proceedings of the National Academy of Sciences 119(7), e2024455119. https://doi.org/10.1073/pnas.2024455119
  26. Fontenele, A. J., de Vasconcelos, N. A. P., Feliciano, T., et al. (2019). Criticality between Cortical States. Physical Review Letters 122, 208101. https://doi.org/10.1103/PhysRevLett.122.208101
  27. O'Byrne, J., and Jerbi, K. (2022). How critical is brain criticality? Trends in Neurosciences 45(11), 820–837. https://doi.org/10.1016/j.tins.2022.08.007
  28. van Gelder, T. (1998). The dynamical hypothesis in cognitive science. Behavioral and Brain Sciences 21(5), 615–628. https://doi.org/10.1017/S0140525X98001733

Primary texts and translations

  1. Vasugupta (or Bhaṭṭa Kallaṭa). Spandakārikā. Sanskrit text on GRETIL, input by Marino Faliero (not proof-read); verse numbers in this paper follow this file. https://gretil.sub.uni-goettingen.de/gretil/corpustei/transformations/html/sa_vasugupta-spandakArikA.htm
  2. Vasugupta. Spandakārikā with the Spandanirṇaya of Kṣemarāja. Sanskrit text on GRETIL, based on the edition by M. K. Shastri (Srinagar, 1925), input by Oliver Hellwig. https://gretil.sub.uni-goettingen.de/gretil/corpustei/sa_vasugupta-spandakArikA-comm.xml

Secondary sources and leads

  1. Dyczkowski, M. S. G. (1992). The Stanzas on Vibration: The SpandaKārikā with Four Commentaries. Albany: State University of New York Press. https://www.jstor.org/stable/jj.18255324
Cite this work

How to cite this paper

Mahapatra, R. R. (2026). Dynamical systems for cognition: attractors, bifurcations and metastability. Mind Frontier, paper 13, version 1.0. Rajesh Ranjan Mahapatra. https://rajeshrmahapatra.com/papers/mind-frontier/U0214-dynamical-systems-for-cognition.html (DOI to be assigned.)

@misc{mahapatra2026u0214,
  author       = {Mahapatra, Rajesh Ranjan},
  title        = {Dynamical systems for cognition: attractors, bifurcations and metastability},
  year         = {2026},
  howpublished = {Mind Frontier, paper 13},
  note         = {Version 1.0},
  doi          = {10.XXXX/placeholder},
  url          = {https://rajeshrmahapatra.com/papers/mind-frontier/U0214-dynamical-systems-for-cognition.html}
}
Top