TBan-III

From Disordered Systems Wiki
Jump to navigation Jump to search

Exercise 2: Edwards-Wilkinson Interface with Stationary Initial Condition

Consider an Edwards-Wilkinson interface in 1+1 dimensions, at temperature T, and of length L with periodic boundary conditions:

∂h(x,t)∂t=ν∇2h(x,t)+η(x,t)

where η(x,t) is a Gaussian white noise with zero mean and variance:

⟨η(x,t)η(x′,t′)⟩=2Tδ(x−x′)δ(t−t′)

The solution can be written in Fourier space as:

h^q(t)=h^q(0)e−νq2t+∫0tdse−νq2(t−s)ηq(s)

with Fourier decomposition:

h^q(t)=1L∫0Leiqxh(x,t),h(x,t)=∑qe−iqxh^q(t),⟨ηq1(t′)ηq2(t)⟩=2TLδq1,−q2δ(t−t′)

where q=2πn/L,n=…,−1,0,1,….

In class, we computed the width of the interface starting from a flat interface at t=0, i.e., h(x,0)=0. The mean square displacement of a point h(x,t) is similar but includes also the contribution of the zero mode. The result is:

⟨Δhflat2⟩=2TLt+{T2tπν,t≪L2,TνL12,t≫L2.

The first term describes the diffusion of the center of mass, while the second comes from the non-zero Fourier modes.

Now consider the case where the initial interface h(x,0) is drawn from the equilibrium distribution at temperature T:

Pstat.[h]∝exp⁡[−ν2T∫0Ldx(∂xh)2]

For simplicity, set the initial center of mass to zero: h^q=0(0)=0. We consider the mean square displacement of the point h(x=0,t). The average is performed over both the thermal noise ⟨⋅⟩ and the initial condition ⋅‾:

⟨Δh2⟩‾=⟨[h(0,t)−h(0,0)]2⟩‾=⟨h2(0,t)⟩‾+⟨h2(0,0)⟩‾−2⟨h(0,t)h(0,0)⟩‾

Questions:

  • Compute the ensemble average of the Gaussian initial condition:
h^q1(0)h^q2(0)‾

Hint: Write the integral in terms of Fourier modes and use ∫0Ldxeiqx=Lδq,0.

  • Show that:
⟨h2(0,0)⟩‾=TνL∑q≠01q2,⟨h(0,t)h(0,0)⟩‾=TνL∑q≠0e−νq2tq2
  • Show that:
⟨h2(0,t)⟩‾=A+⟨Δhflat2⟩

where the term A depends only on the initial condition. Show that:

A(t)=TνL∑q≠0e−2νq2tq2
  • Hence write:
C(t)≡⟨Δh2⟩‾−⟨Δhflat2⟩=2TνL∑n=1∞(1−e−ν(2πn/L)2t)2(2πn/L)2

Estimate C(t) for t≫L2.

  • Estimate C(t) for t≪L2 and large L.

Hint: Write the series as an integral using the continuum variable z=2πn/L. It is helpful to know:

∫0∞ds(1−e−s2)2s2=π(2−2)

Provide the two asymptotic behaviors of ⟨Δh2⟩‾.