diff --git a/src/pair_yukawa.cpp b/src/pair_yukawa.cpp
index 69a88ff0b0..1ad3708b76 100644
--- a/src/pair_yukawa.cpp
+++ b/src/pair_yukawa.cpp
@@ -217,6 +217,7 @@ double PairYukawa::init_one(int i, int j)
     offset[i][j] = a[i][j] * screening / cut[i][j];
   } else offset[i][j] = 0.0;
 
+  a[j][i] = a[i][j];
   offset[j][i] = offset[i][j];
 
   return cut[i][j];