diff --git a/benchmarks/OptimizationFrameworks/Project.toml b/benchmarks/OptimizationFrameworks/Project.toml
index 8967fc3f4..59a7e45cf 100644
--- a/benchmarks/OptimizationFrameworks/Project.toml
+++ b/benchmarks/OptimizationFrameworks/Project.toml
@@ -38,7 +38,7 @@ Enzyme = "0.11, 0.12, 0.13"
 ForwardDiff = "0.10"
 Ipopt = "1"
 JuMP = "1.15"
-ModelingToolkit = "9"
+ModelingToolkit = "9.47"
 NLPModelsIpopt = "0.10"
 Nonconvex = "2"
 NonconvexIpopt = "0.4"
diff --git a/benchmarks/OptimizationFrameworks/optimal_powerflow.jmd b/benchmarks/OptimizationFrameworks/optimal_powerflow.jmd
index b0cc68bfb..6249741cf 100644
--- a/benchmarks/OptimizationFrameworks/optimal_powerflow.jmd
+++ b/benchmarks/OptimizationFrameworks/optimal_powerflow.jmd
@@ -722,7 +722,7 @@ function build_opf_mtk_prob(dataset)
     end
 
     optsys = ModelingToolkit.OptimizationSystem(loss, vars, [], constraints=cons, name=:rosetta)
-    optsys = ModelingToolkit.complete(optsys)
+    optsys = ModelingToolkit.structural_simplify(optsys)
     u0map = [Num(k) => 0.0 for k in collect(unknowns(optsys))]
     ks = collect(Num.(unknowns(optsys)))
     for key in keys(ref[:bus])