diff --git a/Sample/Interface/AutocorrelationModels.cpp b/Sample/Interface/AutocorrelationModels.cpp
index 1fdcf2de8c3d17ba7a11fa5080e6216196761269..aa62772a1a8f97c1e933e0cafe7fdd06796809a1 100644
--- a/Sample/Interface/AutocorrelationModels.cpp
+++ b/Sample/Interface/AutocorrelationModels.cpp
@@ -108,5 +108,7 @@ double K_CorrelationModel::spectralFunction(const R3& q) const
 double K_CorrelationModel::rms() const
 {
     // integration of spectral function: rms^2 = Integrate[PSD(f) * 2*pi*f, {f, 0, max_frequency}]
-    return 0;
+    const double H = m_hurst_parameter;
+    const double val = 2 * pi * m_lateral_corr_length * m_max_spatial_frequency;
+    return m_sigma * std::sqrt(1. - std::pow(1 + val * val, -H));
 }