Search

The gravitational potentials of realistic galaxy models are in general non-integrable, in the sense that they admit orbits that do not have three independent isolating integrals of motion and are therefore chaotic. However, if chaotic orbits are a small minority in a stellar system, it is expected that they have negligible impact on the main dynamical properties of the system. In this paper we address the question of quantifying the importance of chaotic orbits in a stellar system, focusing, for simplicity, on axisymmetric systems. Chaotic orbits have been found in essentially all (non-St\"ackel) axisymmetric gravitational potentials in which they have been looked for. Based on the analysis of the surfaces of section, we add new examples to those in the literature, finding chaotic orbits, as well as resonantly trapped orbits among regular orbits, in Miyamoto-Nagai, flattened logarithmic and shifted Plummer axisymmetric potentials. We define the fractional contributions in mass of chaotic ($\xi_{\rm c}$) and resonantly trapped ($\xi_{\rm t}$) orbits to a stellar system of given distribution function, which are very useful quantities, for instance in the study of the dispersal of stellar streams of galaxy satellites. As a case study, we measure $\xi_{\rm c}$ and $\xi_{\rm t}$ in two axisymmetric stellar systems obtained by populating flattened logarithmic potentials with the Evans ergodic distribution function, finding $\xi_{\rm c}\sim 10^{-4}-10^{-3}$ and $\xi_{\rm t}\sim 10^{-2}-10^{-1}$., Comment: Accepted for publication in MNRAS. 15 pages, 11 figures

Two counterexamples, addressing questions raised in Adamczak (2019) and Poly and Zheng (2019), are provided. Both counterexamples are related to chaoses. Let F n = Y n + Z n , where the random variables Y n and Z n belong to different chaoses of uniformly bounded degree. It may be that F n ⟶ a . s . 0 , F n ⟶ L 2 + δ 0 and E { sup n | F n | δ } ∞ , for some δ > 0 , and yet Y n fails to converge to 0 a.s.

We study the chaotic-like behavior of cosmological simulations by quantifying how minute perturbations grow over time and manifest as macroscopic differences in galaxy properties. When we run pairs of 'shadow' simulations that are identical except for random minute initial displacements to particle positions (e.g. of order 1e-7pc), the results diverge from each other at the individual galaxy level (while the statistical properties of the ensemble of galaxies are unchanged). After cosmological times, the global properties of pairs of 'shadow' galaxies that are matched between the simulations differ from each other generally at a level of ~2-25%, depending on the considered physical quantity. We perform these experiments using cosmological volumes of (25-50Mpc/h)^3 evolved either purely with dark matter, or with baryons and star-formation but no feedback, or using the full feedback model of the IllustrisTNG project. The runs cover four resolution levels spanning a factor of 512 in mass. We find that without feedback the differences between shadow galaxies generally become smaller as the resolution increases, but with the IllustrisTNG model the results are mostly converging towards a 'floor'. This hints at the role of feedback in setting the chaotic properties of galaxy formation. Importantly, we compare the macroscopic differences between shadow galaxies to the overall scatter in various galaxy scaling relations, and conclude that for the star formation-mass and the Tully-Fisher relations the butterfly effect in our simulations contributes significantly to the overall scatter. We find that our results are robust to whether random numbers are used in the sub-grid models or not. We discuss the implications for galaxy formation theory in general and for cosmological simulations in particular., Comment: Key figures: 10 & 11. Accepted for publication in ApJ. This final version includes a new verification of the conclusions in simulations that completely avoid the usage of random numbers

The emergence of cosmic structure is commonly considered one of the most complex phenomena in Nature. However, this complexity has never been defined nor measured in a quantitative and objective way. In this work we propose a method to measure the information content of cosmic structure and to quantify the complexity that emerges from it, based on Information Theory. The emergence of complex evolutionary patterns is studied with a statistical symbolic analysis of the datastream produced by state-of-the-art cosmological simulations of forming galaxy clusters. This powerful approach allows us to measure how many bits of information are necessary to predict the evolution of energy fields in a statistical way, and it offers a simple way to quantify when, where and how the cosmic gas behaves in complex ways. The most complex behaviors are found in the peripheral regions of galaxy clusters, where supersonic flows drive shocks and large energy fluctuations over a few tens of million years. Describing the evolution of magnetic energy requires at least a twice as large amount of bits than for the other energy fields. When radiative cooling and feedback from galaxy formation are considered, the cosmic gas is overall found to double its degree of complexity. In the future, Cosmic Information Theory can significantly increase our understanding of the emergence of cosmic structure as it represents an innovative framework to design and analyze complex simulations of the Universe in a simple, yet powerful way., Comment: 15 pages, 14 figures. MNRAS accepted, in press

Good quality entropy sources are indispensable in most modern cryptographic protocols. Unfortunately, many currently deployed networked devices do not include them and may be vulnerable to Random Number Generator (RNG) attacks. Since most of these systems allow firmware upgrades and have serial communication facilities, the potential for retrofitting them with secure hardware-based entropy sources exists. To this aim, very low-cost, robust, easy to deploy solutions are required. Here, a retrofittable, sub 10$ entropy source based on chaotic dynamics is illustrated, capable of a 32 kbit/s rate or more and offering multiple serial communication options including USB, I2C, SPI or USART. Operation is based on a loop built around the Analog to Digital Converter (ADC) hosted on a standard microcontroller., Comment: 4 pages, 6 figures. Pre-print from conference proceedings; IEEE 21th International Conference on Electronics, Circuits, and Systems (ICECS 2014), pp. 175-178, Dec. 2014

Data assimilation schemes are confronted with the presence of model errors arising from the imperfect description of atmospheric dynamics. These errors are usually modelled on the basis of simple assumptions such as bias, white noise, and first-order Markov process. In the present work, a formulation of the sequential extended Kalman filter is proposed, based on recent findings on the universal deterministic behaviour of model errors in marked contrast with previous approaches. This new scheme is applied in the context of a spatially distributed system proposed by Lorenz. First, it is found that, for short times, the estimation error is accurately approximated by an evolution law in which the variance of the model error (assumed to be a deterministic process) evolves according to a quadratic law, in agreement with the theory. Moreover, the correlation with the initial condition error appears to play a secondary role in the short-time dynamics of the estimation error covariance. Second, the deterministic description of the model error evolution, incorporated into the classical extended Kalman filter equations, reveals that substantial improvements of the filter accuracy can be gained compared with the classical white-noise assumption. The universal short-time quadratic law for the evolution of the model error covariance matrix seems very promising for modelling estimation error dynamics in sequential data assimilation. Copyright © 2008 Royal Meteorological Society.