Brownian motion 18.S995 - L03 & 04

Similar documents
Topics covered so far:

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

Introduction to the mathematical modeling of multi-scale phenomena

For slowly varying probabilities, the continuum form of these equations is. = (r + d)p T (x) (u + l)p D (x) ar x p T(x, t) + a2 r

Anomalous diffusion in cells: experimental data and modelling

Continuous and Discrete random process

Biological motors 18.S995 - L10

Gaussian Random Variables Why we Care

Brownian motion and the Central Limit Theorem

Weak Ergodicity Breaking. Manchester 2016

ROBERT BROWN AND THE POLLEN STUFF

Handbook of Stochastic Methods

Local vs. Nonlocal Diffusions A Tale of Two Laplacians

A short story about FRAP (fluourescence recovery after photobleaching)

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

16. Working with the Langevin and Fokker-Planck equations

Stationary independent increments. 1. Random changes of the form X t+h X t fixed h > 0 are called increments of the process.

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

Cover times of random searches

Gaussian, Markov and stationary processes

Lecture 3: From Random Walks to Continuum Diffusion

The dynamics of small particles whose size is roughly 1 µmt or. smaller, in a fluid at room temperature, is extremely erratic, and is

Weak Ergodicity Breaking

Weak Ergodicity Breaking WCHAOS 2011

Anomalous Transport and Fluctuation Relations: From Theory to Biology

MAE 545: Lecture 2 (9/22) E. coli chemotaxis

Handbook of Stochastic Methods

Anatoly B. Kolomeisky. Department of Chemistry CAN WE UNDERSTAND THE COMPLEX DYNAMICS OF MOTOR PROTEINS USING SIMPLE STOCHASTIC MODELS?

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

Random walks. Marc R. Roussel Department of Chemistry and Biochemistry University of Lethbridge. March 18, 2009

Escape Problems & Stochastic Resonance

Stochastic differential equation models in biology Susanne Ditlevsen

Random walk through Anomalous processes: some applications of Fractional calculus

Properties of the Autocorrelation Function

STOCHASTIC PROCESSES IN PHYSICS AND CHEMISTRY

First-Passage Statistics of Extreme Values

Non equilibrium thermodynamics: foundations, scope, and extension to the meso scale. Miguel Rubi

Lecture 13. Drunk Man Walks

Some Tools From Stochastic Analysis

Lecture 1: Random walk

Diffusion and cellular-level simulation. CS/CME/BioE/Biophys/BMI 279 Nov. 7 and 9, 2017 Ron Dror

t For l = 1 a monomer cannot be destroyed or created from nothing: = b p(2,t) a p(1,t).

CHAPTER V. Brownian motion. V.1 Langevin dynamics

Random Processes Why we Care

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

EQUATION LANGEVIN. Physics, Chemistry and Electrical Engineering. World Scientific. With Applications to Stochastic Problems in. William T.

Smoluchowski Symposia: Why are we here?

Problems in diffusion and absorption: How fast can you hit a target with a random walk?

1 Introduction. 2 Diffusion equation and central limit theorem. The content of these notes is also covered by chapter 3 section B of [1].

Chapter 3. Random Process & Partial Differential Equations

2. Molecules in Motion

Particle displacements by swimming organisms

BioNewMetrics and stochastic analysis of superresolution. David Holcman Ecole Normale Superieure

The Polymers Tug Back

GEM4 Summer School OpenCourseWare

Introduction to Theory of Mesoscopic Systems. Boris Altshuler Columbia University & NEC Laboratories America

Chapter 6 - Random Processes

Electronics Labs - Digital Electronics

Corrections to Theory of Asset Pricing (2008), Pearson, Boston, MA

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

Lecture The Sample Mean and the Sample Variance Under Assumption of Normality

Deterministic. Deterministic data are those can be described by an explicit mathematical relationship

Introduction to Machine Learning CMU-10701

What s more chaotic than chaos itself? Brownian Motion - before, after, and beyond.

Infinitely divisible distributions and the Lévy-Khintchine formula

MD Thermodynamics. Lecture 12 3/26/18. Harvard SEAS AP 275 Atomistic Modeling of Materials Boris Kozinsky

Beyond the Gaussian universality class

3 Biopolymers Uncorrelated chains - Freely jointed chain model

A Hilbert Space for Random Processes

Lecture 3: Central Limit Theorem

Stochastic-Process Limits

Jerry Gollub 7/18/2011

Monte Carlo Integration I [RC] Chapter 3

Exercises. T 2T. e ita φ(t)dt.

Diffusion in a Logarithmic Potential Anomalies, Aging & Ergodicity Breaking

Lecture 9. d N(0, 1). Now we fix n and think of a SRW on [0,1]. We take the k th step at time k n. and our increments are ± 1

Generating and characteristic functions. Generating and Characteristic Functions. Probability generating function. Probability generating function

Theory of Mesoscopic Systems

Scaling exponents for certain 1+1 dimensional directed polymers

Search for Food of Birds, Fish and Insects

Statistical Thermodynamics and Monte-Carlo Evgenii B. Rudnyi and Jan G. Korvink IMTEK Albert Ludwig University Freiburg, Germany

Statistical Mechanics of Active Matter

When is an Integrate-and-fire Neuron like a Poisson Neuron?

Diffusion. CS/CME/BioE/Biophys/BMI 279 Nov. 15 and 20, 2016 Ron Dror

Intelligent Embedded Systems Uncertainty, Information and Learning Mechanisms (Part 1)

Fractional Laplacian

arxiv: v1 [cond-mat.stat-mech] 29 Jun 2015

Stochastic Processes at Single-molecule and Single-cell levels

Continuum Limit of Forward Kolmogorov Equation Friday, March 06, :04 PM

EXPERIMENT 17. To Determine Avogadro s Number by Observations on Brownian Motion. Introduction

Random Walks A&T and F&S 3.1.2

Other properties of M M 1

Direct observation of aggregative nanoparticle growth: Kinetic modeling of the size distribution and growth rate

1.1 Definition of BM and its finite-dimensional distributions

CHE3935. Lecture 10 Brownian Dynamics

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

Stochastic Processes

1 IEOR 6712: Notes on Brownian Motion I

FE 5204 Stochastic Differential Equations

Weak chaos, infinite ergodic theory, and anomalous diffusion

Transcription:

Brownian motion 18.S995 - L03 & 04

Typical length scales http://www2.estrellamountain.edu/faculty/farabee/biobk/biobookcell2.html dunkel@math.mit.edu

Brownian motion

Brownian motion Jan Ingen-Housz (1730-1799) 1784/1785: http://www.physik.uni-augsburg.de/theo1/hanggi/history/bm-history.html

Robert Brown (1773-1858) Linnean society, London 1827: irreguläre Eigenbewegung von Pollen in Flüssigkeit irregular motion of pollen in fluid http://www.brianjford.com/wbbrownc.htm

!"#$%&'()*+,-#./0102/3//4 5"#67,8&(7,#./0932/3114 :"#$;<*%='<>8?7#./09@2/3/94!"#$%&'!((()*+&,-&)%".),#!"#$%&'!(/0/1&2/,)"$-!"#$%&'!(/0/1&2/,)"$- D! 2 RT RT 1 32 mc D! D! 6"#aC N 6" k P 243 "$ R A'7*"#:+B"#!C#90/#./3D14 5,,"#A'E8"#"#C#1F3#./3D14 5,,"#A'E8"#$"C#91G#./3DG4

Jean Baptiste Perrin (1870-1942, Nobelpreis prize 1926) colloidal particles of radius 0.53µm Mouvement brownien et réalité moléculaire, Annales de chimie et de physique VIII 18, 5-114 (1909) successive positions every 30 seconds joined by straight line segments mesh size is 3.2µm Les Atomes, Paris, Alcan (1913) Experimenteller Nachweis der atomistischen Struktur der Materie experimental evidence for atomistic structure of matter

Norbert Wiener (1894-1864) MIT

Relevance in biology intracellular transport intercellular transport microorganisms must beat BM to achieve directed locomotion tracer diffusion = important experimental tool generalized BMs (polymers, membranes, etc.) dunkel@math.mit.edu

Polymers & filaments (D=1) Dogic Lab, Brandeis Physical parameters (e.g. bending rigidity) from fluctuation Drosophila oocyte analysis Goldstein lab, PNAS 2012 dunkel@math.mit.edu

Brownian tracer particles in a bacterial suspension Bacillus subtilis Tracer colloids PRL 2013

http://web.mit.edu/mbuehler/www/research/f103.jpg

Basic idea Split dynamics into deterministic part (drift) random part (diffusion)

Typical problems Determine noise structure transport coefficients first passage (escape) times D http://www.pnas.org/content/104/41/16098/f1.expansion.html

Probability space A B [0, 1] P[;] =0

Expectation values of discrete random variables

Expectation values of continuous random variables

Random walk model X N = x 0 + ` NX S i

1.1 Random walks 1.1.1 Unbiased random walk (RW) Consider the one-dimensional unbiased RW (fixed initial position X 0 length `) = x 0, N steps of X N = x 0 + ` NX S i (1.1) where S i 2 {±1} are iid. random variables (RVs) with P[S i = ±1] = 1/2. Noting that 1 E[S i ] = 1 1 2 +1 1 =0, (1.2) 2 apple E[S i S j ] = ij E[Si 2 ]= ij ( 1) 2 1 2 +(1)2 1 2 = ij, (1.3) we find for the first moment of the RW E[X N ] = x 0 + ` NX E[S i ]=x 0 (1.4) R

1.1 Random walks 1.1.1 Unbiased random walk (RW) Consider the one-dimensional unbiased RW (fixed initial position X 0 length `) = x 0, N steps of X N = x 0 + ` NX S i (1.1) where S i 2 {±1} are iid. random variables (RVs) with P[S i = ±1] = 1/2. Noting that 1 E[S i ] = 1 1 2 +1 1 =0, (1.2) 2 apple E[S i S j ] = ij E[Si 2 ]= ij ( 1) 2 1 2 +(1)2 1 2 = ij, (1.3) we find for the first moment of the RW E[X N ] = x 0 + ` NX E[S i ]=x 0 (1.4) R

1.1 Random walks 1.1.1 Unbiased random walk (RW) Consider the one-dimensional unbiased RW (fixed initial position X 0 length `) = x 0, N steps of X N = x 0 + ` NX S i (1.1) where S i 2 {±1} are iid. random variables (RVs) with P[S i = ±1] = 1/2. Noting that 1 E[S i ] = 1 1 2 +1 1 =0, (1.2) 2 apple E[S i S j ] = ij E[Si 2 ]= ij ( 1) 2 1 2 +(1)2 1 2 = ij, (1.3) we find for the first moment of the RW E[X N ] = x 0 + ` NX E[S i ]=x 0 (1.4) R

Second moment (uncentered) NX E[XN] 2 = E[(x 0 + ` S i ) 2 ] = E[x 2 0 +2x 0` NX NX NX S i + `2 S i S j ] j=1 j=1 NX NX = x 2 0 +2x 0 0+`2 E[S i S j ] NX NX = x 2 0 +2x 0 0+`2 j=1 ij = x 2 0 + `2N. (1.5)

Second moment (uncentered) NX E[XN] 2 = E[(x 0 + ` S i ) 2 ] = E[x 2 0 +2x 0` NX NX NX S i + `2 S i S j ] j=1 j=1 NX NX = x 2 0 +2x 0 0+`2 E[S i S j ] NX NX = x 2 0 +2x 0 0+`2 j=1 ij = x 2 0 + `2N. (1.5)

Second moment (uncentered) NX E[XN] 2 = E[(x 0 + ` S i ) 2 ] = E[x 2 0 +2x 0` NX NX NX S i + `2 S i S j ] j=1 j=1 NX NX = x 2 0 +2x 0 0+`2 E[S i S j ] NX NX = x 2 0 +2x 0 0+`2 j=1 ij = x 2 0 + `2N. (1.5)

Second moment (uncentered) NX E[XN] 2 = E[(x 0 + ` S i ) 2 ] = E[x 2 0 +2x 0` NX NX NX S i + `2 S i S j ] j=1 j=1 NX NX = x 2 0 +2x 0 0+`2 E[S i S j ] NX NX = x 2 0 +2x 0 0+`2 j=1 ij = x 2 0 + `2N. (1.5)

Second moment (uncentered) NX E[XN] 2 = E[(x 0 + ` S i ) 2 ] = E[x 2 0 +2x 0` NX NX NX S i + `2 S i S j ] j=1 j=1 NX NX = x 2 0 +2x 0 0+`2 E[S i S j ] NX NX = x 2 0 +2x 0 0+`2 j=1 ij = x 2 0 + `2N. (1.5)

Continuum limit X N = x 0 + ` NX S i Let

Continuum limit P (N,K) = = N 1 N N 2 1 2 K 2 N N! ((N + K)/2)! ((N K)/2)!. (1.8) The associated probability density function (PDF) can be found by defining p(t, x) := P (N,K) 2` and considering limit, `! 0suchthat = P (t/,x/`) 2` (1.9) D := `2 2 = const, (1.10)

Continuum limit yielding the Gaussian p(t, x) ' r 1 x 2 4 Dt exp 4Dt (1.11) (pset1) Eq. (1.11) is the fundamental solution to the di usion equation, @ t p t = D@ xx p, (1.12) where @ t, @ x, @ xx,... denote partial derivatives. The mean square displacement of the continuous process described by Eq. (1.11) is Z E[X(t) 2 ]= dx x 2 p(t, x) =2Dt, (1.13) in agreement with Eq. (1.7).

Z Remark One often classifies di usion processes by the (asymptotic) power-law growth of the mean square displacement, µ =0:Staticprocesswithnomovement. E[(X(t) X(0)) 2 ] t µ. (1.14) 0 <µ<1:sub-di usion, arises typically when waiting times between subsequent jumps can be long and/or in the presence of a su ciently large number of obstacles (e.g. slow di usion of molecules in crowded cells). µ = 1 : Normal di usion, corresponds to the regime governed by the standard Central Limit Theorem (CLT). 1 <µ<2: Super-di usion, occurs when step-lengths are drawn from distributions with infinite variance (Lévy walks; considered as models of bird or insect movements). µ =2:Ballisticpropagation(deterministicwave-likeprocess).

Brownian motion non-brownian Levy-flight

1.1.2 Biased random walk (BRW) Consider a one-dimensional hopping process on a discrete lattice (spacing `), defined such that during a time-step aparticleatpositionx(t) =`j 2 `Z can either (i) jump a fixed distance ` to the left with probability,or (ii) jump a fixed distance ` to the right with probability, or (iii) remain at its position x with probability (1 ). Assuming that the process is Markovian (does not depend on the past), the evolution of the associated probability vector P (t) =(P (t, x)) = (P j (t)), where x = `j, is governed by the master equation P (t +,x)=(1 ) P (t, x)+ P (t, x `)+ P (t, x + `). (1.15)

Master equations P (t +,x)=(1 ) P (t, x)+ P (t, x `)+ P (t, x + `). (1.15) Technically,, and (1 ) arethenon-zero-elementsofthecorrespondingtransition matrix W =(W ij )withw ij > 0thatgovernstheevolutionofthecolumnprobability vector P (t) =(P j (t)) = (P (t, y)) by P i (t + ) =W ij P j (t) (1.16a) or, more generally, for n steps P (t + n ) =W n P (t). (1.16b) The stationary solutions are the eigenvectors of W with eigenvalue 1. To preserve normalization, one requires P i W ij =1.

Continuum limit Define the density p(t,x) = P (t,x)/`. Assume, ` are small, so that we can Taylor-expand p(t +,x) ' p(t, x)+ @ t p(t, x) (1.17a) p(t, x ± `) ' p(t, x) ± `@ x p(t, x)+`2 2 @ xxp(t, x) (1.17b) Neglecting the higher-order terms, it follows from Eq. (1.15) that p(t, x)+ @ t p(t, x) ' (1 ) p(t, x)+ Dividing by, one obtains the advection-di usion equation with drift velocity u and di usion constant D given by 2 [p(t, x) `@ x p(t, x)+`2 2 @ xxp(t, x)] + [p(t, x)+`@ x p(t, x)+`2 2 @ xxp(t, x)]. (1.18) @ t p = u @ x p + D @ xx p (1.19a) u := ( ) `, D := ( + ) `2 2. (1.19b)

Continuum limit Define the density p(t,x) = P (t,x)/`. Assume, ` are small, so that we can Taylor-expand p(t +,x) ' p(t, x)+ @ t p(t, x) (1.17a) p(t, x ± `) ' p(t, x) ± `@ x p(t, x)+`2 2 @ xxp(t, x) (1.17b) Neglecting the higher-order terms, it follows from Eq. (1.15) that p(t, x)+ @ t p(t, x) ' (1 ) p(t, x)+ Dividing by, one obtains the advection-di usion equation with drift velocity u and di usion constant D given by 2 [p(t, x) `@ x p(t, x)+`2 2 @ xxp(t, x)] + [p(t, x)+`@ x p(t, x)+`2 2 @ xxp(t, x)]. (1.18) @ t p = u @ x p + D @ xx p (1.19a) u := ( ) `, D := ( + ) `2 2. (1.19b)

Continuum limit Define the density p(t,x) = P (t,x)/`. Assume, ` are small, so that we can Taylor-expand p(t +,x) ' p(t, x)+ @ t p(t, x) (1.17a) p(t, x ± `) ' p(t, x) ± `@ x p(t, x)+`2 2 @ xxp(t, x) (1.17b) Neglecting the higher-order terms, it follows from Eq. (1.15) that p(t, x)+ @ t p(t, x) ' (1 ) p(t, x)+ Dividing by, one obtains the advection-di usion equation with drift velocity u and di usion constant D given by 2 [p(t, x) `@ x p(t, x)+`2 2 @ xxp(t, x)] + [p(t, x)+`@ x p(t, x)+`2 2 @ xxp(t, x)]. (1.18) @ t p = u @ x p + D @ xx p (1.19a) u := ( ) `, D := ( + ) `2 2. (1.19b)

Time-dependent solution Dividing by, one obtains the advection-di usion equation with drift velocity u and di usion constant D given by 2 @ t p = u @ x p + D @ xx p (1.19a) u := ( ) `, D := ( + ) `2 2. We recover the classical di usion equation (1.12) from Eq. (1.19a) for = time-dependent fundamental solution of Eq. (1.19a) reads r 1 (x ut) 2 p(t, x) = 4 Dt exp 4Dt (1.19b) = 0.5. The (1.20)

r Remarks Note that Eqs. (1.12) and Eq. (1.19a) can both be written in the current-form @ t p + @ x j x =0 (1.21) with j x = up D@ x p, (1.22) reflecting conservation of probability. Another commonly-used representation is @ t p = Lp, (1.23) where L is a linear di erential operator; in the above example (1.19b) L := u @ x + D @ xx. (1.24) Stationary solutions, if they exist, are eigenfunctions of L with eigenvalue 0. (useful later when discussing Brownian motors)

2024 © sciencedocbox.com Privacy Policy | Terms of Service | Feedback