 Solution

For $q\ne 0$ and $\alpha >0$, each Fourier mode is an Ornstein-Uhlenbeck process with decay rate $r_q=\alpha q^2$. Its explicit Ornstein-Uhlenbeck solution isso the Itô isometry givesTo identify the overdamped Langevin dynamics, choose a friction ${\zeta }_q>0$, a harmonic oscillator potential $V_q={\zeta }_q\alpha q^2{h}_q^2/2$, and force noise $f_q={\zeta }_q{\eta }_q$. The fluctuation-dissipation relation for a Langevin particle then defines $k_B{T}_{\mathrm{eff}}={\zeta }_q{\sigma }^2/2$. Its Boltzmann distribution isThis 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.