diff --git a/src/pybamm/simulation.py b/src/pybamm/simulation.py
index 7bd1e5eea0..24ae029032 100644
--- a/src/pybamm/simulation.py
+++ b/src/pybamm/simulation.py
@@ -1022,8 +1022,8 @@ def step(
     def _get_esoh_solver(self, calc_esoh):
         if (
             calc_esoh is False
-            or isinstance(self._model, pybamm.lead_acid.BaseModel)
-            or isinstance(self._model, pybamm.equivalent_circuit.Thevenin)
+            or not isinstance(self._model, pybamm.lithium_ion.BaseModel)
+            or self._model.options["particle phases"] not in ["1", ("1", "1")]
             or self._model.options["working electrode"] != "both"
         ):
             return None