Brownian Motion: Fokker-Planck Equation
|
|
- Julian Herbert Briggs
- 6 years ago
- Views:
Transcription
1 Chapter 7 Brownian Motion: Fokker-Planck Equation The Fokker-Planck equation is the equation governing the time evolution of the probability density of the Brownian particla. It is a second order differential equation and is exact for the case when the noise acting on the Brownian particle is Gaussian white noise. A general Fokker-Planck equation can be derived from the Chapman-Kolmogorov equation, but we also like to find the Fokker-Planck equation corresponding to the time dependence given by a Langevin equation. The derivation of the Fokker-Planck equation is a two step process. We first derive the equation of motion for the probability density 4/varrho(x, v, t)4 to find the Brownian particle in the interval (x, x+dx) and (v, v+dv) at time t for one realization of the random force ξ(t). We then obtain an equation for P (x, v, t) = ϱ(x, v, t) ξ i.e. the average of ϱ(x, v, t) over many realizations of the random force. The probability density P (x, v, t) is the macroscopically observed probability density for the Brownian particel. 7.1 Probability flow in phase-space Let us obtain the probability to find the Brownian particle in the interval (x, x + dx and (v, v + dv) at time t. We will consider the space of coordiantes x = (x, v). The Brownian particle is located in the infinitesimal ara dxdv with probablity ϱ(x, v, t)dxdv. The velocity of the particle at point (x, v) is given by ẋ = (ẋ, v) and the current density is ẋϱ. Since the Brownian particle must lie somewhere in the phase-space < x <, < v < we have the condition dx dvϱ(x, v, t) = 1 Let us now consider a finite area, or volume, V in this space. Since the Brownian particle cannot be destroyed a change in the probability contained in V must be due to a flow of 89
2 9 Chapter 7 Brownian Motion: Fokker-Planck Equation probability through the surface S surrounding V. Thus d dxdvϱ(x, v, t) = ϱ(x, v, t)ẋ ds V S We can now use Gauss theorem to change the surface integral into a volume integral. dxdv ϱ(x, v, t) = dxdv (ẋϱ(x, v, t)) V t V Since V is fixed and arbitrary we find the continuity equation ϱ(x, v, t) = (ẋϱ(x, v, t)) = t x (ẋϱ(x, v, t)) ( vϱ(x, v, t)) (7.1) v This is the continuity equation in phase-space which just state that probability is conserved. 7.2 Probability flow for Brownian particle In order to write (7.1) explicitly for a Brownian particle we must know the Langevin equation governing the evolution of the particle. For a particle moving in the presence of a potnetial V (x) the Langevin equations are dx dv = v = γ m v + 1 m F (x) + 1 ξ(t) (7.2) m where the force F (x) = V (x). Inserting (7.2) into (7.1) gives ϱ(x, v, t) = t x (vϱ(x, v, t)) + γ m v (vϱ(x, v, t)) 1 m F (x) ϱ(x, v, t) v 1 m ξ(t) v ϱ(x, v, t) = L ϱ(x, v, t) L 1 (t)ϱ(x, v, t) where the differential operators L and L 1 are defined as L = v x γ m γ m v v + 1 m F (x) v L 1 = 1 m ξ(t) v Since ξ(t) is a stochastic variable the time evolution of ϱ will be different for each realization of ξ(t). However when we observe an actual Brownian particle we are observing the average effect of the random force on it. Therefore we introduce an observable probability density P (x, v, t) = ξ(t) ξ. Let ϱ(t) = e L t σ(t) then t σ(x, v, t) = e L t L 1 (t)e L t σ(x, v, t) = V (t)σ(x, v, t)
3 7.2 Probability flow for Brownian particle 91 This equation has the formal solution which follows since formally σ(t) = σ() 1 + ( 1) n 1 = σ(t) = exp 1 V (t 1 )σ(t 1 ) = σ() V (t ) σ() 1 V (t 1 )σ() + n 2 n 1 V (t 1 )V (t 2 )... V (t n 1 )σ() +... ( 1) n n 1 V (t 1 ) σ() = exp n! n= V (t ) σ() V (t 1 )V (t 2 )σ() +... The third step follows since by changing the order of integration and then varibles V (t 1 )V (t 2 ) = 2 t 2 1 V (t 1 )V (t 2 ) = 1 t 1 2 V (t 1 )V (t 2 ) so that Also assume that V (t 1 )V (t 2 ) = V (t 1 ) 1 1 then by taking the derivative n 2 n 1 V (t 1 )V (t 2 )... V (t n 1 ) = 1 n 1 V (t 1 ) (7.3) n! d = V (t) 1 1 n n n V (t 1 )V (t 2 )... V (t n ) n+1 V (t 1 )V (t 2 )... V (t n ) = V (t) 1 n 1 V (t 1 ) = d 1 t n! (n + 1)! n+1 1 V (t 1 ) By induction it therefore follows that (7.3) holds. Taking the average ξ over the Gaussian noice ξ(t) we see that σ(t) ξ is the characteristic function of the random variable X(t) = i 1V (t 1 ). This must again be a Gaussian variable with X(t) ξ = and the variance is X(t) 2 = V (t 1 )V (t 2 ) 2 Since the characteristic function for the Gaussian variable X(t) is exp (ix(t)) = exp ( iµ X X(t) 2 /2 ) we find ( 1 t ) σ(t) ξ = exp 1 2 V (t 1 )V (t 2 ) σ() (7.4) 2
4 92 Chapter 7 Brownian Motion: Fokker-Planck Equation This formula is just a special case of a cumulant expansion. The integral in (7.4) becomes V (t 1 )V (t 2 ) ξ = e L t m ξ(t 1) v e L t 1 e L t 2 1 m ξ(t 2) v e L t 2 ξ g t = 2m 2 1 e L t 1 2 v 2 e L t 1 Then t σ(x, v, t) ξ = g 2m 2 el t 2 v 2 e L t σ(x, v, t) ξ This gives for ϱ(x, v, t) ξ t ϱ(x, v, t) ξ = L ϱ(x, v, t) ξ + g 2 2m 2 v 2 ϱ(x, v, t) ξ and for the probability distribution P (x, v, t) = v t x P (x, v, t) ( γ v m v 1 ) m F (x) P (x, v, t) + g 2 2m 2 P (x, v, t) v2 (7.5) This is the Fokker-Planck equation for the probability P dxdv to find the Brownian particle in the interval (x, x + dx, (v, v + dv) at time t. We can write the Fokker-Planck equation as a continuity equation P (x, v, t) = j t where = e x /x + e v /v and the probability current is ( γ j = e x vp e v m v 1 ) m F (x) P + g 2m 2 v P 7.3 General Fokker-Planck equation We can obtain the Fokker-Planck equation for a quite general Langevin equation for the dynamics of a set of fluctuating variables a = {a 1, a 2,...} We assume a general friction term ν j (a 1, a 2,...) = ν j (a) and assume a Gaussian noise ξ j (t) where The Langevin equation becomes ξ j (t) = ξ i (t 2 )ξ j (t 1 ) = g ij δ(t 2 t 1) d a j(t) = ν j (a) + ξ j (t)
5 7.3 General Fokker-Planck equation 93 or in vector form d a(t) = ν(a) + ξ(t) We ask for the probability distribution Again from conservation of probability P (a, t) = ϱ(a, t) ξ ϱ(a, t) + t a da ϱ(a, t) = Usin the Langevin equation to solve for da/ we find where ϱ(a, t) = (ν(a)ϱ(a, t)) t a a (ξ(t)ϱ(a, t)) = L + L 1 (t) ϱ(a, t) Following the steps as above we find L = ( ) a ν(a) + ν(a) a L 1 (t) = ξ(t) a P (a, t) = t a (ν(a)p (a, t)) Here g is a tensor with elements g ij. Example a g P (a, t) (7.6) a Our previous result can be obtained as a special case of (7.6). The Langevin equations are where Then a = ξ(t) = dx dv = v = γ m v + 1 m F (x) + 1 m ξ(t) ξ(t 2 )ξ(t 1) = 2γk B T δ(t 2 t 1 ) ( x y ( 1 m ξ(t) The Fokker-Planck equation becomes P (x, v, t) = t x vp (x, v, t) v ) (, ν(a) = ), g = ) v γ m v + 1 m ( F ) (x) 2γk B T m 2 ( γm v + 1m F (x) ) P (x, v, t) + γk BT 2 m 2 P (x, v, t) v2
6 94 Chapter 7 Brownian Motion: Fokker-Planck Equation The equilibrium distribution /tp (x, v, t) = is the solution to x vp (x, v, t) ( γm v v + 1m ) F (x) P (x, v, t) + γk BT 2 m 2 P (x, v, t) = v2 The Hamiltonian of the Brownian particle is H = 1 2 mv2 + V (x) i.e we can write the equilibrium condition as H P (x, v, t) + 1 H P (x, v, t) + γ x v m v x m v Assuming that we find x H v f(h) + 1 m v This equation is satisfied if with the solution P = f(h) H v P (x, v, t) + k BT m P (x, v, t) = v H x f(h) + γ H H m v v f(h) + k BT H m v f (H) = 2 H xv = f(h) + k B T f (H) = f(h) = 1 Z e βh, β = 1 k B T The Fokker-Planck equation is therefore consistent with a Boltzmann equilibrium distribution. 7.4 Averages and adjoint operators Sometimes we want the full solution of a Fokker-Planck equation, but sometimes we are interested only in certain averages. These can be found by two distinct but equivalent procedures analogous to the Schrödinger- Heisenberg duality in quantum mechanics. First we can follow the evolution of some initial state P (a, ) by solving the Fokker- Planck equation P (a, t) = LP (a, t) t where the operator L is given by L = a v(a) a D a with v a streaming term and D a diffusion tensor. The Fokker-Planck equation has the formal operator solution P (a, t) = e Lt P (a, )
7 7.5 Green s function 95 This can be used to obtain the average of any dynamical property A(a) A(t) = daa(a)p (a, t) = daa(a)e Lt P (a, ) (7.7) We introduce the adjoint operator L defined by ( ) daφ (a)lψ(a) = daψ (a)l φ(a) or φ Lψ = L φ ψ. By partial integration we find L = v(a) a a D a The average in (7.7) can be obtained by partial integration A(t) = = daa(a)e Lt ( t) n P (a, ) = daa(a)l n P (a, ) n! n= ( t) n da (L ) n A(a) P (a, ) = da e L t A(a) P (a, ) n! n= This defines the time-dependent variable A(a, t) = e L t A(a) The equation of motion for the dynamical variable A(a, t) becomes ( t A(a, t) = L A(a, t) = v(a) a ) a D A(a, t) (7.8) a 7.5 Green s function A formal solution of the Fokker-Planck equation uses the operator e Lt. more explicit solution uses the Green s function G(a, t a ) P (a, t) = da G(a, t)p (a, ) A somewhat The Green s function is the conditional probability to find the systems in state a at time t when is started at a at time t =. It satisfies the equation When the streaming term v(a) is linear in a t G(a, t a ) = LG(a, t a ) G(a, a ) = δ(a a ) v(a) = C a
8 96 Chapter 7 Brownian Motion: Fokker-Planck Equation and the diffusion tensor D is constant independent of a, the Green function can be found explicitly by the Fourier transform Ĝ(k, t a ) = dae ik a G(a, t a ) Then tĝ(k, t a ) = k C Ĝ(k, t a ) k D k kĝ(k, t a ) which follows since a derivatives of a gives a factor ik and multiplication by a gives a derivative with repsect ot. Dividing by Ĝ(k, t a ) gives If we make the Gaussian ansatz t ln Ĝ(k, t a ) = k C k ln Ĝ(k, t a ) k D k ln Ĝ(k, t a ) = im(t) k 1 2 k S(t) k with the unknown vector m(t) and tensor S(t) we find i d m(t) k 1 2 k d S(t) k = ik C m(t) 1 2 k C S(t) k 1 2 k S(t) CT k k D k Identifying equal powers of k this gives d m(t) = C m(t) d S(t) = C S(t) + S(t) CT + 2D (7.9) From the initial value ln Ĝ(k, a ) = ik a we find The solution to (7.9) can then be written m() = a, S() = m(t) = e Ct a (7.1) S(t) = 2 dse C(t s)t D e CT (t s)t (7.11) which one can easily verify by substitution. These quantities have a simple interpretation in terms of averages and mean squared fluctuations since ( ) Ĝ(k, t a ) = daiag(a, t a ) = i a = im(t) k and k= ( ) Ĝ(k, t a ) = k k k= daaag(a, t a ) = aa = m(t)m(t) S(t)
9 7.6 Examples 97 and S(t) = aa a a = a a a a Since G(k, t a ) is a Gaussian the inverse transform also give a Gaussian and follows the same lines as we did before in connection with the multivariate Gaussian distribution and we find 1 1/2 ( G(a, t a ) = 2πdetS(t) = exp 1 ) 2 a m(t) S 1 (t) a m(t) (7.12) 7.6 Examples Free Brownian particle For a free Brownian particle we have the Langevin equation Then with vector notation m dv = γv + ξ(t) ae v, v(a) = γ m ve v, Ξ = 1 m ξ(t)e v, g = 2γk BT m 2 The general Fokker-Planck equation (7.6) P (a, t) = t a (ν(a)p (a, t)) I a D P (a, t) a becomes t P (v, t) = ( γ ) v m vp (v, t) + γk BT 2 m 2 P (v, t) (7.13) v2 This is a linear equation and we can directly obtain the solution for the Green function from (7.12). The moments are obtained from (7.11) with C = γ/mi. Therefore m(t) = e Ct v e v = e t/τ B v e v = m(t)e v, τ B = m γ Also S(t) = 2 dse C(t s) D e CT (t s)t = 2 dse 2(t s)/τ B γk BT m 2 I = k BT m I 1 e 2t/τ B = S(t)I The Green function is then 1 1/2 ) ( (v m(t))2 βm G(v, t v ) = exp ( = 2πS(t) 2S(t) 2π 1 e 2t/τ exp βm B 2 We see that G(v, t v ) i.e a Maxwell-Boltzmann distribution as expected. ( ) βm 1/2 e βmv2 /2, t 2π ( v v e t/τ B) 2 ) 1 e 2t/τ B
10 98 Chapter 7 Brownian Motion: Fokker-Planck Equation Diffusion controlled reactions The rate of chemical reactions in condensed phases is often determined by the rate at which the reacting species get close enough together to react. Once they get within some critical distance of each other, the making and breaking of chemical bonds takes place rapidly. The process by which the reactive particles reach that critical distance is the brownian motion of the particles themselves. Such reactions are called diffusion controlled reactions since the diffusive Brownian motion of the species is the rate determining step. There are two reacting species A and B in solution. At first we consider B to be stationary and species A to diffuse with diffusion constant D A. When an A particle comes within a distance R of a B particle the particles react and effectively disappear from the problem. For simplicity the particles are assumed to be spherical. We also assume that there are no forces between A and B when they are farther apart than R. Let the concentration of A and B particles be n A and n B respectively. The rate at which A particles disappear is then n B times the flux of A particles across the sphere of radius R surrounding a B particle. If j be the current density of A particles we have rate = n B j ds = 4πR 2 n B j r (R, t) where j r is the radial component of j. The integration is over the surface of a sphere of radius R about the stationary B particle, and the current j is supposd not to vary on the surface S. The current is given by Ficks law j = D A n A and the continuity equation for n A leads to the diffusion equation n A t ( = D A 2 1 n A = D A r 2 r 2 n ) A r r In the last step we have expressed the Laplacian in terms of spherical coordinates and used the spherical symmetry in the problem. The initial and boundary conditions are n A (r, t = ) = n A n A (r, t) n A, r n A (R, t) =, r = R The second condition requires that the concentration very far from the sink (B particle) is not perturbed by the sink. The last condition insures that the sphere of radius R is a sink. At r = R the A particle reacts and it is then no longer an A particle. We define a new function u(r, t) n A (r, t) = r then r 2 n A r = u(r, t) + r u(r, t) r
11 7.6 Examples 99 and so u(r, t) satisfies the simple diffusion equation and the boundary conditions are t u(r, t) = D A u(r, t) (7.14) r2 2 u(r, t) n Ar r u(r, t) = The solution to (7.15) with the initial condition U(r, ) = n Ar, r > R is u(r, t) = R dr G(r, t r, )n Ar where the Green functiong satisfies (7.14) with a δ(r r ) initial condition and boundary condition G(R, t r ) = For a homogeneous medium the diffusion equation has the solution G (r, t r, ) = ( ) 1 1/2 e (r r )2 /2Dt 4πDt To satisfy the boundary condition we introduce a mirror source in = s and choose s so that the boundary condition is satisfied, i.e. G(r, t r, ) = G (r, t r, ) G (r, t s, ) Then G(R, t r, ) = provided { (R r ) 2 = (R s) 2, s = r 2R r Therefore the Green-function which satisfies the boundary condition at r = R is G(r, t r, ) = ( ) 1 1/2 e (r r ) 2 /2Dt e (r+r 2R) 2 /2Dt 4πDt The solution for the density of A particles n A (r, t) in the region r > R therefore becomes n A (r, t) = n A r Then by a change of variables n A (r, t) = = n A n A 2πr 1 2R r ( ) 1 1/2 dr e (r r ) 2 /2Dt e (r+r 2R) 2 /2Dt r 4πDt R ( r R 2Dt dξe ξ 2 r ) 2Dtξ 1 2π r R dξe ξ2 = n A 2Dt r R ( ) dξe ξ2 2Dtξ + 2R r 2Dt 1 R r + R 2 r R 2Dt dξe ξ 2 r 2π
12 1 Chapter 7 Brownian Motion: Fokker-Planck Equation The radial component of the current for particle A then becomes R j r (r, t) = D A r n A(r, t) = D A n A r 2 + R 1 e (r R)2 /2Dt R 2 r πdt r 2 2π r R 2Dt dξe ξ 2 and so This gives the rate j r (R, t) D A n A R + 1 πdt rate = D A n A n B 4πR R + 1 = k AB n A n B πdt with the rate constant k AB = D A 4πR R + 1 πdt The time dependent correction is of little interest experimentally. In reality the B particles move around just as do the A particles. We can take account of this motion by replacing D A with D A + D B. Then wth the Stokes-Einstein relation D i = k BT 6πηRi we find k AB = 4πR(D A + D B ) = 2k BT 3η (R A + R B ) 2 R A R B
The Kramers problem and first passage times.
Chapter 8 The Kramers problem and first passage times. The Kramers problem is to find the rate at which a Brownian particle escapes from a potential well over a potential barrier. One method of attack
More informationFokker-Planck Equation with Detailed Balance
Appendix E Fokker-Planck Equation with Detailed Balance A stochastic process is simply a function of two variables, one is the time, the other is a stochastic variable X, defined by specifying: a: the
More informationA path integral approach to the Langevin equation
A path integral approach to the Langevin equation - Ashok Das Reference: A path integral approach to the Langevin equation, A. Das, S. Panda and J. R. L. Santos, arxiv:1411.0256 (to be published in Int.
More informationLangevin Methods. Burkhard Dünweg Max Planck Institute for Polymer Research Ackermannweg 10 D Mainz Germany
Langevin Methods Burkhard Dünweg Max Planck Institute for Polymer Research Ackermannweg 1 D 55128 Mainz Germany Motivation Original idea: Fast and slow degrees of freedom Example: Brownian motion Replace
More informationStatistical Mechanics of Active Matter
Statistical Mechanics of Active Matter Umberto Marini Bettolo Marconi University of Camerino, Italy Naples, 24 May,2017 Umberto Marini Bettolo Marconi (2017) Statistical Mechanics of Active Matter 2017
More informationLANGEVIN THEORY OF BROWNIAN MOTION. Contents. 1 Langevin theory. 1 Langevin theory 1. 2 The Ornstein-Uhlenbeck process 8
Contents LANGEVIN THEORY OF BROWNIAN MOTION 1 Langevin theory 1 2 The Ornstein-Uhlenbeck process 8 1 Langevin theory Einstein (as well as Smoluchowski) was well aware that the theory of Brownian motion
More informationStochastic Particle Methods for Rarefied Gases
CCES Seminar WS 2/3 Stochastic Particle Methods for Rarefied Gases Julian Köllermeier RWTH Aachen University Supervisor: Prof. Dr. Manuel Torrilhon Center for Computational Engineering Science Mathematics
More informationIntroduction to nonequilibrium physics
Introduction to nonequilibrium physics Jae Dong Noh December 18, 2016 Preface This is a note for the lecture given in the 2016 KIAS-SNU Physics Winter Camp which is held at KIAS in December 17 23, 2016.
More informationThe dynamics of small particles whose size is roughly 1 µmt or. smaller, in a fluid at room temperature, is extremely erratic, and is
1 I. BROWNIAN MOTION The dynamics of small particles whose size is roughly 1 µmt or smaller, in a fluid at room temperature, is extremely erratic, and is called Brownian motion. The velocity of such particles
More informationmacroscopic view (phenomenological) microscopic view (atomistic) computing reaction rate rate of reactions experiments thermodynamics
Rate Theory (overview) macroscopic view (phenomenological) rate of reactions experiments thermodynamics Van t Hoff & Arrhenius equation microscopic view (atomistic) statistical mechanics transition state
More informationBrownian 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 informationOpen Quantum Systems and Markov Processes II
Open Quantum Systems and Markov Processes II Theory of Quantum Optics (QIC 895) Sascha Agne sascha.agne@uwaterloo.ca July 20, 2015 Outline 1 1. Introduction to open quantum systems and master equations
More information08. Brownian Motion. University of Rhode Island. Gerhard Müller University of Rhode Island,
University of Rhode Island DigitalCommons@URI Nonequilibrium Statistical Physics Physics Course Materials 1-19-215 8. Brownian Motion Gerhard Müller University of Rhode Island, gmuller@uri.edu Follow this
More informationSTOCHASTIC PROCESSES IN PHYSICS AND CHEMISTRY
STOCHASTIC PROCESSES IN PHYSICS AND CHEMISTRY Third edition N.G. VAN KAMPEN Institute for Theoretical Physics of the University at Utrecht ELSEVIER Amsterdam Boston Heidelberg London New York Oxford Paris
More informationStochastic continuity equation and related processes
Stochastic continuity equation and related processes Gabriele Bassi c Armando Bazzani a Helmut Mais b Giorgio Turchetti a a Dept. of Physics Univ. of Bologna, INFN Sezione di Bologna, ITALY b DESY, Hamburg,
More informationBROWNIAN DYNAMICS SIMULATIONS WITH HYDRODYNAMICS. Abstract
BROWNIAN DYNAMICS SIMULATIONS WITH HYDRODYNAMICS Juan J. Cerdà 1 1 Institut für Computerphysik, Pfaffenwaldring 27, Universität Stuttgart, 70569 Stuttgart, Germany. (Dated: July 21, 2009) Abstract - 1
More informationThis 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 informationSmoluchowski Diffusion Equation
Chapter 4 Smoluchowski Diffusion Equation Contents 4. Derivation of the Smoluchoswki Diffusion Equation for Potential Fields 64 4.2 One-DimensionalDiffusoninaLinearPotential... 67 4.2. Diffusion in an
More information16. Working with the Langevin and Fokker-Planck equations
16. Working with the Langevin and Fokker-Planck equations In the preceding Lecture, we have shown that given a Langevin equation (LE), it is possible to write down an equivalent Fokker-Planck equation
More informationActive Matter Lectures for the 2011 ICTP School on Mathematics and Physics of Soft and Biological Matter Lecture 3: Hydrodynamics of SP Hard Rods
Active Matter Lectures for the 2011 ICTP School on Mathematics and Physics of Soft and Biological Matter Lecture 3: of SP Hard Rods M. Cristina Marchetti Syracuse University Baskaran & MCM, PRE 77 (2008);
More informationExample 4.1 Let X be a random variable and f(t) a given function of time. Then. Y (t) = f(t)x. Y (t) = X sin(ωt + δ)
Chapter 4 Stochastic Processes 4. Definition In the previous chapter we studied random variables as functions on a sample space X(ω), ω Ω, without regard to how these might depend on parameters. We now
More informationPath integrals for classical Markov processes
Path integrals for classical Markov processes Hugo Touchette National Institute for Theoretical Physics (NITheP) Stellenbosch, South Africa Chris Engelbrecht Summer School on Non-Linear Phenomena in Field
More informationMesoscale Simulation Methods. Ronojoy Adhikari The Institute of Mathematical Sciences Chennai
Mesoscale Simulation Methods Ronojoy Adhikari The Institute of Mathematical Sciences Chennai Outline What is mesoscale? Mesoscale statics and dynamics through coarse-graining. Coarse-grained equations
More informationCHAPTER V. Brownian motion. V.1 Langevin dynamics
CHAPTER V Brownian motion In this chapter, we study the very general paradigm provided by Brownian motion. Originally, this motion is that a heavy particle, called Brownian particle, immersed in a fluid
More informationDiffusion 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 informationLangevin Dynamics of a Single Particle
Overview Langevin Dynamics of a Single Particle We consider a spherical particle of radius r immersed in a viscous fluid and suppose that the dynamics of the particle depend on two forces: A drag force
More informationHandbook of Stochastic Methods
C. W. Gardiner Handbook of Stochastic Methods for Physics, Chemistry and the Natural Sciences Third Edition With 30 Figures Springer Contents 1. A Historical Introduction 1 1.1 Motivation I 1.2 Some Historical
More informationNon-equilibrium phenomena and fluctuation relations
Non-equilibrium phenomena and fluctuation relations Lamberto Rondoni Politecnico di Torino Beijing 16 March 2012 http://www.rarenoise.lnl.infn.it/ Outline 1 Background: Local Thermodyamic Equilibrium 2
More informationFluid-Particles Interaction Models Asymptotics, Theory and Numerics I
Fluid-Particles Interaction Models Asymptotics, Theory and Numerics I J. A. Carrillo collaborators: T. Goudon (Lille), P. Lafitte (Lille) and F. Vecil (UAB) (CPDE 2005), (JCP, 2008), (JSC, 2008) ICREA
More informationDiffusion in the cell
Diffusion in the cell Single particle (random walk) Microscopic view Macroscopic view Measuring diffusion Diffusion occurs via Brownian motion (passive) Ex.: D = 100 μm 2 /s for typical protein in water
More informationThe Truth about diffusion (in liquids)
The Truth about diffusion (in liquids) Aleksandar Donev Courant Institute, New York University & Eric Vanden-Eijnden, Courant In honor of Berni Julian Alder LLNL, August 20th 2015 A. Donev (CIMS) Diffusion
More informationExercise Set 4. D s n ds + + V. s dv = V. After using Stokes theorem, the surface integral becomes
Exercise Set Exercise - (a) Let s consider a test volum in the pellet. The substract enters the pellet by diffusion and some is created and disappears due to the chemical reaction. The two contribute to
More informationSolutions of a Focker-Planck Equation
Solutions of a Focker-Planck Equation W. C. Troy March 31, 006 1. Overview.. The Langevin Equation for Brownian Motion. 3. The Focker-Planck equation. 4. The steady state equation. 1 Overview. In these
More informationFrom a Mesoscopic to a Macroscopic Description of Fluid-Particle Interaction
From a Mesoscopic to a Macroscopic Description of Fluid-Particle Interaction Carnegie Mellon University Center for Nonlinear Analysis Working Group, October 2016 Outline 1 Physical Framework 2 3 Free Energy
More informationStochastic equations for thermodynamics
J. Chem. Soc., Faraday Trans. 93 (1997) 1751-1753 [arxiv 1503.09171] Stochastic equations for thermodynamics Roumen Tsekov Department of Physical Chemistry, University of Sofia, 1164 Sofia, ulgaria The
More informationthe renormalization group (RG) idea
the renormalization group (RG) idea Block Spin Partition function Z =Tr s e H. block spin transformation (majority rule) T (s, if s i ; s,...,s 9 )= s i > 0; 0, otherwise. b Block Spin (block-)transformed
More informationSupporting Information
Supporting Information A: Calculation of radial distribution functions To get an effective propagator in one dimension, we first transform 1) into spherical coordinates: x a = ρ sin θ cos φ, y = ρ sin
More informationMathematical Methods for Neurosciences. ENS - Master MVA Paris 6 - Master Maths-Bio ( )
Mathematical Methods for Neurosciences. ENS - Master MVA Paris 6 - Master Maths-Bio (2014-2015) Etienne Tanré - Olivier Faugeras INRIA - Team Tosca November 26th, 2014 E. Tanré (INRIA - Team Tosca) Mathematical
More informationMajor Concepts Kramers Turnover
Major Concepts Kramers Turnover Low/Weak -Friction Limit (Energy Diffusion) Intermediate Regime bounded by TST High/Strong -Friction Limit Smoluchovski (Spatial Diffusion) Fokker-Planck Equation Probability
More informationFluid Equations for Rarefied Gases
1 Fluid Equations for Rarefied Gases Jean-Luc Thiffeault Department of Applied Physics and Applied Mathematics Columbia University http://plasma.ap.columbia.edu/~jeanluc 23 March 2001 with E. A. Spiegel
More informationLecture 12: Detailed balance and Eigenfunction methods
Miranda Holmes-Cerfon Applied Stochastic Analysis, Spring 2015 Lecture 12: Detailed balance and Eigenfunction methods Readings Recommended: Pavliotis [2014] 4.5-4.7 (eigenfunction methods and reversibility),
More information1 Elementary probability
1 Elementary probability Problem 1.1 (*) A coin is thrown several times. Find the probability, that at the n-th experiment: (a) Head appears for the first time (b) Head and Tail have appeared equal number
More informationConfined Brownian Motion:
Master Thesis Department of Physics/Stockholm University and NORDITA Confined Brownian Motion: Fick-Jacobs Equation & Stochastic Thermodynamics Author: Christoph Streißnig Supervisor: Ralf Eichhorn Co-Supervisor:
More informationG : Statistical Mechanics
G25.2651: Statistical Mechanics Notes for Lecture 15 Consider Hamilton s equations in the form I. CLASSICAL LINEAR RESPONSE THEORY q i = H p i ṗ i = H q i We noted early in the course that an ensemble
More information2. Molecules in Motion
2. Molecules in Motion Kinetic Theory of Gases (microscopic viewpoint) assumptions (1) particles of mass m and diameter d; ceaseless random motion (2) dilute gas: d λ, λ = mean free path = average distance
More informationA Detailed Look at a Discrete Randomw Walk with Spatially Dependent Moments and Its Continuum Limit
A Detailed Look at a Discrete Randomw Walk with Spatially Dependent Moments and Its Continuum Limit David Vener Department of Mathematics, MIT May 5, 3 Introduction In 8.366, we discussed the relationship
More informationCosmological Issues. Consider the stress tensor of a fluid in the local orthonormal frame where the metric is η ab
Cosmological Issues Radiation dominated Universe Consider the stress tensor of a fluid in the local orthonormal frame where the metric is η ab ρ 0 0 0 T ab = 0 p 0 0 0 0 p 0 () 0 0 0 p We do not often
More informationLANGEVIN EQUATION AND THERMODYNAMICS
LANGEVIN EQUATION AND THERMODYNAMICS RELATING STOCHASTIC DYNAMICS WITH THERMODYNAMIC LAWS November 10, 2017 1 / 20 MOTIVATION There are at least three levels of description of classical dynamics: thermodynamic,
More informationHandbook of Stochastic Methods
Springer Series in Synergetics 13 Handbook of Stochastic Methods for Physics, Chemistry and the Natural Sciences von Crispin W Gardiner Neuausgabe Handbook of Stochastic Methods Gardiner schnell und portofrei
More informationFluid Equations for Rarefied Gases
1 Fluid Equations for Rarefied Gases Jean-Luc Thiffeault Department of Applied Physics and Applied Mathematics Columbia University http://plasma.ap.columbia.edu/~jeanluc 21 May 2001 with E. A. Spiegel
More informationPath Intergal. 1 Introduction. 2 Derivation From Schrödinger Equation. Shoichi Midorikawa
Path Intergal Shoichi Midorikawa 1 Introduction The Feynman path integral1 is one of the formalism to solve the Schrödinger equation. However this approach is not peculiar to quantum mechanics, and M.
More informationAnomalous 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 information1 Phase Spaces and the Liouville Equation
Phase Spaces and the Liouville Equation emphasize the change of language from deterministic to probablistic description. Under the dynamics: ½ m vi = F i ẋ i = v i with initial data given. What is the
More informationActive and Driven Soft Matter Lecture 3: Self-Propelled Hard Rods
A Tutorial: From Langevin equation to Active and Driven Soft Matter Lecture 3: Self-Propelled Hard Rods M. Cristina Marchetti Syracuse University Boulder School 2009 A Tutorial: From Langevin equation
More informationQuantifying Intermittent Transport in Cell Cytoplasm
Quantifying Intermittent Transport in Cell Cytoplasm Ecole Normale Supérieure, Mathematics and Biology Department. Paris, France. May 19 th 2009 Cellular Transport Introduction Cellular Transport Intermittent
More informationI. QUANTUM MONTE CARLO METHODS: INTRODUCTION AND BASICS
I. QUANTUM MONTE CARLO METHODS: INTRODUCTION AND BASICS Markus Holzmann LPMMC, UJF, Grenoble, and LPTMC, UPMC, Paris markus@lptl.jussieu.fr http://www.lptl.jussieu.fr/users/markus (Dated: January 24, 2012)
More informationGreen-Kubo formula for heat conduction in open systems
Green-Kubo formula for heat conduction in open systems Abhishek Dhar Raman Research Institute, Bangalore, India. Collaborators: Anupam Kundu(Raman Research Institute) Onuttom Narayan(UC Santa Cruz) Green-Kubo
More informationFor one such collision, a new particle, "k", is formed of volume v k = v i + v j. The rate of formation of "k" particles is, N ij
Coagulation Theory and The Smoluchowski Equation: Coagulation is defined as growth of particles by collisions among particles. It is usually associated with dense, 3-d growth, in contrast to aggregation
More information2. Basic assumptions for stellar atmospheres
. Basic assumptions for stellar atmospheres 1. geometry, stationarity. conservation of momentum, mass 3. conservation of energy 4. Local Thermodynamic Equilibrium 1 1. Geometry Stars as gaseous spheres
More information221A Miscellaneous Notes Continuity Equation
221A Miscellaneous Notes Continuity Equation 1 The Continuity Equation As I received questions about the midterm problems, I realized that some of you have a conceptual gap about the continuity equation.
More informationLinear Diffusion and Image Processing. Outline
Outline Linear Diffusion and Image Processing Fourier Transform Convolution Image Restoration: Linear Filtering Diffusion Processes for Noise Filtering linear scale space theory Gauss-Laplace pyramid for
More information3. Stochastic Processes
3. Stochastic Processes We learn in kindergarten about the phenomenon of Brownian motion, the random jittery movement that a particle suffers when it is placed in a liquid. Famously, it is caused by the
More informationNon equilibrium thermodynamic transformations. Giovanni Jona-Lasinio
Non equilibrium thermodynamic transformations Giovanni Jona-Lasinio Kyoto, July 29, 2013 1. PRELIMINARIES 2. RARE FLUCTUATIONS 3. THERMODYNAMIC TRANSFORMATIONS 1. PRELIMINARIES Over the last ten years,
More informationSession 1: Probability and Markov chains
Session 1: Probability and Markov chains 1. Probability distributions and densities. 2. Relevant distributions. 3. Change of variable. 4. Stochastic processes. 5. The Markov property. 6. Markov finite
More informationBAE 820 Physical Principles of Environmental Systems
BAE 820 Physical Principles of Environmental Systems Estimation of diffusion Coefficient Dr. Zifei Liu Diffusion mass transfer Diffusion mass transfer refers to mass in transit due to a species concentration
More informationCosmological Issues. Consider the stress tensor of a fluid in the local orthonormal frame where the metric is η ab
Cosmological Issues 1 Radiation dominated Universe Consider the stress tensor of a fluid in the local orthonormal frame where the metric is η ab ρ 0 0 0 T ab = 0 p 0 0 0 0 p 0 (1) 0 0 0 p We do not often
More informationThe First Principle Calculation of Green Kubo Formula with the Two-Time Ensemble Technique
Commun. Theor. Phys. (Beijing, China 35 (2 pp. 42 46 c International Academic Publishers Vol. 35, No. 4, April 5, 2 The First Principle Calculation of Green Kubo Formula with the Two-Time Ensemble Technique
More information14. Energy transport.
Phys780: Plasma Physics Lecture 14. Energy transport. 1 14. Energy transport. Chapman-Enskog theory. ([8], p.51-75) We derive macroscopic properties of plasma by calculating moments of the kinetic equation
More informationBrownian Motion and Langevin Equations
1 Brownian Motion and Langevin Equations 1.1 Langevin Equation and the Fluctuation- Dissipation Theorem The theory of Brownian motion is perhaps the simplest approximate way to treat the dynamics of nonequilibrium
More informationNon-equilibrium phase transitions
Non-equilibrium phase transitions An Introduction Lecture III Haye Hinrichsen University of Würzburg, Germany March 2006 Third Lecture: Outline 1 Directed Percolation Scaling Theory Langevin Equation 2
More informationSpin-Boson Model. A simple Open Quantum System. M. Miller F. Tschirsich. Quantum Mechanics on Macroscopic Scales Theory of Condensed Matter July 2012
Spin-Boson Model A simple Open Quantum System M. Miller F. Tschirsich Quantum Mechanics on Macroscopic Scales Theory of Condensed Matter July 2012 Outline 1 Bloch-Equations 2 Classical Dissipations 3 Spin-Boson
More informationSolutions: Homework 5
Ex. 5.1: Capacitor Solutions: Homework 5 (a) Consider a parallel plate capacitor with large circular plates, radius a, a distance d apart, with a d. Choose cylindrical coordinates (r,φ,z) and let the z
More information1 Introduction. 2 Diffusion equation and central limit theorem. The content of these notes is also covered by chapter 3 section B of [1].
1 Introduction The content of these notes is also covered by chapter 3 section B of [1]. Diffusion equation and central limit theorem Consider a sequence {ξ i } i=1 i.i.d. ξ i = d ξ with ξ : Ω { Dx, 0,
More informationTemperature Fluctuations and Entropy Formulas
Temperature Fluctuations and Entropy Formulas Does it or does not? T.S. Biró G.G. Barnaföldi P. Ván MTA Heavy Ion Research Group Research Centre for Physics, Budapest February 1, 2014 J.Uffink, J.van Lith:
More informationThe Polymers Tug Back
Tugging at Polymers in Turbulent Flow The Polymers Tug Back Jean-Luc Thiffeault http://plasma.ap.columbia.edu/ jeanluc Department of Applied Physics and Applied Mathematics Columbia University Tugging
More informationSolution. For one question the mean grade is ḡ 1 = 10p = 8 and the standard deviation is 1 = g
Exam 23 -. To see how much of the exam grade spread is pure luck, consider the exam with ten identically difficult problems of ten points each. The exam is multiple-choice that is the correct choice of
More informationLecture 5: Kinetic theory of fluids
Lecture 5: Kinetic theory of fluids September 21, 2015 1 Goal 2 From atoms to probabilities Fluid dynamics descrines fluids as continnum media (fields); however under conditions of strong inhomogeneities
More informationTable of Contents [ntc]
Table of Contents [ntc] 1. Introduction: Contents and Maps Table of contents [ntc] Equilibrium thermodynamics overview [nln6] Thermal equilibrium and nonequilibrium [nln1] Levels of description in statistical
More informationI. Dissipative Dynamics
I. Dissipative Dynamics I.A Brownian Motion of a Particle Observations under a microscope indicate that a dust particle in a liquid drop undergoes a random jittery motion. This is because of the random
More informationGeometric Dyson Brownian motion and May Wigner stability
Geometric Dyson Brownian motion and May Wigner stability Jesper R. Ipsen University of Melbourne Summer school: Randomness in Physics and Mathematics 2 Bielefeld 2016 joint work with Henning Schomerus
More informationPAPER 310 COSMOLOGY. Attempt no more than THREE questions. There are FOUR questions in total. The questions carry equal weight.
MATHEMATICAL TRIPOS Part III Wednesday, 1 June, 2016 9:00 am to 12:00 pm PAPER 310 COSMOLOGY Attempt no more than THREE questions. There are FOUR questions in total. The questions carry equal weight. STATIONERY
More informationPhysics 9 Spring 2011 Midterm 1 Solutions
Physics 9 Spring 2011 Midterm 1 s For the midterm, you may use one sheet of notes with whatever you want to put on it, front and back. Please sit every other seat, and please don t cheat! If something
More informationLecture 12: Detailed balance and Eigenfunction methods
Lecture 12: Detailed balance and Eigenfunction methods Readings Recommended: Pavliotis [2014] 4.5-4.7 (eigenfunction methods and reversibility), 4.2-4.4 (explicit examples of eigenfunction methods) Gardiner
More informationINTRODUCTION TO о JLXJLA Из А lv-/xvj_y JrJrl Y üv_>l3 Second Edition
INTRODUCTION TO о JLXJLA Из А lv-/xvj_y JrJrl Y üv_>l3 Second Edition Kerson Huang CRC Press Taylor & Francis Croup Boca Raton London New York CRC Press is an imprint of the Taylor & Francis Group an Informa
More informationThe Smoluchowski-Kramers Approximation: What model describes a Brownian particle?
The Smoluchowski-Kramers Approximation: What model describes a Brownian particle? Scott Hottovy shottovy@math.arizona.edu University of Arizona Applied Mathematics October 7, 2011 Brown observes a particle
More information5 Applying the Fokker-Planck equation
5 Applying the Fokker-Planck equation We begin with one-dimensional examples, keeping g = constant. Recall: the FPE for the Langevin equation with η(t 1 )η(t ) = κδ(t 1 t ) is = f(x) + g(x)η(t) t = x [f(x)p
More information12. MHD Approximation.
Phys780: Plasma Physics Lecture 12. MHD approximation. 1 12. MHD Approximation. ([3], p. 169-183) The kinetic equation for the distribution function f( v, r, t) provides the most complete and universal
More informationF r (t) = 0, (4.2) F (t) = r= r a (t)
Chapter 4 Stochastic Equations 4.1 Langevin equation To explain the idea of introduction of the Langevin equation, let us consider the concrete example taken from surface physics: motion of an atom (adatom)
More informationNew Physical Principle for Monte-Carlo simulations
EJTP 6, No. 21 (2009) 9 20 Electronic Journal of Theoretical Physics New Physical Principle for Monte-Carlo simulations Michail Zak Jet Propulsion Laboratory California Institute of Technology, Advance
More information2. Basic Assumptions for Stellar Atmospheres
2. Basic Assumptions for Stellar Atmospheres 1. geometry, stationarity 2. conservation of momentum, mass 3. conservation of energy 4. Local Thermodynamic Equilibrium 1 1. Geometry Stars as gaseous spheres!
More informationHomogenization with stochastic differential equations
Homogenization with stochastic differential equations Scott Hottovy shottovy@math.arizona.edu University of Arizona Program in Applied Mathematics October 12, 2011 Modeling with SDE Use SDE to model system
More informationThe Brownian Oscillator
Chapter The Brownian Oscillator Contents.One-DimensionalDiffusioninaHarmonicPotential...38 The one-dimensional Smoluchowski equation, in case that a stationary flux-free equilibrium state p o (x) exists,
More informationVIII.B Equilibrium Dynamics of a Field
VIII.B Equilibrium Dynamics of a Field The next step is to generalize the Langevin formalism to a collection of degrees of freedom, most conveniently described by a continuous field. Let us consider the
More information2. Basic assumptions for stellar atmospheres
. Basic assumptions for stellar atmospheres 1. geometry, stationarity. conservation of momentum, mass 3. conservation of energy 4. Local Thermodynamic Equilibrium 1 1. Geometry Stars as gaseous spheres
More informationIntroduction to nonequilibrium physics
Introduction to nonequilibrium physics Jae Dong Noh December 19, 2016 Preface This is a note for the lecture given in the 2016 KIAS-SNU Physics Winter Camp which is held at KIAS in December 17 23, 2016.
More informationJ10M.1 - Rod on a Rail (M93M.2)
Part I - Mechanics J10M.1 - Rod on a Rail (M93M.2) J10M.1 - Rod on a Rail (M93M.2) s α l θ g z x A uniform rod of length l and mass m moves in the x-z plane. One end of the rod is suspended from a straight
More informationELEMENTS OF PROBABILITY THEORY
ELEMENTS OF PROBABILITY THEORY Elements of Probability Theory A collection of subsets of a set Ω is called a σ algebra if it contains Ω and is closed under the operations of taking complements and countable
More informationVector and Tensor Calculus
Appendices 58 A Vector and Tensor Calculus In relativistic theory one often encounters vector and tensor expressions in both three- and four-dimensional form. The most important of these expressions are
More informationLINEAR RESPONSE THEORY
MIT Department of Chemistry 5.74, Spring 5: Introductory Quantum Mechanics II Instructor: Professor Andrei Tokmakoff p. 8 LINEAR RESPONSE THEORY We have statistically described the time-dependent behavior
More informationGaussian processes for inference in stochastic differential equations
Gaussian processes for inference in stochastic differential equations Manfred Opper, AI group, TU Berlin November 6, 2017 Manfred Opper, AI group, TU Berlin (TU Berlin) inference in SDE November 6, 2017
More information