THE MYSTERY OF STOCHASTIC MECHANICS. Edward Nelson Department of Mathematics Princeton University

Similar documents
Physics 235 Chapter 2. Chapter 2 Newtonian Mechanics Single Particle

An Introduction to Malliavin calculus and its applications

Hamilton- J acobi Equation: Explicit Formulas In this lecture we try to apply the method of characteristics to the Hamilton-Jacobi equation: u t

Linear Response Theory: The connection between QFT and experiments

The expectation value of the field operator.

Lecture 4 Kinetics of a particle Part 3: Impulse and Momentum

Some Basic Information about M-S-D Systems

Math Week 14 April 16-20: sections first order systems of linear differential equations; 7.4 mass-spring systems.

Lecture #11: Wavepacket Dynamics for Harmonic Oscillator

Physics 127b: Statistical Mechanics. Fokker-Planck Equation. Time Evolution

Predator - Prey Model Trajectories and the nonlinear conservation law

2.3 SCHRÖDINGER AND HEISENBERG REPRESENTATIONS

Problem Set on Differential Equations

ψ(t) = V x (0)V x (t)

Homogenization of random Hamilton Jacobi Bellman Equations

14 Autoregressive Moving Average Models

Module 3: The Damped Oscillator-II Lecture 3: The Damped Oscillator-II

Two Coupled Oscillators / Normal Modes

IB Physics Kinematics Worksheet

Major Concepts. Brownian Motion & More. Chemical Kinetics Master equation & Detailed Balance Relaxation rate & Inverse Phenomenological Rate

Fractional Method of Characteristics for Fractional Partial Differential Equations

Applications of the Basic Equations Chapter 3. Paul A. Ullrich

Lecture 10: The Poincaré Inequality in Euclidean space

Lecture 9: Advanced DFT concepts: The Exchange-correlation functional and time-dependent DFT

The Quantum Theory of Atoms and Molecules: The Schrodinger equation. Hilary Term 2008 Dr Grant Ritchie

From Particles to Rigid Bodies

Homework Solution Set # 3. Thursday, September 22, Textbook: Claude Cohen Tannoudji, Bernard Diu and Franck Lalo, Second Volume Complement G X

Lecture 20: Riccati Equations and Least Squares Feedback Control

Class Meeting # 10: Introduction to the Wave Equation

Hamilton Jacobi equations

where the coordinate X (t) describes the system motion. X has its origin at the system static equilibrium position (SEP).

Announcements: Warm-up Exercise:

What Ties Return Volatilities to Price Valuations and Fundamentals? On-Line Appendix

Hamilton- J acobi Equation: Weak S olution We continue the study of the Hamilton-Jacobi equation:

Optimizing heat exchangers

A. Using Newton s second law in one dimension, F net. , write down the differential equation that governs the motion of the block.

6. Stochastic calculus with jump processes

2001 November 15 Exam III Physics 191

KINEMATICS IN ONE DIMENSION

Math 334 Fall 2011 Homework 11 Solutions

An Introduction to Backward Stochastic Differential Equations (BSDEs) PIMS Summer School 2016 in Mathematical Finance.

Backward stochastic dynamics on a filtered probability space

A quantum method to test the existence of consciousness

k 1 k 2 x (1) x 2 = k 1 x 1 = k 2 k 1 +k 2 x (2) x k series x (3) k 2 x 2 = k 1 k 2 = k 1+k 2 = 1 k k 2 k series

Math 333 Problem Set #2 Solution 14 February 2003

System of Linear Differential Equations

Chapter 12: Velocity, acceleration, and forces

Chapter 3 Boundary Value Problem

Summary of shear rate kinematics (part 1)

Chapter 2. First Order Scalar Equations

Cash Flow Valuation Mode Lin Discrete Time

( ) = b n ( t) n " (2.111) or a system with many states to be considered, solving these equations isn t. = k U I ( t,t 0 )! ( t 0 ) (2.

Second quantization and gauge invariance.

Solutions of Sample Problems for Third In-Class Exam Math 246, Spring 2011, Professor David Levermore

In this chapter the model of free motion under gravity is extended to objects projected at an angle. When you have completed it, you should

On Option Pricing by Quantum Mechanics Approach

Undetermined coefficients for local fractional differential equations

Optimal Path Planning for Flexible Redundant Robot Manipulators

PROBLEMS FOR MATH 162 If a problem is starred, all subproblems are due. If only subproblems are starred, only those are due. SLOPES OF TANGENT LINES

Vanishing Viscosity Method. There are another instructive and perhaps more natural discontinuous solutions of the conservation law

2.1: What is physics? Ch02: Motion along a straight line. 2.2: Motion. 2.3: Position, Displacement, Distance

Math 527 Lecture 6: Hamilton-Jacobi Equation: Explicit Formulas

Random Walk on Circle Imagine a Markov process governing the random motion of a particle on a circular

arxiv: v1 [math.pr] 19 Feb 2011

Electrical Circuits. 1. Circuit Laws. Tools Used in Lab 13 Series Circuits Damped Vibrations: Energy Van der Pol Circuit

Week 1 Lecture 2 Problems 2, 5. What if something oscillates with no obvious spring? What is ω? (problem set problem)

Lecture 2-1 Kinematics in One Dimension Displacement, Velocity and Acceleration Everything in the world is moving. Nothing stays still.

Stochastic Model for Cancer Cell Growth through Single Forward Mutation

t is a basis for the solution space to this system, then the matrix having these solutions as columns, t x 1 t, x 2 t,... x n t x 2 t...

Math 10B: Mock Mid II. April 13, 2016

Traveling Waves. Chapter Introduction

Estimation of Poses with Particle Filters

We here collect a few numerical tests, in order to put into evidence the potentialities of HBVMs [4, 6, 7].

Introduction D P. r = constant discount rate, g = Gordon Model (1962): constant dividend growth rate.

Utility maximization in incomplete markets

3, so θ = arccos

Today in Physics 218: radiation reaction

4.5 Constant Acceleration

(1) (2) Differentiation of (1) and then substitution of (3) leads to. Therefore, we will simply consider the second-order linear system given by (4)

EXPLICIT TIME INTEGRATORS FOR NONLINEAR DYNAMICS DERIVED FROM THE MIDPOINT RULE

2. Nonlinear Conservation Law Equations

L07. KALMAN FILTERING FOR NON-LINEAR SYSTEMS. NA568 Mobile Robotics: Methods & Algorithms

Ch1: Introduction and Review

1. VELOCITY AND ACCELERATION

HOMEWORK # 2: MATH 211, SPRING Note: This is the last solution set where I will describe the MATLAB I used to make my pictures.

WEEK-3 Recitation PHYS 131. of the projectile s velocity remains constant throughout the motion, since the acceleration a x

Section 3.8, Mechanical and Electrical Vibrations

Solution: b All the terms must have the dimension of acceleration. We see that, indeed, each term has the units of acceleration

III. Direct evolution of the density: The Liouville Operator

Two Popular Bayesian Estimators: Particle and Kalman Filters. McGill COMP 765 Sept 14 th, 2017

AP Calculus BC Chapter 10 Part 1 AP Exam Problems

Zürich. ETH Master Course: L Autonomous Mobile Robots Localization II

LAB 6: SIMPLE HARMONIC MOTION

Chapter 8 The Complete Response of RL and RC Circuits

Solutions to Assignment 1

MATH 4330/5330, Fourier Analysis Section 6, Proof of Fourier s Theorem for Pointwise Convergence

MOMENTUM CONSERVATION LAW

Optima and Equilibria for Traffic Flow on a Network

u(x) = e x 2 y + 2 ) Integrate and solve for x (1 + x)y + y = cos x Answer: Divide both sides by 1 + x and solve for y. y = x y + cos x

4.6 One Dimensional Kinematics and Integration

Transcription:

THE MYSTERY OF STOCHASTIC MECHANICS Edward Nelson Deparmen of Mahemaics Princeon Universiy 1

Classical Hamilon-Jacobi heory N paricles of various masses on a Euclidean space. Incorporae he masses in he fla Riemannian meric m ij, he mass ensor. Then if v i is a velociy, v i = m ij v j is a momenum. Kineic energy: 1 2 vi v i. Poenial energy: V. Lagrangian: L = 1 2 vi v i V. Posiion a ime of he configuraion: ξ(). Iniial ime:. Final ime: 1. Hamilon s principal funcion: S(x, ) = 1 L ( ξ(s) ) ds. Subsanial derivaive (derivaive along rajecories): Then DS = L. D = + vi i. 2

Vecor field: v wih dξ d = v. Principle of leas acion in Hamilon-Jacobi heory: v is a criical poin for S, for unconsrained variaions. Tha is, le v be anoher vecor field, le δv = v v, and denoe by a prime quaniies wih v replaced by v. Then D(S S) = D S DS + (D D )S = L L δv i i S = L L δv i i S + o(δv). Now L L = v i δv i + o(δv), S S = 1 (v i i S)δv i ds + o(δv). Since his is rue for all variaions, we have he Hamilon-Jacobi condiion: v i = i S. 3

Togeher wih DS = L; i.e., s + vi i S = 1 2 vi v i V, his gives he Hamilon-Jacobi equaion S + 1 2 i i S + V = 0. If we ake he gradien we obain Newon s equaion F = ma. 4

Sochasic Hamilon-Jacobi heory Following Guerra and Morao, consruc a conservaive Markovian dynamics for a Markov process ξ. Wiener process (Brownian moion) w wih Edw i ()dw i () = hd. Kinemaics: dξ() = b ( ξ(), ) d + dw(). Dynamics: δe L d = 0 (heurisically). The rajecories are non-differeniable, so wha is he meaning of 1 2 dξ i d dξ i d in he Lagrangian L? 5

Le d > 0 and dξ() = ξ( + d) ξ(), so dξ d is a quoien, no a derivaive. Compue E 1 dξ i dξ i up o o(1). 2 d d dξ i = +d Noe: dw is of order d 1/2. b i( ξ(r), r ) dr + dw i. dξ i = where +d + dw i b i( ξ() + r = b i d + k b i W k + dw i + O(d 2 ) b ( ) ξ(s), s)ds + w(r) w(), r dr W k = +d [w k (r) w k ()] dr. Therefore 6

1 2 dξi dξ i = 1 2 bi b i d + b i dw i d + k b i W k dw i + 1 2 dwi dw i + o(d 2 ). Now EW k dw i = hδ k i +d (r )d = h 2 δk i, so k b i W k = h 2 ib i. The erm b i dw i d is singular, of order d 3/2, bu is expecaion is 0. Finally, Edw i dw i = h 2 nd. Hence we have he sough-for resul: E 1 2 dξ i d dξ i d = 1 2 bi b i + h 2 ib i + h 2 n d + o(1). 7

The singular erm h n is a consan, no depending on he rajecory, and i drops ou of he 2 d variaion. Le and L + = 1 2 bi b i + h 2 ib i V I = E 1 L + ( ξ(s) ) ds. Le δb be a vecor field, b = b + δb, and as before denoe by primes quaniies wih b replaced by b. Le ξ saisfy dξ () = b ( ξ (), ) d + dw() wih ξ () = ξ() for he iniial ime. Definiion: ξ is criical for L in case I I = o(δb). 8

Sochasic Hamilon s principal funcion: Then S(x, ) = E x, 1 DS = L + L + ( ξ(s), s ) ds. where now D is he mean forward derivaive: ( ) df() Df = lim E. d 0+ d and As before, D(S S) = D S DS + (D D )S Now = L + L + δb i i S = L + L + δb i i S + o(δb). L + L + = b i δb i + h 2 iδb i + o(δb) ( ) I I = E 1 ( b i i S + h ) 2 i δb i ds + o(δb). 9

Digression on Markovian kinemaics: Markov process: given he presen, pas and fuure are independen. The wo direcions of ime are on an equal fooing. In addiion o he forward drif b here is he backward drif b. v = b + b 2 u = b b 2 curren velociy, osmoic velociy. Le ρ be he probabiliy densiy of ξ. Then v = (vρ) u = h ρ 2 ρ curren equaion, osmoic equaion. End of digression. 10

The equaion (*), namely I I = E 1 ( b i i S + h ) 2 i δb i ds + o(δb), is awkward because i involves i δb i. Inegrae by pars: E h 2 iδb i( ξ(s), s ) h ( = i δb i) ρ dx R n 2 = R δbi u i ρ dx. n Since b u = v, I I = E 1 (v i i S)δb i ds + o(δb). This is rue for all variaions, so ξ is criical for L if and only if he sochasic Hamilon-Jacobi condiion holds: v i = i S. 11

Le R = h 2 log ρ, so i R = u i. If we wrie ou DS = L + and express everyhing in erms of R and S, we obain he sochasic Hamilon-Jacobi equaion: (1) S + 1 2 i S i S + V 1 2 i R i R h R = 0. 2 Expressing he curren equaion v = (vρ) in erms of R and S we obain: (2) R + ir i S + h S = 0. 2 These wo equaions are a sysem of coupled nonlinear parial differenial equaions expressing necessary and sufficien condiions for a Markov process o be criical. How can we solve hem? 12

Le ψ = e 1 h (R+iS). Then he sysem is equivalen o he Schrödinger equaion ψ = ī h ] [ h2 + V ψ. 13

magneic fields Riemannian manifolds spin Bose-Einsein and Fermi-Dirac saisics exisence of finie-energy Markov processes momenum inerference 14

Non-localiy Consider wo paricles on R 1 wih iniial wave funcion ψ(0) = 1 2π e 1 4 σ 1 (0)x x where σ 1 (0) = ( 1 ) 1 1 2 and x = ( x 1 x 2 ). The x 1 and x 2 axes are unrelaed; he paricles may be separaed by an arbirarily large amoun a. A ime 0, urn on a linear resoring force (harmonic oscillaor wih circular frequency ω) for he second paricle. Then he paricles are dynamically uncoupled. 15

Since he paricles are very widely separaed and dynamically uncoupled, we should expec ha does no depend on ω. Eξ 1 ()ξ 1 (0) In fac i does no o fourh order in, bu neverheless he rajecory of he firs paricle is immediaely affeced by he choice of ω in a far disan place. For me his is unphysical, especially since he effec does no depend on he separaion a. 16

A wrong predicion Consider wo harmonic oscillaors, abou wo widely separaed poins a 1 and a 2, wih circular frequency 1, and le X i be he Heisenberg posiion operaor, in unis of disance from a i, wih Heisenberg momenum operaor P i. Then for i = 1, 2. X i () = cos X i (0) + sin P i (0) P i () = sin X i (0) + cos P i (0) Le he Heisenberg sae vecor ψ 0 be a real Gaussian cenered a (0, 0), and wrie A = (ψ 0, Aψ 0 ). Then X i () = P i () = 0 since his is rue for = 0. The operaors X 1 () and X 2 (s) commue. Choose ψ 0 so ha he correlaion X 1 (0)X 2 (0) is 90%. Thus he oscillaors are enangled bu dynamically uncoupled. The quanum mechanical correlaion funcion X 1 ()X 2 (0) is periodic of period 2π since X 1 () is. Hence X 1 ()X 2 (0) =.9 whenever is a muliple of 2π. 17

Bu le (ξ 1, ξ 2 ) be he corresponding Markov process of sochasic mechanics. This is a diffusion process and i evenually loses all memory of where i sared. Thus whereas lim E( ξ 1 (2πn)ξ 2 (0) ) = 0 n X 1 (2πn)X 2 (0) =.9. Here we have an empirical difference beween he predicions of quanum mechanics and sochasic mechanics. Measuremens of he posiion of he firs paricle a ime and of he second paricle a ime 0 do no inerfere wih each oher, and he wo heories predic oally differen saisics. Does anyone doub ha quanum mechanics is righ and sochasic mechanics is wrong? 18