diff --git a/CHANGELOG.md b/CHANGELOG.md deleted file mode 100644 index 26cc89842..000000000 --- a/CHANGELOG.md +++ /dev/null @@ -1,33 +0,0 @@ -# Changelog - -All notable changes to this project will be documented in this file. - -The format is based on [Keep a Changelog](https://keepachangelog.com/en/1.1.0/), -and this project adheres to [Semantic Versioning](https://semver.org/spec/v2.0.0.html). - -## [Unreleased] - -### Added - -- `LocalPreferences.toml` file to ensure `TensorOperations` properly precompiles on testing - infrastructure - -### Fixed - -- Dynamic tolerances yielded `NaN` during the initialization stage due to `1 / sqrt(iter)` - where `iter = 0`. -- `InfiniteMPOHamiltonian` environments with low bond dimension and high Krylov dimension now are properly - clamped. - -### Changed - -- The `changebonds(state, ::RandExpand)` algorithm now no longer has to perform a - truncated SVD to obtain the desired spaces, and instead sample the space directly - and then generates a random isometry. This should be slightly more performant, but - otherwise equivalent. - -### Deprecated - -### Removed - -[unreleased]: https://github.com/quantumkithub/pepskit.jl/compare/v0.13.8...HEAD diff --git a/CITATION.cff b/CITATION.cff index ee1d9f18f..bd968cb18 100644 --- a/CITATION.cff +++ b/CITATION.cff @@ -1,18 +1,18 @@ cff-version: 1.2.0 message: "If you use this software, please cite it as below." authors: -- family-names: "Van Damme" - given-names: "Maarten" - orcid: "https://orcid.org/0000-0002-1558-0568" - family-names: "Devos" given-names: "Lukas" orcid: "https://orcid.org/0000-0002-0256-4200" +- family-names: "Van Damme" + given-names: "Maarten" + orcid: "https://orcid.org/0000-0002-1558-0568" - family-names: "Haegeman" given-names: "Jutho" orcid: "https://orcid.org/0000-0002-0858-291X" title: "MPSKit" -version: 0.13.5 +version: 0.13.9 doi: 10.5281/zenodo.10654900 -date-released: 2025-09-05 +date-released: 2026-02-03 url: "https://github.com/QuantumKitHub/MPSKit.jl" diff --git a/Project.toml b/Project.toml index d2125adcc..cda6e601c 100644 --- a/Project.toml +++ b/Project.toml @@ -1,6 +1,6 @@ name = "MPSKit" uuid = "bb1c41ca-d63c-52ed-829e-0820dda26502" -version = "0.13.8" +version = "0.13.9" authors = "Lukas Devos, Maarten Van Damme and contributors" [deps] diff --git a/docs/make.jl b/docs/make.jl index 301cb6230..560e0f61b 100644 --- a/docs/make.jl +++ b/docs/make.jl @@ -70,6 +70,7 @@ makedocs(; "Examples" => "examples/index.md", "Library" => "lib/lib.md", "References" => "references.md", + "Changelog" => "changelog.md", ], checkdocs = :exports, doctest = true, diff --git a/docs/src/changelog.md b/docs/src/changelog.md new file mode 100644 index 000000000..44e916bd4 --- /dev/null +++ b/docs/src/changelog.md @@ -0,0 +1,79 @@ +# Changelog + +All notable changes to this project will be documented in this file. + +The format is based on [Keep a Changelog](https://keepachangelog.com/en/1.1.0/), +and this project adheres to [Semantic Versioning](https://semver.org/spec/v2.0.0.html). + +## Guidelines for updating this changelog + +When making changes to this project, please update the "Unreleased" section with your changes under the appropriate category: + +- **Added** for new features. +- **Changed** for changes in existing functionality. +- **Deprecated** for soon-to-be removed features. +- **Removed** for now removed features. +- **Fixed** for any bug fixes. + +When releasing a new version, move the "Unreleased" changes to a new version section with the release date. + +## [Unreleased](https://github.com/QuantumKitHub/MPSKit.jl/compare/v0.13.9...HEAD) + +### Added + +### Changed + +### Deprecated + +### Removed + +### Fixed + +## [0.13.9](https://github.com/QuantumKitHub/MPSKit.jl/compare/v0.13.8...v0.13.9) - 2026-02-03 + +### Added + +- `LocalPreferences.toml` file to ensure `TensorOperations` properly precompiles on testing + infrastructure +- `GeometryStyle` and `OperatorStyle` traits for dispatching on finite/infinite geometry and + operator types ([#352](https://github.com/QuantumKitHub/MPSKit.jl/pull/352), [#354](https://github.com/QuantumKitHub/MPSKit.jl/pull/354)) +- `Base.isfinite` methods for MPS types ([#347](https://github.com/QuantumKitHub/MPSKit.jl/pull/347)) +- Bose-Hubbard example ([#342](https://github.com/QuantumKitHub/MPSKit.jl/pull/342)) +- WindowMPS example update ([#350](https://github.com/QuantumKitHub/MPSKit.jl/pull/350)) +- Multifusion category compatibility ([#297](https://github.com/QuantumKitHub/MPSKit.jl/pull/297)) + +### Fixed + +- Dynamic tolerances yielded `NaN` during the initialization stage due to `1 / sqrt(iter)` + where `iter = 0` ([#335](https://github.com/QuantumKitHub/MPSKit.jl/pull/335)) +- `InfiniteMPOHamiltonian` environments with low bond dimension and high Krylov dimension now are properly + clamped ([#335](https://github.com/QuantumKitHub/MPSKit.jl/pull/335)) +- Logical operator precedence in `getproperty` function ([#346](https://github.com/QuantumKitHub/MPSKit.jl/pull/346)) +- Typo in `VUMPSSvdCut` ([#361](https://github.com/QuantumKitHub/MPSKit.jl/pull/361)) +- Typo in time formatting for logs ([#336](https://github.com/QuantumKitHub/MPSKit.jl/pull/336)) +- Domain/codomain of `MPODerivativeOperator` ([#370](https://github.com/QuantumKitHub/MPSKit.jl/pull/370)) +- In-place operations handled more carefully ([#337](https://github.com/QuantumKitHub/MPSKit.jl/pull/337)) +- Orthogonalization algorithms now use correct methods ([#373](https://github.com/QuantumKitHub/MPSKit.jl/pull/373)) + +### Changed + +- The `changebonds(state, ::RandExpand)` algorithm now no longer has to perform a + truncated SVD to obtain the desired spaces, and instead sample the space directly + and then generates a random isometry. This should be slightly more performant, but + otherwise equivalent ([#335](https://github.com/QuantumKitHub/MPSKit.jl/pull/335)) +- `IDMRG` refactored to follow the `IterativeSolver` interface and share code between + `IDMRG` and `IDMRG2` ([#348](https://github.com/QuantumKitHub/MPSKit.jl/pull/348)) +- Bumped compatibility for TensorKit 0.16 and MatrixAlgebraKit 0.6 ([#365](https://github.com/QuantumKitHub/MPSKit.jl/pull/365)) +- Removed `_left_orth` and `_right_orth` workarounds in favor of new orthogonalization methods +- Reduced allocation while computing Galerkin error ([#366](https://github.com/QuantumKitHub/MPSKit.jl/pull/366)) +- Updated `show` methods to reflect new TensorKit printing ([#341](https://github.com/QuantumKitHub/MPSKit.jl/pull/341)) +- More informative errors for finite MPS ([#367](https://github.com/QuantumKitHub/MPSKit.jl/pull/367)) +- Minor documentation and docstring improvements ([#363](https://github.com/QuantumKitHub/MPSKit.jl/pull/363), [#372](https://github.com/QuantumKitHub/MPSKit.jl/pull/372), [#371](https://github.com/QuantumKitHub/MPSKit.jl/pull/371)) + +### Deprecated + +### Removed + +## [0.13.8](https://github.com/QuantumKitHub/MPSKit.jl/releases/tag/v0.13.8) - 2024-10-31 + +See full history and previous releases on [GitHub](https://github.com/QuantumKitHub/MPSKit.jl/releases). diff --git a/docs/src/examples/classic2d/1.hard-hexagon/index.md b/docs/src/examples/classic2d/1.hard-hexagon/index.md index c1dafc781..c9575c2e4 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.8839037051703853 S = 1.2807829622000708 ξ = 13.849682582620911 +F = 0.8839037051703852 S = 1.280782962156832 ξ = 13.849682579815157 ```` @@ -100,7 +100,7 @@ c = f.coeffs[2] ```` ```` -0.802524583072294 +0.802523415438401 ```` ````julia @@ -110,7 +110,7 @@ plot!(p, ξ -> f(ξ) / 6; label = "fit") ```` ```@raw html - + ``` --- diff --git a/docs/src/examples/classic2d/1.hard-hexagon/main.ipynb b/docs/src/examples/classic2d/1.hard-hexagon/main.ipynb index 81e443120..4eebe0d8d 100644 --- a/docs/src/examples/classic2d/1.hard-hexagon/main.ipynb +++ b/docs/src/examples/classic2d/1.hard-hexagon/main.ipynb @@ -147,11 +147,11 @@ "file_extension": ".jl", "mimetype": "application/julia", "name": "julia", - "version": "1.12.0" + "version": "1.12.4" }, "kernelspec": { "name": "julia-1.12", - "display_name": "Julia 1.12.0", + "display_name": "Julia 1.12.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 d69b6d38f..2005cddcc 100644 --- a/docs/src/examples/quantum1d/1.ising-cft/index.md +++ b/docs/src/examples/quantum1d/1.ising-cft/index.md @@ -26,19 +26,22 @@ 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))) +12-site FiniteMPOHamiltonian(ComplexF64, TensorKit.ComplexSpace) with maximal dimension 6: +┬─[12]─ ℂ^2 +│ (ℂ^1 ⊞ ℂ^1 ⊞ ℂ^1 ⊞ ⋯ ⊞ ℂ^1 ⊞ ℂ^1 ⊞ ℂ^1) +┼─[11]─ ℂ^2 +│ (ℂ^1 ⊞ ℂ^1 ⊞ ℂ^1 ⊞ ⋯ ⊞ ℂ^1 ⊞ ℂ^1 ⊞ ℂ^1) +┼─[10]─ ℂ^2 +│ (ℂ^1 ⊞ ℂ^1 ⊞ ℂ^1 ⊞ ⋯ ⊞ ℂ^1 ⊞ ℂ^1 ⊞ ℂ^1) +┼─[9]─ ℂ^2 +│ (ℂ^1 ⊞ ℂ^1 ⊞ ℂ^1 ⊞ ⋯ ⊞ ℂ^1 ⊞ ℂ^1 ⊞ ℂ^1) +│ ⋮ +│ (ℂ^1 ⊞ ℂ^1 ⊞ ℂ^1 ⊞ ⋯ ⊞ ℂ^1 ⊞ ℂ^1 ⊞ ℂ^1) +┼─[3]─ ℂ^2 +│ (ℂ^1 ⊞ ℂ^1 ⊞ ℂ^1 ⊞ ⋯ ⊞ ℂ^1 ⊞ ℂ^1 ⊞ ℂ^1) +┼─[2]─ ℂ^2 +│ (ℂ^1 ⊞ ℂ^1 ⊞ ℂ^1 ⊞ ⋯ ⊞ ℂ^1 ⊞ ℂ^1 ⊞ ℂ^1) +┴─[1]─ ℂ^2 ```` @@ -125,23 +128,23 @@ append!(momenta, fix_degeneracies(states[17:18])) ```` 18-element Vector{Float64}: - -5.378354829155499e-17 - -7.106988585426104e-20 - -1.1970611374275794e-17 - -0.5235987755982987 - 0.523598775598299 - -1.0471975511965979 - 1.0471975511965974 - 0.5235987755982988 - -0.5235987755982986 - 1.0471975511965976 - -1.0471975511965979 - 1.705690506440868e-17 + -3.369480642089346e-17 + 3.9966961823263104e-17 + -2.4201642794270328e-17 + 0.5235987755982987 + -0.5235987755982985 + 1.0471975511965979 + -1.047197551196598 + -0.5235987755982991 + 0.5235987755982989 + -1.0471975511965974 + 1.0471975511965979 + -1.4167067287245392e-18 -1.570796326794897 - 1.5707963267948963 - 1.047197551196598 - -1.0471975511965976 1.5707963267948968 + -1.0471975511965979 + 1.0471975511965979 + 1.5707963267948966 -1.5707963267948966 ```` @@ -182,46 +185,34 @@ D = 64 ```` ```` -[ Info: DMRG init: obj = -1.959798006201e+01 err = 3.9034e-01 -[ Info: DMRG 1: obj = -2.548860604005e+01 err = 1.7104122458e-02 time = 2.94 sec -[ Info: DMRG 2: obj = -2.549098968635e+01 err = 4.7769187563e-06 time = 1.19 sec -[ Info: DMRG 3: obj = -2.549098968636e+01 err = 1.4444395774e-07 time = 0.56 sec -[ Info: DMRG 4: obj = -2.549098968636e+01 err = 2.4683334253e-08 time = 0.47 sec -[ Info: DMRG 5: obj = -2.549098968636e+01 err = 9.2494811019e-09 time = 0.47 sec -[ Info: DMRG 6: obj = -2.549098968636e+01 err = 4.2700930707e-09 time = 0.47 sec -[ Info: DMRG 7: obj = -2.549098968636e+01 err = 2.8468477933e-09 time = 0.46 sec -[ Info: DMRG 8: obj = -2.549098968636e+01 err = 2.2214754548e-09 time = 0.48 sec -[ Info: DMRG 9: obj = -2.549098968636e+01 err = 1.7977223058e-09 time = 0.47 sec -[ Info: DMRG 10: obj = -2.549098968636e+01 err = 1.4634778351e-09 time = 0.47 sec -[ Info: DMRG 11: obj = -2.549098968636e+01 err = 1.3630887037e-09 time = 0.48 sec -[ Info: DMRG 12: obj = -2.549098968636e+01 err = 1.5594810847e-09 time = 0.47 sec -[ Info: DMRG 13: obj = -2.549098968636e+01 err = 1.7301607123e-09 time = 0.46 sec -[ Info: DMRG 14: obj = -2.549098968636e+01 err = 1.8601697500e-09 time = 0.46 sec -[ Info: DMRG 15: obj = -2.549098968636e+01 err = 1.9357644569e-09 time = 0.44 sec -[ Info: DMRG 16: obj = -2.549098968636e+01 err = 1.9483292750e-09 time = 1.29 sec -[ Info: DMRG 17: obj = -2.549098968636e+01 err = 1.9124252703e-09 time = 0.40 sec -[ Info: DMRG 18: obj = -2.549098968636e+01 err = 1.8333796427e-09 time = 0.36 sec -[ Info: DMRG 19: obj = -2.549098968636e+01 err = 1.7086486342e-09 time = 0.39 sec -[ Info: DMRG 20: obj = -2.549098968636e+01 err = 1.5546926753e-09 time = 0.40 sec -[ Info: DMRG 21: obj = -2.549098968636e+01 err = 1.3877708971e-09 time = 0.40 sec -[ Info: DMRG 22: obj = -2.549098968636e+01 err = 1.2209468902e-09 time = 0.40 sec -[ Info: DMRG 23: obj = -2.549098968636e+01 err = 1.0630324605e-09 time = 0.41 sec -[ Info: DMRG 24: obj = -2.549098968636e+01 err = 9.1895063136e-10 time = 0.41 sec -[ Info: DMRG 25: obj = -2.549098968636e+01 err = 7.9071631158e-10 time = 0.41 sec -[ Info: DMRG 26: obj = -2.549098968636e+01 err = 6.7846224132e-10 time = 0.41 sec -[ Info: DMRG 27: obj = -2.549098968636e+01 err = 5.8125445739e-10 time = 0.41 sec -[ Info: DMRG 28: obj = -2.549098968636e+01 err = 4.9765022739e-10 time = 0.41 sec -[ Info: DMRG 29: obj = -2.549098968636e+01 err = 4.2604131419e-10 time = 0.41 sec -[ Info: DMRG 30: obj = -2.549098968636e+01 err = 3.6484596128e-10 time = 0.41 sec -[ Info: DMRG 31: obj = -2.549098968636e+01 err = 3.1260439575e-10 time = 0.41 sec -[ Info: DMRG 32: obj = -2.549098968636e+01 err = 2.6801737353e-10 time = 0.41 sec -[ Info: DMRG 33: obj = -2.549098968636e+01 err = 2.2995321952e-10 time = 0.41 sec -[ Info: DMRG 34: obj = -2.549098968636e+01 err = 1.9743900235e-10 time = 0.41 sec -[ Info: DMRG 35: obj = -2.549098968636e+01 err = 1.6964470820e-10 time = 0.43 sec -[ Info: DMRG 36: obj = -2.549098968636e+01 err = 1.4586531713e-10 time = 1.39 sec -[ Info: DMRG 37: obj = -2.549098968636e+01 err = 1.2550327335e-10 time = 0.49 sec -[ Info: DMRG 38: obj = -2.549098968636e+01 err = 1.0805248647e-10 time = 0.41 sec -[ Info: DMRG conv 39: obj = -2.549098968636e+01 err = 9.3084268982e-11 time = 22.02 sec +[ Info: DMRG init: obj = -1.947519345847e+01 err = 7.5797e-02 +[ Info: DMRG 1: obj = -2.549098874118e+01 err = 9.1376055611e-03 time = 3.54 sec +[ Info: DMRG 2: obj = -2.549098968631e+01 err = 1.1687746206e-06 time = 0.99 sec +[ Info: DMRG 3: obj = -2.549098968636e+01 err = 1.5592263716e-07 time = 1.21 sec +[ Info: DMRG 4: obj = -2.549098968636e+01 err = 1.3434579843e-08 time = 0.62 sec +[ Info: DMRG 5: obj = -2.549098968636e+01 err = 5.9955426595e-09 time = 0.69 sec +[ Info: DMRG 6: obj = -2.549098968636e+01 err = 4.4303868458e-09 time = 0.62 sec +[ Info: DMRG 7: obj = -2.549098968636e+01 err = 3.2709219988e-09 time = 0.62 sec +[ Info: DMRG 8: obj = -2.549098968636e+01 err = 2.3764855423e-09 time = 0.63 sec +[ Info: DMRG 9: obj = -2.549098968636e+01 err = 1.7190997617e-09 time = 0.55 sec +[ Info: DMRG 10: obj = -2.549098968636e+01 err = 1.2467309044e-09 time = 0.55 sec +[ Info: DMRG 11: obj = -2.549098968636e+01 err = 9.1015412118e-10 time = 0.59 sec +[ Info: DMRG 12: obj = -2.549098968636e+01 err = 7.7947998405e-10 time = 0.61 sec +[ Info: DMRG 13: obj = -2.549098968636e+01 err = 6.8300434048e-10 time = 0.61 sec +[ Info: DMRG 14: obj = -2.549098968636e+01 err = 5.9682764775e-10 time = 0.67 sec +[ Info: DMRG 15: obj = -2.549098968636e+01 err = 5.2006815587e-10 time = 0.75 sec +[ Info: DMRG 16: obj = -2.549098968636e+01 err = 4.5212003413e-10 time = 1.13 sec +[ Info: DMRG 17: obj = -2.549098968636e+01 err = 3.9234738860e-10 time = 0.67 sec +[ Info: DMRG 18: obj = -2.549098968636e+01 err = 3.4004240572e-10 time = 0.56 sec +[ Info: DMRG 19: obj = -2.549098968636e+01 err = 2.9445818985e-10 time = 0.59 sec +[ Info: DMRG 20: obj = -2.549098968636e+01 err = 2.5485045197e-10 time = 0.62 sec +[ Info: DMRG 21: obj = -2.549098968636e+01 err = 2.2050912081e-10 time = 0.59 sec +[ Info: DMRG 22: obj = -2.549098968636e+01 err = 1.9077747364e-10 time = 0.60 sec +[ Info: DMRG 23: obj = -2.549098968636e+01 err = 1.6506123889e-10 time = 0.59 sec +[ Info: DMRG 24: obj = -2.549098968636e+01 err = 1.4283083882e-10 time = 0.59 sec +[ Info: DMRG 25: obj = -2.549098968636e+01 err = 1.2361940779e-10 time = 0.65 sec +[ Info: DMRG 26: obj = -2.549098968636e+01 err = 1.0701848483e-10 time = 0.59 sec +[ Info: DMRG conv 27: obj = -2.549098968636e+01 err = 9.2672581484e-11 time = 21.17 sec ```` @@ -260,7 +251,7 @@ 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 857d0f4a6..a9142c06d 100644 --- a/docs/src/examples/quantum1d/1.ising-cft/main.ipynb +++ b/docs/src/examples/quantum1d/1.ising-cft/main.ipynb @@ -24,7 +24,7 @@ { "cell_type": "markdown", "source": [ - "The hamiltonian is defined on a finite lattice with periodic boundary conditions,\n", + "The Hamiltonian is defined on a finite lattice with periodic boundary conditions,\n", "which can be implemented as follows:" ], "metadata": {} @@ -45,8 +45,8 @@ "## Exact diagonalisation\n", "\n", "In MPSKit, there is support for exact diagonalisation by leveraging the fact that applying\n", - "the hamiltonian to an untruncated MPS will result in an effective hamiltonian on the center\n", - "site which implements the action of the entire hamiltonian. Thus, optimizing the middle\n", + "the Hamiltonian to an untruncated MPS will result in an effective Hamiltonian on the center\n", + "site which implements the action of the entire Hamiltonian. Thus, optimizing the middle\n", "tensor is equivalent to optimixing a state in the entire Hilbert space, as all other tensors\n", "are just unitary matrices that mix the basis." ], @@ -108,7 +108,7 @@ { "cell_type": "markdown", "source": [ - "We can then calculate the momentum of the groundstate as the expectation value of this\n", + "We can then calculate the momentum of the ground state as the expectation value of this\n", "operator. However, there is a subtlety because of the degeneracies in the energy\n", "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", @@ -203,7 +203,7 @@ { "cell_type": "markdown", "source": [ - "Excitations on top of the groundstate can be found through the use of the quasiparticle\n", + "Excitations on top of the ground state can be found through the use of the quasiparticle\n", "ansatz. This returns quasiparticle states, which can be converted to regular `FiniteMPS`\n", "objects." ], @@ -260,11 +260,11 @@ "file_extension": ".jl", "mimetype": "application/julia", "name": "julia", - "version": "1.12.0" + "version": "1.12.4" }, "kernelspec": { "name": "julia-1.12", - "display_name": "Julia 1.12.0", + "display_name": "Julia 1.12.4", "language": "julia" } }, diff --git a/docs/src/examples/quantum1d/2.haldane/index.md b/docs/src/examples/quantum1d/2.haldane/index.md index 3b3195e15..5cf05ecdc 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.7989253589480493 +0.798925358948053 ```` We can go even further and doublecheck the claim that ``S = 1`` is an edge excitation, by plotting the energy density. @@ -102,7 +102,7 @@ f = fit(Ls .^ (-2), ΔEs, 1) ```` ```` -0.4517340158584577 +0.4517340158584072 ```` ````julia @@ -140,23 +140,23 @@ println("minimum @k = $(kspace[idx]):\t ΔE = $(ΔE)") ```` ```` +[ Info: Found excitations for momentum = 1.8849555921538759 +[ Info: Found excitations for momentum = 1.6755160819145563 +[ Info: Found excitations for momentum = 2.0943951023931953 +[ Info: Found excitations for momentum = 1.4660765716752369 [ 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 = 1.6755160819145563 -[ Info: Found excitations for momentum = 1.8849555921538759 -[ Info: Found excitations for momentum = 2.0943951023931953 +[ Info: Found excitations for momentum = 2.9321531433504737 +[ Info: Found excitations for momentum = 0.41887902047863906 [ Info: Found excitations for momentum = 2.303834612632515 +[ Info: Found excitations for momentum = 0.6283185307179586 [ Info: Found excitations for momentum = 2.5132741228718345 +[ Info: Found excitations for momentum = 1.0471975511965976 [ 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.41047924851920886 +[ Info: Found excitations for momentum = 0.8377580409572781 +minimum @k = 3.141592653589793: ΔE = 0.410479248594856 ```` diff --git a/docs/src/examples/quantum1d/2.haldane/main.ipynb b/docs/src/examples/quantum1d/2.haldane/main.ipynb index 5af3641d7..605d4e4cc 100644 --- a/docs/src/examples/quantum1d/2.haldane/main.ipynb +++ b/docs/src/examples/quantum1d/2.haldane/main.ipynb @@ -22,13 +22,13 @@ { "cell_type": "markdown", "source": [ - "The Heisenberg model is defined by the following hamiltonian:\n", + "The Heisenberg model is defined by the following Hamiltonian:\n", "\n", "$$\n", "H = -J∑_{⟨i,j⟩} (X_i X_j + Y_i Y_j + Z_i Z_j)\n", "$$\n", "\n", - "This hamiltonian has an SU(2) symmetry, which we can enforce by using SU(2)-symmetric tensors:" + "This Hamiltonian has an SU(2) symmetry, which we can enforce by using SU(2)-symmetric tensors:" ], "metadata": {} }, @@ -49,7 +49,7 @@ "## Finite size extrapolation\n", "\n", "We can start the analysis using finite-size methods.\n", - "The groundstate of this model can be approximated using finite MPS through the use of DMRG.\n", + "The ground state of this model can be approximated using finite MPS through the use of DMRG.\n", "\n", "The typical way to find excited states is to minimize the energy while adding an error term\n", "$$λ \\left|gs\\right> \\left< gs\\right|$$\n", @@ -144,7 +144,7 @@ "A much nicer way of obtaining the Haldane gap is by working directly in the thermodynamic limit.\n", "As was already hinted at by the edge modes, this model is in a non-trivial SPT phase.\n", "Thus, care must be taken when selecting the symmetry sectors.\n", - "The groundstate has half-integer edge modes, thus the virtual spaces must also all carry half-integer charges.\n", + "The ground state has half-integer edge modes, thus the virtual spaces must also all carry half-integer charges.\n", "\n", "In contrast with the finite size case, we now should specify a momentum label to the excitations.\n", "This way, it is possible to scan the dispersion relation over the entire momentum space." @@ -187,11 +187,11 @@ "file_extension": ".jl", "mimetype": "application/julia", "name": "julia", - "version": "1.12.0" + "version": "1.12.4" }, "kernelspec": { "name": "julia-1.12", - "display_name": "Julia 1.12.0", + "display_name": "Julia 1.12.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 85b5468bc..114dff78f 100644 --- a/docs/src/examples/quantum1d/3.ising-dqpt/index.md +++ b/docs/src/examples/quantum1d/3.ising-dqpt/index.md @@ -39,13 +39,13 @@ H₀ = transverse_field_ising(FiniteChain(L); g = -0.5) ```` ```` -[ Info: DMRG init: obj = +9.814858589284e+00 err = 6.2050e-01 -[ Info: DMRG 1: obj = -2.040021714839e+01 err = 2.3171232950e-03 time = 0.03 sec -[ Info: DMRG 2: obj = -2.040021715177e+01 err = 2.1683169804e-07 time = 0.02 sec -[ Info: DMRG 3: obj = -2.040021782399e+01 err = 3.6729918367e-05 time = 0.06 sec -[ Info: DMRG 4: obj = -2.040021786693e+01 err = 1.5200955823e-06 time = 0.03 sec -[ Info: DMRG 5: obj = -2.040021786703e+01 err = 2.8646117610e-07 time = 0.03 sec -[ Info: DMRG conv 6: obj = -2.040021786703e+01 err = 7.0484493208e-11 time = 0.18 sec +[ Info: DMRG init: obj = +1.000110728540e+01 err = 1.6485e-01 +[ Info: DMRG 1: obj = -2.040021714732e+01 err = 2.0639142207e-02 time = 0.05 sec +[ Info: DMRG 2: obj = -2.040021715172e+01 err = 4.6522148018e-07 time = 0.02 sec +[ Info: DMRG 3: obj = -2.040021780162e+01 err = 3.6601066271e-05 time = 0.07 sec +[ Info: DMRG 4: obj = -2.040021786698e+01 err = 1.6025039584e-06 time = 0.04 sec +[ Info: DMRG 5: obj = -2.040021786703e+01 err = 1.6324890445e-07 time = 0.03 sec +[ Info: DMRG conv 6: obj = -2.040021786703e+01 err = 7.7060474652e-11 time = 0.23 sec ```` @@ -111,14 +111,13 @@ H₀ = transverse_field_ising(; g = -0.5) ```` ```` -[ Info: VUMPS init: obj = +4.868298549128e-01 err = 3.9079e-01 -[ Info: VUMPS 1: obj = -1.058483590962e+00 err = 7.5407605847e-02 time = 7.02 sec -[ Info: VUMPS 2: obj = -1.063544286238e+00 err = 3.1160989378e-04 time = 0.01 sec -[ Info: VUMPS 3: obj = -1.063544409943e+00 err = 5.9017642250e-06 time = 0.01 sec -[ Info: VUMPS 4: obj = -1.063544409973e+00 err = 1.5639699231e-07 time = 0.00 sec -[ Info: VUMPS 5: obj = -1.063544409973e+00 err = 5.8340710042e-09 time = 0.00 sec -[ Info: VUMPS 6: obj = -1.063544409973e+00 err = 3.3087676594e-10 time = 0.00 sec -[ Info: VUMPS conv 7: obj = -1.063544409973e+00 err = 3.9690136612e-11 time = 7.05 sec +[ Info: VUMPS init: obj = +4.821692686834e-01 err = 3.7170e-01 +[ Info: VUMPS 1: obj = -1.062759988951e+00 err = 2.3248148483e-02 time = 4.48 sec +[ Info: VUMPS 2: obj = -1.063544409807e+00 err = 1.3008047000e-05 time = 0.01 sec +[ Info: VUMPS 3: obj = -1.063544409973e+00 err = 1.7599807614e-07 time = 0.01 sec +[ Info: VUMPS 4: obj = -1.063544409973e+00 err = 9.6477135912e-09 time = 0.01 sec +[ Info: VUMPS 5: obj = -1.063544409973e+00 err = 4.4909944158e-10 time = 0.01 sec +[ Info: VUMPS conv 6: obj = -1.063544409973e+00 err = 4.9916284444e-11 time = 4.52 sec ```` @@ -130,7 +129,7 @@ dot(ψ₀, ψ₀) ```` ```` -1.000000000000001 - 2.1950801504054652e-16im +0.9999999999999987 + 1.749244813012143e-16im ```` so the Loschmidt 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 b0e82357c..195a9171e 100644 --- a/docs/src/examples/quantum1d/3.ising-dqpt/main.ipynb +++ b/docs/src/examples/quantum1d/3.ising-dqpt/main.ipynb @@ -24,7 +24,7 @@ "source": [ "Dynamical quantum phase transitions (DQPT in short) are signatures of equilibrium phase transitions in a dynamical quantity - the Loschmidt echo.\n", "This quantity is given by $L(t) = \\frac{-2}{N} ln(| < \\psi(t) | \\psi(0) > |) $ where $N$ is the system size.\n", - "One typically starts from a groundstate and then quenches the hamiltonian to a different point.\n", + "One typically starts from a ground state and then quenches the Hamiltonian to a different point.\n", "Non analycities in the Loschmidt echo are called 'dynamical quantum phase transitions'.\n", "\n", "In the mentioned paper they work with\n", @@ -36,7 +36,7 @@ "The outline of the tutorial is as follows. We will pick $g₀ = 0.5$, $g₁ = 2.0$, and perform the time evolution at different system sizes and compare with the thermodynamic limit.\n", "For those $g$ we expect non-analyticities to occur at $t_n ≈ 2.35 (n + 1/2)$.\n", "\n", - "First we construct the hamiltonian in mpo form, and obtain the pre-quenched groundstate:" + "First we construct the Hamiltonian in MPO form, and obtain the pre-quenched ground state:" ], "metadata": {} }, @@ -93,7 +93,7 @@ { "cell_type": "markdown", "source": [ - "\"envs\" is a kind of cache object that keeps track of all environments in `ψ`. It is often advantageous to re-use the environment, so that mpskit doesn't need to recalculate everything.\n", + "\"envs\" is a kind of cache object that keeps track of all environments in `ψ`. It is often advantageous to re-use the environment, so that MPSKit doesn't need to recalculate everything.\n", "\n", "Putting it all together, we get" ], @@ -139,7 +139,7 @@ "source": [ "## Infinite MPS quenching\n", "\n", - "Similarly we could start with an initial infinite state and find the pre-quench groundstate:" + "Similarly we could start with an initial infinite state and find the pre-quench ground state:" ], "metadata": {} }, @@ -191,9 +191,8 @@ { "cell_type": "markdown", "source": [ - "This time we cannot use a two-site scheme to grow the bond dimension, as this isn't implemented (yet).\n", - "Instead, we have to make use of the changebonds machinery.\n", - "Multiple algorithms are available, but we will only focus on `OptimalEpand()`.\n", + "We make use of the `changebonds` machinery to grow the bond dimension. This can also be achieved through a two-site scheme.\n", + "Multiple algorithms are available, but we will only focus on `OptimalExpand()`.\n", "Growing the bond dimension by $5$ can be done by calling:" ], "metadata": {} @@ -285,11 +284,11 @@ "file_extension": ".jl", "mimetype": "application/julia", "name": "julia", - "version": "1.12.0" + "version": "1.12.4" }, "kernelspec": { "name": "julia-1.12", - "display_name": "Julia 1.12.0", + "display_name": "Julia 1.12.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 33c25c263..1776379be 100644 --- a/docs/src/examples/quantum1d/4.xxz-heisenberg/index.md +++ b/docs/src/examples/quantum1d/4.xxz-heisenberg/index.md @@ -26,10 +26,12 @@ H = heisenberg_XXX(; spin = 1 // 2) ```` ```` -single site InfiniteMPOHamiltonian{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[1]: 3×1×1×3 JordanMPOTensor(((ℂ^1 ⊞ ℂ^3 ⊞ ℂ^1) ⊗ ⊞(ℂ^2)) ← (⊞(ℂ^2) ⊗ (ℂ^1 ⊞ ℂ^3 ⊞ ℂ^1))) -╵ ⋮ +1-site InfiniteMPOHamiltonian(ComplexF64, TensorKit.ComplexSpace) with maximal dimension 5: +| ⋮ +| (ℂ^1 ⊞ ℂ^3 ⊞ ℂ^1) +┼─[1]─ ℂ^2 +│ (ℂ^1 ⊞ ℂ^3 ⊞ ℂ^1) +| ⋮ ```` @@ -40,11 +42,12 @@ state = InfiniteMPS(2, 20) ```` ```` -single site InfiniteMPS: -│ ⋮ -│ C[1]: TensorMap(ℂ^20 ← ℂ^20) -├── AL[1]: TensorMap((ℂ^20 ⊗ ℂ^2) ← ℂ^20) -│ ⋮ +1-site InfiniteMPS(ComplexF64, TensorKit.ComplexSpace) with maximal dimension 20: +| ⋮ +| ℂ^20 +├─[1]─ ℂ^2 +│ ℂ^20 +| ⋮ ```` @@ -55,208 +58,208 @@ groundstate, cache, delta = find_groundstate(state, H, VUMPS()); ```` ```` -[ Info: VUMPS init: obj = +2.499995549467e-01 err = 1.8944e-03 -[ Info: VUMPS 1: obj = -1.637326006220e-01 err = 3.4380167886e-01 time = 0.02 sec -[ Info: VUMPS 2: obj = +4.022469296605e-02 err = 3.6124904819e-01 time = 0.02 sec -[ Info: VUMPS 3: obj = -1.451781276155e-01 err = 3.7525843588e-01 time = 0.02 sec -[ Info: VUMPS 4: obj = -3.297874086946e-01 err = 3.0974732609e-01 time = 0.03 sec -[ Info: VUMPS 5: obj = +8.480531607776e-02 err = 4.2311149602e-01 time = 0.02 sec -[ Info: VUMPS 6: obj = -4.682900333525e-02 err = 3.9435625697e-01 time = 0.03 sec -[ Info: VUMPS 7: obj = -2.038899540195e-01 err = 3.7392529829e-01 time = 0.02 sec -[ Info: VUMPS 8: obj = -1.567891704735e-01 err = 3.9440419459e-01 time = 0.02 sec -[ Info: VUMPS 9: obj = +1.798128191908e-02 err = 3.8563879385e-01 time = 0.02 sec -[ Info: VUMPS 10: obj = -1.627417051509e-02 err = 4.0179419114e-01 time = 0.02 sec -[ Info: VUMPS 11: obj = +1.971444286808e-01 err = 3.1358642027e-01 time = 0.02 sec -[ Info: VUMPS 12: obj = -1.544732012749e-01 err = 3.9113131940e-01 time = 0.03 sec -[ Info: VUMPS 13: obj = -5.143757541054e-02 err = 3.6158026015e-01 time = 0.02 sec -[ Info: VUMPS 14: obj = -2.378431775937e-01 err = 3.6684775265e-01 time = 0.15 sec -[ Info: VUMPS 15: obj = -2.335972628792e-01 err = 3.5833112567e-01 time = 0.02 sec -[ Info: VUMPS 16: obj = -3.513837130591e-01 err = 3.1447160955e-01 time = 0.03 sec -[ Info: VUMPS 17: obj = -1.269547785447e-01 err = 4.0560571685e-01 time = 0.03 sec -[ Info: VUMPS 18: obj = -5.526651111523e-02 err = 3.9877496947e-01 time = 0.02 sec -[ Info: VUMPS 19: obj = -1.736900085076e-01 err = 3.7193659518e-01 time = 0.02 sec -[ Info: VUMPS 20: obj = -2.041914181357e-01 err = 3.5751923396e-01 time = 0.02 sec -[ Info: VUMPS 21: obj = -3.693165536045e-01 err = 2.7182712254e-01 time = 0.03 sec -[ Info: VUMPS 22: obj = -2.011720710870e-01 err = 3.6676865727e-01 time = 0.02 sec -[ Info: VUMPS 23: obj = -1.633950485920e-01 err = 4.1193224347e-01 time = 0.02 sec -[ Info: VUMPS 24: obj = -4.626612529510e-02 err = 4.0924314691e-01 time = 0.03 sec -[ Info: VUMPS 25: obj = -1.281611280890e-01 err = 3.9044251366e-01 time = 0.02 sec -[ Info: VUMPS 26: obj = -3.368232288231e-01 err = 3.2045176746e-01 time = 0.02 sec -[ Info: VUMPS 27: obj = -5.824695053770e-02 err = 3.9279884760e-01 time = 0.02 sec -[ Info: VUMPS 28: obj = -1.584622209274e-01 err = 3.9976472773e-01 time = 0.11 sec -[ Info: VUMPS 29: obj = -1.002822977314e-01 err = 3.8803754727e-01 time = 0.03 sec -[ Info: VUMPS 30: obj = +7.400348453278e-03 err = 3.9371296126e-01 time = 0.03 sec -[ Info: VUMPS 31: obj = -1.177431901110e-01 err = 3.9678690292e-01 time = 0.02 sec -[ Info: VUMPS 32: obj = -1.648756258411e-01 err = 3.8382308262e-01 time = 0.03 sec -[ Info: VUMPS 33: obj = -1.486770429634e-01 err = 3.6066618432e-01 time = 0.02 sec -[ Info: VUMPS 34: obj = -3.115577559394e-01 err = 3.4055409376e-01 time = 0.03 sec -[ Info: VUMPS 35: obj = -3.431437826768e-01 err = 3.1792959691e-01 time = 0.03 sec -[ Info: VUMPS 36: obj = -4.164170397906e-01 err = 2.0323232703e-01 time = 0.03 sec -[ Info: VUMPS 37: obj = -4.347019355186e-01 err = 1.2218577552e-01 time = 0.04 sec -[ Info: VUMPS 38: obj = +1.126876972824e-02 err = 3.8542612776e-01 time = 0.03 sec -[ Info: VUMPS 39: obj = -1.054351314040e-01 err = 4.0525492793e-01 time = 0.06 sec -[ Info: VUMPS 40: obj = -1.414301602807e-01 err = 4.0916217647e-01 time = 0.03 sec -[ Info: VUMPS 41: obj = -1.893357319314e-01 err = 3.7913998013e-01 time = 0.02 sec -[ Info: VUMPS 42: obj = -2.638139371875e-01 err = 3.5727613224e-01 time = 0.03 sec -[ Info: VUMPS 43: obj = -1.087944796522e-01 err = 3.5097164414e-01 time = 0.02 sec -[ Info: VUMPS 44: obj = -1.601735087524e-01 err = 3.9348691903e-01 time = 0.02 sec -[ Info: VUMPS 45: obj = -3.599764269854e-01 err = 2.9117042431e-01 time = 0.03 sec -[ Info: VUMPS 46: obj = -4.307479018240e-01 err = 1.4590877399e-01 time = 0.03 sec -[ Info: VUMPS 47: obj = +1.409551265290e-01 err = 3.7003839916e-01 time = 0.02 sec -[ Info: VUMPS 48: obj = -1.237917373142e-01 err = 3.8254863769e-01 time = 0.02 sec -[ Info: VUMPS 49: obj = -1.627580205525e-01 err = 3.7367727861e-01 time = 0.02 sec -[ Info: VUMPS 50: obj = -2.769089084247e-01 err = 3.5236198797e-01 time = 0.05 sec -[ Info: VUMPS 51: obj = -2.343888370451e-01 err = 3.6427105579e-01 time = 0.02 sec -[ Info: VUMPS 52: obj = -2.248449151603e-01 err = 3.7977874636e-01 time = 0.03 sec -[ Info: VUMPS 53: obj = -7.550510266759e-02 err = 4.0150761461e-01 time = 0.02 sec -[ Info: VUMPS 54: obj = -1.157957101584e-01 err = 4.1227979674e-01 time = 0.03 sec -[ Info: VUMPS 55: obj = -3.546254055503e-02 err = 3.9957978395e-01 time = 0.03 sec -[ Info: VUMPS 56: obj = -1.129978671789e-01 err = 4.0034284390e-01 time = 0.02 sec -[ Info: VUMPS 57: obj = -2.329800814496e-01 err = 3.7187892563e-01 time = 0.03 sec -[ Info: VUMPS 58: obj = +2.751650323643e-02 err = 4.0247044284e-01 time = 0.03 sec -[ Info: VUMPS 59: obj = +1.917009780601e-01 err = 3.3210927196e-01 time = 0.02 sec -[ Info: VUMPS 60: obj = -1.263549808391e-01 err = 3.9315953588e-01 time = 0.07 sec -[ Info: VUMPS 61: obj = -8.215407597511e-02 err = 3.9533232172e-01 time = 0.02 sec -[ Info: VUMPS 62: obj = -2.064545516441e-01 err = 3.7374649705e-01 time = 0.02 sec -[ Info: VUMPS 63: obj = -2.034027935096e-01 err = 3.9761679200e-01 time = 0.03 sec -[ Info: VUMPS 64: obj = -2.846911108688e-03 err = 3.7662169289e-01 time = 0.03 sec -[ Info: VUMPS 65: obj = +7.553415672703e-02 err = 3.7835378290e-01 time = 0.02 sec -[ Info: VUMPS 66: obj = -1.171790490724e-01 err = 3.5935068701e-01 time = 0.03 sec -[ Info: VUMPS 67: obj = -1.497117506752e-01 err = 3.7395334449e-01 time = 0.03 sec -[ Info: VUMPS 68: obj = -2.039793547678e-01 err = 3.7667227814e-01 time = 0.03 sec -[ Info: VUMPS 69: obj = -2.089140931903e-01 err = 3.6415207135e-01 time = 0.02 sec -[ Info: VUMPS 70: obj = -1.134944918293e-01 err = 3.6436812433e-01 time = 0.06 sec -[ Info: VUMPS 71: obj = -1.157762777756e-01 err = 3.8567331411e-01 time = 0.02 sec -[ Info: VUMPS 72: obj = -2.394261860644e-01 err = 3.6057044079e-01 time = 0.02 sec -[ Info: VUMPS 73: obj = -3.342400580410e-01 err = 3.1317816493e-01 time = 0.02 sec -[ Info: VUMPS 74: obj = -3.904134305503e-01 err = 2.5856894641e-01 time = 0.03 sec -[ Info: VUMPS 75: obj = +8.449256685246e-02 err = 3.8483086625e-01 time = 0.03 sec -[ Info: VUMPS 76: obj = -3.933130622801e-02 err = 4.0208977801e-01 time = 0.03 sec -[ Info: VUMPS 77: obj = +8.653869040795e-02 err = 3.8802124228e-01 time = 0.03 sec -[ Info: VUMPS 78: obj = +9.587848560625e-02 err = 3.5722735528e-01 time = 0.02 sec -[ Info: VUMPS 79: obj = -5.563236275270e-02 err = 4.1092165670e-01 time = 0.03 sec -[ Info: VUMPS 80: obj = -6.486296834251e-03 err = 3.9951995891e-01 time = 0.06 sec -[ Info: VUMPS 81: obj = -1.209081100544e-01 err = 3.6666868713e-01 time = 0.02 sec -[ Info: VUMPS 82: obj = -2.521478349167e-01 err = 3.5309480836e-01 time = 0.02 sec -[ Info: VUMPS 83: obj = -3.517052113063e-01 err = 3.1461835111e-01 time = 0.03 sec -[ Info: VUMPS 84: obj = -3.299898797334e-01 err = 3.1486702189e-01 time = 0.03 sec -[ Info: VUMPS 85: obj = -3.493479897924e-01 err = 3.1357877281e-01 time = 0.03 sec -[ Info: VUMPS 86: obj = -9.705041624306e-02 err = 3.8584904079e-01 time = 0.02 sec -[ Info: VUMPS 87: obj = -6.719503369858e-02 err = 3.8802177181e-01 time = 0.02 sec -[ Info: VUMPS 88: obj = +7.083873868280e-02 err = 3.8329406568e-01 time = 0.02 sec -[ Info: VUMPS 89: obj = +1.427527172763e-01 err = 3.6625276364e-01 time = 0.02 sec -[ Info: VUMPS 90: obj = -2.155800529197e-02 err = 4.3111181536e-01 time = 0.03 sec -[ Info: VUMPS 91: obj = -3.535552416184e-02 err = 4.1031265505e-01 time = 0.05 sec -[ Info: VUMPS 92: obj = +4.038017896823e-02 err = 3.8143678194e-01 time = 0.02 sec -[ Info: VUMPS 93: obj = -1.670680877582e-01 err = 3.8295252948e-01 time = 0.02 sec -[ Info: VUMPS 94: obj = -9.389358669508e-02 err = 4.0264966465e-01 time = 0.04 sec -[ Info: VUMPS 95: obj = -6.539125878566e-02 err = 3.6997274129e-01 time = 0.03 sec -[ Info: VUMPS 96: obj = -2.076204063368e-01 err = 3.8010942600e-01 time = 0.02 sec -[ Info: VUMPS 97: obj = -2.793821365956e-01 err = 3.4935730821e-01 time = 0.02 sec -[ Info: VUMPS 98: obj = -7.458031879980e-02 err = 3.8897663409e-01 time = 0.02 sec -[ Info: VUMPS 99: obj = -2.819439484489e-01 err = 3.6706683600e-01 time = 0.03 sec -[ Info: VUMPS 100: obj = -1.001823219839e-01 err = 3.7624054699e-01 time = 0.02 sec -[ Info: VUMPS 101: obj = -6.394432799715e-02 err = 3.9177673956e-01 time = 0.06 sec -[ Info: VUMPS 102: obj = -1.465694002120e-01 err = 3.8378109832e-01 time = 0.02 sec -[ Info: VUMPS 103: obj = -2.405698991621e-01 err = 3.7543291978e-01 time = 0.02 sec -[ Info: VUMPS 104: obj = -3.513218895760e-01 err = 3.0159411514e-01 time = 0.02 sec -[ Info: VUMPS 105: obj = -1.023723784151e-01 err = 3.6048545911e-01 time = 0.03 sec -[ Info: VUMPS 106: obj = -1.254960641829e-01 err = 3.5705227068e-01 time = 0.03 sec -[ Info: VUMPS 107: obj = -2.439341345938e-01 err = 3.6666868249e-01 time = 0.03 sec -[ Info: VUMPS 108: obj = -1.822440707693e-01 err = 3.9830893422e-01 time = 0.03 sec -[ Info: VUMPS 109: obj = -2.501473062742e-01 err = 3.5974890197e-01 time = 0.02 sec -[ Info: VUMPS 110: obj = -2.907188314950e-01 err = 3.3385674670e-01 time = 0.03 sec -[ Info: VUMPS 111: obj = -3.880803577729e-01 err = 2.6012194034e-01 time = 0.05 sec -[ Info: VUMPS 112: obj = -3.135177544690e-01 err = 3.5859183499e-01 time = 0.03 sec -[ Info: VUMPS 113: obj = -3.239598468980e-01 err = 3.4821622589e-01 time = 0.03 sec -[ Info: VUMPS 114: obj = -4.092283711908e-01 err = 2.1862198400e-01 time = 0.03 sec -[ Info: VUMPS 115: obj = -4.220328070787e-01 err = 1.8801072163e-01 time = 0.03 sec -[ Info: VUMPS 116: obj = +1.547122686700e-02 err = 4.0231050978e-01 time = 0.03 sec -[ Info: VUMPS 117: obj = -7.971332985234e-02 err = 3.7654696853e-01 time = 0.02 sec -[ Info: VUMPS 118: obj = +1.462167544480e-02 err = 3.7464919429e-01 time = 0.02 sec -[ Info: VUMPS 119: obj = -2.507716805320e-01 err = 3.7118480512e-01 time = 0.02 sec -[ Info: VUMPS 120: obj = -3.534787471790e-01 err = 3.1688927832e-01 time = 0.05 sec -[ Info: VUMPS 121: obj = -1.108707448425e-01 err = 4.1451059429e-01 time = 0.02 sec -[ Info: VUMPS 122: obj = -1.842124515165e-01 err = 3.8528792301e-01 time = 0.02 sec -[ Info: VUMPS 123: obj = -2.998975814447e-01 err = 3.4023194006e-01 time = 0.02 sec -[ Info: VUMPS 124: obj = -1.186804385815e-01 err = 4.1899473179e-01 time = 0.04 sec -[ Info: VUMPS 125: obj = -1.948725735956e-01 err = 3.6431703182e-01 time = 0.02 sec -[ Info: VUMPS 126: obj = +3.273407063701e-02 err = 3.4643539632e-01 time = 0.02 sec -[ Info: VUMPS 127: obj = -2.266211833470e-01 err = 3.4825725987e-01 time = 0.03 sec -[ Info: VUMPS 128: obj = -3.116731680576e-01 err = 3.4010142993e-01 time = 0.03 sec -[ Info: VUMPS 129: obj = -2.826698629538e-01 err = 3.6929476876e-01 time = 0.03 sec -[ Info: VUMPS 130: obj = -1.572495813700e-01 err = 3.9971170739e-01 time = 0.06 sec -[ Info: VUMPS 131: obj = -1.732268711440e-02 err = 3.6936890832e-01 time = 0.02 sec -[ Info: VUMPS 132: obj = -4.775357014229e-02 err = 3.9962098427e-01 time = 0.03 sec -[ Info: VUMPS 133: obj = -2.431698382557e-01 err = 3.6404028082e-01 time = 0.02 sec -[ Info: VUMPS 134: obj = -2.821932941813e-01 err = 3.5445708250e-01 time = 0.02 sec -[ Info: VUMPS 135: obj = -9.024883449741e-02 err = 3.7670364392e-01 time = 0.03 sec -[ Info: VUMPS 136: obj = -2.262433195377e-01 err = 3.8532078673e-01 time = 0.03 sec -[ Info: VUMPS 137: obj = -8.163381818312e-02 err = 3.8947952861e-01 time = 0.02 sec -[ Info: VUMPS 138: obj = -7.989801779312e-03 err = 3.8868571809e-01 time = 0.02 sec -[ Info: VUMPS 139: obj = +2.634778860890e-02 err = 3.6599239785e-01 time = 0.02 sec -[ Info: VUMPS 140: obj = -1.307675815427e-01 err = 3.4665317823e-01 time = 0.03 sec -[ Info: VUMPS 141: obj = -3.217088121498e-01 err = 3.2313971501e-01 time = 0.05 sec -[ Info: VUMPS 142: obj = -4.246143494691e-01 err = 1.5885779904e-01 time = 0.03 sec -[ Info: VUMPS 143: obj = -8.774044006206e-02 err = 3.8014515500e-01 time = 0.03 sec -[ Info: VUMPS 144: obj = -8.327043512249e-02 err = 3.7969588944e-01 time = 0.02 sec -[ Info: VUMPS 145: obj = -1.325167744344e-01 err = 3.8919621062e-01 time = 0.02 sec -[ Info: VUMPS 146: obj = -1.662916750343e-01 err = 3.4606874768e-01 time = 0.03 sec -[ Info: VUMPS 147: obj = -1.901394048799e-01 err = 3.8421371861e-01 time = 0.03 sec -[ Info: VUMPS 148: obj = -2.601825005552e-02 err = 3.9297149298e-01 time = 0.03 sec -[ Info: VUMPS 149: obj = -2.376758303565e-01 err = 3.6086607615e-01 time = 0.03 sec -[ Info: VUMPS 150: obj = -2.272499924592e-01 err = 3.7479126222e-01 time = 0.05 sec -[ Info: VUMPS 151: obj = -2.036843938804e-01 err = 3.9925741694e-01 time = 0.03 sec -[ Info: VUMPS 152: obj = -2.313892960788e-01 err = 3.7972643403e-01 time = 0.03 sec -[ Info: VUMPS 153: obj = -1.986998610932e-01 err = 3.9191236797e-01 time = 0.03 sec -[ Info: VUMPS 154: obj = -2.108517047947e-01 err = 3.7580583074e-01 time = 0.03 sec -[ Info: VUMPS 155: obj = -3.176605702081e-01 err = 3.2923368977e-01 time = 0.03 sec -[ Info: VUMPS 156: obj = +9.435660504635e-02 err = 3.8427690796e-01 time = 0.02 sec -[ Info: VUMPS 157: obj = -7.759368642086e-02 err = 3.9643333120e-01 time = 0.02 sec -[ Info: VUMPS 158: obj = -3.957699631929e-02 err = 4.0189632984e-01 time = 0.02 sec -[ Info: VUMPS 159: obj = -1.217059048082e-01 err = 3.7324188583e-01 time = 0.03 sec -[ Info: VUMPS 160: obj = -3.009908896146e-02 err = 3.8322255280e-01 time = 0.05 sec -[ Info: VUMPS 161: obj = -7.237850128759e-02 err = 4.0889926487e-01 time = 0.03 sec -[ Info: VUMPS 162: obj = -1.373382697223e-01 err = 4.2014320380e-01 time = 0.02 sec -[ Info: VUMPS 163: obj = -1.099154797758e-01 err = 4.1417886345e-01 time = 0.03 sec -[ Info: VUMPS 164: obj = -1.374913376578e-01 err = 3.5906146231e-01 time = 0.02 sec -[ Info: VUMPS 165: obj = -1.464572634014e-01 err = 3.8610039285e-01 time = 0.03 sec -[ Info: VUMPS 166: obj = -3.497532054532e-01 err = 3.2149526676e-01 time = 0.03 sec -[ Info: VUMPS 167: obj = -3.538682309356e-01 err = 3.0937161510e-01 time = 0.03 sec -[ Info: VUMPS 168: obj = -1.174426478371e-01 err = 4.0570387888e-01 time = 0.02 sec -[ Info: VUMPS 169: obj = -1.003314337439e-01 err = 3.9330106579e-01 time = 0.03 sec -[ Info: VUMPS 170: obj = -1.397014985899e-01 err = 3.8327251771e-01 time = 0.05 sec -[ Info: VUMPS 171: obj = -1.485590666911e-02 err = 3.6193610363e-01 time = 0.02 sec -[ Info: VUMPS 172: obj = -2.841425862542e-01 err = 3.5517293193e-01 time = 0.03 sec -[ Info: VUMPS 173: obj = -3.175669710557e-01 err = 3.2886842202e-01 time = 0.02 sec -[ Info: VUMPS 174: obj = -4.215876012708e-01 err = 1.8272144933e-01 time = 0.03 sec -[ Info: VUMPS 175: obj = +1.080077314978e-01 err = 3.8650612296e-01 time = 0.03 sec -[ Info: VUMPS 176: obj = -1.105754803548e-01 err = 3.7491324259e-01 time = 0.02 sec -[ Info: VUMPS 177: obj = -1.534253812824e-01 err = 3.6869415832e-01 time = 0.02 sec -[ Info: VUMPS 178: obj = -1.418091290806e-01 err = 3.7395180342e-01 time = 0.02 sec -[ Info: VUMPS 179: obj = -2.897274308719e-01 err = 3.5651563652e-01 time = 0.02 sec -[ Info: VUMPS 180: obj = -2.391169842001e-02 err = 3.7325894327e-01 time = 0.05 sec -[ Info: VUMPS 181: obj = -1.721332943936e-01 err = 3.7821121884e-01 time = 0.03 sec -[ Info: VUMPS 182: obj = -1.052070523338e-01 err = 3.8451236128e-01 time = 0.02 sec -[ Info: VUMPS 183: obj = -6.354898914006e-02 err = 4.3678211247e-01 time = 0.02 sec -[ Info: VUMPS 184: obj = +5.549252056570e-02 err = 3.6090995884e-01 time = 0.03 sec -[ Info: VUMPS 185: obj = +5.565063849145e-02 err = 4.0136597605e-01 time = 0.03 sec -[ Info: VUMPS 186: obj = +1.768269685046e-01 err = 3.2651611381e-01 time = 0.02 sec -[ Info: VUMPS 187: obj = +1.593690655282e-04 err = 3.8841746150e-01 time = 0.02 sec -[ Info: VUMPS 188: obj = +4.130061860625e-02 err = 3.7222674602e-01 time = 0.02 sec -[ Info: VUMPS 189: obj = -1.930146741248e-01 err = 3.8117696997e-01 time = 0.03 sec -[ Info: VUMPS 190: obj = -3.519361182836e-01 err = 3.0263989022e-01 time = 0.02 sec -[ Info: VUMPS 191: obj = -2.601579677807e-01 err = 3.6061622733e-01 time = 0.05 sec -[ Info: VUMPS 192: obj = -1.002975257021e-01 err = 4.1048517065e-01 time = 0.02 sec -[ Info: VUMPS 193: obj = -1.563988174307e-01 err = 3.9315686700e-01 time = 0.02 sec -[ Info: VUMPS 194: obj = -9.911970259716e-02 err = 3.7832099886e-01 time = 0.02 sec -[ Info: VUMPS 195: obj = -1.044710830190e-01 err = 4.0631216492e-01 time = 0.03 sec -[ Info: VUMPS 196: obj = -1.858694684809e-01 err = 3.9336892502e-01 time = 0.02 sec -[ Info: VUMPS 197: obj = -3.090115170517e-01 err = 3.5367444183e-01 time = 0.03 sec -[ Info: VUMPS 198: obj = -1.946258306654e-02 err = 3.6965265025e-01 time = 0.02 sec -[ Info: VUMPS 199: obj = -2.238429310745e-01 err = 3.7577524599e-01 time = 0.03 sec -┌ Warning: VUMPS cancel 200: obj = -3.951359685829e-01 err = 2.5736930769e-01 time = 5.62 sec -└ @ MPSKit ~/Projects/MPSKit.jl/src/algorithms/groundstate/vumps.jl:76 +[ Info: VUMPS init: obj = +2.499992657736e-01 err = 2.3659e-03 +[ Info: VUMPS 1: obj = -1.113720211054e-01 err = 3.4044919149e-01 time = 7.00 sec +[ Info: VUMPS 2: obj = -7.257521722207e-02 err = 3.7501747094e-01 time = 0.29 sec +[ Info: VUMPS 3: obj = -4.083234386194e-02 err = 4.0906387666e-01 time = 0.08 sec +[ Info: VUMPS 4: obj = -1.575841447054e-01 err = 3.6917648883e-01 time = 0.03 sec +[ Info: VUMPS 5: obj = -1.988657390936e-01 err = 3.6791056420e-01 time = 0.03 sec +[ Info: VUMPS 6: obj = -3.197413863086e-01 err = 3.3258645953e-01 time = 0.04 sec +[ Info: VUMPS 7: obj = -7.087435895207e-02 err = 3.8002416604e-01 time = 0.03 sec +[ Info: VUMPS 8: obj = -1.207743137267e-01 err = 3.8655168695e-01 time = 0.02 sec +[ Info: VUMPS 9: obj = -2.086746007253e-01 err = 3.8866997285e-01 time = 0.08 sec +[ Info: VUMPS 10: obj = -2.135149067167e-02 err = 3.9998107520e-01 time = 0.02 sec +[ Info: VUMPS 11: obj = +7.890476682466e-02 err = 3.9662516712e-01 time = 0.02 sec +[ Info: VUMPS 12: obj = +8.179378380698e-02 err = 3.6865902322e-01 time = 0.04 sec +[ Info: VUMPS 13: obj = -3.389058286591e-01 err = 3.3592003962e-01 time = 0.07 sec +[ Info: VUMPS 14: obj = -6.936331877678e-02 err = 4.1384132520e-01 time = 0.02 sec +[ Info: VUMPS 15: obj = -2.095570913077e-01 err = 3.7672687907e-01 time = 0.03 sec +[ Info: VUMPS 16: obj = -8.547010074293e-02 err = 4.0253817589e-01 time = 0.02 sec +[ Info: VUMPS 17: obj = -1.633131236641e-01 err = 3.7085752764e-01 time = 0.01 sec +[ Info: VUMPS 18: obj = -2.440483461804e-01 err = 3.5634709811e-01 time = 0.02 sec +[ Info: VUMPS 19: obj = -2.227348354782e-01 err = 3.6670163212e-01 time = 0.07 sec +[ Info: VUMPS 20: obj = -1.647291995688e-01 err = 3.9547116130e-01 time = 0.03 sec +[ Info: VUMPS 21: obj = -2.272417561286e-01 err = 3.6276621476e-01 time = 0.02 sec +[ Info: VUMPS 22: obj = -3.388919076523e-01 err = 3.1579827529e-01 time = 0.02 sec +[ Info: VUMPS 23: obj = -3.915604574764e-01 err = 2.5208735300e-01 time = 0.07 sec +[ Info: VUMPS 24: obj = -2.716952795096e-01 err = 3.4825332817e-01 time = 0.04 sec +[ Info: VUMPS 25: obj = -2.511746372453e-01 err = 3.6715895791e-01 time = 0.03 sec +[ Info: VUMPS 26: obj = -2.731262108428e-02 err = 3.8195721702e-01 time = 0.04 sec +[ Info: VUMPS 27: obj = -8.401004965257e-02 err = 4.1219830163e-01 time = 0.03 sec +[ Info: VUMPS 28: obj = -5.593112558648e-02 err = 3.7991678397e-01 time = 0.07 sec +[ Info: VUMPS 29: obj = -1.039939899261e-01 err = 4.1791684367e-01 time = 0.04 sec +[ Info: VUMPS 30: obj = +3.646055386088e-02 err = 3.6884640674e-01 time = 0.02 sec +[ Info: VUMPS 31: obj = -1.745961021141e-01 err = 3.7601676764e-01 time = 0.02 sec +[ Info: VUMPS 32: obj = +6.528773402406e-03 err = 3.9555076806e-01 time = 0.02 sec +[ Info: VUMPS 33: obj = -6.925592993288e-02 err = 3.7990669074e-01 time = 0.08 sec +[ Info: VUMPS 34: obj = -1.493711283015e-01 err = 3.9233769652e-01 time = 0.02 sec +[ Info: VUMPS 35: obj = -8.418199500312e-02 err = 3.6340380458e-01 time = 0.02 sec +[ Info: VUMPS 36: obj = -2.276410296255e-02 err = 3.7714183538e-01 time = 0.02 sec +[ Info: VUMPS 37: obj = -1.746306414275e-02 err = 4.3071643087e-01 time = 0.03 sec +[ Info: VUMPS 38: obj = -1.078530946206e-01 err = 4.0554246549e-01 time = 0.08 sec +[ Info: VUMPS 39: obj = -4.201824689181e-03 err = 3.8332228055e-01 time = 0.03 sec +[ Info: VUMPS 40: obj = -1.695196954544e-01 err = 4.0646048881e-01 time = 0.03 sec +[ Info: VUMPS 41: obj = -3.991123357010e-01 err = 2.6404691766e-01 time = 0.03 sec +[ Info: VUMPS 42: obj = -3.731927178774e-04 err = 4.0618023904e-01 time = 0.02 sec +[ Info: VUMPS 43: obj = -1.361859934319e-01 err = 3.8083035475e-01 time = 0.07 sec +[ Info: VUMPS 44: obj = -4.351467677687e-02 err = 3.6067072767e-01 time = 0.01 sec +[ Info: VUMPS 45: obj = +2.036960334903e-02 err = 4.0236374802e-01 time = 0.02 sec +[ Info: VUMPS 46: obj = -2.564786609196e-01 err = 3.6550194687e-01 time = 0.02 sec +[ Info: VUMPS 47: obj = -2.794914053758e-02 err = 3.8554065921e-01 time = 0.03 sec +[ Info: VUMPS 48: obj = -1.456926586076e-01 err = 3.8815324507e-01 time = 0.02 sec +[ Info: VUMPS 49: obj = -1.424477861718e-02 err = 3.7358256943e-01 time = 0.06 sec +[ Info: VUMPS 50: obj = -4.121215827673e-02 err = 3.7267277980e-01 time = 0.02 sec +[ Info: VUMPS 51: obj = -3.082552784588e-01 err = 3.3162188490e-01 time = 0.02 sec +[ Info: VUMPS 52: obj = -2.583717562047e-01 err = 3.6414628453e-01 time = 0.03 sec +[ Info: VUMPS 53: obj = +1.080952855580e-01 err = 3.5691712848e-01 time = 0.02 sec +[ Info: VUMPS 54: obj = -2.388241333550e-01 err = 3.6671504144e-01 time = 0.06 sec +[ Info: VUMPS 55: obj = -2.158989432818e-01 err = 3.5252639302e-01 time = 0.02 sec +[ Info: VUMPS 56: obj = -2.846028253281e-01 err = 3.4280944784e-01 time = 0.02 sec +[ Info: VUMPS 57: obj = -2.844144756572e-01 err = 3.5771158117e-01 time = 0.04 sec +[ Info: VUMPS 58: obj = -3.296915296031e-01 err = 3.4124297389e-01 time = 0.07 sec +[ Info: VUMPS 59: obj = -3.428588110821e-01 err = 3.1055543978e-01 time = 0.05 sec +[ Info: VUMPS 60: obj = -3.689645262669e-01 err = 2.6709254795e-01 time = 0.04 sec +[ Info: VUMPS 61: obj = -1.926296113920e-01 err = 3.9051008708e-01 time = 0.04 sec +[ Info: VUMPS 62: obj = -2.083013559580e-01 err = 3.8360309501e-01 time = 0.09 sec +[ Info: VUMPS 63: obj = -2.330751889106e-01 err = 3.5874442855e-01 time = 0.04 sec +[ Info: VUMPS 64: obj = -1.170346482791e-01 err = 3.8356655434e-01 time = 0.03 sec +[ Info: VUMPS 65: obj = -5.200963578316e-02 err = 3.6833174183e-01 time = 0.03 sec +[ Info: VUMPS 66: obj = -2.197891431095e-02 err = 4.0124794245e-01 time = 0.02 sec +[ Info: VUMPS 67: obj = -1.169084572818e-01 err = 3.9681625944e-01 time = 0.07 sec +[ Info: VUMPS 68: obj = -3.187497224181e-01 err = 3.5389927614e-01 time = 0.03 sec +[ Info: VUMPS 69: obj = -5.774837074080e-02 err = 4.0555259098e-01 time = 0.05 sec +[ Info: VUMPS 70: obj = -1.385925677140e-01 err = 3.5902899521e-01 time = 0.02 sec +[ Info: VUMPS 71: obj = -1.618893406115e-01 err = 3.7624779484e-01 time = 0.02 sec +[ Info: VUMPS 72: obj = -1.501259513086e-01 err = 3.9017543257e-01 time = 0.08 sec +[ Info: VUMPS 73: obj = -2.344578406890e-01 err = 3.7629861390e-01 time = 0.02 sec +[ Info: VUMPS 74: obj = -9.559035249277e-02 err = 3.7722338357e-01 time = 0.02 sec +[ Info: VUMPS 75: obj = -3.559605422071e-01 err = 3.0031220801e-01 time = 0.02 sec +[ Info: VUMPS 76: obj = -1.954426659903e-01 err = 3.7999614851e-01 time = 0.02 sec +[ Info: VUMPS 77: obj = -2.469742715213e-01 err = 3.5932526830e-01 time = 0.07 sec +[ Info: VUMPS 78: obj = +1.378130557437e-02 err = 3.7598886038e-01 time = 0.02 sec +[ Info: VUMPS 79: obj = -1.211793491665e-02 err = 3.6307721624e-01 time = 0.02 sec +[ Info: VUMPS 80: obj = -1.280208560452e-01 err = 4.0740026676e-01 time = 0.03 sec +[ Info: VUMPS 81: obj = -2.377846434769e-01 err = 3.7798079849e-01 time = 0.02 sec +[ Info: VUMPS 82: obj = -3.480513403032e-01 err = 3.1094392423e-01 time = 0.06 sec +[ Info: VUMPS 83: obj = -3.231666531453e-01 err = 3.5795237816e-01 time = 0.03 sec +[ Info: VUMPS 84: obj = -3.230617723979e-01 err = 3.5338137372e-01 time = 0.08 sec +[ Info: VUMPS 85: obj = -1.391043781959e-01 err = 3.7699469945e-01 time = 0.02 sec +[ Info: VUMPS 86: obj = -3.434600465182e-01 err = 3.1162050408e-01 time = 0.08 sec +[ Info: VUMPS 87: obj = -2.896970610423e-01 err = 3.4715325464e-01 time = 0.03 sec +[ Info: VUMPS 88: obj = -3.473500689819e-01 err = 3.0203627103e-01 time = 0.03 sec +[ Info: VUMPS 89: obj = -3.920852749773e-01 err = 2.4261851491e-01 time = 0.04 sec +[ Info: VUMPS 90: obj = -2.323120265885e-03 err = 3.9363699208e-01 time = 0.08 sec +[ Info: VUMPS 91: obj = +6.597719659994e-03 err = 3.6444377682e-01 time = 0.02 sec +[ Info: VUMPS 92: obj = +2.842708237132e-02 err = 3.8036913669e-01 time = 0.02 sec +[ Info: VUMPS 93: obj = +8.866987029294e-02 err = 3.6998433554e-01 time = 0.03 sec +[ Info: VUMPS 94: obj = -5.099903063388e-02 err = 3.6527413192e-01 time = 0.02 sec +[ Info: VUMPS 95: obj = -1.724607412056e-01 err = 3.6271311572e-01 time = 0.06 sec +[ Info: VUMPS 96: obj = -1.359620364521e-01 err = 3.7163521224e-01 time = 0.03 sec +[ Info: VUMPS 97: obj = -2.487229212367e-01 err = 3.5559482532e-01 time = 0.03 sec +[ Info: VUMPS 98: obj = -3.363768011366e-01 err = 3.0691942472e-01 time = 0.02 sec +[ Info: VUMPS 99: obj = -2.349190816289e-01 err = 3.8040603712e-01 time = 0.07 sec +[ Info: VUMPS 100: obj = -2.501766365990e-01 err = 3.3899797906e-01 time = 0.02 sec +[ Info: VUMPS 101: obj = -3.944006951029e-01 err = 2.5365404140e-01 time = 0.02 sec +[ Info: VUMPS 102: obj = -2.995031953737e-02 err = 3.9078970987e-01 time = 0.03 sec +[ Info: VUMPS 103: obj = -1.496727054580e-01 err = 3.8887820541e-01 time = 0.03 sec +[ Info: VUMPS 104: obj = -1.651964515037e-01 err = 3.8482074585e-01 time = 0.06 sec +[ Info: VUMPS 105: obj = -2.136080092702e-01 err = 3.7086717496e-01 time = 0.02 sec +[ Info: VUMPS 106: obj = +2.076780492509e-03 err = 3.7971063787e-01 time = 0.02 sec +[ Info: VUMPS 107: obj = -1.801409094443e-01 err = 3.7203262354e-01 time = 0.02 sec +[ Info: VUMPS 108: obj = -2.824910448641e-01 err = 3.5387251865e-01 time = 0.02 sec +[ Info: VUMPS 109: obj = -1.603436420431e-01 err = 3.9841222680e-01 time = 0.07 sec +[ Info: VUMPS 110: obj = +2.717881863352e-02 err = 3.9583169447e-01 time = 0.02 sec +[ Info: VUMPS 111: obj = +3.649162587405e-02 err = 3.2986973402e-01 time = 0.03 sec +[ Info: VUMPS 112: obj = -1.931191472770e-01 err = 3.8832602682e-01 time = 0.02 sec +[ Info: VUMPS 113: obj = -8.595612263509e-02 err = 3.9456614625e-01 time = 0.02 sec +[ Info: VUMPS 114: obj = -2.052894924418e-01 err = 3.7916182480e-01 time = 0.08 sec +[ Info: VUMPS 115: obj = -1.566623408265e-01 err = 3.7572169212e-01 time = 0.02 sec +[ Info: VUMPS 116: obj = -2.697925193842e-01 err = 3.5224872666e-01 time = 0.03 sec +[ Info: VUMPS 117: obj = -1.653547727749e-01 err = 3.6580445508e-01 time = 0.02 sec +[ Info: VUMPS 118: obj = -3.099498644207e-01 err = 3.3813063847e-01 time = 0.02 sec +[ Info: VUMPS 119: obj = -1.249709134848e-01 err = 3.9297526101e-01 time = 0.07 sec +[ Info: VUMPS 120: obj = -1.759792403391e-01 err = 3.9286117967e-01 time = 0.03 sec +[ Info: VUMPS 121: obj = -1.531238685880e-01 err = 3.7741746395e-01 time = 0.02 sec +[ Info: VUMPS 122: obj = -2.607981777457e-02 err = 4.0319496740e-01 time = 0.03 sec +[ Info: VUMPS 123: obj = -6.592234489034e-02 err = 4.0591425677e-01 time = 0.06 sec +[ Info: VUMPS 124: obj = -6.634520512325e-02 err = 3.5661784954e-01 time = 0.02 sec +[ Info: VUMPS 125: obj = -8.713375234031e-02 err = 3.8753642064e-01 time = 0.03 sec +[ Info: VUMPS 126: obj = -1.125202000113e-01 err = 3.9947548666e-01 time = 0.02 sec +[ Info: VUMPS 127: obj = -7.695667111350e-02 err = 3.8370863324e-01 time = 0.02 sec +[ Info: VUMPS 128: obj = -1.626394492619e-01 err = 3.6111268034e-01 time = 0.06 sec +[ Info: VUMPS 129: obj = -3.049948899402e-01 err = 3.4403324168e-01 time = 0.02 sec +[ Info: VUMPS 130: obj = -1.258602194339e-01 err = 3.9449464476e-01 time = 0.03 sec +[ Info: VUMPS 131: obj = -1.104762253188e-01 err = 3.7940891591e-01 time = 0.03 sec +[ Info: VUMPS 132: obj = +1.961617781097e-01 err = 3.5498814846e-01 time = 0.04 sec +[ Info: VUMPS 133: obj = -5.045117087438e-02 err = 3.9852508192e-01 time = 0.06 sec +[ Info: VUMPS 134: obj = -1.456650620925e-01 err = 3.7944737703e-01 time = 0.02 sec +[ Info: VUMPS 135: obj = -8.263613090518e-02 err = 3.9519353274e-01 time = 0.02 sec +[ Info: VUMPS 136: obj = -1.765342573210e-01 err = 3.8095369610e-01 time = 0.02 sec +[ Info: VUMPS 137: obj = -1.324964043868e-01 err = 3.9563235756e-01 time = 0.02 sec +[ Info: VUMPS 138: obj = -5.801528053563e-03 err = 3.9895803654e-01 time = 0.07 sec +[ Info: VUMPS 139: obj = +3.411995954756e-02 err = 3.7915180484e-01 time = 0.02 sec +[ Info: VUMPS 140: obj = -5.243883424994e-02 err = 3.6892061143e-01 time = 0.02 sec +[ Info: VUMPS 141: obj = -2.240214467858e-01 err = 3.6441762460e-01 time = 0.02 sec +[ Info: VUMPS 142: obj = -2.299190776747e-01 err = 3.6831318301e-01 time = 0.03 sec +[ Info: VUMPS 143: obj = -1.957936901383e-01 err = 3.9355611293e-01 time = 0.06 sec +[ Info: VUMPS 144: obj = -3.444003971032e-01 err = 3.1223261647e-01 time = 0.02 sec +[ Info: VUMPS 145: obj = -4.206622251723e-01 err = 1.7712997133e-01 time = 0.03 sec +[ Info: VUMPS 146: obj = -3.152511465401e-01 err = 3.4248220688e-01 time = 0.04 sec +[ Info: VUMPS 147: obj = -2.523905778375e-01 err = 3.7740870895e-01 time = 0.07 sec +[ Info: VUMPS 148: obj = -4.278492688890e-02 err = 3.5715350579e-01 time = 0.02 sec +[ Info: VUMPS 149: obj = -1.324003340205e-01 err = 4.0013050853e-01 time = 0.02 sec +[ Info: VUMPS 150: obj = -1.068985749883e-01 err = 4.0722208139e-01 time = 0.02 sec +[ Info: VUMPS 151: obj = -2.083862397081e-01 err = 3.6731256501e-01 time = 0.02 sec +[ Info: VUMPS 152: obj = -1.440078983653e-01 err = 3.7426441705e-01 time = 0.06 sec +[ Info: VUMPS 153: obj = -2.520133076021e-01 err = 3.7507045206e-01 time = 0.02 sec +[ Info: VUMPS 154: obj = -1.034831359402e-01 err = 3.3452358556e-01 time = 0.02 sec +[ Info: VUMPS 155: obj = -6.618164126260e-02 err = 3.5135937268e-01 time = 0.03 sec +[ Info: VUMPS 156: obj = -8.742999239375e-02 err = 3.8534304068e-01 time = 0.02 sec +[ Info: VUMPS 157: obj = -1.259714263689e-01 err = 3.9638882368e-01 time = 0.07 sec +[ Info: VUMPS 158: obj = -9.181298201303e-02 err = 3.6918229891e-01 time = 0.02 sec +[ Info: VUMPS 159: obj = -3.681891440506e-01 err = 2.7888162226e-01 time = 0.02 sec +[ Info: VUMPS 160: obj = -3.357053644322e-01 err = 3.2800086170e-01 time = 0.03 sec +[ Info: VUMPS 161: obj = -3.215096951929e-03 err = 3.6664749423e-01 time = 0.07 sec +[ Info: VUMPS 162: obj = +2.646903240101e-03 err = 3.8676335475e-01 time = 0.02 sec +[ Info: VUMPS 163: obj = +7.709386278480e-02 err = 3.5056948614e-01 time = 0.01 sec +[ Info: VUMPS 164: obj = -1.446953476068e-01 err = 3.4180912676e-01 time = 0.02 sec +[ Info: VUMPS 165: obj = -3.581100223225e-01 err = 3.1025051513e-01 time = 0.02 sec +[ Info: VUMPS 166: obj = -2.406814929481e-01 err = 3.8072841015e-01 time = 0.07 sec +[ Info: VUMPS 167: obj = -3.531848910098e-01 err = 3.1053005603e-01 time = 0.02 sec +[ Info: VUMPS 168: obj = -3.297669874782e-01 err = 3.2648804299e-01 time = 0.03 sec +[ Info: VUMPS 169: obj = -1.800883981792e-01 err = 3.7239998435e-01 time = 0.03 sec +[ Info: VUMPS 170: obj = -4.056383454871e-02 err = 3.8301203586e-01 time = 0.03 sec +[ Info: VUMPS 171: obj = -4.448126315315e-02 err = 3.9634205142e-01 time = 0.07 sec +[ Info: VUMPS 172: obj = -3.109186702795e-01 err = 3.4936775901e-01 time = 0.02 sec +[ Info: VUMPS 173: obj = -1.687932549901e-01 err = 3.9713516798e-01 time = 0.02 sec +[ Info: VUMPS 174: obj = -1.703811557590e-01 err = 3.7182890670e-01 time = 0.02 sec +[ Info: VUMPS 175: obj = -1.360221590893e-01 err = 3.7750383659e-01 time = 0.02 sec +[ Info: VUMPS 176: obj = -1.746982448560e-01 err = 3.8728329904e-01 time = 0.07 sec +[ Info: VUMPS 177: obj = -9.378022642433e-02 err = 3.9508912967e-01 time = 0.02 sec +[ Info: VUMPS 178: obj = -1.283340267182e-01 err = 3.8629991015e-01 time = 0.02 sec +[ Info: VUMPS 179: obj = -8.976213980219e-02 err = 4.0495071985e-01 time = 0.02 sec +[ Info: VUMPS 180: obj = -2.669351107647e-01 err = 3.3526515364e-01 time = 0.07 sec +[ Info: VUMPS 181: obj = -3.341709529585e-02 err = 3.8011147844e-01 time = 0.02 sec +[ Info: VUMPS 182: obj = -3.109125460864e-01 err = 3.4567156611e-01 time = 0.03 sec +[ Info: VUMPS 183: obj = -2.183171915980e-01 err = 3.7747243244e-01 time = 0.02 sec +[ Info: VUMPS 184: obj = -2.372842132147e-01 err = 3.7374609557e-01 time = 0.03 sec +[ Info: VUMPS 185: obj = -2.543753578738e-01 err = 3.5816792692e-01 time = 0.07 sec +[ Info: VUMPS 186: obj = -2.469096804733e-01 err = 3.7131984627e-01 time = 0.03 sec +[ Info: VUMPS 187: obj = -1.830902139527e-01 err = 3.9043241343e-01 time = 0.02 sec +[ Info: VUMPS 188: obj = -1.574766929467e-01 err = 3.7319600051e-01 time = 0.03 sec +[ Info: VUMPS 189: obj = -6.991335970984e-02 err = 3.9222107611e-01 time = 0.07 sec +[ Info: VUMPS 190: obj = -1.792846142089e-01 err = 3.6356547438e-01 time = 0.02 sec +[ Info: VUMPS 191: obj = -1.225042448135e-01 err = 3.7360289247e-01 time = 0.02 sec +[ Info: VUMPS 192: obj = -2.129551247072e-01 err = 3.8299381445e-01 time = 0.02 sec +[ Info: VUMPS 193: obj = -8.283800637018e-02 err = 3.8500403241e-01 time = 0.02 sec +[ Info: VUMPS 194: obj = -2.897858859525e-01 err = 3.5007431093e-01 time = 0.06 sec +[ Info: VUMPS 195: obj = -5.369820439974e-02 err = 4.1837283102e-01 time = 0.02 sec +[ Info: VUMPS 196: obj = -2.800880298694e-01 err = 3.7189978262e-01 time = 0.03 sec +[ Info: VUMPS 197: obj = -1.669323399699e-01 err = 3.9481332800e-01 time = 0.03 sec +[ Info: VUMPS 198: obj = -6.212304176188e-02 err = 4.2526435790e-01 time = 0.03 sec +[ Info: VUMPS 199: obj = -4.204616169557e-02 err = 3.7541303170e-01 time = 0.07 sec +┌ Warning: VUMPS cancel 200: obj = -1.383333715979e-01 err = 3.7699700566e-01 time = 14.12 sec +└ @ MPSKit /home/ldevos/LocalProjects/MPSKit.jl/src/algorithms/groundstate/vumps.jl:76 ```` @@ -269,9 +272,9 @@ groundstate, cache, delta = find_groundstate(state, H, GradientGrassmann(; maxit ```` ```` -[ Info: CG: initializing with f = 0.249999554947, ‖∇f‖ = 1.3395e-03 -┌ Warning: CG: not converged to requested tol after 20 iterations and time 6.30 s: f = -0.441802293068, ‖∇f‖ = 1.0800e-02 -└ @ OptimKit ~/.julia/packages/OptimKit/G6i79/src/cg.jl:172 +[ Info: CG: initializing with f = 2.499992657736e-01, ‖∇f‖ = 1.6729e-03 +┌ Warning: CG: not converged to requested tol after 20 iterations and time 8.29 s: f = -4.426111048892e-01, ‖∇f‖ = 5.5480e-03 +└ @ OptimKit ~/.julia/packages/OptimKit/OEwMx/src/cg.jl:174 ```` @@ -283,7 +286,7 @@ transferplot(groundstate, groundstate) ```` ```@raw html - + ``` We can clearly see multiple eigenvalues close to the unit circle. @@ -300,12 +303,14 @@ state = InfiniteMPS(fill(2, 2), fill(20, 2)) ```` ```` -2-site InfiniteMPS: -│ ⋮ -│ C[2]: TensorMap(ℂ^20 ← ℂ^20) -├── AL[2]: TensorMap((ℂ^20 ⊗ ℂ^2) ← ℂ^20) -├── AL[1]: TensorMap((ℂ^20 ⊗ ℂ^2) ← ℂ^20) -│ ⋮ +2-site InfiniteMPS(ComplexF64, TensorKit.ComplexSpace) with maximal dimension 20: +| ⋮ +| ℂ^20 +├─[2]─ ℂ^2 +│ ℂ^20 +├─[1]─ ℂ^2 +│ ℂ^20 +| ⋮ ```` @@ -322,108 +327,108 @@ groundstate, envs, delta = find_groundstate( ```` ```` -[ Info: VUMPS init: obj = +4.988938408716e-01 err = 6.6329e-02 -[ Info: VUMPS 1: obj = -4.739248140491e-01 err = 3.3816850792e-01 time = 0.03 sec -[ Info: VUMPS 2: obj = -8.743783163834e-01 err = 8.7782281077e-02 time = 0.03 sec -[ Info: VUMPS 3: obj = -8.852289662056e-01 err = 1.1261717077e-02 time = 0.03 sec -[ Info: VUMPS 4: obj = -8.859341540347e-01 err = 6.0581249357e-03 time = 0.12 sec -[ Info: VUMPS 5: obj = -8.861158625188e-01 err = 4.1548303785e-03 time = 0.03 sec -[ Info: VUMPS 6: obj = -8.861822061550e-01 err = 2.8484872369e-03 time = 0.03 sec -[ Info: VUMPS 7: obj = -8.862103323128e-01 err = 2.1781142069e-03 time = 0.03 sec -[ Info: VUMPS 8: obj = -8.862234788867e-01 err = 1.6559901778e-03 time = 0.03 sec -[ Info: VUMPS 9: obj = -8.862298897521e-01 err = 1.3551274661e-03 time = 0.03 sec -[ Info: VUMPS 10: obj = -8.862331118301e-01 err = 1.1002219058e-03 time = 0.04 sec -[ Info: VUMPS 11: obj = -8.862347480843e-01 err = 9.4695562130e-04 time = 0.09 sec -[ Info: VUMPS 12: obj = -8.862355951169e-01 err = 8.0634081582e-04 time = 0.04 sec -[ Info: VUMPS 13: obj = -8.862360364673e-01 err = 7.1458186619e-04 time = 0.03 sec -[ Info: VUMPS 14: obj = -8.862362718663e-01 err = 6.3113350829e-04 time = 0.04 sec -[ Info: VUMPS 15: obj = -8.862363996891e-01 err = 5.6139764728e-04 time = 0.04 sec -[ Info: VUMPS 16: obj = -8.862364715964e-01 err = 5.0511718667e-04 time = 0.04 sec -[ Info: VUMPS 17: obj = -8.862365136698e-01 err = 4.5105971790e-04 time = 0.06 sec -[ Info: VUMPS 18: obj = -8.862365394867e-01 err = 4.0773238562e-04 time = 0.03 sec -[ Info: VUMPS 19: obj = -8.862365562591e-01 err = 3.6503807001e-04 time = 0.04 sec -[ Info: VUMPS 20: obj = -8.862365676678e-01 err = 3.3014036465e-04 time = 0.03 sec -[ Info: VUMPS 21: obj = -8.862365758645e-01 err = 2.9598570227e-04 time = 0.04 sec -[ Info: VUMPS 22: obj = -8.862365819255e-01 err = 2.6762287043e-04 time = 0.07 sec -[ Info: VUMPS 23: obj = -8.862365865819e-01 err = 2.4021341209e-04 time = 0.03 sec -[ Info: VUMPS 24: obj = -8.862365902034e-01 err = 2.1707752879e-04 time = 0.03 sec -[ Info: VUMPS 25: obj = -8.862365930830e-01 err = 1.9507319179e-04 time = 0.03 sec -[ Info: VUMPS 26: obj = -8.862365953820e-01 err = 1.7616283550e-04 time = 0.03 sec -[ Info: VUMPS 27: obj = -8.862365972387e-01 err = 1.5851786128e-04 time = 0.03 sec -[ Info: VUMPS 28: obj = -8.862365987416e-01 err = 1.4302884556e-04 time = 0.07 sec -[ Info: VUMPS 29: obj = -8.862365999644e-01 err = 1.2889703895e-04 time = 0.03 sec -[ Info: VUMPS 30: obj = -8.862366009623e-01 err = 1.1617835138e-04 time = 0.03 sec -[ Info: VUMPS 31: obj = -8.862366017775e-01 err = 1.0488359393e-04 time = 0.03 sec -[ Info: VUMPS 32: obj = -8.862366024472e-01 err = 9.4408469738e-05 time = 0.03 sec -[ Info: VUMPS 33: obj = -8.862366029961e-01 err = 8.5404044738e-05 time = 0.03 sec -[ Info: VUMPS 34: obj = -8.862366034498e-01 err = 7.6748331848e-05 time = 0.03 sec -[ Info: VUMPS 35: obj = -8.862366038231e-01 err = 6.9594627182e-05 time = 0.06 sec -[ Info: VUMPS 36: obj = -8.862366041338e-01 err = 6.2415914486e-05 time = 0.03 sec -[ Info: VUMPS 37: obj = -8.862366043907e-01 err = 5.6758538061e-05 time = 0.03 sec -[ Info: VUMPS 38: obj = -8.862366046064e-01 err = 5.0780009127e-05 time = 0.03 sec -[ Info: VUMPS 39: obj = -8.862366047860e-01 err = 4.6333072643e-05 time = 0.03 sec -[ Info: VUMPS 40: obj = -8.862366049382e-01 err = 4.1331207549e-05 time = 0.03 sec -[ Info: VUMPS 41: obj = -8.862366050662e-01 err = 3.7863899058e-05 time = 0.06 sec -[ Info: VUMPS 42: obj = -8.862366051761e-01 err = 3.3658336653e-05 time = 0.03 sec -[ Info: VUMPS 43: obj = -8.862366052696e-01 err = 3.0983965257e-05 time = 0.03 sec -[ Info: VUMPS 44: obj = -8.862366053511e-01 err = 2.7523802341e-05 time = 0.04 sec -[ Info: VUMPS 45: obj = -8.862366054215e-01 err = 2.5397047125e-05 time = 0.03 sec -[ Info: VUMPS 46: obj = -8.862366054839e-01 err = 2.2595519406e-05 time = 0.03 sec -[ Info: VUMPS 47: obj = -8.862366055388e-01 err = 2.0862977631e-05 time = 0.06 sec -[ Info: VUMPS 48: obj = -8.862366055883e-01 err = 1.8599136674e-05 time = 0.03 sec -[ Info: VUMPS 49: obj = -8.862366056327e-01 err = 1.7188042263e-05 time = 0.03 sec -[ Info: VUMPS 50: obj = -8.862366056736e-01 err = 1.5364204857e-05 time = 0.04 sec -[ Info: VUMPS 51: obj = -8.862366057110e-01 err = 1.4215375363e-05 time = 0.03 sec -[ Info: VUMPS 52: obj = -8.862366057460e-01 err = 1.2752194421e-05 time = 0.04 sec -[ Info: VUMPS 53: obj = -8.862366057787e-01 err = 1.1817843549e-05 time = 0.06 sec -[ Info: VUMPS 54: obj = -8.862366058097e-01 err = 1.0651271784e-05 time = 0.03 sec -[ Info: VUMPS 55: obj = -8.862366058392e-01 err = 9.8929892467e-06 time = 0.04 sec -[ Info: VUMPS 56: obj = -8.862366058675e-01 err = 8.9706820132e-06 time = 0.03 sec -[ Info: VUMPS 57: obj = -8.862366058948e-01 err = 8.3567048894e-06 time = 0.04 sec -[ Info: VUMPS 58: obj = -8.862366059214e-01 err = 7.6361746845e-06 time = 0.04 sec -[ Info: VUMPS 59: obj = -8.862366059472e-01 err = 7.1407423616e-06 time = 0.07 sec -[ Info: VUMPS 60: obj = -8.862366059725e-01 err = 6.5866324690e-06 time = 0.04 sec -[ Info: VUMPS 61: obj = -8.862366059973e-01 err = 6.1886994474e-06 time = 0.04 sec -[ Info: VUMPS 62: obj = -8.862366060219e-01 err = 5.7708691344e-06 time = 0.04 sec -[ Info: VUMPS 63: obj = -8.862366060460e-01 err = 5.4526696713e-06 time = 0.04 sec -[ Info: VUMPS 64: obj = -8.862366060700e-01 err = 5.1458563859e-06 time = 0.07 sec -[ Info: VUMPS 65: obj = -8.862366060938e-01 err = 4.8920859232e-06 time = 0.03 sec -[ Info: VUMPS 66: obj = -8.862366061174e-01 err = 4.6741015849e-06 time = 0.03 sec -[ Info: VUMPS 67: obj = -8.862366061409e-01 err = 4.4720328824e-06 time = 0.03 sec -[ Info: VUMPS 68: obj = -8.862366061643e-01 err = 4.3236853227e-06 time = 0.04 sec -[ Info: VUMPS 69: obj = -8.862366061876e-01 err = 4.1624035800e-06 time = 0.04 sec -[ Info: VUMPS 70: obj = -8.862366062109e-01 err = 4.0670654185e-06 time = 0.06 sec -[ Info: VUMPS 71: obj = -8.862366062341e-01 err = 3.9372522985e-06 time = 0.04 sec -[ Info: VUMPS 72: obj = -8.862366062573e-01 err = 3.8814470100e-06 time = 0.04 sec -[ Info: VUMPS 73: obj = -8.862366062805e-01 err = 3.7755570661e-06 time = 0.04 sec -[ Info: VUMPS 74: obj = -8.862366063037e-01 err = 3.7485088165e-06 time = 0.04 sec -[ Info: VUMPS 75: obj = -8.862366063269e-01 err = 3.6605911620e-06 time = 0.04 sec -[ Info: VUMPS 76: obj = -8.862366063501e-01 err = 3.6540367594e-06 time = 0.06 sec -[ Info: VUMPS 77: obj = -8.862366063733e-01 err = 3.5795056881e-06 time = 0.04 sec -[ Info: VUMPS 78: obj = -8.862366063965e-01 err = 3.5873001441e-06 time = 0.04 sec -[ Info: VUMPS 79: obj = -8.862366064197e-01 err = 3.5226970394e-06 time = 0.04 sec -[ Info: VUMPS 80: obj = -8.862366064429e-01 err = 3.5403874329e-06 time = 0.04 sec -[ Info: VUMPS 81: obj = -8.862366064662e-01 err = 3.4831436727e-06 time = 0.07 sec -[ Info: VUMPS 82: obj = -8.862366064895e-01 err = 3.5075667400e-06 time = 0.03 sec -[ Info: VUMPS 83: obj = -8.862366065128e-01 err = 3.4557926414e-06 time = 0.04 sec -[ Info: VUMPS 84: obj = -8.862366065361e-01 err = 3.4847339467e-06 time = 0.04 sec -[ Info: VUMPS 85: obj = -8.862366065595e-01 err = 3.4370450912e-06 time = 0.04 sec -[ Info: VUMPS 86: obj = -8.862366065829e-01 err = 3.4689709047e-06 time = 0.04 sec -[ Info: VUMPS 87: obj = -8.862366066063e-01 err = 3.4243527295e-06 time = 0.06 sec -[ Info: VUMPS 88: obj = -8.862366066297e-01 err = 3.4582089166e-06 time = 0.04 sec -[ Info: VUMPS 89: obj = -8.862366066532e-01 err = 3.4159150211e-06 time = 0.03 sec -[ Info: VUMPS 90: obj = -8.862366066767e-01 err = 3.4509811256e-06 time = 0.04 sec -[ Info: VUMPS 91: obj = -8.862366067002e-01 err = 3.4104587388e-06 time = 0.04 sec -[ Info: VUMPS 92: obj = -8.862366067238e-01 err = 3.4462444730e-06 time = 0.07 sec -[ Info: VUMPS 93: obj = -8.862366067473e-01 err = 3.4093962434e-06 time = 0.03 sec -[ Info: VUMPS 94: obj = -8.862366067709e-01 err = 3.4432535858e-06 time = 0.03 sec -[ Info: VUMPS 95: obj = -8.862366067945e-01 err = 3.4135991515e-06 time = 0.03 sec -[ Info: VUMPS 96: obj = -8.862366068182e-01 err = 3.4414722655e-06 time = 0.04 sec -[ Info: VUMPS 97: obj = -8.862366068418e-01 err = 3.4175821732e-06 time = 0.04 sec -[ Info: VUMPS 98: obj = -8.862366068655e-01 err = 3.4405116799e-06 time = 0.06 sec -[ Info: VUMPS 99: obj = -8.862366068892e-01 err = 3.4212410568e-06 time = 0.04 sec -┌ Warning: VUMPS cancel 100: obj = -8.862366069129e-01 err = 3.4400872838e-06 time = 4.02 sec -└ @ MPSKit ~/Projects/MPSKit.jl/src/algorithms/groundstate/vumps.jl:76 +[ Info: VUMPS init: obj = +4.987085387825e-01 err = 7.4815e-02 +[ Info: VUMPS 1: obj = -4.070483296483e-02 err = 3.8421207912e-01 time = 0.42 sec +[ Info: VUMPS 2: obj = -8.588956105873e-01 err = 1.4311389089e-01 time = 0.04 sec +[ Info: VUMPS 3: obj = -8.846913958875e-01 err = 1.6288696714e-02 time = 0.03 sec +[ Info: VUMPS 4: obj = -8.858742904256e-01 err = 6.3271162473e-03 time = 0.03 sec +[ Info: VUMPS 5: obj = -8.860994898111e-01 err = 4.7506122055e-03 time = 0.07 sec +[ Info: VUMPS 6: obj = -8.861800643367e-01 err = 2.9027017342e-03 time = 0.04 sec +[ Info: VUMPS 7: obj = -8.862112267309e-01 err = 2.2025348252e-03 time = 0.05 sec +[ Info: VUMPS 8: obj = -8.862256494201e-01 err = 1.5437281695e-03 time = 0.08 sec +[ Info: VUMPS 9: obj = -8.862322270149e-01 err = 1.2285599184e-03 time = 0.05 sec +[ Info: VUMPS 10: obj = -8.862354505034e-01 err = 9.5416020728e-04 time = 0.09 sec +[ Info: VUMPS 11: obj = -8.862369964136e-01 err = 7.5140552439e-04 time = 0.06 sec +[ Info: VUMPS 12: obj = -8.862377649153e-01 err = 6.2947047843e-04 time = 0.06 sec +[ Info: VUMPS 13: obj = -8.862381482190e-01 err = 5.3567155770e-04 time = 0.08 sec +[ Info: VUMPS 14: obj = -8.862383501305e-01 err = 4.9726432030e-04 time = 0.05 sec +[ Info: VUMPS 15: obj = -8.862384635741e-01 err = 4.7373205998e-04 time = 0.04 sec +[ Info: VUMPS 16: obj = -8.862385365416e-01 err = 4.7834381662e-04 time = 0.08 sec +[ Info: VUMPS 17: obj = -8.862385914623e-01 err = 4.9010937476e-04 time = 0.06 sec +[ Info: VUMPS 18: obj = -8.862386397747e-01 err = 5.1562261789e-04 time = 0.09 sec +[ Info: VUMPS 19: obj = -8.862386887080e-01 err = 5.4370412689e-04 time = 0.07 sec +[ Info: VUMPS 20: obj = -8.862387403517e-01 err = 5.8200261169e-04 time = 0.06 sec +[ Info: VUMPS 21: obj = -8.862388009252e-01 err = 6.1775277631e-04 time = 0.11 sec +[ Info: VUMPS 22: obj = -8.862388688746e-01 err = 6.6580030857e-04 time = 0.06 sec +[ Info: VUMPS 23: obj = -8.862389551333e-01 err = 7.0184437921e-04 time = 0.06 sec +[ Info: VUMPS 24: obj = -8.862390538865e-01 err = 7.5533855251e-04 time = 0.09 sec +[ Info: VUMPS 25: obj = -8.862391870753e-01 err = 7.7812725098e-04 time = 0.08 sec +[ Info: VUMPS 26: obj = -8.862393378986e-01 err = 8.2535209132e-04 time = 0.10 sec +[ Info: VUMPS 27: obj = -8.862395466874e-01 err = 8.0853992627e-04 time = 0.05 sec +[ Info: VUMPS 28: obj = -8.862397673707e-01 err = 8.2479835919e-04 time = 0.06 sec +[ Info: VUMPS 29: obj = -8.862400528237e-01 err = 7.3793510189e-04 time = 0.07 sec +[ Info: VUMPS 30: obj = -8.862403146804e-01 err = 6.9781267895e-04 time = 0.05 sec +[ Info: VUMPS 31: obj = -8.862405935152e-01 err = 5.6033545341e-04 time = 0.06 sec +[ Info: VUMPS 32: obj = -8.862408144890e-01 err = 4.7977846041e-04 time = 0.09 sec +[ Info: VUMPS 33: obj = -8.862410054255e-01 err = 3.7069947399e-04 time = 0.06 sec +[ Info: VUMPS 34: obj = -8.862411519178e-01 err = 3.0069863398e-04 time = 0.08 sec +[ Info: VUMPS 35: obj = -8.862412707892e-01 err = 2.4560744504e-04 time = 0.06 sec +[ Info: VUMPS 36: obj = -8.862413665434e-01 err = 2.0272830003e-04 time = 0.06 sec +[ Info: VUMPS 37: obj = -8.862414453558e-01 err = 1.7669114346e-04 time = 0.09 sec +[ Info: VUMPS 38: obj = -8.862415101392e-01 err = 1.4966363352e-04 time = 0.06 sec +[ Info: VUMPS 39: obj = -8.862415632632e-01 err = 1.3429888391e-04 time = 0.09 sec +[ Info: VUMPS 40: obj = -8.862416064083e-01 err = 1.1453939487e-04 time = 0.06 sec +[ Info: VUMPS 41: obj = -8.862416411106e-01 err = 1.0327940645e-04 time = 0.09 sec +[ Info: VUMPS 42: obj = -8.862416687158e-01 err = 8.8267214320e-05 time = 0.06 sec +[ Info: VUMPS 43: obj = -8.862416904652e-01 err = 7.8860063320e-05 time = 0.06 sec +[ Info: VUMPS 44: obj = -8.862417074438e-01 err = 6.7541649023e-05 time = 0.08 sec +[ Info: VUMPS 45: obj = -8.862417205990e-01 err = 5.9561454565e-05 time = 0.04 sec +[ Info: VUMPS 46: obj = -8.862417307234e-01 err = 5.0969707803e-05 time = 0.08 sec +[ Info: VUMPS 47: obj = -8.862417384749e-01 err = 4.4565392480e-05 time = 0.05 sec +[ Info: VUMPS 48: obj = -8.862417443833e-01 err = 3.8093480996e-05 time = 0.06 sec +[ Info: VUMPS 49: obj = -8.862417488723e-01 err = 3.3125644419e-05 time = 0.07 sec +[ Info: VUMPS 50: obj = -8.862417522738e-01 err = 2.8298971982e-05 time = 0.06 sec +[ Info: VUMPS 51: obj = -8.862417548472e-01 err = 2.4528466837e-05 time = 0.10 sec +[ Info: VUMPS 52: obj = -8.862417567915e-01 err = 2.0959047274e-05 time = 0.04 sec +[ Info: VUMPS 53: obj = -8.862417582599e-01 err = 1.8137306177e-05 time = 0.03 sec +[ Info: VUMPS 54: obj = -8.862417593686e-01 err = 1.5513404897e-05 time = 0.06 sec +[ Info: VUMPS 55: obj = -8.862417602060e-01 err = 1.3420069417e-05 time = 0.05 sec +[ Info: VUMPS 56: obj = -8.862417608389e-01 err = 1.1498110058e-05 time = 0.08 sec +[ Info: VUMPS 57: obj = -8.862417613176e-01 err = 9.9527379200e-06 time = 0.05 sec +[ Info: VUMPS 58: obj = -8.862417616801e-01 err = 8.5468069101e-06 time = 0.09 sec +[ Info: VUMPS 59: obj = -8.862417619550e-01 err = 7.4079343001e-06 time = 0.04 sec +[ Info: VUMPS 60: obj = -8.862417621636e-01 err = 6.3789249797e-06 time = 0.05 sec +[ Info: VUMPS 61: obj = -8.862417623223e-01 err = 5.5392243963e-06 time = 0.08 sec +[ Info: VUMPS 62: obj = -8.862417624431e-01 err = 4.7843861025e-06 time = 0.05 sec +[ Info: VUMPS 63: obj = -8.862417625353e-01 err = 4.1638877981e-06 time = 0.09 sec +[ Info: VUMPS 64: obj = -8.862417626056e-01 err = 3.6081493661e-06 time = 0.06 sec +[ Info: VUMPS 65: obj = -8.862417626595e-01 err = 3.1478809644e-06 time = 0.06 sec +[ Info: VUMPS 66: obj = -8.862417627008e-01 err = 2.7367561575e-06 time = 0.09 sec +[ Info: VUMPS 67: obj = -8.862417627324e-01 err = 2.3936907479e-06 time = 0.06 sec +[ Info: VUMPS 68: obj = -8.862417627568e-01 err = 2.0878342869e-06 time = 0.09 sec +[ Info: VUMPS 69: obj = -8.862417627755e-01 err = 1.8307028023e-06 time = 0.07 sec +[ Info: VUMPS 70: obj = -8.862417627899e-01 err = 1.6017548673e-06 time = 0.04 sec +[ Info: VUMPS 71: obj = -8.862417628011e-01 err = 1.4078813339e-06 time = 0.08 sec +[ Info: VUMPS 72: obj = -8.862417628097e-01 err = 1.2353984722e-06 time = 0.05 sec +[ Info: VUMPS 73: obj = -8.862417628163e-01 err = 1.0883299189e-06 time = 0.09 sec +[ Info: VUMPS 74: obj = -8.862417628215e-01 err = 9.5754695615e-07 time = 0.06 sec +[ Info: VUMPS 75: obj = -8.862417628254e-01 err = 8.4531618087e-07 time = 0.06 sec +[ Info: VUMPS 76: obj = -8.862417628285e-01 err = 7.4553090808e-07 time = 0.08 sec +[ Info: VUMPS 77: obj = -8.862417628309e-01 err = 6.5939726315e-07 time = 0.05 sec +[ Info: VUMPS 78: obj = -8.862417628328e-01 err = 5.8281301292e-07 time = 0.06 sec +[ Info: VUMPS 79: obj = -8.862417628343e-01 err = 5.1635751307e-07 time = 0.05 sec +[ Info: VUMPS 80: obj = -8.862417628354e-01 err = 4.5725918727e-07 time = 0.05 sec +[ Info: VUMPS 81: obj = -8.862417628363e-01 err = 4.0573779114e-07 time = 0.10 sec +[ Info: VUMPS 82: obj = -8.862417628370e-01 err = 3.5990703959e-07 time = 0.05 sec +[ Info: VUMPS 83: obj = -8.862417628375e-01 err = 3.1978939846e-07 time = 0.09 sec +[ Info: VUMPS 84: obj = -8.862417628379e-01 err = 2.8408974871e-07 time = 0.06 sec +[ Info: VUMPS 85: obj = -8.862417628382e-01 err = 2.5273008172e-07 time = 0.06 sec +[ Info: VUMPS 86: obj = -8.862417628385e-01 err = 2.2481213666e-07 time = 0.08 sec +[ Info: VUMPS 87: obj = -8.862417628387e-01 err = 2.0021375141e-07 time = 0.06 sec +[ Info: VUMPS 88: obj = -8.862417628389e-01 err = 1.7830492076e-07 time = 0.08 sec +[ Info: VUMPS 89: obj = -8.862417628390e-01 err = 1.5895093159e-07 time = 0.06 sec +[ Info: VUMPS 90: obj = -8.862417628391e-01 err = 1.4170464909e-07 time = 0.05 sec +[ Info: VUMPS 91: obj = -8.862417628391e-01 err = 1.2643555469e-07 time = 0.08 sec +[ Info: VUMPS 92: obj = -8.862417628392e-01 err = 1.1282240133e-07 time = 0.06 sec +[ Info: VUMPS 93: obj = -8.862417628393e-01 err = 1.0074688630e-07 time = 0.08 sec +[ Info: VUMPS 94: obj = -8.862417628393e-01 err = 8.9975370500e-08 time = 0.04 sec +[ Info: VUMPS 95: obj = -8.862417628393e-01 err = 8.0404798877e-08 time = 0.06 sec +[ Info: VUMPS 96: obj = -8.862417628393e-01 err = 7.1863230951e-08 time = 0.08 sec +[ Info: VUMPS 97: obj = -8.862417628394e-01 err = 6.4263171589e-08 time = 0.05 sec +[ Info: VUMPS 98: obj = -8.862417628394e-01 err = 5.7476665744e-08 time = 0.07 sec +[ Info: VUMPS 99: obj = -8.862417628394e-01 err = 5.1430713814e-08 time = 0.04 sec +┌ Warning: VUMPS cancel 100: obj = -8.862417628394e-01 err = 4.6029104045e-08 time = 6.81 sec +└ @ MPSKit /home/ldevos/LocalProjects/MPSKit.jl/src/algorithms/groundstate/vumps.jl:76 ```` @@ -438,7 +443,7 @@ entanglementplot(groundstate) ```` ```@raw html - + ``` We see that some eigenvalues clearly belong to a group, and are almost degenerate. @@ -475,7 +480,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:160 +└ @ MPSKit /home/ldevos/LocalProjects/MPSKit.jl/src/states/infinitemps.jl:160 ```` @@ -490,80 +495,74 @@ groundstate, cache, delta = find_groundstate(state, H2, VUMPS(; maxiter = 400, t ```` 52 70 -[ Info: VUMPS init: obj = -1.141245853330e-02 err = 4.0215e-01 -[ Info: VUMPS 1: obj = -8.788989897232e-01 err = 8.3126786202e-02 time = 0.13 sec -[ Info: VUMPS 2: obj = -8.857995903945e-01 err = 6.7432032291e-03 time = 0.03 sec -[ Info: VUMPS 3: obj = -8.861329058794e-01 err = 3.4413406908e-03 time = 0.03 sec -[ Info: VUMPS 4: obj = -8.862249298781e-01 err = 2.0932139560e-03 time = 0.03 sec -[ Info: VUMPS 5: obj = -8.862609030986e-01 err = 1.0275091935e-03 time = 0.03 sec -[ Info: VUMPS 6: obj = -8.862754866598e-01 err = 7.4461572015e-04 time = 0.04 sec -[ Info: VUMPS 7: obj = -8.862819270462e-01 err = 5.8510697988e-04 time = 0.05 sec -[ Info: VUMPS 8: obj = -8.862849475005e-01 err = 4.8578532406e-04 time = 0.04 sec -[ Info: VUMPS 9: obj = -8.862864073427e-01 err = 4.0315541644e-04 time = 0.05 sec -[ Info: VUMPS 10: obj = -8.862871279225e-01 err = 3.5445010854e-04 time = 0.05 sec -[ Info: VUMPS 11: obj = -8.862874891630e-01 err = 3.1822823463e-04 time = 0.05 sec -[ Info: VUMPS 12: obj = -8.862876773439e-01 err = 2.7443998839e-04 time = 0.05 sec -[ Info: VUMPS 13: obj = -8.862877793817e-01 err = 2.2054254501e-04 time = 0.12 sec -[ Info: VUMPS 14: obj = -8.862878353269e-01 err = 1.6605215489e-04 time = 0.04 sec -[ Info: VUMPS 15: obj = -8.862878654264e-01 err = 1.1979575296e-04 time = 0.04 sec -[ Info: VUMPS 16: obj = -8.862878812303e-01 err = 8.4769125422e-05 time = 0.05 sec -[ Info: VUMPS 17: obj = -8.862878894026e-01 err = 5.9778171991e-05 time = 0.05 sec -[ Info: VUMPS 18: obj = -8.862878936100e-01 err = 4.2285093952e-05 time = 0.05 sec -[ Info: VUMPS 19: obj = -8.862878957743e-01 err = 3.0076620714e-05 time = 0.05 sec -[ Info: VUMPS 20: obj = -8.862878968901e-01 err = 2.1497977400e-05 time = 0.05 sec -[ Info: VUMPS 21: obj = -8.862878974666e-01 err = 1.5431141918e-05 time = 0.05 sec -[ Info: VUMPS 22: obj = -8.862878977650e-01 err = 1.1107227296e-05 time = 0.05 sec -[ Info: VUMPS 23: obj = -8.862878979198e-01 err = 8.0121943412e-06 time = 0.05 sec -[ Info: VUMPS 24: obj = -8.862878980003e-01 err = 5.7884597378e-06 time = 0.05 sec -[ Info: VUMPS 25: obj = -8.862878980421e-01 err = 4.1862569480e-06 time = 0.10 sec -[ Info: VUMPS 26: obj = -8.862878980639e-01 err = 3.0295966217e-06 time = 0.05 sec -[ Info: VUMPS 27: obj = -8.862878980752e-01 err = 2.1934508934e-06 time = 0.05 sec -[ Info: VUMPS 28: obj = -8.862878980811e-01 err = 1.5884579833e-06 time = 0.05 sec -[ Info: VUMPS 29: obj = -8.862878980842e-01 err = 1.1504627834e-06 time = 0.05 sec -[ Info: VUMPS 30: obj = -8.862878980858e-01 err = 8.3325912425e-07 time = 0.05 sec -[ Info: VUMPS 31: obj = -8.862878980867e-01 err = 6.0349319385e-07 time = 0.05 sec -[ Info: VUMPS 32: obj = -8.862878980871e-01 err = 4.3705144218e-07 time = 0.05 sec -[ Info: VUMPS 33: obj = -8.862878980874e-01 err = 3.1648278136e-07 time = 0.09 sec -[ Info: VUMPS 34: obj = -8.862878980875e-01 err = 2.2914950477e-07 time = 0.04 sec -[ Info: VUMPS 35: obj = -8.862878980876e-01 err = 1.6589609004e-07 time = 0.04 sec -[ Info: VUMPS 36: obj = -8.862878980876e-01 err = 1.2008839249e-07 time = 0.04 sec -[ Info: VUMPS 37: obj = -8.862878980876e-01 err = 8.6918909626e-08 time = 0.05 sec -[ Info: VUMPS 38: obj = -8.862878980876e-01 err = 6.2903889958e-08 time = 0.05 sec -[ Info: VUMPS 39: obj = -8.862878980876e-01 err = 4.5519032022e-08 time = 0.05 sec -[ Info: VUMPS 40: obj = -8.862878980877e-01 err = 3.2940004401e-08 time = 0.05 sec -[ Info: VUMPS 41: obj = -8.862878980877e-01 err = 2.3832479903e-08 time = 0.09 sec -[ Info: VUMPS 42: obj = -8.862878980877e-01 err = 1.7241026505e-08 time = 0.05 sec -[ Info: VUMPS 43: obj = -8.862878980877e-01 err = 1.2471412997e-08 time = 0.05 sec -[ Info: VUMPS 44: obj = -8.862878980877e-01 err = 9.0205419450e-09 time = 0.05 sec -[ Info: VUMPS 45: obj = -8.862878980877e-01 err = 6.5240508358e-09 time = 0.05 sec -[ Info: VUMPS 46: obj = -8.862878980877e-01 err = 4.7181601548e-09 time = 0.05 sec -[ Info: VUMPS 47: obj = -8.862878980877e-01 err = 3.4119389042e-09 time = 0.05 sec -[ Info: VUMPS 48: obj = -8.862878980877e-01 err = 2.4672054622e-09 time = 0.05 sec -[ Info: VUMPS 49: obj = -8.862878980877e-01 err = 1.7839684924e-09 time = 0.05 sec -[ Info: VUMPS 50: obj = -8.862878980877e-01 err = 1.2898741454e-09 time = 0.09 sec -[ Info: VUMPS 51: obj = -8.862878980877e-01 err = 9.3258970310e-10 time = 0.05 sec -[ Info: VUMPS 52: obj = -8.862878980877e-01 err = 6.7423602771e-10 time = 0.04 sec -[ Info: VUMPS 53: obj = -8.862878980877e-01 err = 4.8743673761e-10 time = 0.04 sec -[ Info: VUMPS 54: obj = -8.862878980877e-01 err = 3.5238100286e-10 time = 0.04 sec -[ Info: VUMPS 55: obj = -8.862878980877e-01 err = 2.5473951258e-10 time = 0.05 sec -[ Info: VUMPS 56: obj = -8.862878980878e-01 err = 1.8414828185e-10 time = 0.05 sec -[ Info: VUMPS 57: obj = -8.862878980878e-01 err = 1.3311459508e-10 time = 0.05 sec -[ Info: VUMPS 58: obj = -8.862878980878e-01 err = 9.6221597943e-11 time = 0.05 sec -[ Info: VUMPS 59: obj = -8.862878980878e-01 err = 6.9553567834e-11 time = 0.09 sec -[ Info: VUMPS 60: obj = -8.862878980878e-01 err = 5.0274813406e-11 time = 0.05 sec -[ Info: VUMPS 61: obj = -8.862878980878e-01 err = 3.6335883698e-11 time = 0.05 sec -[ Info: VUMPS 62: obj = -8.862878980878e-01 err = 2.6254468114e-11 time = 0.04 sec -[ Info: VUMPS 63: obj = -8.862878980878e-01 err = 1.8980171184e-11 time = 0.04 sec -[ Info: VUMPS 64: obj = -8.862878980878e-01 err = 1.3721774296e-11 time = 0.04 sec -[ Info: VUMPS 65: obj = -8.862878980878e-01 err = 9.9185056725e-12 time = 0.04 sec -[ Info: VUMPS 66: obj = -8.862878980878e-01 err = 7.1729564828e-12 time = 0.04 sec -[ Info: VUMPS 67: obj = -8.862878980878e-01 err = 5.1845841611e-12 time = 0.04 sec -[ Info: VUMPS 68: obj = -8.862878980878e-01 err = 3.7451968982e-12 time = 0.04 sec -[ Info: VUMPS 69: obj = -8.862878980878e-01 err = 2.7085301699e-12 time = 0.08 sec -[ Info: VUMPS 70: obj = -8.862878980878e-01 err = 1.9547335304e-12 time = 0.04 sec -[ Info: VUMPS 71: obj = -8.862878980878e-01 err = 1.4119541453e-12 time = 0.04 sec -[ Info: VUMPS 72: obj = -8.862878980878e-01 err = 1.0222255848e-12 time = 0.03 sec -[ Info: VUMPS conv 73: obj = -8.862878980878e-01 err = 7.3834786497e-13 time = 3.71 sec +[ Info: VUMPS init: obj = +8.454690130663e-02 err = 3.6812e-01 +[ Info: VUMPS 1: obj = -8.807747096663e-01 err = 7.4524923622e-02 time = 3.77 sec +[ Info: VUMPS 2: obj = -8.858788324414e-01 err = 6.9171953600e-03 time = 0.05 sec +[ Info: VUMPS 3: obj = -8.861621536444e-01 err = 2.6767683452e-03 time = 0.04 sec +[ Info: VUMPS 4: obj = -8.862392626495e-01 err = 1.6032192901e-03 time = 0.04 sec +[ Info: VUMPS 5: obj = -8.862672547653e-01 err = 9.5323528320e-04 time = 0.05 sec +[ Info: VUMPS 6: obj = -8.862784830480e-01 err = 7.0061763044e-04 time = 0.05 sec +[ Info: VUMPS 7: obj = -8.862834114803e-01 err = 5.7030493713e-04 time = 0.06 sec +[ Info: VUMPS 8: obj = -8.862857129161e-01 err = 4.5154675114e-04 time = 0.15 sec +[ Info: VUMPS 9: obj = -8.862868209497e-01 err = 3.5140725914e-04 time = 0.05 sec +[ Info: VUMPS 10: obj = -8.862873648329e-01 err = 2.6806862728e-04 time = 0.13 sec +[ Info: VUMPS 11: obj = -8.862876338388e-01 err = 2.0070574932e-04 time = 0.11 sec +[ Info: VUMPS 12: obj = -8.862877672196e-01 err = 1.4816530119e-04 time = 0.07 sec +[ Info: VUMPS 13: obj = -8.862878333408e-01 err = 1.0821757653e-04 time = 0.15 sec +[ Info: VUMPS 14: obj = -8.862878660733e-01 err = 7.8417469213e-05 time = 0.06 sec +[ Info: VUMPS 15: obj = -8.862878822601e-01 err = 5.6494124158e-05 time = 0.07 sec +[ Info: VUMPS 16: obj = -8.862878902612e-01 err = 4.0537829957e-05 time = 0.06 sec +[ Info: VUMPS 17: obj = -8.862878942156e-01 err = 2.9004225089e-05 time = 0.07 sec +[ Info: VUMPS 18: obj = -8.862878961706e-01 err = 2.0708366147e-05 time = 0.07 sec +[ Info: VUMPS 19: obj = -8.862878971378e-01 err = 1.4762413368e-05 time = 0.13 sec +[ Info: VUMPS 20: obj = -8.862878976166e-01 err = 1.0511055800e-05 time = 0.03 sec +[ Info: VUMPS 21: obj = -8.862878978539e-01 err = 7.4778223881e-06 time = 0.03 sec +[ Info: VUMPS 22: obj = -8.862878979715e-01 err = 5.3158051331e-06 time = 0.07 sec +[ Info: VUMPS 23: obj = -8.862878980299e-01 err = 3.7764425487e-06 time = 0.07 sec +[ Info: VUMPS 24: obj = -8.862878980589e-01 err = 2.6814072095e-06 time = 0.14 sec +[ Info: VUMPS 25: obj = -8.862878980733e-01 err = 1.9030014616e-06 time = 0.06 sec +[ Info: VUMPS 26: obj = -8.862878980805e-01 err = 1.3500577199e-06 time = 0.06 sec +[ Info: VUMPS 27: obj = -8.862878980841e-01 err = 9.5735794398e-07 time = 0.06 sec +[ Info: VUMPS 28: obj = -8.862878980859e-01 err = 6.7863772480e-07 time = 0.06 sec +[ Info: VUMPS 29: obj = -8.862878980867e-01 err = 4.8090393886e-07 time = 0.14 sec +[ Info: VUMPS 30: obj = -8.862878980872e-01 err = 3.4067956729e-07 time = 0.06 sec +[ Info: VUMPS 31: obj = -8.862878980874e-01 err = 2.4127441738e-07 time = 0.05 sec +[ Info: VUMPS 32: obj = -8.862878980875e-01 err = 1.7082422697e-07 time = 0.05 sec +[ Info: VUMPS 33: obj = -8.862878980876e-01 err = 1.2091935762e-07 time = 0.06 sec +[ Info: VUMPS 34: obj = -8.862878980876e-01 err = 8.5574898934e-08 time = 0.05 sec +[ Info: VUMPS 35: obj = -8.862878980876e-01 err = 6.0549094354e-08 time = 0.13 sec +[ Info: VUMPS 36: obj = -8.862878980877e-01 err = 4.2833729305e-08 time = 0.05 sec +[ Info: VUMPS 37: obj = -8.862878980877e-01 err = 3.0296142094e-08 time = 0.06 sec +[ Info: VUMPS 38: obj = -8.862878980877e-01 err = 2.1424846144e-08 time = 0.06 sec +[ Info: VUMPS 39: obj = -8.862878980877e-01 err = 1.5148957750e-08 time = 0.06 sec +[ Info: VUMPS 40: obj = -8.862878980877e-01 err = 1.0709953737e-08 time = 0.12 sec +[ Info: VUMPS 41: obj = -8.862878980877e-01 err = 7.5707891268e-09 time = 0.03 sec +[ Info: VUMPS 42: obj = -8.862878980877e-01 err = 5.3510570418e-09 time = 0.04 sec +[ Info: VUMPS 43: obj = -8.862878980877e-01 err = 3.7817943502e-09 time = 0.07 sec +[ Info: VUMPS 44: obj = -8.862878980877e-01 err = 2.6724193960e-09 time = 0.06 sec +[ Info: VUMPS 45: obj = -8.862878980877e-01 err = 1.8882939261e-09 time = 0.07 sec +[ Info: VUMPS 46: obj = -8.862878980877e-01 err = 1.3341314150e-09 time = 0.11 sec +[ Info: VUMPS 47: obj = -8.862878980877e-01 err = 9.4252990942e-10 time = 0.06 sec +[ Info: VUMPS 48: obj = -8.862878980877e-01 err = 6.6582278922e-10 time = 0.06 sec +[ Info: VUMPS 49: obj = -8.862878980877e-01 err = 4.7032251942e-10 time = 0.06 sec +[ Info: VUMPS 50: obj = -8.862878980877e-01 err = 3.3220639408e-10 time = 0.06 sec +[ Info: VUMPS 51: obj = -8.862878980878e-01 err = 2.3463502230e-10 time = 0.12 sec +[ Info: VUMPS 52: obj = -8.862878980878e-01 err = 1.6571486183e-10 time = 0.05 sec +[ Info: VUMPS 53: obj = -8.862878980878e-01 err = 1.1703442495e-10 time = 0.05 sec +[ Info: VUMPS 54: obj = -8.862878980878e-01 err = 8.2650109655e-11 time = 0.06 sec +[ Info: VUMPS 55: obj = -8.862878980878e-01 err = 5.8367474734e-11 time = 0.05 sec +[ Info: VUMPS 56: obj = -8.862878980878e-01 err = 4.1213417017e-11 time = 0.12 sec +[ Info: VUMPS 57: obj = -8.862878980878e-01 err = 2.9101697547e-11 time = 0.04 sec +[ Info: VUMPS 58: obj = -8.862878980878e-01 err = 2.0551179926e-11 time = 0.06 sec +[ Info: VUMPS 59: obj = -8.862878980878e-01 err = 1.4510549999e-11 time = 0.06 sec +[ Info: VUMPS 60: obj = -8.862878980878e-01 err = 1.0245548104e-11 time = 0.05 sec +[ Info: VUMPS 61: obj = -8.862878980878e-01 err = 7.2325909689e-12 time = 0.06 sec +[ Info: VUMPS 62: obj = -8.862878980878e-01 err = 5.1092588216e-12 time = 0.12 sec +[ Info: VUMPS 63: obj = -8.862878980878e-01 err = 3.6043616497e-12 time = 0.02 sec +[ Info: VUMPS 64: obj = -8.862878980878e-01 err = 2.5462748087e-12 time = 0.04 sec +[ Info: VUMPS 65: obj = -8.862878980878e-01 err = 1.7984804673e-12 time = 0.03 sec +[ Info: VUMPS 66: obj = -8.862878980878e-01 err = 1.2696913652e-12 time = 0.02 sec +[ Info: VUMPS conv 67: obj = -8.862878980879e-01 err = 8.9456922075e-13 time = 8.25 sec ```` diff --git a/docs/src/examples/quantum1d/4.xxz-heisenberg/main.ipynb b/docs/src/examples/quantum1d/4.xxz-heisenberg/main.ipynb index 9d851880a..39316fad3 100644 --- a/docs/src/examples/quantum1d/4.xxz-heisenberg/main.ipynb +++ b/docs/src/examples/quantum1d/4.xxz-heisenberg/main.ipynb @@ -24,7 +24,7 @@ "source": [ "## Failure\n", "\n", - "First we should define the hamiltonian we want to work with.\n", + "First we should define the Hamiltonian we want to work with.\n", "Then we specify an initial guess, which we then further optimize.\n", "Working directly in the thermodynamic limit, this is achieved as follows:" ], @@ -42,7 +42,7 @@ { "cell_type": "markdown", "source": [ - "We then need an intial state, which we shall later optimize. In this example we work directly in the thermodynamic limit." + "We then need an initial state, which we shall later optimize. In this example we work directly in the thermodynamic limit." ], "metadata": {} }, @@ -58,7 +58,7 @@ { "cell_type": "markdown", "source": [ - "The groundstate can then be found by calling `find_groundstate`." + "The ground state can then be found by calling `find_groundstate`." ], "metadata": {} }, @@ -75,7 +75,7 @@ "cell_type": "markdown", "source": [ "As you can see, VUMPS struggles to converge.\n", - "On it's own, that is already quite curious.\n", + "On its own, that is already quite curious.\n", "Maybe we can do better using another algorithm, such as gradient descent." ], "metadata": {} @@ -137,9 +137,9 @@ { "cell_type": "markdown", "source": [ - "In MPSKit, we require that the periodicity of the hamiltonian equals that of the state it is applied to.\n", - "This is not a big obstacle, you can simply repeat the original hamiltonian.\n", - "Alternatively, the hamiltonian can be constructed directly on a two-site unitcell by making use of MPSKitModels.jl's `@mpoham`." + "In MPSKit, we require that the periodicity of the Hamiltonian equals that of the state it is applied to.\n", + "This is not a big obstacle, you can simply repeat the original Hamiltonian.\n", + "Alternatively, the Hamiltonian can be constructed directly on a two-site unit cell by making use of MPSKitModels.jl's `@mpoham`." ], "metadata": {} }, @@ -193,9 +193,9 @@ "source": [ "## Symmetries\n", "\n", - "The XXZ Heisenberg hamiltonian is SU(2) symmetric and we can exploit this to greatly speed up the simulation.\n", + "The XXZ Heisenberg Hamiltonian is SU(2) symmetric and we can exploit this to greatly speed up the simulation.\n", "\n", - "It is cumbersome to construct symmetric hamiltonians, but luckily su(2) symmetric XXZ is already implemented:" + "It is cumbersome to construct symmetric Hamiltonians, but luckily SU(2) symmetric XXZ is already implemented:" ], "metadata": {} }, @@ -267,11 +267,11 @@ "file_extension": ".jl", "mimetype": "application/julia", "name": "julia", - "version": "1.12.0" + "version": "1.12.4" }, "kernelspec": { "name": "julia-1.12", - "display_name": "Julia 1.12.0", + "display_name": "Julia 1.12.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 cb5eeef18..903e39db1 100644 --- a/docs/src/examples/quantum1d/5.haldane-spt/index.md +++ b/docs/src/examples/quantum1d/5.haldane-spt/index.md @@ -54,10 +54,12 @@ H = heisenberg_hamiltonian() ```` ```` -single site InfiniteMPOHamiltonian{MPSKit.JordanMPOTensor{ComplexF64, TensorKit.GradedSpace{TensorKitSectors.SU2Irrep, TensorKit.SortedVectorDict{TensorKitSectors.SU2Irrep, Int64}}, Union{TensorKit.BraidingTensor{ComplexF64, TensorKit.GradedSpace{TensorKitSectors.SU2Irrep, TensorKit.SortedVectorDict{TensorKitSectors.SU2Irrep, Int64}}}, TensorKit.TensorMap{ComplexF64, TensorKit.GradedSpace{TensorKitSectors.SU2Irrep, TensorKit.SortedVectorDict{TensorKitSectors.SU2Irrep, Int64}}, 2, 2, Vector{ComplexF64}}}, TensorKit.TensorMap{ComplexF64, TensorKit.GradedSpace{TensorKitSectors.SU2Irrep, TensorKit.SortedVectorDict{TensorKitSectors.SU2Irrep, Int64}}, 2, 1, Vector{ComplexF64}}, TensorKit.TensorMap{ComplexF64, TensorKit.GradedSpace{TensorKitSectors.SU2Irrep, TensorKit.SortedVectorDict{TensorKitSectors.SU2Irrep, Int64}}, 1, 2, Vector{ComplexF64}}, TensorKit.TensorMap{ComplexF64, TensorKit.GradedSpace{TensorKitSectors.SU2Irrep, TensorKit.SortedVectorDict{TensorKitSectors.SU2Irrep, Int64}}, 1, 1, Vector{ComplexF64}}}}: -╷ ⋮ -┼ W[1]: 3×1×1×3 JordanMPOTensor(((Rep[SU₂](0=>1) ⊞ Rep[SU₂](1=>1) ⊞ Rep[SU₂](0=>1)) ⊗ ⊞(Rep[SU₂](1=>1))) ← (⊞(Rep[SU₂](1=>1)) ⊗ (Rep[SU₂](0=>1) ⊞ Rep[SU₂](1=>1) ⊞ Rep[SU₂](0=>1)))) -╵ ⋮ +1-site InfiniteMPOHamiltonian(ComplexF64, Rep[SU₂]) with maximal dimension 5: +| ⋮ +| ((0=>1) ⊞ (1=>1) ⊞ (0=>1)) +┼─[1]─ (1 => 1) +│ ((0=>1) ⊞ (1=>1) ⊞ (0=>1)) +| ⋮ ```` @@ -101,7 +103,7 @@ transferplot(ψ; sectors, title = "Transfer matrix spectrum", legend = :outertop ```` ```@raw html - + ``` Nevertheless, using the symmetry, this can be remedied rather easily, by imposing the @@ -128,7 +130,7 @@ E_plus = expectation_value(ψ_plus, H) ```` ```` --1.401419331339302 - 1.2835881733092454e-17im +-1.4014193313393002 - 1.1731597242433267e-16im ```` ````julia @@ -139,7 +141,7 @@ E_minus = expectation_value(ψ_minus, H) ```` ```` --1.401483973963085 + 2.188327525534786e-16im +-1.4014839739630833 + 2.1061363600641873e-16im ```` ````julia @@ -155,7 +157,7 @@ plot(transferp_plus, transferp_minus; layout = (1, 2), size = (800, 400)) ```` ```@raw html - + ``` ````julia @@ -165,7 +167,7 @@ plot(entanglementp_plus, entanglementp_minus; layout = (1, 2), size = (800, 400) ```` ```@raw html - + ``` As we can see, the ground state can be found in the non-trivial SPT phase, $\ket{\psi_-}$. We @@ -192,8 +194,8 @@ println("S_plus = $S_plus") ```` ```` -S_minus + log(2) = 1.5486227235438952 -S_plus = 1.54503235302521 +S_minus + log(2) = 1.5486227235412815 +S_plus = 1.545032353058148 ```` diff --git a/docs/src/examples/quantum1d/5.haldane-spt/main.ipynb b/docs/src/examples/quantum1d/5.haldane-spt/main.ipynb index 6c03ba77a..57e53f288 100644 --- a/docs/src/examples/quantum1d/5.haldane-spt/main.ipynb +++ b/docs/src/examples/quantum1d/5.haldane-spt/main.ipynb @@ -110,8 +110,8 @@ "cell_type": "markdown", "source": [ "Nevertheless, using the symmetry, this can be remedied rather easily, by imposing the\n", - "groundstate to belong to a single class, and comparing the results. We can readily obtain 3\n", - "different criteria for determining the SPT phase of the groundstate.\n", + "ground state to belong to a single class, and comparing the results. We can readily obtain 3\n", + "different criteria for determining the SPT phase of the ground state.\n", "\n", "Firstly, we can compare variational energies for states of similar bond dimensions. As we\n", "expect the state of the wrong SPT phase to have to expend some of its expressiveness in\n", @@ -158,19 +158,19 @@ { "cell_type": "markdown", "source": [ - "As we can see, the groundstate can be found in the non-trivial SPT phase, $\\ket{\\psi_-}$. We\n", + "As we can see, the ground state can be found in the non-trivial SPT phase, $\\ket{\\psi_-}$. We\n", "can obtain an intuitive understanding of $\\ket{\\psi_+}$ by considering the following\n", - "diagram. If we denote the MPS tensors that make up the groundstate as $A_-$, we can\n", - "construct a state in the trivial SPT phase that approximates the groundstate as follows:\n", + "diagram. If we denote the MPS tensors that make up the ground state as $A_-$, we can\n", + "construct a state in the trivial SPT phase that approximates the ground state as follows:\n", "\n", "\n", "\"SPT\n", "\n", "\n", "In other words, we can factorize a purely virtual isomorphism of $S = 1/2$ in order to\n", - "obtain the groundstate. This then also explains the degeneracies in the entanglement\n", + "obtain the ground state. This then also explains the degeneracies in the entanglement\n", "spectrum as well as in the transfer matrix spectrum. Finally, we can further confirm this\n", - "intuition by looking at the entanglement entropy of the groundstate. As we can see, the\n", + "intuition by looking at the entanglement entropy of the ground state. As we can see, the\n", "entanglement entropy of the state in the wrong SPT phase is exactly $log(2)$ higher than the\n", "one in the right SPT phase, which is exactly what we would expect from the diagram above." ], @@ -204,11 +204,11 @@ "file_extension": ".jl", "mimetype": "application/julia", "name": "julia", - "version": "1.12.0" + "version": "1.12.4" }, "kernelspec": { "name": "julia-1.12", - "display_name": "Julia 1.12.0", + "display_name": "Julia 1.12.4", "language": "julia" } }, diff --git a/docs/src/examples/quantum1d/6.hubbard/index.md b/docs/src/examples/quantum1d/6.hubbard/index.md index f032b7506..91c9c424b 100644 --- a/docs/src/examples/quantum1d/6.hubbard/index.md +++ b/docs/src/examples/quantum1d/6.hubbard/index.md @@ -45,6 +45,7 @@ using Plots using Interpolations using Optim + const t = 1.0 const mu = 0.0 const U = 3.0 @@ -114,54 +115,54 @@ Groundstate energy: ```` ```` -[ Info: VUMPS init: obj = -1.579081573143e+00 err = 5.2411e-01 -[ Info: VUMPS conv 7: obj = -4.377030027572e+00 err = 9.5143497631e-03 time = 8.04 sec -[ Info: VUMPS init: obj = -4.377030027572e+00 err = 1.6737e-02 -[ Info: VUMPS conv 6: obj = -4.378747216681e+00 err = 1.5036098687e-04 time = 0.37 sec -[ Info: VUMPS init: obj = -4.378747216681e+00 err = 8.0072e-03 -[ Info: VUMPS conv 6: obj = -4.379161050350e+00 err = 1.6595379794e-04 time = 0.38 sec -[ Info: VUMPS init: obj = -4.379161050350e+00 err = 6.1111e-03 -[ Info: VUMPS conv 5: obj = -4.379452139513e+00 err = 1.7635676107e-04 time = 0.40 sec -[ Info: VUMPS init: obj = -4.379452139513e+00 err = 5.6960e-03 -[ Info: VUMPS conv 4: obj = -4.379651657382e+00 err = 1.7362696368e-04 time = 0.38 sec -[ Info: VUMPS init: obj = -4.379651657382e+00 err = 4.1057e-03 -[ Info: VUMPS conv 3: obj = -4.379734893707e+00 err = 1.9418536588e-04 time = 0.32 sec -[ Info: VUMPS init: obj = -4.379734893707e+00 err = 3.5576e-03 -[ Info: VUMPS conv 3: obj = -4.379797734073e+00 err = 1.3530833964e-04 time = 0.50 sec -[ Info: VUMPS init: obj = -4.379797734073e+00 err = 2.7583e-03 -[ Info: VUMPS conv 2: obj = -4.379838443218e+00 err = 1.7548586831e-04 time = 0.34 sec -[ Info: VUMPS init: obj = -4.379838443218e+00 err = 2.7244e-03 -[ Info: VUMPS conv 3: obj = -4.379878819613e+00 err = 1.9910900463e-04 time = 0.68 sec -[ Info: VUMPS init: obj = -4.379878819613e+00 err = 2.6903e-03 -[ Info: VUMPS conv 3: obj = -4.379929209976e+00 err = 1.7353866489e-04 time = 0.89 sec -[ Info: VUMPS init: obj = -4.379929209976e+00 err = 2.5534e-03 -[ Info: VUMPS conv 3: obj = -4.379968033957e+00 err = 1.8030717590e-04 time = 1.15 sec -[ Info: VUMPS init: obj = -4.379968033957e+00 err = 1.7669e-03 -[ Info: VUMPS conv 2: obj = -4.379986875953e+00 err = 1.8529970888e-04 time = 0.75 sec -[ Info: VUMPS init: obj = -4.379986875953e+00 err = 1.5829e-03 -[ Info: VUMPS conv 2: obj = -4.380001233738e+00 err = 1.8857698446e-04 time = 0.86 sec -[ Info: VUMPS init: obj = -4.380001233739e+00 err = 1.4982e-03 -[ Info: VUMPS conv 2: obj = -4.380013190309e+00 err = 1.4999356908e-04 time = 1.14 sec -[ Info: VUMPS init: obj = -4.380013190309e+00 err = 1.4212e-03 -[ Info: VUMPS conv 2: obj = -4.380024391046e+00 err = 1.7666048459e-04 time = 1.45 sec -[ Info: VUMPS init: obj = -4.380024391046e+00 err = 1.3326e-03 -[ Info: VUMPS conv 2: obj = -4.380038154943e+00 err = 1.5589147881e-04 time = 1.71 sec -[ Info: VUMPS init: obj = -4.380038154943e+00 err = 1.0028e-03 -[ Info: VUMPS conv 1: obj = -4.380043680273e+00 err = 1.6697843620e-04 time = 0.74 sec -[ Info: VUMPS init: obj = -4.380043680277e+00 err = 9.0960e-04 -[ Info: VUMPS conv 1: obj = -4.380048637666e+00 err = 1.8527886431e-04 time = 0.83 sec -[ Info: VUMPS init: obj = -4.380048637678e+00 err = 8.3048e-04 -[ Info: VUMPS conv 1: obj = -4.380053195003e+00 err = 1.8076356558e-04 time = 0.94 sec -[ Info: VUMPS init: obj = -4.380053195036e+00 err = 6.8162e-04 -[ Info: VUMPS conv 1: obj = -4.380057137143e+00 err = 1.8928785082e-04 time = 1.17 sec -[ Info: VUMPS init: obj = -4.380057137148e+00 err = 6.0309e-04 -[ Info: VUMPS conv 1: obj = -4.380060540455e+00 err = 1.8334294673e-04 time = 1.33 sec -[ Info: VUMPS init: obj = -4.379609563543e+00 err = 4.0950e-03 -[ Info: VUMPS conv 19: obj = -4.379763089580e+00 err = 9.9477192904e-06 time = 8.38 sec -[ Info: CG: initializing with f = -4.379763089225, ‖∇f‖ = 3.1536e-05 -[ Info: CG: converged after 170 iterations and time 89.57 s: f = -4.379763311396, ‖∇f‖ = 9.4220e-07 +[ Info: VUMPS init: obj = -1.472548758147e+00 err = 5.1344e-01 +[ Info: VUMPS conv 6: obj = -4.376626743591e+00 err = 9.9886660196e-03 time = 4.42 sec +[ Info: VUMPS init: obj = -4.376626743590e+00 err = 2.1408e-02 +[ Info: VUMPS conv 8: obj = -4.378746003112e+00 err = 1.5955449521e-04 time = 0.42 sec +[ Info: VUMPS init: obj = -4.378746003112e+00 err = 8.1588e-03 +[ Info: VUMPS conv 6: obj = -4.379161057076e+00 err = 1.5466178028e-04 time = 0.43 sec +[ Info: VUMPS init: obj = -4.379161057076e+00 err = 6.1185e-03 +[ Info: VUMPS conv 5: obj = -4.379452137265e+00 err = 1.6793132081e-04 time = 0.42 sec +[ Info: VUMPS init: obj = -4.379452137265e+00 err = 5.6949e-03 +[ Info: VUMPS conv 4: obj = -4.379651664636e+00 err = 1.6997715633e-04 time = 0.51 sec +[ Info: VUMPS init: obj = -4.379651664636e+00 err = 4.1040e-03 +[ Info: VUMPS conv 3: obj = -4.379734890151e+00 err = 1.9458349625e-04 time = 0.43 sec +[ Info: VUMPS init: obj = -4.379734890151e+00 err = 3.5568e-03 +[ Info: VUMPS conv 3: obj = -4.379797732181e+00 err = 1.3874790123e-04 time = 0.49 sec +[ Info: VUMPS init: obj = -4.379797732181e+00 err = 2.7580e-03 +[ Info: VUMPS conv 2: obj = -4.379838437704e+00 err = 1.7903427058e-04 time = 0.37 sec +[ Info: VUMPS init: obj = -4.379838437704e+00 err = 2.7247e-03 +[ Info: VUMPS conv 3: obj = -4.379878817727e+00 err = 1.9890160987e-04 time = 0.86 sec +[ Info: VUMPS init: obj = -4.379878817727e+00 err = 2.6892e-03 +[ Info: VUMPS conv 3: obj = -4.379929192662e+00 err = 1.7735649634e-04 time = 1.10 sec +[ Info: VUMPS init: obj = -4.379929192662e+00 err = 2.5531e-03 +[ Info: VUMPS conv 3: obj = -4.379968013458e+00 err = 1.8131482590e-04 time = 2.36 sec +[ Info: VUMPS init: obj = -4.379968013458e+00 err = 1.7666e-03 +[ Info: VUMPS conv 2: obj = -4.379986848440e+00 err = 1.8697865996e-04 time = 1.01 sec +[ Info: VUMPS init: obj = -4.379986848440e+00 err = 1.5825e-03 +[ Info: VUMPS conv 2: obj = -4.380000756174e+00 err = 1.9154925356e-04 time = 1.54 sec +[ Info: VUMPS init: obj = -4.380000756174e+00 err = 1.5051e-03 +[ Info: VUMPS conv 2: obj = -4.380013136300e+00 err = 1.5515231683e-04 time = 1.88 sec +[ Info: VUMPS init: obj = -4.380013136300e+00 err = 1.4208e-03 +[ Info: VUMPS conv 2: obj = -4.380024381179e+00 err = 1.7768123800e-04 time = 2.03 sec +[ Info: VUMPS init: obj = -4.380024381179e+00 err = 1.3308e-03 +[ Info: VUMPS conv 2: obj = -4.380038135198e+00 err = 1.5610887648e-04 time = 2.47 sec +[ Info: VUMPS init: obj = -4.380038135198e+00 err = 1.0025e-03 +[ Info: VUMPS conv 1: obj = -4.380043663927e+00 err = 1.6686329927e-04 time = 1.29 sec +[ Info: VUMPS init: obj = -4.380043663928e+00 err = 9.0871e-04 +[ Info: VUMPS conv 1: obj = -4.380048620164e+00 err = 1.8550404960e-04 time = 1.73 sec +[ Info: VUMPS init: obj = -4.380048620164e+00 err = 8.3016e-04 +[ Info: VUMPS conv 1: obj = -4.380053179097e+00 err = 1.8078418705e-04 time = 1.88 sec +[ Info: VUMPS init: obj = -4.380053179097e+00 err = 6.8103e-04 +[ Info: VUMPS conv 1: obj = -4.380057119334e+00 err = 1.9088736457e-04 time = 3.48 sec +[ Info: VUMPS init: obj = -4.380057119334e+00 err = 6.0308e-04 +[ Info: VUMPS conv 1: obj = -4.380060522532e+00 err = 1.8437058341e-04 time = 3.48 sec +[ Info: VUMPS init: obj = -4.379610294430e+00 err = 4.0923e-03 +[ Info: VUMPS conv 19: obj = -4.379763567395e+00 err = 9.9309403645e-06 time = 11.96 sec +[ Info: CG: initializing with f = -4.379763567040e+00, ‖∇f‖ = 3.1464e-05 +[ Info: CG: converged after 154 iterations and time 1.64 m: f = -4.379763577686e+00, ‖∇f‖ = 9.2421e-07 ┌ Info: Groundstate energy: -│ * numerical: -2.1899960609614597 +│ * numerical: -2.1899960610249782 └ * analytic: -2.190038374277775 ```` @@ -195,56 +196,56 @@ Groundstate energy: ```` ```` -[ Info: VUMPS init: obj = +8.497854905331e-01 err = 6.4635e-01 -[ Info: VUMPS conv 1: obj = -4.000000000000e+00 err = 2.0411791683e-15 time = 4.67 sec +[ Info: VUMPS init: obj = -5.588042682004e-02 err = 9.2053e-01 +[ Info: VUMPS conv 1: obj = -4.000000000000e+00 err = 2.1240270466e-15 time = 2.96 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.06 sec -[ Info: VUMPS init: obj = -4.289650419703e+00 err = 1.1203e-01 -[ Info: VUMPS conv 6: obj = -4.359865567633e+00 err = 1.0046406773e-04 time = 0.16 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 = 2.56 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.09 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.29 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.18 sec -[ Info: VUMPS init: obj = -4.379272965037e+00 err = 9.9127e-03 -[ Info: VUMPS conv 4: obj = -4.379592225819e+00 err = 1.5571559713e-04 time = 0.19 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.27 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 = 0.88 sec -└ @ MPSKit ~/Projects/MPSKit.jl/src/algorithms/groundstate/vumps.jl:76 -[ 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.37 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.53 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.42 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.63 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.35 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.49 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.70 sec -[ Info: VUMPS init: obj = -4.380072130335e+00 err = 1.3827e-04 -[ Info: VUMPS conv 1: obj = -4.380073468947e+00 err = 8.7912179067e-05 time = 1.13 sec -[ Info: VUMPS init: obj = -4.380073469160e+00 err = 9.9636e-05 -[ Info: VUMPS conv 1: obj = -4.380074456012e+00 err = 6.8484929131e-05 time = 3.50 sec -[ Info: VUMPS init: obj = -4.380074456584e+00 err = 7.7488e-05 -[ Info: VUMPS conv 1: obj = -4.380075158440e+00 err = 6.5951776337e-05 time = 3.49 sec -[ Info: VUMPS init: obj = -4.380075159046e+00 err = 6.3891e-05 -[ Info: VUMPS conv 1: obj = -4.380075656659e+00 err = 4.7588615231e-05 time = 8.17 sec -[ Info: VUMPS init: obj = -4.379308902467e+00 err = 8.0000e-03 -┌ Warning: VUMPS cancel 100: obj = -4.379693243679e+00 err = 1.5955101978e-05 time = 15.41 sec -└ @ MPSKit ~/Projects/MPSKit.jl/src/algorithms/groundstate/vumps.jl:76 -[ Info: CG: initializing with f = -4.379693243679, ‖∇f‖ = 5.7833e-05 -[ Info: CG: converged after 13 iterations and time 6.25 s: f = -4.379693244608, ‖∇f‖ = 6.3799e-07 +[ Info: VUMPS conv 4: obj = -4.289650419749e+00 err = 1.8514003381e-04 time = 0.05 sec +[ Info: VUMPS init: obj = -4.289650419749e+00 err = 1.1203e-01 +[ Info: VUMPS conv 6: obj = -4.359865567620e+00 err = 1.0046942911e-04 time = 0.14 sec +[ Info: VUMPS init: obj = -4.359865567619e+00 err = 4.3643e-02 +[ Info: VUMPS conv 6: obj = -4.372880928482e+00 err = 1.3025843115e-04 time = 2.31 sec +[ Info: VUMPS init: obj = -4.372880928482e+00 err = 3.2693e-02 +[ Info: VUMPS conv 4: obj = -4.375236954488e+00 err = 1.1814239608e-04 time = 0.17 sec +[ Info: VUMPS init: obj = -4.375236954488e+00 err = 2.9487e-02 +[ Info: VUMPS conv 7: obj = -4.378159084364e+00 err = 1.1896740056e-04 time = 0.51 sec +[ Info: VUMPS init: obj = -4.378159084364e+00 err = 1.9312e-02 +[ Info: VUMPS conv 5: obj = -4.379272966040e+00 err = 1.5785413165e-04 time = 0.34 sec +[ Info: VUMPS init: obj = -4.379272966040e+00 err = 9.9128e-03 +[ Info: VUMPS conv 4: obj = -4.379592229143e+00 err = 1.5550378745e-04 time = 0.21 sec +[ Info: VUMPS init: obj = -4.379592229143e+00 err = 6.4841e-03 +[ Info: VUMPS conv 4: obj = -4.379819377264e+00 err = 1.7492038571e-04 time = 0.27 sec +[ Info: VUMPS init: obj = -4.379819377264e+00 err = 3.8754e-03 +┌ Warning: VUMPS cancel 10: obj = -4.379964033305e+00 err = 2.1228930049e-04 time = 1.29 sec +└ @ MPSKit /home/ldevos/LocalProjects/MPSKit.jl/src/algorithms/groundstate/vumps.jl:76 +[ Info: VUMPS init: obj = -4.379964033305e+00 err = 2.8978e-03 +[ Info: VUMPS conv 3: obj = -4.380010384710e+00 err = 1.4775284542e-04 time = 0.52 sec +[ Info: VUMPS init: obj = -4.380010384710e+00 err = 2.0609e-03 +[ Info: VUMPS conv 3: obj = -4.380041751503e+00 err = 1.6327798118e-04 time = 0.86 sec +[ Info: VUMPS init: obj = -4.380041751502e+00 err = 1.2364e-03 +[ Info: VUMPS conv 2: obj = -4.380055778759e+00 err = 1.8366845284e-04 time = 0.53 sec +[ Info: VUMPS init: obj = -4.380055778759e+00 err = 8.5857e-04 +[ Info: VUMPS conv 2: obj = -4.380064749427e+00 err = 1.3905442267e-04 time = 0.96 sec +[ Info: VUMPS init: obj = -4.380064749427e+00 err = 5.2502e-04 +[ Info: VUMPS conv 1: obj = -4.380067974777e+00 err = 1.5646700070e-04 time = 0.51 sec +[ Info: VUMPS init: obj = -4.380067974777e+00 err = 3.3275e-04 +[ Info: VUMPS conv 1: obj = -4.380070351418e+00 err = 1.3123916502e-04 time = 1.05 sec +[ Info: VUMPS init: obj = -4.380070351418e+00 err = 2.0348e-04 +[ Info: VUMPS conv 1: obj = -4.380072125256e+00 err = 1.1119707628e-04 time = 1.50 sec +[ Info: VUMPS init: obj = -4.380072125256e+00 err = 1.3635e-04 +[ Info: VUMPS conv 1: obj = -4.380073467830e+00 err = 8.5045032312e-05 time = 2.38 sec +[ Info: VUMPS init: obj = -4.380073467830e+00 err = 9.7226e-05 +[ Info: VUMPS conv 1: obj = -4.380074455763e+00 err = 6.4430026630e-05 time = 3.97 sec +[ Info: VUMPS init: obj = -4.380074455763e+00 err = 7.3787e-05 +[ Info: VUMPS conv 1: obj = -4.380075159887e+00 err = 6.2144398837e-05 time = 6.74 sec +[ Info: VUMPS init: obj = -4.380075159887e+00 err = 5.9899e-05 +[ Info: VUMPS conv 1: obj = -4.380075661721e+00 err = 4.2515939995e-05 time = 11.38 sec +[ Info: VUMPS init: obj = -4.379308795201e+00 err = 7.9930e-03 +┌ Warning: VUMPS cancel 100: obj = -4.379692711472e+00 err = 1.5979764630e-05 time = 17.26 sec +└ @ MPSKit /home/ldevos/LocalProjects/MPSKit.jl/src/algorithms/groundstate/vumps.jl:76 +[ Info: CG: initializing with f = -4.379692711471e+00, ‖∇f‖ = 5.7923e-05 +[ Info: CG: converged after 13 iterations and time 5.38 s: f = -4.379692712392e+00, ‖∇f‖ = 6.2087e-07 ┌ Info: Groundstate energy: -│ * numerical: -2.1900153475149375 +│ * numerical: -2.190015347514449 └ * analytic: -2.190038374277775 ```` @@ -281,71 +282,71 @@ E_holon, ϕ_holon = excitations( ```` ```` +[ Info: Found excitations for momentum = 2.945243112740431 [ 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.356194490192345 -[ Info: Found excitations for momentum = -2.1598449493429825 -[ 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 -[ Info: Found excitations for momentum = -0.5890486225480862 -[ 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.19634954084936207 +[ Info: Found excitations for momentum = -2.552544031041707 +[ Info: Found excitations for momentum = -2.748893571891069 +[ Info: Found excitations for momentum = -2.945243112740431 +[ Info: Found excitations for momentum = -0.39269908169872414 [ Info: Found excitations for momentum = 0.39269908169872414 +[ Info: Found excitations for momentum = 2.552544031041707 +[ Info: Found excitations for momentum = 2.748893571891069 +[ Info: Found excitations for momentum = -0.5890486225480862 [ Info: Found excitations for momentum = 0.5890486225480862 +[ Info: Found excitations for momentum = -0.7853981633974483 [ Info: Found excitations for momentum = 0.7853981633974483 +[ Info: Found excitations for momentum = 2.356194490192345 +[ Info: Found excitations for momentum = -2.356194490192345 +[ Info: Found excitations for momentum = 1.5707963267948966 +[ Info: Found excitations for momentum = -1.5707963267948966 [ Info: Found excitations for momentum = 0.9817477042468103 -[ Info: Found excitations for momentum = 1.1780972450961724 +[ Info: Found excitations for momentum = 2.1598449493429825 +[ Info: Found excitations for momentum = -0.9817477042468103 +[ Info: Found excitations for momentum = -2.1598449493429825 +[ Info: Found excitations for momentum = -1.1780972450961724 +[ Info: Found excitations for momentum = 1.7671458676442586 +[ Info: Found excitations for momentum = 1.9634954084936207 [ Info: Found excitations for momentum = 1.3744467859455345 +[ Info: Found excitations for momentum = -1.7671458676442586 +[ Info: Found excitations for momentum = -1.3744467859455345 +[ Info: Found excitations for momentum = -1.9634954084936207 +[ Info: Found excitations for momentum = 3.141592653589793 +[ Info: Found excitations for momentum = 1.1780972450961724 +[ 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.3744467859455345 [ Info: Found excitations for momentum = 1.7671458676442586 [ Info: Found excitations for momentum = 1.9634954084936207 +[ Info: Found excitations for momentum = -1.5707963267948966 +[ Info: Found excitations for momentum = -1.9634954084936207 +[ Info: Found excitations for momentum = 1.1780972450961724 [ 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.748893571891069 -[ Info: Found excitations for momentum = 2.945243112740431 -[ Info: Found excitations for momentum = 3.141592653589793 -[ 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.356194490192345 +[ Info: Found excitations for momentum = 0.0 [ Info: Found excitations for momentum = -2.1598449493429825 -[ 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 = -3.141592653589793 +[ Info: Found excitations for momentum = -2.356194490192345 +[ Info: Found excitations for momentum = 2.552544031041707 +[ Info: Found excitations for momentum = 0.7853981633974483 +[ Info: Found excitations for momentum = 2.356194490192345 [ 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 = 2.748893571891069 +[ Info: Found excitations for momentum = -2.748893571891069 +[ Info: Found excitations for momentum = 0.5890486225480862 +[ Info: Found excitations for momentum = 0.39269908169872414 +[ Info: Found excitations for momentum = -2.552544031041707 [ 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 = 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 -[ Info: Found excitations for momentum = 1.9634954084936207 -[ 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.748893571891069 [ Info: Found excitations for momentum = 2.945243112740431 +[ Info: Found excitations for momentum = 0.19634954084936207 +[ Info: Found excitations for momentum = -2.945243112740431 [ Info: Found excitations for momentum = 3.141592653589793 ```` @@ -400,7 +401,7 @@ end ```` ```@raw html - + ``` The plot shows some discrepancies between the numerical and analytic results. @@ -422,7 +423,7 @@ end ```` ```@raw html - + ``` The second discrepancy is that while the spinon dispersion is well-reproduced, the holon dispersion is not. @@ -490,7 +491,7 @@ end ```` ```@raw html - + ``` --- diff --git a/docs/src/examples/quantum1d/6.hubbard/main.ipynb b/docs/src/examples/quantum1d/6.hubbard/main.ipynb index 932573445..1524ffac0 100644 --- a/docs/src/examples/quantum1d/6.hubbard/main.ipynb +++ b/docs/src/examples/quantum1d/6.hubbard/main.ipynb @@ -52,6 +52,7 @@ "using Interpolations\n", "using Optim\n", "\n", + "\n", "const t = 1.0\n", "const mu = 0.0\n", "const U = 3.0" @@ -62,7 +63,7 @@ { "cell_type": "markdown", "source": [ - "For this case, the groundstate energy has an analytic solution, which can be used to benchmark the numerical results.\n", + "For this case, the ground state energy has an analytic solution, which can be used to benchmark the numerical results.\n", "It follows from Eq. (6.82) in []().\n", "\n", "$$\n", @@ -384,11 +385,11 @@ "file_extension": ".jl", "mimetype": "application/julia", "name": "julia", - "version": "1.12.0" + "version": "1.12.4" }, "kernelspec": { "name": "julia-1.12", - "display_name": "Julia 1.12.0", + "display_name": "Julia 1.12.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 4b01212f7..983ad0858 100644 --- a/docs/src/examples/quantum1d/7.xy-finiteT/index.md +++ b/docs/src/examples/quantum1d/7.xy-finiteT/index.md @@ -125,11 +125,11 @@ println("Exact (N=Inf):\t", groundstate_energy(J, Inf)) ```` ```` -[ Info: DMRG2 1: obj = -5.004084850990e+00 err = 9.7717634223e-01 time = 1.69 min -[ Info: DMRG2 2: obj = -5.004096939452e+00 err = 1.1654525568e-06 time = 0.49 sec -[ Info: DMRG2 3: obj = -5.004096975044e+00 err = 2.4907168550e-09 time = 0.61 sec -[ Info: DMRG2 conv 4: obj = -5.004096975044e+00 err = 1.1812772982e-13 time = 1.71 min -Numerical: -0.15637803047010954 +[ Info: DMRG2 1: obj = -5.004084801485e+00 err = 9.7485774328e-01 time = 1.43 min +[ Info: DMRG2 2: obj = -5.004096940647e+00 err = 1.1899230994e-06 time = 1.27 sec +[ Info: DMRG2 3: obj = -5.004096975044e+00 err = 2.2262868216e-09 time = 0.80 sec +[ Info: DMRG2 conv 4: obj = -5.004096975044e+00 err = 1.1612932838e-13 time = 1.47 min +Numerical: -0.15637803047010942 Exact (N=32): -0.15637803047254015 Exact (N=Inf): -0.15915494309189535 diff --git a/docs/src/examples/quantum1d/7.xy-finiteT/main.ipynb b/docs/src/examples/quantum1d/7.xy-finiteT/main.ipynb index 8936d8bde..c6a7df3ce 100644 --- a/docs/src/examples/quantum1d/7.xy-finiteT/main.ipynb +++ b/docs/src/examples/quantum1d/7.xy-finiteT/main.ipynb @@ -696,11 +696,11 @@ "file_extension": ".jl", "mimetype": "application/julia", "name": "julia", - "version": "1.12.0" + "version": "1.12.4" }, "kernelspec": { "name": "julia-1.12", - "display_name": "Julia 1.12.0", + "display_name": "Julia 1.12.4", "language": "julia" } }, diff --git a/docs/src/examples/quantum1d/8.bose-hubbard/index.md b/docs/src/examples/quantum1d/8.bose-hubbard/index.md index 476461d7a..7c1ff7ee0 100644 --- a/docs/src/examples/quantum1d/8.bose-hubbard/index.md +++ b/docs/src/examples/quantum1d/8.bose-hubbard/index.md @@ -11,6 +11,7 @@ using Markdown using MPSKit, MPSKitModels, TensorKit using Plots, LaTeXStrings + theme(:wong) default(fontfamily = "Computer Modern", label = nothing, dpi = 100, framestyle = :box) ```` @@ -153,9 +154,9 @@ println("Energy: ", expectation_value(ground_state, hamiltonian)) ```` ```` -[ Info: VUMPS init: obj = +5.102229926686e-01 err = 6.1730e-01 -[ Info: VUMPS conv 44: obj = -6.757777651150e-01 err = 9.9483095620e-07 time = 13.31 sec -Energy: -0.6757777651149829 + 1.3734202787398213e-17im +[ Info: VUMPS init: obj = +4.962471958690e-01 err = 5.8876e-01 +[ Info: VUMPS conv 50: obj = -6.757777651150e-01 err = 9.9482617748e-07 time = 6.72 sec +Energy: -0.6757777651149999 - 6.246015340721835e-17im ```` @@ -197,7 +198,7 @@ hline!([abs2(expectation_value(ground_state, (0,) => a_op))], ls = :dash, c = :b ```` ```@raw html - + ``` We see that the correlations drop off exponentially, indicating the existence of a gapped @@ -211,7 +212,7 @@ hline!([abs2(expectation_value(ground_state, (0,) => a_op))], ls = :dash, c = :b ```` ```@raw html - + ``` In this case, the correlation function drops off algebraically and eventually saturates as @@ -275,7 +276,7 @@ scatter!( ```` ```@raw html - + ``` This shows that any finite bond dimension MPS necessarily breaks the symmetry of the system, @@ -290,13 +291,13 @@ quasicondensate_density = map(state -> abs2(expectation_value(state, (0,) => a_o ```` 7-element Vector{Float64}: - 0.30974277868294353 - 0.2881478518741761 - 0.2702001385414079 - 0.2571272814147897 - 0.24685385675266389 - 0.23539778700973993 - 0.2279965507292309 + 0.30974277207656425 + 0.28814775930068737 + 0.2702000951730164 + 0.25712728156142256 + 0.2468538617129866 + 0.2353979140629328 + 0.2279965552408022 ```` We may now also visualize the momentum distribution function, which is obtained as the @@ -350,7 +351,7 @@ plot(ks, momentum_distribution, lab = "D = " .* string.(permutedims(Ds)), lw = 1 ```` ```@raw html - + ``` We see that the density seems to peak around $k=0$, this time seemingly becoming more @@ -434,7 +435,7 @@ superfluid_stiffness_profile(0.01, 0.3, 5, 4) # mott insulator ```` ```@raw html - + ``` Now that we know what phases to expect, we can plot the phase diagram by scanning over a @@ -465,7 +466,7 @@ heatmap(ts, mus, order_parameters, xlabel = L"t/U", ylabel = L"\mu/U", title = L ```` ```@raw html - + ``` Although the bond dimension here is quite low, we already see the deformation of the Mott diff --git a/docs/src/examples/quantum1d/8.bose-hubbard/main.ipynb b/docs/src/examples/quantum1d/8.bose-hubbard/main.ipynb index 83681551a..d4ace3b04 100644 --- a/docs/src/examples/quantum1d/8.bose-hubbard/main.ipynb +++ b/docs/src/examples/quantum1d/8.bose-hubbard/main.ipynb @@ -8,6 +8,7 @@ "using MPSKit, MPSKitModels, TensorKit\n", "using Plots, LaTeXStrings\n", "\n", + "\n", "theme(:wong)\n", "default(fontfamily = \"Computer Modern\", label = nothing, dpi = 100, framestyle = :box)" ], @@ -542,11 +543,11 @@ "file_extension": ".jl", "mimetype": "application/julia", "name": "julia", - "version": "1.11.6" + "version": "1.12.4" }, "kernelspec": { - "name": "julia-1.11", - "display_name": "Julia 1.11.6", + "name": "julia-1.12", + "display_name": "Julia 1.12.4", "language": "julia" } }, diff --git a/examples/Cache.toml b/examples/Cache.toml index 0c2fdb60b..b41aab9b5 100644 --- a/examples/Cache.toml +++ b/examples/Cache.toml @@ -1,12 +1,12 @@ [classic2d] -"1.hard-hexagon" = "2c01966230ba56faa19a6751a734e53b8db744e7522b03b6a554662ccc76c550" +"1.hard-hexagon" = "b78b475bee20d0366789a01277b6242089c9a610290df012bd66ad880127553b" [quantum1d] -"2.haldane" = "804433b690faa1ce268a430edb4185af77a38de7a9f53e097795688a53c30c5e" -"6.hubbard" = "af14e1df11b2392a69260ce061a716370a2ce50b5c907f0a7d2548e796d543fe" -"7.xy-finiteT" = "7afb26fb9ff8ef84722fee3442eaed2ddffda4a5f70a7844b61ce5178093706c" -"8.bose-hubbard" = "4ca414d4b76bbffbd88a7c153dce245c9fab5b7cee0fe0b706b9a10ed87e29e9" -"3.ising-dqpt" = "2d314fd05a75c5c91ff81391c7bf166171b35a7b32933e9490756dc584fe8437" -"5.haldane-spt" = "3de1b3baa1e4c5dc2252bbe8688d0c6ac8e5d5265deeb6ab66818bbfb02dd4aa" -"4.xxz-heisenberg" = "1a2093a182f3ce070b53488484d1024e45496b291c1b74e3f370201d4c056be2" -"1.ising-cft" = "fc99e02aed275ca8f588a780a9de09ab09e16901681fcb72103e1e7dc4be157c" +"2.haldane" = "b81cbd6a4b9347be45acb779fa6e91b3c64c127d39e8a299aea64198ca24fba7" +"6.hubbard" = "f5710bddbb5bc7dc95b69978d45d70fe0a3d3d65dfbad66ad985d8977914c710" +"7.xy-finiteT" = "0f330a157bea739a43a82a937791c680b2fa6e5e479171ee5ef318d1fdae7bcf" +"8.bose-hubbard" = "a060ac3e29a9169e9420f7125cdb020c83b26a380508ac85f108bff1562bc70c" +"3.ising-dqpt" = "735882054a689df10ce0123c56fcc131d5d1156b8d12e9237d544e81b086250d" +"5.haldane-spt" = "9fce2c2d34c9e3a3b91dff3c3cba6404db6880fd92e66dcdfbbe7a6f3067495b" +"4.xxz-heisenberg" = "9033fb104c3f658ccb1b8e335b986b7abda59e0b5359be086d3b62e1cc32ebd6" +"1.ising-cft" = "ee6512b8e92aaaf42142b80024b905e9396aff75d4d60f77a29c18b234d1f65c" diff --git a/examples/classic2d/1.hard-hexagon/main.jl b/examples/classic2d/1.hard-hexagon/main.jl index 600703c4a..bdd6dc41a 100644 --- a/examples/classic2d/1.hard-hexagon/main.jl +++ b/examples/classic2d/1.hard-hexagon/main.jl @@ -12,6 +12,10 @@ To follow the tutorial you need the following packages. using MPSKit, MPSKitModels, TensorKit, Plots, Polynomials +#src # for reproducibility: +#src using Random +#src Random.seed!(123) + md""" The [hard hexagon model](https://en.wikipedia.org/wiki/Hard_hexagon_model) is a 2-dimensional lattice model of a gas, where particles are allowed to be on the vertices of a triangular lattice, but no two particles may be adjacent. This can be encoded in a transfer matrix with a local MPO tensor using anyonic symmetries, and the resulting MPO has been implemented in MPSKitModels. diff --git a/examples/quantum1d/4.xxz-heisenberg/main.jl b/examples/quantum1d/4.xxz-heisenberg/main.jl index 52930bdc0..8c684e911 100644 --- a/examples/quantum1d/4.xxz-heisenberg/main.jl +++ b/examples/quantum1d/4.xxz-heisenberg/main.jl @@ -7,6 +7,10 @@ The necessary packages to follow this tutorial are: using MPSKit, MPSKitModels, TensorKit, Plots +#src # for reproducibility: +#src using Random +#src Random.seed!(123) + md""" ## Failure diff --git a/examples/quantum1d/6.hubbard/main.jl b/examples/quantum1d/6.hubbard/main.jl index 50eabfe8b..155816c85 100644 --- a/examples/quantum1d/6.hubbard/main.jl +++ b/examples/quantum1d/6.hubbard/main.jl @@ -36,6 +36,10 @@ using Plots using Interpolations using Optim +#src # for reproducibility: +#src using Random +#src Random.seed!(123) + const t = 1.0 const mu = 0.0 const U = 3.0 diff --git a/examples/quantum1d/7.xy-finiteT/main.jl b/examples/quantum1d/7.xy-finiteT/main.jl index fd154438f..03ad41999 100644 --- a/examples/quantum1d/7.xy-finiteT/main.jl +++ b/examples/quantum1d/7.xy-finiteT/main.jl @@ -9,6 +9,10 @@ using Plots using LinearAlgebra using BenchmarkFreeFermions +#src # for reproducibility: +#src using Random +#src Random.seed!(123) + md""" # Finite temperature XY model diff --git a/examples/quantum1d/8.bose-hubbard/main.jl b/examples/quantum1d/8.bose-hubbard/main.jl index 293a23822..571107698 100644 --- a/examples/quantum1d/8.bose-hubbard/main.jl +++ b/examples/quantum1d/8.bose-hubbard/main.jl @@ -2,6 +2,10 @@ using Markdown using MPSKit, MPSKitModels, TensorKit using Plots, LaTeXStrings +#src # for reproducibility: +#src using Random +#src Random.seed!(123) + theme(:wong) default(fontfamily = "Computer Modern", label = nothing, dpi = 100, framestyle = :box)