# inordinatum

Physics and Mathematics of Disordered Systems

## A very brief introduction to directed polymers

A classical problem in the field of disordered systems is an elastic string (or elastic polymer) stretched between two points in a random environment. One is interested in knowing how the polymer geometry (roughness), ground state energy, and excitations change due to the disordered environment. A fundamental question is if the polymer becomes rough for arbitrarily weak disorder, or if there is a phase transition between a weakly disordered (smooth) and a strongly disordered (rough) phase.

In the following I will try explain some basic results on the simple case of a two-dimensional directed polymer (i.e. a “stretched” polymer which does not contain loops or overhangs). I will explain why for large polymer lengths $L$, moments of the partition sum $Z$ grow as $\ln \overline{Z^n} \propto n^3 L$. I will show why this indicates a rough polymer with roughness exponent $\zeta=\frac{2}{3}$. This analysis holds for arbitrarily weak disorder, meaning that in two dimensions disorder is always relevant, and the polymer is always in the rough phase.

Although none of these observations are new, the existing results are often buried in technical literature (there is not even a Wikipedia page describing the directed polymer problem). I hope this self-contained summary will be useful for people interested in a quick introduction to the subject. Some original references where more details can be found are cited at the end of the post.

## The directed polymer partition sum, and its moments

We can consider a polymer configuration as a trajectory for a particle moving between two points. We call the polymer directed if the particle always moves from the initial towards the final point, i.e. if the trajectory does not contain loops or overhangs.

Formally, this means that the polymer configuration in $1+d$ dimensions can be parametrized by giving $d$ transversal coordinates $(x_1...x_d)=\vec{x}$ as a function of the longitudinal coordinate $t$ (distance along the path). The partition sum is then given by a path integral

$Z = \int \mathcal{D}[\vec{x}]e^{-\int_0^L \mathrm{d}t \, \left[(\dot{\vec{x}})^2 + \eta(\vec{x},t)\right]}$

$(\dot{\vec{x}})^2$ is the elastic energy term, which is minimized for straight paths. $\eta(\vec{x},t)$ is the energy it costs the polymer to pass through the site $(\vec{x},t)$ (and which is minimized for rough paths picking out the minimal energy sites). We model a disordered environment by taking $\eta$ to be uncorrelated Gaussian white noise:

$\displaystyle \overline{\eta(\vec{x}_1,t_1)\eta(\vec{x}_2,t_2)} = 2c\delta^{(d)}(\vec{x}_1-\vec{x}_2)\delta(t_1-t_2)$

Applying the Feynman-Kac formula, the partition sum $Z(\vec{x},t)$ satisfies the following stochastic differential equation:

$\displaystyle \partial_t Z(\vec{x},t) = \nabla^2 Z(\vec{x},t) +\eta(\vec{x},t) Z(\vec{x},t)$

$Z$ is a random quantity, because $\eta$ is random. Since computing the distribution of $Z$ is difficult, we will consider its moments, defined as

$\displaystyle \Psi_n(\vec{x}_1...\vec{x}_n,t) := \overline{Z(\vec{x}_1,t)\cdots Z(\vec{x}_n,t)}$

The time evolution of $\Psi_n$ is determined by applying the Itô formula to the equation for $Z$:

$\displaystyle \partial_t\Psi_n(\vec{x}_1...\vec{x}_n,t) = \left[\sum_{j=1}^n \nabla_{\vec{x}_j}^2 + 2c\sum_{j

This is, up to an overall sign, the Schrödinger equation for $n$ bosons interacting through pairwise pointlike attraction of strength $c$. In one dimension, this is called the Lieb-Liniger model. It is easy to determine the ground state of this system exactly, which will allow us to obtain the leading behaviour of $Z^n$ for large $t$.

## Solution in one spatial dimension: Bethe ansatz ground state

In $d=1$, the ground state of the Hamiltonian $H_n$ defined above has energy $E_n = \frac{c^2}{12}n(n^2-1)$ and is given by

$\displaystyle \Psi_n(\vec{x}_1...\vec{x}_n) \propto \exp\left(-\frac{c}{2}\sum_{j

This is a special case of the general Bethe ansatz ubiquitious in modern theoretical physics. Its application to the directed polymer problem, as discussed here, was recognized by Kardar (see references below).

To check that this ansatz for $\Psi$ gives, indeed, an eigenstate, let us take two derivatives with respect to $x_j$:

$\displaystyle \begin{array}{lcl} \partial_{x_j}^2 \Psi_n & = & \partial_{x_j} \left[-\frac{c}{2}\sum_{j \neq k} \text{sgn}(x_j-x_k) \Psi_n\right] \\ & = & \left[-c\sum_{j \neq k} \delta(x_j-x_k)\right] \Psi_n + \left[-\frac{c}{2}\sum_{k \neq j} \text{sgn}(x_j-x_k)\right]^2 \Psi_n \end{array}$

Taking the sum over all $j$ we get

$\displaystyle \sum_{j=1}^n \partial_{x_j}^2 \Psi_n = -2c\sum_{j < k} \delta(x_j-x_k) \Psi_n + \frac{c^2}{4}\sum_{j=1}^n \left[\sum_{j \neq k} \text{sgn}(x_j-x_k)\right]^2 \Psi_n$

To compute the second term, note that it is symmetric under permutations of the $x_j$. Let us assume (without loss of generality) that $x_1 < x_2 < ... < x_n$. Then

$\displaystyle \begin{array}{lcl} \sum_{j=1}^n \left[\sum_{j \neq k} \text{sgn}(x_j-x_k)\right]^2 & = & \sum_{j=1}^n \left[\sum_{k>j} (-1)+\sum_{k

Thus, we see that as claimed,

$\displaystyle H_n \Psi_n = E_n \Psi_n$

with the energy $E_n = \frac{c^2}{12}n(n^2-1)$. To see that this is, indeed, the ground state, is beyond the scope of this post (essentially it would require one to take a finite system volume $V$, construct all excited states using the Bethe ansatz, and check that they all have higher energy. Since one knows the number of eigenstates in a finite volume, this would exclude any lower-energy state).

Believing that this is, indeed, the ground state our argument shows that for large polymer length $L$, the moments of the partition sum $Z$ grow as:

$\displaystyle \overline{Z^n (L)} \propto e^{\frac{c^2}{12}n(n^2-1) L}$

## Implications for the string geometry

The asymptotics of $\overline{Z^n}$ derived above is consistent with a free energy that scales with polymer length $L$ as $\ln Z = F =: f L^\frac{1}{3}$. Furthermore, the distribution $P(f)$ of the rescaled free energy should have a tail decaying as

$\displaystyle P(f) \propto e^{-f^{\frac{3}{2}}} \quad \text{for } f \gg 1$

To see this, let us compute $\overline{Z^n}$ assuming this form for $P(f)$:

$\displaystyle \overline{Z^n} = \int \mathrm{d}f e^{n f L^{\frac{1}{3}}-f^{\frac{3}{2}}}$

For large $n$, we expect $f$ to be large; thus the integral can be approximated by the maximum of the exponent. It occurs at

$f \propto \left(n L^{\frac{1}{3}}\right)^2$

and thus we get

$\displaystyle \ln \overline{Z^n} \propto L n^3$,

consistent with the large-$n$ behaviour of the explicit formula for $\overline{Z^n}$ derived above.
Thus, the asymptotics of $\overline{Z^n}$ allowed us to determine the tail of the free energy distribution $P(f)$ and the scaling $F = f L^\frac{1}{3}$. The scaling of the typical displacement $x$ with polymer length $L$ follows:

$\displaystyle L^\frac{1}{3} \sim F \sim \int_0^L \mathrm{d}t (\dot{x})^2 \sim L \left(\frac{x}{L}\right)^2 \Rightarrow x \sim L^\frac{2}{3}$

Thus, the polymer is rough in the sense that the typical fluctuations don’t remain bounded, but increase with polymer length as

$\displaystyle \overline{\left[x(L)-x(0)\right]^2}\sim L^{2\zeta}$

with the so-called roughness exponent $\zeta= \frac{2}{3}$.

This whole analysis remains valid at an arbitrarily weak coupling strength $c$. Thus, the polymer is rough for arbitrarily weak disorder (or, equivalently, arbitrarily high temperatures). In order to observe a true phase transition between weak-disorder and strong-disorder behaviour one needs to go to higher dimensions.

## Higher dimensions and phase transitions

As we have now seen, in two dimensions (one longitudinal + one transversal) the directed polymer is always in a rough phase. It turns out (but it would take too long to explain here) that the same holds in three (1+2) dimensions. However, starting from one longitudinal + three transversal dimensions, weak disorder is insufficient to make the polymer rough. In other words, there is a high-temperature phase where the moments of $Z$ scale as $\ln \overline{Z^n} \propto n^2$, indicating Gaussian fluctuations, i.e. purely thermal roughness.

## References

The original argument for obtaining the moments of $Z$ from the Bethe ansatz solution of the Lieb-Liniger model is due to Kardar (see e.g. Nucl. Phys. B 290, 1989 and the references therein). The saddle-point argument for the tail of the free energy distribution is due to Zhang (see e.g. PRB 42, 1990). If you feel some other reference should be included please let me know!

Written by inordinatum

July 29, 2012 at 9:42 pm

### 10 Responses

1. […] few posts ago I briefly talked about the problem of a directed polymer in a random medium. In that post, I discussed the scaling behaviour of the polymer fluctuations due to the disordered […]

2. A time derivative, $\partial_t \Psi_n$, seems to be missing in front of the evolution equation for $\Psi$, the Schrodinger equation. Great web notes, congratulations!

Juan M Lopez

September 16, 2016 at 3:45 pm

• Thanks, I just corrected that! Glad you enjoyed it 🙂

inordinatum

September 17, 2016 at 9:34 am

3. […] a cartesian lattice in 1+1 dimensions, the roughness exponent is known exactly (see e.g. my earlier introduction to directed polymers) and, in some cases, even the distribution of the free energy can be computed (see my earlier post […]

4. Great Introduction again!
I missed exactly how one gets the ground energy value. Is their a simple counting argument?
Thanks a lot

GUILLAUME A

January 31, 2018 at 4:27 pm

• Thanks!
The ground energy value is calculated in section “Solution in one spatial dimension: Bethe ansatz ground state” – the calculations there show that when you apply the Hamiltonian (including the Laplacian and the delta-functions sum) to the wave function $\Psi_n$ you obtain again $\Psi_n$, multiplied by $E_n$… Let me know if something is unclear there!

inordinatum

February 3, 2018 at 4:03 pm

• Oh yes you are right! I was wondering how does one finds the eigen vectors/values together by some general method… but of course it is all enclosed in the Bethe ansatz itself ! My bad, I never learned that properly, now I know at least that bit 😉
Thanks again!

GUILLAUME A

February 3, 2018 at 7:00 pm

• Let me try to ask a more ‘clever’ question this time:
So as it looks, Bethe ansatz works pretty well for delta correlated potentials, I guess we can work out generalizations ? And as a corollary, could you indicate maybe a few cases where the roughness exponent is modified in specific ways ? I am interested to delve deeper into it.
Again : Great blog!

GUILLAUME A

February 3, 2018 at 7:18 pm

• Thanks 🙂
The roughness exponent is actually pretty robust, on large scales it appears independently of the small-scale details of the system. I.e., even if you smoothen the delta-correlated potential slightly, or add small thermal fluctuations, it will remain the same.
A different roughness exponent is obtained if you modify the universality class of your model, e.g. by changing the spatial dimension or by introducing long-range correlations (e.g. long-range elastic coupling or long-range correlated disorder). See e.g. the following papers and the references therein:
– Peng, Havlin, Schwartz, Eugene Stanley 1991, “Directed-polymer and ballistic-deposition growth with correlated noise”, available on ResearchGate
– Rosso, Krauth 2002, “Roughness at the depinning threshold for a long-range elastic string”, arXiv: cond-mat/0107527
– Fedorenko, Le Doussal, Wiese 2006, “Statics and dynamics of elastic manifolds in media with long-range correlated disorder”, arxiv:cond-mat/0609234

inordinatum

February 5, 2018 at 9:16 pm

• Fascinating!

GUILLAUME A

February 6, 2018 at 10:04 am