diff --git a/src/qibo/backends/abstract.py b/src/qibo/backends/abstract.py
index 8f84a39532..9ab4b62cf6 100644
--- a/src/qibo/backends/abstract.py
+++ b/src/qibo/backends/abstract.py
@@ -378,7 +378,7 @@ def calculate_singular_value_decomposition(self, matrix):  # pragma: no cover
 
     @abc.abstractmethod
     def calculate_jacobian_matrix(
-        self, circuit, parameters, return_complex: bool = True
+        self, circuit, parameters, initial_state=None, return_complex: bool = True
     ):  # pragma: no cover
         """Calculate the Jacobian matrix of ``circuit`` with respect to varables ``params``."""
         raise_error(NotImplementedError)
diff --git a/src/qibo/backends/numpy.py b/src/qibo/backends/numpy.py
index be1b8cefdf..80601ac2a9 100644
--- a/src/qibo/backends/numpy.py
+++ b/src/qibo/backends/numpy.py
@@ -800,7 +800,7 @@ def calculate_singular_value_decomposition(self, matrix):
         return self.np.linalg.svd(matrix)
 
     def calculate_jacobian_matrix(
-        self, circuit, parameters=None, return_complex: bool = True
+        self, circuit, parameters=None, initial_state=None, return_complex: bool = True
     ):
         raise_error(
             NotImplementedError,