diff --git a/Project.toml b/Project.toml
index d21d44c..7cf0447 100644
--- a/Project.toml
+++ b/Project.toml
@@ -1,7 +1,7 @@
 name = "MPSKitModels"
 uuid = "ca635005-6f8c-4cd1-b51d-8491250ef2ab"
+version = "0.4.5"
 authors = ["Maarten Van Damme", "Lukas Devos", "Gertian Roose", "Klaas Gunst"]
-version = "0.4.4"
 
 [deps]
 LinearAlgebra = "37e2e46d-f89d-539d-b4ee-838fcccc9c8e"
@@ -17,7 +17,7 @@ MPSKit = "0.13.1"
 MacroTools = "0.5"
 PrecompileTools = "1"
 QuadGK = "2.11.1"
-TensorKit = "0.13,0.14, 0.15"
+TensorKit = "0.15, 0.16"
 TensorOperations = "5"
 TupleTools = "1"
 julia = "1.10"
diff --git a/test/hubbardoperators.jl b/test/hubbardoperators.jl
index ea67ebb..4072b40 100644
--- a/test/hubbardoperators.jl
+++ b/test/hubbardoperators.jl
@@ -72,7 +72,7 @@ end
     mu = randn()
 
     H_triv = hubbard_hamiltonian(Trivial, Trivial; t, U, mu, L)
-    vals_triv = mapreduce(vcat, eigvals(H_triv)) do (c, v)
+    vals_triv = mapreduce(vcat, pairs(eigvals(H_triv))) do (c, v)
         return repeat(real.(v), dim(c))
     end
     sort!(vals_triv)
@@ -82,7 +82,7 @@ end
             continue
         end
         H_symm = hubbard_hamiltonian(particle_symmetry, spin_symmetry; t, U, mu, L)
-        vals_symm = mapreduce(vcat, eigvals(H_symm)) do (c, v)
+        vals_symm = mapreduce(vcat, pairs(eigvals(H_symm))) do (c, v)
             return repeat(real.(v), dim(c))
         end
         sort!(vals_symm)
diff --git a/test/tjoperators.jl b/test/tjoperators.jl
index a092e1e..2268c0a 100644
--- a/test/tjoperators.jl
+++ b/test/tjoperators.jl
@@ -105,7 +105,7 @@ end
 
     for slave_fermion in (false, true)
         H_triv = tjhamiltonian(Trivial, Trivial; t, J, mu, L, slave_fermion)
-        vals_triv = mapreduce(vcat, eigvals(H_triv)) do (c, v)
+        vals_triv = mapreduce(vcat, pairs(eigvals(H_triv))) do (c, v)
             return repeat(real.(v), dim(c))
         end
         sort!(vals_triv)
@@ -117,7 +117,7 @@ end
             H_symm = tjhamiltonian(
                 particle_symmetry, spin_symmetry; t, J, mu, L, slave_fermion
             )
-            vals_symm = mapreduce(vcat, eigvals(H_symm)) do (c, v)
+            vals_symm = mapreduce(vcat, pairs(eigvals(H_symm))) do (c, v)
                 return repeat(real.(v), dim(c))
             end
             sort!(vals_symm)