L-3

Goal: This lecture is dedicated to a classical model in disordered systems: the directed polymer in random media. It has been introduced to model vortices in superconductur or domain wall in magnetic film. We will focus here on the algorithms that identify the ground state or compute the free energy at temperature T, as well as, on the Cole-Hopf transformation that map this model on the KPZ equation.

Polymers, interfaces and manifolds in random media

We consider the following potential energy

$E_{pot}=\int dr{\frac {1}{2}}(\nabla h)^{2}+V(h(r),r)$

The first term represents the elasticity of the manifold and the second term is the quenched disorder, due to the impurities. In general, the medium is D-dimensional, the internal coordinate of the manifold is d-dimensional and the height filed is N-dimensional. Hence,the following equations always holds:

$D=d+N$

In practice, we will study two cases:

• Directed Polymers ($d=1$), $D=1+N$. Examples are vortices, fronts...
• Elastic interfaces ($N=1$), $D=d+1$. Examples are domain walls...

Today we restrict to polymers. Note that they are directed because their configuration $h(r)$ is uni-valuated. It is useful to study the model using the following change of variable

$h\to x,\quad r\to t$

Directed polymers

Dijkstra Algorithm and transfer matrix

Sketch of the discrete Directed Polymer model. At each time the polymer grows either one step left either one step right. A random energy $V(\tau ,x)$ is associated at each node and the total energy is simply $E[x(\tau )]=\sum _{\tau =0}^{t}V(\tau ,x)$.

We introduce a lattice model for the directed polymer (see figure). In a companion notebook we provide the implementation of the powerful Dijkstra algorithm.

Dijkstra allows to identify the minimal energy among the exponential number of configurations $x(\tau )$

$E_{\min }=\min _{x(\tau )}E[x(\tau )].$

We are also interested in the ground state configuration $x_{\min }(\tau )$. For both quantities we expect scale invariance with two exponents $\theta ,\zeta$ for the energy and for the roughness

$E_{\min }=c_{\infty }t+b_{\infty }t^{\theta }\chi ,\quad x_{\min }(t/2))\sim a_{\infty }t^{\zeta }{\tilde {\chi }}$

Universal exponents: Both $\theta ,\zeta$ are Independent of the lattice, the disorder distribution, the elastic constants, or the boudanry conditions. Note that $\omega =\theta$, while for an interface $\omega =d\theta$.

Non-universal constants: $c_{\infty },b_{\infty },a_{\infty }$ are of order 1 and depend on the lattice, the disorder distribution, the elastic constants... However $c_{\infty }$ is independent on the boudanry conditions!

Universal distributions: $\chi ,{\tilde {\chi }}$ are instead universal, but depends on the boundary condtions. Starting from 2000 a magic connection has been revealed between this model and the smallest eigenvalues of random matrices. In particular I discuss two different boundary conditions:

• Droplet: $x(\tau =0)=x(\tau =t)=0$. In this case, up to rescaling, $\chi$ is distributed as the smallest eigenvalue of a GUE random matrix (Tracy Widom distribution $F_{2}(\chi )$)

• Flat: $x(\tau =0)=0$ while the other end $x(\tau =t)$ is free. In this case, up to rescaling, $\chi$ is distributed as the smallest eigenvalue of a GOE random matrix (Tracy Widom distribution $F_{1}(\chi )$)

Entropy and scaling relation

It is useful to compute the entropy

${\text{Entropy}}=\ln {\binom {t}{\frac {t-x}{2}}}\approx t\ln 2-{\frac {x^{2}}{t}}+O(x^{4})$

From which one could guess from dimensional analysis

$\theta =2\zeta -1$

We will see that this relation is actually exact.

Back to the continuum model

Let us consider polymers $x(\tau )$ of length $t$, starting in $0$ and ending in $x$ and at thermal equlibrium at temperature $T$. The partition function of the model writes as

$Z(x,t)=\int _{x(0)=0}^{x(t)=x}{\cal {D}}x(\tau )\exp \left[-{\frac {1}{T}}\int _{0}^{t}d\tau {\frac {1}{2}}(\partial _{\tau }x)^{2}+V(x(\tau ),\tau )\right]$

For simplicity, we assume a white noise, ${\overline {V(x,\tau )V(x',\tau ')}}=D\delta (x-x')\delta (\tau -\tau ')$. Here, the partition function is written as a sum over all possible paths, corresponding to all possible polymer configurations that start at $0$ and end at $x$, weighted by the appropriate Boltzmann factor.

Polymer partition function and propagator of a quantum particle

Let's perform the following change of variables: $\tau =it'$. We also identifies $T$ with $\hbar$ and ${\tilde {t}}=-it$ as the fine time.

$Z(x,{\tilde {t}})=\int _{x(0)=0}^{x({\tilde {t}})=x}{\cal {D}}x(t')\exp \left[{\frac {i}{\hbar }}\int _{0}^{\tilde {t}}dt'{\frac {1}{2}}(\partial _{t'}x)^{2}-V(x(t'),t')\right]$

Note that $S[x]=\int _{0}^{\tilde {t}}dt'{\frac {1}{2}}(\partial _{t'}x)^{2}-V(x(t'),t')$ is the classical action of a particle with kinetic energy ${\frac {1}{2}}(\partial _{\tau }x)^{2}$ and time dependent potential $V(x(\tau ),\tau )$, evolving from time zero to time ${\tilde {t}}$. From the Feymann path integral formulation, Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle Z[x,\tilde t]} is the propagator of the quantum particle.

In absence of disorder, one can find the propagator of the free particle, that, in the original variables, writes:

Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle Z_{\text{free}}(x,t)=\frac{e^{-x^2/(2Tt)}}{\sqrt{2 Tt}} }

Feynman-Kac formula

Let's derive the Feyman Kac formula for Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle Z(x,t)} in the general case:

• First, focus on free paths and introduce the following probability

Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle P[A,x,t] =\int_{x(0)=0}^{x(t)=x} {\cal D} x(\tau) \exp\left[- \frac{1}{T} \int_0^t d \tau \frac{1}2(\partial_\tau x)^2 \right] \delta\left( \int_0^t d \tau V(x(\tau),\tau)-A \right) }

• Second, the moments generating function

Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle Z_p(x ,t) = \int_{-\infty}^\infty d A e^{-p A} P[A,x,t] =\int_{x(0)=0}^{x(t)=x} {\cal D} x(\tau) e^{-\frac{1}{T} \int_0^t d \tau \frac{1}2(\partial_\tau x)^2 -p \int_0^t d \tau V(x(\tau),\tau)} }

• Third, the backward approach. Consider free paths evolving up to Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle t+dt} and reaching Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle x}  :

Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle Z_p(x,t+dt)= \left\langle e^{-p \int_0^{t+dt} d \tau V(x(\tau),\tau)}\right\rangle= \left\langle e^{-p \int_0^{t} d \tau V(x(\tau),\tau)}\right\rangle e^{-p V(x,t) d t } =[1 -p V(x,t) d t +\dots]\left\langle Z_p(x-\Delta x,t) \right\rangle_{\Delta x} }

Here $\langle \ldots \rangle$ is the average over all free paths, while $\langle \ldots \rangle _{\Delta x}$ is the average over the last jump, namely Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle \langle \Delta x \rangle=0 } and Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle \langle \Delta x^2 \rangle=T d t } .

• At the lowest order we have

Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle Z_p(x,t+dt)= Z_p(x,t) +dt \left[ \frac{T}{2} \partial_x^2 Z_p -p V(x,t) Z_p \right] +O(dt^2) }

Replacing Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle p=1/T} we obtain the partition function is the solution of the Schrodinger-like equation:

Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle \partial_t Z(x,t) =- \hat H Z = - \left[ -\frac{T}{2}\frac{d^2 }{d x^2} + \frac{V(x,\tau)}{T}\right] Z(x,t) }

Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle Z[x,t=0]=\delta(x) }

Remarks

Remark 1:

This equation is a diffusive equation with multiplicative noise. Edwards Wilkinson is instead a diffusive equation with additive noise. The Cole Hopf transformation allows to map the diffusive equation with multiplicative noise in a non-linear equation with additive noise. We will apply this tranformation and have a surprise.

Remark 1: This hamiltonian is time dependent because of the multiplicative noise Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle V(x,\tau)/T} . For a time independent hamiltonian, we can use the spectrum of the operator. In general we will have to parts:

• A discrete set of eigenvalues Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle E_n} with the eigenstates Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle \psi_n(x)}
• A continuum part where the states Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle \psi_E(x)} have energy Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle E} . We define the density of states (DOS) Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle \rho(E)} , such that the number of states with energy in Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle (E, E+dE)} is Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle \rho(E) dE } .

In this case Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle Z[x,t] } can be written has the sum of two contributions:

Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle Z[x,t] = \left( e^{- \hat H t} \right)_{0 \to x}= \sum_n \psi_n(0) \psi_n^*(x) e^{- E_n t} + \int_0^\infty dE \, \rho(E) \psi_E(0) \psi_E^*(x) e^{- E t}. }

Exercise: free particle in 1D

• Step 1: The spectrum For a free particle there a no discrete eigenvalue, but only a continuum spectrum. In 1D the hamiltonian is Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle \hat H =-(T/2) \partial_x^2} .

Shows that the continuum spectrum has the form

Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle \psi_k(x) =\frac{1}{\sqrt{2 \pi}} e^{i k x} ,\quad \text{with} \; E_k =\frac{ T k^2}{2} }

Here Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle k } is a real number, while Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle E_k} is a positive number. Note that the states are delocalized on the entire real axis and they cannot be normalized.

• Step 2: From Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle k} to Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle E} For a given energy Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle E > } there are two states

Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle \psi_E^{L}(x) =\frac{1}{\sqrt{2 \pi}} e^{-i \sqrt{2 E/T} x} ,\quad \text{and} \psi_E^{R}(x) =\frac{1}{\sqrt{2 \pi}} e^{i \sqrt{2 E/T} x} }

Where the supercripts stay for left and right. Using Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle \rho(E) \int_\{-\infty}^\infty dk \delta(E-E_k) } shows that for both left and right :

Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle \rho(E) = \frac{1}{\sqrt{2 E T } }

• Use the spectral decomposition of the propagator to recover

Cole Hopf Transformation

Replacing

• Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle T =2 \nu }
• Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle x = r }
• Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle Z(x,t) = \exp\left(\frac{\lambda}{2 \nu} h(r,t) \right) }
• Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle - V(x,t)=\lambda \eta(r,t) }

You get

Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle \partial_t h(r,t)= \nu \nabla^2 h(r,t)+ \frac{\lambda}{2} (\nabla h)^2 + \eta(r,t) }

The KPZ equation!

We can establish a KPZ/Directed polymer dictionary, valid in any dimension. Let us remark that the free energy of the polymer is

Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle F= - T \ln Z(x,t) = \frac{-1}{\lambda} h(r,t) }

At low temperature, the free energy approaches the ground state energy, Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle E_{\min}} .

Dictionary

KPZ	KPZ exponents	Directed polymer	Directed polymer exponents
Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle r }	Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle r \sim t^{1/z}}	Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle x }	Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle x\sim t^{\zeta}}
Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle t}	Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle h(r,t) \sim t^{\alpha/z}}	$t$
Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle h}	Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle h(r,t) \sim r^{\alpha}}	Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle F, E_{\min}}	Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle \overline{(E_{\min} - \overline{E_{\min}})^2} \sim t^{2\theta} }

We conclude that

Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle \theta =\alpha/z, \quad \zeta=1/z }

Moreover, the scaling relation Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle \theta =2 \zeta- 1 } is a reincarnation of the Galilean invariance Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://wikimedia.org/api/rest_v1/":): {\displaystyle \alpha +z =2 } .