The starting point is a practice that has become standard in climate modeling over the past two decades. Large ensembles of simulations, typically tens to hundreds of members, are run with identical forcing scenarios but slightly different initial conditions. Because the climate system is chaotic, tiny differences in where a simulation starts grow exponentially, so that after some time the ensemble members scatter across the range of states the system permits. The ensemble mean is then treated as the climatological mean, the ensemble spread as internal variability, and the time evolution of these statistics as the forced response to whatever forcing scenario was imposed. Projects such as the CESM Large Ensemble and the Max Planck Institute Grand Ensemble have made this approach a pillar of modern climate science, and its popularity is only expected to grow.

But there is an assumption buried inside this procedure, and Drótos and Bódai set out to expose it. Implicitly, researchers assume that the distribution of ensemble members eventually forgets how the members were initialized and converges to something objective, a probability distribution that belongs to the system and its forcing rather than to the personal choices of whoever set up the experiment. If that assumption fails, every modeling group running the same model under the same scenario could end up defining its own private climate, a subjective object that would make comparisons between studies ambiguous. The authors call this the uniqueness criterion: a satisfactory definition of climate must rest on a probability measure that is unique, depending only on the system, the forcing, and other objective factors.

The mathematical machinery that delivers uniqueness comes from dynamical systems theory. The Earth system, and the global models that describe it, can be viewed as nonautonomous dissipative deterministic dynamical systems, in which energy is dissipated and the equations depend explicitly on time through the forcing. Under a time-dependent forcing, the appropriate geometric object is the snapshot or pullback attractor, the set of states the system can occupy at each instant, equipped with a natural probability measure. Earlier work by Drótos, Bódai, and collaborators, notably a landmark 2015 paper, showed that an ensemble of trajectories converges to this natural measure in an approximately exponential fashion, so that after a finite convergence time the ensemble represents the same distribution regardless of how it was initialized. In an intermediate-complexity climate model, that convergence time turned out to be a few decades. The naive proposal, then, is to define climate at any instant as the natural probability measure of the snapshot attractor at that instant, with climate change being the evolution of this measure under the forcing.

Elegant as it is, this naive definition runs into a serious caveat, and dealing with it is the central novelty of the new paper. The convergence process itself has multiple characteristic time scales, because the loss of memory about initial conditions decomposes into a sum of exponentially decaying contributions. In the spectral theory of transfer operators, these decay rates are given by the real parts of the eigenvalues of the Ruelle–Perron–Frobenius operator, while the imaginary parts describe predictable, oscillatory evolution. Some modes of the system forget their initial conditions within years, others within decades, and some, such as those tied to the deep ocean, may take centuries or even up to a thousand years. If a study targets the coming century, waiting for the slowest modes to converge is pointless: their unpredictable variability would swamp the analysis in what the authors, citing David Stainforth, connect to the so-called Edinburgh paradox, the explosion of uncertainty when all possible long-term behaviors are included.

The resolution proposed by Drótos and Bódai is a conditional definition of climate. If there is a sufficiently large gap in the spectrum of convergence time scales, one can demand complete convergence in the faster-converging modes while conditioning on the actual state of the slower-converging ones. Climate is then a unique probability measure, but one that is conditional on a realization of the slow variables, which the authors call the predictable context. Crucially, the slow variables need not even be identified explicitly: initializing an ensemble by small perturbations of a model state automatically yields the desired conditional measure once the fast modes have decayed. The context itself is objective, because by definition the slow modes remain predictable over the time span of interest, so their state can in principle be learned from observations. This is what rescues uniqueness without pretending that the deep ocean’s memory can be ignored or averaged away.

The definition has consequences that reach beyond philosophy. One is that climate change and forced response, usually treated as synonyms, can come apart. If the slow modes evolve predictably and influence the converged statistics of the fast ones, the conditional climate changes even without any explicit forcing in the equations of motion. Such a change is genuine climate change under the conditional definition, but it is not, or not entirely, a forced response, since it does not originate from external time dependence. Disentangling the two requires comparing against an unforced evolution, in a spirit analogous to the removal of spurious model drift when estimating forced trends. Another consequence concerns initialization: for climate projections, the authors argue, the state of the slower-converging modes should be taken from observations, whereas much current practice samples arbitrary time instants of a long control run, which may be problematic regardless of how climate is defined.

Does the real Earth system cooperate? The authors’ preliminary assessment is cautiously optimistic. Ocean response studies suggest a possible separation by roughly a factor of ten between the time scales of the mixed layer and the deeper layers, and recent work indicates that convergence associated with the Atlantic Meridional Overturning Circulation takes up to about forty years. Taken together, these findings hint that for investigations spanning around a century, climate might be meaningfully defined through a probability measure obtained after a convergence time of a few decades, perhaps up to four. But open questions remain. Observed scaling behavior in climate time series, the possible role of multidecadal oscillations, regime transitions, and intertwined basins of attraction could all complicate the picture. The authors illustrate the regime problem with a stochastic slow-fast toy model: an ensemble initialized during a transition between regimes fails to converge on the fast time scale, destroying uniqueness, whereas initialization away from transitions preserves it.

To make the framework testable, the paper proposes a concrete initialization scheme for Earth system models. Pairs of ensembles are generated, the second initialized from a member of the first after a controlled delay, and the delays are increased across successive pairs sampled from different epochs of a control run. By checking whether the perturbed ensemble converges to its parent within the delay, researchers can determine the longest time span over which a practically unique, and therefore objectively definable, climate exists for a given variable, without ever needing to compute the full spectrum of convergence time scales. The authors are careful about limits: their conclusions are qualitative, the exponential-like convergence they confirmed numerically leaves room for model-specific detail, and if the required separation of time scales proves too small, the notion of climate may have to remain subjective, treated more like probabilistic ensemble weather forecasting. Even then, the framework offers a recipe: evaluate any statistical quantifier with respect to the converged ensemble, rather than inventing ensemble-based statistics one by one, a point connected to the violation of Birkhoff’s ergodic theorem in systems with explicit time dependence. What emerges is not merely a technical fix but a conceptual foundation, one that could change how the next generation of large ensemble experiments is designed, initialized, and interpreted.

Subject of Research: A dynamical-systems framework for defining climate via ensemble convergence and time scales of convergence

Article Title: Can we define climate by means of an ensemble? A tale of time scales of convergence

Article References: Drótos, G., & Bódai, T. (2026). Can we define climate by means of an ensemble? A tale of time scales of convergence. Earth System Dynamics, 17(5), 1529-1549. https://doi.org/10.5194/esd-17-1529-2026

Image Credits: AI Generated

DOI: 10.5194/esd-17-1529-2026

Keywords: climate definition, ensemble simulations, chaos theory, snapshot attractor, pullback attractor, internal variability, convergence time scales, Ruelle–Perron–Frobenius operator, Earth system models, forced response, climate change, dynamical systems

To understand why this matters, it helps to know what gerophysics is. The term describes an emerging interdisciplinary field that applies concepts from physics, mathematics, and complex systems science to the study of aging. Rather than focusing solely on molecular damage or cellular senescence, gerophysics borrows tools such as dynamical systems theory, entropy, network theory, and stochastic processes to describe how biological systems change over time. A recent report from the Global Conference on Gerophysics, discussed by Unfried and colleagues in the work that prompted this commentary, outlined a framework intended to connect aging processes across biological scales, from individual cells to the whole organism. Temprado and Sleimen-Malkoun argue that this framework, ambitious as it is, remains incomplete without an explicit account of behavior.

The reasoning behind that argument rests on a decades-long insight from the study of living systems: healthy physiological and behavioral systems do not produce perfectly regular signals, nor do they produce pure noise. Instead, they generate complex fluctuations that unfold across multiple time scales. A healthy heartbeat, for example, is not metronomic; it varies in subtle, structured ways. The same is true of posture, gait, and even the timing of daily activity. This structured variability is not a defect. It reflects the flexibility an organism needs to respond to an ever-changing environment. When complexity declines, the authors suggest, patterns become either more predictable or more random, and either shift may signal a reduced capacity for adaptation.

There is already substantial evidence that such complexity measurements are relevant to aging. Studies of brain activity, muscular signals, and movement variability have all shown that the richness of these signals tends to change as people grow older. Previous research has even explored whether alterations in movement patterns can help predict future falls, one of the most consequential health risks in later life. The commentary’s authors argue that these scattered findings should be viewed as interconnected aspects of a single phenomenon rather than isolated observations. If aging erodes the adaptive flexibility of the body, that erosion should leave fingerprints everywhere, from neural oscillations to the way a person walks across a room.

One of the most intriguing points in the commentary concerns how complexity changes across different biological levels. The authors discuss preliminary, unpublished findings suggesting that network entropy may decrease with age even as molecular entropy increases. If that pattern is confirmed by future research, it would carry a major implication: complexity at one biological level cannot necessarily be inferred from measurements at another. A blood test that captures molecular entropy might say nothing about the dynamical richness of a person’s movement or brain rhythms. Behavioral complexity, in other words, may carry information that molecular biomarkers alone cannot provide, which is precisely why the authors want it built into the gerophysics framework from the start.

From this line of thinking emerges the commentary’s most distinctive proposal: the concept of “complexity reserve.” Temprado and Sleimen-Malkoun define it as the residual buffer of dynamical flexibility available before a biological system crosses into irreversible functional decline. The idea extends existing notions of cognitive reserve and physical reserve, which describe how some individuals tolerate brain or body damage better than others, but it shifts the focus to how behavior changes dynamically in response to challenges. Rather than asking how much function a person has at a single moment, complexity reserve asks how much adaptive capacity remains in the system before it tips past a critical threshold.

How might such a reserve actually be measured? The authors sketch two complementary approaches. The first would examine how well individuals maintain organized behavioral patterns when exposed to progressively demanding tasks or environmental changes, effectively stress-testing the system’s flexibility in a controlled way. The second would estimate how close a person’s current functional state is to a transition point at which even a relatively small stressor could disrupt stability. This second approach echoes ideas from dynamical systems theory, in which systems approaching a critical transition often show characteristic warning signs. Either method, if validated, could turn an abstract theoretical concept into something a clinician could assess.

The potential implications for frailty are what make the proposal clinically exciting. Frailty is a condition characterized by reduced physiological reserve and increased vulnerability to stress, and it often precedes falls, hospitalization, and loss of independence. Current assessments typically rely on isolated health measurements taken at a single point in time. The commentary suggests an alternative: continuous behavioral monitoring that could identify declining adaptability before more obvious functional deterioration occurs. Instead of waiting for a person to become visibly frail, clinicians might one day detect the erosion of behavioral complexity early enough to intervene. The authors also raise the possibility that past experiences shape resilience, and that exposure to stimulating environments and varied physical or cognitive challenges could help maintain behavioral flexibility. Complexity reserve might even help explain why individuals respond so differently to exercise programs or environmental enrichment, though the authors are careful to note that these possibilities remain hypotheses requiring empirical validation.

Another forward-looking direction involves a phenomenon known as “complexity matching,” the idea that interactions between systems can be influenced by the similarity of their temporal patterns. When two people converse, for instance, the rhythmic structure of their speech and gestures can become coupled. The authors propose investigating whether exposure to stimuli with fractal-like structures, patterns that repeat similarly across many scales, could help restore aspects of age-related behavioral complexity. They are explicit, however, that this remains a theoretical possibility, not an established intervention for reversing aging-related decline. The distinction matters in a field where hype often outruns evidence, and the commentary’s authors are unusually candid about the limits of their own proposal.

Indeed, intellectual honesty is a defining feature of the piece. The commentary is conceptual rather than experimental. It presents no new clinical trial results and does not demonstrate that behavioral complexity measurements can reliably predict frailty or improve health outcomes. The proposed relationships between behavioral complexity, biological resilience, and responses to interventions will require carefully designed studies, including longitudinal measurements that track the same individuals over time and standardized assessments of adaptive behavior. The authors also caution against assuming that different measures of complexity represent the same underlying biological process. Establishing how molecular, network, and behavioral dynamics relate to one another will be essential before these concepts can be translated into practical assessments or interventions. In science, complexity is an easy word to use and a hard thing to pin down.

Even with those caveats, the commentary represents a potentially valuable addition to the growing field of gerophysics. By integrating measurements of how organisms respond to environmental challenges with molecular and physiological models of aging, researchers may gain a more complete picture of resilience and functional decline than either approach can offer alone. The framework opens new, testable directions for studying how adaptability changes with age, and it suggests that the road to healthier aging may run not only through blood tests and brain scans, but through the rich, fluctuating patterns of everyday behavior. If the authors are right, the way we move through the world may quietly announce, long before symptoms appear, how much adaptive life we have left.

Subject of Research: Behavioral complexity and gerophysics as markers of aging and frailty

Article Title: Behavioral complexity could offer new insights into aging and frailty

Article References: Behavioral complexity could offer new insights into aging and frailty. (n.d.). Original publication

Image Credits: AI Generated

DOI: Not provided

Keywords: aging, frailty, behavioral complexity, gerophysics, dynamical systems, entropy, complexity reserve, movement variability, biomarkers, resilience, physiological rhythms, complexity matching

The central problem the researchers set out to address is a gap in how climate tipping points are usually studied. Earth System Models, the vast numerical frameworks that simulate the climate at fine spatial scales, do project tipping behaviour in future scenarios. Yet these models are computationally so demanding that they have been run over only a narrow range of emissions pathways. Very few simulations exist in which warming is deliberately reversed, the so-called overshoot scenarios in which temperatures rise past a threshold and then decline. That scarcity leaves scientists with a limited understanding of hysteresis, the phenomenon in which a system, once tipped, refuses to return to its original state even after the forcing that triggered the change has been undone. Without such understanding, policymakers hoping to stabilise the climate after a temporary overshoot cannot know which damages would be reversible and which would be locked in.

Huntingford and colleagues turned to the mathematics of nonlinear dynamical systems, which have been refined over decades and are well suited to describing abrupt changes of state. Their starting point is a cubic equation of a form long used to describe large-scale environmental systems, including the Atlantic Meridional Overturning Circulation. The equation contains a bifurcation parameter, a quantity that in this case is set by the amount of global warming since pre-industrial times. As warming increases, the parameter moves toward a fold bifurcation, the mathematical point at which the stable state the system has occupied disappears and the system is forced to jump to an alternative equilibrium. Crucially, if the parameter only briefly exceeds that point, tipping may be avoided, a behaviour that depends on the inertia of the system, its tendency to respond slowly over long timescales.

The genuine novelty of the new paper lies in the algebra that maps real, measurable attributes of a climate component onto the equation’s abstract parameters. The framework requires just five quantities: the level of global warming at which tipping occurs, the extent or magnitude of the system at the moment of tipping, the lower temperature at which hysteresis ends and the system can return to its earlier state, the system’s extent at that lower temperature, and a single parameter describing inertia. From these five user-defined values, the authors derive closed-form expressions for the four coefficients of the governing equation. This means that anyone with knowledge of a system’s tipping threshold and its hysteresis behaviour, whether drawn from palaeoclimate records or from complex model output, can calibrate the simple equation directly, without any fitting procedure.

To drive the model, the team constructed a warming trajectory that begins with the smoothed historical record of global temperatures from the NASA-GISS dataset, spanning 1880 to 2024, and then extends it smoothly into the future with a quadratic overshoot profile. The coefficients of that quadratic are constrained by the final value and the rate of change of the historical record, ensuring a seamless transition from observed warming to the idealised future, and by a user-chosen peak warming level. In their numerical example, the authors set a peak of 2.8 degrees Celsius above pre-industrial levels, well beyond the tipping threshold in their illustrative configuration, and then let temperatures decline back down. This construction allows the equation to be tested across the full arc of an overshoot: the approach to the threshold, the passage beyond it, and the long return.

The simulations reveal how decisively inertia shapes the outcome. With low inertia, the system tips as warming peaks, jumping to the alternative state and then tracing a full hysteresis loop as temperatures fall, only returning to its original branch once cooling drops below the lower fold. At intermediate inertia, the system makes an extensive excursion toward the new state but ultimately recovers, sliding back toward its initial condition without ever completing the jump. At very high inertia, the state variable barely moves at all. The physical intuition is straightforward: a sluggish system such as a massive ice sheet may simply not have time to respond before the forcing recedes, whereas a fast-responding system such as a coral reef, with little inertia, is far more vulnerable to a transient overshoot.

Going beyond numerical experiments, the authors performed a scale analysis by rewriting the equation in non-dimensional form. This transformation collapses the problem into a compact expression governed by three dimensionless parameter clusters that combine the system’s tipping attributes, its inertia, and the amplitude and curvature of the warming overshoot. From this form, the team recovered an inequality, consistent with earlier theoretical work by Ritchie and colleagues, that cleanly separates the cases in which an overshoot triggers full tipping and hysteresis from those in which the system escapes unscathed. For their illustrative parameters, the analysis shows that avoiding tipping requires the inertia parameter to exceed a critical value of roughly 27.8, in close agreement with the numerical simulations. The agreement between the analytical threshold and the computed trajectories is a satisfying validation of the framework.

The authors are careful about the limitations of their approach. The cubic structure of the equation, to some extent, predetermines the shape of the hysteresis loop, and real climate components may require perturbation terms or asymmetric basins of attraction. The model tracks a single state variable, whereas scientists worry that the activation of one tipping element could alter the timing of others, creating cascades of the kind explored in coupled network models. The framework also addresses only fold bifurcations, while tipping can also arise from oscillatory instabilities associated with Hopf bifurcations, or even without any bifurcation at all, through rate-induced tipping of the kind implicated in peatland fires. In cases where a simple equation cannot capture the dominant qualitative behaviour, the authors note that more complex models from the climate modelling hierarchy remain the appropriate tools.

Even with those caveats, the potential applications are considerable. Because the equation is computationally trivial to run, it can be forced across a far wider ensemble of warming pathways than any Earth System Model could ever explore, including the overshoot trajectories that are becoming central to climate policy debates. If the same framework is fitted to multiple Earth System Models, the resulting spread in parameter values offers a concise way to quantify why the models disagree about when tipping occurs, covering the warming threshold, the size of the jump, and the amount of cooling needed to escape hysteresis. There is also a path toward better early warning systems: adding high-frequency noise to the calibrated equation and studying how the system’s variability changes as tipping approaches could sharpen the statistical signals that researchers monitor in real climate data.

The timing of this work gives it particular resonance. Recent analyses suggest the planet may already be at or near the 1.5 degree threshold, meaning that any eventual stabilisation at that level is likely to occur only after a temporary overshoot, potentially enabled by carbon removal technologies. Whether such an overshoot is a survivable detour or a one-way door depends on the inertia and hysteresis characteristics of each tipping element, quantities that this new framework is designed to quantify from the evidence that already exists. By providing a complete, reproducible manual for mapping climate components onto a single dynamical equation, Huntingford and his colleagues have given the tipping points research community a tool that is simple enough to be used widely, yet rich enough to capture the full drama of a system that jumps, locks itself into a new state, and only lets go when the world has cooled far more than it warmed in the first place.

Subject of Research: A simple nonlinear dynamical system for representing climate tipping points, hysteresis and inertia under global warming overshoot scenarios

Article Title: A simple dynamical system for representing climate tipping points with hysteresis

Article References: Huntingford, C., Ritchie, P. D. L., & Clarke, J. (2026). A simple dynamical system for representing climate tipping points with hysteresis. Nonlinear Processes in Geophysics, 33(3), 385-399. https://doi.org/10.5194/npg-33-385-2026

Image Credits: AI Generated

DOI: 10.5194/npg-33-385-2026

Keywords: climate tipping points, hysteresis, dynamical systems, bifurcation, Earth System Models, global warming overshoot, AMOC, system inertia, nonlinear processes, early warning signals, Paris Agreement, fold bifurcation

The difficulty stems from the sheer number of factors that shape collective dynamics. The behavior of any single component in a network depends not only on its own internal dynamics but also on how many connections it has, how strong those connections are, what kinds of interactions its neighbors exert upon it, and how long it takes for information or influence to propagate through the network. Structural complexity and time delays are known to matter individually, but how they combine to determine the boundary between stability and oscillation has remained poorly understood, particularly for large systems whose architecture resists straightforward analysis.

A new study published in National Science Research addresses this gap with an analytical framework that brings these ingredients together in a single, tractable description. The work shows how network complexity, propagation delays and interaction types jointly determine the transition between two contrasting dynamical regimes: amplitude death, in which coupling suppresses oscillations and the network settles into a steady state, and sustained oscillations, in which rhythmic activity persists indefinitely across the system. By deriving explicit relationships among these quantities, the framework turns a question that previously could only be probed case by case through simulation into one that can be answered analytically.

The mathematical backbone of the study is the network of coupled Stuart-Landau oscillators, a widely used and well-understood model for oscillatory dynamics. Each oscillator in such a network behaves, in isolation, like a simple limit-cycle system whose amplitude and phase evolve according to well-characterized equations. When many such oscillators are coupled together, the collective behavior depends delicately on the coupling structure. This makes the Stuart-Landau framework an ideal testing ground: it is simple enough to permit rigorous analysis, yet rich enough to display the full range of collective phenomena observed in real networks, from complete synchronization to the complete suppression of activity known as amplitude death.

Working within this setting, the researchers derived a relationship between the effective complexity of the network and the critical delay at which a previously stable system begins to oscillate. The resulting picture is striking. Increasing network complexity generally reduces the amount of propagation delay needed to trigger oscillations, meaning that richer, more highly connected architectures are intrinsically closer to the oscillatory regime. In sufficiently complex networks, the analysis shows, sustained oscillations can emerge even in the complete absence of propagation delay. Complexity alone, in other words, can be enough to destabilize a steady state that would remain perfectly stable in a simpler network with identical components and coupling strengths.

The framework goes further by examining how the nature of the interactions between nodes reshapes this transition boundary. The researchers considered four interaction types: cooperative, in which connections reinforce one another; competitive, in which connections oppose one another; mixed, combining both; and random. Each type was found to modify the transition boundary in the plane spanned by complexity and delay. The ordering that emerges is clear and systematic. Cooperative interactions promote the onset of sustained oscillations most readily, requiring the lowest critical complexity level for rhythmic activity to appear. Competitive, mixed and random interactions demand progressively higher critical complexity levels before sustained oscillations can emerge, so the same network architecture that oscillates readily under cooperative coupling may remain steady under competitive coupling.

This interaction-dependent hierarchy carries practical implications for any field in which network design matters. In engineering contexts such as power grids or communication networks, where oscillations can be destructive, the results suggest that the sign and structure of interactions should be treated as a design parameter on equal footing with topology and delay management. In biological contexts, where cooperative interactions are common, the findings offer a possible explanation for why rhythmic behavior arises so readily in neural circuits, gene regulatory networks and ecological communities: the interaction structure itself may lower the barrier to oscillation, allowing rhythms to emerge without requiring long propagation delays or extreme architectural complexity.

A central strength of the study lies in its effort to verify that the predicted transitions survive contact with the physical world rather than existing only in idealized computation. To this end, the researchers constructed a digital-analog hardware emulation platform. A microcontroller updated the network dynamics in real time, while external electronic circuits converted selected network states into measurable voltage signals. This hardware-in-the-loop approach deliberately exposed the theoretical predictions to the imperfections of real instrumentation, including finite sampling rates, signal quantization, transistor switching and other implementation artifacts that are absent from purely numerical experiments.

The observed hardware transitions matched the predicted critical delays, providing evidence that the framework captures a robust physical phenomenon rather than a fragile artifact of floating-point arithmetic. The agreement across three independent levels of scrutiny, namely analytical derivation, numerical simulation and hardware emulation, indicates that the complexity-delay transition remains observable in systems subject to the sampling, quantization and noise inherent in practical implementations. For researchers who wish to apply these results to real engineered or biological systems, this robustness is arguably as important as the analytical results themselves, because real networks never satisfy idealized assumptions exactly.

To demonstrate applicability beyond synthetic architectures, the researchers further applied their framework to the connectome of the nematode worm Caenorhabditis elegans, one of the most completely mapped neural networks in biology. The analysis of this biologically derived topology illustrated how the theoretical tools can be carried over to empirical network data, opening a path toward assessing whether the oscillatory tendencies of real neural systems can be anticipated from their structural complexity, interaction types and signal propagation delays alone. Taken together, the study offers a unified lens on a question that touches neuroscience, ecology, epidemiology and engineering alike: when does a network hold steady, and when does it begin to sing? By showing that complexity, delay and interaction type jointly draw the boundary, and by validating that boundary in hardware and in a real connectome, the work provides both a conceptual map and a practical toolkit for navigating the transition between silence and rhythm in complex systems.

Subject of Research: Transitions between amplitude death and sustained oscillations in complex networks of coupled oscillators

Article Title: What shape the oscillatory transitions in complex networks?

Article References: What shape the oscillatory transitions in complex networks?. (n.d.). Original publication

Image Credits: AI Generated

DOI: Not provided

Keywords: complex networks, oscillations, amplitude death, Stuart-Landau oscillators, propagation delays, network complexity, cooperative interactions, competitive interactions, coupled oscillators, hardware emulation, C. elegans, dynamical systems

Unimodular gravity has long intrigued theorists precisely because of how it handles the cosmological constant problem, arguably the most notorious discrepancy between theory and observation in physics. Quantum field theory predicts a vacuum energy vastly larger than what astronomers measure, yet in unimodular gravity the situation changes fundamentally. By restricting the symmetry of Einstein’s theory so that transformations must preserve the four-dimensional volume element, the theory yields only the trace-free part of Einstein’s field equations. The cosmological constant then ceases to be a parameter fixed by the microscopic physics and instead appears as a constant of integration, much like the energy of a pendulum, whose value is determined by initial conditions rather than by quantum loops.

The price of this elegant restructuring is a modified conservation law. In standard general relativity, the energy-momentum tensor of matter is covariantly conserved, which forces the cosmological constant to be truly constant. In unimodular gravity, that constraint is relaxed, and the divergence of the energy-momentum tensor is balanced by the gradient of the Lagrange multiplier that enforces the volume restriction. Solving the resulting equation reveals that the cosmological term splits into two pieces: a fixed integration constant and an arbitrary function of time, which the authors call the energy diffusion function. This function quantifies exactly how much local energy conservation fails, and it acts as a continuous source or sink of energy for the cosmic fluids filling the universe.

Physically motivated guesses for this diffusion function have been proposed before. One striking example draws on the continuous spontaneous localization model of quantum collapse, in which energy is genuinely created during wavefunction collapse, suggesting a diffusion term proportional to the energy density of the fluid itself. Earlier work showed that when only dark matter diffuses, the resulting cosmology closely reproduces the standard Lambda-CDM model, and related diffusion scenarios have even been explored as a way to ease the Hubble tension, the persistent disagreement between different measurements of the universe’s expansion rate. The trouble, as the Chilean team emphasizes, is that no widely accepted diffusion function follows from an established physical process, so most proposals rest on phenomenological convenience or mathematical simplicity rather than principle.

The new study attacks this problem from the opposite direction. Rather than postulating a diffusion law and studying its consequences, the authors recast the cosmological equations as an autonomous dynamical system and ask what diffusion functions the phase-space structure itself demands. The key move is to introduce a dimensionless variable called the diffusion slope, defined as minus the logarithmic derivative of the diffusion function with respect to the e-fold number, a natural clock for cosmic expansion. For the system to close, this slope must be invertible along the trajectories — a condition satisfied whenever it evolves monotonically. Remarkably, when the slope is held constant, the formalism automatically spits out the power-law diffusion functions that had previously been introduced by hand, suggesting that the dynamical systems perspective can rediscover and organize known models rather than merely accommodate them.

With the autonomous system in hand, the team mapped out the fixed points of cosmic evolution and their stability. The familiar matter-dominated era appears as a saddle point, not an attractor, meaning the universe inevitably passes through it rather than lingering. The de Sitter solution, where expansion accelerates at a constant rate under the influence of the integration constant, emerges as a genuine late-time attractor for positive diffusion slopes. More intriguingly, the analysis uncovered a novel matter-diffusion scaling solution in which the diffusion term tracks the dark matter density, maintaining a constant fractional share of the cosmic energy budget. Within a specific range of parameters, this scaling regime itself drives accelerated expansion — entirely without a cosmological constant — although it represents a transient stage rather than a final destiny.

Perhaps the most provocative result concerns the purely diffusion-dominated configuration. When the diffusion function becomes constant, it behaves exactly like a cosmological constant, producing exponential expansion with an effective equation-of-state parameter of minus one. In that regime, the diffusion sector alone can power late-time acceleration, with no fundamental vacuum energy required at all. The stability analysis revealed subtlety here: linear theory alone cannot settle the fate of this point because one eigenvalue vanishes, so the researchers supplemented their analytic work with numerical integration of the phase-space flow, showing that trajectories are attracted along one direction but repelled along another. The diffusion-dominated solution is therefore a saddle, a waystation the cosmos may visit but not a permanent home — the true endpoint remains the de Sitter state governed by the integration constant.

Beyond the asymptotic regimes, the authors built a full reconstruction machinery analogous to potential reconstruction in scalar-field cosmology. By specifying a curvature function that controls how the diffusion slope bends in logarithmic space, one can integrate the slope evolution and then reconstruct the diffusion function itself through a simple exponential integral. The constant-curvature family already displays rich behavior: curvature equal to one recovers power-law diffusion, curvature greater than one drives the diffusion term smoothly to zero as the universe expands, and curvature below one produces a finite-time singularity in the slope, signaling a breakdown of the description. More elaborate curvature functions that cross unity allow the slope to station at multiple values, enabling trajectories that interpolate between a rapidly decaying diffusion contribution in the early universe and an asymptotically constant component at late times — precisely the kind of transition a viable dark energy candidate might need.

The framework also clarifies what thermodynamics demands. Using the Gibbs relation, the team showed that the second law of thermodynamics requires the diffusion function to be non-increasing as the universe expands, guaranteeing positive entropy production. Notably, this condition constrains only the trend, not the sign, of the diffusion term, which can push the effective cosmological constant either up or down. The reconstruction formalism respects this automatically: because the diffusion function is built from a strictly positive exponential factor, it cannot change sign dynamically, keeping every reconstructed model consistent with the thermodynamic constraint from the outset.

The implications reach well beyond formal elegance. Because the diffusion sector exchanges energy with dark matter, it should alter the growth of cosmic structure, leaving fingerprints in the matter power spectrum and in the growth rate of galaxies that could distinguish diffusion cosmologies from Lambda-CDM. The authors identify this perturbative analysis, together with a direct confrontation of the reconstructed models with observational data, as the natural next step of their program. If the coming generation of surveys continues to hint that dark energy is not perfectly constant — as recent measurements have suggested — then a framework that generates and classifies diffusion models from first principles, rather than by trial and error, may prove exactly what cosmologists need to make sense of a universe that refuses to sit still.

Subject of Research: Cosmological diffusion models and energy non-conservation in unimodular gravity, analyzed through dynamical systems reconstruction

Article Title: Unimodular gravity with arbitrary diffusion function: a dynamical system reconstruction approach

Article References: Unimodular gravity with arbitrary diffusion function: a dynamical system reconstruction approach. (n.d.). https://doi.org/10.1140/epjc/s10052-026-16363-y

Image Credits: AI Generated

DOI: 10.1140/epjc/s10052-026-16363-y

Keywords: unimodular gravity, cosmological constant problem, dark energy, energy diffusion function, dynamical systems, phase-space fixed points, de Sitter expansion, scaling solutions, entropy production, Lambda-CDM, Hubble tension, theoretical cosmology

The research, led by Yutong Cai and Zhaoyue Zhong, with Zefeng Zhang, Bo-Wei Qin and senior author Wei Lin at Fudan University in Shanghai, introduces a rigorous mathematical treatment of what the authors call modulatability. Rather than asking how to push an oscillator toward a target behavior through trial and error, the team asked a more fundamental question about the geometry of the problem itself. If you want to impose several properties on an oscillating system at once, say a specific frequency, a specific amplitude, and a specific phase relationship across parts of a network, how many adjustable knobs do you need? Their central theoretical result is strikingly elegant: the local uniqueness of a solution is generally guaranteed when the number of free parameters equals the number of properties being targeted. Too few knobs and the problem is overconstrained, with most target combinations simply unreachable. Too many and solutions proliferate, making it hard to know which one a solver will find.

This dimensional matching principle may sound abstract, but it has immediate practical consequences. In the biological sciences, researchers routinely try to reprogram oscillators whose dynamics involve dozens of interacting variables. Synthetic biologists, for example, have spent two decades engineering genetic circuits such as the repressilator, a ring of three genes that repress one another in turn and thereby generate oscillating protein concentrations. Tuning such circuits to display both a desired period and a desired amplitude simultaneously has historically required laborious computational redesign or experimental screening. The Fudan team’s framework reframes the whole enterprise: instead of simulating trajectories forward and hoping to stumble on good parameters, the modulation task is treated as a property-based inverse problem, in which the desired properties are fixed and the equations are solved backward for the parameters that produce them.

The technical machinery behind the framework combines several ideas from dynamical systems and numerical analysis. The oscillatory trajectories of the system are represented through a Fourier expansion, allowing periodic solutions to be described by a finite set of coefficients once a truncation order is chosen. The desired properties, such as minima of particular protein activities or the period of the collective rhythm, are expressed as functions of these coefficients and of the system’s kinetic parameters. Solving for the parameters then becomes a root-finding problem in a space whose dimension is chosen deliberately to match the number of targets. The authors deploy Newton-type iterations to converge on solutions, and they use continuation strategies, stepping gradually through parameter space, to follow modulation paths even when the targets are far from the system’s natural behavior. Sweeps over system dimension, Fourier truncation order, number of targets, iteration tolerance, and continuation steps, reported in the paper’s extended data, map out how computational cost scales with each factor.

Crucially, the framework is not merely a simulation exercise. The team validated it on electronic analogs of genetic circuits, building hardware implementations of the repressilator concept in which voltages play the role of protein concentrations. From recorded voltage time courses alone, the researchers inferred the underlying parameters of the oscillating circuit and then modulated them to hit orthogonal targets, meaning goals that do not interfere with one another, such as independently shifting the period while holding amplitude fixed. The experiments included orthogonal period modulation tasks with two and with three simultaneous targets, and the inferred parameter distributions were tight and unimodal, suggesting that the inverse problem is well-posed in precisely the sense the theory predicts. This hardware demonstration matters because electronic circuits are a standard proving ground for ideas intended eventually for living cells, capturing the nonlinear feedback structure of genetic networks while remaining far easier to measure and perturb.

The data-driven side of the work is where the framework arguably delivers its most compelling numbers. In benchmark comparisons against baseline inference methods, the property-based approach achieved shorter runtimes and higher accuracy, a combination that is rare in inverse problems, where speed is usually bought at the cost of precision. The authors also demonstrated the framework on a model of an engineered Sir2-HAP negative feedback loop, a genetic oscillator linked to cellular longevity in yeast. In that system, low activity of the SIR2 and HAP proteins corresponds to a detrimental aging zone that the oscillating state variables pass through. By identifying a parameter set that up-regulates the minimal activities of both proteins, the framework showed in stochastic simulations, which added Ornstein-Uhlenbeck noise to eight kinetic terms, that the modulated oscillator avoids the aging zone more effectively, with distributions of minimum protein levels shifted upward across one hundred independent runs. It is a concrete illustration of how abstract parameter identification could translate into a design principle for extending cellular lifespan.

Networks add another layer of complexity, and the paper addresses it head-on. In one extended demonstration, the team applied their method to a FitzHugh-Nagumo neuronal network of one hundred nodes arranged on a randomly generated Barabasi-Albert topology, a structure mimicking the heterogeneous, hub-rich connectivity seen in many real systems. With fifty target nodes whose amplitudes needed adjustment, the researchers modulated random subsets of edges and measured how often the modulation succeeded while the network remained oscillatory. The results, summarized across many repeated trials, revealed how success rates depend on the balance between the number of targeted nodes and the number of edges available for modulation, with the equal-dimension condition again emerging as the natural dividing line. For anyone trying to control synchronization patterns in power grids, neural tissue, or coupled laser arrays, this kind of systematic map of when control is feasible is exactly the missing piece.

What makes the study resonate beyond applied mathematics is the sheer breadth of oscillatory phenomena it touches. The references span circadian medicine, the bioelectrical phase transition that patterns the first vertebrate heartbeats, collective oscillations emerging in massive human crowds, amplitude and frequency modulation of subthalamic beta oscillations in Parkinson’s disease, and synthetic oscillators designed to slow cellular aging. In each of these domains, the practical question is the same: which properties of the rhythm can be moved, and by touching which parameters? The modulatability framework offers a principled answer, replacing intuition with a dimension-counting criterion and replacing forward simulation with inverse problem solving. It also connects to a broader theoretical conversation about controlling complex networks, complementing earlier work on controllability, functional control of oscillator networks, and inverse problems for dynamic patterns in coupled systems.

The authors have made their work unusually accessible. The source code is publicly available on GitHub, with a frozen release archived on Zenodo, and a figure-to-script mapping document allows readers to reproduce individual results directly. Source data, including network topologies, modulation parameters, oscillatory trajectories, stochastic realizations, and the experimental voltage recordings from the electronic circuit, are deposited without restrictions. For a field where reproducibility of nonlinear dynamics can be notoriously difficult, this level of transparency lowers the barrier for other groups to test the framework on their own oscillators, whether biological, neural, mechanical, or photonic.

The larger significance of the work lies in its reframing. For years, the modulation of oscillators has been pursued through two parallel traditions, one statistical and one rooted in dynamical systems theory, each producing partial answers. By proving a general condition for when modulation parameters exist uniquely and packaging the insight into an efficient computational pipeline, the Fudan team has given both theorists and experimentalists a common language for a problem that cuts across biology, neuroscience, and engineering. If rhythms govern health, and disrupted oscillations are implicated in conditions from diabetes to psychiatric illness, then knowing exactly when and how a rhythm can be retuned is more than a mathematical curiosity. It is a step toward making oscillators, whether in a cell, a circuit board, or a brain, into things we can genuinely engineer.

Subject of Research: Mathematical theory and computational framework for the modulatability of complex oscillators

Article Title: Modulatability of complex oscillators

Article References: Modulatability of complex oscillators. (n.d.). https://doi.org/10.1038/s43588-026-01050-5

Image Credits: AI Generated

DOI: 10.1038/s43588-026-01050-5

Keywords: oscillators, dynamical systems, inverse problems, complex networks, synchronization, genetic circuits, repressilator, circadian rhythms, parameter identification, nonlinear dynamics, Nature Computational Science, synthetic biology

The result matters because integrability is the dividing line between order and chaos in classical mechanics. An integrable system possesses enough constants of motion, one for each degree of freedom, to pin down every trajectory analytically. The two-body problem under gravity is the textbook example: Kepler’s laws follow from conserved energy, angular momentum and the direction of the orbital plane. Add a third body, even confined to a line, and the supply of conserved quantities runs dry. Poincaré recognized more than a century ago that the general three-body problem resists solution, and subsequent work has mapped out precisely which special cases escape this fate. The new study extends that map into the territory of charged particles, where electrostatic attraction and repulsion join gravity in shaping the dynamics.

In the charged three-body problem, the particles interact through a potential that combines inverse-square force terms, generalizing the purely gravitational case studied for centuries. When the charges and masses are chosen so that some particles attract and others repel, the motion along the line can be far richer than in the gravitational version, where all interactions pull inward. The authors focus on the equal-mass case, a natural benchmark that nonetheless captures the essential difficulty. Their central claim is rigorous: for equal masses, the system admits no additional analytic first integral beyond the obvious conserved energy, and therefore cannot be solved by quadrature.

The proof rests on a powerful modern tool known as Morales-Ramis theory, developed by Jesús Morales-Ruiz and Jean-Pierre Ramis at the turn of the millennium. The theory forges a link between the classical question of integrability and the differential Galois theory of linear equations. The strategy is to examine how nearby trajectories behave in the neighborhood of a particular solution of the system. If the full nonlinear system is integrable, then the linearized equations governing small perturbations around that special solution must themselves be solvable in a precise sense: their differential Galois group must be virtually Abelian. If the Galois group turns out to be too large, typically the full symplectic group, integrability is impossible.

The special solutions at the heart of the analysis are called central configurations. These are arrangements of the point masses, here specific positions along the line, at which the force on each particle is exactly proportional to its distance from the center of mass. Central configurations act as the skeleton of the three-body problem: they generate the famous Euler and Lagrange solutions and organize the topology of the phase space. Crucially for the authors’ purposes, the Hessian matrix of the potential function, its matrix of second derivatives, can be evaluated at each central configuration, and the eigenvalues of that matrix feed directly into the Morales-Ramis integrability criteria.

Here the analysis collides with a notorious algebraic obstacle. Finding the central configurations of the three-body problem on a line requires solving the Euler quintic, a fifth-degree polynomial whose roots correspond to the allowed arrangements. Fifth-degree polynomials generally cannot be solved by radicals, and the roots of the Euler quintic are unwieldy expressions that resist direct manipulation. Earlier approaches to non-integrability proofs often stalled at exactly this point, forced either to handle the roots symbolically in full generality or to restrict attention to special parameter values where the quintic factors nicely.

Shibayama and Takeguchi sidestep the problem with an elegant algebraic maneuver. Rather than computing the roots of the Euler quintic, they exploit the classical relations between roots and coefficients of a polynomial. The Morales-Ramis criteria do not actually require the eigenvalues themselves; they constrain the eigenvalues through their elementary symmetric polynomials, the building blocks that appear in Vieta’s formulas. The authors show that these symmetric polynomials, evaluated at the eigenvalues associated with each root of the quintic, can be expressed purely in terms of the physical parameters of the system, namely the masses and the coefficients governing the interaction terms in the potential. This reformulation converts an intractable root-finding problem into a tractable calculation with polynomial expressions in the parameters, and it allows the non-integrability conditions to be checked across the equal-mass case without ever writing down a single root.

The payoff is a proof that the eigenvalue configurations demanded by integrability cannot occur. The Morales-Ramis theory imposes strong arithmetic restrictions: the eigenvalues of the Hessian at a central configuration must satisfy rigid algebraic relations if an additional analytic integral is to exist. By expressing the relevant symmetric functions in terms of masses and interaction parameters, the authors demonstrate that for equal masses these relations fail, so the differential Galois group of the variational equations is large enough to rule out integrability. The conclusion is that the one-dimensional charged three-body problem with equal masses is genuinely non-integrable, no matter how the remaining interaction parameters are tuned within the family considered.

As a companion result, the paper offers an alternative proof of non-integrability for the purely gravitational case, revisiting territory that Shibayama had explored in earlier work on the collinear three-body problem. What distinguishes the new treatment is its use of symbolic computation. The authors employ the computer algebra system Maple to carry out the lengthy manipulations of the variational equations and the verification of the Galois-theoretic obstructions, providing a machine-checked pathway through calculations that would be error-prone by hand. The approach signals a broader trend in celestial mechanics, where computer-assisted algebra is increasingly deployed to certify results that were once the exclusive province of painstaking manual derivation, following a line of work by researchers such as Alin Bostan, Thierry Combot and Mohab Safey El Din on computing integrability conditions for parametrized potentials.

The study builds on and complements recent advances in the field. In 2025, Maria Przybylska and Andrzej Maciejewski published a non-integrability result for the charged three-body problem in the same journal, and the new work sharpens the picture by handling the equal-mass one-dimensional case with a method that tames the Euler quintic. Hiroshi Yoshida’s classic 1987 criterion for homogeneous potentials, along with the Morales-Ruiz and Simó analysis of n-body non-integrability, forms the theoretical backbone of the enterprise. Taken together, these results chart a program: identify the exact boundary between solvable and unsolvable regimes of few-body dynamics. The Kyoto team’s contribution shows that even when the geometry is reduced to a single line and the masses are made identical, the charged three-body problem remains on the unsolvable side of that boundary, a reminder that deterministic simplicity in the laws of motion does not guarantee solvability in the mathematics that follows from them.

Subject of Research: Non-integrability of the one-dimensional charged three-body problem via Morales-Ramis theory

Article Title: Non-integrability of some one-dimensional charged three-body problems

Article References: Shibayama, M., & Takeguchi, Y. (2026). Non-integrability of some one-dimensional charged three-body problems. Celestial Mechanics and Dynamical Astronomy, 138(5), Article 57. https://doi.org/10.1007/s10569-026-10332-z

Image Credits: AI Generated

DOI: 10.1007/s10569-026-10332-z

Keywords: three-body problem, non-integrability, Morales-Ramis theory, celestial mechanics, central configurations, Euler quintic, differential Galois theory, Hamiltonian systems, charged particles, symbolic computation, Kyoto University, dynamical systems