diff --git a/include/teqp/models/multifluid.hpp b/include/teqp/models/multifluid.hpp
index cf72dff1e875b6d35a1b54a615ca29461c77714d..daec32a065e4b35f4f0f9025e03e7393ca00fd7d 100644
--- a/include/teqp/models/multifluid.hpp
+++ b/include/teqp/models/multifluid.hpp
@@ -42,6 +42,7 @@ public:
         resulttype alphar = 0.0;
         auto N = molefracs.size();
         for (auto i = 0; i < N; ++i) {
+            if (molefracs[i] == 0.0) { continue; }
             alphar = alphar + molefracs[i] * EOSs[i].alphar(tau, delta);
         }
         return alphar;