diff --git a/qiskit_addon_sqd/counts.py b/qiskit_addon_sqd/counts.py
index a996c87..45d4513 100644
--- a/qiskit_addon_sqd/counts.py
+++ b/qiskit_addon_sqd/counts.py
@@ -25,11 +25,10 @@ def counts_to_arrays(counts: dict[str, float | int]) -> tuple[np.ndarray, np.nda
         counts: The counts dictionary to convert
 
     Returns:
-        A tuple containing:
-            - A 2D array representing the sampled bitstrings. Each row represents a
-              bitstring, and each element is a ``bool`` representation of the
-              bit's value
-            - A 1D array containing the probability with which each bitstring was sampled
+        - A 2D array representing the sampled bitstrings. Each row represents a
+          bitstring, and each element is a ``bool`` representation of the
+          bit's value
+        - A 1D array containing the probability with which each bitstring was sampled
 
     """
     if not counts:
diff --git a/qiskit_addon_sqd/fermion.py b/qiskit_addon_sqd/fermion.py
index d14f721..b62a7cb 100644
--- a/qiskit_addon_sqd/fermion.py
+++ b/qiskit_addon_sqd/fermion.py
@@ -64,11 +64,10 @@ def solve_fermion(
         verbose: A verbosity level between 0 and 10
 
     Returns:
-            A tuple containing:
-                - Minimum energy from SCI calculation
-                - SCI coefficients
-                - Average orbital occupancy
-                - Expectation value of spin-squared
+        - Minimum energy from SCI calculation
+        - SCI coefficients
+        - Average orbital occupancy
+        - Expectation value of spin-squared
 
     """
     if isinstance(bitstring_matrix, tuple):
@@ -168,10 +167,9 @@ def optimize_orbitals(
         learning_rate: The learning rate to use during gradient descent
 
     Returns:
-        A tuple containing:
-            - The groundstate energy found during the last optimization iteration
-            - An optimized 1D array defining the orbital transform
-            - Average orbital occupancy
+        - The groundstate energy found during the last optimization iteration
+        - An optimized 1D array defining the orbital transform
+        - Average orbital occupancy
 
     """
     if isinstance(bitstring_matrix, tuple):