diff --git a/euphonic/force_constants.py b/euphonic/force_constants.py
index 5e0555cbb..b1e29a91a 100644
--- a/euphonic/force_constants.py
+++ b/euphonic/force_constants.py
@@ -372,7 +372,8 @@ def calculate_qpoint_phonon_modes(
 
         .. math::
 
-         \\frac{d\\omega_{\\mathbf{q}\\nu}}{d\\mathbf{q}} = \\frac{1}{2\\omega_{\\mathbf{q}\\nu}}{(E^{-1}\\frac{dD}{d\\mathbf{q}}E)}
+          \\frac{d\\omega_{\\mathbf{q}\\nu}}{d\\mathbf{q}} = \\frac{1}{2\\omega_{\\mathbf{q}\\nu}}{(E^{-1}\\frac{dD}{d\\mathbf{q}}E)}
+
         :math:`\\frac{dD}{d\\mathbf{q}}` can be obtained by differentiating the Fourier equation above:
 
         .. math::
diff --git a/euphonic/structure_factor.py b/euphonic/structure_factor.py
index 599c48888..97d2ddbd1 100644
--- a/euphonic/structure_factor.py
+++ b/euphonic/structure_factor.py
@@ -198,7 +198,7 @@ def calculate_sqw_map(self,
         StructureFactor.structure_factors is defined as the mode-resolved
         :math:`S(Q, \\omega_{\\mathbf{q}\\nu})` per atom of sample. To create
         an :math:`S(Q,\\omega)` map, it is binned in :math:`\\omega`
-        and the Bose population factor is applied [1]:
+        and the Bose population factor is applied [1]_:
 
         .. math::