# 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] Emanuel-Ciprian Cismas. Euler-Poincaré-Arnold equations on semi-direct products II. Discrete & Continuous Dynamical Systems - A, 2016, 36 (11) : 5993-6022. doi: 10.3934/dcds.2016063 [2] Jeffrey K. Lawson, Tanya Schmah, Cristina Stoica. Euler-Poincaré reduction for systems with configuration space isotropy. Journal of Geometric Mechanics, 2011, 3 (2) : 261-275. doi: 10.3934/jgm.2011.3.261 [3] Joachim Escher, Rossen Ivanov, Boris Kolev. Euler equations on a semi-direct product of the diffeomorphisms group by itself. Journal of Geometric Mechanics, 2011, 3 (3) : 313-322. doi: 10.3934/jgm.2011.3.313 [4] Marco Castrillón López, Pablo M. Chacón, Pedro L. García. Lagrange-Poincaré reduction in affine principal bundles. Journal of Geometric Mechanics, 2013, 5 (4) : 399-414. doi: 10.3934/jgm.2013.5.399 [5] Pilar Bayer, Dionís Remón. A reduction point algorithm for cocompact Fuchsian groups and applications. Advances in Mathematics of Communications, 2014, 8 (2) : 223-239. doi: 10.3934/amc.2014.8.223 [6] David Blázquez-Sanz, Juan J. Morales-Ruiz. Lie's reduction method and differential Galois theory in the complex analytic context. Discrete & Continuous Dynamical Systems - A, 2012, 32 (2) : 353-379. doi: 10.3934/dcds.2012.32.353 [7] Luciana A. Alves, Luiz A. B. San Martin. Multiplicative ergodic theorem on flag bundles of semi-simple Lie groups. Discrete & Continuous Dynamical Systems - A, 2013, 33 (4) : 1247-1273. doi: 10.3934/dcds.2013.33.1247 [8] Pierluigi Colli, Antonio Segatti. Uniform attractors for a phase transition model coupling momentum balance and phase dynamics. Discrete & Continuous Dynamical Systems - A, 2008, 22 (4) : 909-932. doi: 10.3934/dcds.2008.22.909 [9] Tian Ma, Shouhong Wang. Cahn-Hilliard equations and phase transition dynamics for binary systems. Discrete & Continuous Dynamical Systems - B, 2009, 11 (3) : 741-784. doi: 10.3934/dcdsb.2009.11.741 [10] Juan Carlos Marrero, D. Martín de Diego, Diana Sosa. Variational constrained mechanics on Lie affgebroids. Discrete & Continuous Dynamical Systems - S, 2010, 3 (1) : 105-128. doi: 10.3934/dcdss.2010.3.105 [11] Katarzyna Grabowska, Marcin Zając. The Tulczyjew triple in mechanics on a Lie group. Journal of Geometric Mechanics, 2016, 8 (4) : 413-435. doi: 10.3934/jgm.2016014 [12] Hiroaki Yoshimura, Jerrold E. Marsden. Dirac cotangent bundle reduction. Journal of Geometric Mechanics, 2009, 1 (1) : 87-158. doi: 10.3934/jgm.2009.1.87 [13] Inês Cruz, M. Esmeralda Sousa-Dias. Reduction of cluster iteration maps. Journal of Geometric Mechanics, 2014, 6 (3) : 297-318. doi: 10.3934/jgm.2014.6.297 [14] Laura Luzzi, Ghaya Rekaya-Ben Othman, Jean-Claude Belfiore. Algebraic reduction for the Golden Code. Advances in Mathematics of Communications, 2012, 6 (1) : 1-26. doi: 10.3934/amc.2012.6.1 [15] Jean-François Babadjian, Francesca Prinari, Elvira Zappale. Dimensional reduction for supremal functionals. Discrete & Continuous Dynamical Systems - A, 2012, 32 (5) : 1503-1535. doi: 10.3934/dcds.2012.32.1503 [16] Martins Bruveris, David C. P. Ellis, Darryl D. Holm, François Gay-Balmaz. Un-reduction. Journal of Geometric Mechanics, 2011, 3 (4) : 363-387. doi: 10.3934/jgm.2011.3.363 [17] Katarzyna Grabowska, Paweƚ Urbański. Geometry of Routh reduction. Journal of Geometric Mechanics, 2019, 11 (1) : 23-44. doi: 10.3934/jgm.2019002 [18] Thomas Espitau, Antoine Joux. Certified lattice reduction. Advances in Mathematics of Communications, 2020, 14 (1) : 137-159. doi: 10.3934/amc.2020011 [19] André Caldas, Mauro Patrão. Entropy of endomorphisms of Lie groups. Discrete & Continuous Dynamical Systems - A, 2013, 33 (4) : 1351-1363. doi: 10.3934/dcds.2013.33.1351 [20] Gerard Thompson. Invariant metrics on Lie groups. Journal of Geometric Mechanics, 2015, 7 (4) : 517-526. doi: 10.3934/jgm.2015.7.517

2018 Impact Factor: 0.525