diff --git a/include/teqp/algorithms/critical_tracing.hpp b/include/teqp/algorithms/critical_tracing.hpp
index b3b8a8a4585debbd875f344f5c7fd6abe42b0cec..2d8c96d5baee82c782f58b60b4aa48298e367d96 100644
--- a/include/teqp/algorithms/critical_tracing.hpp
+++ b/include/teqp/algorithms/critical_tracing.hpp
@@ -233,8 +233,8 @@ struct CriticalTracing {
 
         // The temperature derivative of total Psi w.r.t.T from a centered finite difference in T
         auto dT = 1e-7;
-        auto plusT = get_derivs(model, T + dT, rhovec).tot;
-        auto minusT = get_derivs(model, T - dT, rhovec).tot;
+        auto plusT = get_derivs(model, T + dT, rhovec, all_derivs.ei.v0).tot;
+        auto minusT = get_derivs(model, T - dT, rhovec, all_derivs.ei.v0).tot;
         auto derivT = (plusT - minusT) / (2.0 * dT);
 
         // Solve the eigenvalue problem for the given T & rho