Ehrenfest regularization of Hamiltonian systems

Imagine a freely rotating rigid body. The body has three principal axes of rotation. It follows from mathematical analysis of the evolution equations that pure rotations around the major and minor axes are stable while rotation around the middle axis is unstable. However, only rotation around the major axis (with highest moment of inertia) is stable in physical reality (as demonstrated by the unexpected change of rotation of the Explorer 1 probe). We propose a general method of Ehrenfest regularization of Hamiltonian equations by which the reversible Hamiltonian equations are equipped with irreversible terms constructed from the Hamiltonian dynamics itself. The method is demonstrated on harmonic oscillator, rigid body motion (solving the problem of stable minor axis rotation), ideal fluid mechanics and kinetic theory. In particular, the regularization can be seen as a birth of irreversibility and dissipation. In addition, we discuss and propose discretizations of the Ehrenfest regularized evolution equations such that key model characteristics (behavior of energy and entropy) are valid in the numerical scheme as well

Entropy and Entropy Production in Multiscale Dynamics

Heat conduction is investigated on three levels: equilibrium, Fourier, and Cattaneo. The Fourier level is either the point of departure for investigating the approach to equilibrium or the final stage in the investigation of the approach from the Cattaneo level. Both investigations bring to the Fourier level an entropy and a thermodynamics. In the absence of external and internal influences preventing the approach to equilibrium the entropy that arises in the latter investigation is the production of the classical entropy that arises in the former investigation. If the approach to equilibrium is prevented, then the entropy that arises in the investigation of the approach from the Cattaneo level to the Fourier level still brings to the Fourier level the entropy and the thermodynamics even if the classical entropy and the classical thermodynamics is absent. We also note that vanishing total entropy production as a characterization of equilibrium state is insufficient.Comment: Submitted to the Journal of Non-equilibrium Thermodynamic

Heterogeneity induces spatiotemporal oscillations in reaction-diffusions systems

We report on a novel instability arising in activator-inhibitor reaction-diffusion (RD) systems with a simple spatial heterogeneity. This instability gives rise to periodic creation, translation, and destruction of spike solutions that are commonly formed due to Turing instabilities. While this behavior is oscillatory in nature, it occurs purely within the Turing space such that no region of the domain would give rise to a Hopf bifurcation for the homogeneous equilibrium. We use the shadow limit of the Gierer-Meinhardt system to show that the speed of spike movement can be predicted from well-known asymptotic theory, but that this theory is unable to explain the emergence of these spatiotemporal oscillations. Instead, we numerically explore this system and show that the oscillatory behavior is caused by the destabilization of a steady spike pattern due to the creation of a new spike arising from endogeneous activator production. We demonstrate that on the edge of this instability, the period of the oscillations goes to infinity, although it does not fit the profile of any well known bifurcation of a limit cycle. We show that nearby stationary states are either Turing unstable, or undergo saddle-node bifurcations near the onset of the oscillatory instability, suggesting that the periodic motion does not emerge from a local equilibrium. We demonstrate the robustness of this spatiotemporal oscillation by exploring small localized heterogeneity, and showing that this behavior also occurs in the Schnakenberg RD model. Our results suggest that this phenomenon is ubiquitous in spatially heterogeneous RD systems, but that current tools, such as stability of spike solutions and shadow-limit asymptotics, do not elucidate understanding. This opens several avenues for further mathematical analysis and highlights difficulties in explaining how robust patterning emerges from Turing's mechanism in the presence of even small spatial heterogeneity

From one pattern into another: analysis of Turing patterns in heterogeneous domains via WKBJ

Pattern formation from homogeneity is well-studied, but less is known concerning symmetry-breaking instabilities in heterogeneous media. It is nontrivial to separate observed spatial patterning due to inherent spatial heterogeneity from emergent patterning due to nonlinear instability. We employ WKBJ asymptotics to investigate Turing instabilities for a spatially heterogeneous reaction-diffusion system, and derive conditions for instability which are local versions of the classical Turing conditions We find that the structure of unstable modes differs substantially from the typical trigonometric functions seen in the spatially homogeneous setting. Modes of different growth rates are localized to different spatial regions. This localization helps explain common amplitude modulations observed in simulations of Turing systems in heterogeneous settings. We numerically demonstrate this theory, giving an illustrative example of the emergent instabilities and the striking complexity arising from spatially heterogeneous reaction-diffusion systems. Our results give insight both into systems driven by exogenous heterogeneity, as well as successive pattern forming processes, noting that most scenarios in biology do not involve symmetry breaking from homogeneity, but instead consist of sequential evolutions of heterogeneous states. The instability mechanism reported here precisely captures such evolution, and extends Turing’s original thesis to a far wider and more realistic class of systems

Hierarchical patterning modes orchestrate hair follicle morphogenesis

Two theories address the origin of repeating patterns, such as hair follicles, limb digits, and intestinal villi, during development. The Turing reaction–diffusion system posits that interacting diffusible signals produced by static cells first define a prepattern that then induces cell rearrangements to produce an anatomical structure. The second theory, that of mesenchymal self-organisation, proposes that mobile cells can form periodic patterns of cell aggregates directly, without reference to any prepattern. Early hair follicle development is characterised by the rapid appearance of periodic arrangements of altered gene expression in the epidermis and prominent clustering of the adjacent dermal mesenchymal cells. We assess the contributions and interplay between reaction–diffusion and mesenchymal self-organisation processes in hair follicle patterning, identifying a network of fibroblast growth factor (FGF), wingless-related integration site (WNT), and bone morphogenetic protein (BMP) signalling interactions capable of spontaneously producing a periodic pattern. Using time-lapse imaging, we find that mesenchymal cell condensation at hair follicles is locally directed by an epidermal prepattern. However, imposing this prepattern’s condition of high FGF and low BMP activity across the entire skin reveals a latent dermal capacity to undergo spatially patterned self-organisation in the absence of epithelial direction. This mesenchymal self-organisation relies on restricted transforming growth factor (TGF) β signalling, which serves to drive chemotactic mesenchymal patterning when reaction–diffusion patterning is suppressed, but, in normal conditions, facilitates cell movement to locally prepatterned sources of FGF. This work illustrates a hierarchy of periodic patterning modes operating in organogenesis