diff --git a/Project.toml b/Project.toml
index 20c97db..b33f05e 100644
--- a/Project.toml
+++ b/Project.toml
@@ -1,7 +1,7 @@
 name = "ModelOrderReduction"
 uuid = "207801d6-6cee-43a9-ad0c-f0c64933efa0"
 authors = ["Bowen S. Zhu <bowenzhu@mit.edu> and contributors"]
-version = "0.1.1"
+version = "0.1.2"
 
 [deps]
 DocStringExtensions = "ffbed154-4ef7-542d-bbb7-c09d3a79fcae"
@@ -15,7 +15,7 @@ TSVD = "9449cd9e-2762-5aa3-a617-5413e99d722e"
 [compat]
 DocStringExtensions = "0.8, 0.9"
 LinearAlgebra = "1"
-ModelingToolkit = "9"
+ModelingToolkit = "10"
 RandomizedLinAlg = "0.1"
 Setfield = "0.8, 1"
 SparseArrays = "1"
diff --git a/docs/Project.toml b/docs/Project.toml
index 702b756..2370b58 100644
--- a/docs/Project.toml
+++ b/docs/Project.toml
@@ -14,5 +14,5 @@ Documenter = "1"
 LaTeXStrings = "1.3"
 MethodOfLines = "0.11"
 ModelOrderReduction = "0.1"
-ModelingToolkit = "9"
+ModelingToolkit = "10"
 Plots = "1.40"
diff --git a/src/deim.jl b/src/deim.jl
index 9334413..a2e72eb 100644
--- a/src/deim.jl
+++ b/src/deim.jl
@@ -112,7 +112,7 @@ The LHS of equations in `sys` are all assumed to be 1st order derivatives. Use
 `ModelingToolkit.ode_order_lowering` to transform higher order ODEs before applying DEIM.
 
 `sys` is assumed to have no internal systems. End users are encouraged to call
-`ModelingToolkit.structural_simplify` beforehand.
+`ModelingToolkit.mtkcompile` beforehand.
 
 The POD basis used for DEIM interpolation is obtained from the snapshot matrix of the
 nonlinear terms, which is computed by executing the runtime-generated function for
diff --git a/src/utils.jl b/src/utils.jl
index 24d8ef5..2e6bfca 100644
--- a/src/utils.jl
+++ b/src/utils.jl
@@ -7,7 +7,7 @@ $(TYPEDSIGNATURES)
 Return the differential equations and other non-differential equations.
 
 For differential equations, this function assume the derivative is the only single term at
-the left hand side, which is typically the result of `ModelingToolkit.structural_simplify`.
+the left hand side, which is typically the result of `ModelingToolkit.mtkcompile`.
 
 Equations from subsystems are not included.
 """
diff --git a/test/Project.toml b/test/Project.toml
index bfc7f97..5e29555 100644
--- a/test/Project.toml
+++ b/test/Project.toml
@@ -10,5 +10,5 @@ Test = "8dfed614-e22c-5e08-85e1-65c5234f0b40"
 Aqua = "0.8"
 OrdinaryDiffEq = "6"
 MethodOfLines = "0.11"
-ModelingToolkit = "9"
+ModelingToolkit = "10"
 SafeTestsets = "0.1"
diff --git a/test/deim.jl b/test/deim.jl
index 13f0d33..365e59b 100644
--- a/test/deim.jl
+++ b/test/deim.jl
@@ -31,7 +31,7 @@ dxs = [x => dx]
 order = 2
 discretization = MOLFiniteDifference(dxs, t; approx_order = order)
 ode_sys, tspan = symbolic_discretize(pde_sys, discretization)
-simp_sys = structural_simplify(ode_sys) # field substitutions is non-empty
+simp_sys = mtkcompile(ode_sys) # field substitutions is non-empty
 ode_prob = ODEProblem(simp_sys, nothing, tspan)
 sol = solve(ode_prob, Tsit5(), saveat = 1.0)