diff --git a/src/util/util_get_energy_momentum.f90 b/src/util/util_get_energy_momentum.f90
index b935a680c..ae59c158a 100644
--- a/src/util/util_get_energy_momentum.f90
+++ b/src/util/util_get_energy_momentum.f90
@@ -199,7 +199,7 @@ subroutine util_get_energy_potential_triangular(npl, lmask, GMcb, Gmass, mass, x
          do concurrent(j = i+1:npl, lmask(i) .and. lmask(j))
             pepl(j) = - (Gmass(i) * mass(j)) / norm2(xb(:, i) - xb(:, j))
          end do
-         pe = pe + sum(pepl(i+1:npl), lmask(i) .and. lmask(j))
+         pe = pe + sum(pepl(i+1:npl), lmask(i))
       end do
       !$omp end parallel do
       pe = pe + sum(pecb(1:npl), lmask(1:npl))