L-4: Difference between revisions

From Disordered Systems Wiki
Jump to navigation Jump to search
Line 4: Line 4:
= KPZ : from 1d to the Cayley tree=
= KPZ : from 1d to the Cayley tree=


We know a lot about KPZ, but still we have much to understand:
We know a lot about KPZ, but there is still much to understand:


* In <math>d=1</math> we found <math>\theta=1/3</math> and a glassy regime present at all temperatures. The stationary solution of the KPZ equation describes, at long times, the fluctions of quantities like  <math>E_{\min}[x] - E_{\min}[x']</math>. However it does not identify the actual distribution of <math> E_{\min}</math> for a given <math>x</math>. In particular we have no idea from where Tracy Widom comes from.
* In <math>d=1</math>, we have found <math>\theta=1/3</math> and a glassy regime present at all temperatures. The stationary solution of the KPZ equation describes, at long times, the fluctuations of quantities such as <math>E_{\min}[x] - E_{\min}[x']</math>. However, it does not determine the actual distribution of <math>E_{\min}</math> for a given <math>x</math>. In particular, we have no clear understanding of the origin of the Tracy-Widom distribution.


* In <math>d=\infty</math>, there is an exact solution for the Cayley tree that predicts a freezing transition to an 1RSB phase (<math>\theta=0</math>).
* In <math>d=\infty</math>, an exact solution exists for the Cayley tree, predicting a freezing transition to a 1RSB phase (<math>\theta=0</math>).


* In finite dimension, but larger than 1, there are no exact solutions. Numerical simulations find <math>\theta >0</math> in <math>d=2</math>. The case <math>d>2</math> is very interesting.
* In finite dimensions greater than one, no exact solutions are available. Numerical simulations indicate <math>\theta > 0</math> in <math>d=2</math>. The case <math>d > 2</math> remains particularly intriguing.


==Let's do replica!==
==Let's do replica!==
To make progress in disordered systems we have to go through the moments of the partition function. For simplicity we consider polymers starting in <math>0</math> and ending in  <math>x</math>. We recall that  
To make progress in disordered systems, we need to analyze the moments of the partition function. For simplicity, we consider polymers starting at <math>0</math> and ending at <math>x</math>. We recall that
* <math>V(x,\tau)</math> is a Gaussian field with
 
* <math>V(x,\tau)</math> is a Gaussian field with
<center> <math>
<center> <math>
\overline{V(x,\tau)}=0, \quad  \overline{V(x,\tau) V(x',\tau')} = D \delta^d(x-x') \delta(\tau-\tau')  
\overline{V(x,\tau)}=0, \quad  \overline{V(x,\tau) V(x',\tau')} = D \delta^d(x-x') \delta(\tau-\tau')  
</math></center>
</math></center>
* From the Wick theorem, for a generic Gaussian <math> W </math> field we have
 
* From Wick's theorem, for a generic Gaussian field <math> W </math>, we have
<center><math>
<center><math>
\overline{\exp(W)} = \exp\left[\overline{W} +\frac{1}{2} (\overline{W^2}-\overline{W}^2)\right] </math></center>
\overline{\exp(W)} = \exp\left[\overline{W} +\frac{1}{2} \left(\overline{W^2}-\overline{W}^2\right)\right]
</math></center>
 
===The first moment===
===The first moment===
The first moment of the partition function is simple to compute and corresponds to a single replica
The first moment of the partition function is straightforward to compute and corresponds to a single replica
<center> <math>
 
\overline{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\right]  \overline{\exp\left[- \frac{1}{T} \int d \tau V(x(\tau),\tau ) \right]}
<center>
</math></center>
<math>
Note that the term <math> T^2 \overline{W^2} = \int d \tau_1 d\tau_2 \overline{V(x,\tau_1)V(x,\tau_2)}= D t \delta_0</math> has a short distance divergence due to the delta-function. Hence we can write:
\overline{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\right]  \overline{\exp\left[- \frac{1}{T} \int d \tau V(x(\tau),\tau ) \right]}
<center> <math>
</math>
\overline{Z(x,t) } = \frac{1}{(2 \pi t T)^{d/2}}\exp\left[ -\frac{1}{2} \frac{ x^2}{t T} \right]  \exp\left[ \frac{D  t \delta_0}{2T^2}  \right]
</center>
</math></center>
 
Note that the term <math> T^2 \overline{W^2} = \int d \tau_1 d\tau_2 \overline{V(x,\tau_1)V(x,\tau_2)}= D t \delta_0</math> exhibits a short-distance divergence due to the delta function. Hence, we can write:
 
<center>
<math>
\overline{Z(x,t) } = \frac{1}{(2 \pi t T)^{d/2}}\exp\left[ -\frac{1}{2} \frac{ x^2}{t T} \right]  \exp\left[ \frac{D  t \delta_0}{2T^2}  \right]
</math> 
</center> 
 
=== The second moment === 
 
For the second moment, there are two replicas: 
 
* **Step 1: The second moment is** 
 
<center> 
<math> 
\overline{Z(x,t)^2 } =\int {\cal D} x_1\int  {\cal D} x_2 \exp\left[-  \int_0^t d \tau  \frac{1}{2T}[(\partial_\tau x_1)^2+ (\partial_\tau x_2)^2] \right]  \overline{\exp\left[- \frac{1}{T} \int_0^t d \tau_1 V(x_1(\tau_1),\tau_1 ) - \frac{1}{T} \int_0^t d \tau_2 V(x_2(\tau_2),\tau_2 )\right]} 
</math>
</center>


=== The second moment ===
* **Step 2: Using Wick's theorem, we obtain** 
For the second moment instead the are two replica


* Step 1: The second moment is
<center>
<center> <math>
<math>
\overline{Z(x,t)^2 } =\int {\cal D} x_1\int  {\cal D} x_2 \exp\left[-  \int_0^t d \tau  \frac{1}{2T}[(\partial_\tau x_1)^2+ (\partial_\tau x_2)^2 \right]  \overline{\exp\left[- \frac{1}{T} \int_0^t d \tau_1 V(x_1(\tau_1),\tau_1 ) - \frac{1}{T} \int_0^t d \tau_2 V(x_2(\tau_2),\tau_2 )\right]}
\overline{Z(x,t)^2 } = \exp\left[ \frac{D  t \delta_0}{T^2}  \right]\int {\cal D} x_1\int  {\cal D} x_2 \exp\left[-  \int_0^t d \tau  \frac{1}{2T}[(\partial_\tau x_1)^2+ (\partial_\tau x_2)^2 - \frac{D}{T^2} \delta^d[x_1(\tau)-x_2(\tau)]\right]
</math></center>
</math>
</center>


* Step 2: Use Wick and derive:
* **Step 3: Changing coordinates** <math>X=(x_1+x_2)/2; \; u=x_1-x_2</math>, we get 
<center> <math>
\overline{Z(x,t)^2 } = \exp\left[ \frac{D  t \delta_0}{T^2}  \right]\int {\cal D} x_1\int  {\cal D} x_2 \exp\left[-  \int_0^t d \tau  \frac{1}{2T}[(\partial_\tau x_1)^2+ (\partial_\tau x_2)^2 - \frac{D}{T^2} \delta^d[x_1(\tau)-x_2(\tau)]\right]
</math></center>


* Step 3: Now  change coordinate <math>X=(x_1+x_2)/2; \; u=x_1-x_2</math> and get:
<center>
<center> <math>
<math>
\overline{Z(x,t)^2} = (\overline{Z(x,t)})^2 \frac{\int_{u(0)=0}^{u(t)=0} {\cal D} u  \exp\left[-  \int_0^t d \tau  \frac{1}{4T}(\partial_\tau u)^2- \frac{D}{T^2} \delta^d[u(\tau)]\right]}{\int_{u(0)=0}^{u(t)=0} {\cal D} u  \exp\left[-  \int_0^t d \tau  \frac{1}{4T}(\partial_\tau u)^2\right]}
\overline{Z(x,t)^2} = (\overline{Z(x,t)})^2 \frac{\int_{u(0)=0}^{u(t)=0} {\cal D} u  \exp\left[-  \int_0^t d \tau  \frac{1}{4T}(\partial_\tau u)^2- \frac{D}{T^2} \delta^d[u(\tau)]\right]}{\int_{u(0)=0}^{u(t)=0} {\cal D} u  \exp\left[-  \int_0^t d \tau  \frac{1}{4T}(\partial_\tau u)^2\right]}
</math></center>
</math>
</center>


==Discussion==
==Discussion==

Revision as of 10:18, 2 February 2025

Goal : final lecture on KPZ and directed polymers at finite dimension. We will show that for d>2 a "glass transition" takes place.


KPZ : from 1d to the Cayley tree

We know a lot about KPZ, but there is still much to understand:

  • In d=1, we have found θ=1/3 and a glassy regime present at all temperatures. The stationary solution of the KPZ equation describes, at long times, the fluctuations of quantities such as Emin[x]Emin[x]. However, it does not determine the actual distribution of Emin for a given x. In particular, we have no clear understanding of the origin of the Tracy-Widom distribution.
  • In d=, an exact solution exists for the Cayley tree, predicting a freezing transition to a 1RSB phase (θ=0).
  • In finite dimensions greater than one, no exact solutions are available. Numerical simulations indicate θ>0 in d=2. The case d>2 remains particularly intriguing.

Let's do replica!

To make progress in disordered systems, we need to analyze the moments of the partition function. For simplicity, we consider polymers starting at 0 and ending at x. We recall that:

  • V(x,τ) is a Gaussian field with
V(x,τ)=0,V(x,τ)V(x,τ)=Dδd(xx)δ(ττ)
  • From Wick's theorem, for a generic Gaussian field W, we have
exp(W)=exp[W+12(W2W2)]

The first moment

The first moment of the partition function is straightforward to compute and corresponds to a single replica:

Z(x,t)=x(0)=0x(t)=x𝒟x(τ)exp[1T0tdτ12(τx)2]exp[1TdτV(x(τ),τ)]

Note that the term T2W2=dτ1dτ2V(x,τ1)V(x,τ2)=Dtδ0 exhibits a short-distance divergence due to the delta function. Hence, we can write:

Z(x,t)=1(2πtT)d/2exp[12x2tT]exp[Dtδ02T2]

The second moment

For the second moment, there are two replicas:

  • **Step 1: The second moment is**

Z(x,t)2=𝒟x1𝒟x2exp[0tdτ12T[(τx1)2+(τx2)2]]exp[1T0tdτ1V(x1(τ1),τ1)1T0tdτ2V(x2(τ2),τ2)]

  • **Step 2: Using Wick's theorem, we obtain**

Z(x,t)2=exp[Dtδ0T2]𝒟x1𝒟x2exp[0tdτ12T[(τx1)2+(τx2)2DT2δd[x1(τ)x2(τ)]]

  • **Step 3: Changing coordinates** X=(x1+x2)/2;u=x1x2, we get

Z(x,t)2=(Z(x,t))2u(0)=0u(t)=0𝒟uexp[0tdτ14T(τu)2DT2δd[u(τ)]]u(0)=0u(t)=0𝒟uexp[0tdτ14T(τu)2]

Discussion

Hence, the quantity Z(x,t)2/(Z(x,t))2 can be computed.

  • The denominator u(0)=0u(t)=0𝒟uexp[0tdτ14T(τu)2] is the free propagator and gives a contribution (4Tt)d/2 .
  • Let us define the numerator
W(0,t)=u(0)=0u(t)=0𝒟uexp[0tdτ14T(τu)2DT2δd[u(τ)]]

Remark 1: From Valentina's lecture, remember that if

Z(x,t)2(Z(x,t))2=1

the partition function is self-averaging and lnZ(x,t)=lnZ(x,t). The condition above is sufficient but not necessary. It is enough that Z(x,t)2/(Z(x,t))2<const, when t, to have the equivalence between annealed and quenched averages.

Remark 2: From Feynman-Kac we can write the following equation

tW(x,t)=H^W(x,t)

Here the Hamiltonian reads:

H^=2T2DT2δd[u]

The single particle potential is time independent and actractive .

W(x,t)=x|exp(H^t)|0

At large times the behaviour is dominatated by the low energy part of the spectrum.

  • In d2 an actractive potential always gives a bound state. In particular the ground state has a negative energy E0<0. Hence at large times
W(x,t)=e|E0|t

grows exponentially. This means that at all temperature, when t

lnZ(x,t)lnZ(x,t)
  • For d>2 the low part of the spectrum is controlled by the strength of the prefactor DT2. At high temperature we have a continuum positive spectrum, at low temperature bound states exist. Hence, when t
{lnZ(x,t)=lnZ(x,t)forT>TclnZ(x,t)lnZ(x,t)forT<Tc

This transition, in d=3, is between a high temeprature, θ=0 phase and a low temeprature θ>0 no RSB phase.