diff --git a/src/ekore/operator_matrix_elements/polarized/space_like/as2.py b/src/ekore/operator_matrix_elements/polarized/space_like/as2.py
index 402274a0c..eab9b1b00 100644
--- a/src/ekore/operator_matrix_elements/polarized/space_like/as2.py
+++ b/src/ekore/operator_matrix_elements/polarized/space_like/as2.py
@@ -116,7 +116,7 @@ def A_hq_ps(n, cache, L, nf):
         * ((n + 2) / (n**2 * (n + 1) ** 2))
         * (
             (n - 1) * (2 * S2 + zeta2)
-            - ((4 * n**3 - 4 * n**2 - 3 * n - 1) / (n**2 * (n + 1) ** 2))
+            + ((4 * n**3 - 4 * n**2 - 3 * n - 1) / (n**2 * (n + 1) ** 2))
         )
     )
     # term that differentiates between M scheme and Larin scheme,
diff --git a/tests/ekore/operator_matrix_elements/polarized/space_like/test_nnlo.py b/tests/ekore/operator_matrix_elements/polarized/space_like/test_nnlo.py
index c107aff89..697f90fd3 100644
--- a/tests/ekore/operator_matrix_elements/polarized/space_like/test_nnlo.py
+++ b/tests/ekore/operator_matrix_elements/polarized/space_like/test_nnlo.py
@@ -53,18 +53,18 @@ def test_hg():
 def test_hq():
     refs = {
         0: [
-            -0.4642489711934156,
-            -0.23400164609053498,
-            -0.13602229539775093,
-            -0.08807148818984994,
-            -0.061448213442336426,
+            -0.5439814814814814,
+            -0.2517004115226338,
+            -0.14165471418159098,
+            -0.09036414682928498,
+            -0.06254577811991613,
         ],
         10: [
-            -24.229681069958847,
-            -10.410791769547325,
-            -5.706357572365681,
-            -3.587147848500778,
-            -2.459209109825822,
+            -24.309413580246915,
+            -10.428490534979424,
+            -5.7119899911495216,
+            -3.589440507140213,
+            -2.460306674503402,
         ],
     }
     zqq_shift = [