diff --git a/include/amici/forwardproblem.h b/include/amici/forwardproblem.h
index 2871ab28f8..9b5252c243 100644
--- a/include/amici/forwardproblem.h
+++ b/include/amici/forwardproblem.h
@@ -7,7 +7,6 @@
 #include "amici/vector.h"
 #include <amici/amici.h>
 
-#include <memory>
 #include <sundials/sundials_direct.h>
 #include <vector>
 
diff --git a/src/forwardproblem.cpp b/src/forwardproblem.cpp
index e3dcecd1a9..67ccd455f3 100644
--- a/src/forwardproblem.cpp
+++ b/src/forwardproblem.cpp
@@ -130,6 +130,8 @@ void ForwardProblem::workForwardProblem() {
 
     /* store initial state and sensitivity*/
     initial_state_ = getSimulationState();
+    // store root information at t0
+    model->froot(t_, x_, dx_, rootvals_);
 
     // get list of trigger timepoints for fixed-time triggered events
     auto trigger_timepoints = model->get_trigger_timepoints();