arxiv: v1 [q-bio.nc] 5 Dec 2015

Size: px
Start display at page:

Download "arxiv: v1 [q-bio.nc] 5 Dec 2015"

Transcription

1 Modelling axon growing using CTRW arxiv: v1 [q-bio.nc] 5 Dec 2015 Xavier Descombes, Inria Sophia Antipolis, France Elena Zhizhina, Sergey Komech Institute for Information Transmission Problems, Russian Academy of Sciences, Russia 1 Description and analysis of the model The main goal of this study is to propose a mathematical model describing paths of the axon growth cones and differences in the behavior of normal and mutant axons. We introduce a probabilistic model for axon growing, such that each family of axons is described as an ensemble of trajectories of a continuous time random walk (CTRW) model under different parameters in the case of normal and mutant axons. We describe different regimes in the model and conclude how the behavior of axons depends on the parameters of the model. Biological observations of the axonal growth process say us that the guiding development of axons to their targets is operated by chemical signals from the cellular environment. To simulate this control mechanism we propose the CTRW model, where a random waiting time reflects a reaction time of the growth cones on the neighboring chemical environment. Continuous-time random walks (CTRW) are a natural generalization of the usual random walk. The mathematical analysis can be dated back to the pioneering work of Montroll and Weiss in the sixties [1]. At present they are extensively used in applications to physics, chemistry, and other sciences, see e.g. [2], [3]. We consider three families of (static) trajectories associated with axon growth cones. One of the family is for normal axons, the other two are for two different axon mutations. We don t have dynamical data, and any trajectory presents positions of an axon tip during the observation time. The observation time is the same for all axons, but we don t have information about the stopping time for development of each axon. On figures 1-3 we present three sets of trajectories of axon growth cones. In fact, axons on real data images have different starting positions, but we shifted all of them to have the same starting point. The work is partially supported by RAS-CNRS grant EDC26091 and RFBR grants and

2 Figure 1: Axonal trajectories On fig. 4 we present the values of the lengths for the normal and mutant axons. One can see that the average length for the mutant axons is less then for the normal ones. The second observation concerns the value of the average deviation of trajectories from the origin. In the case of mutant 2 trajectories it is much less then in the normal case. To explain these facts we propose a model of axon growing based on a CTRW model. A CTRW model is a symmetrical random walk subordinated to a renewal process. It is defined by two probability distributions: a distribution p(x y), x, y Z d for a spacial random walk and a distribution f(t), t 0 for a waiting time. Then CTRW can be described as follows: a particle wait at position x a random time distributed by f(t), then it jumps from x to y according to the distribution p(x y), and so on. Then a more long waiting time implies reducing of trajectories lengths. Since we don t have any dynamical data for axon growing, we couldn t control the absolute time and should introduce in our model an artificial time. We study here three cases. The first case (A uniform partition). Assuming an equal absolute time T for growing of all axons and the existence of the average speed for any trajectory of the CTRW X(t) Z d, t [0, T ] (depending on the length of the trajectory) we can consider uniform partitions of trajectories proportional to the length L of axons: X(0) = 0, X(T ) = x end, X( k 10 T ) = x k, k = 0, 1,..., 10, where x end is the position of the end of the trajectory of the length L, and x k, k = 0, 1,..., 10 are points on the trajectory corresponding to the length k 10L, k = 0, 1,..., 10. 2

3 Figure 2: Length distributions The second case (A uniform partition after elongation). In this case we assume different time for axon growing but the same evolving time T for all axons. For long axons we consider again an equal absolute growing/evolving time as above, assuming different average speeds of the growth dependent on the length L. For short axons we say that the growing process stops at some moment of time, and after that the short axons just keep their final position in evolution process. Let L n be an average length of trajectories of normal axons, then we put L n = T as the absolute evolving time and consider the following elongation in time for all short axons with a length L less than L n = T : we take the last part of the trajectory (from time L to time T ) coinciding with the end point of axons. Thus we obtain that all short trajectories evolve during the same time T = L n with a constant speed 1, and we assume that it is the same absolute growing/evolving time T = L n for long trajectories. Finally, we construct the uniform partitions depending on the length of axons L, which is the same as in the first case. The third case (A non uniform partition). Assuming again the same absolute time for growing of all axons and using observations from biological experiments that a speed of axon growing is decreasing in time (especially for mutant axons), we consider a non uniform partitions of trajectories corresponding to the following lengths l(k): l(k) = αl tk 1 e α e αt dt, t k = k, k = 0, 1, with α = 2. Since we have a permanent drift along the X-axis, then we will consider the deviation only along the Y-axis, considering X direction as a direction of an artificial time. We will use the following result for the mean squared deviation Y 2 (t), where Y (t) is the position of the Y-coordinate of the CTRW when time t is large enough. 3

4 Proposition. Let p(u) = p( u) is a distribution of a symmetrical random walk on Z 1, and we consider two cases: 1) diffusive behaviour: τ = tf(t)dt <, 0 2) anomalous diffusion: τ = and f(t) 1 t as t with α (0, 1). α+1 Then for large enough t we have for the mean squared displacement and Y 2 (t) 2a τ t in the first case (1) Y 2 (t) t α in the second case with α (0, 1). (2) Here a is the dispersion of the symmetrical random walk. The proof of the Proposition is based on the Fourier-Laplace transform of the characteristic function of the random walk and on the formula Ŷ 2 (s) = where ˆθ(λ, s) = 0 e st e iλy (t) dt. Y 2 (t) e st dt = 2 ˆθ(λ, s) λ 2 λ=0, The main idea of our analysis is to compare statistical characteristics of the trajectories from different families (for normal and mutant axons separately) with the mean squared displacements for CTRW given by (1)-(2). Using the empirical law for Y 2 (t) as a function of t we can conclude which parameters of the CTRW model imply the similar law for Y 2 (t). On Fig we present our calculations of Y 2 (t) and (Y (t) < Y (t) >) 2 using the statistical data and construct corresponding interpolation curves in each of three cases under our consideration. Here is the empirical average. 2 Conclusions We propose a continuous time random walk model as a model of axon growing. We can observe two different regimes during evolution of axon growth cones. The second regime, which is started approximately after 0.6 L, is related with, see formula (1). That can be the result of decreasing averaged waiting time τ or increasing dispersion of the spacial random walk a. The distributions of the waiting time for normal and mutant axons are different, and by the formula (1) in the proposition the average waiting time for mutant axons is longer (greater) than for normal axons. Using our results on the shape of Y 2 (t) for mutant axons we can formulate a hypothesis that for normal axons the time between renewals has a finite mean, whereas for mutant axons the time between renewals is greater in average or even can have an infinite mean. In the latter case, the scaling limit for the CTRW is an operator increasing coefficient 2a τ 4

5 Levy motion subordinated to the hitting time process of a classical stable subordinator, see e.g. [4]. In this case, for large enough t we have Y 2 (t) t α with 0 < α < 1. That means that the hitting time for this CTRW to reach distant compact sets is much more greater for mutant axons than for normal ones. References [1] E.W. Montroll and G.H.: Weiss Random walks on lattices. II. J. Math. Phys. 6, (1965) MR [2] Kolomeisky, A.B.: Continuous-time random walks at all times. J. Chem. Phys. 131, (2009). doi: / [3] Rudnick, J., Gaspari, G.: Elements of the Random Walk. Cambridge University Press (2004) [4] M. Meerschaert, H.-P. Scheffler Limit theorems for continuous-time random walks with infinite mean waiting time, J. Appl. Prob. 41, 623 (2004). 5

6 Figure 3: The first case: second moments as a function of time(% of length) Figure 4: The first case: variances as a function of time(% of length) 6

7 Figure 5: The second case (elongation of short axons): second moments as a function of time(% of length) Figure 6: The second case (elongation of short axons): variances as a function of time(% of length) 7

8 Figure 7: The third case : second moments as a function of time(% of length with weights) Figure 8: The third case : variances as a function of time(% of length with weights) 8

Anomalous diffusion in cells: experimental data and modelling

Anomalous diffusion in cells: experimental data and modelling Anomalous diffusion in cells: experimental data and modelling Hugues BERRY INRIA Lyon, France http://www.inrialpes.fr/berry/ CIMPA School Mathematical models in biology and medicine H. BERRY - Anomalous

More information

Lecture 25: Large Steps and Long Waiting Times

Lecture 25: Large Steps and Long Waiting Times Lecture 25: Large Steps and Long Waiting Times Scribe: Geraint Jones (and Martin Z. Bazant) Department of Economics, MIT Proofreader: Sahand Jamal Rahi Department of Physics, MIT scribed: May 10, 2005,

More information

Functional Limit theorems for the quadratic variation of a continuous time random walk and for certain stochastic integrals

Functional Limit theorems for the quadratic variation of a continuous time random walk and for certain stochastic integrals Functional Limit theorems for the quadratic variation of a continuous time random walk and for certain stochastic integrals Noèlia Viles Cuadros BCAM- Basque Center of Applied Mathematics with Prof. Enrico

More information

Weak Ergodicity Breaking WCHAOS 2011

Weak Ergodicity Breaking WCHAOS 2011 Weak Ergodicity Breaking Eli Barkai Bar-Ilan University Bel, Burov, Korabel, Margolin, Rebenshtok WCHAOS 211 Outline Single molecule experiments exhibit weak ergodicity breaking. Blinking quantum dots,

More information

Continuous and Discrete random process

Continuous and Discrete random process Continuous and Discrete random and Discrete stochastic es. Continuous stochastic taking values in R. Many real data falls into the continuous category: Meteorological data, molecular motion, traffic data...

More information

Weak chaos, infinite ergodic theory, and anomalous diffusion

Weak chaos, infinite ergodic theory, and anomalous diffusion Weak chaos, infinite ergodic theory, and anomalous diffusion Rainer Klages Queen Mary University of London, School of Mathematical Sciences Marseille, CCT11, 24 May 2011 Weak chaos, infinite ergodic theory,

More information

Weak Ergodicity Breaking

Weak Ergodicity Breaking Weak Ergodicity Breaking Eli Barkai Bar-Ilan University 215 Ergodicity Ergodicity: time averages = ensemble averages. x = lim t t x(τ)dτ t. x = xp eq (x)dx. Ergodicity out of equilibrium δ 2 (, t) = t

More information

Volatility and Returns in Korean Futures Exchange Markets

Volatility and Returns in Korean Futures Exchange Markets Volatility and Returns in Korean Futures Exchange Markets Kyungsik Kim*,, Seong-Min Yoon and Jum Soo Choi Department of Physics, Pukyong National University, Pusan 608-737, Korea Division of Economics,

More information

Random walk through Anomalous processes: some applications of Fractional calculus

Random walk through Anomalous processes: some applications of Fractional calculus Random walk through Anomalous processes: some applications of Fractional calculus Nickolay Korabel Fractional Calculus Day @ UCMerced 12 June 2013 Fractional Calculus Mathematics Engineering Physics Levy

More information

Weak Ergodicity Breaking. Manchester 2016

Weak Ergodicity Breaking. Manchester 2016 Weak Ergodicity Breaking Eli Barkai Bar-Ilan University Burov, Froemberg, Garini, Metzler PCCP 16 (44), 24128 (2014) Akimoto, Saito Manchester 2016 Outline Experiments: anomalous diffusion of single molecules

More information

Covariance structure of continuous time random walk limit processes

Covariance structure of continuous time random walk limit processes Covariance structure of continuous time random walk limit processes Alla Sikorskii Department of Statistics and Probability Michigan State University Fractional Calculus, Probability and Non-local Operators:

More information

Lecture notes for /12.586, Modeling Environmental Complexity. D. H. Rothman, MIT September 24, Anomalous diffusion

Lecture notes for /12.586, Modeling Environmental Complexity. D. H. Rothman, MIT September 24, Anomalous diffusion Lecture notes for 12.086/12.586, Modeling Environmental Complexity D. H. Rothman, MIT September 24, 2014 Contents 1 Anomalous diffusion 1 1.1 Beyond the central limit theorem................ 2 1.2 Large

More information

The Portevin-Le Chatelier effect in the Continuous Time Random Walk. framework

The Portevin-Le Chatelier effect in the Continuous Time Random Walk. framework The Portevin-Le Chatelier effect in the Continuous Time Random Walk framework A. Sarkar * and P. Barat Variable Energy Cyclotron Centre 1/AF Bidhan Nagar, Kolkata 764, India PACS numbers: 62.2.Fe, 5.4.Fb,

More information

SUBDYNAMICS OF FINANCIAL DATA FROM FRACTIONAL FOKKER PLANCK EQUATION

SUBDYNAMICS OF FINANCIAL DATA FROM FRACTIONAL FOKKER PLANCK EQUATION Vol. 4 (29) ACTA PHYSICA POLONICA B No 5 SUBDYNAMICS OF FINANCIAL DATA FROM FRACTIONAL FOKKER PLANCK EQUATION Joanna Janczura, Agnieszka Wyłomańska Institute of Mathematics and Computer Science Hugo Steinhaus

More information

Cover times of random searches

Cover times of random searches Cover times of random searches by Chupeau et al. NATURE PHYSICS www.nature.com/naturephysics 1 I. DEFINITIONS AND DERIVATION OF THE COVER TIME DISTRIBUTION We consider a random walker moving on a network

More information

arxiv: v1 [cond-mat.stat-mech] 23 Feb 2017

arxiv: v1 [cond-mat.stat-mech] 23 Feb 2017 First passage properties in a crowded environment Vincent Tejedor 1 1 Cour des comptes, 13, rue Cambon, 75 001 Paris, France February 27, 2017 arxiv:1702.07394v1 [cond-mat.stat-mech] 23 Feb 2017 Introduction

More information

MORTAL CREEPERS SEARCHING FOR A TARGET

MORTAL CREEPERS SEARCHING FOR A TARGET MORTAL CREEPERS SEARCHING FOR A TARGET D. Campos, E. Abad, V. Méndez, S. B. Yuste, K. Lindenberg Barcelona, Extremadura, San Diego We present a simple paradigm for detection of an immobile target by a

More information

arxiv: v1 [cond-mat.stat-mech] 23 Apr 2010

arxiv: v1 [cond-mat.stat-mech] 23 Apr 2010 Fractional Fokker-Planck Equations for Subdiffusion with Space-and-Time-Dependent Forces. B.. Henry Department of Applied Mathematics, School of Mathematics and Statistics, University of New South Wales,

More information

arxiv:cond-mat/ v2 [cond-mat.stat-mech] 23 Feb 1998

arxiv:cond-mat/ v2 [cond-mat.stat-mech] 23 Feb 1998 Multiplicative processes and power laws arxiv:cond-mat/978231v2 [cond-mat.stat-mech] 23 Feb 1998 Didier Sornette 1,2 1 Laboratoire de Physique de la Matière Condensée, CNRS UMR6632 Université des Sciences,

More information

Non-homogeneous random walks on a semi-infinite strip

Non-homogeneous random walks on a semi-infinite strip Non-homogeneous random walks on a semi-infinite strip Chak Hei Lo Joint work with Andrew R. Wade World Congress in Probability and Statistics 11th July, 2016 Outline Motivation: Lamperti s problem Our

More information

CONTINUOUS TIME RANDOM WALK WITH CORRELATED WAITING TIMES

CONTINUOUS TIME RANDOM WALK WITH CORRELATED WAITING TIMES CONTINUOUS TIME RANDOM WALK WITH CORRELATED WAITING TIMES Aleksei V. Chechkin, 1,2 Michael Hofmann 3, and Igor M. Sokolov 3 1 School of Chemistry, Tel Aviv University, Ramat Aviv, 69978 Tel Aviv, Israel

More information

Anomalous diffusion in biology: fractional Brownian motion, Lévy flights

Anomalous diffusion in biology: fractional Brownian motion, Lévy flights Anomalous diffusion in biology: fractional Brownian motion, Lévy flights Jan Korbel Faculty of Nuclear Sciences and Physical Engineering, CTU, Prague Minisymposium on fundamental aspects behind cell systems

More information

Semi-Markov approach to continuous time random walk limit processes

Semi-Markov approach to continuous time random walk limit processes Semi-Markov approach to continuous time random walk limit processes Joint Mathematics Meetings San Antonio, TX January 13, 2015 Mark M. Meerschaert Michigan State University East Lansing, Michigan, USA

More information

Revisiting the microfield formulation using Continuous. Time Random. Walk

Revisiting the microfield formulation using Continuous. Time Random. Walk Revisiting the microfield formulation using Continuous Time Random Walk Theory D. Boland, A. Bouzaher, H. Capes, F. Catoire, M. Christova a, L. Godbert-Mouret, M. Koubiti, S. Mekkaoui, Y. Marandet, J.

More information

The continuous time random walk (CTRW) was introduced by Montroll and Weiss 1.

The continuous time random walk (CTRW) was introduced by Montroll and Weiss 1. 1 I. CONTINUOUS TIME RANDOM WALK The continuou time random walk (CTRW) wa introduced by Montroll and Wei 1. Unlike dicrete time random walk treated o far, in the CTRW the number of jump n made by the walker

More information

Random Walk Model with Waiting Times Depending on the Preceding Jump Length

Random Walk Model with Waiting Times Depending on the Preceding Jump Length Journal of Statistical Physics, Vol. 123, No. 4, May 2006 ( C 2006 ) DOI: 10.1007/s10955-006-9104-z Random Walk Model with Waiting Times Depending on the Preceding Jump Length Vasily Yu. Zaburdaev 1 Received

More information

Lecture 1: Random walk

Lecture 1: Random walk Lecture : Random walk Paul C Bressloff (Spring 209). D random walk q p r- r r+ Figure 2: A random walk on a D lattice. Consider a particle that hops at discrete times between neighboring sites on a one

More information

Continuous time random walks in finite domains and general boundary conditions: some formal considerations

Continuous time random walks in finite domains and general boundary conditions: some formal considerations IOP PUBISHING JOURNA OF PHYSICS A: MATHEMATICA AND THEORETICA J. Phys. A: Math. Theor. 41 (008) 15004 (11pp) doi:10.1088/1751-8113/41/1/15004 Continuous time random walks in finite domains and general

More information

Tailing of the breakthrough curve in aquifer contaminant transport: equivalent longitudinal macrodispersivity and occurrence of anomalous transport

Tailing of the breakthrough curve in aquifer contaminant transport: equivalent longitudinal macrodispersivity and occurrence of anomalous transport 5 GQ7: Securing Groundwater Quality in Urban and Industrial Environments (Proc. 6th International Groundwater Quality Conference held in Fremantle, Western Australia, 2 7 December 27). IAHS Publ. no. XXX,

More information

Chapter 1. An Introduction to Fractional Diffusion

Chapter 1. An Introduction to Fractional Diffusion Chapter An Introduction to Fractional Diffusion B.I. Henry, T.A.M. Langlands and P. Straka School of Mathematics and Statistics, University of New South Wales, Sydney, NSW 252 Australia b.henry@unsw.edu.au

More information

arxiv: v2 [cond-mat.stat-mech] 12 Apr 2011

arxiv: v2 [cond-mat.stat-mech] 12 Apr 2011 epl draft arxiv:8.762v2 cond-mat.stat-mech] 2 Apr 2 Raman Research Institute, Bangalore 568, India PACS 2.5.-r Probability theory, stochastic processes, and statistics PACS 5.4.-a Fluctuation phenomena,

More information

Russia. St. Petersburg, , Russia

Russia. St. Petersburg, , Russia PACS 81.07.Bc, 61.72.Cc An asymptotic formula for displacement field in triangular lattice with vacancy V.A. Tsaplin 1, V.A. Kuzkin 1,2 vtsaplin@yandex.ru, kuzkinva@gmail.com 1 Institute for Problems in

More information

Remember... Average rate of change slope of a secant (between two points)

Remember... Average rate of change slope of a secant (between two points) 3.7 Rates of Change in the Natural and Social Sciences Remember... Average rate of change slope of a secant (between two points) Instantaneous rate of change slope of a tangent derivative We will assume

More information

Remember... Average rate of change slope of a secant (between two points)

Remember... Average rate of change slope of a secant (between two points) 3.7 Rates of Change in the Natural and Social Sciences Remember... Average rate of change slope of a secant (between two points) Instantaneous rate of change slope of a tangent derivative We will assume

More information

Random Times and Their Properties

Random Times and Their Properties Chapter 6 Random Times and Their Properties Section 6.1 recalls the definition of a filtration (a growing collection of σ-fields) and of stopping times (basically, measurable random times). Section 6.2

More information

arxiv:cond-mat/ v1 2 Sep 2003

arxiv:cond-mat/ v1 2 Sep 2003 Bulk Mediated Surface Diffusion: The Infinite System Case Jorge A. Revelli a, Carlos. E. Budde b, Domingo Prato b and Horacio S. Wio a,c (a) Grupo de Física Estadística, Centro Atómico Bariloche and Instituto

More information

FRACTIONAL BROWNIAN MOTION AND DYNAMIC APPROACH TO COMPLEXITY. Rasit Cakir. Dissertation Prepared for the Degree of DOCTOR OF PHILOSOPHY

FRACTIONAL BROWNIAN MOTION AND DYNAMIC APPROACH TO COMPLEXITY. Rasit Cakir. Dissertation Prepared for the Degree of DOCTOR OF PHILOSOPHY FRACTIONAL BROWNIAN MOTION AND DYNAMIC APPROACH TO COMPLEXITY Rasit Cakir Dissertation Prepared for the Degree of DOCTOR OF PHILOSOPHY UNIVERSITY OF NORTH TEXAS August 27 APPROVED: Paolo Grigolini, Major

More information

Noisy continuous time random walks

Noisy continuous time random walks Noisy continuous time random walks Jae-Hyung Jeon, Eli Barkai, and Ralf Metzler Citation: J. Chem. Phys. 139, 121916 (2013); doi: 10.1063/1.4816635 View online: http://dx.doi.org/10.1063/1.4816635 View

More information

Theory of fractional Lévy diffusion of cold atoms in optical lattices

Theory of fractional Lévy diffusion of cold atoms in optical lattices Theory of fractional Lévy diffusion of cold atoms in optical lattices, Erez Aghion, David Kessler Bar-Ilan Univ. PRL, 108 230602 (2012) PRX, 4 011022 (2014) Fractional Calculus, Leibniz (1695) L Hospital:

More information

An Introduction to Probability Theory and Its Applications

An Introduction to Probability Theory and Its Applications An Introduction to Probability Theory and Its Applications WILLIAM FELLER (1906-1970) Eugene Higgins Professor of Mathematics Princeton University VOLUME II SECOND EDITION JOHN WILEY & SONS Contents I

More information

Introduction: the plague. The Scaling Laws Of Human Travel. Introduction: SARS (severe acute respiratory syndrom)

Introduction: the plague. The Scaling Laws Of Human Travel. Introduction: SARS (severe acute respiratory syndrom) The Scaling Laws Of Human Travel Verena Zuber Department of Statistics University of Munich 7. July 2006 Introduction: the plague highly infectious disease several pandemics have influenced the human history

More information

Algebraic geometry for shallow capillary-gravity waves

Algebraic geometry for shallow capillary-gravity waves Algebraic geometry for shallow capillary-gravity waves DENYS DUTYKH 1 Chargé de Recherche CNRS 1 Université de Savoie Laboratoire de Mathématiques (LAMA) 73376 Le Bourget-du-Lac France Seminar of Computational

More information

A short story about FRAP (fluourescence recovery after photobleaching)

A short story about FRAP (fluourescence recovery after photobleaching) 1 A short story about FRAP (fluourescence recovery after photobleaching) Bob Pego (Carnegie Mellon) work involving team leader: James McNally (NCI Lab of Receptor Biology & Gene Expression) Biophysicists:

More information

ON THE FIRST TIME THAT AN ITO PROCESS HITS A BARRIER

ON THE FIRST TIME THAT AN ITO PROCESS HITS A BARRIER ON THE FIRST TIME THAT AN ITO PROCESS HITS A BARRIER GERARDO HERNANDEZ-DEL-VALLE arxiv:1209.2411v1 [math.pr] 10 Sep 2012 Abstract. This work deals with first hitting time densities of Ito processes whose

More information

A Stochastic Model for TCP with Stationary Random Losses

A Stochastic Model for TCP with Stationary Random Losses A Stochastic Model for TCP with Stationary Random Losses Eitan Altman, Kostya Avrachenkov Chadi Barakat INRIA Sophia Antipolis - France ACM SIGCOMM August 31, 2000 Stockholm, Sweden Introduction Outline

More information

Latent voter model on random regular graphs

Latent voter model on random regular graphs Latent voter model on random regular graphs Shirshendu Chatterjee Cornell University (visiting Duke U.) Work in progress with Rick Durrett April 25, 2011 Outline Definition of voter model and duality with

More information

Brownian motion and the Central Limit Theorem

Brownian motion and the Central Limit Theorem Brownian motion and the Central Limit Theorem Amir Bar January 4, 3 Based on Shang-Keng Ma, Statistical Mechanics, sections.,.7 and the course s notes section 6. Introduction In this tutorial we shall

More information

Contents. Preface. Notation

Contents. Preface. Notation Contents Preface Notation xi xv 1 The fractional Laplacian in one dimension 1 1.1 Random walkers with constant steps.............. 1 1.1.1 Particle number density distribution.......... 2 1.1.2 Numerical

More information

CTRW Limits: Governing Equations and Fractal Dimensions

CTRW Limits: Governing Equations and Fractal Dimensions CTRW Limits: Governing Equations and Fractal Dimensions Erkan Nane DEPARTMENT OF MATHEMATICS AND STATISTICS AUBURN UNIVERSITY August 19-23, 2013 Joint work with Z-Q. Chen, M. D Ovidio, M.M. Meerschaert,

More information

Cauchy Problems in Bounded Domains and Iterated Processes

Cauchy Problems in Bounded Domains and Iterated Processes Cauchy Problems in Bounded Domains and Iterated Processes Erkan Nane DEPARTMENT OF MATHEMATICS AND STATISTICS AUBURN UNIVERSITY August 18, 21 Continuous time random walks R d X 2 X 3 X 1 J 1 J 2 J 3 J

More information

Chain differentials with an application to the mathematical fear operator

Chain differentials with an application to the mathematical fear operator Chain differentials with an application to the mathematical fear operator Pierre Bernhard I3S, University of Nice Sophia Antipolis and CNRS, ESSI, B.P. 145, 06903 Sophia Antipolis cedex, France January

More information

Jump-type Levy Processes

Jump-type Levy Processes Jump-type Levy Processes Ernst Eberlein Handbook of Financial Time Series Outline Table of contents Probabilistic Structure of Levy Processes Levy process Levy-Ito decomposition Jump part Probabilistic

More information

Upper and lower bounds for ruin probability

Upper and lower bounds for ruin probability Upper and lower bounds for ruin probability E. Pancheva,Z.Volkovich and L.Morozensky 3 Institute of Mathematics and Informatics, the Bulgarian Academy of Sciences, 3 Sofia, Bulgaria pancheva@math.bas.bg

More information

Anomalous Lévy diffusion: From the flight of an albatross to optical lattices. Eric Lutz Abteilung für Quantenphysik, Universität Ulm

Anomalous Lévy diffusion: From the flight of an albatross to optical lattices. Eric Lutz Abteilung für Quantenphysik, Universität Ulm Anomalous Lévy diffusion: From the flight of an albatross to optical lattices Eric Lutz Abteilung für Quantenphysik, Universität Ulm Outline 1 Lévy distributions Broad distributions Central limit theorem

More information

Winning quick and dirty: the greedy random walk

Winning quick and dirty: the greedy random walk Winning quick and dirty: the greedy random walk E. Ben-Naim and S. Redner Theory Division and Center for Nonlinear Studies, Los Alamos National Laboratory, Los Alamos, New Meico 87545 USA As a strategy

More information

A Note on the Probabilistic Evolution for P Systems

A Note on the Probabilistic Evolution for P Systems A Note on the Probabilistic Evolution for P Systems Sergey Verlan Laboratoire d Algorithmique, Complexité et Logique, Université Paris Est Créteil Val de Marne, 61, av. gén. de Gaulle, 94010, Créteil,

More information

The wave model of metric spaces

The wave model of metric spaces arxiv:1901.04317v1 [math.fa] 10 Jan 019 The wave model of metric spaces M. I. Belishev, S. A. Simonov Abstract Let Ω be a metric space, A t denote the metric neighborhood of the set A Ω of the radius t;

More information

This is a Gaussian probability centered around m = 0 (the most probable and mean position is the origin) and the mean square displacement m 2 = n,or

This is a Gaussian probability centered around m = 0 (the most probable and mean position is the origin) and the mean square displacement m 2 = n,or Physics 7b: Statistical Mechanics Brownian Motion Brownian motion is the motion of a particle due to the buffeting by the molecules in a gas or liquid. The particle must be small enough that the effects

More information

From discrete time random walks to numerical methods for fractional order differential equations

From discrete time random walks to numerical methods for fractional order differential equations From discrete time random walks to numerical methods for fractional order differential equations Click to edit Present s Name Dr. Christopher Angstmann SANUM 2016 Never Stand Still Faculty of Science School

More information

PERSISTENT RANDOM WALK EFFECT IN THE SUBDIFFUSION-REACTION PROCESS

PERSISTENT RANDOM WALK EFFECT IN THE SUBDIFFUSION-REACTION PROCESS Vol. 45 (214) ACTA PHYSICA POLONICA B No 9 PERSISTENT RANDOM WALK EFFECT IN THE SUBDIFFUSION-REACTION PROCESS Tadeusz Kosztołowicz, Mateusz Piwnik Institute of Physics, Jan Kochanowski University Świętokrzyska

More information

Local time path integrals and their application to Lévy random walks

Local time path integrals and their application to Lévy random walks Local time path integrals and their application to Lévy random walks Václav Zatloukal (www.zatlovac.eu) Faculty of Nuclear Sciences and Physical Engineering Czech Technical University in Prague talk given

More information

COMPLEX HERMITE POLYNOMIALS: FROM THE SEMI-CIRCULAR LAW TO THE CIRCULAR LAW

COMPLEX HERMITE POLYNOMIALS: FROM THE SEMI-CIRCULAR LAW TO THE CIRCULAR LAW Serials Publications www.serialspublications.com OMPLEX HERMITE POLYOMIALS: FROM THE SEMI-IRULAR LAW TO THE IRULAR LAW MIHEL LEDOUX Abstract. We study asymptotics of orthogonal polynomial measures of the

More information

variability and their application for environment protection

variability and their application for environment protection Main factors of climate variability and their application for environment protection problems in Siberia Vladimir i Penenko & Elena Tsvetova Institute of Computational Mathematics and Mathematical Geophysics

More information

Kinetics of the scavenger reaction

Kinetics of the scavenger reaction J. Phys. A: Math. Gen. 17 (1984) L451-L455. Printed in Great Britain LETTER TO THE EDITOR Kinetics of the scavenger reaction S Redner and K Kang Center for Polymer Studiest and Department of Physics, Boston

More information

Anomalous Transport and Fluctuation Relations: From Theory to Biology

Anomalous Transport and Fluctuation Relations: From Theory to Biology Anomalous Transport and Fluctuation Relations: From Theory to Biology Aleksei V. Chechkin 1, Peter Dieterich 2, Rainer Klages 3 1 Institute for Theoretical Physics, Kharkov, Ukraine 2 Institute for Physiology,

More information

How to beat a random walk with a clock?

How to beat a random walk with a clock? Introduction Moody walk analysis R/A Conclusion How to beat a random walk with a clock? Locating a Target With an Agent Guided by Unreliable Local Advice Nicolas Hanusse 1 David Ilcinkas 1 Adrian Kosowski

More information

A Feynman-Kac Path-Integral Implementation for Poisson s Equation

A Feynman-Kac Path-Integral Implementation for Poisson s Equation for Poisson s Equation Chi-Ok Hwang and Michael Mascagni Department of Computer Science, Florida State University, 203 Love Building Tallahassee, FL 32306-4530 Abstract. This study presents a Feynman-Kac

More information

Random walks in a 1D Levy random environment

Random walks in a 1D Levy random environment Random Combinatorial Structures and Statistical Mechanics, University of Warwick, Venice, 6-10 May, 2013 Random walks in a 1D Levy random environment Alessandra Bianchi Mathematics Department, University

More information

BRANCHING PROCESSES AND THEIR APPLICATIONS: Lecture 15: Crump-Mode-Jagers processes and queueing systems with processor sharing

BRANCHING PROCESSES AND THEIR APPLICATIONS: Lecture 15: Crump-Mode-Jagers processes and queueing systems with processor sharing BRANCHING PROCESSES AND THEIR APPLICATIONS: Lecture 5: Crump-Mode-Jagers processes and queueing systems with processor sharing June 7, 5 Crump-Mode-Jagers process counted by random characteristics We give

More information

q P (T b < T a Z 0 = z, X 1 = 1) = p P (T b < T a Z 0 = z + 1) + q P (T b < T a Z 0 = z 1)

q P (T b < T a Z 0 = z, X 1 = 1) = p P (T b < T a Z 0 = z + 1) + q P (T b < T a Z 0 = z 1) Random Walks Suppose X 0 is fixed at some value, and X 1, X 2,..., are iid Bernoulli trials with P (X i = 1) = p and P (X i = 1) = q = 1 p. Let Z n = X 1 + + X n (the n th partial sum of the X i ). The

More information

Outline of the talk How to describe restricted diffusion? How to monitor restricted diffusion? Laplacian eigenfunctions in NMR Other applications Loca

Outline of the talk How to describe restricted diffusion? How to monitor restricted diffusion? Laplacian eigenfunctions in NMR Other applications Loca Laplacian Eigenfunctions in NMR Denis S. Grebenkov Laboratoire de Physique de la Matière Condensée CNRS Ecole Polytechnique, Palaiseau, France IPAM Workshop «Laplacian Eigenvalues and Eigenfunctions» February

More information

Fig 1: Stationary and Non Stationary Time Series

Fig 1: Stationary and Non Stationary Time Series Module 23 Independence and Stationarity Objective: To introduce the concepts of Statistical Independence, Stationarity and its types w.r.to random processes. This module also presents the concept of Ergodicity.

More information

Moment analysis of the telegraph random process

Moment analysis of the telegraph random process BULETINUL ACADEMIEI DE ŞTIINŢE A REPUBLICII MOLDOVA. MATEMATICA Number 68), 0, Pages 90 07 ISSN 04 7696 Moment analysis of the telegraph random process Alexander D.Kolesnik Abstract. We consider the Goldstein-Kac

More information

Re-entrant transition in a one-dimensional growth model

Re-entrant transition in a one-dimensional growth model Re-entrant transition in a one-dimensional growth model Ostap Hryniv Dept. of Mathematical Sciences Durham University Ostap.Hryniv@durham.ac.uk X International Conference Stochastic and Analytic Methods

More information

Local vs. Nonlocal Diffusions A Tale of Two Laplacians

Local vs. Nonlocal Diffusions A Tale of Two Laplacians Local vs. Nonlocal Diffusions A Tale of Two Laplacians Jinqiao Duan Dept of Applied Mathematics Illinois Institute of Technology Chicago duan@iit.edu Outline 1 Einstein & Wiener: The Local diffusion 2

More information

CRRAO Advanced Institute of Mathematics, Statistics and Computer Science (AIMSCS) Research Report. B. L. S. Prakasa Rao

CRRAO Advanced Institute of Mathematics, Statistics and Computer Science (AIMSCS) Research Report. B. L. S. Prakasa Rao CRRAO Advanced Institute of Mathematics, Statistics and Computer Science (AIMSCS) Research Report Author (s): B. L. S. Prakasa Rao Title of the Report: FILTERED FRACTIONAL POISSON ROCESSES Research Report

More information

arxiv: v1 [cond-mat.stat-mech] 6 Mar 2008

arxiv: v1 [cond-mat.stat-mech] 6 Mar 2008 CD2dBS-v2 Convergence dynamics of 2-dimensional isotropic and anisotropic Bak-Sneppen models Burhan Bakar and Ugur Tirnakli Department of Physics, Faculty of Science, Ege University, 35100 Izmir, Turkey

More information

Spectral functions of subordinated Brownian motion

Spectral functions of subordinated Brownian motion Spectral functions of subordinated Brownian motion M.A. Fahrenwaldt 12 1 Institut für Mathematische Stochastik Leibniz Universität Hannover, Germany 2 EBZ Business School, Bochum, Germany Berlin, 23 October

More information

Numerical Algorithm for Optimal Control of Continuity Equations

Numerical Algorithm for Optimal Control of Continuity Equations Numerical Algorithm for Optimal Control of Continuity Equations Nikolay Pogodaev Matrosov Institute for System Dynamics and Control Theory Lermontov str., 134 664033 Irkutsk, Russia Krasovskii Institute

More information

Lévy Walks and scaling in quenched disordered media

Lévy Walks and scaling in quenched disordered media Lévy Walks and scaling in quenched disordered media A. Vezzani, - CNR Modena - Parma L. Caniparoli - Sissa Trieste S.Lepri - Firenze CNR - ISC Raffaella Burioni - Parma Sperlonga - September 2010 Work

More information

Normal and Anomalous Diffusion: A Tutorial

Normal and Anomalous Diffusion: A Tutorial Normal and Anomalous Diffusion: A Tutorial Loukas Vlahos, Heinz Isliker Department of Physics, University of Thessaloniki, 54124 Thessaloniki, Greece Yannis Kominis, and Kyriakos Hizanidis School of Electrical

More information

Biased Random Walks in Biology

Biased Random Walks in Biology Biased Random Walks in Biology Edward Alexander Codling Submitted in accordance with the requirements for the degree of Doctor of Philosophy The University of Leeds, Department of Applied Mathematics.

More information

Controlling chaotic transport in Hamiltonian systems

Controlling chaotic transport in Hamiltonian systems Controlling chaotic transport in Hamiltonian systems Guido Ciraolo Facoltà di Ingegneria, Università di Firenze via S. Marta, I-50129 Firenze, Italy Cristel Chandre, Ricardo Lima, Michel Vittot CPT-CNRS,

More information

Asymptotics for posterior hazards

Asymptotics for posterior hazards Asymptotics for posterior hazards Pierpaolo De Blasi University of Turin 10th August 2007, BNR Workshop, Isaac Newton Intitute, Cambridge, UK Joint work with Giovanni Peccati (Université Paris VI) and

More information

Cauchy Problems Solved by Running Subordinate Processes

Cauchy Problems Solved by Running Subordinate Processes Cauchy Problems Solved by Running Subordinate Processes Erkan Nane DEPARTMENT OF MATHEMATICS AND STATISTICS AUBURN UNIVERSITY August 18, 21 Continuous time random walks R d X 2 X 3 X 1 J 1 J 2 J 3 J 4

More information

Random Walks and Diffusion. APC 514 December 5, 2002 Cox&Shvartsman

Random Walks and Diffusion. APC 514 December 5, 2002 Cox&Shvartsman Random Walks and Diffusion APC 514 December 5, 2002 Cox&Shvartsman Cell motility over adhesive substrates: a periodic phenomenon Simple model for linear motion: DIMILLA PA, BARBEE K, LAUFFENBURGER DA MATHEMATICAL-MODEL

More information

On angular momentum balance for particle systems with periodic boundary conditions

On angular momentum balance for particle systems with periodic boundary conditions On angular momentum balance for particle systems with periodic boundary conditions Vitaly A. Kuzkin 1 Institute for Problems in Mechanical Engineering RAS Saint Petersburg State Polytechnical University

More information

A mathematical model for a copolymer in an emulsion

A mathematical model for a copolymer in an emulsion J Math Chem (2010) 48:83 94 DOI 10.1007/s10910-009-9564-y ORIGINAL PAPER A mathematical model for a copolymer in an emulsion F. den Hollander N. Pétrélis Received: 3 June 2007 / Accepted: 22 April 2009

More information

Diffusion in a Logarithmic Potential Anomalies, Aging & Ergodicity Breaking

Diffusion in a Logarithmic Potential Anomalies, Aging & Ergodicity Breaking Diffusion in a Logarithmic Potential Anomalies, Aging & Ergodicity Breaking David Kessler Bar-Ilan Univ. E. Barkai (Bar-Ilan) A. Dechant (Augsburg) E. Lutz (Augsburg) PRL, 105 120602 (2010) arxiv:1105.5496

More information

A State Space Model for Wind Forecast Correction

A State Space Model for Wind Forecast Correction A State Space Model for Wind Forecast Correction Valrie Monbe, Pierre Ailliot 2, and Anne Cuzol 1 1 Lab-STICC, Université Européenne de Bretagne, France (e-mail: valerie.monbet@univ-ubs.fr, anne.cuzol@univ-ubs.fr)

More information

Matrix-valued stochastic processes

Matrix-valued stochastic processes Matrix-valued stochastic processes and applications Małgorzata Snarska (Cracow University of Economics) Grodek, February 2017 M. Snarska (UEK) Matrix Diffusion Grodek, February 2017 1 / 18 Outline 1 Introduction

More information

Diffusion of a density in a static fluid

Diffusion of a density in a static fluid Diffusion of a density in a static fluid u(x, y, z, t), density (M/L 3 ) of a substance (dye). Diffusion: motion of particles from places where the density is higher to places where it is lower, due to

More information

Lie Symmetry of Ito Stochastic Differential Equation Driven by Poisson Process

Lie Symmetry of Ito Stochastic Differential Equation Driven by Poisson Process American Review of Mathematics Statistics June 016, Vol. 4, No. 1, pp. 17-30 ISSN: 374-348 (Print), 374-356 (Online) Copyright The Author(s).All Rights Reserved. Published by American Research Institute

More information

AGAT 2016, Cargèse a point-vortex toy model

AGAT 2016, Cargèse a point-vortex toy model AGAT 2016, Cargèse a point-vortex toy model Jean-François Pinton CNRS & ENS de Lyon M.P. Rast, JFP, PRE 79 (2009) M.P. Rast, JFP, PRL 107 (2011) M.P. Rast, JFP, P.D. Mininni, PRE 93 (2009) Motivations

More information

Solution for Problem 7.1. We argue by contradiction. If the limit were not infinite, then since τ M (ω) is nondecreasing we would have

Solution for Problem 7.1. We argue by contradiction. If the limit were not infinite, then since τ M (ω) is nondecreasing we would have 362 Problem Hints and Solutions sup g n (ω, t) g(ω, t) sup g(ω, s) g(ω, t) µ n (ω). t T s,t: s t 1/n By the uniform continuity of t g(ω, t) on [, T], one has for each ω that µ n (ω) as n. Two applications

More information

Convergence of Banach valued stochastic processes of Pettis and McShane integrable functions

Convergence of Banach valued stochastic processes of Pettis and McShane integrable functions Convergence of Banach valued stochastic processes of Pettis and McShane integrable functions V. Marraffa Department of Mathematics, Via rchirafi 34, 90123 Palermo, Italy bstract It is shown that if (X

More information

Motion in Space. MATH 311, Calculus III. J. Robert Buchanan. Fall Department of Mathematics. J. Robert Buchanan Motion in Space

Motion in Space. MATH 311, Calculus III. J. Robert Buchanan. Fall Department of Mathematics. J. Robert Buchanan Motion in Space Motion in Space MATH 311, Calculus III J. Robert Buchanan Department of Mathematics Fall 2011 Background Suppose the position vector of a moving object is given by r(t) = f (t), g(t), h(t), Background

More information

arxiv:cond-mat/ v2 [cond-mat.dis-nn] 11 Oct 2002

arxiv:cond-mat/ v2 [cond-mat.dis-nn] 11 Oct 2002 arxiv:cond-mat/21166v2 [cond-mat.dis-nn] 11 Oct 22 REVISITING THE DERIVATION OF THE FRACTIONAL DIFFUSION EQUATION ENRICO SCALAS Dipartimento di Scienze e Tecnologie Avanzate, Università del Piemonte Orientale

More information

Branching Random Walks with Immigration

Branching Random Walks with Immigration UDC 519.21 Branching Random Walks with Immigration D. Han, Yu. Makarova, S. Molchanov, E. Yarovaya Department of Mathematics and Statistics, University of North Carolina at Charlotte, Charlotte, NC 28223,

More information

Anomalous diffusion of volcanic earthquakes

Anomalous diffusion of volcanic earthquakes Anomalous diffusion of volcanic earthquakes SUMIYOSHI ABE 1,2 and NORIKAZU SUZUKI 3 1 Department of Physical Engineering, Mie University, Mie 514-8507, Japan 2 Institute of Physics, Kazan Federal University,

More information