diff --git a/docs/src/examples/classic2d/1.hard-hexagon/index.md b/docs/src/examples/classic2d/1.hard-hexagon/index.md index 2981f1758..e94b03bb6 100644 --- a/docs/src/examples/classic2d/1.hard-hexagon/index.md +++ b/docs/src/examples/classic2d/1.hard-hexagon/index.md @@ -58,7 +58,7 @@ println("F = $F\tS = $S\tξ = $ξ") ```` ```` -F = 0.8839037051703845 S = 1.2807829622316382 ξ = 13.849682584604567 +F = 0.8839037051703852 S = 1.2807829622066444 ξ = 13.849682583029384 ```` @@ -98,7 +98,7 @@ c = f.coeffs[2] ```` ```` -0.8025265211777018 +0.8025249301530266 ```` ````julia @@ -111,57 +111,57 @@ plot!(p, ξ -> f(ξ) / 6; label="fit") - + - + - + - + - - + + - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + ``` diff --git a/docs/src/examples/classic2d/1.hard-hexagon/main.ipynb b/docs/src/examples/classic2d/1.hard-hexagon/main.ipynb index 4639509e0..fd7e7be84 100644 --- a/docs/src/examples/classic2d/1.hard-hexagon/main.ipynb +++ b/docs/src/examples/classic2d/1.hard-hexagon/main.ipynb @@ -145,11 +145,11 @@ "file_extension": ".jl", "mimetype": "application/julia", "name": "julia", - "version": "1.11.4" + "version": "1.10.4" }, "kernelspec": { - "name": "julia-1.11", - "display_name": "Julia 1.11.4", + "name": "julia-1.10", + "display_name": "Julia 1.10.4", "language": "julia" } }, diff --git a/docs/src/examples/quantum1d/1.ising-cft/index.md b/docs/src/examples/quantum1d/1.ising-cft/index.md index be3efa005..cd40011ed 100644 --- a/docs/src/examples/quantum1d/1.ising-cft/index.md +++ b/docs/src/examples/quantum1d/1.ising-cft/index.md @@ -27,18 +27,18 @@ H = periodic_boundary_conditions(transverse_field_ising(), L) ```` 12-site FiniteMPOHamiltonian{MPSKit.JordanMPOTensor{ComplexF64, TensorKit.ComplexSpace, Union{TensorKit.BraidingTensor{ComplexF64, TensorKit.ComplexSpace}, TensorKit.TensorMap{ComplexF64, TensorKit.ComplexSpace, 2, 2, Vector{ComplexF64}}}, TensorKit.TensorMap{ComplexF64, TensorKit.ComplexSpace, 2, 1, Vector{ComplexF64}}, TensorKit.TensorMap{ComplexF64, TensorKit.ComplexSpace, 1, 2, Vector{ComplexF64}}, TensorKit.TensorMap{ComplexF64, TensorKit.ComplexSpace, 1, 1, Vector{ComplexF64}}}}: -┬ W[12]: 6×1×1×1 JordanMPOTensor(((ℂ^1 ⊕ ℂ^1 ⊕ ℂ^1 ⊕ ℂ^1 ⊕ ℂ^1 ⊕ ℂ^1) ⊗ ⊕(ℂ^2)) ← (⊕(ℂ^2) ⊗ ⊕(ℂ^1))) -┼ W[11]: 6×1×1×6 JordanMPOTensor(((ℂ^1 ⊕ ℂ^1 ⊕ ℂ^1 ⊕ ℂ^1 ⊕ ℂ^1 ⊕ ℂ^1) ⊗ ⊕(ℂ^2)) ← (⊕(ℂ^2) ⊗ (ℂ^1 ⊕ ℂ^1 ⊕ ℂ^1 ⊕ ℂ^1 ⊕ ℂ^1 ⊕ ℂ^1))) -┼ W[10]: 6×1×1×6 JordanMPOTensor(((ℂ^1 ⊕ ℂ^1 ⊕ ℂ^1 ⊕ ℂ^1 ⊕ ℂ^1 ⊕ ℂ^1) ⊗ ⊕(ℂ^2)) ← (⊕(ℂ^2) ⊗ (ℂ^1 ⊕ ℂ^1 ⊕ ℂ^1 ⊕ ℂ^1 ⊕ ℂ^1 ⊕ ℂ^1))) -┼ W[9]: 6×1×1×6 JordanMPOTensor(((ℂ^1 ⊕ ℂ^1 ⊕ ℂ^1 ⊕ ℂ^1 ⊕ ℂ^1 ⊕ ℂ^1) ⊗ ⊕(ℂ^2)) ← (⊕(ℂ^2) ⊗ (ℂ^1 ⊕ ℂ^1 ⊕ ℂ^1 ⊕ ℂ^1 ⊕ ℂ^1 ⊕ ℂ^1))) -┼ W[8]: 6×1×1×6 JordanMPOTensor(((ℂ^1 ⊕ ℂ^1 ⊕ ℂ^1 ⊕ ℂ^1 ⊕ ℂ^1 ⊕ ℂ^1) ⊗ ⊕(ℂ^2)) ← (⊕(ℂ^2) ⊗ (ℂ^1 ⊕ ℂ^1 ⊕ ℂ^1 ⊕ ℂ^1 ⊕ ℂ^1 ⊕ ℂ^1))) -┼ W[7]: 6×1×1×6 JordanMPOTensor(((ℂ^1 ⊕ ℂ^1 ⊕ ℂ^1 ⊕ ℂ^1 ⊕ ℂ^1 ⊕ ℂ^1) ⊗ ⊕(ℂ^2)) ← (⊕(ℂ^2) ⊗ (ℂ^1 ⊕ ℂ^1 ⊕ ℂ^1 ⊕ ℂ^1 ⊕ ℂ^1 ⊕ ℂ^1))) -┼ W[6]: 6×1×1×6 JordanMPOTensor(((ℂ^1 ⊕ ℂ^1 ⊕ ℂ^1 ⊕ ℂ^1 ⊕ ℂ^1 ⊕ ℂ^1) ⊗ ⊕(ℂ^2)) ← (⊕(ℂ^2) ⊗ (ℂ^1 ⊕ ℂ^1 ⊕ ℂ^1 ⊕ ℂ^1 ⊕ ℂ^1 ⊕ ℂ^1))) -┼ W[5]: 6×1×1×6 JordanMPOTensor(((ℂ^1 ⊕ ℂ^1 ⊕ ℂ^1 ⊕ ℂ^1 ⊕ ℂ^1 ⊕ ℂ^1) ⊗ ⊕(ℂ^2)) ← (⊕(ℂ^2) ⊗ (ℂ^1 ⊕ ℂ^1 ⊕ ℂ^1 ⊕ ℂ^1 ⊕ ℂ^1 ⊕ ℂ^1))) -┼ W[4]: 6×1×1×6 JordanMPOTensor(((ℂ^1 ⊕ ℂ^1 ⊕ ℂ^1 ⊕ ℂ^1 ⊕ ℂ^1 ⊕ ℂ^1) ⊗ ⊕(ℂ^2)) ← (⊕(ℂ^2) ⊗ (ℂ^1 ⊕ ℂ^1 ⊕ ℂ^1 ⊕ ℂ^1 ⊕ ℂ^1 ⊕ ℂ^1))) -┼ W[3]: 6×1×1×6 JordanMPOTensor(((ℂ^1 ⊕ ℂ^1 ⊕ ℂ^1 ⊕ ℂ^1 ⊕ ℂ^1 ⊕ ℂ^1) ⊗ ⊕(ℂ^2)) ← (⊕(ℂ^2) ⊗ (ℂ^1 ⊕ ℂ^1 ⊕ ℂ^1 ⊕ ℂ^1 ⊕ ℂ^1 ⊕ ℂ^1))) -┼ W[2]: 6×1×1×6 JordanMPOTensor(((ℂ^1 ⊕ ℂ^1 ⊕ ℂ^1 ⊕ ℂ^1 ⊕ ℂ^1 ⊕ ℂ^1) ⊗ ⊕(ℂ^2)) ← (⊕(ℂ^2) ⊗ (ℂ^1 ⊕ ℂ^1 ⊕ ℂ^1 ⊕ ℂ^1 ⊕ ℂ^1 ⊕ ℂ^1))) -┴ W[1]: 1×1×1×6 JordanMPOTensor((⊕(ℂ^1) ⊗ ⊕(ℂ^2)) ← (⊕(ℂ^2) ⊗ (ℂ^1 ⊕ ℂ^1 ⊕ ℂ^1 ⊕ ℂ^1 ⊕ ℂ^1 ⊕ ℂ^1))) +┬ W[12]: 6×1×1×1 JordanMPOTensor(((ℂ^1 ⊕ ℂ^1 ⊕ ℂ^1 ⊕ ⋯ ⊕ ℂ^1 ⊗ ℂ^1 ⊗ ℂ^1) ⊗ ⊕(ℂ^2)) ← (⊕(ℂ^2) ⊗ ⊕(ℂ^1))) +┼ W[11]: 6×1×1×6 JordanMPOTensor(((ℂ^1 ⊕ ℂ^1 ⊕ ℂ^1 ⊕ ⋯ ⊕ ℂ^1 ⊗ ℂ^1 ⊗ ℂ^1) ⊗ ⊕(ℂ^2)) ← (⊕(ℂ^2) ⊗ (ℂ^1 ⊕ ℂ^1 ⊕ ℂ^1 ⊕ ⋯ ⊕ ℂ^1 ⊗ ℂ^1 ⊗ ℂ^1))) +┼ W[10]: 6×1×1×6 JordanMPOTensor(((ℂ^1 ⊕ ℂ^1 ⊕ ℂ^1 ⊕ ⋯ ⊕ ℂ^1 ⊗ ℂ^1 ⊗ ℂ^1) ⊗ ⊕(ℂ^2)) ← (⊕(ℂ^2) ⊗ (ℂ^1 ⊕ ℂ^1 ⊕ ℂ^1 ⊕ ⋯ ⊕ ℂ^1 ⊗ ℂ^1 ⊗ ℂ^1))) +┼ W[9]: 6×1×1×6 JordanMPOTensor(((ℂ^1 ⊕ ℂ^1 ⊕ ℂ^1 ⊕ ⋯ ⊕ ℂ^1 ⊗ ℂ^1 ⊗ ℂ^1) ⊗ ⊕(ℂ^2)) ← (⊕(ℂ^2) ⊗ (ℂ^1 ⊕ ℂ^1 ⊕ ℂ^1 ⊕ ⋯ ⊕ ℂ^1 ⊗ ℂ^1 ⊗ ℂ^1))) +┼ W[8]: 6×1×1×6 JordanMPOTensor(((ℂ^1 ⊕ ℂ^1 ⊕ ℂ^1 ⊕ ⋯ ⊕ ℂ^1 ⊗ ℂ^1 ⊗ ℂ^1) ⊗ ⊕(ℂ^2)) ← (⊕(ℂ^2) ⊗ (ℂ^1 ⊕ ℂ^1 ⊕ ℂ^1 ⊕ ⋯ ⊕ ℂ^1 ⊗ ℂ^1 ⊗ ℂ^1))) +┼ W[7]: 6×1×1×6 JordanMPOTensor(((ℂ^1 ⊕ ℂ^1 ⊕ ℂ^1 ⊕ ⋯ ⊕ ℂ^1 ⊗ ℂ^1 ⊗ ℂ^1) ⊗ ⊕(ℂ^2)) ← (⊕(ℂ^2) ⊗ (ℂ^1 ⊕ ℂ^1 ⊕ ℂ^1 ⊕ ⋯ ⊕ ℂ^1 ⊗ ℂ^1 ⊗ ℂ^1))) +┼ W[6]: 6×1×1×6 JordanMPOTensor(((ℂ^1 ⊕ ℂ^1 ⊕ ℂ^1 ⊕ ⋯ ⊕ ℂ^1 ⊗ ℂ^1 ⊗ ℂ^1) ⊗ ⊕(ℂ^2)) ← (⊕(ℂ^2) ⊗ (ℂ^1 ⊕ ℂ^1 ⊕ ℂ^1 ⊕ ⋯ ⊕ ℂ^1 ⊗ ℂ^1 ⊗ ℂ^1))) +┼ W[5]: 6×1×1×6 JordanMPOTensor(((ℂ^1 ⊕ ℂ^1 ⊕ ℂ^1 ⊕ ⋯ ⊕ ℂ^1 ⊗ ℂ^1 ⊗ ℂ^1) ⊗ ⊕(ℂ^2)) ← (⊕(ℂ^2) ⊗ (ℂ^1 ⊕ ℂ^1 ⊕ ℂ^1 ⊕ ⋯ ⊕ ℂ^1 ⊗ ℂ^1 ⊗ ℂ^1))) +┼ W[4]: 6×1×1×6 JordanMPOTensor(((ℂ^1 ⊕ ℂ^1 ⊕ ℂ^1 ⊕ ⋯ ⊕ ℂ^1 ⊗ ℂ^1 ⊗ ℂ^1) ⊗ ⊕(ℂ^2)) ← (⊕(ℂ^2) ⊗ (ℂ^1 ⊕ ℂ^1 ⊕ ℂ^1 ⊕ ⋯ ⊕ ℂ^1 ⊗ ℂ^1 ⊗ ℂ^1))) +┼ W[3]: 6×1×1×6 JordanMPOTensor(((ℂ^1 ⊕ ℂ^1 ⊕ ℂ^1 ⊕ ⋯ ⊕ ℂ^1 ⊗ ℂ^1 ⊗ ℂ^1) ⊗ ⊕(ℂ^2)) ← (⊕(ℂ^2) ⊗ (ℂ^1 ⊕ ℂ^1 ⊕ ℂ^1 ⊕ ⋯ ⊕ ℂ^1 ⊗ ℂ^1 ⊗ ℂ^1))) +┼ W[2]: 6×1×1×6 JordanMPOTensor(((ℂ^1 ⊕ ℂ^1 ⊕ ℂ^1 ⊕ ⋯ ⊕ ℂ^1 ⊗ ℂ^1 ⊗ ℂ^1) ⊗ ⊕(ℂ^2)) ← (⊕(ℂ^2) ⊗ (ℂ^1 ⊕ ℂ^1 ⊕ ℂ^1 ⊕ ⋯ ⊕ ℂ^1 ⊗ ℂ^1 ⊗ ℂ^1))) +┴ W[1]: 1×1×1×6 JordanMPOTensor((⊕(ℂ^1) ⊗ ⊕(ℂ^2)) ← (⊕(ℂ^2) ⊗ (ℂ^1 ⊕ ℂ^1 ⊕ ℂ^1 ⊕ ⋯ ⊕ ℂ^1 ⊗ ℂ^1 ⊗ ℂ^1))) ```` @@ -63,69 +63,67 @@ plot(real.(energies); - + - + - + - + - - + + - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + ``` @@ -146,8 +144,8 @@ or in the code as: ````julia function O_shift(L) - τ = BraidingTensor{ComplexF64}(ℂ^2, ℂ^2) - O = TensorMap(τ) + I = id(ComplexF64, ℂ^2) + @tensor O[W S; N E] := I[W; N] * I[S; E] return periodic_boundary_conditions(InfiniteMPO([O]), L) end ```` @@ -161,7 +159,8 @@ operator. However, there is a subtlety because of the degeneracies in the energy eigenvalues. The eigensolver will find an orthonormal basis within each energy subspace, but this basis is not necessarily a basis of eigenstates of the translation operator. In order to fix this, we diagonalize the translation operator within each energy subspace. -The resulting energy levels have one-to-one correspondence to the operators in CFT, where the momentum is related to their conformal spin as $P_n = \frac{2\pi}{L}S_n$. +The resulting energy levels have one-to-one correspondence to the operators in CFT, where +the momentum is related to their conformal spin as $P_n = \frac{2\pi}{L}S_n$. ````julia function fix_degeneracies(basis) @@ -190,34 +189,40 @@ append!(momenta, fix_degeneracies(states[17:18])) ```` 18-element Vector{Float64}: - 3.1857542530442574e-30 - 1.4051172205135508e-30 - -6.311400512248983e-31 - -1.110223024625158e-16 - -1.3602461870053588e-18 - 1.0408340855860866e-16 - -4.1470733097570545e-18 - 3.371727499927069e-18 - 1.5959455978986586e-16 - 7.055770862123108e-18 - 5.5511151231257815e-17 - -2.973502681954427e-31 - -6.331740687315362e-17 - 5.359415723744804e-18 - 1.301042606982604e-18 - -9.920449878242352e-17 - 0.0 - -6.247725059346555e-19 + -3.9672804251717156e-17 + -3.7809066376698395e-17 + 1.9364434176550726e-17 + -0.5235987755982991 + 0.5235987755982987 + 1.047197551196598 + -1.047197551196597 + -0.5235987755982985 + 0.5235987755982987 + -1.0471975511965976 + 1.0471975511965976 + -1.1073637563936691e-17 + -1.5707963267948968 + 1.5707963267948968 + -1.0471975511965979 + 1.0471975511965972 + 1.5707963267948968 + -1.5707963267948963 ```` -Calculating scaling dimensions from the energy gap +We can compute the scaling dimensions $\Delta_n$ of the operators in the CFT from the +energy gap of the corresponding excitations as $E_n - E_0 = \frac{2\pi v}{L} \Delta_n$, +where $v = 2$. If we plot these scaling dimensions against the conformal spin $S_n$ from +above, we retrieve the familiar spectrum of the Ising CFT. ````julia v = 2.0 Δ = real.(energies[1:18] .- energies[1]) ./ (2π * v / L) -p = plot(momenta, real.(Δ); - seriestype=:scatter, xlabel="momentum", ylabel="energy", legend=false) -vline!(p, [2π / L * i for i in -3:3]; color="gray", linestyle=:dash) +S = momenta ./ (2π / L) + +p = plot(S, real.(Δ); + seriestype=:scatter, xlabel="conformal spin (S)", ylabel="scaling dimension (Δ)", + legend=false) +vline!(p, -3:3; color="gray", linestyle=:dash) hline!(p, [0, 1 / 8, 1, 9 / 8, 2, 17 / 8]; color="gray", linestyle=:dash) p ```` @@ -226,81 +231,81 @@ p - + - + - + - + - - + + - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + ``` @@ -318,27 +323,36 @@ D = 64 ```` ```` -[ Info: DMRG init: obj = -1.971321475690e+01 err = 7.4003e-02 -[ Info: DMRG 1: obj = -2.549098964429e+01 err = 6.9440494815e-03 time = 2.67 sec -[ Info: DMRG 2: obj = -2.549098968636e+01 err = 9.0032288881e-07 time = 0.66 sec -[ Info: DMRG 3: obj = -2.549098968636e+01 err = 1.3353746489e-07 time = 0.57 sec -[ Info: DMRG 4: obj = -2.549098968636e+01 err = 1.0102919353e-08 time = 0.35 sec -[ Info: DMRG 5: obj = -2.549098968636e+01 err = 3.6147044323e-09 time = 0.40 sec -[ Info: DMRG 6: obj = -2.549098968636e+01 err = 2.3487622576e-09 time = 0.39 sec -[ Info: DMRG 7: obj = -2.549098968636e+01 err = 1.6417769963e-09 time = 0.47 sec -[ Info: DMRG 8: obj = -2.549098968636e+01 err = 1.1602119470e-09 time = 0.64 sec -[ Info: DMRG 9: obj = -2.549098968636e+01 err = 8.1567961851e-10 time = 0.39 sec -[ Info: DMRG 10: obj = -2.549098968636e+01 err = 5.7049959097e-10 time = 0.39 sec -[ Info: DMRG 11: obj = -2.549098968636e+01 err = 4.6429052737e-10 time = 0.40 sec -[ Info: DMRG 12: obj = -2.549098968636e+01 err = 3.8662080912e-10 time = 0.40 sec -[ Info: DMRG 13: obj = -2.549098968636e+01 err = 3.2276986394e-10 time = 0.40 sec -[ Info: DMRG 14: obj = -2.549098968636e+01 err = 2.7010926832e-10 time = 0.40 sec -[ Info: DMRG 15: obj = -2.549098968636e+01 err = 2.2653819206e-10 time = 0.38 sec -[ Info: DMRG 16: obj = -2.549098968636e+01 err = 1.9037561223e-10 time = 0.40 sec -[ Info: DMRG 17: obj = -2.549098968636e+01 err = 1.6027409710e-10 time = 0.49 sec -[ Info: DMRG 18: obj = -2.549098968636e+01 err = 1.3515030264e-10 time = 0.54 sec -[ Info: DMRG 19: obj = -2.549098968636e+01 err = 1.1412981266e-10 time = 0.39 sec -[ Info: DMRG conv 20: obj = -2.549098968636e+01 err = 9.6503774976e-11 time = 11.14 sec +[ Info: DMRG init: obj = -1.990218724629e+01 err = 7.4846e-02 +[ Info: DMRG 1: obj = -2.549098967874e+01 err = 4.1389989697e-03 time = 1.98 sec +[ Info: DMRG 2: obj = -2.549098968636e+01 err = 8.3286537771e-07 time = 0.43 sec +[ Info: DMRG 3: obj = -2.549098968636e+01 err = 7.8836932549e-08 time = 0.78 sec +[ Info: DMRG 4: obj = -2.549098968636e+01 err = 1.2157986090e-08 time = 0.52 sec +[ Info: DMRG 5: obj = -2.549098968636e+01 err = 4.5963842433e-09 time = 0.56 sec +[ Info: DMRG 6: obj = -2.549098968636e+01 err = 2.3229977031e-09 time = 0.55 sec +[ Info: DMRG 7: obj = -2.549098968636e+01 err = 1.9807447155e-09 time = 0.79 sec +[ Info: DMRG 8: obj = -2.549098968636e+01 err = 1.8101666690e-09 time = 0.70 sec +[ Info: DMRG 9: obj = -2.549098968636e+01 err = 1.6596660041e-09 time = 1.00 sec +[ Info: DMRG 10: obj = -2.549098968636e+01 err = 1.5028466015e-09 time = 1.00 sec +[ Info: DMRG 11: obj = -2.549098968636e+01 err = 1.3419670975e-09 time = 0.99 sec +[ Info: DMRG 12: obj = -2.549098968636e+01 err = 1.1842249707e-09 time = 1.00 sec +[ Info: DMRG 13: obj = -2.549098968636e+01 err = 1.0356798368e-09 time = 0.99 sec +[ Info: DMRG 14: obj = -2.549098968636e+01 err = 9.0001540544e-10 time = 0.99 sec +[ Info: DMRG 15: obj = -2.549098968636e+01 err = 7.7880124636e-10 time = 1.01 sec +[ Info: DMRG 16: obj = -2.549098968636e+01 err = 6.7212700978e-10 time = 1.00 sec +[ Info: DMRG 17: obj = -2.549098968636e+01 err = 5.7919723678e-10 time = 1.02 sec +[ Info: DMRG 18: obj = -2.549098968636e+01 err = 4.9877203200e-10 time = 1.03 sec +[ Info: DMRG 19: obj = -2.549098968636e+01 err = 4.2945107431e-10 time = 1.02 sec +[ Info: DMRG 20: obj = -2.549098968636e+01 err = 3.6984011575e-10 time = 1.03 sec +[ Info: DMRG 21: obj = -2.549098968636e+01 err = 3.1863801222e-10 time = 1.05 sec +[ Info: DMRG 22: obj = -2.549098968636e+01 err = 2.7467488610e-10 time = 1.06 sec +[ Info: DMRG 23: obj = -2.549098968636e+01 err = 2.3692226590e-10 time = 1.05 sec +[ Info: DMRG 24: obj = -2.549098968636e+01 err = 2.0448836836e-10 time = 1.06 sec +[ Info: DMRG 25: obj = -2.549098968636e+01 err = 1.7660632136e-10 time = 1.06 sec +[ Info: DMRG 26: obj = -2.549098968636e+01 err = 1.5262000256e-10 time = 1.05 sec +[ Info: DMRG 27: obj = -2.549098968636e+01 err = 1.3196928368e-10 time = 1.13 sec +[ Info: DMRG 28: obj = -2.549098968636e+01 err = 1.1417667953e-10 time = 1.35 sec +[ Info: DMRG conv 29: obj = -2.549098968636e+01 err = 9.8835132816e-11 time = 28.86 sec ```` @@ -347,121 +361,108 @@ ansatz. This returns quasiparticle states, which can be converted to regular `Fi objects. ````julia -E_ex, qps = excitations(H_mps, QuasiparticleAnsatz(), ψ, envs; num=16) +E_ex, qps = excitations(H_mps, QuasiparticleAnsatz(), ψ, envs; num=18) states_mps = vcat(ψ, map(qp -> convert(FiniteMPS, qp), qps)) -E_mps = map(x -> expectation_value(x, H_mps), states_mps) +energies_mps = map(x -> expectation_value(x, H_mps), states_mps) momenta_mps = Float64[] -append!(momenta_mps, fix_degeneracies(states[1:1])) -append!(momenta_mps, fix_degeneracies(states[2:2])) -append!(momenta_mps, fix_degeneracies(states[3:3])) -append!(momenta_mps, fix_degeneracies(states[4:5])) -append!(momenta_mps, fix_degeneracies(states[6:9])) -append!(momenta_mps, fix_degeneracies(states[10:11])) -append!(momenta_mps, fix_degeneracies(states[12:12])) -append!(momenta_mps, fix_degeneracies(states[13:16])) +append!(momenta_mps, fix_degeneracies(states_mps[1:1])) +append!(momenta_mps, fix_degeneracies(states_mps[2:2])) +append!(momenta_mps, fix_degeneracies(states_mps[3:3])) +append!(momenta_mps, fix_degeneracies(states_mps[4:5])) +append!(momenta_mps, fix_degeneracies(states_mps[6:9])) +append!(momenta_mps, fix_degeneracies(states_mps[10:11])) +append!(momenta_mps, fix_degeneracies(states_mps[12:12])) +append!(momenta_mps, fix_degeneracies(states_mps[13:16])) +append!(momenta_mps, fix_degeneracies(states_mps[17:18])) v = 2.0 -Δ = real.(energies[1:18] .- energies[1]) ./ (2π * v / L) -plot(momenta_mps, Δ; - seriestype=:scatter, xlabel="momentum", ylabel="energy", legend=false) -vline!(p, [2π / L * i for i in -3:3]; color="gray", linestyle=:dash) -hline!(p, [0, 1 / 8, 1, 9 / 8, 2, 17 / 8]; color="gray", linestyle=:dash) -p +Δ_mps = real.(energies_mps[1:18] .- energies_mps[1]) ./ (2π * v / L_mps) +S_mps = momenta_mps ./ (2π / L_mps) + +p_mps = plot(S_mps, real.(Δ_mps); + seriestype=:scatter, xlabel="conformal spin (S)", + ylabel="scaling dimension (Δ)", legend=false) +vline!(p_mps, -3:3; color="gray", linestyle=:dash) +hline!(p_mps, [0, 1 / 8, 1, 9 / 8, 2, 17 / 8]; color="gray", linestyle=:dash) +p_mps ```` ```@raw html - + - + - + - + - - + + - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + ``` diff --git a/docs/src/examples/quantum1d/1.ising-cft/main.ipynb b/docs/src/examples/quantum1d/1.ising-cft/main.ipynb index fec792848..11419eaee 100644 --- a/docs/src/examples/quantum1d/1.ising-cft/main.ipynb +++ b/docs/src/examples/quantum1d/1.ising-cft/main.ipynb @@ -96,8 +96,8 @@ "cell_type": "code", "source": [ "function O_shift(L)\n", - " τ = BraidingTensor{ComplexF64}(ℂ^2, ℂ^2)\n", - " O = TensorMap(τ)\n", + " I = id(ComplexF64, ℂ^2)\n", + " @tensor O[W S; N E] := I[W; N] * I[S; E]\n", " return periodic_boundary_conditions(InfiniteMPO([O]), L)\n", "end" ], @@ -112,7 +112,8 @@ "eigenvalues. The eigensolver will find an orthonormal basis within each energy subspace, but\n", "this basis is not necessarily a basis of eigenstates of the translation operator. In order\n", "to fix this, we diagonalize the translation operator within each energy subspace.\n", - "The resulting energy levels have one-to-one correspondence to the operators in CFT, where the momentum is related to their conformal spin as $P_n = \\frac{2\\pi}{L}S_n$." + "The resulting energy levels have one-to-one correspondence to the operators in CFT, where\n", + "the momentum is related to their conformal spin as $P_n = \\frac{2\\pi}{L}S_n$." ], "metadata": {} }, @@ -149,7 +150,10 @@ { "cell_type": "markdown", "source": [ - "Calculating scaling dimensions from the energy gap" + "We can compute the scaling dimensions $\\Delta_n$ of the operators in the CFT from the\n", + "energy gap of the corresponding excitations as $E_n - E_0 = \\frac{2\\pi v}{L} \\Delta_n$,\n", + "where $v = 2$. If we plot these scaling dimensions against the conformal spin $S_n$ from\n", + "above, we retrieve the familiar spectrum of the Ising CFT." ], "metadata": {} }, @@ -159,9 +163,12 @@ "source": [ "v = 2.0\n", "Δ = real.(energies[1:18] .- energies[1]) ./ (2π * v / L)\n", - "p = plot(momenta, real.(Δ);\n", - " seriestype=:scatter, xlabel=\"momentum\", ylabel=\"energy\", legend=false)\n", - "vline!(p, [2π / L * i for i in -3:3]; color=\"gray\", linestyle=:dash)\n", + "S = momenta ./ (2π / L)\n", + "\n", + "p = plot(S, real.(Δ);\n", + " seriestype=:scatter, xlabel=\"conformal spin (S)\", ylabel=\"scaling dimension (Δ)\",\n", + " legend=false)\n", + "vline!(p, -3:3; color=\"gray\", linestyle=:dash)\n", "hline!(p, [0, 1 / 8, 1, 9 / 8, 2, 17 / 8]; color=\"gray\", linestyle=:dash)\n", "p" ], @@ -203,27 +210,31 @@ "outputs": [], "cell_type": "code", "source": [ - "E_ex, qps = excitations(H_mps, QuasiparticleAnsatz(), ψ, envs; num=16)\n", + "E_ex, qps = excitations(H_mps, QuasiparticleAnsatz(), ψ, envs; num=18)\n", "states_mps = vcat(ψ, map(qp -> convert(FiniteMPS, qp), qps))\n", - "E_mps = map(x -> expectation_value(x, H_mps), states_mps)\n", + "energies_mps = map(x -> expectation_value(x, H_mps), states_mps)\n", "\n", "momenta_mps = Float64[]\n", - "append!(momenta_mps, fix_degeneracies(states[1:1]))\n", - "append!(momenta_mps, fix_degeneracies(states[2:2]))\n", - "append!(momenta_mps, fix_degeneracies(states[3:3]))\n", - "append!(momenta_mps, fix_degeneracies(states[4:5]))\n", - "append!(momenta_mps, fix_degeneracies(states[6:9]))\n", - "append!(momenta_mps, fix_degeneracies(states[10:11]))\n", - "append!(momenta_mps, fix_degeneracies(states[12:12]))\n", - "append!(momenta_mps, fix_degeneracies(states[13:16]))\n", + "append!(momenta_mps, fix_degeneracies(states_mps[1:1]))\n", + "append!(momenta_mps, fix_degeneracies(states_mps[2:2]))\n", + "append!(momenta_mps, fix_degeneracies(states_mps[3:3]))\n", + "append!(momenta_mps, fix_degeneracies(states_mps[4:5]))\n", + "append!(momenta_mps, fix_degeneracies(states_mps[6:9]))\n", + "append!(momenta_mps, fix_degeneracies(states_mps[10:11]))\n", + "append!(momenta_mps, fix_degeneracies(states_mps[12:12]))\n", + "append!(momenta_mps, fix_degeneracies(states_mps[13:16]))\n", + "append!(momenta_mps, fix_degeneracies(states_mps[17:18]))\n", "\n", "v = 2.0\n", - "Δ = real.(energies[1:18] .- energies[1]) ./ (2π * v / L)\n", - "plot(momenta_mps, Δ;\n", - " seriestype=:scatter, xlabel=\"momentum\", ylabel=\"energy\", legend=false)\n", - "vline!(p, [2π / L * i for i in -3:3]; color=\"gray\", linestyle=:dash)\n", - "hline!(p, [0, 1 / 8, 1, 9 / 8, 2, 17 / 8]; color=\"gray\", linestyle=:dash)\n", - "p" + "Δ_mps = real.(energies_mps[1:18] .- energies_mps[1]) ./ (2π * v / L_mps)\n", + "S_mps = momenta_mps ./ (2π / L_mps)\n", + "\n", + "p_mps = plot(S_mps, real.(Δ_mps);\n", + " seriestype=:scatter, xlabel=\"conformal spin (S)\",\n", + " ylabel=\"scaling dimension (Δ)\", legend=false)\n", + "vline!(p_mps, -3:3; color=\"gray\", linestyle=:dash)\n", + "hline!(p_mps, [0, 1 / 8, 1, 9 / 8, 2, 17 / 8]; color=\"gray\", linestyle=:dash)\n", + "p_mps" ], "metadata": {}, "execution_count": null @@ -244,11 +255,11 @@ "file_extension": ".jl", "mimetype": "application/julia", "name": "julia", - "version": "1.11.5" + "version": "1.10.4" }, "kernelspec": { - "name": "julia-1.11", - "display_name": "Julia 1.11.5", + "name": "julia-1.10", + "display_name": "Julia 1.10.4", "language": "julia" } }, diff --git a/docs/src/examples/quantum1d/2.haldane/index.md b/docs/src/examples/quantum1d/2.haldane/index.md index 52fe3c4fb..5a1d86df5 100644 --- a/docs/src/examples/quantum1d/2.haldane/index.md +++ b/docs/src/examples/quantum1d/2.haldane/index.md @@ -61,7 +61,7 @@ En_2, st_2 = excitations(H, QuasiparticleAnsatz(), ψ, envs; sector=SU2Irrep(2)) ```` ```` -0.7989253589480537 +0.7989253589480468 ```` We can go even further and doublecheck the claim that ``S = 1`` is an edge excitation, by plotting the energy density. @@ -83,51 +83,51 @@ plot!(p_density, e₂; label="S = 2") - + - + - + - + - - + + - - - - - - - - - - - - - - - - - - - - - - - - - - - - - + + + + + + + + + + + + + + + + + + + + + + + + + + + + + ``` @@ -150,7 +150,7 @@ f = fit(Ls .^ (-2), ΔEs, 1) ```` ```` -0.45173401585844547 +0.4517340158582977 ```` ````julia @@ -163,53 +163,53 @@ plot!(p_size_extrapolation, x -> f(x); label="fit") - + - + - + - + - - + + - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + ``` @@ -238,23 +238,25 @@ println("minimum @k = $(kspace[idx]):\t ΔE = $(ΔE)") ```` ```` -[ Info: Found excitations for momentum = 0.20943951023931953 +┌ Warning: resorting to η +└ @ OptimKit ~/.julia/packages/OptimKit/G6i79/src/cg.jl:207 [ Info: Found excitations for momentum = 0.0 +[ Info: Found excitations for momentum = 0.20943951023931953 [ Info: Found excitations for momentum = 0.41887902047863906 [ Info: Found excitations for momentum = 0.6283185307179586 +[ Info: Found excitations for momentum = 0.8377580409572781 [ Info: Found excitations for momentum = 1.0471975511965976 [ Info: Found excitations for momentum = 1.2566370614359172 [ Info: Found excitations for momentum = 1.4660765716752369 -[ Info: Found excitations for momentum = 0.8377580409572781 [ Info: Found excitations for momentum = 1.6755160819145563 [ Info: Found excitations for momentum = 1.8849555921538759 [ Info: Found excitations for momentum = 2.0943951023931953 [ Info: Found excitations for momentum = 2.303834612632515 -[ Info: Found excitations for momentum = 2.9321531433504737 -[ Info: Found excitations for momentum = 2.722713633111154 [ Info: Found excitations for momentum = 2.5132741228718345 +[ Info: Found excitations for momentum = 2.722713633111154 +[ Info: Found excitations for momentum = 2.9321531433504737 [ Info: Found excitations for momentum = 3.141592653589793 -minimum @k = 3.141592653589793: ΔE = 0.4104792487438299 +minimum @k = 3.141592653589793: ΔE = 0.41047924862708485 ```` @@ -266,47 +268,47 @@ plot(kspace, real.(Es); xaxis="momentum", yaxis="ΔE", label="S = 1") - + - + - + - + - - + + - - - - - - - - - - - - - - - - - - - - - - - - - + + + + + + + + + + + + + + + + + + + + + + + + + ``` diff --git a/docs/src/examples/quantum1d/2.haldane/main.ipynb b/docs/src/examples/quantum1d/2.haldane/main.ipynb index 0a2bc7798..a3ec1bc5c 100644 --- a/docs/src/examples/quantum1d/2.haldane/main.ipynb +++ b/docs/src/examples/quantum1d/2.haldane/main.ipynb @@ -187,11 +187,11 @@ "file_extension": ".jl", "mimetype": "application/julia", "name": "julia", - "version": "1.11.4" + "version": "1.10.4" }, "kernelspec": { - "name": "julia-1.11", - "display_name": "Julia 1.11.4", + "name": "julia-1.10", + "display_name": "Julia 1.10.4", "language": "julia" } }, diff --git a/docs/src/examples/quantum1d/3.ising-dqpt/index.md b/docs/src/examples/quantum1d/3.ising-dqpt/index.md index 1f62c4a88..f72756afe 100644 --- a/docs/src/examples/quantum1d/3.ising-dqpt/index.md +++ b/docs/src/examples/quantum1d/3.ising-dqpt/index.md @@ -39,13 +39,14 @@ H₀ = transverse_field_ising(FiniteChain(L); g=-0.5) ```` ```` -[ Info: DMRG init: obj = +9.627138678740e+00 err = 1.5287e-01 -[ Info: DMRG 1: obj = -2.040021714898e+01 err = 1.5582585194e-02 time = 0.10 sec -[ Info: DMRG 2: obj = -2.040021715175e+01 err = 3.4655395009e-07 time = 0.02 sec -[ Info: DMRG 3: obj = -2.040021763026e+01 err = 3.7385206468e-05 time = 0.06 sec -[ Info: DMRG 4: obj = -2.040021786702e+01 err = 2.1782648849e-06 time = 0.04 sec -[ Info: DMRG 5: obj = -2.040021786703e+01 err = 1.3998187485e-07 time = 0.02 sec -[ Info: DMRG conv 6: obj = -2.040021786703e+01 err = 9.0384912069e-11 time = 0.26 sec +[ Info: DMRG init: obj = +9.906929661608e+00 err = 1.4654e-01 +[ Info: DMRG 1: obj = -2.040021714938e+01 err = 9.3181641986e-04 time = 0.05 sec +[ Info: DMRG 2: obj = -2.040021715179e+01 err = 2.4688856530e-07 time = 0.03 sec +[ Info: DMRG 3: obj = -2.040021786221e+01 err = 4.2747944525e-05 time = 0.08 sec +[ Info: DMRG 4: obj = -2.040021786699e+01 err = 1.6446674043e-06 time = 0.04 sec +[ Info: DMRG 5: obj = -2.040021786703e+01 err = 2.4678293656e-07 time = 0.03 sec +[ Info: DMRG 6: obj = -2.040021786703e+01 err = 2.3749087526e-10 time = 0.03 sec +[ Info: DMRG conv 7: obj = -2.040021786703e+01 err = 4.3310784899e-12 time = 0.29 sec ```` @@ -111,14 +112,13 @@ H₀ = transverse_field_ising(; g=-0.5) ```` ```` -[ Info: VUMPS init: obj = +4.981979532800e-01 err = 3.9530e-01 -[ Info: VUMPS 1: obj = -1.030701827098e+00 err = 1.0340065851e-01 time = 6.72 sec -[ Info: VUMPS 2: obj = -1.063544326030e+00 err = 2.9155624606e-04 time = 0.01 sec -[ Info: VUMPS 3: obj = -1.063544409971e+00 err = 1.8767724186e-06 time = 0.01 sec -[ Info: VUMPS 4: obj = -1.063544409973e+00 err = 2.0959736412e-08 time = 0.00 sec -[ Info: VUMPS 5: obj = -1.063544409973e+00 err = 2.7005924256e-09 time = 0.00 sec -[ Info: VUMPS 6: obj = -1.063544409973e+00 err = 3.0434524435e-10 time = 0.00 sec -[ Info: VUMPS conv 7: obj = -1.063544409973e+00 err = 2.4566755222e-11 time = 6.75 sec +[ Info: VUMPS init: obj = +4.937592959715e-01 err = 3.8640e-01 +[ Info: VUMPS 1: obj = -1.041532472680e+00 err = 1.0207037612e-01 time = 7.91 sec +[ Info: VUMPS 2: obj = -1.063544395012e+00 err = 1.0570052184e-04 time = 0.01 sec +[ Info: VUMPS 3: obj = -1.063544409972e+00 err = 1.3451291600e-06 time = 0.01 sec +[ Info: VUMPS 4: obj = -1.063544409973e+00 err = 2.1982905922e-08 time = 0.00 sec +[ Info: VUMPS 5: obj = -1.063544409973e+00 err = 8.6218866571e-10 time = 0.00 sec +[ Info: VUMPS conv 6: obj = -1.063544409973e+00 err = 7.1345139950e-11 time = 7.94 sec ```` @@ -130,7 +130,7 @@ dot(ψ₀, ψ₀) ```` ```` -1.0000000000000053 - 8.519001927175606e-17im +0.9999999999999991 + 5.43344267830843e-16im ```` so the loschmidth echo takes on the pleasant form diff --git a/docs/src/examples/quantum1d/3.ising-dqpt/main.ipynb b/docs/src/examples/quantum1d/3.ising-dqpt/main.ipynb index 005397e62..983169527 100644 --- a/docs/src/examples/quantum1d/3.ising-dqpt/main.ipynb +++ b/docs/src/examples/quantum1d/3.ising-dqpt/main.ipynb @@ -285,11 +285,11 @@ "file_extension": ".jl", "mimetype": "application/julia", "name": "julia", - "version": "1.11.4" + "version": "1.10.4" }, "kernelspec": { - "name": "julia-1.11", - "display_name": "Julia 1.11.4", + "name": "julia-1.10", + "display_name": "Julia 1.10.4", "language": "julia" } }, diff --git a/docs/src/examples/quantum1d/4.xxz-heisenberg/index.md b/docs/src/examples/quantum1d/4.xxz-heisenberg/index.md index 2e43e0cc4..01541310e 100644 --- a/docs/src/examples/quantum1d/4.xxz-heisenberg/index.md +++ b/docs/src/examples/quantum1d/4.xxz-heisenberg/index.md @@ -55,208 +55,208 @@ groundstate, cache, delta = find_groundstate(state, H, VUMPS()); ```` ```` -[ Info: VUMPS init: obj = +2.499997151608e-01 err = 1.5076e-03 -[ Info: VUMPS 1: obj = -1.102936647693e-01 err = 3.8393870916e-01 time = 0.02 sec -[ Info: VUMPS 2: obj = -2.518819686637e-03 err = 3.8247870042e-01 time = 0.02 sec -[ Info: VUMPS 3: obj = -1.283698632877e-01 err = 4.5231528318e-01 time = 0.02 sec -[ Info: VUMPS 4: obj = -1.318247939270e-01 err = 3.9194798602e-01 time = 0.02 sec -[ Info: VUMPS 5: obj = -1.563893197526e-01 err = 3.6030425585e-01 time = 0.02 sec -[ Info: VUMPS 6: obj = -1.275381556576e-01 err = 3.7241063143e-01 time = 0.09 sec -[ Info: VUMPS 7: obj = -2.001728737220e-01 err = 3.5694925283e-01 time = 0.09 sec -[ Info: VUMPS 8: obj = -4.148514795033e-01 err = 1.9978091568e-01 time = 0.02 sec -[ Info: VUMPS 9: obj = -8.948544839910e-02 err = 3.8051752551e-01 time = 0.02 sec -[ Info: VUMPS 10: obj = +3.907102087904e-03 err = 4.0778106525e-01 time = 0.02 sec -[ Info: VUMPS 11: obj = -4.210510632748e-02 err = 3.9552473178e-01 time = 0.02 sec -[ Info: VUMPS 12: obj = -1.143304888714e-01 err = 3.8847677538e-01 time = 0.02 sec -[ Info: VUMPS 13: obj = -2.684102122959e-01 err = 3.5839329032e-01 time = 0.02 sec -[ Info: VUMPS 14: obj = -3.140375987930e-01 err = 3.5220148478e-01 time = 0.07 sec -[ Info: VUMPS 15: obj = -3.494719807888e-01 err = 3.0216614383e-01 time = 0.03 sec -[ Info: VUMPS 16: obj = -2.600823546903e-01 err = 3.4812671114e-01 time = 0.02 sec -[ Info: VUMPS 17: obj = -2.378681434544e-01 err = 3.7837254224e-01 time = 0.02 sec -[ Info: VUMPS 18: obj = -1.165503147031e-01 err = 3.9507067038e-01 time = 0.02 sec -[ Info: VUMPS 19: obj = -3.371493327212e-01 err = 3.2991102430e-01 time = 0.02 sec -[ Info: VUMPS 20: obj = -2.694806272840e-01 err = 3.6205665137e-01 time = 0.05 sec -[ Info: VUMPS 21: obj = -2.381146966165e-01 err = 3.6628398837e-01 time = 0.03 sec -[ Info: VUMPS 22: obj = -3.533056802218e-01 err = 3.0627687902e-01 time = 0.02 sec -[ Info: VUMPS 23: obj = -3.323562393144e-01 err = 3.2310539990e-01 time = 0.03 sec -[ Info: VUMPS 24: obj = -2.930878664518e-01 err = 3.5261734042e-01 time = 0.03 sec -[ Info: VUMPS 25: obj = -3.195260603801e-01 err = 3.4213237114e-01 time = 0.02 sec -[ Info: VUMPS 26: obj = -3.009445857195e-01 err = 3.6267895477e-01 time = 0.05 sec -[ Info: VUMPS 27: obj = -1.945482035700e-01 err = 3.8036687114e-01 time = 0.02 sec -[ Info: VUMPS 28: obj = -3.129608668588e-01 err = 3.3946739811e-01 time = 0.02 sec -[ Info: VUMPS 29: obj = -2.685329650730e-01 err = 3.6719306482e-01 time = 0.02 sec -[ Info: VUMPS 30: obj = -1.255622850827e-01 err = 3.9830069704e-01 time = 0.02 sec -[ Info: VUMPS 31: obj = -2.259105658052e-01 err = 3.5495588739e-01 time = 0.02 sec -[ Info: VUMPS 32: obj = -2.961189169588e-01 err = 3.4982483151e-01 time = 0.02 sec -[ Info: VUMPS 33: obj = -3.982754765049e-01 err = 2.4610327930e-01 time = 0.05 sec -[ Info: VUMPS 34: obj = -3.868968141119e-01 err = 2.6461805850e-01 time = 0.08 sec -[ Info: VUMPS 35: obj = -1.455777740261e-01 err = 4.0452794718e-01 time = 0.02 sec -[ Info: VUMPS 36: obj = -2.384731273529e-01 err = 3.7152623297e-01 time = 0.02 sec -[ Info: VUMPS 37: obj = -2.655148561526e-01 err = 3.5351438817e-01 time = 0.02 sec -[ Info: VUMPS 38: obj = -1.824684314123e-01 err = 3.8362950466e-01 time = 0.02 sec -[ Info: VUMPS 39: obj = -4.340710641979e-02 err = 3.6536524706e-01 time = 0.06 sec -[ Info: VUMPS 40: obj = -7.432642869747e-02 err = 4.1533450888e-01 time = 0.02 sec -[ Info: VUMPS 41: obj = -1.277746081261e-01 err = 4.0325440215e-01 time = 0.02 sec -[ Info: VUMPS 42: obj = -1.637433441229e-01 err = 3.9644561915e-01 time = 0.02 sec -[ Info: VUMPS 43: obj = -2.494122089360e-02 err = 3.9811228669e-01 time = 0.01 sec -[ Info: VUMPS 44: obj = -5.437621624074e-02 err = 4.1086160744e-01 time = 0.02 sec -[ Info: VUMPS 45: obj = -9.499156814866e-02 err = 3.9556960135e-01 time = 0.02 sec -[ Info: VUMPS 46: obj = -2.664160174156e-01 err = 3.4236577100e-01 time = 0.06 sec -[ Info: VUMPS 47: obj = -1.712858760526e-01 err = 3.9929081714e-01 time = 0.02 sec -[ Info: VUMPS 48: obj = -2.278081372696e-01 err = 3.6618881136e-01 time = 0.02 sec -[ Info: VUMPS 49: obj = -3.229675684963e-01 err = 3.1862991140e-01 time = 0.02 sec -[ Info: VUMPS 50: obj = -2.110632456326e-01 err = 3.8859303819e-01 time = 0.02 sec -[ Info: VUMPS 51: obj = -1.101060326910e-01 err = 3.3901651620e-01 time = 0.03 sec -[ Info: VUMPS 52: obj = +7.333338049546e-02 err = 3.6898787075e-01 time = 0.04 sec -[ Info: VUMPS 53: obj = -2.159728827455e-01 err = 3.4788128440e-01 time = 0.02 sec -[ Info: VUMPS 54: obj = -3.781744951044e-01 err = 2.8734162232e-01 time = 0.02 sec -[ Info: VUMPS 55: obj = -2.498878569904e-01 err = 3.7437122967e-01 time = 0.02 sec -[ Info: VUMPS 56: obj = -9.182753559612e-02 err = 4.1119786611e-01 time = 0.02 sec -[ Info: VUMPS 57: obj = -9.896252196567e-02 err = 3.9430375805e-01 time = 0.01 sec -[ Info: VUMPS 58: obj = +9.467977464259e-02 err = 3.8040243513e-01 time = 0.02 sec -[ Info: VUMPS 59: obj = -4.749770024744e-02 err = 3.6638080888e-01 time = 0.05 sec -[ Info: VUMPS 60: obj = -2.530337073792e-01 err = 3.5986314085e-01 time = 0.01 sec -[ Info: VUMPS 61: obj = -1.569761362938e-01 err = 3.8385844092e-01 time = 0.02 sec -[ Info: VUMPS 62: obj = -2.259009259309e-01 err = 3.7057427131e-01 time = 0.02 sec -[ Info: VUMPS 63: obj = +6.333173021014e-02 err = 4.1117754784e-01 time = 0.02 sec -[ Info: VUMPS 64: obj = -1.342383442709e-01 err = 3.6454520924e-01 time = 0.01 sec -[ Info: VUMPS 65: obj = -2.731276613159e-01 err = 3.5374290954e-01 time = 0.02 sec -[ Info: VUMPS 66: obj = -1.467205296305e-02 err = 3.7407661717e-01 time = 0.05 sec -[ Info: VUMPS 67: obj = -7.117743661345e-02 err = 3.5898879533e-01 time = 0.02 sec -[ Info: VUMPS 68: obj = -2.109524282916e-01 err = 3.7190622381e-01 time = 0.02 sec -[ Info: VUMPS 69: obj = -2.429823175649e-01 err = 3.5944789711e-01 time = 0.02 sec -[ Info: VUMPS 70: obj = -2.401276795251e-01 err = 3.5294125102e-01 time = 0.02 sec -[ Info: VUMPS 71: obj = -3.770033736284e-01 err = 2.7476469956e-01 time = 0.02 sec -[ Info: VUMPS 72: obj = +6.405996956936e-02 err = 3.5202436851e-01 time = 0.05 sec -[ Info: VUMPS 73: obj = -2.333733284165e-01 err = 3.6882340641e-01 time = 0.02 sec -[ Info: VUMPS 74: obj = -8.840535012858e-02 err = 3.8308266948e-01 time = 0.02 sec -[ Info: VUMPS 75: obj = -3.397593157665e-01 err = 3.2422050821e-01 time = 0.02 sec -[ Info: VUMPS 76: obj = -1.577057931895e-01 err = 3.7953340218e-01 time = 0.02 sec -[ Info: VUMPS 77: obj = -1.283119074775e-01 err = 3.8701219784e-01 time = 0.02 sec -[ Info: VUMPS 78: obj = -2.838119695877e-01 err = 3.5153675863e-01 time = 0.02 sec -[ Info: VUMPS 79: obj = -3.773014424183e-01 err = 2.8096984820e-01 time = 0.05 sec -[ Info: VUMPS 80: obj = -4.248680541774e-01 err = 1.6995008124e-01 time = 0.02 sec -[ Info: VUMPS 81: obj = -3.898082089496e-01 err = 2.6625454039e-01 time = 0.03 sec -[ Info: VUMPS 82: obj = -9.519451168265e-03 err = 4.2221971320e-01 time = 0.02 sec -[ Info: VUMPS 83: obj = -5.231292682682e-02 err = 3.6959419670e-01 time = 0.02 sec -[ Info: VUMPS 84: obj = -1.866826486661e-01 err = 3.8258698851e-01 time = 0.01 sec -[ Info: VUMPS 85: obj = -1.692405363751e-01 err = 3.8770899838e-01 time = 0.05 sec -[ Info: VUMPS 86: obj = -2.083324644893e-01 err = 3.5431521899e-01 time = 0.02 sec -[ Info: VUMPS 87: obj = -2.939979120876e-01 err = 3.5532871233e-01 time = 0.02 sec -[ Info: VUMPS 88: obj = -2.828666406513e-01 err = 3.4873331725e-01 time = 0.02 sec -[ Info: VUMPS 89: obj = -1.692982578604e-01 err = 3.8168691441e-01 time = 0.02 sec -[ Info: VUMPS 90: obj = -2.363440425785e-01 err = 3.4723691699e-01 time = 0.02 sec -[ Info: VUMPS 91: obj = -1.414396078859e-01 err = 4.0303952514e-01 time = 0.05 sec -[ Info: VUMPS 92: obj = -9.111717266883e-02 err = 3.9613002109e-01 time = 0.02 sec -[ Info: VUMPS 93: obj = +3.874548352332e-02 err = 3.6730004996e-01 time = 0.02 sec -[ Info: VUMPS 94: obj = -7.194785728285e-02 err = 3.9110410465e-01 time = 0.02 sec -[ Info: VUMPS 95: obj = -1.043822044301e-01 err = 3.9084477906e-01 time = 0.02 sec -[ Info: VUMPS 96: obj = -1.029644744149e-01 err = 4.1397652279e-01 time = 0.02 sec -[ Info: VUMPS 97: obj = -2.861705477721e-01 err = 3.5209989794e-01 time = 0.01 sec -[ Info: VUMPS 98: obj = -3.470187676038e-01 err = 3.0920621601e-01 time = 0.05 sec -[ Info: VUMPS 99: obj = -5.822199548527e-02 err = 3.5959218485e-01 time = 0.02 sec -[ Info: VUMPS 100: obj = -1.847189630830e-01 err = 3.6268585762e-01 time = 0.02 sec -[ Info: VUMPS 101: obj = -1.998247793394e-01 err = 3.7026241762e-01 time = 0.02 sec -[ Info: VUMPS 102: obj = -2.845861059390e-01 err = 3.5015115387e-01 time = 0.02 sec -[ Info: VUMPS 103: obj = +1.001145171552e-01 err = 3.7897389825e-01 time = 0.02 sec -[ Info: VUMPS 104: obj = +5.388940615771e-03 err = 3.8837861849e-01 time = 0.02 sec -[ Info: VUMPS 105: obj = +6.174310736557e-02 err = 3.6429568298e-01 time = 0.05 sec -[ Info: VUMPS 106: obj = -1.953458747936e-02 err = 4.1105781305e-01 time = 0.02 sec -[ Info: VUMPS 107: obj = -1.227404516478e-01 err = 3.9527799477e-01 time = 0.02 sec -[ Info: VUMPS 108: obj = -2.249890803810e-01 err = 3.7797056287e-01 time = 0.02 sec -[ Info: VUMPS 109: obj = -6.728916323045e-02 err = 3.6671401448e-01 time = 0.03 sec -[ Info: VUMPS 110: obj = -1.503938692777e-01 err = 3.7432233900e-01 time = 0.02 sec -[ Info: VUMPS 111: obj = -5.369841176091e-02 err = 4.2078297478e-01 time = 0.02 sec -[ Info: VUMPS 112: obj = +3.503395811604e-02 err = 4.0241955891e-01 time = 0.05 sec -[ Info: VUMPS 113: obj = -1.641757316788e-01 err = 3.8361173302e-01 time = 0.01 sec -[ Info: VUMPS 114: obj = -2.285888795528e-01 err = 3.6399759156e-01 time = 0.02 sec -[ Info: VUMPS 115: obj = -3.184239125421e-01 err = 3.5620987147e-01 time = 0.03 sec -[ Info: VUMPS 116: obj = -4.004343201842e-01 err = 2.3427592725e-01 time = 0.02 sec -[ Info: VUMPS 117: obj = +1.250815154369e-01 err = 3.8171291086e-01 time = 0.02 sec -[ Info: VUMPS 118: obj = -8.382309493919e-02 err = 3.9885688763e-01 time = 0.05 sec -[ Info: VUMPS 119: obj = -1.480032425391e-01 err = 3.7961317585e-01 time = 0.01 sec -[ Info: VUMPS 120: obj = -2.603417820356e-01 err = 3.5989698304e-01 time = 0.02 sec -[ Info: VUMPS 121: obj = -3.231635462838e-01 err = 3.3803029927e-01 time = 0.02 sec -[ Info: VUMPS 122: obj = -3.960892923915e-01 err = 2.5894024392e-01 time = 0.02 sec -[ Info: VUMPS 123: obj = +1.383960399671e-01 err = 3.2560793126e-01 time = 0.02 sec -[ Info: VUMPS 124: obj = +2.266107848362e-03 err = 4.1080709820e-01 time = 0.02 sec -[ Info: VUMPS 125: obj = -2.283058118711e-01 err = 3.7356440605e-01 time = 0.05 sec -[ Info: VUMPS 126: obj = -3.418107299565e-01 err = 3.1899399020e-01 time = 0.02 sec -[ Info: VUMPS 127: obj = -2.802812290403e-01 err = 3.5828911261e-01 time = 0.02 sec -[ Info: VUMPS 128: obj = -1.362964889519e-01 err = 3.9593008281e-01 time = 0.02 sec -[ Info: VUMPS 129: obj = -6.715310983028e-02 err = 3.8272455592e-01 time = 0.02 sec -[ Info: VUMPS 130: obj = -2.592620383388e-01 err = 3.5195567454e-01 time = 0.02 sec -[ Info: VUMPS 131: obj = -4.252097482878e-01 err = 1.6873141588e-01 time = 0.05 sec -[ Info: VUMPS 132: obj = +3.863457166648e-02 err = 3.2307300405e-01 time = 0.02 sec -[ Info: VUMPS 133: obj = +2.300685431813e-03 err = 3.8827309741e-01 time = 0.02 sec -[ Info: VUMPS 134: obj = -1.800971988749e-01 err = 3.9325052328e-01 time = 0.02 sec -[ Info: VUMPS 135: obj = +2.655829501197e-03 err = 3.8180278781e-01 time = 0.02 sec -[ Info: VUMPS 136: obj = -1.960012423358e-01 err = 3.7736036559e-01 time = 0.02 sec -[ Info: VUMPS 137: obj = -3.945976063716e-01 err = 2.5776349595e-01 time = 0.05 sec -[ Info: VUMPS 138: obj = +9.835127345556e-02 err = 3.9127467384e-01 time = 0.02 sec -[ Info: VUMPS 139: obj = -1.960607007481e-01 err = 3.9818909217e-01 time = 0.02 sec -[ Info: VUMPS 140: obj = -2.076438228078e-01 err = 3.6858724855e-01 time = 0.01 sec -[ Info: VUMPS 141: obj = -2.752417182485e-01 err = 3.5409227604e-01 time = 0.02 sec -[ Info: VUMPS 142: obj = -2.323570171993e-01 err = 3.8479850876e-01 time = 0.03 sec -[ Info: VUMPS 143: obj = -2.478431862183e-01 err = 3.6601199783e-01 time = 0.02 sec -[ Info: VUMPS 144: obj = -1.303161187004e-01 err = 3.9306929113e-01 time = 0.05 sec -[ Info: VUMPS 145: obj = -1.753964431937e-01 err = 3.8947518565e-01 time = 0.02 sec -[ Info: VUMPS 146: obj = +4.823037317739e-02 err = 3.7010938139e-01 time = 0.02 sec -[ Info: VUMPS 147: obj = -1.720633754678e-01 err = 3.8439582822e-01 time = 0.02 sec -[ Info: VUMPS 148: obj = -1.497805800054e-01 err = 3.8708211808e-01 time = 0.02 sec -[ Info: VUMPS 149: obj = -1.574371937148e-01 err = 3.9621587227e-01 time = 0.02 sec -[ Info: VUMPS 150: obj = -4.188596075056e-02 err = 4.0489025033e-01 time = 0.05 sec -[ Info: VUMPS 151: obj = -2.931751248451e-01 err = 3.3701305583e-01 time = 0.02 sec -[ Info: VUMPS 152: obj = -3.542385510115e-01 err = 3.0922083054e-01 time = 0.02 sec -[ Info: VUMPS 153: obj = -3.153636373307e-01 err = 3.2780719969e-01 time = 0.02 sec -[ Info: VUMPS 154: obj = -2.099485771080e-01 err = 3.7583374896e-01 time = 0.03 sec -[ Info: VUMPS 155: obj = -2.325036131152e-01 err = 3.6785157612e-01 time = 0.02 sec -[ Info: VUMPS 156: obj = -3.306527885138e-01 err = 3.1451420802e-01 time = 0.05 sec -[ Info: VUMPS 157: obj = -3.683042109253e-01 err = 3.0098744465e-01 time = 0.03 sec -[ Info: VUMPS 158: obj = -2.899970230712e-01 err = 3.4912426816e-01 time = 0.03 sec -[ Info: VUMPS 159: obj = -1.548982258892e-01 err = 3.9634701519e-01 time = 0.02 sec -[ Info: VUMPS 160: obj = -1.807853460354e-01 err = 4.0347363337e-01 time = 0.02 sec -[ Info: VUMPS 161: obj = -9.508228951537e-02 err = 3.7574701420e-01 time = 0.02 sec -[ Info: VUMPS 162: obj = -1.993563335295e-01 err = 3.5850713945e-01 time = 0.05 sec -[ Info: VUMPS 163: obj = +9.436575352874e-02 err = 3.6482840550e-01 time = 0.02 sec -[ Info: VUMPS 164: obj = +7.903937847467e-02 err = 3.8855241138e-01 time = 0.02 sec -[ Info: VUMPS 165: obj = -5.495752329074e-02 err = 3.8048454765e-01 time = 0.01 sec -[ Info: VUMPS 166: obj = -2.977368905239e-01 err = 3.5841635043e-01 time = 0.02 sec -[ Info: VUMPS 167: obj = -2.285161464463e-01 err = 3.6155361498e-01 time = 0.02 sec -[ Info: VUMPS 168: obj = -2.969068468462e-01 err = 3.3763604978e-01 time = 0.02 sec -[ Info: VUMPS 169: obj = -3.372311890964e-01 err = 3.1700106891e-01 time = 0.05 sec -[ Info: VUMPS 170: obj = -3.045611491085e-01 err = 3.0628902149e-01 time = 0.03 sec -[ Info: VUMPS 171: obj = -3.017393279199e-01 err = 3.6408415132e-01 time = 0.02 sec -[ Info: VUMPS 172: obj = -1.249057059582e-01 err = 3.8003500781e-01 time = 0.02 sec -[ Info: VUMPS 173: obj = -1.804070758425e-01 err = 3.8048914783e-01 time = 0.02 sec -[ Info: VUMPS 174: obj = -1.587740528232e-01 err = 3.8521882028e-01 time = 0.02 sec -[ Info: VUMPS 175: obj = -2.851791605376e-01 err = 3.3829945088e-01 time = 0.05 sec -[ Info: VUMPS 176: obj = -6.876942688524e-02 err = 4.0447417286e-01 time = 0.02 sec -[ Info: VUMPS 177: obj = +4.088424378843e-02 err = 3.9727430522e-01 time = 0.02 sec -[ Info: VUMPS 178: obj = +4.349001316215e-02 err = 3.6572321243e-01 time = 0.02 sec -[ Info: VUMPS 179: obj = -2.273704295370e-01 err = 3.8731200334e-01 time = 0.02 sec -[ Info: VUMPS 180: obj = -2.616628489407e-01 err = 3.4696217488e-01 time = 0.02 sec -[ Info: VUMPS 181: obj = -2.694623492030e-01 err = 3.4874627834e-01 time = 0.02 sec -[ Info: VUMPS 182: obj = -2.783882521737e-01 err = 3.5076502054e-01 time = 0.05 sec -[ Info: VUMPS 183: obj = -2.136144754439e-01 err = 3.8565666056e-01 time = 0.02 sec -[ Info: VUMPS 184: obj = -2.172901485715e-01 err = 3.8072126413e-01 time = 0.02 sec -[ Info: VUMPS 185: obj = +8.912737400891e-02 err = 3.5445943853e-01 time = 0.03 sec -[ Info: VUMPS 186: obj = -4.405574529414e-02 err = 4.0903072415e-01 time = 0.02 sec -[ Info: VUMPS 187: obj = -1.606066852805e-01 err = 3.8640765878e-01 time = 0.02 sec -[ Info: VUMPS 188: obj = +2.554966161756e-02 err = 3.7067877383e-01 time = 0.05 sec -[ Info: VUMPS 189: obj = -1.520562341146e-01 err = 3.7246997774e-01 time = 0.01 sec -[ Info: VUMPS 190: obj = -5.455120357931e-02 err = 3.7880083957e-01 time = 0.02 sec -[ Info: VUMPS 191: obj = -1.821662375375e-01 err = 3.9326693657e-01 time = 0.02 sec -[ Info: VUMPS 192: obj = -1.577721407955e-01 err = 4.1429770418e-01 time = 0.03 sec -[ Info: VUMPS 193: obj = +2.586895357233e-02 err = 3.5544800306e-01 time = 0.02 sec -[ Info: VUMPS 194: obj = -2.583240143294e-01 err = 3.7987694829e-01 time = 0.02 sec -[ Info: VUMPS 195: obj = -2.224656299792e-01 err = 3.6394963764e-01 time = 0.05 sec -[ Info: VUMPS 196: obj = -1.254311298138e-01 err = 3.7548601279e-01 time = 0.02 sec -[ Info: VUMPS 197: obj = -1.723763414272e-01 err = 3.5494740522e-01 time = 0.02 sec -[ Info: VUMPS 198: obj = -2.762110763835e-01 err = 3.4218268756e-01 time = 0.02 sec -[ Info: VUMPS 199: obj = -3.239089102881e-01 err = 3.2528270621e-01 time = 0.02 sec -┌ Warning: VUMPS cancel 200: obj = -2.870558938604e-01 err = 3.3497651674e-01 time = 5.07 sec -└ @ MPSKit ~/Projects/MPSKit.jl/src/algorithms/groundstate/vumps.jl:73 +[ Info: VUMPS init: obj = +2.499952940580e-01 err = 6.3441e-03 +[ Info: VUMPS 1: obj = -1.289444927303e-01 err = 3.6505534112e-01 time = 0.02 sec +[ Info: VUMPS 2: obj = -3.046604131170e-01 err = 3.5438451443e-01 time = 0.03 sec +[ Info: VUMPS 3: obj = -1.151380226837e-01 err = 3.7396849153e-01 time = 0.03 sec +[ Info: VUMPS 4: obj = -7.859906223027e-02 err = 3.8101931668e-01 time = 0.02 sec +[ Info: VUMPS 5: obj = -8.570311273681e-02 err = 4.5415600150e-01 time = 0.04 sec +[ Info: VUMPS 6: obj = -1.607866955582e-01 err = 3.7328533114e-01 time = 0.02 sec +[ Info: VUMPS 7: obj = -2.590055398932e-01 err = 3.6315039435e-01 time = 0.02 sec +[ Info: VUMPS 8: obj = -3.985401345056e-01 err = 2.4797713082e-01 time = 0.03 sec +[ Info: VUMPS 9: obj = -2.980269171709e-01 err = 3.5009326244e-01 time = 0.04 sec +[ Info: VUMPS 10: obj = -7.766811941836e-02 err = 3.7146004460e-01 time = 0.04 sec +[ Info: VUMPS 11: obj = -1.513433486022e-01 err = 3.9298119907e-01 time = 0.04 sec +[ Info: VUMPS 12: obj = -1.636153457041e-01 err = 3.8567245083e-01 time = 0.02 sec +[ Info: VUMPS 13: obj = -8.850847675954e-02 err = 3.8567080613e-01 time = 0.03 sec +[ Info: VUMPS 14: obj = -1.933682984363e-01 err = 3.8578605637e-01 time = 0.08 sec +[ Info: VUMPS 15: obj = -4.668791737320e-02 err = 4.0507831578e-01 time = 0.03 sec +[ Info: VUMPS 16: obj = -6.685207964181e-02 err = 4.4181563160e-01 time = 0.03 sec +[ Info: VUMPS 17: obj = -2.672203890897e-01 err = 3.3906178990e-01 time = 0.03 sec +[ Info: VUMPS 18: obj = -3.428995365711e-01 err = 3.2686747553e-01 time = 0.03 sec +[ Info: VUMPS 19: obj = -2.863734743969e-01 err = 3.5064087667e-01 time = 0.04 sec +[ Info: VUMPS 20: obj = -1.067319756533e-01 err = 3.8187305789e-01 time = 0.03 sec +[ Info: VUMPS 21: obj = -1.318548228297e-01 err = 3.6108954481e-01 time = 0.02 sec +[ Info: VUMPS 22: obj = -4.044279153596e-02 err = 3.7520252599e-01 time = 0.02 sec +[ Info: VUMPS 23: obj = -3.945104968675e-02 err = 3.9224571991e-01 time = 0.03 sec +[ Info: VUMPS 24: obj = -2.671412452672e-01 err = 3.6585228539e-01 time = 0.03 sec +[ Info: VUMPS 25: obj = -3.310347250642e-01 err = 3.1899906877e-01 time = 0.03 sec +[ Info: VUMPS 26: obj = -2.168383952062e-01 err = 3.7166034087e-01 time = 0.03 sec +[ Info: VUMPS 27: obj = -3.507015235665e-01 err = 3.0330946464e-01 time = 0.04 sec +[ Info: VUMPS 28: obj = -2.810889490868e-01 err = 3.5384589143e-01 time = 0.04 sec +[ Info: VUMPS 29: obj = -1.491679715343e-01 err = 3.7712897712e-01 time = 0.03 sec +[ Info: VUMPS 30: obj = -2.103310942986e-01 err = 3.6804474254e-01 time = 0.02 sec +[ Info: VUMPS 31: obj = -2.570407107876e-01 err = 3.5647467556e-01 time = 0.03 sec +[ Info: VUMPS 32: obj = -2.622474733096e-01 err = 3.6359049960e-01 time = 0.03 sec +[ Info: VUMPS 33: obj = -2.747408111088e-01 err = 3.5217438909e-01 time = 0.02 sec +[ Info: VUMPS 34: obj = -2.795011381225e-01 err = 3.5577849785e-01 time = 0.03 sec +[ Info: VUMPS 35: obj = +6.109686520960e-02 err = 3.7478423163e-01 time = 0.03 sec +[ Info: VUMPS 36: obj = -3.235177003488e-01 err = 3.2501656960e-01 time = 0.03 sec +[ Info: VUMPS 37: obj = -3.607718277418e-01 err = 3.0399210370e-01 time = 0.03 sec +[ Info: VUMPS 38: obj = -1.522405023809e-01 err = 3.8460611720e-01 time = 0.03 sec +[ Info: VUMPS 39: obj = -1.401486230679e-01 err = 4.2359144803e-01 time = 0.04 sec +[ Info: VUMPS 40: obj = -3.016337945443e-01 err = 3.3760785136e-01 time = 0.02 sec +[ Info: VUMPS 41: obj = -3.592630330471e-01 err = 2.8277667197e-01 time = 0.04 sec +[ Info: VUMPS 42: obj = -9.648261609327e-02 err = 3.9800363224e-01 time = 0.03 sec +[ Info: VUMPS 43: obj = -2.068120169738e-01 err = 3.4795790120e-01 time = 0.07 sec +[ Info: VUMPS 44: obj = -3.320461487347e-01 err = 3.3031447114e-01 time = 0.03 sec +[ Info: VUMPS 45: obj = -2.748266113721e-01 err = 3.3719047760e-01 time = 0.03 sec +[ Info: VUMPS 46: obj = +5.825197478980e-02 err = 3.7299269867e-01 time = 0.03 sec +[ Info: VUMPS 47: obj = -2.209737585465e-01 err = 3.6224489048e-01 time = 0.03 sec +[ Info: VUMPS 48: obj = -3.495730972764e-01 err = 3.1370483925e-01 time = 0.03 sec +[ Info: VUMPS 49: obj = -1.184295942341e-02 err = 3.8132184336e-01 time = 0.03 sec +[ Info: VUMPS 50: obj = -2.253818587474e-01 err = 3.6351801471e-01 time = 0.02 sec +[ Info: VUMPS 51: obj = -2.611723114106e-01 err = 3.6263568517e-01 time = 0.04 sec +[ Info: VUMPS 52: obj = -2.391735319717e-01 err = 3.4034883778e-01 time = 0.03 sec +[ Info: VUMPS 53: obj = -2.076637304355e-01 err = 3.7785371429e-01 time = 0.03 sec +[ Info: VUMPS 54: obj = +1.233357556406e-01 err = 3.6817687888e-01 time = 0.03 sec +[ Info: VUMPS 55: obj = -1.064547626837e-02 err = 3.6875802199e-01 time = 0.03 sec +[ Info: VUMPS 56: obj = -1.051495232666e-01 err = 4.0131031658e-01 time = 0.03 sec +[ Info: VUMPS 57: obj = -1.412941986925e-01 err = 3.8676020133e-01 time = 0.02 sec +[ Info: VUMPS 58: obj = -5.105594856383e-02 err = 4.1918593555e-01 time = 0.03 sec +[ Info: VUMPS 59: obj = -9.429180669034e-02 err = 3.5663676393e-01 time = 0.02 sec +[ Info: VUMPS 60: obj = -3.211645595234e-01 err = 3.3556696008e-01 time = 0.03 sec +[ Info: VUMPS 61: obj = -4.997897659120e-02 err = 3.6911454941e-01 time = 0.03 sec +[ Info: VUMPS 62: obj = +5.061815735283e-02 err = 3.9240870402e-01 time = 0.03 sec +[ Info: VUMPS 63: obj = -2.858216622151e-02 err = 3.9031720030e-01 time = 0.03 sec +[ Info: VUMPS 64: obj = -3.066897558850e-01 err = 3.4219120029e-01 time = 0.02 sec +[ Info: VUMPS 65: obj = -1.035726948011e-01 err = 3.9992774981e-01 time = 0.03 sec +[ Info: VUMPS 66: obj = +9.359087106400e-02 err = 3.2551512023e-01 time = 0.03 sec +[ Info: VUMPS 67: obj = -1.506323639493e-01 err = 3.8851758838e-01 time = 0.03 sec +[ Info: VUMPS 68: obj = -2.694024194911e-01 err = 3.4908885270e-01 time = 0.02 sec +[ Info: VUMPS 69: obj = -1.666459363907e-01 err = 3.7191916042e-01 time = 0.03 sec +[ Info: VUMPS 70: obj = -3.295401946651e-01 err = 3.2485997665e-01 time = 0.02 sec +[ Info: VUMPS 71: obj = -1.716492900593e-01 err = 3.4774685002e-01 time = 0.05 sec +[ Info: VUMPS 72: obj = -2.062243894572e-01 err = 3.4367048189e-01 time = 0.03 sec +[ Info: VUMPS 73: obj = -2.270916009662e-01 err = 3.4883289649e-01 time = 0.04 sec +[ Info: VUMPS 74: obj = +6.972996650156e-03 err = 3.6664495072e-01 time = 0.03 sec +[ Info: VUMPS 75: obj = -2.140897386871e-01 err = 3.6896873908e-01 time = 0.02 sec +[ Info: VUMPS 76: obj = -2.776215003655e-01 err = 3.5162263908e-01 time = 0.03 sec +[ Info: VUMPS 77: obj = -2.548517060346e-01 err = 3.5985886200e-01 time = 0.04 sec +[ Info: VUMPS 78: obj = -3.627093129660e-01 err = 2.8522543561e-01 time = 0.04 sec +[ Info: VUMPS 79: obj = -4.246979337066e-01 err = 1.8433888114e-01 time = 0.04 sec +[ Info: VUMPS 80: obj = -2.768088632211e-01 err = 3.6889157455e-01 time = 0.04 sec +[ Info: VUMPS 81: obj = -3.333529535371e-01 err = 3.2461308190e-01 time = 0.04 sec +[ Info: VUMPS 82: obj = +9.739742280028e-02 err = 3.6140350665e-01 time = 0.04 sec +[ Info: VUMPS 83: obj = -8.001219240663e-02 err = 3.3065484133e-01 time = 0.03 sec +[ Info: VUMPS 84: obj = -1.494867091384e-01 err = 3.6515333400e-01 time = 0.02 sec +[ Info: VUMPS 85: obj = -3.357585607378e-01 err = 3.2062204533e-01 time = 0.03 sec +[ Info: VUMPS 86: obj = -1.214350882454e-01 err = 3.8756516414e-01 time = 0.03 sec +[ Info: VUMPS 87: obj = -1.720890136641e-01 err = 3.7209498473e-01 time = 0.03 sec +[ Info: VUMPS 88: obj = -2.513176623582e-01 err = 3.6428215490e-01 time = 0.03 sec +[ Info: VUMPS 89: obj = -1.167825455975e-01 err = 3.7816700280e-01 time = 0.03 sec +[ Info: VUMPS 90: obj = -1.487085246354e-01 err = 3.8452112070e-01 time = 0.03 sec +[ Info: VUMPS 91: obj = -1.492511193594e-01 err = 3.6888008833e-01 time = 0.04 sec +[ Info: VUMPS 92: obj = -2.942886788237e-01 err = 3.3926362090e-01 time = 0.02 sec +[ Info: VUMPS 93: obj = -2.513546381842e-01 err = 3.5669021055e-01 time = 0.04 sec +[ Info: VUMPS 94: obj = -2.592696224289e-01 err = 3.8742412819e-01 time = 0.04 sec +[ Info: VUMPS 95: obj = -1.670969129733e-01 err = 4.0396565424e-01 time = 0.04 sec +[ Info: VUMPS 96: obj = -8.982884784032e-02 err = 4.0673356787e-01 time = 0.04 sec +[ Info: VUMPS 97: obj = -1.162574886965e-01 err = 4.0014599231e-01 time = 0.03 sec +[ Info: VUMPS 98: obj = -9.485766256473e-02 err = 3.6567908225e-01 time = 0.02 sec +[ Info: VUMPS 99: obj = -1.626897465348e-01 err = 3.9836027057e-01 time = 0.06 sec +[ Info: VUMPS 100: obj = -2.820414385240e-01 err = 3.3860149095e-01 time = 0.02 sec +[ Info: VUMPS 101: obj = +3.653164331064e-03 err = 3.8214630304e-01 time = 0.04 sec +[ Info: VUMPS 102: obj = -1.488426870678e-01 err = 3.5921878796e-01 time = 0.02 sec +[ Info: VUMPS 103: obj = -3.107550441980e-01 err = 3.3205507908e-01 time = 0.03 sec +[ Info: VUMPS 104: obj = -2.991475468045e-01 err = 3.4083808929e-01 time = 0.04 sec +[ Info: VUMPS 105: obj = -1.521003309052e-01 err = 4.0772463716e-01 time = 0.05 sec +[ Info: VUMPS 106: obj = -8.123652136498e-02 err = 3.7546849297e-01 time = 0.03 sec +[ Info: VUMPS 107: obj = -2.596284651709e-01 err = 3.4976630049e-01 time = 0.03 sec +[ Info: VUMPS 108: obj = -2.935877769698e-01 err = 3.3783749995e-01 time = 0.04 sec +[ Info: VUMPS 109: obj = -2.402370295781e-01 err = 3.6810295199e-01 time = 0.04 sec +[ Info: VUMPS 110: obj = -3.150592845722e-01 err = 3.3233905128e-01 time = 0.04 sec +[ Info: VUMPS 111: obj = -3.266942906048e-01 err = 3.3265129842e-01 time = 0.05 sec +[ Info: VUMPS 112: obj = -3.089623364833e-02 err = 3.8736213578e-01 time = 0.03 sec +[ Info: VUMPS 113: obj = +1.187921857755e-01 err = 3.5846597520e-01 time = 0.03 sec +[ Info: VUMPS 114: obj = -6.812831791961e-02 err = 3.6601110780e-01 time = 0.03 sec +[ Info: VUMPS 115: obj = -2.169545739645e-02 err = 3.7328578080e-01 time = 0.02 sec +[ Info: VUMPS 116: obj = -6.651891606433e-02 err = 3.7799893055e-01 time = 0.03 sec +[ Info: VUMPS 117: obj = -7.188213532364e-02 err = 3.9041626270e-01 time = 0.02 sec +[ Info: VUMPS 118: obj = -1.858708677528e-01 err = 3.8364886262e-01 time = 0.02 sec +[ Info: VUMPS 119: obj = -1.426092420792e-01 err = 4.0022614912e-01 time = 0.04 sec +[ Info: VUMPS 120: obj = -5.966800560939e-02 err = 3.7475162153e-01 time = 0.03 sec +[ Info: VUMPS 121: obj = -1.927011968813e-01 err = 3.7213373940e-01 time = 0.02 sec +[ Info: VUMPS 122: obj = -1.332482420009e-01 err = 3.9823332335e-01 time = 0.03 sec +[ Info: VUMPS 123: obj = -5.583604028721e-02 err = 3.8811845128e-01 time = 0.02 sec +[ Info: VUMPS 124: obj = -1.173769997886e-01 err = 3.9758982598e-01 time = 0.03 sec +[ Info: VUMPS 125: obj = -1.623302909431e-01 err = 3.9006196969e-01 time = 0.03 sec +[ Info: VUMPS 126: obj = -1.945385790084e-01 err = 3.7801724719e-01 time = 0.06 sec +[ Info: VUMPS 127: obj = -1.198866351263e-01 err = 3.8641158608e-01 time = 0.04 sec +[ Info: VUMPS 128: obj = -4.982327612008e-02 err = 3.8983422860e-01 time = 0.02 sec +[ Info: VUMPS 129: obj = -1.566708742911e-01 err = 3.6570433006e-01 time = 0.03 sec +[ Info: VUMPS 130: obj = -3.746845720087e-01 err = 2.7790316765e-01 time = 0.03 sec +[ Info: VUMPS 131: obj = -4.113990855397e-01 err = 2.1938734166e-01 time = 0.04 sec +[ Info: VUMPS 132: obj = +2.987037930122e-02 err = 3.9820435991e-01 time = 0.03 sec +[ Info: VUMPS 133: obj = -1.361451053253e-01 err = 3.5853656852e-01 time = 0.04 sec +[ Info: VUMPS 134: obj = +1.459826659542e-02 err = 3.6344259860e-01 time = 0.02 sec +[ Info: VUMPS 135: obj = -1.625462338965e-01 err = 3.8750540849e-01 time = 0.03 sec +[ Info: VUMPS 136: obj = -1.945260498995e-01 err = 3.7322465363e-01 time = 0.03 sec +[ Info: VUMPS 137: obj = -2.378874547137e-01 err = 3.5236618024e-01 time = 0.02 sec +[ Info: VUMPS 138: obj = -3.397752713889e-01 err = 3.1681390396e-01 time = 0.04 sec +[ Info: VUMPS 139: obj = +3.872021456876e-02 err = 3.9044746200e-01 time = 0.03 sec +[ Info: VUMPS 140: obj = -2.931409673729e-02 err = 4.0519520800e-01 time = 0.03 sec +[ Info: VUMPS 141: obj = -3.558928281290e-02 err = 3.9408216725e-01 time = 0.03 sec +[ Info: VUMPS 142: obj = -7.419482234602e-02 err = 3.5450836272e-01 time = 0.03 sec +[ Info: VUMPS 143: obj = -1.320472388492e-01 err = 3.8629858112e-01 time = 0.02 sec +[ Info: VUMPS 144: obj = -9.540887842838e-02 err = 3.6612197907e-01 time = 0.03 sec +[ Info: VUMPS 145: obj = -1.425512039944e-01 err = 3.8151598515e-01 time = 0.02 sec +[ Info: VUMPS 146: obj = -1.889146812574e-01 err = 3.7124283218e-01 time = 0.03 sec +[ Info: VUMPS 147: obj = -2.150431811539e-01 err = 3.5676583668e-01 time = 0.02 sec +[ Info: VUMPS 148: obj = -5.190628673026e-02 err = 3.8934562140e-01 time = 0.02 sec +[ Info: VUMPS 149: obj = -1.783612875181e-01 err = 3.6973710837e-01 time = 0.03 sec +[ Info: VUMPS 150: obj = -1.973875305247e-01 err = 3.7126107705e-01 time = 0.02 sec +[ Info: VUMPS 151: obj = -2.899275740844e-01 err = 3.6876702526e-01 time = 0.03 sec +[ Info: VUMPS 152: obj = -1.958917623537e-01 err = 3.6026243367e-01 time = 0.04 sec +[ Info: VUMPS 153: obj = -1.658053486005e-01 err = 3.7110180704e-01 time = 0.03 sec +[ Info: VUMPS 154: obj = -1.826177780294e-01 err = 4.1264125075e-01 time = 0.06 sec +[ Info: VUMPS 155: obj = -1.059120939251e-02 err = 3.9019355474e-01 time = 0.04 sec +[ Info: VUMPS 156: obj = -1.197063315403e-01 err = 3.8266860150e-01 time = 0.03 sec +[ Info: VUMPS 157: obj = -1.518053401692e-01 err = 3.9773151275e-01 time = 0.03 sec +[ Info: VUMPS 158: obj = -3.627713872160e-01 err = 3.0031072551e-01 time = 0.03 sec +[ Info: VUMPS 159: obj = -3.817082298432e-01 err = 2.7493910948e-01 time = 0.04 sec +[ Info: VUMPS 160: obj = -3.859619689587e-01 err = 2.7048867198e-01 time = 0.04 sec +[ Info: VUMPS 161: obj = -4.307827862211e-01 err = 1.2679659771e-01 time = 0.05 sec +[ Info: VUMPS 162: obj = -4.409832181424e-01 err = 5.1946083987e-02 time = 0.05 sec +[ Info: VUMPS 163: obj = -3.974418237321e-01 err = 2.6311596995e-01 time = 0.06 sec +[ Info: VUMPS 164: obj = -1.705693106104e-01 err = 3.7489488650e-01 time = 0.04 sec +[ Info: VUMPS 165: obj = -2.811377309143e-01 err = 3.5373076558e-01 time = 0.02 sec +[ Info: VUMPS 166: obj = -3.580031548549e-01 err = 3.0296648457e-01 time = 0.03 sec +[ Info: VUMPS 167: obj = -3.498256310994e-01 err = 3.0138171612e-01 time = 0.04 sec +[ Info: VUMPS 168: obj = -3.961930709656e-01 err = 2.3106926668e-01 time = 0.05 sec +[ Info: VUMPS 169: obj = -3.012322573839e-01 err = 3.4871517230e-01 time = 0.05 sec +[ Info: VUMPS 170: obj = -3.236118252223e-01 err = 3.2683561230e-01 time = 0.03 sec +[ Info: VUMPS 171: obj = -4.207960751090e-01 err = 1.9439910804e-01 time = 0.04 sec +[ Info: VUMPS 172: obj = +6.486820083537e-02 err = 3.7997567109e-01 time = 0.04 sec +[ Info: VUMPS 173: obj = -2.205989590288e-01 err = 3.5392631042e-01 time = 0.02 sec +[ Info: VUMPS 174: obj = -3.040541044421e-01 err = 3.4795284083e-01 time = 0.02 sec +[ Info: VUMPS 175: obj = -2.700316483862e-01 err = 3.5371268695e-01 time = 0.04 sec +[ Info: VUMPS 176: obj = -1.451879669536e-01 err = 3.7441207912e-01 time = 0.03 sec +[ Info: VUMPS 177: obj = -1.675392520421e-01 err = 3.5821939134e-01 time = 0.03 sec +[ Info: VUMPS 178: obj = -1.067646782592e-01 err = 3.5902159365e-01 time = 0.02 sec +[ Info: VUMPS 179: obj = +5.087477936752e-02 err = 3.8379007723e-01 time = 0.06 sec +[ Info: VUMPS 180: obj = +2.700648055200e-02 err = 3.5560967889e-01 time = 0.03 sec +[ Info: VUMPS 181: obj = -4.040983881654e-02 err = 3.1955366201e-01 time = 0.03 sec +[ Info: VUMPS 182: obj = -2.974889920215e-01 err = 3.4443915929e-01 time = 0.03 sec +[ Info: VUMPS 183: obj = -3.470958294750e-01 err = 3.1478705957e-01 time = 0.03 sec +[ Info: VUMPS 184: obj = -2.257989223447e-01 err = 3.8360838997e-01 time = 0.04 sec +[ Info: VUMPS 185: obj = -1.892095032963e-01 err = 3.9861950100e-01 time = 0.04 sec +[ Info: VUMPS 186: obj = -1.178703030512e-01 err = 3.9409430451e-01 time = 0.03 sec +[ Info: VUMPS 187: obj = -1.803058178275e-01 err = 3.8179418477e-01 time = 0.02 sec +[ Info: VUMPS 188: obj = +2.026839055691e-02 err = 3.9195067790e-01 time = 0.02 sec +[ Info: VUMPS 189: obj = -2.213025303557e-01 err = 3.6496697752e-01 time = 0.03 sec +[ Info: VUMPS 190: obj = -3.308949018046e-01 err = 3.0825188354e-01 time = 0.03 sec +[ Info: VUMPS 191: obj = -3.358030111876e-01 err = 3.2213471244e-01 time = 0.04 sec +[ Info: VUMPS 192: obj = -3.879783512495e-01 err = 2.6554883101e-01 time = 0.04 sec +[ Info: VUMPS 193: obj = -4.235762853319e-01 err = 1.6965146723e-01 time = 0.04 sec +[ Info: VUMPS 194: obj = -2.904118238494e-01 err = 3.5667684109e-01 time = 0.06 sec +[ Info: VUMPS 195: obj = -2.142437785647e-01 err = 3.8141868514e-01 time = 0.04 sec +[ Info: VUMPS 196: obj = -7.074078690463e-02 err = 3.4787569806e-01 time = 0.04 sec +[ Info: VUMPS 197: obj = -1.976927769100e-01 err = 3.7553265539e-01 time = 0.02 sec +[ Info: VUMPS 198: obj = -2.828727739426e-01 err = 3.5757076377e-01 time = 0.03 sec +[ Info: VUMPS 199: obj = -2.312961866796e-01 err = 3.7739801742e-01 time = 0.04 sec +┌ Warning: VUMPS cancel 200: obj = -2.166543836900e-01 err = 3.7105307386e-01 time = 6.55 sec +└ @ MPSKit ~/git/MPSKit.jl/src/algorithms/groundstate/vumps.jl:73 ```` @@ -269,8 +269,8 @@ groundstate, cache, delta = find_groundstate(state, H, GradientGrassmann(; maxit ```` ```` -[ Info: CG: initializing with f = 0.249999715161, ‖∇f‖ = 1.0660e-03 -┌ Warning: CG: not converged to requested tol after 20 iterations and time 4.69 s: f = -0.442407384926, ‖∇f‖ = 8.5132e-03 +[ Info: CG: initializing with f = 0.249995294058, ‖∇f‖ = 4.4866e-03 +┌ Warning: CG: not converged to requested tol after 20 iterations and time 8.67 s: f = -0.441284034826, ‖∇f‖ = 1.1624e-02 └ @ OptimKit ~/.julia/packages/OptimKit/G6i79/src/cg.jl:172 ```` @@ -286,151 +286,143 @@ transferplot(groundstate, groundstate) - + - + - + - + - - + + - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + ``` @@ -469,108 +461,108 @@ groundstate, envs, delta = find_groundstate(state, H2, ```` ```` -[ Info: VUMPS init: obj = +4.991752933537e-01 err = 6.3268e-02 -[ Info: VUMPS 1: obj = -3.177585801487e-01 err = 3.4457002483e-01 time = 0.26 sec -[ Info: VUMPS 2: obj = -8.783277426699e-01 err = 7.2859563338e-02 time = 0.02 sec -[ Info: VUMPS 3: obj = -8.852412984334e-01 err = 1.2557144867e-02 time = 0.01 sec -[ Info: VUMPS 4: obj = -8.859163712655e-01 err = 6.9408088858e-03 time = 0.02 sec -[ Info: VUMPS 5: obj = -8.861199316380e-01 err = 3.8910203290e-03 time = 0.09 sec -[ Info: VUMPS 6: obj = -8.861865553641e-01 err = 2.8480322641e-03 time = 0.02 sec -[ Info: VUMPS 7: obj = -8.862145191872e-01 err = 2.0944450809e-03 time = 0.02 sec -[ Info: VUMPS 8: obj = -8.862266281943e-01 err = 1.7300547846e-03 time = 0.02 sec -[ Info: VUMPS 9: obj = -8.862324798521e-01 err = 1.4474994746e-03 time = 0.04 sec -[ Info: VUMPS 10: obj = -8.862352100866e-01 err = 1.4187844649e-03 time = 0.10 sec -[ Info: VUMPS 11: obj = -8.862366322409e-01 err = 1.3017638159e-03 time = 0.05 sec -[ Info: VUMPS 12: obj = -8.862372899823e-01 err = 1.4040262285e-03 time = 0.05 sec -[ Info: VUMPS 13: obj = -8.862377268299e-01 err = 1.3307584406e-03 time = 0.05 sec -[ Info: VUMPS 14: obj = -8.862379533217e-01 err = 1.4146443881e-03 time = 0.05 sec -[ Info: VUMPS 15: obj = -8.862382038030e-01 err = 1.3272177810e-03 time = 0.09 sec -[ Info: VUMPS 16: obj = -8.862383890237e-01 err = 1.3377148402e-03 time = 0.04 sec -[ Info: VUMPS 17: obj = -8.862386470981e-01 err = 1.1956124852e-03 time = 0.03 sec -[ Info: VUMPS 18: obj = -8.862388507153e-01 err = 1.1254536444e-03 time = 0.03 sec -[ Info: VUMPS 19: obj = -8.862391024178e-01 err = 9.3067565334e-04 time = 0.06 sec -[ Info: VUMPS 20: obj = -8.862392792626e-01 err = 8.0592942148e-04 time = 0.03 sec -[ Info: VUMPS 21: obj = -8.862394408212e-01 err = 6.3751094232e-04 time = 0.03 sec -[ Info: VUMPS 22: obj = -8.862395447014e-01 err = 5.2641943438e-04 time = 0.08 sec -[ Info: VUMPS 23: obj = -8.862396202592e-01 err = 4.2357063364e-04 time = 0.05 sec -[ Info: VUMPS 24: obj = -8.862396677820e-01 err = 3.5346208123e-04 time = 0.05 sec -[ Info: VUMPS 25: obj = -8.862396999787e-01 err = 2.9779376807e-04 time = 0.05 sec -[ Info: VUMPS 26: obj = -8.862397209167e-01 err = 2.5591990952e-04 time = 0.09 sec -[ Info: VUMPS 27: obj = -8.862397356027e-01 err = 2.2322838649e-04 time = 0.05 sec -[ Info: VUMPS 28: obj = -8.862397458518e-01 err = 1.9539015992e-04 time = 0.05 sec -[ Info: VUMPS 29: obj = -8.862397536040e-01 err = 1.7307555129e-04 time = 0.09 sec -[ Info: VUMPS 30: obj = -8.862397594558e-01 err = 1.5260211513e-04 time = 0.05 sec -[ Info: VUMPS 31: obj = -8.862397642025e-01 err = 1.3574708036e-04 time = 0.05 sec -[ Info: VUMPS 32: obj = -8.862397680164e-01 err = 1.2004395690e-04 time = 0.05 sec -[ Info: VUMPS 33: obj = -8.862397712590e-01 err = 1.0674955980e-04 time = 0.11 sec -[ Info: VUMPS 34: obj = -8.862397739766e-01 err = 9.4662033971e-05 time = 0.05 sec -[ Info: VUMPS 35: obj = -8.862397763530e-01 err = 8.4025399031e-05 time = 0.05 sec -[ Info: VUMPS 36: obj = -8.862397784016e-01 err = 7.4857517793e-05 time = 0.09 sec -[ Info: VUMPS 37: obj = -8.862397802247e-01 err = 6.6308156905e-05 time = 0.03 sec -[ Info: VUMPS 38: obj = -8.862397818290e-01 err = 5.9532927397e-05 time = 0.02 sec -[ Info: VUMPS 39: obj = -8.862397832750e-01 err = 5.2635281822e-05 time = 0.03 sec -[ Info: VUMPS 40: obj = -8.862397845683e-01 err = 4.7812478221e-05 time = 0.05 sec -[ Info: VUMPS 41: obj = -8.862397857460e-01 err = 4.2209834688e-05 time = 0.05 sec -[ Info: VUMPS 42: obj = -8.862397868137e-01 err = 3.8962795037e-05 time = 0.04 sec -[ Info: VUMPS 43: obj = -8.862397877944e-01 err = 3.4360884968e-05 time = 0.05 sec -[ Info: VUMPS 44: obj = -8.862397886936e-01 err = 3.2362735927e-05 time = 0.08 sec -[ Info: VUMPS 45: obj = -8.862397895258e-01 err = 2.8521727631e-05 time = 0.05 sec -[ Info: VUMPS 46: obj = -8.862397902957e-01 err = 2.7488241281e-05 time = 0.05 sec -[ Info: VUMPS 47: obj = -8.862397910129e-01 err = 2.4218467594e-05 time = 0.09 sec -[ Info: VUMPS 48: obj = -8.862397916814e-01 err = 2.3904114480e-05 time = 0.05 sec -[ Info: VUMPS 49: obj = -8.862397923073e-01 err = 2.1059493401e-05 time = 0.05 sec -[ Info: VUMPS 50: obj = -8.862397928943e-01 err = 2.1258345276e-05 time = 0.05 sec -[ Info: VUMPS 51: obj = -8.862397934463e-01 err = 1.8738335007e-05 time = 0.07 sec -[ Info: VUMPS 52: obj = -8.862397939664e-01 err = 1.9277045662e-05 time = 0.02 sec -[ Info: VUMPS 53: obj = -8.862397944573e-01 err = 1.7514627759e-05 time = 0.04 sec -[ Info: VUMPS 54: obj = -8.862397949216e-01 err = 1.7756990313e-05 time = 0.08 sec -[ Info: VUMPS 55: obj = -8.862397953612e-01 err = 1.6507778424e-05 time = 0.05 sec -[ Info: VUMPS 56: obj = -8.862397957782e-01 err = 1.6553773913e-05 time = 0.05 sec -[ Info: VUMPS 57: obj = -8.862397961741e-01 err = 1.5646923962e-05 time = 0.05 sec -[ Info: VUMPS 58: obj = -8.862397965506e-01 err = 1.5568425529e-05 time = 0.09 sec -[ Info: VUMPS 59: obj = -8.862397969088e-01 err = 1.4888553209e-05 time = 0.05 sec -[ Info: VUMPS 60: obj = -8.862397972502e-01 err = 1.4734924565e-05 time = 0.05 sec -[ Info: VUMPS 61: obj = -8.862397975756e-01 err = 1.4205866057e-05 time = 0.09 sec -[ Info: VUMPS 62: obj = -8.862397978863e-01 err = 1.4009885396e-05 time = 0.05 sec -[ Info: VUMPS 63: obj = -8.862397981830e-01 err = 1.3582081159e-05 time = 0.05 sec -[ Info: VUMPS 64: obj = -8.862397984666e-01 err = 1.3364902771e-05 time = 0.05 sec -[ Info: VUMPS 65: obj = -8.862397987379e-01 err = 1.3006347133e-05 time = 0.11 sec -[ Info: VUMPS 66: obj = -8.862397989976e-01 err = 1.2781264172e-05 time = 0.04 sec -[ Info: VUMPS 67: obj = -8.862397992464e-01 err = 1.2471321790e-05 time = 0.05 sec -[ Info: VUMPS 68: obj = -8.862397994848e-01 err = 1.2246441659e-05 time = 0.09 sec -[ Info: VUMPS 69: obj = -8.862397997134e-01 err = 1.1971757049e-05 time = 0.03 sec -[ Info: VUMPS 70: obj = -8.862397999328e-01 err = 1.1751836617e-05 time = 0.04 sec -[ Info: VUMPS 71: obj = -8.862398001434e-01 err = 1.1503687228e-05 time = 0.05 sec -[ Info: VUMPS 72: obj = -8.862398003457e-01 err = 1.1291356616e-05 time = 0.08 sec -[ Info: VUMPS 73: obj = -8.862398005401e-01 err = 1.1063958815e-05 time = 0.05 sec -[ Info: VUMPS 74: obj = -8.862398007271e-01 err = 1.0860522804e-05 time = 0.05 sec -[ Info: VUMPS 75: obj = -8.862398009068e-01 err = 1.0649959354e-05 time = 0.05 sec -[ Info: VUMPS 76: obj = -8.862398010799e-01 err = 1.0455912668e-05 time = 0.09 sec -[ Info: VUMPS 77: obj = -8.862398012464e-01 err = 1.0259455807e-05 time = 0.05 sec -[ Info: VUMPS 78: obj = -8.862398014068e-01 err = 1.0074813285e-05 time = 0.04 sec -[ Info: VUMPS 79: obj = -8.862398015614e-01 err = 9.8905015821e-06 time = 0.09 sec -[ Info: VUMPS 80: obj = -8.862398017103e-01 err = 9.7150036409e-06 time = 0.05 sec -[ Info: VUMPS 81: obj = -8.862398018540e-01 err = 9.5413743747e-06 time = 0.05 sec -[ Info: VUMPS 82: obj = -8.862398019925e-01 err = 9.3746148885e-06 time = 0.05 sec -[ Info: VUMPS 83: obj = -8.862398021261e-01 err = 9.2105357844e-06 time = 0.07 sec -[ Info: VUMPS 84: obj = -8.862398022551e-01 err = 9.0520382611e-06 time = 0.04 sec -[ Info: VUMPS 85: obj = -8.862398023797e-01 err = 8.8966029204e-06 time = 0.05 sec -[ Info: VUMPS 86: obj = -8.862398025000e-01 err = 8.7458737419e-06 time = 0.09 sec -[ Info: VUMPS 87: obj = -8.862398026162e-01 err = 8.5983314572e-06 time = 0.05 sec -[ Info: VUMPS 88: obj = -8.862398027285e-01 err = 8.4548765300e-06 time = 0.04 sec -[ Info: VUMPS 89: obj = -8.862398028370e-01 err = 8.3145903294e-06 time = 0.05 sec -[ Info: VUMPS 90: obj = -8.862398029420e-01 err = 8.1779367252e-06 time = 0.11 sec -[ Info: VUMPS 91: obj = -8.862398030435e-01 err = 8.0443547655e-06 time = 0.05 sec -[ Info: VUMPS 92: obj = -8.862398031417e-01 err = 7.9140564546e-06 time = 0.05 sec -[ Info: VUMPS 93: obj = -8.862398032367e-01 err = 7.7866927647e-06 time = 0.09 sec -[ Info: VUMPS 94: obj = -8.862398033287e-01 err = 7.6623336924e-06 time = 0.04 sec -[ Info: VUMPS 95: obj = -8.862398034177e-01 err = 7.5407551359e-06 time = 0.05 sec -[ Info: VUMPS 96: obj = -8.862398035039e-01 err = 7.4219497671e-06 time = 0.05 sec -[ Info: VUMPS 97: obj = -8.862398035875e-01 err = 7.3057668246e-06 time = 0.08 sec -[ Info: VUMPS 98: obj = -8.862398036684e-01 err = 7.1921592390e-06 time = 0.05 sec -[ Info: VUMPS 99: obj = -8.862398037467e-01 err = 7.0810192783e-06 time = 0.05 sec -┌ Warning: VUMPS cancel 100: obj = -8.862398038227e-01 err = 6.9722815663e-06 time = 5.64 sec -└ @ MPSKit ~/Projects/MPSKit.jl/src/algorithms/groundstate/vumps.jl:73 +[ Info: VUMPS init: obj = +4.984019701133e-01 err = 8.5802e-02 +[ Info: VUMPS 1: obj = -1.587916149322e-01 err = 3.0275893671e-01 time = 0.04 sec +[ Info: VUMPS 2: obj = -8.673455076862e-01 err = 1.2709823901e-01 time = 0.03 sec +[ Info: VUMPS 3: obj = -8.851574658168e-01 err = 1.3169797608e-02 time = 0.03 sec +[ Info: VUMPS 4: obj = -8.859091848496e-01 err = 6.4392637469e-03 time = 0.10 sec +[ Info: VUMPS 5: obj = -8.861109088796e-01 err = 4.0495543537e-03 time = 0.03 sec +[ Info: VUMPS 6: obj = -8.861799998365e-01 err = 3.2423747666e-03 time = 0.04 sec +[ Info: VUMPS 7: obj = -8.862098845534e-01 err = 2.5878808961e-03 time = 0.04 sec +[ Info: VUMPS 8: obj = -8.862231663400e-01 err = 2.2980643190e-03 time = 0.04 sec +[ Info: VUMPS 9: obj = -8.862302410094e-01 err = 1.9030792239e-03 time = 0.04 sec +[ Info: VUMPS 10: obj = -8.862337476758e-01 err = 1.5771139900e-03 time = 0.04 sec +[ Info: VUMPS 11: obj = -8.862357477670e-01 err = 1.3386319635e-03 time = 0.05 sec +[ Info: VUMPS 12: obj = -8.862368077969e-01 err = 1.1118556991e-03 time = 0.04 sec +[ Info: VUMPS 13: obj = -8.862374143689e-01 err = 9.6368190823e-04 time = 0.05 sec +[ Info: VUMPS 14: obj = -8.862377772951e-01 err = 8.9527621802e-04 time = 0.04 sec +[ Info: VUMPS 15: obj = -8.862380175272e-01 err = 8.0363531590e-04 time = 0.05 sec +[ Info: VUMPS 16: obj = -8.862381809910e-01 err = 8.1310379974e-04 time = 0.05 sec +[ Info: VUMPS 17: obj = -8.862383289606e-01 err = 7.6372744918e-04 time = 0.04 sec +[ Info: VUMPS 18: obj = -8.862384512228e-01 err = 7.7168418781e-04 time = 0.05 sec +[ Info: VUMPS 19: obj = -8.862385939943e-01 err = 7.2169619957e-04 time = 0.05 sec +[ Info: VUMPS 20: obj = -8.862387296086e-01 err = 6.9847750143e-04 time = 0.14 sec +[ Info: VUMPS 21: obj = -8.862388823501e-01 err = 6.4519851323e-04 time = 0.08 sec +[ Info: VUMPS 22: obj = -8.862390299430e-01 err = 5.9543430924e-04 time = 0.08 sec +[ Info: VUMPS 23: obj = -8.862391796049e-01 err = 5.2965230331e-04 time = 0.08 sec +[ Info: VUMPS 24: obj = -8.862393105628e-01 err = 4.6746973021e-04 time = 0.08 sec +[ Info: VUMPS 25: obj = -8.862394256169e-01 err = 4.0012165330e-04 time = 0.08 sec +[ Info: VUMPS 26: obj = -8.862395146213e-01 err = 3.4348612692e-04 time = 0.05 sec +[ Info: VUMPS 27: obj = -8.862395831043e-01 err = 2.9133403018e-04 time = 0.05 sec +[ Info: VUMPS 28: obj = -8.862396319943e-01 err = 2.5002069641e-04 time = 0.05 sec +[ Info: VUMPS 29: obj = -8.862396671228e-01 err = 2.1555353820e-04 time = 0.05 sec +[ Info: VUMPS 30: obj = -8.862396917938e-01 err = 1.8787657127e-04 time = 0.05 sec +[ Info: VUMPS 31: obj = -8.862397096146e-01 err = 1.6535795978e-04 time = 0.05 sec +[ Info: VUMPS 32: obj = -8.862397226287e-01 err = 1.4647959939e-04 time = 0.04 sec +[ Info: VUMPS 33: obj = -8.862397325704e-01 err = 1.3058029797e-04 time = 0.04 sec +[ Info: VUMPS 34: obj = -8.862397403210e-01 err = 1.1703644593e-04 time = 0.05 sec +[ Info: VUMPS 35: obj = -8.862397466535e-01 err = 1.0479443543e-04 time = 0.05 sec +[ Info: VUMPS 36: obj = -8.862397518883e-01 err = 9.4757697496e-05 time = 0.07 sec +[ Info: VUMPS 37: obj = -8.862397563857e-01 err = 8.4793760071e-05 time = 0.05 sec +[ Info: VUMPS 38: obj = -8.862397602492e-01 err = 7.7330476046e-05 time = 0.05 sec +[ Info: VUMPS 39: obj = -8.862397636664e-01 err = 6.9007080884e-05 time = 0.05 sec +[ Info: VUMPS 40: obj = -8.862397666688e-01 err = 6.3574783777e-05 time = 0.05 sec +[ Info: VUMPS 41: obj = -8.862397693669e-01 err = 5.6541518218e-05 time = 0.05 sec +[ Info: VUMPS 42: obj = -8.862397717706e-01 err = 5.2754788593e-05 time = 0.05 sec +[ Info: VUMPS 43: obj = -8.862397739510e-01 err = 4.6775683210e-05 time = 0.05 sec +[ Info: VUMPS 44: obj = -8.862397759130e-01 err = 4.4320904414e-05 time = 0.05 sec +[ Info: VUMPS 45: obj = -8.862397777047e-01 err = 3.9206640444e-05 time = 0.05 sec +[ Info: VUMPS 46: obj = -8.862397793305e-01 err = 3.7808257666e-05 time = 0.05 sec +[ Info: VUMPS 47: obj = -8.862397808237e-01 err = 3.3398078106e-05 time = 0.05 sec +[ Info: VUMPS 48: obj = -8.862397821887e-01 err = 3.2810027183e-05 time = 0.05 sec +[ Info: VUMPS 49: obj = -8.862397834489e-01 err = 2.8968297396e-05 time = 0.05 sec +[ Info: VUMPS 50: obj = -8.862397846088e-01 err = 2.8975455638e-05 time = 0.07 sec +[ Info: VUMPS 51: obj = -8.862397856847e-01 err = 2.5588978527e-05 time = 0.05 sec +[ Info: VUMPS 52: obj = -8.862397866809e-01 err = 2.6012564811e-05 time = 0.05 sec +[ Info: VUMPS 53: obj = -8.862397876092e-01 err = 2.3457901855e-05 time = 0.05 sec +[ Info: VUMPS 54: obj = -8.862397884732e-01 err = 2.3689215433e-05 time = 0.05 sec +[ Info: VUMPS 55: obj = -8.862397892815e-01 err = 2.1824090291e-05 time = 0.05 sec +[ Info: VUMPS 56: obj = -8.862397900373e-01 err = 2.1828940191e-05 time = 0.05 sec +[ Info: VUMPS 57: obj = -8.862397907469e-01 err = 2.0446827272e-05 time = 0.05 sec +[ Info: VUMPS 58: obj = -8.862397914131e-01 err = 2.0302723176e-05 time = 0.05 sec +[ Info: VUMPS 59: obj = -8.862397920405e-01 err = 1.9254981490e-05 time = 0.05 sec +[ Info: VUMPS 60: obj = -8.862397926316e-01 err = 1.9019089435e-05 time = 0.05 sec +[ Info: VUMPS 61: obj = -8.862397931898e-01 err = 1.8202179915e-05 time = 0.05 sec +[ Info: VUMPS 62: obj = -8.862397937172e-01 err = 1.7914246472e-05 time = 0.05 sec +[ Info: VUMPS 63: obj = -8.862397942165e-01 err = 1.7257569224e-05 time = 0.05 sec +[ Info: VUMPS 64: obj = -8.862397946895e-01 err = 1.6944227459e-05 time = 0.05 sec +[ Info: VUMPS 65: obj = -8.862397951382e-01 err = 1.6400243806e-05 time = 0.07 sec +[ Info: VUMPS 66: obj = -8.862397955643e-01 err = 1.6078941938e-05 time = 0.05 sec +[ Info: VUMPS 67: obj = -8.862397959693e-01 err = 1.5615639391e-05 time = 0.05 sec +[ Info: VUMPS 68: obj = -8.862397963546e-01 err = 1.5297227747e-05 time = 0.05 sec +[ Info: VUMPS 69: obj = -8.862397967215e-01 err = 1.4893178345e-05 time = 0.05 sec +[ Info: VUMPS 70: obj = -8.862397970712e-01 err = 1.4584085106e-05 time = 0.05 sec +[ Info: VUMPS 71: obj = -8.862397974047e-01 err = 1.4224856670e-05 time = 0.05 sec +[ Info: VUMPS 72: obj = -8.862397977230e-01 err = 1.3928611175e-05 time = 0.05 sec +[ Info: VUMPS 73: obj = -8.862397980271e-01 err = 1.3604370840e-05 time = 0.05 sec +[ Info: VUMPS 74: obj = -8.862397983177e-01 err = 1.3322648086e-05 time = 0.05 sec +[ Info: VUMPS 75: obj = -8.862397985958e-01 err = 1.3026579806e-05 time = 0.05 sec +[ Info: VUMPS 76: obj = -8.862397988619e-01 err = 1.2759903586e-05 time = 0.05 sec +[ Info: VUMPS 77: obj = -8.862397991167e-01 err = 1.2487175232e-05 time = 0.05 sec +[ Info: VUMPS 78: obj = -8.862397993609e-01 err = 1.2235379085e-05 time = 0.05 sec +[ Info: VUMPS 79: obj = -8.862397995950e-01 err = 1.1982470825e-05 time = 0.05 sec +[ Info: VUMPS 80: obj = -8.862397998196e-01 err = 1.1744996531e-05 time = 0.07 sec +[ Info: VUMPS 81: obj = -8.862398000351e-01 err = 1.1509266452e-05 time = 0.05 sec +[ Info: VUMPS 82: obj = -8.862398002421e-01 err = 1.1285351075e-05 time = 0.05 sec +[ Info: VUMPS 83: obj = -8.862398004410e-01 err = 1.1064753788e-05 time = 0.05 sec +[ Info: VUMPS 84: obj = -8.862398006321e-01 err = 1.0853545836e-05 time = 0.05 sec +[ Info: VUMPS 85: obj = -8.862398008159e-01 err = 1.0646449681e-05 time = 0.05 sec +[ Info: VUMPS 86: obj = -8.862398009928e-01 err = 1.0447078851e-05 time = 0.05 sec +[ Info: VUMPS 87: obj = -8.862398011630e-01 err = 1.0252146279e-05 time = 0.05 sec +[ Info: VUMPS 88: obj = -8.862398013268e-01 err = 1.0063763685e-05 time = 0.05 sec +[ Info: VUMPS 89: obj = -8.862398014847e-01 err = 9.8798722021e-06 time = 0.05 sec +[ Info: VUMPS 90: obj = -8.862398016367e-01 err = 9.7016721333e-06 time = 0.04 sec +[ Info: VUMPS 91: obj = -8.862398017834e-01 err = 9.5278612850e-06 time = 0.04 sec +[ Info: VUMPS 92: obj = -8.862398019247e-01 err = 9.3590916174e-06 time = 0.04 sec +[ Info: VUMPS 93: obj = -8.862398020611e-01 err = 9.1945267654e-06 time = 0.04 sec +[ Info: VUMPS 94: obj = -8.862398021927e-01 err = 9.0344925557e-06 time = 0.05 sec +[ Info: VUMPS 95: obj = -8.862398023197e-01 err = 8.8784395519e-06 time = 0.07 sec +[ Info: VUMPS 96: obj = -8.862398024423e-01 err = 8.7265026811e-06 time = 0.05 sec +[ Info: VUMPS 97: obj = -8.862398025608e-01 err = 8.5783096813e-06 time = 0.05 sec +[ Info: VUMPS 98: obj = -8.862398026752e-01 err = 8.4338863245e-06 time = 0.05 sec +[ Info: VUMPS 99: obj = -8.862398027858e-01 err = 8.2929703487e-06 time = 0.05 sec +┌ Warning: VUMPS cancel 100: obj = -8.862398028927e-01 err = 8.1555273601e-06 time = 5.09 sec +└ @ MPSKit ~/git/MPSKit.jl/src/algorithms/groundstate/vumps.jl:73 ```` @@ -588,86 +580,80 @@ entanglementplot(groundstate) - + - + - + - + - - + + - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + ``` @@ -706,7 +692,7 @@ state = InfiniteMPS([P, P], [V1, V2]); ```` ┌ Warning: Constructing an MPS from tensors that are not full rank -└ @ MPSKit ~/Projects/MPSKit.jl/src/states/infinitemps.jl:149 +└ @ MPSKit ~/git/MPSKit.jl/src/states/infinitemps.jl:149 ```` @@ -722,74 +708,76 @@ groundstate, cache, delta = find_groundstate(state, H2, ```` 52 70 -[ Info: VUMPS init: obj = +1.757765749207e-02 err = 3.4286e-01 -[ Info: VUMPS 1: obj = -8.696634007330e-01 err = 1.3207614531e-01 time = 0.21 sec -[ Info: VUMPS 2: obj = -8.857455218686e-01 err = 7.9112939847e-03 time = 0.12 sec -[ Info: VUMPS 3: obj = -8.861028581917e-01 err = 3.6033827124e-03 time = 0.03 sec -[ Info: VUMPS 4: obj = -8.862158212696e-01 err = 1.6926177755e-03 time = 0.03 sec -[ Info: VUMPS 5: obj = -8.862577312672e-01 err = 1.0611468264e-03 time = 0.03 sec -[ Info: VUMPS 6: obj = -8.862746027064e-01 err = 7.7042891777e-04 time = 0.03 sec -[ Info: VUMPS 7: obj = -8.862817300286e-01 err = 6.2709123174e-04 time = 0.04 sec -[ Info: VUMPS 8: obj = -8.862849391575e-01 err = 5.0493868777e-04 time = 0.05 sec -[ Info: VUMPS 9: obj = -8.862864498235e-01 err = 3.9930999383e-04 time = 0.03 sec -[ Info: VUMPS 10: obj = -8.862871831010e-01 err = 3.0800857369e-04 time = 0.03 sec -[ Info: VUMPS 11: obj = -8.862875441594e-01 err = 2.3263846933e-04 time = 0.05 sec -[ Info: VUMPS 12: obj = -8.862877229110e-01 err = 1.7276003533e-04 time = 0.10 sec -[ Info: VUMPS 13: obj = -8.862878114835e-01 err = 1.2669078880e-04 time = 0.05 sec -[ Info: VUMPS 14: obj = -8.862878553110e-01 err = 9.2023000727e-05 time = 0.13 sec -[ Info: VUMPS 15: obj = -8.862878769711e-01 err = 6.6386757063e-05 time = 0.05 sec -[ Info: VUMPS 16: obj = -8.862878876657e-01 err = 4.7658401917e-05 time = 0.05 sec -[ Info: VUMPS 17: obj = -8.862878929424e-01 err = 3.4092229469e-05 time = 0.05 sec -[ Info: VUMPS 18: obj = -8.862878955462e-01 err = 2.4325841163e-05 time = 0.05 sec -[ Info: VUMPS 19: obj = -8.862878968315e-01 err = 1.7324877324e-05 time = 0.06 sec -[ Info: VUMPS 20: obj = -8.862878974663e-01 err = 1.2323307122e-05 time = 0.06 sec -[ Info: VUMPS 21: obj = -8.862878977801e-01 err = 8.7563035932e-06 time = 0.05 sec -[ Info: VUMPS 22: obj = -8.862878979353e-01 err = 6.2166861423e-06 time = 0.04 sec -[ Info: VUMPS 23: obj = -8.862878980121e-01 err = 4.4107365461e-06 time = 0.04 sec -[ Info: VUMPS 24: obj = -8.862878980502e-01 err = 3.1277152106e-06 time = 0.04 sec -[ Info: VUMPS 25: obj = -8.862878980690e-01 err = 2.2170142160e-06 time = 0.11 sec -[ Info: VUMPS 26: obj = -8.862878980784e-01 err = 1.5707437234e-06 time = 0.05 sec -[ Info: VUMPS 27: obj = -8.862878980830e-01 err = 1.1124442525e-06 time = 0.04 sec -[ Info: VUMPS 28: obj = -8.862878980853e-01 err = 7.8759729582e-07 time = 0.05 sec -[ Info: VUMPS 29: obj = -8.862878980865e-01 err = 5.5743614600e-07 time = 0.06 sec -[ Info: VUMPS 30: obj = -8.862878980870e-01 err = 3.9442227022e-07 time = 0.05 sec -[ Info: VUMPS 31: obj = -8.862878980873e-01 err = 2.7900505779e-07 time = 0.04 sec -[ Info: VUMPS 32: obj = -8.862878980875e-01 err = 1.9731276183e-07 time = 0.06 sec -[ Info: VUMPS 33: obj = -8.862878980875e-01 err = 1.3950769360e-07 time = 0.06 sec -[ Info: VUMPS 34: obj = -8.862878980876e-01 err = 9.8611263682e-08 time = 0.11 sec -[ Info: VUMPS 35: obj = -8.862878980876e-01 err = 6.9693039590e-08 time = 0.05 sec -[ Info: VUMPS 36: obj = -8.862878980876e-01 err = 4.9246283392e-08 time = 0.05 sec -[ Info: VUMPS 37: obj = -8.862878980876e-01 err = 3.4792297123e-08 time = 0.05 sec -[ Info: VUMPS 38: obj = -8.862878980876e-01 err = 2.4576688585e-08 time = 0.04 sec -[ Info: VUMPS 39: obj = -8.862878980876e-01 err = 1.7357975423e-08 time = 0.05 sec -[ Info: VUMPS 40: obj = -8.862878980877e-01 err = 1.2257865533e-08 time = 0.04 sec -[ Info: VUMPS 41: obj = -8.862878980877e-01 err = 8.6551624916e-09 time = 0.06 sec -[ Info: VUMPS 42: obj = -8.862878980877e-01 err = 6.1106026828e-09 time = 0.08 sec -[ Info: VUMPS 43: obj = -8.862878980877e-01 err = 4.3136802051e-09 time = 0.06 sec -[ Info: VUMPS 44: obj = -8.862878980877e-01 err = 3.0449401934e-09 time = 0.05 sec -[ Info: VUMPS 45: obj = -8.862878980877e-01 err = 2.1490894708e-09 time = 0.06 sec -[ Info: VUMPS 46: obj = -8.862878980877e-01 err = 1.5166673693e-09 time = 0.06 sec -[ Info: VUMPS 47: obj = -8.862878980877e-01 err = 1.0702667096e-09 time = 0.06 sec -[ Info: VUMPS 48: obj = -8.862878980877e-01 err = 7.5519801704e-10 time = 0.05 sec -[ Info: VUMPS 49: obj = -8.862878980877e-01 err = 5.3284409175e-10 time = 0.08 sec -[ Info: VUMPS 50: obj = -8.862878980877e-01 err = 3.7593665457e-10 time = 0.05 sec -[ Info: VUMPS 51: obj = -8.862878980877e-01 err = 2.6521797987e-10 time = 0.04 sec -[ Info: VUMPS 52: obj = -8.862878980877e-01 err = 1.8709489707e-10 time = 0.05 sec -[ Info: VUMPS 53: obj = -8.862878980877e-01 err = 1.3197998674e-10 time = 0.05 sec -[ Info: VUMPS 54: obj = -8.862878980877e-01 err = 9.3098181427e-11 time = 0.05 sec -[ Info: VUMPS 55: obj = -8.862878980877e-01 err = 6.5667849471e-11 time = 0.05 sec -[ Info: VUMPS 56: obj = -8.862878980877e-01 err = 4.6315648967e-11 time = 0.04 sec -[ Info: VUMPS 57: obj = -8.862878980877e-01 err = 3.2665823891e-11 time = 0.10 sec -[ Info: VUMPS 58: obj = -8.862878980878e-01 err = 2.3038786394e-11 time = 0.03 sec -[ Info: VUMPS 59: obj = -8.862878980878e-01 err = 1.6249291109e-11 time = 0.05 sec -[ Info: VUMPS 60: obj = -8.862878980878e-01 err = 1.1461841113e-11 time = 0.04 sec -[ Info: VUMPS 61: obj = -8.862878980878e-01 err = 8.0811426926e-12 time = 0.03 sec -[ Info: VUMPS 62: obj = -8.862878980878e-01 err = 5.6989473760e-12 time = 0.03 sec -[ Info: VUMPS 63: obj = -8.862878980878e-01 err = 4.0199684636e-12 time = 0.02 sec -[ Info: VUMPS 64: obj = -8.862878980878e-01 err = 2.8323491405e-12 time = 0.03 sec -[ Info: VUMPS 65: obj = -8.862878980878e-01 err = 1.9985430493e-12 time = 0.03 sec -[ Info: VUMPS 66: obj = -8.862878980878e-01 err = 1.4098856438e-12 time = 0.07 sec -[ Info: VUMPS conv 67: obj = -8.862878980878e-01 err = 9.9583426980e-13 time = 3.67 sec +[ Info: VUMPS init: obj = +2.876758033543e-02 err = 4.0644e-01 +[ Info: VUMPS 1: obj = -8.737906665472e-01 err = 1.2045339353e-01 time = 0.02 sec +[ Info: VUMPS 2: obj = -8.857246450469e-01 err = 7.4507510196e-03 time = 0.02 sec +[ Info: VUMPS 3: obj = -8.861168843924e-01 err = 3.8806594310e-03 time = 0.02 sec +[ Info: VUMPS 4: obj = -8.862224229828e-01 err = 1.8866990312e-03 time = 0.02 sec +[ Info: VUMPS 5: obj = -8.862604919935e-01 err = 9.5632297117e-04 time = 0.03 sec +[ Info: VUMPS 6: obj = -8.862756622261e-01 err = 6.9689351881e-04 time = 0.03 sec +[ Info: VUMPS 7: obj = -8.862821568604e-01 err = 5.4384343431e-04 time = 0.03 sec +[ Info: VUMPS 8: obj = -8.862851144799e-01 err = 4.3498937171e-04 time = 0.03 sec +[ Info: VUMPS 9: obj = -8.862865208255e-01 err = 3.4022284480e-04 time = 0.04 sec +[ Info: VUMPS 10: obj = -8.862872092893e-01 err = 2.6155639792e-04 time = 0.04 sec +[ Info: VUMPS 11: obj = -8.862875517233e-01 err = 1.9816403910e-04 time = 0.04 sec +[ Info: VUMPS 12: obj = -8.862877233984e-01 err = 1.4833091166e-04 time = 0.04 sec +[ Info: VUMPS 13: obj = -8.862878098491e-01 err = 1.1002989377e-04 time = 0.04 sec +[ Info: VUMPS 14: obj = -8.862878534672e-01 err = 8.1031399491e-05 time = 0.04 sec +[ Info: VUMPS 15: obj = -8.862878754996e-01 err = 5.9348893102e-05 time = 0.04 sec +[ Info: VUMPS 16: obj = -8.862878866404e-01 err = 4.3288174232e-05 time = 0.05 sec +[ Info: VUMPS 17: obj = -8.862878922801e-01 err = 3.1476890398e-05 time = 0.05 sec +[ Info: VUMPS 18: obj = -8.862878951379e-01 err = 2.2833727293e-05 time = 0.04 sec +[ Info: VUMPS 19: obj = -8.862878965878e-01 err = 1.6534345712e-05 time = 0.04 sec +[ Info: VUMPS 20: obj = -8.862878973242e-01 err = 1.1956033801e-05 time = 0.04 sec +[ Info: VUMPS 21: obj = -8.862878976987e-01 err = 8.6354344118e-06 time = 0.04 sec +[ Info: VUMPS 22: obj = -8.862878978893e-01 err = 6.2316368830e-06 time = 0.14 sec +[ Info: VUMPS 23: obj = -8.862878979864e-01 err = 4.4936852018e-06 time = 0.04 sec +[ Info: VUMPS 24: obj = -8.862878980359e-01 err = 3.2380407694e-06 time = 0.04 sec +[ Info: VUMPS 25: obj = -8.862878980612e-01 err = 2.3317908791e-06 time = 0.05 sec +[ Info: VUMPS 26: obj = -8.862878980741e-01 err = 1.6782796290e-06 time = 0.04 sec +[ Info: VUMPS 27: obj = -8.862878980807e-01 err = 1.2073259818e-06 time = 0.04 sec +[ Info: VUMPS 28: obj = -8.862878980841e-01 err = 8.6818060956e-07 time = 0.04 sec +[ Info: VUMPS 29: obj = -8.862878980858e-01 err = 6.2400887691e-07 time = 0.04 sec +[ Info: VUMPS 30: obj = -8.862878980867e-01 err = 4.4833124713e-07 time = 0.04 sec +[ Info: VUMPS 31: obj = -8.862878980871e-01 err = 3.2198399757e-07 time = 0.04 sec +[ Info: VUMPS 32: obj = -8.862878980874e-01 err = 2.3117466228e-07 time = 0.04 sec +[ Info: VUMPS 33: obj = -8.862878980875e-01 err = 1.6592578229e-07 time = 0.04 sec +[ Info: VUMPS 34: obj = -8.862878980876e-01 err = 1.1905960048e-07 time = 0.04 sec +[ Info: VUMPS 35: obj = -8.862878980876e-01 err = 8.5408410847e-08 time = 0.04 sec +[ Info: VUMPS 36: obj = -8.862878980876e-01 err = 6.1253530320e-08 time = 0.04 sec +[ Info: VUMPS 37: obj = -8.862878980876e-01 err = 4.3920166462e-08 time = 0.04 sec +[ Info: VUMPS 38: obj = -8.862878980877e-01 err = 3.1485556009e-08 time = 0.04 sec +[ Info: VUMPS 39: obj = -8.862878980877e-01 err = 2.2566832411e-08 time = 0.04 sec +[ Info: VUMPS 40: obj = -8.862878980877e-01 err = 1.6171599796e-08 time = 0.04 sec +[ Info: VUMPS 41: obj = -8.862878980877e-01 err = 1.1586823278e-08 time = 0.04 sec +[ Info: VUMPS 42: obj = -8.862878980877e-01 err = 8.3006991909e-09 time = 0.04 sec +[ Info: VUMPS 43: obj = -8.862878980877e-01 err = 5.9456786414e-09 time = 0.04 sec +[ Info: VUMPS 44: obj = -8.862878980877e-01 err = 4.2582593636e-09 time = 0.10 sec +[ Info: VUMPS 45: obj = -8.862878980877e-01 err = 3.0493882526e-09 time = 0.04 sec +[ Info: VUMPS 46: obj = -8.862878980877e-01 err = 2.1834676831e-09 time = 0.04 sec +[ Info: VUMPS 47: obj = -8.862878980877e-01 err = 1.5632717223e-09 time = 0.04 sec +[ Info: VUMPS 48: obj = -8.862878980877e-01 err = 1.1191534669e-09 time = 0.04 sec +[ Info: VUMPS 49: obj = -8.862878980877e-01 err = 8.0114053028e-10 time = 0.04 sec +[ Info: VUMPS 50: obj = -8.862878980877e-01 err = 5.7344926629e-10 time = 0.04 sec +[ Info: VUMPS 51: obj = -8.862878980877e-01 err = 4.1044557276e-10 time = 0.04 sec +[ Info: VUMPS 52: obj = -8.862878980877e-01 err = 2.9375533030e-10 time = 0.04 sec +[ Info: VUMPS 53: obj = -8.862878980877e-01 err = 2.1022890604e-10 time = 0.04 sec +[ Info: VUMPS 54: obj = -8.862878980877e-01 err = 1.5044585215e-10 time = 0.04 sec +[ Info: VUMPS 55: obj = -8.862878980877e-01 err = 1.0765471631e-10 time = 0.04 sec +[ Info: VUMPS 56: obj = -8.862878980878e-01 err = 7.7031624881e-11 time = 0.04 sec +[ Info: VUMPS 57: obj = -8.862878980878e-01 err = 5.5116540056e-11 time = 0.04 sec +[ Info: VUMPS 58: obj = -8.862878980878e-01 err = 3.9433877868e-11 time = 0.04 sec +[ Info: VUMPS 59: obj = -8.862878980878e-01 err = 2.8212244668e-11 time = 0.04 sec +[ Info: VUMPS 60: obj = -8.862878980878e-01 err = 2.0184013518e-11 time = 0.03 sec +[ Info: VUMPS 61: obj = -8.862878980878e-01 err = 1.4441829469e-11 time = 0.03 sec +[ Info: VUMPS 62: obj = -8.862878980878e-01 err = 1.0334072641e-11 time = 0.03 sec +[ Info: VUMPS 63: obj = -8.862878980878e-01 err = 7.3919424370e-12 time = 0.03 sec +[ Info: VUMPS 64: obj = -8.862878980878e-01 err = 5.2875140115e-12 time = 0.03 sec +[ Info: VUMPS 65: obj = -8.862878980878e-01 err = 3.7817587214e-12 time = 0.03 sec +[ Info: VUMPS 66: obj = -8.862878980878e-01 err = 2.7051025623e-12 time = 0.03 sec +[ Info: VUMPS 67: obj = -8.862878980878e-01 err = 1.9384032618e-12 time = 0.03 sec +[ Info: VUMPS 68: obj = -8.862878980878e-01 err = 1.3868042929e-12 time = 0.06 sec +[ Info: VUMPS conv 69: obj = -8.862878980878e-01 err = 9.8790427735e-13 time = 2.76 sec ```` diff --git a/docs/src/examples/quantum1d/4.xxz-heisenberg/main.ipynb b/docs/src/examples/quantum1d/4.xxz-heisenberg/main.ipynb index 35adb4c49..0fecfb03c 100644 --- a/docs/src/examples/quantum1d/4.xxz-heisenberg/main.ipynb +++ b/docs/src/examples/quantum1d/4.xxz-heisenberg/main.ipynb @@ -267,11 +267,11 @@ "file_extension": ".jl", "mimetype": "application/julia", "name": "julia", - "version": "1.11.4" + "version": "1.10.4" }, "kernelspec": { - "name": "julia-1.11", - "display_name": "Julia 1.11.4", + "name": "julia-1.10", + "display_name": "Julia 1.10.4", "language": "julia" } }, diff --git a/docs/src/examples/quantum1d/5.haldane-spt/index.md b/docs/src/examples/quantum1d/5.haldane-spt/index.md index 2d653e563..c54642d41 100644 --- a/docs/src/examples/quantum1d/5.haldane-spt/index.md +++ b/docs/src/examples/quantum1d/5.haldane-spt/index.md @@ -104,276 +104,274 @@ transferplot(ψ; sectors, title="Transfer matrix spectrum", legend=:outertop) - + - + - + - + - - + + - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + ``` @@ -401,7 +399,7 @@ E_plus = expectation_value(ψ_plus, H) ```` ```` --1.4014193313393013 + 2.816079042392764e-17im +-1.4014193313393015 - 2.5829324229205153e-17im ```` ````julia @@ -412,7 +410,7 @@ E_minus = expectation_value(ψ_minus, H) ```` ```` --1.4014839739630858 - 3.028001148219207e-17im +-1.4014839739630842 - 3.3657493938637046e-16im ```` ````julia @@ -427,468 +425,464 @@ plot(transferp_plus, transferp_minus; layout=(1, 2), size=(800, 400)) - + - + - + - + - - + + - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + - - + + - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + ``` @@ -902,104 +896,104 @@ plot(entanglementp_plus, entanglementp_minus; layout=(1, 2), size=(800, 400)) - + - + - + - + - - + + - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + - - + + - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + ``` @@ -1025,8 +1019,8 @@ println("S_plus = $S_plus") ```` ```` -S_minus + log(2) = 1.5486227235400487 -S_plus = 1.545032353059228 +S_minus + log(2) = 1.5486227235443117 +S_plus = 1.545032353057881 ```` diff --git a/docs/src/examples/quantum1d/5.haldane-spt/main.ipynb b/docs/src/examples/quantum1d/5.haldane-spt/main.ipynb index 6ee2e5fc6..567c463e6 100644 --- a/docs/src/examples/quantum1d/5.haldane-spt/main.ipynb +++ b/docs/src/examples/quantum1d/5.haldane-spt/main.ipynb @@ -198,11 +198,11 @@ "file_extension": ".jl", "mimetype": "application/julia", "name": "julia", - "version": "1.11.4" + "version": "1.10.4" }, "kernelspec": { - "name": "julia-1.11", - "display_name": "Julia 1.11.4", + "name": "julia-1.10", + "display_name": "Julia 1.10.4", "language": "julia" } }, diff --git a/docs/src/examples/quantum1d/6.hubbard/index.md b/docs/src/examples/quantum1d/6.hubbard/index.md index d503d15c6..d05b363bd 100644 --- a/docs/src/examples/quantum1d/6.hubbard/index.md +++ b/docs/src/examples/quantum1d/6.hubbard/index.md @@ -108,57 +108,54 @@ Groundstate energy: ```` ```` -[ Info: VUMPS init: obj = -1.459840683100e+00 err = 5.6981e-01 -[ Info: VUMPS conv 12: obj = -4.377110142621e+00 err = 9.8483052508e-03 time = 10.30 sec -[ Info: VUMPS init: obj = -4.377110142621e+00 err = 1.6381e-02 -┌ Warning: VUMPS cancel 10: obj = -4.378360604966e+00 err = 1.3577487013e-02 time = 0.32 sec -└ @ MPSKit ~/Projects/MPSKit.jl/src/algorithms/groundstate/vumps.jl:73 -[ Info: VUMPS init: obj = -4.378360604966e+00 err = 1.4940e-02 -┌ Warning: VUMPS cancel 10: obj = -4.378841582201e+00 err = 1.2373568873e-02 time = 0.41 sec -└ @ MPSKit ~/Projects/MPSKit.jl/src/algorithms/groundstate/vumps.jl:73 -[ Info: VUMPS init: obj = -4.378841582201e+00 err = 1.4876e-02 -[ Info: VUMPS conv 5: obj = -4.379480020439e+00 err = 1.8699990778e-04 time = 0.31 sec -[ Info: VUMPS init: obj = -4.379480020439e+00 err = 5.1669e-03 -[ Info: VUMPS conv 5: obj = -4.379639700173e+00 err = 1.9060279760e-04 time = 0.38 sec -[ Info: VUMPS init: obj = -4.379639700173e+00 err = 4.7624e-03 -[ Info: VUMPS conv 4: obj = -4.379745769831e+00 err = 1.8370454603e-04 time = 0.69 sec -[ Info: VUMPS init: obj = -4.379745769831e+00 err = 3.5050e-03 -[ Info: VUMPS conv 3: obj = -4.379799852985e+00 err = 1.4079134809e-04 time = 0.57 sec -[ Info: VUMPS init: obj = -4.379799852985e+00 err = 2.8030e-03 -[ Info: VUMPS conv 2: obj = -4.379838075265e+00 err = 1.7047756792e-04 time = 0.43 sec -[ Info: VUMPS init: obj = -4.379838075265e+00 err = 2.7778e-03 -┌ Warning: VUMPS cancel 10: obj = -4.379884982067e+00 err = 2.6174773183e-04 time = 3.89 sec -└ @ MPSKit ~/Projects/MPSKit.jl/src/algorithms/groundstate/vumps.jl:73 -[ Info: VUMPS init: obj = -4.379884982067e+00 err = 2.8617e-03 -[ Info: VUMPS conv 2: obj = -4.379925411329e+00 err = 1.7224103315e-04 time = 0.79 sec -[ Info: VUMPS init: obj = -4.379925411329e+00 err = 2.7213e-03 -[ Info: VUMPS conv 3: obj = -4.379967808299e+00 err = 1.8697039936e-04 time = 1.64 sec -[ Info: VUMPS init: obj = -4.379967808299e+00 err = 1.7885e-03 -[ Info: VUMPS conv 2: obj = -4.379986798935e+00 err = 1.8550344471e-04 time = 1.02 sec -[ Info: VUMPS init: obj = -4.379986798935e+00 err = 1.5831e-03 -[ Info: VUMPS conv 2: obj = -4.380000986838e+00 err = 1.8747241571e-04 time = 1.23 sec -[ Info: VUMPS init: obj = -4.380000986838e+00 err = 1.5082e-03 -[ Info: VUMPS conv 2: obj = -4.380013146139e+00 err = 1.5036056628e-04 time = 2.00 sec -[ Info: VUMPS init: obj = -4.380013146139e+00 err = 1.4240e-03 -[ Info: VUMPS conv 2: obj = -4.380024374442e+00 err = 1.8015023821e-04 time = 2.02 sec -[ Info: VUMPS init: obj = -4.380024374442e+00 err = 1.3383e-03 -[ Info: VUMPS conv 2: obj = -4.380038125427e+00 err = 1.6836301386e-04 time = 2.50 sec -[ Info: VUMPS init: obj = -4.380038125427e+00 err = 1.0049e-03 -[ Info: VUMPS conv 1: obj = -4.380043660553e+00 err = 1.7391973790e-04 time = 0.74 sec -[ Info: VUMPS init: obj = -4.380043660559e+00 err = 9.1190e-04 -[ Info: VUMPS conv 1: obj = -4.380048624312e+00 err = 1.9138539579e-04 time = 1.45 sec -[ Info: VUMPS init: obj = -4.380048624328e+00 err = 8.3215e-04 -[ Info: VUMPS conv 1: obj = -4.380053188785e+00 err = 1.8374077529e-04 time = 1.52 sec -[ Info: VUMPS init: obj = -4.380053188823e+00 err = 6.8338e-04 -[ Info: VUMPS conv 1: obj = -4.380057133613e+00 err = 1.9170795863e-04 time = 2.81 sec -[ Info: VUMPS init: obj = -4.380057133620e+00 err = 6.0431e-04 -[ Info: VUMPS conv 1: obj = -4.380060541780e+00 err = 1.8351019737e-04 time = 2.56 sec -[ Info: VUMPS init: obj = -4.379609492466e+00 err = 4.0941e-03 -[ Info: VUMPS conv 19: obj = -4.379763054684e+00 err = 9.9541478983e-06 time = 14.21 sec -[ Info: CG: initializing with f = -4.379763054330, ‖∇f‖ = 3.1546e-05 -[ Info: CG: converged after 164 iterations and time 120.42 s: f = -4.379763208139, ‖∇f‖ = 9.5378e-07 +[ Info: VUMPS init: obj = -1.418138980906e+00 err = 5.4975e-01 +[ Info: VUMPS conv 7: obj = -4.376662745131e+00 err = 8.7994117461e-03 time = 5.77 sec +[ Info: VUMPS init: obj = -4.376662745131e+00 err = 2.1471e-02 +[ Info: VUMPS conv 8: obj = -4.378745586158e+00 err = 1.5392590844e-04 time = 0.67 sec +[ Info: VUMPS init: obj = -4.378745586158e+00 err = 8.1759e-03 +[ Info: VUMPS conv 6: obj = -4.379161067061e+00 err = 1.5555838326e-04 time = 0.44 sec +[ Info: VUMPS init: obj = -4.379161067061e+00 err = 6.1229e-03 +[ Info: VUMPS conv 5: obj = -4.379452165287e+00 err = 1.7041880980e-04 time = 0.61 sec +[ Info: VUMPS init: obj = -4.379452165287e+00 err = 5.6967e-03 +[ Info: VUMPS conv 4: obj = -4.379651685338e+00 err = 1.9086298845e-04 time = 0.47 sec +[ Info: VUMPS init: obj = -4.379651685338e+00 err = 4.1042e-03 +[ Info: VUMPS conv 4: obj = -4.379735588923e+00 err = 1.4107256275e-04 time = 0.59 sec +[ Info: VUMPS init: obj = -4.379735588923e+00 err = 3.5761e-03 +[ Info: VUMPS conv 3: obj = -4.379797862274e+00 err = 1.3847431310e-04 time = 0.47 sec +[ Info: VUMPS init: obj = -4.379797862274e+00 err = 2.7695e-03 +[ Info: VUMPS conv 2: obj = -4.379838498353e+00 err = 1.8068332663e-04 time = 0.33 sec +[ Info: VUMPS init: obj = -4.379838498352e+00 err = 2.7281e-03 +[ Info: VUMPS conv 3: obj = -4.379878843329e+00 err = 1.9308669443e-04 time = 0.73 sec +[ Info: VUMPS init: obj = -4.379878843329e+00 err = 2.6902e-03 +[ Info: VUMPS conv 3: obj = -4.379929213991e+00 err = 1.7246879281e-04 time = 0.97 sec +[ Info: VUMPS init: obj = -4.379929213991e+00 err = 2.5537e-03 +[ Info: VUMPS conv 3: obj = -4.379968050470e+00 err = 1.7633765416e-04 time = 1.23 sec +[ Info: VUMPS init: obj = -4.379968050470e+00 err = 1.7676e-03 +[ Info: VUMPS conv 2: obj = -4.379986897845e+00 err = 1.7982149771e-04 time = 0.93 sec +[ Info: VUMPS init: obj = -4.379986897845e+00 err = 1.5832e-03 +[ Info: VUMPS conv 2: obj = -4.380001017729e+00 err = 1.8565198322e-04 time = 1.22 sec +[ Info: VUMPS init: obj = -4.380001017729e+00 err = 1.5102e-03 +[ Info: VUMPS conv 2: obj = -4.380013168324e+00 err = 1.5019893736e-04 time = 1.55 sec +[ Info: VUMPS init: obj = -4.380013168324e+00 err = 1.4228e-03 +[ Info: VUMPS conv 2: obj = -4.380024392465e+00 err = 1.7662069399e-04 time = 2.36 sec +[ Info: VUMPS init: obj = -4.380024392465e+00 err = 1.3320e-03 +[ Info: VUMPS conv 2: obj = -4.380038151479e+00 err = 1.5891915641e-04 time = 2.25 sec +[ Info: VUMPS init: obj = -4.380038151479e+00 err = 1.0034e-03 +[ Info: VUMPS conv 1: obj = -4.380043679645e+00 err = 1.6778281273e-04 time = 0.83 sec +[ Info: VUMPS init: obj = -4.380043679649e+00 err = 9.1040e-04 +[ Info: VUMPS conv 1: obj = -4.380048640213e+00 err = 1.8625092812e-04 time = 1.14 sec +[ Info: VUMPS init: obj = -4.380048640225e+00 err = 8.3097e-04 +[ Info: VUMPS conv 1: obj = -4.380053201812e+00 err = 1.8243182942e-04 time = 1.01 sec +[ Info: VUMPS init: obj = -4.380053201846e+00 err = 6.8228e-04 +[ Info: VUMPS conv 1: obj = -4.380057143072e+00 err = 1.9139721436e-04 time = 1.85 sec +[ Info: VUMPS init: obj = -4.380057143076e+00 err = 6.0387e-04 +[ Info: VUMPS conv 1: obj = -4.380060547201e+00 err = 1.8483235697e-04 time = 1.31 sec +[ Info: VUMPS init: obj = -4.379609329367e+00 err = 4.0967e-03 +[ Info: VUMPS conv 19: obj = -4.379762979233e+00 err = 9.9502287579e-06 time = 10.30 sec +[ Info: CG: initializing with f = -4.379762978878, ‖∇f‖ = 3.1550e-05 +[ Info: CG: converged after 166 iterations and time 130.43 s: f = -4.379763024867, ‖∇f‖ = 9.5155e-07 ┌ Info: Groundstate energy: -│ * numerical: -2.1899960609691504 +│ * numerical: -2.1899960610082476 └ * analytic: -2.190038374277775 ```` @@ -192,56 +189,56 @@ Groundstate energy: ```` ```` -[ Info: VUMPS init: obj = +2.937065302538e-01 err = 8.7934e-01 -[ Info: VUMPS conv 1: obj = -4.000000000000e+00 err = 1.5900820202e-15 time = 3.70 sec +[ Info: VUMPS init: obj = -3.310849814614e-03 err = 9.6889e-01 +[ Info: VUMPS conv 1: obj = -4.000000000000e+00 err = 1.7910070597e-15 time = 4.40 sec [ Info: VUMPS init: obj = -4.000000000000e+00 err = 3.3634e-01 [ Info: VUMPS conv 4: obj = -4.289650419703e+00 err = 1.8515305877e-04 time = 0.04 sec [ Info: VUMPS init: obj = -4.289650419703e+00 err = 1.1203e-01 -[ Info: VUMPS conv 6: obj = -4.359777339623e+00 err = 1.0862346095e-04 time = 0.63 sec -[ Info: VUMPS init: obj = -4.359777339623e+00 err = 4.4013e-02 -[ Info: VUMPS conv 6: obj = -4.372880960838e+00 err = 1.2171760223e-04 time = 1.72 sec -[ Info: VUMPS init: obj = -4.372880960838e+00 err = 3.2691e-02 -[ Info: VUMPS conv 4: obj = -4.375236954439e+00 err = 1.1813015096e-04 time = 0.25 sec -[ Info: VUMPS init: obj = -4.375236954439e+00 err = 2.9487e-02 -[ Info: VUMPS conv 7: obj = -4.378159083962e+00 err = 1.1902464282e-04 time = 0.51 sec -[ Info: VUMPS init: obj = -4.378159083962e+00 err = 1.9312e-02 -[ Info: VUMPS conv 5: obj = -4.379272964964e+00 err = 1.5792771069e-04 time = 0.43 sec -[ Info: VUMPS init: obj = -4.379272964963e+00 err = 9.9127e-03 -[ Info: VUMPS conv 4: obj = -4.379592225772e+00 err = 1.5571650676e-04 time = 0.43 sec -[ Info: VUMPS init: obj = -4.379592225772e+00 err = 6.4839e-03 -[ Info: VUMPS conv 4: obj = -4.379819373751e+00 err = 1.7517281564e-04 time = 0.58 sec -[ Info: VUMPS init: obj = -4.379819373751e+00 err = 3.8752e-03 -┌ Warning: VUMPS cancel 10: obj = -4.379964036343e+00 err = 2.1162274451e-04 time = 2.01 sec -└ @ MPSKit ~/Projects/MPSKit.jl/src/algorithms/groundstate/vumps.jl:73 -[ Info: VUMPS init: obj = -4.379964036342e+00 err = 2.8980e-03 -[ Info: VUMPS conv 3: obj = -4.380010370438e+00 err = 1.4802687624e-04 time = 0.98 sec -[ Info: VUMPS init: obj = -4.380010370438e+00 err = 2.0599e-03 -[ Info: VUMPS conv 3: obj = -4.380041745658e+00 err = 1.6343274067e-04 time = 1.23 sec -[ Info: VUMPS init: obj = -4.380041745658e+00 err = 1.2361e-03 -[ Info: VUMPS conv 2: obj = -4.380055779868e+00 err = 1.8388289197e-04 time = 0.89 sec -[ Info: VUMPS init: obj = -4.380055779868e+00 err = 8.5853e-04 -[ Info: VUMPS conv 2: obj = -4.380064751094e+00 err = 1.3920397153e-04 time = 1.08 sec -[ Info: VUMPS init: obj = -4.380064751094e+00 err = 5.2218e-04 -[ Info: VUMPS conv 1: obj = -4.380067968749e+00 err = 1.5790331564e-04 time = 0.52 sec -[ Info: VUMPS init: obj = -4.380067968843e+00 err = 2.8321e-04 -[ Info: VUMPS conv 1: obj = -4.380069973425e+00 err = 1.1574874270e-04 time = 0.71 sec -[ Info: VUMPS init: obj = -4.380069973478e+00 err = 2.6036e-04 -[ Info: VUMPS conv 1: obj = -4.380071401187e+00 err = 1.0406988217e-04 time = 0.78 sec -[ Info: VUMPS init: obj = -4.380071401378e+00 err = 1.8008e-04 -[ Info: VUMPS conv 1: obj = -4.380072413626e+00 err = 9.4811643803e-05 time = 0.74 sec -[ Info: VUMPS init: obj = -4.380072413986e+00 err = 1.6570e-04 -[ Info: VUMPS conv 1: obj = -4.380073148781e+00 err = 8.3106483694e-05 time = 0.84 sec -[ Info: VUMPS init: obj = -4.380073149081e+00 err = 1.8963e-04 -[ Info: VUMPS conv 1: obj = -4.380073703995e+00 err = 7.4441883673e-05 time = 0.96 sec -[ Info: VUMPS init: obj = -4.380073704254e+00 err = 1.7788e-04 -[ Info: VUMPS conv 1: obj = -4.380074139011e+00 err = 6.9099006045e-05 time = 0.95 sec -[ Info: VUMPS init: obj = -4.379297063145e+00 err = 7.6292e-03 -┌ Warning: VUMPS cancel 100: obj = -4.379633680570e+00 err = 1.6837508787e-05 time = 24.07 sec -└ @ MPSKit ~/Projects/MPSKit.jl/src/algorithms/groundstate/vumps.jl:73 -[ Info: CG: initializing with f = -4.379633680570, ‖∇f‖ = 6.3970e-05 -[ Info: CG: converged after 20 iterations and time 10.12 s: f = -4.379633682318, ‖∇f‖ = 9.3330e-07 +[ Info: VUMPS conv 6: obj = -4.359865567633e+00 err = 1.0046406773e-04 time = 0.10 sec +[ Info: VUMPS init: obj = -4.359865567633e+00 err = 4.3643e-02 +[ Info: VUMPS conv 6: obj = -4.372880928577e+00 err = 1.3023373974e-04 time = 1.98 sec +[ Info: VUMPS init: obj = -4.372880928577e+00 err = 3.2693e-02 +[ Info: VUMPS conv 4: obj = -4.375236954734e+00 err = 1.1803697903e-04 time = 0.11 sec +[ Info: VUMPS init: obj = -4.375236954734e+00 err = 2.9487e-02 +[ Info: VUMPS conv 7: obj = -4.378159083989e+00 err = 1.1902548728e-04 time = 0.28 sec +[ Info: VUMPS init: obj = -4.378159083989e+00 err = 1.9312e-02 +[ Info: VUMPS conv 5: obj = -4.379272965038e+00 err = 1.5792669082e-04 time = 0.26 sec +[ Info: VUMPS init: obj = -4.379272965037e+00 err = 9.9127e-03 +[ Info: VUMPS conv 4: obj = -4.379592225819e+00 err = 1.5571559712e-04 time = 0.64 sec +[ Info: VUMPS init: obj = -4.379592225819e+00 err = 6.4839e-03 +[ Info: VUMPS conv 4: obj = -4.379819373805e+00 err = 1.7517067710e-04 time = 0.29 sec +[ Info: VUMPS init: obj = -4.379819373805e+00 err = 3.8752e-03 +┌ Warning: VUMPS cancel 10: obj = -4.379964036336e+00 err = 2.1162213354e-04 time = 1.09 sec +└ @ MPSKit ~/git/MPSKit.jl/src/algorithms/groundstate/vumps.jl:73 +[ Info: VUMPS init: obj = -4.379964036335e+00 err = 2.8980e-03 +[ Info: VUMPS conv 3: obj = -4.380010370460e+00 err = 1.4802655346e-04 time = 0.43 sec +[ Info: VUMPS init: obj = -4.380010370460e+00 err = 2.0599e-03 +[ Info: VUMPS conv 3: obj = -4.380041745689e+00 err = 1.6343265011e-04 time = 0.61 sec +[ Info: VUMPS init: obj = -4.380041745689e+00 err = 1.2361e-03 +[ Info: VUMPS conv 2: obj = -4.380055779898e+00 err = 1.8387821190e-04 time = 0.54 sec +[ Info: VUMPS init: obj = -4.380055779898e+00 err = 8.5853e-04 +[ Info: VUMPS conv 2: obj = -4.380064751120e+00 err = 1.3920368100e-04 time = 0.82 sec +[ Info: VUMPS init: obj = -4.380064751120e+00 err = 5.2503e-04 +[ Info: VUMPS conv 1: obj = -4.380067975496e+00 err = 1.5780382243e-04 time = 0.45 sec +[ Info: VUMPS init: obj = -4.380067975589e+00 err = 3.3360e-04 +[ Info: VUMPS conv 1: obj = -4.380070355793e+00 err = 1.3292468609e-04 time = 0.61 sec +[ Info: VUMPS init: obj = -4.380070355905e+00 err = 2.0509e-04 +[ Info: VUMPS conv 1: obj = -4.380072130176e+00 err = 1.1399026715e-04 time = 0.89 sec +[ Info: VUMPS init: obj = -4.380072130335e+00 err = 1.3827e-04 +[ Info: VUMPS conv 1: obj = -4.380073468947e+00 err = 8.7912179069e-05 time = 1.30 sec +[ Info: VUMPS init: obj = -4.380073469160e+00 err = 9.9636e-05 +[ Info: VUMPS conv 1: obj = -4.380074456012e+00 err = 6.8484929133e-05 time = 2.04 sec +[ Info: VUMPS init: obj = -4.380074456584e+00 err = 7.7488e-05 +[ Info: VUMPS conv 1: obj = -4.380075158440e+00 err = 6.5951776335e-05 time = 4.26 sec +[ Info: VUMPS init: obj = -4.380075159046e+00 err = 6.3891e-05 +[ Info: VUMPS conv 1: obj = -4.380075656659e+00 err = 4.7588615250e-05 time = 6.33 sec +[ Info: VUMPS init: obj = -4.379308902466e+00 err = 8.0000e-03 +┌ Warning: VUMPS cancel 100: obj = -4.379693243679e+00 err = 1.5955101981e-05 time = 17.79 sec +└ @ MPSKit ~/git/MPSKit.jl/src/algorithms/groundstate/vumps.jl:73 +[ Info: CG: initializing with f = -4.379693243679, ‖∇f‖ = 5.7833e-05 +[ Info: CG: converged after 13 iterations and time 4.96 s: f = -4.379693244608, ‖∇f‖ = 6.3799e-07 ┌ Info: Groundstate energy: -│ * numerical: -2.1900011720804304 +│ * numerical: -2.1900153475149438 └ * analytic: -2.190038374277775 ```` @@ -277,28 +274,28 @@ E_holon, ϕ_holon = excitations(H_u1_su2, alg, momenta, ```` [ Info: Found excitations for momentum = -3.141592653589793 +[ Info: Found excitations for momentum = -2.945243112740431 [ Info: Found excitations for momentum = -2.748893571891069 [ Info: Found excitations for momentum = -2.552544031041707 -[ Info: Found excitations for momentum = -2.945243112740431 [ Info: Found excitations for momentum = -2.356194490192345 [ Info: Found excitations for momentum = -2.1598449493429825 -[ Info: Found excitations for momentum = -1.7671458676442586 [ Info: Found excitations for momentum = -1.9634954084936207 +[ Info: Found excitations for momentum = -1.7671458676442586 [ Info: Found excitations for momentum = -1.5707963267948966 [ Info: Found excitations for momentum = -1.3744467859455345 -[ Info: Found excitations for momentum = -0.9817477042468103 [ Info: Found excitations for momentum = -1.1780972450961724 +[ Info: Found excitations for momentum = -0.9817477042468103 [ Info: Found excitations for momentum = -0.7853981633974483 [ Info: Found excitations for momentum = -0.5890486225480862 -[ Info: Found excitations for momentum = -0.19634954084936207 [ Info: Found excitations for momentum = -0.39269908169872414 +[ Info: Found excitations for momentum = -0.19634954084936207 [ Info: Found excitations for momentum = 0.0 [ Info: Found excitations for momentum = 0.19634954084936207 [ Info: Found excitations for momentum = 0.39269908169872414 [ Info: Found excitations for momentum = 0.5890486225480862 [ Info: Found excitations for momentum = 0.7853981633974483 -[ Info: Found excitations for momentum = 1.1780972450961724 [ Info: Found excitations for momentum = 0.9817477042468103 +[ Info: Found excitations for momentum = 1.1780972450961724 [ Info: Found excitations for momentum = 1.3744467859455345 [ Info: Found excitations for momentum = 1.5707963267948966 [ Info: Found excitations for momentum = 1.7671458676442586 @@ -306,19 +303,19 @@ E_holon, ϕ_holon = excitations(H_u1_su2, alg, momenta, [ Info: Found excitations for momentum = 2.1598449493429825 [ Info: Found excitations for momentum = 2.356194490192345 [ Info: Found excitations for momentum = 2.552544031041707 -[ Info: Found excitations for momentum = 2.945243112740431 [ Info: Found excitations for momentum = 2.748893571891069 +[ Info: Found excitations for momentum = 2.945243112740431 [ Info: Found excitations for momentum = 3.141592653589793 -[ Info: Found excitations for momentum = -2.748893571891069 [ Info: Found excitations for momentum = -3.141592653589793 -[ Info: Found excitations for momentum = -2.552544031041707 [ Info: Found excitations for momentum = -2.945243112740431 -[ Info: Found excitations for momentum = -1.9634954084936207 +[ Info: Found excitations for momentum = -2.748893571891069 +[ Info: Found excitations for momentum = -2.552544031041707 [ Info: Found excitations for momentum = -2.356194490192345 -[ Info: Found excitations for momentum = -1.7671458676442586 [ Info: Found excitations for momentum = -2.1598449493429825 -[ Info: Found excitations for momentum = -1.3744467859455345 +[ Info: Found excitations for momentum = -1.9634954084936207 +[ Info: Found excitations for momentum = -1.7671458676442586 [ Info: Found excitations for momentum = -1.5707963267948966 +[ Info: Found excitations for momentum = -1.3744467859455345 [ Info: Found excitations for momentum = -1.1780972450961724 [ Info: Found excitations for momentum = -0.9817477042468103 [ Info: Found excitations for momentum = -0.7853981633974483 @@ -330,8 +327,8 @@ E_holon, ϕ_holon = excitations(H_u1_su2, alg, momenta, [ Info: Found excitations for momentum = 0.39269908169872414 [ Info: Found excitations for momentum = 0.5890486225480862 [ Info: Found excitations for momentum = 0.7853981633974483 -[ Info: Found excitations for momentum = 1.1780972450961724 [ Info: Found excitations for momentum = 0.9817477042468103 +[ Info: Found excitations for momentum = 1.1780972450961724 [ Info: Found excitations for momentum = 1.3744467859455345 [ Info: Found excitations for momentum = 1.5707963267948966 [ Info: Found excitations for momentum = 1.7671458676442586 @@ -398,123 +395,123 @@ end - + - + - + - + - - + + - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + ``` @@ -540,123 +537,121 @@ end - + - + - + - + - - + + - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + ``` @@ -721,137 +716,135 @@ end - + - + - + - + - - + + - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + - - - - - - - - - + + + + + + + + + - - - + + + ``` diff --git a/docs/src/examples/quantum1d/6.hubbard/main.ipynb b/docs/src/examples/quantum1d/6.hubbard/main.ipynb index 5c0b7885b..e28d0fa35 100644 --- a/docs/src/examples/quantum1d/6.hubbard/main.ipynb +++ b/docs/src/examples/quantum1d/6.hubbard/main.ipynb @@ -369,11 +369,11 @@ "file_extension": ".jl", "mimetype": "application/julia", "name": "julia", - "version": "1.11.4" + "version": "1.10.4" }, "kernelspec": { - "name": "julia-1.11", - "display_name": "Julia 1.11.4", + "name": "julia-1.10", + "display_name": "Julia 1.10.4", "language": "julia" } }, diff --git a/docs/src/examples/quantum1d/7.xy-finiteT/index.md b/docs/src/examples/quantum1d/7.xy-finiteT/index.md index f080e52b9..5cc714959 100644 --- a/docs/src/examples/quantum1d/7.xy-finiteT/index.md +++ b/docs/src/examples/quantum1d/7.xy-finiteT/index.md @@ -124,11 +124,11 @@ println("Exact (N=Inf):\t", groundstate_energy(J, Inf)) ```` ```` -[ Info: DMRG2 1: obj = -5.004084850344e+00 err = 9.8257238519e-01 time = 1.50 min -[ Info: DMRG2 2: obj = -5.004096941404e+00 err = 1.1498230845e-06 time = 0.46 sec -[ Info: DMRG2 3: obj = -5.004096975044e+00 err = 1.9792426587e-09 time = 0.61 sec -[ Info: DMRG2 conv 4: obj = -5.004096975044e+00 err = 1.1590728377e-13 time = 1.52 min -Numerical: -0.15637803047010948 +[ Info: DMRG2 1: obj = -5.004084861350e+00 err = 9.8293025834e-01 time = 1.27 min +[ Info: DMRG2 2: obj = -5.004096939909e+00 err = 1.1403678304e-06 time = 0.55 sec +[ Info: DMRG2 3: obj = -5.004096975044e+00 err = 2.3606914201e-09 time = 0.91 sec +[ Info: DMRG2 conv 4: obj = -5.004096975044e+00 err = 1.1624035068e-13 time = 1.31 min +Numerical: -0.15637803047010954 Exact (N=32): -0.15637803047254015 Exact (N=Inf): -0.15915494309189535 @@ -237,89 +237,89 @@ end - + - + - + - + - - + + - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + - - + + - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + ``` @@ -372,87 +372,87 @@ end - + - + - + - + - - + + - - - - - - - - - - - - - - - - - - - - - - - - - - - - - + + + + + + + + + + + + + + + + + + + + + + + + + + + + + - - + + - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + ``` @@ -506,89 +506,89 @@ end - + - + - + - + - - + + - - - - - - - - - - - - - - - - - - - - - - - - - - - - - + + + + + + + + + + + + + + + + + + + + + + + + + + + + + - - + + - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + ``` @@ -660,73 +660,73 @@ end - + - + - + - + - - + + - - - - - - - - - - - - - - - - - - - - - - - - - - - + + + + + + + + + + + + + + + + + + + + + + + + + + + - - + + - - - - - - - - - - - - - - - - - - - + + + + + + + + + + + + + + + + + + + ``` @@ -806,76 +806,76 @@ end - + - + - + - + - - + + - - - - - - - - - - - - - - - - - - - - - - - - - - - - - + + + + + + + + + + + + + + + + + + + + + + + + + + + + + - - + + - - - - - - - - - - - - - - - - - - - - + + + + + + + + + + + + + + + + + + + + ``` @@ -944,79 +944,79 @@ end - + - + - + - + - - + + - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + - - + + - - - - - - - - - - - - - - - - - - - - - + + + + + + + + + + + + + + + + + + + + + ``` diff --git a/docs/src/examples/quantum1d/7.xy-finiteT/main.ipynb b/docs/src/examples/quantum1d/7.xy-finiteT/main.ipynb index a19eabae6..7a14ede2e 100644 --- a/docs/src/examples/quantum1d/7.xy-finiteT/main.ipynb +++ b/docs/src/examples/quantum1d/7.xy-finiteT/main.ipynb @@ -653,11 +653,11 @@ "file_extension": ".jl", "mimetype": "application/julia", "name": "julia", - "version": "1.11.4" + "version": "1.10.4" }, "kernelspec": { - "name": "julia-1.11", - "display_name": "Julia 1.11.4", + "name": "julia-1.10", + "display_name": "Julia 1.10.4", "language": "julia" } }, diff --git a/examples/Cache.toml b/examples/Cache.toml index 76f8dad0c..faa708cac 100644 --- a/examples/Cache.toml +++ b/examples/Cache.toml @@ -1,11 +1,11 @@ [classic2d] -"1.hard-hexagon" = "b380d99250b1908c03248a569ccb2729891f94ccad2a840dc9866af63ce98d2e" +"1.hard-hexagon" = "6a351c0f8ee413bb8e3ab4ee09de9493856bd1c308c060b6a7e052c5d38b7bc9" [quantum1d] -"2.haldane" = "4f3e177a30df907cf7b9a47958ddf97741e4d0398ae1876e0f4977ccb93ae651" -"6.hubbard" = "a0f80e68edf8b12b90b1a7c8fe8e87d14be2f3dbea4fa0491ef53cf5eea0dbd7" -"7.xy-finiteT" = "0243b44634d38d04ddd96b180fd1fc24e20485e97135cb7c805c9d60baab9ccb" -"3.ising-dqpt" = "f437ac4560db51f67e3997867c02bb8844b1077731b1dd1a06c61ab7c4846e97" -"5.haldane-spt" = "6443696b20d416652951c646cf78b604a7fd0d79804995588491ea53ca210d57" -"4.xxz-heisenberg" = "f8359ad29a028b01d41913014c0a442b99c3bd5cd403c23bd2a59f646a7c78cf" -"1.ising-cft" = "bd1253451973aa73c4fee6b64d81ee52cc2feba5b9632f87ce332eea3a05f460" +"2.haldane" = "769a4ef9d7a878d011f2e5bbcc790f2db23c63fd673db9a4729a8501392c1fdf" +"6.hubbard" = "708c865dd9681d36d5d1b3dc43914d3df59d3f97ae7bcc6e3a24092890b8d44f" +"7.xy-finiteT" = "8ba5bd93407263202d3f07aa5e36d8387308a9107202ea06df9755ee128058c1" +"3.ising-dqpt" = "bfa10659a7c002fc7237948367474064f81071d773d493f30f222c0842fc416f" +"5.haldane-spt" = "1f4a054ff1528f1f273fbeec6258a7bad21fb789bfffb073dbb01e697d1b275c" +"4.xxz-heisenberg" = "8d27d51408b6dbb5550afd7d2d80307e134f8a32a163923e12e74ad3ce8c5041" +"1.ising-cft" = "b0b2106eb3b6ba5b4c505adfcae7a4ed49a788e4545df08d053416feccec1797" diff --git a/examples/make.jl b/examples/make.jl index be471a00c..6cdab5780 100644 --- a/examples/make.jl +++ b/examples/make.jl @@ -36,7 +36,15 @@ function setcached(root, name) return open(f -> TOML.print(f, cache), CACHEFILE, "w") end -checksum(root, name) = bytes2hex(sha256(joinpath(@__DIR__, root, name, "main.jl"))) +# generate checksum based on path relative to ~/.../MPSKit.jl +# such that different users do not have to rerun already cached examples +function checksum(root, name) + example_path = joinpath(@__DIR__, root, name, "main.jl") + @assert isfile(example_path) + return open(example_path, "r") do io + return bytes2hex(sha256(io)) + end +end # ---------------------------------------------------------------------------------------- # # Building diff --git a/examples/quantum1d/1.ising-cft/main.jl b/examples/quantum1d/1.ising-cft/main.jl index 27c9594f0..8bb4cd781 100644 --- a/examples/quantum1d/1.ising-cft/main.jl +++ b/examples/quantum1d/1.ising-cft/main.jl @@ -53,8 +53,8 @@ or in the code as: """ function O_shift(L) - τ = BraidingTensor{ComplexF64}(ℂ^2, ℂ^2) - O = TensorMap(τ) + I = id(ComplexF64, ℂ^2) + @tensor O[W S; N E] := I[W; N] * I[S; E] return periodic_boundary_conditions(InfiniteMPO([O]), L) end @@ -64,7 +64,8 @@ operator. However, there is a subtlety because of the degeneracies in the energy eigenvalues. The eigensolver will find an orthonormal basis within each energy subspace, but this basis is not necessarily a basis of eigenstates of the translation operator. In order to fix this, we diagonalize the translation operator within each energy subspace. -The resulting energy levels have one-to-one correspondence to the operators in CFT, where the momentum is related to their conformal spin as $P_n = \frac{2\pi}{L}S_n$. +The resulting energy levels have one-to-one correspondence to the operators in CFT, where +the momentum is related to their conformal spin as $P_n = \frac{2\pi}{L}S_n$. """ function fix_degeneracies(basis) @@ -90,12 +91,21 @@ append!(momenta, fix_degeneracies(states[12:12])) append!(momenta, fix_degeneracies(states[13:16])) append!(momenta, fix_degeneracies(states[17:18])) -# Calculating scaling dimensions from the energy gap +md""" +We can compute the scaling dimensions $\Delta_n$ of the operators in the CFT from the +energy gap of the corresponding excitations as $E_n - E_0 = \frac{2\pi v}{L} \Delta_n$, +where $v = 2$. If we plot these scaling dimensions against the conformal spin $S_n$ from +above, we retrieve the familiar spectrum of the Ising CFT. +""" + v = 2.0 Δ = real.(energies[1:18] .- energies[1]) ./ (2π * v / L) -p = plot(momenta, real.(Δ); - seriestype=:scatter, xlabel="momentum", ylabel="energy", legend=false) -vline!(p, [2π / L * i for i in -3:3]; color="gray", linestyle=:dash) +S = momenta ./ (2π / L) + +p = plot(S, real.(Δ); + seriestype=:scatter, xlabel="conformal spin (S)", ylabel="scaling dimension (Δ)", + legend=false) +vline!(p, -3:3; color="gray", linestyle=:dash) hline!(p, [0, 1 / 8, 1, 9 / 8, 2, 17 / 8]; color="gray", linestyle=:dash) p @@ -117,24 +127,28 @@ ansatz. This returns quasiparticle states, which can be converted to regular `Fi objects. """ -E_ex, qps = excitations(H_mps, QuasiparticleAnsatz(), ψ, envs; num=16) +E_ex, qps = excitations(H_mps, QuasiparticleAnsatz(), ψ, envs; num=18) states_mps = vcat(ψ, map(qp -> convert(FiniteMPS, qp), qps)) -E_mps = map(x -> expectation_value(x, H_mps), states_mps) +energies_mps = map(x -> expectation_value(x, H_mps), states_mps) momenta_mps = Float64[] -append!(momenta_mps, fix_degeneracies(states[1:1])) -append!(momenta_mps, fix_degeneracies(states[2:2])) -append!(momenta_mps, fix_degeneracies(states[3:3])) -append!(momenta_mps, fix_degeneracies(states[4:5])) -append!(momenta_mps, fix_degeneracies(states[6:9])) -append!(momenta_mps, fix_degeneracies(states[10:11])) -append!(momenta_mps, fix_degeneracies(states[12:12])) -append!(momenta_mps, fix_degeneracies(states[13:16])) +append!(momenta_mps, fix_degeneracies(states_mps[1:1])) +append!(momenta_mps, fix_degeneracies(states_mps[2:2])) +append!(momenta_mps, fix_degeneracies(states_mps[3:3])) +append!(momenta_mps, fix_degeneracies(states_mps[4:5])) +append!(momenta_mps, fix_degeneracies(states_mps[6:9])) +append!(momenta_mps, fix_degeneracies(states_mps[10:11])) +append!(momenta_mps, fix_degeneracies(states_mps[12:12])) +append!(momenta_mps, fix_degeneracies(states_mps[13:16])) +append!(momenta_mps, fix_degeneracies(states_mps[17:18])) v = 2.0 -Δ = real.(energies[1:18] .- energies[1]) ./ (2π * v / L) -plot(momenta_mps, Δ; - seriestype=:scatter, xlabel="momentum", ylabel="energy", legend=false) -vline!(p, [2π / L * i for i in -3:3]; color="gray", linestyle=:dash) -hline!(p, [0, 1 / 8, 1, 9 / 8, 2, 17 / 8]; color="gray", linestyle=:dash) -p +Δ_mps = real.(energies_mps[1:18] .- energies_mps[1]) ./ (2π * v / L_mps) +S_mps = momenta_mps ./ (2π / L_mps) + +p_mps = plot(S_mps, real.(Δ_mps); + seriestype=:scatter, xlabel="conformal spin (S)", + ylabel="scaling dimension (Δ)", legend=false) +vline!(p_mps, -3:3; color="gray", linestyle=:dash) +hline!(p_mps, [0, 1 / 8, 1, 9 / 8, 2, 17 / 8]; color="gray", linestyle=:dash) +p_mps