Chaotic properties of systems with Markov dynamics

We present a general approach for computing the dynamic partition function of a continuous-time Markov process. The Ruelle topological pressure is identified with the large deviation function of a physical observable. We construct for the first time a corresponding finite Kolmogorov-Sinai entropy for these processes. Then, as an example, the latter is computed for a symmetric exclusion process. We further present the first exact calculation of the topological pressure for an N-body stochastic interacting system, namely an infinite-range Ising model endowed with spin-flip dynamics. Expressions for the Kolmogorov-Sinai and the topological entropies follow.Comment: 4 pages, to appear in the Physical Review Letter

Waves and Instabilities in Accretion Disks: MHD Spectroscopic Analysis

A complete analytical and numerical treatment of all magnetohydrodynamic waves and instabilities for radially stratified, magnetized accretion disks is presented. The instabilities are a possible source of anomalous transport. While recovering results on known hydrodynamicand both weak- and strong-field magnetohydrodynamic perturbations, the full magnetohydrodynamic spectra for a realistic accretion disk model demonstrates a much richer variety of instabilities accessible to the plasma than previously realized. We show that both weakly and strongly magnetized accretion disks are prone to strong non-axisymmetric instabilities.The ability to characterize all waves arising in accretion disks holds great promise for magnetohydrodynamic spectroscopic analysis.Comment: FOM-Institute for plasma physics "Rijnhuizen", Nieuwegein, the Netherlands 12 pages, 3 figures, Accepted for publication in ApJ

Observations of Toroidal Coupling for Low-N Alfven Modes in the Tca Tokamak

The antenna structure in the TCA tokamak is phased to excite preferentially Alfven waves with known toroidal and poloidal wave numbers. Surprisingly, the loading spectrum includes both discrete and continuum modes with poloidal wave numbers incompatible with the antenna phasing. These additional modes, which are important for our heating experiments, can be attributed to linear mode coupling induced by the toroidicity of the plasma column, when we take into account ion-cyclotron effects

Chaotic Properties of Dilute Two and Three Dimensional Random Lorentz Gases I: Equilibrium Systems

We compute the Lyapunov spectrum and the Kolmogorov-Sinai entropy for a moving particle placed in a dilute, random array of hard disk or hard sphere scatterers - i.e. the dilute Lorentz gas model. This is carried out in two ways: First we use simple kinetic theory arguments to compute the Lyapunov spectrum for both two and three dimensional systems. In order to provide a method that can easily be generalized to non-uniform systems we then use a method based upon extensions of the Lorentz-Boltzmann (LB) equation to include variables that characterize the chaotic behavior of the system. The extended LB equations depend upon the number of dimensions and on whether one is computing positive or negative Lyapunov exponents. In the latter case the extended LB equation is closely related to an "anti-Lorentz-Boltzmann equation" where the collision operator has the opposite sign from the ordinary LB equation. Finally we compare our results with computer simulations of Dellago and Posch and find very good agreement.Comment: 48 pages, 3 ps fig

Simulating Three-Dimensional Hydrodynamics on a Cellular-Automata Machine

We demonstrate how three-dimensional fluid flow simulations can be carried out on the Cellular Automata Machine 8 (CAM-8), a special-purpose computer for cellular-automata computations. The principal algorithmic innovation is the use of a lattice-gas model with a 16-bit collision operator that is specially adapted to the machine architecture. It is shown how the collision rules can be optimized to obtain a low viscosity of the fluid. Predictions of the viscosity based on a Boltzmann approximation agree well with measurements of the viscosity made on CAM-8. Several test simulations of flows in simple geometries -- channels, pipes, and a cubic array of spheres -- are carried out. Measurements of average flux in these geometries compare well with theoretical predictions.Comment: 19 pages, REVTeX and epsf macros require

Three-dimensional lattice-Boltzmann simulations of critical spinodal decomposition in binary immiscible fluids

We use a modified Shan-Chen, noiseless lattice-BGK model for binary immiscible, incompressible, athermal fluids in three dimensions to simulate the coarsening of domains following a deep quench below the spinodal point from a symmetric and homogeneous mixture into a two-phase configuration. We find the average domain size growing with time as $t^\gamma$, where $\gamma$ increases in the range $0.545 < \gamma < 0.717$, consistent with a crossover between diffusive $t^{1/3}$ and hydrodynamic viscous, $t^{1.0}$, behaviour. We find good collapse onto a single scaling function, yet the domain growth exponents differ from others' works' for similar values of the unique characteristic length and time that can be constructed out of the fluid's parameters. This rebuts claims of universality for the dynamical scaling hypothesis. At early times, we also find a crossover from $q^2$ to $q^4$ in the scaled structure function, which disappears when the dynamical scaling reasonably improves at later times. This excludes noise as the cause for a $q^2$ behaviour, as proposed by others. We also observe exponential temporal growth of the structure function during the initial stages of the dynamics and for wavenumbers less than a threshold value.Comment: 45 pages, 18 figures. Accepted for publication in Physical Review

Persistence exponents in a 3D symmetric binary fluid mixture

The persistence exponent, theta, is defined by N_F sim t^theta, where t is the time since the start of the coarsening process and the "no-flip fraction", N_F, is the number of points that have not seen a change of "color" since t=0. Here we investigate numerically the persistence exponent for a binary fluid system where the coarsening is dominated by hydrodynamic transport. We find that N_F follows a power law decay (as opposed to exponential) with the value of theta somewhat dependent on the domain growth rate (L sim t^alpha, where L is the average domain size), in the range theta=1.23 +-0.1 (alpha = 2/3) to theta=1.37 +-0.2 (alpha=1). These alpha values correspond to the inertial and viscous hydrodynamic regimes respectively.Comment: 9 pages RevTex, 9 figures included as eps files on last 3 pages, submitted to Phys Rev

Sensitivity and performance of the Advanced LIGO detectors in the third observing run

On April 1st, 2019, the Advanced Laser Interferometer Gravitational-Wave Observatory (aLIGO), joined by the Advanced Virgo detector, began the third observing run, a year-long dedicated search for gravitational radiation. The LIGO detectors have achieved a higher duty cycle and greater sensitivity to gravitational waves than ever before, with LIGO Hanford achieving angle-averaged sensitivity to binary neutron star coalescences to a distance of 111 Mpc, and LIGO Livingston to 134 Mpc with duty factors of 74.6% and 77.0% respectively. The improvement in sensitivity and stability is a result of several upgrades to the detectors, including doubled intracavity power, the addition of an in-vacuum optical parametric oscillator for squeezed-light injection, replacement of core optics and end reaction masses, and installation of acoustic mode dampers. This paper explores the purposes behind these upgrades, and explains to the best of our knowledge the noise currently limiting the sensitivity of each detector. © 2020 authors. Published by the American Physical Society

Chromatin accessibility dynamics across C. elegans development and ageing

An essential step for understanding the transcriptional circuits that control development and physiology is the global identification and characterization of regulatory elements. Here we present the first map of regulatory elements across the development and ageing of an animal, identifying 42,245 elements accessible in at least one C. elegans stage. Based on nuclear transcription profiles, we define 15,714 protein-coding promoters and 19,231 putative enhancers, and find that both types of element can drive orientation-independent transcription. Additionally, more than 1000 promoters produce transcripts antisense to protein coding genes, suggesting involvement in a widespread regulatory mechanism. We find that the accessibility of most elements changes during development and/or ageing and that patterns of accessibility change are linked to specific developmental or physiological processes. The map and characterization of regulatory elements across C. elegans life provides a platform for understanding how transcription controls development and ageing