990 likes | 1.35k Views
Lecture 6: Langevin equations. Outline: linear/nonlinear, additive and multiplicative noise soluble linear example w/ additive noise: Ornstein-Uhlenbeck process general 1-d nonlinear equation with multiplicative noise relation to Fokker-Planck equation
E N D
Lecture 6: Langevin equations • Outline: • linear/nonlinear, additive and multiplicative noise • soluble linear example w/ additive noise: Ornstein-Uhlenbeck process • general 1-d nonlinear equation with multiplicative noise • relation to Fokker-Planck equation • Ito formulation, relation between Ito & Stratonovich approaches
Stochastic differential equations Differential equations which contain (“are driven by”) random functions
Stochastic differential equations Differential equations which contain (“are driven by”) random functions Langevin equations: the random function is Gaussian white noiseξ(t):
Stochastic differential equations Differential equations which contain (“are driven by”) random functions Langevin equations: the random function is Gaussian white noiseξ(t):
Stochastic differential equations Differential equations which contain (“are driven by”) random functions Langevin equations: the random function is Gaussian white noiseξ(t): etc.
Stochastic differential equations Differential equations which contain (“are driven by”) random functions Langevin equations: the random function is Gaussian white noiseξ(t): etc. Simple example (Brownian motion/ Ornstein-Uhlenbeck process):
Stochastic differential equations Differential equations which contain (“are driven by”) random functions Langevin equations: the random function is Gaussian white noiseξ(t): etc. Simple example (Brownian motion/ Ornstein-Uhlenbeck process): Solution v(t) is random (because it depends on ξ(t))
Stochastic differential equations Differential equations which contain (“are driven by”) random functions Langevin equations: the random function is Gaussian white noiseξ(t): etc. Simple example (Brownian motion/ Ornstein-Uhlenbeck process): Solution v(t) is random (because it depends on ξ(t)) Want to know P[v], averages over distribution of ξ(t)
More generally, multivariate:
More generally, multivariate: higher-order:
More generally, multivariate: higher-order: nonlinear:
More generally, multivariate: higher-order: nonlinear: multiplicative noise:
Brownian motion solution (with m = 1):
Brownian motion solution (with m = 1): averages:
Brownian motion solution (with m = 1): averages:
Brownian motion solution (with m = 1): averages:
Brownian motion solution (with m = 1): averages:
Brownian motion solution (with m = 1): averages:
Brownian motion solution (with m = 1): averages:
Brown (2) equal-time correlation:
Brown (2) equal-time correlation: but from equilibrium stat mech:
Brown (2) equal-time correlation: but from equilibrium stat mech:
Brown (2) equal-time correlation: but from equilibrium stat mech: (another Einstein relation)
Brown (2) equal-time correlation: but from equilibrium stat mech: (another Einstein relation) Note: OU model also applies with v -> x (position) to overdamped motion in a parabolic potential
Solution using Fourier transform solution:
Solution using Fourier transform solution:
Solution using Fourier transform solution:
Solution using Fourier transform solution:
Solution using Fourier transform solution: inverse FT:
Solution using Fourier transform solution: inverse FT:
Solution using Fourier transform solution: inverse FT:
Solution using Fourier transform solution: inverse FT:
Solution using Fourier transform solution: inverse FT: (as in direct calculation)
Damped oscillator FT: inverse FT:
General OU process damped oscillator:
General OU process damped oscillator:
General OU process damped oscillator: Is a 2-d OU process with
General OU process damped oscillator: Is a 2-d OU process with x(t) is not a Markov process (2nd order equation), but (x(t),p(t)) is (1st order equation).
Formal solution by FT damped oscillator case: