Solution (source code)

= Solution

For $q\ne0$ and $\alpha>0$, each <Fourier mode> is an <Ornstein-Uhlenbeck process> with decay rate $r_q=\alpha q^2$. Its <explicit Ornstein-Uhlenbeck solution> is
$$
h_q(t)=e^{-r_qt}h_q(0)+\sigma\int_0^t e^{-r_q(t-s)}\,dW_q(s),
$$
so the <Itô isometry> gives
$$
\langle h_q(t)^2\rangle
=e^{-2r_qt}\langle h_q(0)^2\rangle
+\frac{\sigma^2}{2r_q}(1-e^{-2r_qt}).
$$
To identify the <overdamped Langevin dynamics>, choose a friction $\zeta_q>0$, a <harmonic oscillator> potential $V_q=\zeta_q\alpha q^2h_q^2/2$, and force noise $f_q=\zeta_q\eta_q$. The <fluctuation-dissipation relation for a Langevin particle> then defines $k_BT_{\rm eff}=\zeta_q\sigma^2/2$. Its <Boltzmann distribution> is
$$
P_q(h)=\sqrt{\frac{\alpha q^2}{\pi\sigma^2}}
\exp\!\left(-\frac{\alpha q^2h^2}{\sigma^2}\right),
\qquad \boxed{\langle h_q^2\rangle_{\rm st}=\frac{\sigma^2}{2\alpha q^2}\propto q^{-2}}.
$$
This is the <stationary spectrum of a linear fluctuating interface>, also following from the <equipartition theorem>. The numerical prefactor uses exactly the mode-noise normalization given in the paper. The nonzero modes can therefore have a stationary <Gaussian distribution> even though the unpinned zero mode cannot.