Update XXZ and Haldane examples

Author: lkdvos

docs/src/examples/quantum1d/4.xxz-heisenberg/index.md

@@ -34,7 +34,7 @@
    "outputs": [],
    "cell_type": "code",
    "source": [
-    "H = heisenberg_XXX(; spin=1 // 2);"
+    "H = heisenberg_XXX(; spin=1 // 2)"
    ],
    "metadata": {},
    "execution_count": null
@@ -50,8 +50,7 @@
    "outputs": [],
    "cell_type": "code",
    "source": [
-    "random_data = TensorMap(rand, ComplexF64, ℂ^20 * ℂ^2, ℂ^20);\n",
-    "state = InfiniteMPS([random_data]);"
+    "state = InfiniteMPS(2, 20)"
    ],
    "metadata": {},
    "execution_count": null
@@ -130,9 +129,7 @@
    "outputs": [],
    "cell_type": "code",
    "source": [
-    "A = TensorMap(rand, ComplexF64, ℂ^20 * ℂ^2, ℂ^20);\n",
-    "B = TensorMap(rand, ComplexF64, ℂ^20 * ℂ^2, ℂ^20);\n",
-    "state = InfiniteMPS([A, B]);"
+    "state = InfiniteMPS(fill(2, 2), fill(20, 2))"
    ],
    "metadata": {},
    "execution_count": null
@@ -152,7 +149,7 @@
    "source": [
     "# H2 = repeat(H, 2); -- copies the one-site version\n",
     "H2 = heisenberg_XXX(ComplexF64, Trivial, InfiniteChain(2); spin=1 // 2)\n",
-    "groundstate, cache, delta = find_groundstate(state, H2,\n",
+    "groundstate, envs, delta = find_groundstate(state, H2,\n",
     "                                             VUMPS(; maxiter=100, tol=1e-12));"
    ],
    "metadata": {},
@@ -170,9 +167,9 @@
    "outputs": [],
    "cell_type": "code",
    "source": [
-    "groundstate, cache, delta = find_groundstate(state, H2,\n",
+    "groundstate, envs, delta = find_groundstate(state, H2,\n",
     "                                             IDMRG2(; trscheme=truncdim(50), maxiter=20,\n",
-    "                                                    tol=1e-12))\n",
+    "                                                    tol=1e-12));\n",
     "entanglementplot(groundstate)"
    ],
    "metadata": {},
@@ -270,11 +267,11 @@
    "file_extension": ".jl",
    "mimetype": "application/julia",
    "name": "julia",
-   "version": "1.11.1"
+   "version": "1.11.2"
   },
   "kernelspec": {
    "name": "julia-1.11",
-   "display_name": "Julia 1.11.1",
+   "display_name": "Julia 1.11.2",
    "language": "julia"
   }
  },

docs/src/examples/quantum1d/4.xxz-heisenberg/main.ipynb

@@ -46,13 +46,13 @@
     "function heisenberg_hamiltonian(; J=-1.0)\n",
     "    s = SU2Irrep(1)\n",
     "    ℋ = SU2Space(1 => 1)\n",
-    "    SS = TensorMap(zeros, ComplexF64, ℋ ⊗ ℋ ← ℋ ⊗ ℋ)\n",
+    "    SS = zeros(ComplexF64, ℋ ⊗ ℋ ← ℋ ⊗ ℋ)\n",
     "    for (S, data) in blocks(SS)\n",
     "        data .= -0.5J * (casimir(S) - casimir(s) - casimir(s))\n",
     "    end\n",
-    "    return MPOHamiltonian(SS)\n",
+    "    return InfiniteMPOHamiltonian(SS)\n",
     "end\n",
-    "H = heisenberg_hamiltonian();"
+    "H = heisenberg_hamiltonian()"
    ],
    "metadata": {},
    "execution_count": null
@@ -198,11 +198,11 @@
    "file_extension": ".jl",
    "mimetype": "application/julia",
    "name": "julia",
-   "version": "1.11.1"
+   "version": "1.11.2"
   },
   "kernelspec": {
    "name": "julia-1.11",
-   "display_name": "Julia 1.11.1",
+   "display_name": "Julia 1.11.2",
    "language": "julia"
   }
  },

docs/src/examples/quantum1d/5.haldane-spt/index.md

@@ -4,4 +4,6 @@
 [quantum1d]
 "2.haldane" = "4f3e177a30df907cf7b9a47958ddf97741e4d0398ae1876e0f4977ccb93ae651"
 "3.ising-dqpt" = "f437ac4560db51f67e3997867c02bb8844b1077731b1dd1a06c61ab7c4846e97"
+"5.haldane-spt" = "6443696b20d416652951c646cf78b604a7fd0d79804995588491ea53ca210d57"
 "1.ising-cft" = "bd1253451973aa73c4fee6b64d81ee52cc2feba5b9632f87ce332eea3a05f460"
+"4.xxz-heisenberg" = "f8359ad29a028b01d41913014c0a442b99c3bd5cd403c23bd2a59f646a7c78cf"

docs/src/examples/quantum1d/5.haldane-spt/main.ipynb

@@ -15,14 +15,13 @@ Then we specify an initial guess, which we then further optimize.
 Working directly in the thermodynamic limit, this is achieved as follows:
 """
 
-H = heisenberg_XXX(; spin=1 // 2);
+H = heisenberg_XXX(; spin=1 // 2)
 
 md"""
 We then need an intial state, which we shall later optimize. In this example we work directly in the thermodynamic limit.
 """
 
-random_data = TensorMap(rand, ComplexF64, ℂ^20 * ℂ^2, ℂ^20);
-state = InfiniteMPS([random_data]);
+state = InfiniteMPS(2, 20)
 
 md"""
 The groundstate can then be found by calling `find_groundstate`.
@@ -58,9 +57,7 @@ md"""
 Let's initialize a different initial state, this time with a 2-site unit cell:
 """
 
-A = TensorMap(rand, ComplexF64, ℂ^20 * ℂ^2, ℂ^20);
-B = TensorMap(rand, ComplexF64, ℂ^20 * ℂ^2, ℂ^20);
-state = InfiniteMPS([A, B]);
+state = InfiniteMPS(fill(2, 2), fill(20, 2))
 
 md"""
 In MPSKit, we require that the periodicity of the hamiltonian equals that of the state it is applied to.
@@ -70,17 +67,17 @@ Alternatively, the hamiltonian can be constructed directly on a two-site unitcel
 
 ## H2 = repeat(H, 2); -- copies the one-site version
 H2 = heisenberg_XXX(ComplexF64, Trivial, InfiniteChain(2); spin=1 // 2)
-groundstate, cache, delta = find_groundstate(state, H2,
+groundstate, envs, delta = find_groundstate(state, H2,
                                              VUMPS(; maxiter=100, tol=1e-12));
 
 md"""
 We get convergence, but it takes an enormous amount of iterations.
 The reason behind this becomes more obvious at higher bond dimensions:
 """
 
-groundstate, cache, delta = find_groundstate(state, H2,
+groundstate, envs, delta = find_groundstate(state, H2,
                                              IDMRG2(; trscheme=truncdim(50), maxiter=20,
-                                                    tol=1e-12))
+                                                    tol=1e-12));
 entanglementplot(groundstate)
 
 md"""

examples/Cache.toml

@@ -37,13 +37,13 @@ casimir(s::SU2Irrep) = s.j * (s.j + 1)
 function heisenberg_hamiltonian(; J=-1.0)
     s = SU2Irrep(1)
     ℋ = SU2Space(1 => 1)
-    SS = TensorMap(zeros, ComplexF64, ℋ ⊗ ℋ ← ℋ ⊗ ℋ)
+    SS = zeros(ComplexF64, ℋ ⊗ ℋ ← ℋ ⊗ ℋ)
     for (S, data) in blocks(SS)
         data .= -0.5J * (casimir(S) - casimir(s) - casimir(s))
     end
-    return MPOHamiltonian(SS)
+    return InfiniteMPOHamiltonian(SS)
 end
-H = heisenberg_hamiltonian();
+H = heisenberg_hamiltonian()
 
 md"""
 ## Symmetry-Protected Topological Order