diff --git a/src/DIELECTRIC/pppm_dielectric.cpp b/src/DIELECTRIC/pppm_dielectric.cpp
index e02cc33162..8eac7a9f78 100644
--- a/src/DIELECTRIC/pppm_dielectric.cpp
+++ b/src/DIELECTRIC/pppm_dielectric.cpp
@@ -292,7 +292,7 @@ void PPPMDielectric::qsum_qsq(int warning_flag)
   double qsqsume;
 
 #if defined(_OPENMP)
-#pragma omp parallel for default(shared) reduction(+:qsum_local,qsqsum_local)
+#pragma omp parallel for default(shared) reduction(+:qsum_local,qsqsum_local,qsqsume_local)
 #endif
   for (int i = 0; i < nlocal; i++) {
     qsum_local += q[i];
diff --git a/src/DIELECTRIC/pppm_disp_dielectric.cpp b/src/DIELECTRIC/pppm_disp_dielectric.cpp
index e5149ae427..8a2d786770 100644
--- a/src/DIELECTRIC/pppm_disp_dielectric.cpp
+++ b/src/DIELECTRIC/pppm_disp_dielectric.cpp
@@ -572,7 +572,7 @@ void PPPMDispDielectric::qsum_qsq(int warning_flag)
   double qsqsume;
 
 #if defined(_OPENMP)
-#pragma omp parallel for default(shared) reduction(+:qsum_local,qsqsum_local)
+#pragma omp parallel for default(shared) reduction(+:qsum_local,qsqsum_local,qsqsume_local)
 #endif
   for (int i = 0; i < nlocal; i++) {
     qsum_local += q[i];