diff --git a/src/apps/molresponse/FrequencyResponse.cpp b/src/apps/molresponse/FrequencyResponse.cpp
index eb7b09e25d1..76b77f24b67 100644
--- a/src/apps/molresponse/FrequencyResponse.cpp
+++ b/src/apps/molresponse/FrequencyResponse.cpp
@@ -272,7 +272,7 @@ auto FrequencyResponse::update_response(World &world, X_space &chi, XCOperator<d
 
     auto [new_res, bsh_norms] = update_residual(world, chi, new_chi, r_params.calc_type(), old_residuals, xres_old);
     inner_to_json(world, "r_x", response_context.inner(new_res, new_res), iter_function_data);
-    if (iteration >= 0) {// & (iteration % 3 == 0)) {
+    if (iteration >= 0 && iteration <r_params.maxsub()+5) {// & (iteration % 3 == 0)) {
         new_chi = kain_x_space_update(world, chi, new_res, kain_x_space);
     }
     chi_norm=new_chi.norm2s();
diff --git a/src/apps/molresponse/testing/runners.hpp b/src/apps/molresponse/testing/runners.hpp
index 47690670063..0e98c1f92b0 100644
--- a/src/apps/molresponse/testing/runners.hpp
+++ b/src/apps/molresponse/testing/runners.hpp
@@ -555,7 +555,7 @@ void set_frequency_response_parameters(World &world, ResponseParameters &r_param
         }
         //r_params.set_user_defined_value("archive", std::string("../restartdata"));
         r_params.set_user_defined_value("maxiter", size_t(35));
-        r_params.set_user_defined_value("maxsub", size_t(3));
+        r_params.set_user_defined_value("maxsub", size_t(10));
         r_params.set_user_defined_value("kain", true);
         r_params.set_user_defined_value("omega", frequency);
         r_params.set_user_defined_value("first_order", true);