Detection of Time Reversibility in Time Series by Ordinal Patterns Analysis J

Total Page:16

File Type:pdf, Size:1020Kb

Detection of Time Reversibility in Time Series by Ordinal Patterns Analysis J Detection of time reversibility in time series by ordinal patterns analysis J. H. Mart´ınez,1, a) J. L. Herrera-Diestra,2 and M. Chavez3 1)INSERM-UM1127, Sorbonne Universit´e,Institut du Cerveau et de la Moelle Epini`ere. France 2)ICTP South American Institute for Fundamental Research, IFT-UNESP. Brazil 3)CNRS UMR7225, H^opitalPiti´eSalp^etri`ere. France (Dated: 13 September 2018) Time irreversibility is a common signature of nonlinear processes, and a fundamental property of non- equilibrium systems driven by non-conservative forces. A time series is said to be reversible if its statistical properties are invariant regardless of the direction of time. Here we propose the Time Reversibility from Ordinal Patterns method (TiROP) to assess time-reversibility from an observed finite time series. TiROP captures the information of scalar observations in time forward, as well as its time-reversed counterpart by means of ordinal patterns. The method compares both underlying information contents by quantifying its (dis)-similarity via Jensen-Shannon divergence. The statistic is contrasted with a population of divergences coming from a set of surrogates to unveil the temporal nature and its involved time scales. We tested TiROP in different synthetic and real, linear and non linear time series, juxtaposed with results from the classical Ramsey's time reversibility test. Our results depict a novel, fast-computation, and fully data-driven method- ology to assess time-reversibility at different time scales with no further assumptions over data. This approach adds new insights about the current non-linear analysis techniques, and also could shed light on determining new physiological biomarkers of high reliability and computational efficiency. PACS numbers: 05.45.Tp, 05.70.Ln, 89.75.Kd, 87.23.-n, 87.19.le, 89.65.Gh, 87.10.Vg Keywords: Time reversibility, Time series, Ordinal patterns analysis, Nonlinearity, Surrogate data Most time series observed from real systems are non-conservative forces (memory)5, therefore, it is ex- inherently nonlinear, thus detecting this property pected to be present in the scalar observation of different is of full interest in natural or social sciences. One biological and physical systems. Indeed, time irreversibil- feature that ensures the nonlinear character of a ity has been reported in ecological and epidemiologi- system is the time irreversibiliity. A time series is cal time series6,7, in tremor time series of patients with said to be reversible if its statistical properties are Parkinson's disease8, in electroencephalographic (EEG) invariant regardless of the direction of time. Here recordings of epileptic patients9{11, or in cardiac inter- we propose the Time Reversibility from Ordinal beat interval time series extracted from patients and Patterns (TiROP) method to assess the tempo- healthy subjects under different cardiac conditions12{15. ral symmetry of linear and nonlinear time series Any time series that is a realisation of a stationary, at different scales. Our approach is based on a linear Gaussian process is time reversible, because of the fast-computing symbolic representation of the ob- symmetry of their covariance functions16{18. Neverthe- served data. Here, TiROP is compared with a less, a non-Gaussian amplitude distribution could be due classical time-reversibility test in a rich variety to a static nonlinear transformation of a stationary lin- of synthetic and real time series from different ear Gaussian process, and by itself is no proof of tempo- systems, including ecology, epidemiology, econ- ral irreversibility. Furthermore, non-Gaussian processes omy and neuroscience. Our results confirm that modeled as outputs of linear systems are reversible19. TiROP has a remarkable performance at unveil- In contrast, the output of a non-linear system excited ing the time scales involved in the temporal irre- by non-Gaussian noises is time irreversible20. Non-linear versibility of a broad range of processes. and non-Gaussian linear models typically have temporal directionality as a property of their higher-order depen- dency18. The study of time reversibility properties of time series might therefore provide meaningful insights I. INTRODUCTION into the underlying nonlinear mechanisms of the observed arXiv:1809.04377v1 [physics.data-an] 12 Sep 2018 data. A time series is said to be reversible if its statisti- Classical time reversibility tests require higher-order 21{23 cal properties are invariant regardless of the direction of moments of the studied signal Xt to be finite . Other time. Time irreversibility is a fundamental property of tests have been devised by directly comparing the dis- 1{3 4,9,24 non-equilibrium systems and dynamics resulting from tribution of vectors fXt;Xt+1; ··· ;Xt+Dg and its time-reversed version fXt+D;Xt+D−1; ··· ;Xtg, or from the projection of dynamics onto a finite number of planes14,25. In the last years, some works have proposed a)Electronic mail: [email protected] statistical tests for irreversibility based on the so-called visibility graphs26, i.e., the mutual visibility relationships 2 between points in a one-dimensional landscape represent- (a) 27{29 ing Xt . These works show that irreversible dynamics results in an asymmetry between the probability distri- ´ butions of graph properties (e.g. links or paths-based Xt Xt characteristics). Recently, this approach has been ex- tended for the study of non-stationary processes30,31. original time series time-reversed signal For real-valued time series, some studies have proposed (b) time-reversibility tests based on different symbolization procedures to characterize the dimensional phase spaces 12,32,33 � of Xt and its time-reversed version . These sym- bolic transformations are generally done by defining a (c) ´ quantization procedure to transform the time series into �(P(�),P(�)) a discrete sequence of unique patterns or symbols13,15,34. P(�) P(´�) Some of these reversibility tests use a priori binomial statistics to assess statistical significance of findings32,33. Nevertheless, such tests assume independence of the ob- served symbols, which is unlikely to occur in real data �1 �2 �3 �4 �5 �6 �1 �2 �3 �4 �5 �6 with temporal correlations. In case of such serial correla- tions, a rigorous theoretical framework cannot be derived FIG. 1. Main steps of the TiROP algorithm for evaluating the time-reversibility of a time series Xt. (a) (Left) Original and Monte Carlo simulations (e.g. parametric or non- time series represented in blue. (Right) The time reversed sig- parametric re-sampling) must be performed to estimate 0 nal Xt represented in orange. (b) Patterns π's extracted from the significance level of time reversibility tests12,13,15. 0 Xt and Xt for D = 3. (c) Probability distributions P (π) and 0 0 In this work we propose a novel procedure, the Time P (π) extracted from Xt and Xt, respectively. The Jensen- Reversibility from Ordinal Patterns method (TiROP), Shannon δ captures the dissimilarity between the information that compares the empirical distributions of the forward content in both distributions and backward statistics of a time series. To estimate the asymmetry between both probability distributions we 40 use the ordinal symbolic representations35,36. In contrast proposals , we considered here the dynamical transfor- 35 with other approaches based on symbolic analysis, the or- mation by Band and Pompe . This method maps a dinal patterns analysis used here is fully data-driven, i.e., time series Xt with t = 1;:::;T to a finite number of pat- the symbolic transformation does not require any a pri- terns that encode the relative amplitudes observed in the ori threshold, or any knowledge about the data sequence. D-dimensional vector Xt = fXt;Xt+τ ;:::;Xt+(D−1)τ g. We complete our time reversibility test with surrogate The elements of the vector Xt are mapped uniquely onto data analysis without making assumptions on the under- the permutation π = (π0; π1; : : : ; πD−1) of (0; 1;:::;D − 37,38 lying generating process . 1) that fulfills Xt+π0τ 6 Xt+π1τ ; 6 ::: 6 Xt+πD−1τ . The proposed framework is validated on synthetic data Each order pattern (permutation) represents thus a sub- simulated with linear, nonlinear, non-Gaussian stochas- set of the whole embedding state space. tic and deterministic processes. The method is also illus- The set of all possible ordinal patterns derived from trated on a collection of different real time series. The a time series is noted as St, whose cardinality is D! at reliability and performances of our method are also com- most. The whole sequence of ordinal patterns extracted pared with those obtained by a classical moment-based from Xt is known as the symbolic representation of the method. The remainder of the paper is organized as time series. The information content of Xt is captured by follows: Section II describes the proposed framework, the probability density P (π) of finding a particular pat- as well as the comparative method used to benchmark tern of order D in St. The higher the order is, the more our solution. Experimental results and evaluation of the information is captured from the time series. To sam- method in synthetic time series are in Section III; while ple the empirical distribution of ordinal patterns densely the evaluation of the test on real data is provided in Sec- enough for a reliable estimation of its probability distri- 41 tion IV. Finally, we conclude the paper with a discussion bution we follow the condition T > (D + 1)! in Section V. The analysis of ordinal representations has some prac- tical advantages36: i) it is computationally efficient, ii) it is fully data-driven with no further assumptions about II. METHODS the data range to find appropriate partitions and, iii) a small D is generally useful in descriptive data analysis35. Capturing information dynamics from time series Furthermore, this symbolisation method is known to be relatively robust against noise, and useful for time series with weak stationarity39,42{46.
Recommended publications
  • Irreversibility of Financial Time Series: a Graph-Theoretical Approach
    Irreversibility of financial time series: a graph-theoretical approach Ryan Flanagan and Lucas Lacasa School of Mathematical Sciences, Queen Mary University of London The relation between time series irreversibility and entropy production has been recently inves- tigated in thermodynamic systems operating away from equilibrium. In this work we explore this concept in the context of financial time series. We make use of visibility algorithms to quantify in graph-theoretical terms time irreversibility of 35 financial indices evolving over the period 1998- 2012. We show that this metric is complementary to standard measures based on volatility and exploit it to both classify periods of financial stress and to rank companies accordingly. We then validate this approach by finding that a projection in principal components space of financial years based on time irreversibility features clusters together periods of financial stress from stable periods. Relations between irreversibility, efficiency and predictability are briefly discussed. I. INTRODUCTION The quantitative analysis of financial time series [1] is a classical field in econometrics that has received in the last decades valuable inputs from statistical physics, nonlinear dynamics and complex systems communities (see [2{4] and references therein for seminal contributions). In particular, the presence of long-range dependence, detected via multifractal measures has been used to quantify the level of development of a given market [5]. This approach was subsequently extended to address the problem of quantifying the degree of market inefficiency [6{10]. Here the adjective efficient refers to a market who is capable of integrating, at any given time, all available information in the price of an asset (the so-called weak-form of the Efficient Market Hypothesis).
    [Show full text]
  • PHYSICA Time-Reversal Symmetry in Dynamical Systems
    PHYSICA ELSEVIER Physica D 112 (1998) 1-39 Time-reversal symmetry in dynamical systems: A survey Jeroen S.W. Lamb a,., John A.G. Roberts b a Mathematics Institute, University of Warwick, Coventry CV4 7AL, UK b School of Mathematics, La Trobe University, Bundoora, Vic. 3083, Australia Abstract In this paper we survey the topic of time-reversal symmetry in dynamical systems. We begin with a brief discussion of the position of time-reversal symmetry in physics. After defining time-reversal symmetry as it applies to dynamical systems, we then introduce a major theme of our survey, namely the relation of time-reversible dynamical sytems to equivariant and Hamiltonian dynamical systems. We follow with a survey of the state of the art on the theory of reversible dynamical systems, including results on symmetric periodic orbits, local bifurcation theory, homoclinic orbits, and renormalization and scaling. Some areas of physics and mathematics in which reversible dynamical systems arise are discussed. In an appendix, we provide an extensive bibliography on the topic of time-reversal symmetry in dynamical systems. 1991 MSC: 58Fxx Keywords: Dynamical systems; Time-reversal symmetry; Reversibility 1. Introduction scope of our survey. Our survey also does not include a discussion on reversible cellular automata. For fur- Time-reversal symmetry is one of the fundamen- ther reading in these areas we recommend the books tal symmetries discussed in natural science. Con- by Brush [6], Sachs [22] and Hawking [14], and the sequently, it arises in many physically motivated survey paper by Toffoli and Margolus [23]. dynamical systems, in particular in classical and Our survey is largely self-contained and accompa- quantum mechanics.
    [Show full text]
  • Time, Symmetry and Structure: a Study in the Foundations of Quantum Theory
    View metadata, citation and similar papers at core.ac.uk brought to you by CORE provided by D-Scholarship@Pitt TIME, SYMMETRY AND STRUCTURE: A STUDY IN THE FOUNDATIONS OF QUANTUM THEORY by Bryan W. Roberts Submitted to the Graduate Faculty of the Department of History & Philosophy of Science in partial fulfillment of the requirements for the degree of Doctor of Philosophy University of Pittsburgh 2012 UNIVERSITY OF PITTSBURGH HISTORY & PHILOSOPHY OF SCIENCE This dissertation was presented by Bryan W. Roberts It was defended on May 30th, 2012 and approved by John Earman, PhD, History & Philosophy of Science John D. Norton, PhD, History & Philosophy of Science Robert Batterman, PhD, Philosophy Laura Ruetsche, PhD, Philosophy (University of Michigan) Giovanni Valente, PhD, Philosophy Dissertation Advisors: John Earman, PhD, History & Philosophy of Science, John D. Norton, PhD, History & Philosophy of Science ii Copyright c by Bryan W. Roberts 2012 iii TIME, SYMMETRY AND STRUCTURE: A STUDY IN THE FOUNDATIONS OF QUANTUM THEORY Bryan W. Roberts, PhD University of Pittsburgh, 2012 This dissertation is about the sense in which the laws of quantum theory distinguish between the past and the future. I begin with an account of what it means for quantum theory to make such a distinction, by providing a novel derivation of the meaning of \time reversal." I then show that if Galilei invariant quantum theory does distinguish a preferred direction in time, then this has consequences for the ontology of the theory. In particular, it requires matter to admit \internal" degrees of freedom, in that the position observable generates a maximal abelian algebra.
    [Show full text]
  • A Note on Time-Reversibility of Multivariate Linear Processes
    A NOTE ON TIME-REVERSIBILITY OF MULTIVARIATE LINEAR PROCESSES KUNG-SIK CHAN Department of Statistics & Actuarial Science, University of Iowa, Iowa City, IA 52242, U.S.A. [email protected] LOP-HING HO Department of Mathematics & Statistics, Wichita State University, Wichita, KS 67260 [email protected] HOWELL TONG Department of Statistics, London School of Economics, WC2A 2AE, U.K. [email protected] March, 2005 SUMMARY We solve an important open problem by deriving some readily verifiable necessary and sufficient conditions for a multivariate non-Gaussian linear process to be time-reversible, under two sets of regularity conditions on the contemporaneous dependence structure of the innovations. One set of regularity conditions concerns the case of independent-component innovations, in which case a multivariate non-Gaussian linear process is time-reversible if and only if the coefficients consist of essentially symmetric columns with column-specific origins of symmetry or symmetric pairs of columns with pair-specific origins of symmetry. On the other hand, for dependent-component innovations plus other regularity conditions, a multivariate non-Gaussian linear process is time- reversible if and only if the coefficients are essentially symmetric about some origin. 1 Some Key Words: Cumulants; Distributional equivalence; Non-Gaussian distribution; Time series; t-distribution; Symmetry. 1 INTRODUCTION & LITERATURE REVIEW Symmetry is all pervasive in art and science. Dawid (1988) has emphasised its role in statis- tics. In his discussion of the paper, Tong (op. cit., p. 30) has identified time-reversibility as an additional symmetry for stationary time series, stationarity itself being an expression of sym- metry.
    [Show full text]
  • Probability Approximations Via the Poisson Clumping Heuristic
    Preface If you place a large number of points randomly in the unit square, what is the distribution of the radius of the largest circle containing no points? Of the smallest circle containing 4 points? Why do Brownian sample paths have local maxima but not points of increase, and how nearly do they have points of increase? Given two long strings of letters drawn i.i.d. from a finite alphabet, how long is the longest consecutive (resp. non-consecutive) substring appearing in both strings? If an imaginary particle performs a simple random walk on the vertices of a high-dimensional cube, how long does it take to visit every vertex? If a particle moves under the influence of a potential field and random perturbations of velocity, how long does it take to escape from a deep potential well? If cars on a freeway move with constant speed (random from car to car), what is the longest stretch of empty road you will see during a long journey? If you take a large i.i.d. sample from a 2-dimensional rotationally-invariant distribution, what is the maximum over all half-spaces of the deviation between the empirical and true distributions? These questions cover a wide cross-section of theoretical and applied probability. The common theme is that they all deal with maxima or min- ima, in some sense. The purpose of this book is to explain a simple idea which enables one to write down, with little effort, approximate solutions to such questions. Let us try to say this idea in one paragraph.
    [Show full text]
  • Time-Reversibility, Causality and Compression-Complexity
    entropy Article Time-Reversibility, Causality and Compression-Complexity Aditi Kathpalia 1,2,* and Nithin Nagaraj 2 1 Department of Complex Systems, Institute of Computer Science of the Czech Academy of Sciences, Czech Academy of Sciences, Pod Vodárenskou vˇeží271/2, 182 07 Prague, Czech Republic 2 Consciousness Studies Programme, National Institute of Advanced Studies (NIAS), Indian Institute of Science Campus, Bengaluru 560012, India; [email protected] * Correspondence: [email protected] Abstract: Detection of the temporal reversibility of a given process is an interesting time series analysis scheme that enables the useful characterisation of processes and offers an insight into the underlying processes generating the time series. Reversibility detection measures have been widely employed in the study of ecological, epidemiological and physiological time series. Further, the time reversal of given data provides a promising tool for analysis of causality measures as well as studying the causal properties of processes. In this work, the recently proposed Compression-Complexity Causality (CCC) measure (by the authors) is shown to be free of the assumption that the "cause precedes the effect", making it a promising tool for causal analysis of reversible processes. CCC is a data-driven interventional measure of causality (second rung on the Ladder of Causation) that is based on Effort-to-Compress (ETC), a well-established robust method to characterize the complexity of time series for analysis and classification. For the detection of the temporal reversibility of processes, we propose a novel measure called the Compressive Potential based Asymmetry Measure. This asymmetry measure compares the probability of the occurrence of patterns at different scales between the forward-time and time-reversed process using ETC.
    [Show full text]
  • On the Analysis of a Random Walk-Jump Chain with Tree-Based Transitions, and Its Applications to Faulty Dichotomous Search∗
    CORE Metadata, citation and similar papers at core.ac.uk Provided by Agder University Research Archive On the Analysis of a Random Walk-Jump Chain with Tree-based Transitions, and its Applications to Faulty Dichotomous Search∗ Anis Yazidiy and B. John Oommenz Abstract Abstract: Random Walks (RWs) have been extensively studied for more than a century [1]. These walks have traditionally been on a line, and the generalizations for two and three dimensions, have been by extending the random steps to the corresponding neighboring positions in one or many of the dimensions. Among the most popular RWs on a line are the various models for birth and death processes, renewal processes and the gambler’s ruin problem. All of these RWs operate “on a discretized line”, and the walk is achieved by performing small steps to the current-state’s neighbor states. Indeed, it is this neighbor-step motion that renders their analyses tractable. When some of the transitions are to non-neighbour states, a formal analysis is, typically, impossible because the difference equations of the steady-state probabilities are not solvable. One endeavor on such an analysis is found in [2]. The problem is far more complex when the transitions of the walk follow an underlying tree-like structure. The analysis of RWs on a tree have received little attention, even though it is an important topic since a tree is a counter-part space representation of a line whenever there is some ordering on the nodes on the line. Nevertheless, RWs on a tree entail moving to non-neighbor states in the space, which makes the analysis involved, and in many cases, impossible.
    [Show full text]
  • Time-Reversal Symmetry Breaking in Quantum Billiards
    Time-Reversal Symmetry Breaking in Quantum Billiards Vom Fachbereich Physik der Technischen Universit¨at Darmstadt zur Erlangung des Grades eines Doktors der Naturwissenschaften (Dr. rer. nat.) genehmigte Dissertation angefertigt von Dipl.-Phys. Florian Sch¨afer aus Dieburg Darmstadt 2009 D 17 Referent: Professor Dr. rer. nat. Dr. h.c. mult. A. Richter Korreferent: Professor Dr. rer. nat. J. Wambach Tag der Einreichung: 1. Dezember 2008 Tag der Prufung:¨ 26. Januar 2009 Abstract The present doctoral thesis describes experimentally measured properties of the resonance spectra of flat microwave billiards with partially broken time- reversal invariance induced by an embedded magnetized ferrite. A vector net- work analyzer determines the complex scattering matrix elements. The data is interpreted in terms of the scattering formalism developed in nuclear physics. At low excitation frequencies the scattering matrix displays isolated reso- nances. At these the effect of the ferrite on isolated resonances (singlets) and pairs of nearly degenerate resonances (doublets) is investigated. The hallmark of time-reversal symmetry breaking is the violation of reciprocity, i.e. of the sym- metry of the scattering matrix. One finds that reciprocity holds in singlets; it is violated in doublets. This is modeled by an effective Hamiltonian of the res- onator. A comparison of the model to the data yields time-reversal symmetry breaking matrix elements in the order of the level spacing. Their dependence on the magnetization of the ferrite is understood in terms of its magnetic properties. At higher excitation frequencies the resonances overlap and the scattering ma- trix elements fluctuate irregularly (Ericson fluctuations). They are analyzed in terms of correlation functions.
    [Show full text]
  • Heat, Temperature and Clausius Inequality in a Model for Active
    www.nature.com/scientificreports OPEN Heat, temperature and Clausius inequality in a model for active Brownian particles Received: 05 December 2016 Umberto Marini Bettolo Marconi1, Andrea Puglisi2 & Claudio Maggi3 Accepted: 17 March 2017 Methods of stochastic thermodynamics and hydrodynamics are applied to a recently introduced Published: 21 April 2017 model of active particles. The model consists of an overdamped particle subject to Gaussian coloured noise. Inspired by stochastic thermodynamics, we derive from the system’s Fokker-Planck equation the average exchanges of heat and work with the active bath and the associated entropy production. We show that a Clausius inequality holds, with the local (non-uniform) temperature of the active bath replacing the uniform temperature usually encountered in equilibrium systems. Furthermore, by restricting the dynamical space to the first velocity moments of the local distribution function we derive a hydrodynamic description where local pressure, kinetic temperature and internal heat fluxes appear and are consistent with the previous thermodynamic analysis. The procedure also shows under which conditions one obtains the unified coloured noise approximation (UCNA): such an approximation neglects the fast relaxation to the active bath and therefore yields detailed balance and zero entropy production. In the last part, by using multiple time-scale analysis, we provide a constructive method (alternative to UCNA) to determine the solution of the Kramers equation and go beyond the detailed balance
    [Show full text]
  • Arxiv:Math-Ph/0007009V2 27 Jun 2001
    MATHEMATICAL FORMALISM FOR ISOTHERMAL LINEAR IRREVERSIBILITY Hong Qian Department of Applied Mathematics, Universityi of Washington, Seattle, WA 98195-2420 E-mail: [email protected] We prove the equivalence among symmetricity, time reversibility, and zero entropy production of the station- ary solutions of linear stochastic differential equations. A sufficient and necessary reversibility condition expressed in terms of the coefficients of the equations is given. The existence of a linear stationary irreversible process is es- tablished. Concerning reversibility, we show that there is a contradistinction between any 1-dimensional station- ary Gaussian process and stationary Gaussian process of dimension n > 1. A concrete criterion for differentiat- ing stationarity and sweeping behavior is also obtained. The mathematical result is a natural generalization of Einstein’s fluctuation-dissipation relation, and provides a arXiv:math-ph/0007009v2 27 Jun 2001 rigorous basis for the isothermal irreversibility in a linear regime which is the basis for applying Onsager’s theory to macromolecules in aqueous solution. Key Words: Brownian motion, Entropy Production Rate, Fluctuation-Dissipation Relation, Reversibility, Stochastic Macromolecular Mechanics, Sweeping 1. INTRODUCTION The stochastic differential equation dx = b(x)+Γξ(t), x Rn (1) dt ∈ with Γ being a nonsingular matrix and ξ(t) being the “derivative” of a n-dimensional Wiener process, has wide applications in science and engi- neering [9]. One standard method for attacking this equation is by finding the fundamental solution to its corresponding Fokker-Planck (Kolmogorov forward) equation ∂P 1 = A P b(x)P , (A = ΓΓT ) (2) ∂t ∇ · 2 ∇ − a parabolic equation for the transition probability P (x, t x0) defined on the entire Rn with the integrability condition | P (x, t) dx < + .
    [Show full text]
  • Lecture-19: Queueing Networks
    Lecture-19: Queueing Networks 1 Migration Processes Corollary 1.1. Consider an M/M/s queue with Poisson(l) arrivals and each server having exponential service time exp(m) service. If l > sm, then the output process in steady state is Poisson(l). Proof. Let X(t) denote the number of customers in the system at time t. Since M/M/s process is a birth and death process, it follows from the previous proposition that fX(t); t ≥ 0g is time reversible. Now going forward in time, the time instants at which X(t) increases by unity are the arrival instants of a Poisson process. Hence, by time reversibility, the time points at which X(t) increases by unity when we go backwards in time also constitutes a Poisson process. But these instants are exactly the departure instants of the forward process. Hence the result. Lemma 1.2. For an ergodic M/M/1 queue in steady state, the following are true. 1. The number of customers present in the system at time t is independent of the sequence of past departures. 2. For FCFS discipline, the waiting time spent in the system (waiting in the queue plus the service time) by a customer is independent of the departure process prior to its departure. Proof. Proofs follow by looking at reversed process. 1. Since the arrival process is Poisson, the future arrivals are independent of the number of customers in the system at the current instant. Looking backwards in time, future arrivals are the past departures. Hence by time reversibility, the number of customers currently in the system are also independent of the past departures.
    [Show full text]
  • Dynamical Systems Theory of Irreversibility
    P. Collet, M. Courbage, S. M´etens, A. Neishtadt, and S. Zaslavsky, Eds., Chaotic Dynamics and Transport in Classical and Quantum Systems, (Kluwer Academic Publishers, Dordrecht, 2005) pp. 107-157. DYNAMICAL SYSTEMS THEORY OF IRREVERSIBILITY Pierre Gaspard Center for Nonlinear Phenomena and Complex Systems, Facult´edes Sciences, Universit´eLibre de Bruxelles, Campus Plaine, Code Postal 231, B-1050 Brussels, Belgium Recent work on the connections between dynamical systems theory and nonequilibrium statisti- cal mechanics is reviewed with emphasis on results which are compatible with Liouville’s theorem. Starting from a general discussion of time-reversal symmetry in the Newtonian scheme, it is shown that the Liouvillian eigenstates associated with the Pollicott-Ruelle resonances spontaneously break the time-reversal symmetry. We explain that such a feature is compatible with the time reversibility of Newton’s equations because of a selection of trajectories which are not time-reversal symmetric. The Pollicott-Ruelle resonances and their associated eigenstates can be constructed not only for decay processes but also for transport processes such as diffusion or viscosity, as well as for reaction- diffusion processes. The Pollicott-Ruelle resonances thus describe the relaxation toward the ther- modynamic equilibrium. The entropy production of these relaxation processes can be calculated and shown to take the value expected from nonequilibrium thermodynamics. In nonequilibrium steady states, an identity is obtained which shows that the entropy production directly character- izes the breaking of time-reversal symmetry by nonequilibrium boundary conditions. The extension to quantum systems is also discussed. I. INTRODUCTION The methods of dynamical systems theory are ubiquitous in nonequilibrium statistical mechanics since the pio- neering work by Sinai and coworkers on the ergodic theory of hard-ball systems [1–4].
    [Show full text]