diff --git a/openmc/deplete/abc.py b/openmc/deplete/abc.py
index 7bf7f1086a..784023f26f 100644
--- a/openmc/deplete/abc.py
+++ b/openmc/deplete/abc.py
@@ -754,8 +754,9 @@ class Integrator(ABC):
         rates = res.rates[0]
         k = ufloat(res.k[0, 0], res.k[0, 1])
 
-        # Scale reaction rates by ratio of source rates
-        rates *= source_rate / res.source_rate
+        if res.source_rate != 0.0:
+            # Scale reaction rates by ratio of source rates
+            rates *= source_rate / res.source_rate
         return bos_conc, OperatorResult(k, rates)
 
     def _get_start_data(self):