diff --git a/src/eko/io/runcards.py b/src/eko/io/runcards.py
index 9de351c19..b83cf0db8 100644
--- a/src/eko/io/runcards.py
+++ b/src/eko/io/runcards.py
@@ -49,7 +49,7 @@ class TheoryCard(DictLike):
     n3lo_ad_variation: N3LOAdVariation
     """|N3LO| anomalous dimension variation: ``(gg, gq, qg, qq, nsp, nsm,
     nsv)``."""
-    use_fhmruvv: Optional[bool] = False
+    use_fhmruvv: Optional[bool] = True
     """If True use the |FHMRUVV| |N3LO| anomalous dimensions."""
     matching_order: Optional[Order] = None
     """Matching conditions perturbative order tuple, ``(QCD, QED)``."""
@@ -213,7 +213,7 @@ def new_theory(self):
         new["n3lo_ad_variation"] = old.get("n3lo_ad_variation", (0, 0, 0, 0, 0, 0, 0))
         # here PTO: 0 means truly LO, no QED matching is available so far.
         new["matching_order"] = old.get("PTO_matching", [old["PTO"], 0])
-        new["use_fhmruvv"] = old.get("use_fhmruvv", False)
+        new["use_fhmruvv"] = old.get("use_fhmruvv", True)
         return TheoryCard.from_dict(new)
 
     @property
diff --git a/src/ekobox/cards.py b/src/ekobox/cards.py
index 295993634..7e9774c7a 100644
--- a/src/ekobox/cards.py
+++ b/src/ekobox/cards.py
@@ -24,7 +24,7 @@
     xif=1.0,
     n3lo_ad_variation=(0, 0, 0, 0, 0, 0, 0),
     matching_order=[0, 0],
-    use_fhmruvv=False,
+    use_fhmruvv=True,
 )
 
 _operator = dict(
diff --git a/src/ekore/anomalous_dimensions/unpolarized/space_like/__init__.py b/src/ekore/anomalous_dimensions/unpolarized/space_like/__init__.py
index cf569d9f7..167110115 100644
--- a/src/ekore/anomalous_dimensions/unpolarized/space_like/__init__.py
+++ b/src/ekore/anomalous_dimensions/unpolarized/space_like/__init__.py
@@ -25,7 +25,7 @@
 
 
 @nb.njit(cache=True)
-def gamma_ns(order, mode, n, nf, n3lo_ad_variation, use_fhmruvv=False):
+def gamma_ns(order, mode, n, nf, n3lo_ad_variation, use_fhmruvv=True):
     r"""Compute the tower of the non-singlet anomalous dimensions.
 
     Parameters
@@ -100,7 +100,7 @@ def gamma_ns(order, mode, n, nf, n3lo_ad_variation, use_fhmruvv=False):
 
 
 @nb.njit(cache=True)
-def gamma_singlet(order, n, nf, n3lo_ad_variation, use_fhmruvv=False):
+def gamma_singlet(order, n, nf, n3lo_ad_variation, use_fhmruvv=True):
     r"""Compute the tower of the singlet anomalous dimensions matrices.
 
     Parameters
@@ -137,7 +137,7 @@ def gamma_singlet(order, n, nf, n3lo_ad_variation, use_fhmruvv=False):
 
 
 @nb.njit(cache=True)
-def gamma_ns_qed(order, mode, n, nf, n3lo_ad_variation, use_fhmruvv=False):
+def gamma_ns_qed(order, mode, n, nf, n3lo_ad_variation, use_fhmruvv=True):
     r"""Compute the grid of the QED non-singlet anomalous dimensions.
 
     Parameters
@@ -288,7 +288,7 @@ def choose_ns_ad_aem2(mode, n, nf, cache):
 
 
 @nb.njit(cache=True)
-def gamma_singlet_qed(order, n, nf, n3lo_ad_variation, use_fhmruvv=False):
+def gamma_singlet_qed(order, n, nf, n3lo_ad_variation, use_fhmruvv=True):
     r"""Compute the grid of the QED singlet anomalous dimensions matrices.
 
     Parameters
@@ -331,7 +331,7 @@ def gamma_singlet_qed(order, n, nf, n3lo_ad_variation, use_fhmruvv=False):
 
 
 @nb.njit(cache=True)
-def gamma_valence_qed(order, n, nf, n3lo_ad_variation, use_fhmruvv=False):
+def gamma_valence_qed(order, n, nf, n3lo_ad_variation, use_fhmruvv=True):
     r"""Compute the grid of the QED valence anomalous dimensions matrices.
 
     Parameters
diff --git a/tests/eko/evolution_operator/test_init.py b/tests/eko/evolution_operator/test_init.py
index ed003f925..69cc29e5a 100644
--- a/tests/eko/evolution_operator/test_init.py
+++ b/tests/eko/evolution_operator/test_init.py
@@ -45,7 +45,7 @@ def test_quad_ker_errors():
                 is_polarized=True,
                 is_time_like=True,
                 n3lo_ad_variation=(0, 0, 0, 0, 0, 0, 0),
-                use_fhmruvv=False,
+                use_fhmruvv=True,
             )
 
 
@@ -112,7 +112,7 @@ def test_quad_ker(monkeypatch):
                     is_polarized=p,
                     is_time_like=t,
                     n3lo_ad_variation=(0, 0, 0, 0, 0, 0, 0),
-                    use_fhmruvv=False,
+                    use_fhmruvv=True,
                 )
                 np.testing.assert_allclose(res_ns, res)
     for label in [(br.non_singlet_pids_map["ns+"], 0), (100, 100)]:
@@ -141,7 +141,7 @@ def test_quad_ker(monkeypatch):
                     is_polarized=polarized,
                     is_time_like=False,
                     n3lo_ad_variation=(0, 0, 0, 0, 0, 0, 0),
-                    use_fhmruvv=False,
+                    use_fhmruvv=True,
                 )
                 np.testing.assert_allclose(res_sv, 1.0)
     for label in [
@@ -177,7 +177,7 @@ def test_quad_ker(monkeypatch):
                 n3lo_ad_variation=(0, 0, 0, 0, 0, 0, 0),
                 is_polarized=False,
                 is_time_like=False,
-                use_fhmruvv=False,
+                use_fhmruvv=True,
             )
             np.testing.assert_allclose(res_sv, 1.0)
 
@@ -205,7 +205,7 @@ def test_quad_ker(monkeypatch):
         n3lo_ad_variation=(0, 0, 0, 0, 0, 0, 0),
         is_polarized=False,
         is_time_like=False,
-        use_fhmruvv=False,
+        use_fhmruvv=True,
     )
     np.testing.assert_allclose(res_ns, 0.0)