Local time path integrals and their application to Lévy random walks
|
|
- Aubrey Turner
- 6 years ago
- Views:
Transcription
1 Local time path integrals and their application to Lévy random walks Václav Zatloukal ( Faculty of Nuclear Sciences and Physical Engineering Czech Technical University in Prague talk given at Path Integration in Complex Dynamical Systems 8 February 27 Leiden V. Zatloukal Local time path integrals /33
2 Motivation TIME ZONES: V. Zatloukal Local time path integrals 2/33
3 Motivation TIME ZONES: NOT INCLUDED IN THIS PRESENTATION V. Zatloukal Local time path integrals 2/33
4 Motivation Local time of a stochastic process: xhtl l=2 l=.5 t - l= -2 V. Zatloukal Local time path integrals 3/33
5 Motivation Local time of a stochastic process: xhtl x x x l=2 l=.5 t LHxL LHxL LHxL l= l=2 l=.5 l= Stochastic trajectory x( )! Local time profile L(x) V. Zatloukal Local time path integrals 3/33
6 Outline Preliminaries: di usion equation and path integral Introduction of local time of a stochastic process Correlation functions and functionals of the local time Time-independent systems and the resolvent method Local time of Lévy random walks (Lévy flights) Local-time representation for Gaussian path integrals In collaboration with Petr Jizba (Czech Technical University in Prague) and Angel Alastuey (ENS Lyon). V. Zatloukal Local time path integrals 4/33
7 Di usion equation and path integral Di usion (or heat, or Fokker-Planck) equation t + H( i@ x, x, t)] P(x, t) = () where H(p, x, t) is a generic Hamiltonian, and x 2 R. [~ =] Solution for initial condition P(x, t a )= (x x a ) represented by phase-space path integral (x b t b x a t a ) P(x b, t b )= x(tz b )=x b x(t a )=x a Dx Z Dp 2 ea[p,x] (2) where the action A[p, x] = Z tb t a d ip( )ẋ( ) H(p( ), x( ), ) (3) V. Zatloukal Local time path integrals 5/33
8 Di usion equation and path integral We will gradually specify the Hamiltonian: H(p, x, t) generic, possibly time-dependent [ H(p, x) time-independent [ H(p) time and position-independent [ H (p) =D (p 2 ) /2 Lévy Hamiltonian In the end, special attention to H(p, x) = p2 2M + V (x) usual quantum Hamiltonian V. Zatloukal Local time path integrals 6/33
9 Di usion equation and path integral Significance of the heat equation: Di usion processes (or continuous random walks) P(x, t) is the probability of transition to point x after time t Quantum mechanics: by rotation to imaginary times t! it P are transition amplitudes Quantum statistical physics: identifying t = ~, where = k B T P are matrix elements of the equilibrium density matrix... In addition, corresponding path integrals describe: Polymers or other line-like objects Stock prices on financial markets... V. Zatloukal Local time path integrals 7/33
10 Local time of a stochastic process xhtl x x x l=2 l=.5 t LHxL LHxL LHxL l= l=2 l=.5 l= For stochastic trajectory x( ) define L(x; t a, t b, x( )) = Z tb t a d (x x( )) (4) P. Lévy, Sur certains processus stochastiques homogènes, Compos. Math. 7, 283(939). V. Zatloukal Local time path integrals 8/33
11 Local time of a stochastic process L is functional of x( ) and function of x: L(x) R R L(x)dx = t b t a L(x) has compact support L(x) is continuous L(x)... trajectories of a new stochastic process V. Zatloukal Local time path integrals 9/33
12 Correlation functions of local time hl(x )...L(x n )i = Z x(t b )=x b x(t a )=x a Dx = X 2S n Z t a <t <...<t n <t b dt...dt n Z Dp 2 ea[p,x] L(x )...L(x n ) ny (x (k+) t k+ x (k) t k ) (5) where are permutations of indices {,...,n}, with () = and (n + ) = n +, and we have denoted x = x a, t = t a, x n+ = x b, t n+ = t b. k= V. Zatloukal Local time path integrals / 33
13 Generic functionals of local time L(x )...L(x n ) is an example of functional of the local time F [L(x)]. In general, hf [L(x)]i = = F x(tz b )=x b x(t a )=x a Dx apple U(x) Z Dp 2 ea[p,x] F [L(x)] x(tz b )=x b x(t a )=x a Dx Z Dp 2 ea U[p,x] U= (6) where A U [p, x] corresponds to H U (p, x, t) =H(p, x, t)+u(x). We have employed the key identity Z tb t a U(x( )) d = Z U(x)L(x) dx (7) V. Zatloukal Local time path integrals / 33
14 Time-independent systems and the resolvent method Assume H(p, x) [ independent of time t ] ) (x b t b x a t a )=hx b e (t b t a )Ĥ x a i! take t a =, t b t from now on. Laplace representation hf [L(x)]i E = R dt e te hf [L(x)]i Correlation functions: hl(x )...L(x n )i E = X 2S n ny k= R(x (k), x (k+), E) (8) [ e.g. hl(x)i E = R(x a, x)r(x, x b ) ] Here R(x, x, E) hx (Ĥ + E) xi (9) denotes matrix elements of the resolvent operator ˆR(E) =(Ĥ E). V. Zatloukal Local time path integrals 2 / 33
15 Time-independent systems and the resolvent method Generic functionals: apple hf [L(x)]i E = F R U (x a, x b, E) U(x) U= where ˆR U (E) =(ĤU E) is the resolvent corresponding to the extended Hamiltonian H U = H + U., () One-point distribution: take U(x )=u (x x)! ) R U (x a, x b )=R(x a, x b ) u R(x a, x)r(x, x b ) +ur(x, x) (L L(x)) E = (L)R(x a, x)r(x, x b ) R(x, x) 2 e apple + (L) R(x a, x b ) L R(x,x) R(x a, x)r(x, x b ) R(x, x) () (2) V. Zatloukal Local time path integrals 3 / 33
16 *Multipointresolvents Resolvent identity for reads (ĤU + E) =(Ĥ + E) (ĤU + E) U(ˆx)(Ĥ + E) (3) U(x) = R U (x a, x b )=R(x a, x b ) Setting x a = x k for k =,...,n, where M jk R U (x j, x b )= nx u j (x x j ) (4) j= nx k= nx u j R(x a, x j )R U (x j, x b ) (5) j= are elements of the inverse of the matrix M jk R(x k, x b ) (6) M jk = jk + u k R(x j, x k ) (7) V. Zatloukal Local time path integrals 4 / 33
17 Distribution functions and moments of local time Invert the Laplace transform to return to the time domain: h...i E!h...i Moments of local time: µ(x,...,x n )= hl(x )...L(x n )i hi (8) n-point distribution functions of local time: Q n j= W (L,...,L n ; x,...,x n )= (L j L(x j )) hi (9) The normalization factor hi equals (x b t b x a t a ) for averaging over paths with fixed initial point x a and final point x b. Averaging over paths with arbitrary x b : h...i R dx b h...i!hi,µ, W,... V. Zatloukal Local time path integrals 5 / 33
18 Lévy random walks Hamiltonians H(p) only dependent on momentum form now on: P(x, t) (x b t x a ) = Z dp 2 e th(p) e ipx, x = x b x a (2) For the Lévy Hamiltonian H (p) =D (p 2 ) /2 [ we consider 2 [, 2] ]! symmetric Lévy stable distribution Heavy tails for For =2Gaussian: P =2 (x, t) = e x 2 4D 2 t p 4 D2 t <2: P (x, t) x! c x P l HxL l=2 l=.5 l= For =Cauchy(-Lorentz): P = (x, t) = D t (D t) 2 + x x V. Zatloukal Local time path integrals 6 / 33
19 Lévy stable distributions Stability: X Levy, X 2 Levy ) X + X 2 Levy Characteristic function: h exp ipµ cp i sgn(p) i, = ( tan( /2) 6= (2/ ) log p =...tail power ( x ), c...width, µ...shift of origin,...asymmetry Generalized central limit theorem: i.i.d. X,...,X n distribution with possibly infinite variance lim n! X X n N n Levy Applications in financial markets (fluctuation of stock prices), physics (e.g. François Bardou et al., Lévy Statistics and Laser Cooling: How Rare Events Bring Atoms to Rest, CUP(22).,...),... V. Zatloukal Local time path integrals 7 / 33
20 Lévy random walks xhtl x x x l=2 l=.5 t LHxL LHxL LHxL l= l=2 l=.5 l= For decreasing : Longer jumps (Lévy flights) Separation of local time L(x) into peaks V. Zatloukal Local time path integrals 8 / 33
21 Resolvent of Lévy Hamiltonian R (x, x, E) = Z dp 2... Linnik (or geometric stable) distribution e ip(x x) D (p 2 ) /2 + E (2) For =2 Gaussian resolvent: R l HxL.8 R 2 (x, E) = e p E/D2 x p 4D2 E.6 l=2 l=.5 For = Cauchy resolvent: apple sin E x D R (x, E) = D + cos E x D si E x D ci E x D l= divergent for = at x = x V. Zatloukal Local time path integrals 9 / 33
22 One-point distributions of local time at the origin ) x b = x a : (L L(x a )) E = e L R () =exp h LE i, D / sin W (L; x a )= = (D t)/ ( ) (L Z hi L(xa )) de e Et exp h LE cos i sin hle sin i (22) Gaussian case: For W =2 (L; x a )= 2D 2L t =(Cauchy case): singularity R () = + e D 2L 2 /t! flattening of the distribution W l HLL l=2 l=.5 l= L V. Zatloukal Local time path integrals 2 / 33
23 One-point distributions of local time at the origin 2) free x b : (L L(x a )) E = ER () e L R () = E / exp h LE i (L W (L; x a )= Z = de e Et E / L(xa )) hi (23) h exp LE cos i sin hle sin + i W l * HLL Gaussian case: W =2 (L; x a)= r 4D2 t e D 2L 2 /t l=2 l=.5 l= L V. Zatloukal Local time path integrals 2 / 33
24 *One-pointdistributionofBrownianlocaltime Gaussian case =2(Brownian motion): generic x b and x W 2 (L; x) = (L) x a x + x x b +2D 2 L apple t ( xa x + x x b +2D 2 L) 2 (x b x a ) 2 exp 4D 2 t apple ( xa x + x x b ) 2 (x b x a ) 2 + (L) exp 4D 2 t (24) r apple W2 4D2 (L; x) = (L) t exp ( xa x +2D 2 L) 2 4D 2 t + (L) erf apple xa x p 4D2 t (25) The -term vanishes for x between x a and x b. V. Zatloukal Local time path integrals 22 / 33
25 First moments of local time [ x b = x a ] [free x b ] m l HxL m l * HxL.5 l=2.5 l=2 l=.5 l=.5. l=. l= x x µ =2 (x) = s (x t b xa) 2 e 4D 2 t 4D 2 erfc " xa x + x x b p 4D2 t # µ =2 (x) = s t D 2 e (xa x) 2 4D 2 t x a x 2D 2 " # xa x erfc p 4D2 t µ = (x) = n 2 D 2 h t (D t) 2 +(xa+x b 2x) 2i h 2(x x a ) h +2(x b x) (D t) 2 +(x b x) 2 (x x a ) 2i arctan D t x b x h +D t (D t) 2 +(x x a ) 2 +(x b x) 2i ln (D t) 2 +(x x a ) 2 (x b x) 2i arctan D t x xa (x xa) 2 (x b x) h 2 (D t) 2 +(x xa) 2ih (D t) 2 +(x b x) 2 i ) ) µ = (x) = " ln + 2 D x D t x a 2 # V. Zatloukal Local time path integrals 23 / 33
26 *SecondmomentsofBrownianlocaltime ) fixed x b µ 2 (x, x 2 )= t 2D 2 +( 2 $ 2 ) apple 2 exp 2 (x b x a ) 2 4D 2 t p apple 2 p 4D2 t erfc 2 p 4D2 t apple (xb x a ) 2 exp 4D 2 t where 2 x a x + x x 2 + x 2 x b, 2 x a x 2 + x 2 x + x x b 2) free x b µ 2 (x, x 2 )= 4D 2 p apple t p 2 ( exp 2 ) 2 D2 4D 2 t + apple ( 2 ) 2 2D 2 + t erfc apple 2 p 4D2 t +( 2 $ 2) where 2 x a x + x x 2, 2 x a x 2 + x 2 x V. Zatloukal Local time path integrals 24 / 33
27 Path integrals in quantum statistical physics From now on H(p, x) = p2 2M + V (x) (26) Motivation in quantum statistical physics: Gibbs operator: e! partition function: Tr (e! thermal density matrix: e Ĥ Ĥ ), =/k B T Ĥ /Tr (e Ĥ ) (x a, x b, ) hx b e Ĥ x a i... Heat kernel for di usion generated by H, with time variable + H( i~@x, x) (x a, x, )= (27) V. Zatloukal Local time path integrals 25 / 33
28 Path integrals in quantum statistical physics ) Feynman path-integral representation For (x a, x b, )= Z x( ~)=x b x()=x a Dx( )exp!... high temperatures: ~ Z ~ d apple M 2 ẋ 2 ( )+V (x( )) (28) V (x( )) = V (x a )+V (x a )(x( ) x a )+...) e V (x a ) p 2 ~ 2 /M P. Jizba and V. Zatloukal, Path-integral approach to the Wigner-Kirkwood expansion, Phys. Rev. E 89, 235(24) 2) Spectral representation (x a, x b, )= X n e E n n(x a ) n (x b ) where Ĥ n i = E n n i (29) For!... low temperatures: e E gs gs(x a ) gs (x b ) V. Zatloukal Local time path integrals 26 / 33
29 Local-time path integral 3) Local-time representation? path integral over local time profiles! s = / ~ L(x) = Z ~ d (x x( )) (3) small large V. Zatloukal Local time path integrals 27 / 33
30 Local-time path integral (x a, x b, )= Z x( ~)=xb x()=x a Dx( )exp ( ~ Z ~ M 2 ẋ 2 ( )d ~ Z R V (X )L(X )dx ) + change of variables: x( )! L(x) (x a, x b, )= Z DL(x)W [L(x); + ~, x a, x b ]exp Z R L(x)V (x)dx where L(x) and W [L(x)] = if R R L(x)dx 6= ~! identify the weights W [L(x)] V. Zatloukal Local time path integrals 28 / 33
31 Local-time path integral: glimpses of the derivation Di usion equation in Laplace picture: E + H( i~@x, x) e (x a, x, E) = (x a x) Quantum-field-theoretic representation: e (x a, x b, E) =2 R D (x) (x a ) (x b )e h E+Ĥ i / R D (x)e h E+Ĥ i Replica trick: a/b =lim D! ab D R PD o 2 e =lim D! D D D (x) (x a ) (x b )expn = h E + Ĥ i Spherical coordinates in -space: radial part = q ~ ~ Inverse Laplace transform ) Local-time representation of (x a, x b, ) V. Zatloukal Local time path integrals 29 / 33
32 Local-time path integral (x a, x b, ) hx b e Z (X+ )= + (X )= Ĥ x a i = D (x) lim lim X ±!± Z X+ X D! 2 D 2 2 (x)dx lim ( +) ±!! D 2 (x a ) (x b )exp{ A [ (x)]} where (x) 2 (x) $ L(x)/~ The action of a radial harmonic oscillator ($ Bessel process indexed by x) A [ (x)] Z X+ X dx apple ~ 2 2M (x) 2 + V (x) 2 (x)+ M ~ 2 (x) 2 2 (x) (3) V. Zatloukal Local time path integrals 3 / 33
33 Local-time path integral at low temperatures Rescaling! p : A [ p (x)] = Z X+ X dx apple ~ 2 2M (x) 2 + V (x) 2 (x)+ M ~ 2 (x) (x)!: Saddle-point approximation of the path integral (! neglect the last term) $ Minimization of the functional: Z X+ X dx apple ~ 2 2M (x) 2 + V (x) 2 (x) = h Ĥ i under the constraint h i = ) Rayleigh-Ritz variational principle for the ground state V. Zatloukal Local time path integrals 3 / 33
34 Generic functionals of the local time Average value of generic functional F [L(x)]: Z x( ~)=xb x()=x a = lim Dx( )F [L(X )] exp lim X ±!± D! Z (X+ )= + (X )= 2 D 2 ( ~ Z ~ lim ( +) ±! D (x)f [~ 2 (x)] d D 2 apple M 2 ẋ 2 ( )+V (x( )) Z X+ X 2 (x)dx! ) (x a ) (x b )e A [ (x)] V. Zatloukal Local time path integrals 32 / 33
35 Summary We discussed local times of stochastic processes for generic Hamiltonian Calculated moments and distributions of local time In particular, for Lévy random walks For Gaussian path integrals derived local-time path-integral representation xhtl x x x l=2 l=.5 t LHxL LHxL LHxL l= -2-2 l=2-2 l=.5-2 l= P. Jizba and VZ, Local-time representation of path integrals, Phys. Rev. E 92, 6237(25)[arXiv:56.888]. VZ, Local time of Lévy random walks: a path integral approach, [arxiv: ]. V. Zatloukal Local time path integrals 33 / 33
Green Function of the Double-Fractional Fokker-Planck Equation
Green Function of the Double-Fractional Fokker-Planck Equation Hagen Kleinert, Václav Zatloukal Max Planck Institute for the History of Science, Berlin and Faculty of Nuclear Sciences and Physical Engineering,
More informationSeries Solution of Linear Ordinary Differential Equations
Series Solution of Linear Ordinary Differential Equations Department of Mathematics IIT Guwahati Aim: To study methods for determining series expansions for solutions to linear ODE with variable coefficients.
More informationAnomalous 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 informationWeak 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 informationIntroduction to Path Integrals
Introduction to Path Integrals Consider ordinary quantum mechanics of a single particle in one space dimension. Let s work in the coordinate space and study the evolution kernel Ut B, x B ; T A, x A )
More informationTheory 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 informationFrom 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 informationWHITE NOISE APPROACH TO FEYNMAN INTEGRALS. Takeyuki Hida
J. Korean Math. Soc. 38 (21), No. 2, pp. 275 281 WHITE NOISE APPROACH TO FEYNMAN INTEGRALS Takeyuki Hida Abstract. The trajectory of a classical dynamics is detrmined by the least action principle. As
More informationAdvanced Quantum Mechanics
Advanced Quantum Mechanics Rajdeep Sensarma! sensarma@theory.tifr.res.in Lecture #22 Path Integrals and QM Recap of Last Class Statistical Mechanics and path integrals in imaginary time Imaginary time
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 information1 Solution to Homework 4
Solution to Homework Section. 5. The characteristic equation is r r + = (r )(r ) = 0 r = or r =. y(t) = c e t + c e t y = c e t + c e t. y(0) =, y (0) = c + c =, c + c = c =, c =. To find the maximum value
More informationPath integrals in quantum mechanics
Path integrals in quantum mechanics Phys V3500/G8099 handout #1 References: there s a nice discussion of this material in the first chapter of L.S. Schulman, Techniques and applications of path integration.
More informationDiffusion 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 informationMath 4263 Homework Set 1
Homework Set 1 1. Solve the following PDE/BVP 2. Solve the following PDE/BVP 2u t + 3u x = 0 u (x, 0) = sin (x) u x + e x u y = 0 u (0, y) = y 2 3. (a) Find the curves γ : t (x (t), y (t)) such that that
More informationEcon 508B: Lecture 5
Econ 508B: Lecture 5 Expectation, MGF and CGF Hongyi Liu Washington University in St. Louis July 31, 2017 Hongyi Liu (Washington University in St. Louis) Math Camp 2017 Stats July 31, 2017 1 / 23 Outline
More informationMultiple Random Variables
Multiple Random Variables Joint Probability Density Let X and Y be two random variables. Their joint distribution function is F ( XY x, y) P X x Y y. F XY ( ) 1, < x
More informationFrom random matrices to free groups, through non-crossing partitions. Michael Anshelevich
From random matrices to free groups, through non-crossing partitions Michael Anshelevich March 4, 22 RANDOM MATRICES For each N, A (N), B (N) = independent N N symmetric Gaussian random matrices, i.e.
More informationFractional Quantum Mechanics and Lévy Path Integrals
arxiv:hep-ph/9910419v2 22 Oct 1999 Fractional Quantum Mechanics and Lévy Path Integrals Nikolai Laskin Isotrace Laboratory, University of Toronto 60 St. George Street, Toronto, ON M5S 1A7 Canada Abstract
More informationSome Contra-Arguments for the Use of Stable Distributions in Financial Modeling
arxiv:1602.00256v1 [stat.ap] 31 Jan 2016 Some Contra-Arguments for the Use of Stable Distributions in Financial Modeling Lev B Klebanov, Gregory Temnov, Ashot V. Kakosyan Abstract In the present paper,
More informationEuclidean path integral formalism: from quantum mechanics to quantum field theory
: from quantum mechanics to quantum field theory Dr. Philippe de Forcrand Tutor: Dr. Marco Panero ETH Zürich 30th March, 2009 Introduction Real time Euclidean time Vacuum s expectation values Euclidean
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 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 informationmacroscopic view (phenomenological) microscopic view (atomistic) computing reaction rate rate of reactions experiments thermodynamics
Rate heory (overview) macroscopic view (phenomenological) rate of reactions experiments thermodynamics Van t Hoff & Arrhenius equation microscopic view (atomistic) statistical mechanics transition state
More informationQUANTUM MARKOVIAN KINETIC EQUATION FOR HARMONIC OSCILLATOR. Boris V. Bondarev
QUANTUM MARKOVIAN KINETIC EQUATION FOR HARMONIC OSCILLATOR Boris V. Bondarev Moscow Aviation Institute, Volokolamsk road, 4, 15871, Moscow, Russia E-mail: bondarev.b@mail.ru arxiv:130.0303v1 [physics.gen-ph]
More informationLocal 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 informationMath 162: Calculus IIA
Math 62: Calculus IIA Final Exam ANSWERS December 9, 26 Part A. (5 points) Evaluate the integral x 4 x 2 dx Substitute x 2 cos θ: x 8 cos dx θ ( 2 sin θ) dθ 4 x 2 2 sin θ 8 cos θ dθ 8 cos 2 θ cos θ dθ
More informationWeak 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 information221A Lecture Notes Path Integral
A Lecture Notes Path Integral Feynman s Path Integral Formulation Feynman s formulation of quantum mechanics using the so-called path integral is arguably the most elegant. It can be stated in a single
More informationE[X n ]= dn dt n M X(t). ). What is the mgf? Solution. Found this the other day in the Kernel matching exercise: 1 M X (t) =
Chapter 7 Generating functions Definition 7.. Let X be a random variable. The moment generating function is given by M X (t) =E[e tx ], provided that the expectation exists for t in some neighborhood of
More informationC. Show your answer in part B agrees with your answer in part A in the limit that the constant c 0.
Problem #1 A. A projectile of mass m is shot vertically in the gravitational field. Its initial velocity is v o. Assuming there is no air resistance, how high does m go? B. Now assume the projectile is
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 informationGaussian integrals and Feynman diagrams. February 28
Gaussian integrals and Feynman diagrams February 28 Introduction A mathematician is one to whom the equality e x2 2 dx = 2π is as obvious as that twice two makes four is to you. Lord W.T. Kelvin to his
More informationHJB equations. Seminar in Stochastic Modelling in Economics and Finance January 10, 2011
Department of Probability and Mathematical Statistics Faculty of Mathematics and Physics, Charles University in Prague petrasek@karlin.mff.cuni.cz Seminar in Stochastic Modelling in Economics and Finance
More informationMATH Solutions to Probability Exercises
MATH 5 9 MATH 5 9 Problem. Suppose we flip a fair coin once and observe either T for tails or H for heads. Let X denote the random variable that equals when we observe tails and equals when we observe
More informationPhysics 351 Wednesday, April 22, 2015
Physics 351 Wednesday, April 22, 2015 HW13 due Friday. The last one! You read Taylor s Chapter 16 this week (waves, stress, strain, fluids), most of which is Phys 230 review. Next weekend, you ll read
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 informationGRAPH QUANTUM MECHANICS
GRAPH QUANTUM MECHANICS PAVEL MNEV Abstract. We discuss the problem of counting paths going along the edges of a graph as a toy model for Feynman s path integral in quantum mechanics. Let Γ be a graph.
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 informationMath Camp II. Calculus. Yiqing Xu. August 27, 2014 MIT
Math Camp II Calculus Yiqing Xu MIT August 27, 2014 1 Sequence and Limit 2 Derivatives 3 OLS Asymptotics 4 Integrals Sequence Definition A sequence {y n } = {y 1, y 2, y 3,..., y n } is an ordered set
More informationFinal Exam May 4, 2016
1 Math 425 / AMCS 525 Dr. DeTurck Final Exam May 4, 2016 You may use your book and notes on this exam. Show your work in the exam book. Work only the problems that correspond to the section that you prepared.
More information7. Kinetics controlled by fluctuations: Kramers theory of activated processes
7. Kinetics controlled by fluctuations: Kramers theory of activated processes Macroscopic kinetic processes (time dependent concentrations) Elementary kinetic process Reaction mechanism Unimolecular processes
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 informationPhysics 741 Graduate Quantum Mechanics 1 Solutions to Midterm Exam, Fall x i x dx i x i x x i x dx
Physics 74 Graduate Quantum Mechanics Solutions to Midterm Exam, Fall 4. [ points] Consider the wave function x Nexp x ix (a) [6] What is the correct normaliation N? The normaliation condition is. exp,
More information4. Supplementary Notes on Time and Space Evolution of a Neutrino Beam
Lecture Notes for Quantum Physics II & III 8.05 & 8.059 Academic Year 1996/1997 4. Supplementary Notes on Time and Space Evolution of a Neutrino Beam c D. Stelitano 1996 As an example of a two-state system
More informationarxiv: v1 [quant-ph] 15 Dec 2011
Sharp and Infinite Boundaries in the Path Integral Formalism Phillip Dluhy and Asim Gangopadhyaya Loyola University Chicago, Department of Physics, Chicago, IL 666 Abstract arxiv:.3674v [quant-ph 5 Dec
More informationChemistry 432 Problem Set 4 Spring 2018 Solutions
Chemistry 4 Problem Set 4 Spring 18 Solutions 1. V I II III a b c A one-dimensional particle of mass m is confined to move under the influence of the potential x a V V (x) = a < x b b x c elsewhere and
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 informationEntropy Gradient Maximization for Systems with Discoupled Time
Entropy Gradient Maximization for Systems with Discoupled Time I. V. Drozdov Abstract The conventional definition of time based on the period of oscillations or fluctuations around equilibrium states is
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 informationHamilton-Jacobi theory
Hamilton-Jacobi theory November 9, 04 We conclude with the crowning theorem of Hamiltonian dynamics: a proof that for any Hamiltonian dynamical system there exists a canonical transformation to a set of
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 informationwhere r n = dn+1 x(t)
Random Variables Overview Probability Random variables Transforms of pdfs Moments and cumulants Useful distributions Random vectors Linear transformations of random vectors The multivariate normal distribution
More informationRandom Correlation Matrices, Top Eigenvalue with Heavy Tails and Financial Applications
Random Correlation Matrices, Top Eigenvalue with Heavy Tails and Financial Applications J.P Bouchaud with: M. Potters, G. Biroli, L. Laloux, M. A. Miceli http://www.cfm.fr Portfolio theory: Basics Portfolio
More informationUNIVERSITY OF SURREY FACULTY OF ENGINEERING AND PHYSICAL SCIENCES DEPARTMENT OF PHYSICS. BSc and MPhys Undergraduate Programmes in Physics LEVEL HE2
Phys/Level /1/9/Semester, 009-10 (1 handout) UNIVERSITY OF SURREY FACULTY OF ENGINEERING AND PHYSICAL SCIENCES DEPARTMENT OF PHYSICS BSc and MPhys Undergraduate Programmes in Physics LEVEL HE PAPER 1 MATHEMATICAL,
More information1/f Fluctuations from the Microscopic Herding Model
1/f Fluctuations from the Microscopic Herding Model Bronislovas Kaulakys with Vygintas Gontis and Julius Ruseckas Institute of Theoretical Physics and Astronomy Vilnius University, Lithuania www.itpa.lt/kaulakys
More informationWinter 2019 Math 106 Topics in Applied Mathematics. Lecture 8: Importance Sampling
Winter 2019 Math 106 Topics in Applied Mathematics Data-driven Uncertainty Quantification Yoonsang Lee (yoonsang.lee@dartmouth.edu) Lecture 8: Importance Sampling 8.1 Importance Sampling Importance sampling
More informationRecursive calculation of effective potential and variational resummation
JOURNAL OF MATHEMATICAL PHYSICS 46, 032101 2005 Recursive calculation of effective potential and variational resummation Sebastian F. Brandt and Hagen Kleinert Freie Universität Berlin, Institut für Theoretische
More informationTopics for the Qualifying Examination
Topics for the Qualifying Examination Quantum Mechanics I and II 1. Quantum kinematics and dynamics 1.1 Postulates of Quantum Mechanics. 1.2 Configuration space vs. Hilbert space, wave function vs. state
More informationp-adic Feynman s path integrals
p-adic Feynman s path integrals G.S. Djordjević, B. Dragovich and Lj. Nešić Abstract The Feynman path integral method plays even more important role in p-adic and adelic quantum mechanics than in ordinary
More informationarxiv: v7 [quant-ph] 22 Aug 2017
Quantum Mechanics with a non-zero quantum correlation time Jean-Philippe Bouchaud 1 1 Capital Fund Management, rue de l Université, 75007 Paris, France. (Dated: October 8, 018) arxiv:170.00771v7 [quant-ph]
More informationTopic 5: The Difference Equation
Topic 5: The Difference Equation Yulei Luo Economics, HKU October 30, 2017 Luo, Y. (Economics, HKU) ME October 30, 2017 1 / 42 Discrete-time, Differences, and Difference Equations When time is taken to
More informationMATH1231 CALCULUS. Session II Dr John Roberts (based on notes of A./Prof. Bruce Henry) Red Center Room 3065
MATH1231 CALCULUS Session II 2007. Dr John Roberts (based on notes of A./Prof. Bruce Henry) Red Center Room 3065 Jag.Roberts@unsw.edu.au MATH1231 CALCULUS p.1/66 Overview Systematic Integration Techniques
More informationExercises. T 2T. e ita φ(t)dt.
Exercises. Set #. Construct an example of a sequence of probability measures P n on R which converge weakly to a probability measure P but so that the first moments m,n = xdp n do not converge to m = xdp.
More informationMultivariate Normal-Laplace Distribution and Processes
CHAPTER 4 Multivariate Normal-Laplace Distribution and Processes The normal-laplace distribution, which results from the convolution of independent normal and Laplace random variables is introduced by
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 informationA Toy Model. Viscosity
A Toy Model for Sponsoring Viscosity Jean BELLISSARD Westfälische Wilhelms-Universität, Münster Department of Mathematics SFB 878, Münster, Germany Georgia Institute of Technology, Atlanta School of Mathematics
More informationMATH115. Indeterminate Forms and Improper Integrals. Paolo Lorenzo Bautista. June 24, De La Salle University
MATH115 Indeterminate Forms and Improper Integrals Paolo Lorenzo Bautista De La Salle University June 24, 2014 PLBautista (DLSU) MATH115 June 24, 2014 1 / 25 Theorem (Mean-Value Theorem) Let f be a function
More informationMath 250B Final Exam Review Session Spring 2015 SOLUTIONS
Math 5B Final Exam Review Session Spring 5 SOLUTIONS Problem Solve x x + y + 54te 3t and y x + 4y + 9e 3t λ SOLUTION: We have det(a λi) if and only if if and 4 λ only if λ 3λ This means that the eigenvalues
More informationSolutions of Math 53 Midterm Exam I
Solutions of Math 53 Midterm Exam I Problem 1: (1) [8 points] Draw a direction field for the given differential equation y 0 = t + y. (2) [8 points] Based on the direction field, determine the behavior
More informationMATH115. Infinite Series. Paolo Lorenzo Bautista. July 17, De La Salle University. PLBautista (DLSU) MATH115 July 17, / 43
MATH115 Infinite Series Paolo Lorenzo Bautista De La Salle University July 17, 2014 PLBautista (DLSU) MATH115 July 17, 2014 1 / 43 Infinite Series Definition If {u n } is a sequence and s n = u 1 + u 2
More informationFeynman s path integral approach to quantum physics and its relativistic generalization
Feynman s path integral approach to quantum physics and its relativistic generalization Jürgen Struckmeier j.struckmeier@gsi.de, www.gsi.de/ struck Vortrag im Rahmen des Winterseminars Aktuelle Probleme
More informationChemistry 532 Problem Set 7 Spring 2012 Solutions
Chemistry 53 Problem Set 7 Spring 01 Solutions 1. The study of the time-independent Schrödinger equation for a one-dimensional particle subject to the potential function leads to the differential equation
More informationAdvanced Quantum Mechanics
Advanced Quantum Mechanics Rajdeep Sensarma sensarma@theory.tifr.res.in Quantum Dynamics Lecture #2 Recap of Last Class Schrodinger and Heisenberg Picture Time Evolution operator/ Propagator : Retarded
More informationMATH 23 Exam 2 Review Solutions
MATH 23 Exam 2 Review Solutions Problem 1. Use the method of reduction of order to find a second solution of the given differential equation x 2 y (x 0.1875)y = 0, x > 0, y 1 (x) = x 1/4 e 2 x Solution
More informationMicroscopic Hamiltonian dynamics perturbed by a conservative noise
Microscopic Hamiltonian dynamics perturbed by a conservative noise CNRS, Ens Lyon April 2008 Introduction Introduction Fourier s law : Consider a macroscopic system in contact with two heat baths with
More informationTransient Phenomena in Quantum Bound States Subjected to a Sudden Perturbation
Symmetry, Integrability and Geometry: Methods and Applications Vol. (5), Paper 3, 9 pages Transient Phenomena in Quantum Bound States Subjected to a Sudden Perturbation Marcos MOSHINSKY and Emerson SADURNÍ
More informationMinimum Bias Events at ATLAS
Camille Bélanger-Champagne Lehman McGill College University City University of ew York Thermodynamics Charged Particle and Correlations Statistical Mechanics in Minimum Bias Events at ATLAS Statistical
More informationt L(q, q)dt (11.1) ,..., S ) + S q 1
Chapter 11 WKB and the path integral In this chapter we discuss two reformulations of the Schrödinger equations that can be used to study the transition from quantum mechanics to classical mechanics. They
More informationLecture 6 Quantum Mechanical Systems and Measurements
Lecture 6 Quantum Mechanical Systems and Measurements Today s Program: 1. Simple Harmonic Oscillator (SHO). Principle of spectral decomposition. 3. Predicting the results of measurements, fourth postulate
More informationMathematics for Chemistry: Exam Set 1
Mathematics for Chemistry: Exam Set 1 June 18, 017 1 mark Questions 1. The minimum value of the rank of any 5 3 matrix is 0 1 3. The trace of an identity n n matrix is equal to 1-1 0 n 3. A square matrix
More informationLEGENDRE POLYNOMIALS AND APPLICATIONS. We construct Legendre polynomials and apply them to solve Dirichlet problems in spherical coordinates.
LEGENDRE POLYNOMIALS AND APPLICATIONS We construct Legendre polynomials and apply them to solve Dirichlet problems in spherical coordinates.. Legendre equation: series solutions The Legendre equation is
More informationJUST THE MATHS UNIT NUMBER DIFFERENTIATION APPLICATIONS 5 (Maclaurin s and Taylor s series) A.J.Hobson
JUST THE MATHS UNIT NUMBER.5 DIFFERENTIATION APPLICATIONS 5 (Maclaurin s and Taylor s series) by A.J.Hobson.5. Maclaurin s series.5. Standard series.5.3 Taylor s series.5.4 Exercises.5.5 Answers to exercises
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 informationUniversality of distribution functions in random matrix theory Arno Kuijlaars Katholieke Universiteit Leuven, Belgium
Universality of distribution functions in random matrix theory Arno Kuijlaars Katholieke Universiteit Leuven, Belgium SEA 06@MIT, Workshop on Stochastic Eigen-Analysis and its Applications, MIT, Cambridge,
More informationconditional cdf, conditional pdf, total probability theorem?
6 Multiple Random Variables 6.0 INTRODUCTION scalar vs. random variable cdf, pdf transformation of a random variable conditional cdf, conditional pdf, total probability theorem expectation of a random
More informationAGAT 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 informationRegular Variation and Extreme Events for Stochastic Processes
1 Regular Variation and Extreme Events for Stochastic Processes FILIP LINDSKOG Royal Institute of Technology, Stockholm 2005 based on joint work with Henrik Hult www.math.kth.se/ lindskog 2 Extremes for
More informationTOP EIGENVALUE OF CAUCHY RANDOM MATRICES
TOP EIGENVALUE OF CAUCHY RANDOM MATRICES with Satya N. Majumdar, Gregory Schehr and Dario Villamaina Pierpaolo Vivo (LPTMS - CNRS - Paris XI) Gaussian Ensembles N = 5 Semicircle Law LARGEST EIGENVALUE
More informationPHYS-454 The position and momentum representations
PHYS-454 The position and momentum representations 1 Τhe continuous spectrum-a n So far we have seen problems where the involved operators have a discrete spectrum of eigenfunctions and eigenvalues.! n
More informationJump-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 informationPath Integrals. Andreas Wipf Theoretisch-Physikalisches-Institut Friedrich-Schiller-Universität, Max Wien Platz Jena
Path Integrals Andreas Wipf Theoretisch-Physikalisches-Institut Friedrich-Schiller-Universität, Max Wien Platz 1 07743 Jena 5. Auflage WS 2008/09 1. Auflage, SS 1991 (ETH-Zürich) I ask readers to report
More information04. Random Variables: Concepts
University of Rhode Island DigitalCommons@URI Nonequilibrium Statistical Physics Physics Course Materials 215 4. Random Variables: Concepts Gerhard Müller University of Rhode Island, gmuller@uri.edu Creative
More informationFunctional 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 informationMath 5588 Final Exam Solutions
Math 5588 Final Exam Solutions Prof. Jeff Calder May 9, 2017 1. Find the function u : [0, 1] R that minimizes I(u) = subject to u(0) = 0 and u(1) = 1. 1 0 e u(x) u (x) + u (x) 2 dx, Solution. Since the
More informationLecture 3 Dynamics 29
Lecture 3 Dynamics 29 30 LECTURE 3. DYNAMICS 3.1 Introduction Having described the states and the observables of a quantum system, we shall now introduce the rules that determine their time evolution.
More informationQuantum Mechanics on Heisenberg group. Ovidiu Calin Der-Chen Chang Peter Greiner
Quantum Mechanics on Heisenberg group Ovidiu Calin Der-Chen Chang Peter Greiner I think I can safely say that no one undestands quantum mechanics Richard Feynman A little bit of History... 1901 Max Plank:the
More informationPhys 622 Problems Chapter 6
1 Problem 1 Elastic scattering Phys 622 Problems Chapter 6 A heavy scatterer interacts with a fast electron with a potential V (r) = V e r/r. (a) Find the differential cross section dσ dω = f(θ) 2 in the
More informationNon-perturbative effects in ABJM theory
October 9, 2015 Non-perturbative effects in ABJM theory 1 / 43 Outline 1. Non-perturbative effects 1.1 General aspects 1.2 Non-perturbative aspects of string theory 1.3 Non-perturbative effects in M-theory
More information1 Equal-time and Time-ordered Green Functions
1 Equal-time and Time-ordered Green Functions Predictions for observables in quantum field theories are made by computing expectation values of products of field operators, which are called Green functions
More informationStochastic (intermittent) Spikes and Strong Noise Limit of SDEs.
Stochastic (intermittent) Spikes and Strong Noise Limit of SDEs. D. Bernard in collaboration with M. Bauer and (in part) A. Tilloy. IPAM-UCLA, Jan 2017. 1 Strong Noise Limit of (some) SDEs Stochastic differential
More information