# American Institute of Mathematical Sciences

December  2019, 11(4): 447-485. doi: 10.3934/jgm.2019023

## Networks of coadjoint orbits: From geometric to statistical mechanics

 Department of Mathematics, Imperial College, London SW7 2AZ, UK

* Corresponding author: So Takao

To Darryl Holm, for his 70th birthday

Received  April 2018 Revised  May 2019 Published  November 2019

Fund Project: AA acknowledges EPSRC funding through award EP/N014529/1 via the EPSRC Centre for Mathematics of Precision Healthcare. ST acknowledges funding through Schrödinger scholarship scheme.

A class of network models with symmetry group $G$ that evolve as a Lie-Poisson system is derived from the framework of geometric mechanics, which generalises the classical Heisenberg model studied in statistical mechanics. We considered two ways of coupling the spins: one via the momentum and the other via the position and studied in details the equilibrium solutions and their corresponding nonlinear stability properties using the energy-Casimir method. We then took the example $G = SO(3)$ and saw that the momentum-coupled system reduces to the classical Heisenberg model with massive spins and the position-coupled case reduces to a new system that has a broken symmetry group $SO(3)/SO(2)$ similar to the heavy top. In the latter system, we numerically observed an interesting synchronisation-like phenomenon for a certain class of initial conditions. Adding a type of noise and dissipation that preserves the coadjoint orbit of the network model, we found that the invariant measure is given by the Gibbs measure, from which the notion of temperature is defined. We then observed a surprising 'triple-humped' phase transition in the heavy top-like lattice model, where the spins switched from one equilibrium position to another before losing magnetisation as we increased the temperature. This work is only a first step towards connecting geometric mechanics with statistical mechanics and several interesting problems are open for further investigation.

Citation: Alexis Arnaudon, So Takao. Networks of coadjoint orbits: From geometric to statistical mechanics. Journal of Geometric Mechanics, 2019, 11 (4) : 447-485. doi: 10.3934/jgm.2019023
##### References:

show all references

##### References:
In this figure, we plot the maximum real part of the eigenvalues $\lambda_S$ of the linearised system (58) corresponding to equilibrium solutions ${\boldsymbol \Pi}_e$ (eigenvectors of $\mathbb L$) with energy $\lambda_e$ (eigenvalue of $\mathbb L$ corresponding to eigenvector ${\boldsymbol \Pi}_e$). The blue crosses indicate the maximum real eigenvalues that are non-zero, corresponding to unstable solutions, and the red dots indicate those that are zero, corresponding to linearly stable solutions. The purple dots indicate the algebraic multiplicity $n$ of $\lambda_e$, whose values are to be read from the right axis. As expected, the equilibria corresponding to the highest and lowest energy are stable in both cases
In this figure, we display four non-trivial examples of linearly stable equilibrium solutions of the $20\times 20$ rigid body lattice. All the solutions are parallel to one of the axes (either $\Pi_1$ or $\Pi_3$) and have an 'argyle' pattern of aligned spins (2(a), 2(d)) or anti-aligned spins (2(b), 2(c)). The former pattern corresponds to configurations with low energy (negative values of $\lambda$), and the latter corresponds to configurations with high energy (positive values of $\lambda$). All four solutions given here correspond to the red dots in figure 1
This figure shows the possible $(\lambda_1, \lambda_2)$-pairs that one can have from solving equation (86). Fixing $\lambda_1$, the possible values of $\lambda_2$ that solve (86) are the eigenvalues of the graph Laplacian $\mathbb L(\lambda_1)$. We see that the eigenvalues $\lambda_2$ in the left panel 3(a) collapse as $\lambda_1 \rightarrow 0$, which does not occur in the right panel 3(b)
The two panels in this figure display the maximum, real part of the eigenvalues $\lambda_S$ of the linearised system (98) around the equilibrium configuration ${\boldsymbol \Gamma}_e$ (eigenvectors of $\mathbb L(\lambda_1)$) plotted against $\lambda_2$ (minus the eigenvalue of $\mathbb L(\lambda_1)$ corresponding to ${\boldsymbol \Gamma}_e$) for $\lambda_1 = 0.5$. Again, the blue crosses correspond to unstable equilibria, the red dots correspond to linearly stable equilibria and the purple dots correspond to the multiplicity of $\lambda_2$, whose values are to be read from the right axis. As expected, the configuration corresponding to the highest $\lambda_2$ (or, the lowest eigenvalue of $\mathbb L(\lambda_1)$) is stable in both cases
This figure shows the time evolution of the averaged spins (5(a), 5(c), 5(e)) and energy (5(b), 5(d), 5(f)) in a $20\times 20$ heavy top lattice, starting from different initial conditions, with dissipation (dashed lines) or without dissipation (bold lines). The blue, orange and green lines in the left panel correspond to the $\Gamma_1$, $\Gamma_2$ and $\Gamma_3$-components of the spins averaged over the lattice respectively and the same colours in the right panel correspond to the total, kinetic and potential energy. We observe metastability of the $\Gamma_1$ and $\Gamma_2$ axes in panels 5(c) and 5(a) in the dissipative case, and a partial synchronisation around the $\Gamma_3$-axis in the non-dissipative case in 5(a), 5(c) and 5(e)
The two panels in this figure show the temperature phase transition of the rigid body lattice for two different cases. The blue dots correspond to averaged magnetisation computed via direct simulation of the full stochastic equation (47) and the black line corresponds to the mean field approximation (103). Both cases exhibit a second order phase transition, characterised by a sudden loss of magnetisation (magnitude of the order parameter) at the critical temperature
The four panels in this figure show phase transitions in the heavy top lattice. In figure 7(a), the blue dots represent the $\Gamma_3$-components of the averaged position $\left\langle {\bf{\Gamma }} \right\rangle$ and in 7(b)-7(d), the blue, green and red dots represent the $\Gamma_1$, $\Gamma_2$ and $\Gamma_3$-components of the averaged position respectively. The black line in all four figures correspond to the mean field approximation (106). The top left panel 7(a) shows a phase transition similar to the rigid body phase transition, and the other panels 7(b)-7(d) display a 'triple-humped' phase transition. The corresponding critical temperatures are indicated by the vertical dashed lines. In 7(b)-7(d), we also varied the Casimirs (where $c_1 = {\boldsymbol \Pi} \cdot {\boldsymbol \Gamma}$ and $c_1 = \|{\boldsymbol \Gamma}\|^2$) to show that the three critical temperatures depend on their values, and may even disappear for large enough $c_1/c_2$. Due to the high dimension and non-compactness of the phase space, the mean field approximation requires intensive computation to be precise, explaining the visible errors, especially at low temperatures
 [1] Buddhadev Pal, Pankaj Kumar. A family of multiply warped product semi-riemannian einstein metrics. Journal of Geometric Mechanics, 2020, 12 (4) : 553-562. doi: 10.3934/jgm.2020017 [2] Javier Fernández, Cora Tori, Marcela Zuccalli. Lagrangian reduction of nonholonomic discrete mechanical systems by stages. Journal of Geometric Mechanics, 2020, 12 (4) : 607-639. doi: 10.3934/jgm.2020029 [3] Tian Ma, Shouhong Wang. Topological phase transition III: Solar surface eruptions and sunspots. Discrete & Continuous Dynamical Systems - B, 2020  doi: 10.3934/dcdsb.2020350 [4] Shuang Chen, Jinqiao Duan, Ji Li. Effective reduction of a three-dimensional circadian oscillator model. Discrete & Continuous Dynamical Systems - B, 2020  doi: 10.3934/dcdsb.2020349 [5] Annegret Glitzky, Matthias Liero, Grigor Nika. Dimension reduction of thermistor models for large-area organic light-emitting diodes. Discrete & Continuous Dynamical Systems - S, 2020  doi: 10.3934/dcdss.2020460 [6] Andrew D. Lewis. Erratum for "nonholonomic and constrained variational mechanics". Journal of Geometric Mechanics, 2020, 12 (4) : 671-675. doi: 10.3934/jgm.2020033 [7] Yifan Chen, Thomas Y. Hou. Function approximation via the subsampled Poincaré inequality. Discrete & Continuous Dynamical Systems - A, 2021, 41 (1) : 169-199. doi: 10.3934/dcds.2020296 [8] Alberto Bressan, Wen Shen. A posteriori error estimates for self-similar solutions to the Euler equations. Discrete & Continuous Dynamical Systems - A, 2021, 41 (1) : 113-130. doi: 10.3934/dcds.2020168 [9] Yuxia Guo, Shaolong Peng. A direct method of moving planes for fully nonlinear nonlocal operators and applications. Discrete & Continuous Dynamical Systems - S, 2020  doi: 10.3934/dcdss.2020462 [10] Helmut Abels, Andreas Marquardt. On a linearized Mullins-Sekerka/Stokes system for two-phase flows. Discrete & Continuous Dynamical Systems - S, 2020  doi: 10.3934/dcdss.2020467 [11] Wenjun Liu, Yukun Xiao, Xiaoqing Yue. Classification of finite irreducible conformal modules over Lie conformal algebra $\mathcal{W}(a, b, r)$. Electronic Research Archive, , () : -. doi: 10.3934/era.2020123 [12] Mingjun Zhou, Jingxue Yin. Continuous subsonic-sonic flows in a two-dimensional semi-infinitely long nozzle. Electronic Research Archive, , () : -. doi: 10.3934/era.2020122 [13] Chao Xing, Jiaojiao Pan, Hong Luo. Stability and dynamic transition of a toxin-producing phytoplankton-zooplankton model with additional food. Communications on Pure & Applied Analysis, 2021, 20 (1) : 427-448. doi: 10.3934/cpaa.2020275 [14] Weiwei Liu, Jinliang Wang, Yuming Chen. Threshold dynamics of a delayed nonlocal reaction-diffusion cholera model. Discrete & Continuous Dynamical Systems - B, 2020  doi: 10.3934/dcdsb.2020316 [15] Manil T. Mohan. First order necessary conditions of optimality for the two dimensional tidal dynamics system. Mathematical Control & Related Fields, 2020  doi: 10.3934/mcrf.2020045 [16] Cuicui Li, Lin Zhou, Zhidong Teng, Buyu Wen. The threshold dynamics of a discrete-time echinococcosis transmission model. Discrete & Continuous Dynamical Systems - B, 2020  doi: 10.3934/dcdsb.2020339 [17] Shao-Xia Qiao, Li-Jun Du. Propagation dynamics of nonlocal dispersal equations with inhomogeneous bistable nonlinearity. Electronic Research Archive, , () : -. doi: 10.3934/era.2020116 [18] Ebraheem O. Alzahrani, Muhammad Altaf Khan. Androgen driven evolutionary population dynamics in prostate cancer growth. Discrete & Continuous Dynamical Systems - S, 2020  doi: 10.3934/dcdss.2020426 [19] A. M. Elaiw, N. H. AlShamrani, A. Abdel-Aty, H. Dutta. Stability analysis of a general HIV dynamics model with multi-stages of infected cells and two routes of infection. Discrete & Continuous Dynamical Systems - S, 2020  doi: 10.3934/dcdss.2020441 [20] José Luis López. A quantum approach to Keller-Segel dynamics via a dissipative nonlinear Schrödinger equation. Discrete & Continuous Dynamical Systems - A, 2020  doi: 10.3934/dcds.2020376

2019 Impact Factor: 0.649