diff --git a/hoomd/hpmc/external/external.py b/hoomd/hpmc/external/external.py
index 02af8f8e9a..2d2e0de6cb 100644
--- a/hoomd/hpmc/external/external.py
+++ b/hoomd/hpmc/external/external.py
@@ -76,7 +76,7 @@ def energy(self):
 
         .. math::
 
-            U = \\sum_{i=0}^\\mathrm{N_particles-1}
+            U = \\sum_{i=0}^{N_\\mathrm{particles}-1}
             U_{\\mathrm{external},i}
 
         See `hoomd.hpmc.integrate` for the full expression which includes
diff --git a/hoomd/hpmc/integrate.py b/hoomd/hpmc/integrate.py
index 1735b5facb..e842e8847c 100644
--- a/hoomd/hpmc/integrate.py
+++ b/hoomd/hpmc/integrate.py
@@ -157,7 +157,7 @@
 .. math::
 
     U_\mathrm{external} =
-        \sum_{i=0}^\mathrm{N_particles-1} U_{\mathrm{external},i}(\vec{r}_i,
+        \sum_{i=0}^\mathrm{N_{particles}-1} U_{\mathrm{external},i}(\vec{r}_i,
                                                                  \mathbf{q}_i)
 
 Potential classes in `hoomd.hpmc.pair` evaluate
diff --git a/hoomd/hpmc/pair/pair.py b/hoomd/hpmc/pair/pair.py
index 96a6ff75a3..600ed956e7 100644
--- a/hoomd/hpmc/pair/pair.py
+++ b/hoomd/hpmc/pair/pair.py
@@ -83,8 +83,8 @@ def energy(self):
 
         .. math::
 
-            U = \\sum_{i=0}^\\mathrm{N_particles-1}
-            \\sum_{j=i+1}^\\mathrm{N_particles-1}
+            U = \\sum_{i=0}^{N_\\mathrm{particles}-1}
+            \\sum_{j=i+1}^{N_\\mathrm{particles}-1}
             U_{\\mathrm{pair},ij}
 
         See `hoomd.hpmc.integrate` for the full expression which includes
diff --git a/hoomd/md/external/wall.py b/hoomd/md/external/wall.py
index 572ab6991f..9cddf0518d 100644
--- a/hoomd/md/external/wall.py
+++ b/hoomd/md/external/wall.py
@@ -6,7 +6,7 @@
 
 .. math::
 
-    U_\mathrm{wall} = \sum_{i=0}^{\mathrm{N_particles-1}}
+    U_\mathrm{wall} = \sum_{i=0}^{N_\mathrm{particles}-1}
                       \sum_{w \in walls} U_w(d_i),
 
 where :math:`d_i` is the signed distance between particle :math:`i` and the wall
diff --git a/hoomd/md/integrate.py b/hoomd/md/integrate.py
index ed4a95c18a..63231bb120 100644
--- a/hoomd/md/integrate.py
+++ b/hoomd/md/integrate.py
@@ -379,7 +379,7 @@ def linear_momentum(self):
 
         .. math::
 
-            \\vec{p} = \\sum_{i=0}^\\mathrm{N_particles-1} m_i \\vec{v}_i
+            \\vec{p} = \\sum_{i=0}^{N_\\mathrm{particles}-1} m_i \\vec{v}_i
         """
         v = self._cpp_obj.computeLinearMomentum()
         return (v.x, v.y, v.z)
diff --git a/hoomd/md/many_body.py b/hoomd/md/many_body.py
index 3f085468ca..c735e8c652 100644
--- a/hoomd/md/many_body.py
+++ b/hoomd/md/many_body.py
@@ -6,7 +6,7 @@
 
 .. math::
 
-    U_\mathrm{many-body} = \frac{1}{2} \sum_{i=0}^\mathrm{N_particles-1}
+    U_\mathrm{many-body} = \frac{1}{2} \sum_{i=0}^{N_\mathrm{particles}-1}
                            \sum_{j \ne i}
                            \sum_{j \ne k} U(\vec{r}_{ij}, \vec{r}_{ik})
 
diff --git a/hoomd/md/pair/__init__.py b/hoomd/md/pair/__init__.py
index 6344f6cf90..ee09159bce 100644
--- a/hoomd/md/pair/__init__.py
+++ b/hoomd/md/pair/__init__.py
@@ -6,7 +6,7 @@
 
 .. math::
 
-    U_\mathrm{pair,total} = \frac{1}{2} \sum_{i=0}^\mathrm{N_particles-1}
+    U_\mathrm{pair,total} = \frac{1}{2} \sum_{i=0}^{N_\mathrm{particles}-1}
                       \sum_{j \ne i, (i,j) \notin \mathrm{exclusions}}
                       U_\mathrm{pair}(r_{ij})
 
diff --git a/hoomd/md/pair/aniso.py b/hoomd/md/pair/aniso.py
index bd08565d7f..3a6d1e5c62 100644
--- a/hoomd/md/pair/aniso.py
+++ b/hoomd/md/pair/aniso.py
@@ -6,7 +6,7 @@
 
 .. math::
 
-    U_\mathrm{pair,total} = \frac{1}{2} \sum_{i=0}^\mathrm{N_particles-1}
+    U_\mathrm{pair,total} = \frac{1}{2} \sum_{i=0}^{N_\mathrm{particles}-1}
                       \sum_{j \ne i, (i,j) \notin \mathrm{exclusions}}
                       U_\mathrm{pair}(r_{ij}, \mathbf{q}_i, \mathbf{q}_j)