Draft Forbes Group Website (Build by Nikola). The official site is hosted at:
License: GPL3
ubuntu2004
Path Integrals
Based on a discussion with Fred Gittes, we compute the propagator for small systems using a path integral.
One of Feynman great insights was that the propagator for quantum systems can be constructed through a path integral. He expressed this by saying that the probability for a particle starting at point and time to be later at a point and time could be obtained by averaging the phases ParseError: KaTeX parse error: Undefined control sequence: \I at position 4: e^{\̲I̲ ̲S[x]} over all possible trajectories that the particle might take from to where is the classical action:
ParseError: KaTeX parse error: Undefined control sequence: \d at position 42: …x, \dot{x}, t) \̲d̲{t}, \qquad L…Formally, this can be expressed in terms of a new type of integral:
ParseError: KaTeX parse error: Undefined control sequence: \op at position 79: … = \braket{x_f|\̲o̲p̲{U}(t_f,t_i|x_i…where ParseError: KaTeX parse error: Undefined control sequence: \op at position 1: \̲o̲p̲{U}(t_f, t_i) is the quantum mechanical propagator:
ParseError: KaTeX parse error: Undefined control sequence: \op at position 19: …t{\psi(t_f)} = \̲o̲p̲{U}(t_f, t_i)\k…In the expression for ParseError: KaTeX parse error: Undefined control sequence: \op at position 1: \̲o̲p̲{U}(t_f, t_i), the integral must be time-ordered as signified by the operator which means that in every term, operators must appear in descending time order. Formally this can be understood by Taylor expanding the exponential and then manually rearranging all terms so that they are in the correct order.
For example, expanding only to quadratic order, we would have:
ParseError: KaTeX parse error: Undefined control sequence: \op at position 1: \̲o̲p̲{U}(t_f, t_i) =…Of course, if the Hamiltonian commutes at different times ParseError: KaTeX parse error: Undefined control sequence: \op at position 2: [\̲o̲p̲{H}(t_1), \op{H… – for example, if ParseError: KaTeX parse error: Undefined control sequence: \op at position 3: V(\̲o̲p̲{X}, t) = (\op{… is time independent – then there is no need to worry about the time-ordering and one can simply ignore this complication.
Consider the path integral representation of the propagator corresponding to a time-independent Hamiltonian:
ParseError: KaTeX parse error: Undefined control sequence: \I at position 75: …= \int \D{x}e^{\̲I̲ ̲S/\hbar}, \qqua…Discretizing in time we have
ParseError: KaTeX parse error: Undefined control sequence: \d at position 24: …um_{i=0}^{N-1} \̲d̲{t}\;\left[ …The matrices ParseError: KaTeX parse error: Undefined control sequence: \mat at position 2: A\̲m̲a̲t̲{M} are unitary ParseError: KaTeX parse error: Undefined control sequence: \abs at position 1: \̲a̲b̲s̲{A}^2\mat{M}^\d… if properly normalized:
ParseError: KaTeX parse error: Undefined control sequence: \mat at position 2: [\̲m̲a̲t̲{M}^\dagger\cdo…Hence (the factor of ParseError: KaTeX parse error: Undefined control sequence: \I at position 1: \̲I̲ included in will help later):
ParseError: KaTeX parse error: Undefined control sequence: \I at position 25: …{\frac{m}{2\pi \̲I̲\hbar\d{t}}}, \…Free Particle
If we neglect the potential , then we can explicitly perform the matrix multiplication, completing the square in and absorbing the integrating factor into :
ParseError: KaTeX parse error: Undefined control sequence: \mat at position 4: [(A\̲m̲a̲t̲{M})^2]_{xy} = …Note: the integral over is a little dubious since the factor of ParseError: KaTeX parse error: Undefined control sequence: \I at position 1: \̲I̲ makes it non-convergent.
Likewise:
ParseError: KaTeX parse error: Undefined control sequence: \mat at position 4: [(A\̲m̲a̲t̲{M})^N]_{xy} …Hence (including the normalization factor),
ParseError: KaTeX parse error: Undefined control sequence: \I at position 41: …{m}{2\pi \hbar \̲I̲ ̲T}}\exp\left\{\…Discretization
To make the path integral concrete (and well defined) we can express it in terms of a discrete lattice in space and time. A useful approach is to use periodic functions on a lattice of points in a box of length . With this discretization both the positions and wave-numbers become discrete with running over a finite range .
ParseError: KaTeX parse error: Undefined control sequence: \d at position 8: x_n = n\̲d̲{x} = n\frac{L}…Two sets of normalizations are useful and follow from the following lattice relationship
ParseError: KaTeX parse error: Undefined control sequence: \I at position 12: \sum_{l}e^{\̲I̲ ̲k_l(x_m - x_n)}…Physical:
ParseError: KaTeX parse error: Undefined control sequence: \I at position 23: …{x_m|k_n} = e^{\̲I̲ ̲k_n x_m}, \qqua…In this normalization, continuum equations can be directly translated by replacing integrals as follows:
ParseError: KaTeX parse error: Undefined control sequence: \d at position 6: \int \̲d̲{x} \rightarrow…Numerical:
The Numerical normalization is obtained from the Physical normalization by rescaling:
ParseError: KaTeX parse error: Undefined control sequence: \d at position 46: …et{x_n}}{\sqrt{\̲d̲{x}}}, \qquad …In this normalization is unitary.
Discrete Unitarity
Fred Gittes points out that in discrete form, there is a special relationship between ParseError: KaTeX parse error: Undefined control sequence: \d at position 1: \̲d̲{t} and ParseError: KaTeX parse error: Undefined control sequence: \d at position 1: \̲d̲{x} such that ParseError: KaTeX parse error: Undefined control sequence: \mat at position 1: \̲m̲a̲t̲{M}^\dagger\cdo… is unitary:
ParseError: KaTeX parse error: Undefined control sequence: \mat at position 4: [(A\̲m̲a̲t̲{M})^\dagger \c…where the last relationship is true for the lattice iff:
ParseError: KaTeX parse error: Undefined control sequence: \d at position 15: \frac{m}{\hbar\̲d̲{t}}\d{x}^2 = \…To see this, note that with the special ParseError: KaTeX parse error: Undefined control sequence: \d at position 1: \̲d̲{t}, the last sum is . This leaves the first factor, whose exponent vanishes when leaving a factor of unity.
The discrete matrix ParseError: KaTeX parse error: Undefined control sequence: \mat at position 1: \̲m̲a̲t̲{M} is thus unitary with this special value of ParseError: KaTeX parse error: Undefined control sequence: \d at position 1: \̲d̲{t} if the normalization factor of at this special value. As we shall see below,
$$A = \sqrt{\frac{m}{2\pi \I\hbar\d{t}}} \\$$Using the numerical basis, we can diagonalize . ParseError: KaTeX parse error: Undefined control sequence: \I at position 120: …\left\{ -\̲I̲ ̲(k_mx - k_ny) …
Free Particle on a Periodic Lattice
The last sum can be simplified if we can shift ParseError: KaTeX parse error: Undefined control sequence: \d at position 20: …ightarrow q + l\̲d̲{x} to complete the square which requires
ParseError: KaTeX parse error: Undefined control sequence: \d at position 30: …r k_n}{m}\frac{\̲d̲{t}}{\d{x}}to be an integer. This is the case if ParseError: KaTeX parse error: Undefined control sequence: \d at position 1: \̲d̲{t} is chosen in terms of an integer :
ParseError: KaTeX parse error: Undefined control sequence: \d at position 1: \̲d̲{t} = j\frac{m …If this is true, then we can do this shift and we have:
ParseError: KaTeX parse error: Undefined control sequence: \I at position 54: …ft\{ \frac{\̲I̲ ̲m l^2 \d{x}^2}{…Harmonic Oscillator
Harmonic Oscillator
---------------------------------------------------------------------------
AssertionError Traceback (most recent call last)
<ipython-input-8-1846013a711c> in <module>()
21 M = dx*np.exp(1j/h*(m*(x_-y_)**2/2.0/dt - V(x_+y_)))/B
22 I = M.conj().T.dot(M)
---> 23 assert np.allclose(np.eye(Nx), I)
AssertionError: