We study a numerical instability of direct simulations with truncated equation chains for the circular cumulant representation and two approaches to its suppression. The approaches are tested for a chimera-bearing hierarchical population of coupled oscillators. The stabilization techniques can be efficiently applied without significant effect on the natural system dynamics within a finite vicinity of the Ott-Antonsen manifold for direct numerical simulations with up to 20 cumulants; with increasing deviation from the Ott-Antonsen manifold the stabilization becomes more problematic.
We demonstrate the application of the circular cumulant approach for thermodynamically large populations of phase elements, where the Ott-Antonsen properties are violated by a multiplicative intrinsic noise. The infinite cumulant equation chain is de
rived for the case of a sinusoidal sensitivity of the phase to noise. For inhomogeneous populations, a Lorentzian distribution of natural frequencies is adopted. Two-cumulant model reductions, which serve as a generalization of the Ott-Antonsen ansatz, are reported. The accuracy of these model reductions and the macroscopic collective dynamics of the system are explored for the case of a Kuramototype global coupling. The Ott-Antonsen ansatz and the Gaussian approximation are found to be not uniformly accurate for non-high frequencies.
The cumulant representation is common in classical statistical physics for variables on the real line and the issue of closures of cumulant expansions is well elaborated. The case of phase variables significantly differs from the case of linear ones;
the relevant order parameters are the Kuramoto-Daido ones but not the conventional moments. One can formally introduce `circular cumulants for Kuramoto-Daido order parameters, similar to the conventional cumulants for moments. The circular cumulant expansions allow to advance beyond the Ott-Antonsen theory and consider populations of real oscillators. First, we show that truncation of circular cumulant expansions, except for the Ott-Antonsen case, is forbidden. Second, we compare this situation to the case of the Gaussian distribution of a linear variable, where the second cumulant is nonzero and all the higher cumulants are zero, and elucidate why keeping up to the second cumulant is admissible for a linear variable, but forbidden for circular cumulants. Third, we discuss the implication of this truncation issue to populations of quadratic integrate-and-fire neurons [E. Montbrio, D. Pazo, A. Roxin, Phys. Rev. X, vol. 5, 021028 (2015)], where within the framework of macroscopic description, the firing rate diverges for any finite truncation of the cumulant series, and discuss how one should handle these situations. Fourth, we consider the cumulant-based low-dimensional reductions for macroscopic population dynamics in the context of this truncation issue. These reductions are applicable, where the cumulant series exponentially decay with the cumulant order, i.e., they form a geometric progression hierarchy. Fifth, we demonstrate the formation of this hierarchy for generic distributions on the circle and experimental data for coupled biological and electrochemical oscillators.
We study synchronization dynamics of a population of pulse-coupled oscillators. In particular, we focus our attention in the interplay between networks topological disorder and its synchronization features. Firstly, we analyze synchronization time $T
$ in random networks, and find a scaling law which relates $T$ to networks connectivity. Then, we carry on comparing synchronization time for several other topological configurations, characterized by a different degree of randomness. The analysis shows that regular lattices perform better than any other disordered network. The fact can be understood by considering the variability in the number of links between two adjacent neighbors. This phenomenon is equivalent to have a non-random topology with a distribution of interactions and it can be removed by an adequate local normalization of the couplings.
A finite element program is presented to simulate the process of packing and coiling elastic wires in two- and three-dimensional confining cavities. The wire is represented by third order beam elements and embedded into a corotational formulation to
capture the geometric nonlinearity resulting from large rotations and deformations. The hyperbolic equations of motion are integrated in time using two different integration methods from the Newmark family: an implicit iterative Newton-Raphson line search solver, and an explicit predictor-corrector scheme, both with adaptive time stepping. These two approaches reveal fundamentally different suitability for the problem of strongly self-interacting bodies found in densely packed cavities. Generalizing the spherical confinement symmetry investigated in recent studies, the packing of a wire in hard ellipsoidal cavities is simulated in the frictionless elastic limit. Evidence is given that packings in oblate spheroids and scalene ellipsoids are energetically preferred to spheres.
Supercritical fluids play a significant role in elucidating fundamental aspects of liquid matter under extreme conditions. They have been extensively studied at pressures and temperatures relevant to various industrial applications. However, much les
s is known about the structural behaviour of supercritical fluids and no structural crossovers have been observed in static compression experiments in any temperature and pressure ranges beyond the critical point. The structure of supercritical state is currently perceived to be uniform everywhere on the pressure-temperature phase diagram, and to change only in a monotonic way even moving around the critical point, not only along isotherms or isobars. Conversely, we observe structural crossovers for the first time in a deeply supercritical sample through diffraction measurements in a diamond anvil cell and discover a new thermodynamic boundary on the pressure-temperature diagram. We explain the existence of these crossovers in the framework of the phonon theory of liquids using molecular dynamics simulations. The obtained results are of prime importance since it implies a global reconsideration of the mere essence of the supercritical phase. Furthermore, this discovery may pave the way to new unexpected applications and to the exploration of exotic behaviour of confined fluids relevant to geo- and planetary sciences.
Irina V. Tyulkina
,Denis S. Goldobin
,Arkady Pikovsky
.
(2019)
.
"Stabilization of direct numerical simulation for finite truncations of circular cumulant expansions"
.
Denis Goldobin
هل ترغب بارسال اشعارات عن اخر التحديثات في شمرا-اكاديميا