added spin projection in Kondo Hamiltonian

abhirup-m · Nov 8, 2024 · 1185f24 · 1185f24
1 parent 21313d8
commit 1185f24
Showing 1 changed file with 8 additions and 0 deletions.
diff --git a/src/modelHamiltonians.jl b/src/modelHamiltonians.jl
@@ -6,6 +6,10 @@ function KondoModel(
     )
     hamiltonian = Tuple{String, Vector{Int64}, Float64}[]
 
+    # project onto spin sector of impurity
+    largestScale = maximum(abs.[kondoJ, hop_t, globalField])
+    append!(hamiltonian, [("n", [1], -1e2 * largestScale), ("n", [2], -1e2 * largestScale), ("nn", [1, 2], 2e2 * largestScale)])
+
     # intra-bath hopping
     for site in 1:(numBathSites-1)
         push!(hamiltonian, ("+-",  [1 + 2 * site, 3 + 2 * site], -hop_t)) # c^†_{j,up} c_{j+1,up}
@@ -44,6 +48,10 @@ function KondoModel(
     numBathSites = length(dispersion)
     hamiltonian = Tuple{String, Vector{Int64}, Float64}[]
 
+    # project onto spin sector of impurity
+    largestScale = maximum(abs.[kondoJ, dispersion..., globalField])
+    append!(hamiltonian, [("n", [1], -1e2 * largestScale), ("n", [2], -1e2 * largestScale), ("nn", [1, 2], 2e2 * largestScale)])
+
     # kinetic energy
     for site in 1:numBathSites
         push!(hamiltonian, ("n",  [1 + 2 * site], dispersion[site])) # up spin