diff --git a/docs/source/hubbard_1d.md b/docs/source/hubbard_1d.md
index d21c43f78a03adcd6d7a9009378b217633867657..8f1a7ea0f8989e4f18aa7d79ff83e9281b9b2469 100644
--- a/docs/source/hubbard_1d.md
+++ b/docs/source/hubbard_1d.md
@@ -91,7 +91,7 @@ def compute_sol(U, h_0, nk, filling=2):
     }
     guess = pymf.generate_guess(frozenset(h_int), len(list(h_0.values())[0]))
     full_model = pymf.Model(h_0, h_int, filling)
-    mf_sol = pymf.solver(full_model, guess, nk=nk)
+    mf_sol = pymf.solver(full_model, guess, nk=nk, optimizer_kwargs={"M":0})
     return pymf.add_tb(h_0, mf_sol)
 
 
diff --git a/examples/1d_hubbard_totalenergy.ipynb b/examples/1d_hubbard_totalenergy.ipynb
index 30492b2d4c8e78f3dffc8b85b5c5a5ed28143566..5464fdc76f2902005f31695001f450b175d71446 100644
--- a/examples/1d_hubbard_totalenergy.ipynb
+++ b/examples/1d_hubbard_totalenergy.ipynb
@@ -124,7 +124,7 @@
     "    }\n",
     "    guess = pymf.generate_guess(frozenset(h_int), len(list(h_0.values())[0]))\n",
     "    full_model = pymf.Model(h_0, h_int, filling)\n",
-    "    mf_sol = pymf.solver(full_model, guess, nk=nk)\n",
+    "    mf_sol = pymf.solver(full_model, guess, nk=nk, optimizer_kwargs={\"M\":0})\n",
     "    return pymf.add_tb(h_0, mf_sol)\n",
     "\n",
     "\n",