diff --git a/Project.toml b/Project.toml
index 1f5e4234..bc8e3e0d 100644
--- a/Project.toml
+++ b/Project.toml
@@ -25,7 +25,6 @@ NamedGraphs = "678767b0-92e7-4007-89e4-4527a8725b19"
 PackageExtensionCompat = "65ce6f38-6b18-4e1d-a461-8949797d7930"
 Printf = "de0858da-6303-5e67-8744-51eddeeeb8d7"
 Random = "9a3f8284-a2c9-5f02-9a11-845980a1fd5c"
-Revise = "295af30f-e4ad-537b-8983-00126c2a3abe"
 SerializedElementArrays = "d3ce8812-9567-47e9-a7b5-65a6d70a3065"
 SimpleTraits = "699a6c99-e7fa-54fc-8d76-47d257e15c1d"
 SparseArrayKit = "a9a3c162-d163-4c15-8926-b8794fbefed2"
diff --git a/src/lib/ModelHamiltonians/src/ModelHamiltonians.jl b/src/lib/ModelHamiltonians/src/ModelHamiltonians.jl
index 53bf6adf..f54faea3 100644
--- a/src/lib/ModelHamiltonians/src/ModelHamiltonians.jl
+++ b/src/lib/ModelHamiltonians/src/ModelHamiltonians.jl
@@ -108,8 +108,8 @@ heisenberg(N::Integer; kwargs...) = heisenberg(path_graph(N); kwargs...)
 """
 Next-to-nearest-neighbor Ising model (ZZX) on a general graph
 """
-function ising(g::AbstractGraph; J1=-1, J2=0, h=0, hl=0)
-  (; J1, J2, h, hl) = map(to_callable, (; J1, J2, h, hl))
+function ising(g::AbstractGraph; J1=-1, J2=0, h=0)
+  (; J1, J2, h) = map(to_callable, (; J1, J2, h))
   ℋ = OpSum()
   for e in edges(g)
     ℋ += J1(e), "Sz", src(e), "Sz", dst(e)
@@ -129,7 +129,6 @@ function ising(g::AbstractGraph; J1=-1, J2=0, h=0, hl=0)
   end
   for v in vertices(g)
     ℋ += h(v), "Sx", v
-    ℋ += hl(v), "Sz", v
   end
   return ℋ
 end
diff --git a/src/solvers/dmrg.jl b/src/solvers/dmrg.jl
index f9718c0c..19511505 100644
--- a/src/solvers/dmrg.jl
+++ b/src/solvers/dmrg.jl
@@ -20,7 +20,7 @@ function ITensorMPS.dmrg(
   state = alternating_update(
     operator, init_state; nsweeps, nsites, updater, region_observer!, kwargs...
   )
-  eigval = first(eigvals_ref[])
+  eigval = only(eigvals_ref[])
   return eigval, state
 end
 
diff --git a/src/solvers/insert/insert.jl b/src/solvers/insert/insert.jl
index 29df13e6..11aed223 100644
--- a/src/solvers/insert/insert.jl
+++ b/src/solvers/insert/insert.jl
@@ -23,7 +23,7 @@ function default_inserter(
     indsTe = inds(state[ortho_vert])
     L, phi, spec = factorize(phi, indsTe; tags=tags(state, e), maxdim, mindim, cutoff)
     state[ortho_vert] = L
-    
+
   else
     v = ortho_vert
   end