diff --git a/R/run.R b/R/run.R
index 66756f8..ce4813e 100644
--- a/R/run.R
+++ b/R/run.R
@@ -16,6 +16,9 @@
 #' Cannot be greater than `n_runs`.
 #' @param verbose (\code{logical}) Whether to print a run summary on completion.
 #' Default is `TRUE`.
+#' @param timeout (\code{integer}) control the maximum time (in second) of a command execution.
+#' Can be use for certain inference methods (e.g SMC-APF), will bias the result.
+#' Default is 0 (no timeout).
 #' @param ... See [treepplr::tp_runtime_options()] for all supported arguments.
 #'
 #' @details
@@ -64,6 +67,7 @@ tp_run <- function(
   n_runs = 1,
   n_processes = 3,
   verbose = TRUE,
+  timeout = 0,
   ...
 ) {
   # start time
@@ -104,7 +108,7 @@ tp_run <- function(
     sampler$exe_path,
     data,
     options_to_string(tpplc_options[["runtime"]]),
-    paste(">", output_path)
+    ">"
   )
 
   if (n_runs > 1) {
@@ -117,18 +121,11 @@ tp_run <- function(
       1:n_runs,
       FUN = function(i) {
         output_path_i <- file.path(dir_path, paste0(out_file_name, i, ".json"))
-        command_i <- paste(
-          "LD_LIBRARY_PATH= ",
-          sampler$exe_path,
-          data,
-          options_to_string(tpplc_options[["runtime"]]),
-          paste(">", output_path_i)
-        )
-        system(command_i)
+        system(paste(command, output_path_i), timeout = timeout)
       }
     )
   } else {
-    system(command)
+    system(paste(command, output_path), timeout = timeout)
   }
 
   # the output (JSON) files
diff --git a/man/tp_run.Rd b/man/tp_run.Rd
index 84436c4..6bbd28a 100644
--- a/man/tp_run.Rd
+++ b/man/tp_run.Rd
@@ -12,6 +12,7 @@ tp_run(
   n_runs = 1,
   n_processes = 3,
   verbose = TRUE,
+  timeout = 0,
   ...
 )
 }
@@ -34,6 +35,10 @@ Cannot be greater than \code{n_runs}.}
 \item{verbose}{(\code{logical}) Whether to print a run summary on completion.
 Default is \code{TRUE}.}
 
+\item{timeout}{(\code{integer}) control the maximum time (in second) of a command execution.
+Can be use for certain inference methods (e.g SMC-APF), will bias the result.
+Default is 0 (no timeout).}
+
 \item{...}{See \code{\link[=tp_runtime_options]{tp_runtime_options()}} for all supported arguments.}
 }
 \value{
diff --git a/tests/testthat/test-run.R b/tests/testthat/test-run.R
index ccde83c..8483b65 100644
--- a/tests/testthat/test-run.R
+++ b/tests/testthat/test-run.R
@@ -59,3 +59,13 @@ test_that("Test-run_1e : tp_run no_parser", {
   )
   expect_true(is.character(run_smc))
 })
+
+test_that("Test-run_1f : tp_run with timeout", {
+  cat("\tTest-run_1f : tp_run with timeout \n")
+  expect_error(suppressWarnings(tp_run(
+    sampler = tp_compile(model = "crbd", method = "smc-apf", sweeps = 10, particles = 10000),
+    data = tp_data(data_input = "crbd"),
+    timeout = 1,
+    verbose = FALSE
+  )))
+})